A FINITE DIFFERENCE METHOD AT ARBITRARY MESHES FOR THE BENDING OF PLATES WITH VARIABLE THICKNESS

A finite difference method at arbitrary meshes for the bending of plates with variable thickness is presented in this paper. The method is completely general with respect to various boundary conditions, load cases and shapes of plates. This difference scheme is simple and the numerical results agree well with those obtained by other methods.

RATIONAL FINITE ELEMENT METHOD FOR ELASTIC BENDING OF REISSNER PLATES

In this paper; some deformation patterns defined by differential equations of the elastic system are introduced into the revised functional for the incompatible elements. And therefore the rational FEM, which is perfect combination of the analytic methods and numeric methods, has been presented. This new approach satisfies not only the mechanical requirement of the elements but also the geometric requirement of the complex structures. What's more the quadrilateral element obtained accordingly for the elastic bending of thick plates demonstrates such advantages as high precision for computation and availability of accurate integration for stiffness matrices.

Finite element analysis of elasto-plastic plate bending problems using transition rectangular plate elements

In this work, the finite element analysis of the elasto-plastic plate bending problems is carried out using transition rectangular plate elements. The shape functions of the transition plate elements are derived based on a practical rule. The transition plate elements are all quadrilateral and can be used to obtain efficient finite element models using minimum number of elements. The mesh convergence rates of the models including the transition elements are compared with the regular element models. To verify the developed elements, simple tests are demonstrated and various elasto-plastic problems are solved. Their results are compared with ANSYS results.

SPECIAL HIGH-ORDER BENDING CRACK TIP FINITE ELEMENT FOR THE REISSNER PLATE

Based on the crack tip field expansion of the Reissner plate, a special high order bending crack tip element is developed, and the element stiffness matrix is given in the explicit form, which is especially convenient for engineering analyses. A numerical example is presented and compared with previous results to demonstrate the efficiency and accuracy of the special element.

A Comparative Study of Stress Recovery Method and Error Estimation of Plate Bending Problem Using DKMQ Element

Recovery by Equilibrium in Patches (REP) is a recovery method introduced by B. Boroomand. This method is using patch as recovery media as is used by Superconvergent Patch Recovery (SPR) which is well known as a good recovery method. In this research, a numerical study of REP implementation is held to estimate error in finite element analysis using DKMQ element. The numerical study is performed with both uniform and adaptive h-type mesh refinement. The result is compared with three other recovery method, i.e. SPR method, averaging method, and projection method.

EXPERIMENT AND SIMULATION FOR THICK-PLATE BENDINGBY HIGH FREQUENCY INDUCTOR

In order to produce thick plates with complicated curved surface, a prototype bending machine by the use of high frequency inductor was developed. The bending mechanism is based on the localized stresses which are induced from the difference of temperature in thickness by the high frequency inductor. The operating speed and the thickness of plate were examined from the experiment, and the variation of the temperature was measured. Finite element analysis was carried out in the second part based on the experimentally obtained temperature distribution. The so-called Mindlin plate element was used in order to perform the simulation efficiently. The strategy to produce such curved surface in the practical process was discussed and further perspective of the production system was described. (Edited author abstract) 6 Refs.

P-FEMOL ANALYSIS OF PLATE BENDING PROBLEMS

In this paper, the p- version of the finite element method of lines (FEMOL) for the analysis of the Mindlin-Reissner plate bending problems is presented and a class of p-FEMOL elements with polynomial degrees as high as nine is developed. Numerical examples given in this paper show tremendous performance of the present method: namely, rapid convergence rate, high accuracy for both displacements and stress resultants, removal of shear-locking trouble, capability of dealing with difficult problems such as the boundary layer behavior near a free edge and stress concentration around a hole.

Static and free vibration analyses of functionally graded porous variable-thickness plates using an edge-based smoothed finite element method

The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES-FEM) associate with the mixed interpolation of tensorial components technique for the three-node triangular element(MITC3), so-called ES-MITC3. This ES-MITC3 element is performed to eliminate the shear locking problem and to enhance the accuracy of the existing MITC3 element. In the ES-MITC3 element, the stiffness matrices are obtained by using the strain smoothing technique over the smoothing domains formed by two adjacent MITC3 triangular elements sharing an edge. Materials of the plate are FGP with a power-law index(k) and maximum porosity distributions(U) in the forms of cosine functions. The influences of some geometric parameters, material properties on static bending, and natural frequency of the FGP variable-thickness plates are examined in detail.

B-Spline Wavelet on Interval Finite Element Method for Static and Vibration Analysis of Stiffened Flexible Thin Plate

A new wavelet finite element method(WFEM)is constructed in this paper and two elements for bending and free vibration problems of a stiffened plate are analyzed.By means of generalized potential energy function and virtual work principle,the formulations of the bending and free vibration problems of the stiffened plate are derived separately.Then,the scaling functions of the B-spline wavelet on the interval(BSWI)are introduced to discrete the solving field variables instead of conventional polynomial interpolation.Finally,the corresponding two problems can be resolved following the traditional finite element frame.There are some advantages of the constructed elements in structural analysis.Due to the excellent features of the wavelet,such as multi-scale and localization characteristics,and the excellent numerical approximation property of the BSWI,the precise and efficient analysis can be achieved.Besides,transformation matrix is used to translate the meaningless wavelet coefficients into physical space,thus the resolving process is simplified.In order to verify the superiority of the constructed method in stiffened plate analysis,several numerical examples are given in the end.

ANALYSIS OF INTER-LAMINAR STRESSES FOR COMPOSITE LAMINATED PLATE WITH INTERFACIAL DAMAGE

A constitutive model for composite laminated plates with the damage effect of the intra-layers and inter-laminar interface is presented. The model is based on the general six-degrees-of-freedom plate theory, the discontinuity of displacement on the interfaces are depicted by three shape functions, which are formulated according to solutions satisfying three equilibrium equations, By using the variation principle, the three-dimensional non-linear equilibrium differential equations of the laminated plates with two different damage models are derived. Then, considering a simply supported laminated plate with damage, an analytical solution is presented using finite difference method to obtain the inter-laminar stresses.

TRAPEZOIDAL PLATE BENDING ELEMENT WITH DOUBLE SET PARAMETERS

Using double set parameter method, a 12-parameter trapezoidal plate bending element is presented. The first set of degrees of freedom, which make the element convergent, are the values at the four vertices and the middle points of the four sides together with the mean values of the outer normal derivatives along four sides. The second set of degree of freedom, which make the number of unknowns in the resulting discrete system small and computation convenient are values and the first derivatives at the four vertices of the element. The convergence of the element is proved.

Study on Bending Conditions of Plate in Coil Box

In a coil box between the roughing and finishing stands on a hot strip mill,a problem has been encountered that the entry region of the plate touches the bending rolls and deforms.As a result,the defective coil occurs.The condition of plate bending,which forms a new deformation feature in coiling,is analyzed.In this paper,the authors focus on the research of the effects of coiling parameters,such as the thickness of plate,roll speed and feeding speed of plate in coil box,and on specific plate bending.A finite element method is developed to simulate this coiling process.Based on numerical simulation,the effects of the coiling parameters on the mechanics and deformation of the bending plate are obtained.Numerical simulation tests have verified the validity of the developed model.

BUCKLING ANALYSIS UNDER COMBINED LOADING OF THIN-WALLED PLATE ASSEMBLIES USING BUBBLE FUNCTIONS

Bubble functions are finite element modes that are zero on the boundary of the element but nonzero at the other point. The present paper adds bubble functions to the ordinary Complex Finite Strip Method(CFSM) to calculate the elastic local buckling stress of plates and plate assemblies. The results indicate that the use of bubble functions greatly improves the convergence of the Finite Strip Method(FSM) in terms of strip subdivision, and leads to much smaller storage required for the structure stiffness and stability matrices. Numerical examples are given, including plates and plate structures subjected to a combination of longitudinal and transverse compression, bending and shear. This study illustrates the power of bubble functions in solving stability problems of plates and plate structures.

BENDING PROBLEM OF A FINITE COMPOSITE LAMINATED PLATE WEAKENED BY MULTIPLE ELLIPTICAL HOLES

Based on the classical composite laminate theory,the bending problem of a finite composite plate weakened by multiple elliptical holes is studied by means of the complex variable method.The present work is intended to express the complex potentials in the form of Faber series aided by the use of the least squares boundary collocation techniques on the finite boundaries.As a result,concise and high accuracy solutions are presented for the stress distribution around the holes.Finally,numerical examples are presented to discuss the effects of some parameters on the stress concentration around the holes.

Study of Longitudinal Oscillations of a Building on the Basis of Continuum Plate Model

Continuum plate model in the form of a cantilever anisotropic plate developed in the framework of the bimoment theory of plates describing seismic oscillations of buildings is proposed in this paper as a dynamic model of a building. Formulas for the reduced moduli of elasticity, shear and density of the plate model of a building are given. Longitudinal oscillations of a building are studied using the continuum plate and box-like models of the building with Finite Element Model. Numerical results are obtained in the form of graphs, followed by their analysis.

Analysis of Critical Load of the Orthotropic Plate Structures

Advanced design based on the concept of orthotropic structure includes better use of materials, less weight compared to the equivalent isotropic construction and controlled effectively reserve resistance in all its segments. In this case a calculation of critical load is exposed using the FDM （Finite Difference Method） concept of thin plates subjected to complex loads due to forces in the middle-plane. Results of calculation model, discussed in this paper, are given in graphic form. Presented results should serve as an indicator of the expansion of theoretical base of similar models, which can be reasonably use by researchers and engineers in their practices, and by students for educational purposes.

Steady Natural Convection from a Vertical Hot Plate with Variable Radiation

The natural convection from a vertical hot plate with radiation and constant flux is studied numerically to know the velocity and temperature distribution characteristics over a vertical hot plate.The governing equations of the natural convection in two-dimension are solved with the implicit finite difference method,whereas the discretized equations are solved with the iterative relaxation method.The results show that the velocity and the temperature increase along the vertical wall.The influence of the radiation parameter in the boundary layer is significant in increasing the velocity and temperature profiles.The velocity profiles increase with the increase of the radiation parameter.The temperature profiles near the wall plate parallel each other due to the constant heat flux applied to the wall.The influence of the radiation parameter is significant either in velocity or temperature characteristics.At the same time,the effect of the Prandtl number greater than 0.71 is not sensitive to the velocity and temperature variations elsewhere.

黏弹性纳米流体在垂直板上的自然对流与传热分析

基于分数阶Maxwell模型和分数阶Fourier定律构建黏弹性纳米流体在垂直板上的非定常二维边界层自然对流与传热控制方程,利用有限差分和L1算法获得数值稳定解,对不同物理参数下的速度、温度、平均表面摩擦系数和平均Nusselt数的变化趋势进行图形化分析。结果显示,速度和温度边界层均表现出短暂记忆和延迟特性;速度分数导数参数削弱了自然对流,而速度松弛时间的影响却相反;温度分数导数参数削弱了自然对流和热传导,而温度松弛时间的影响却相反。

平行板微通道中一类不可压缩微极性流体在高Zeta势下的时间周期电渗流

在高Zeta势下,研究平行板微通道中一类不可压缩微极性流体的时间周期电渗流.在不使用DebyeHüickel线性近似条件下,利用有限差分法数值求解非线性Poisson-Boltzmann方程和不可压缩微极性流体的连续性方程、动量方程、角动量方程及本构方程,在低Zeta势下将所得结果与使用Debye-Hückel线性近似得到的解析解比较,证明本文数值方法是可行的;讨论高Zeta势下电动宽度m、电振荡频率Ω、微极性参数k1等无量纲参数对不可压缩微极性流体的速度和微旋转效应的影响.研究表明:1)随着Zeta势的增大,微极性流体的速度、微旋转、体积流量、微旋强度以及剪切应力增大,说明与低Zeta势相比,高Zeta势对微极性流体电渗流有显著的促进作用.2)在高Zeta势下,随着微极性参数的增大,微极性流体的速度减小,但是对微旋转效应呈现先增强后减弱的趋势.3)在高Zeta势下,当电振荡频率较低(小于1)时,电动宽度的增大促进微极性流体的流动,但抑制其微旋转;当电振荡频率较高(大于1)时,电动宽度的增大抑制微极性流体的流动及微旋转,但促进体积流量快速增大并趋于恒定.4)在高Zeta势下,当电振荡频率较低(小于1)时,微极性流体电渗流速度和微旋转随着电振荡频率的变化呈现明显的振荡变化趋势,但是速度和微旋转的峰值、体积流量及微旋强度均保持不变;当电振荡频率较高(大于1)时,随着电振荡频率的增大,微极性流体电渗流速度和微旋转的幅值减小,体积流量及微旋强度减小直至趋于零.5)在高Zeta势下,壁面剪切应力σ21及σ12的幅值随电动宽度的增大而增大;当电振荡频率较低(小于1)时,壁面剪切应力σ21与σ12不随电振荡频率的增大而变化,均取恒定值,且微极性参数的取值不影响壁面剪切应力σ21的幅值;当电振荡频率较高(大于1)时,壁面剪切应力σ21及σ12的幅值随电振荡频率的增大而减小,且壁面剪切应力σ21的幅值随着微极性参数的增大而减小,而壁面剪切应力σ12的振幅随着微极性参数的增大而线性减小.

Winkler地基上四边自由正交各向异性矩形中厚板弯曲的有限积分变换解

本文应用有限积分变换法研究Winkler地基上四边自由正交各向异性矩形中厚板的弯曲问题.具体由正交各向异性矩形中厚板弯曲的基本方程组和边界条件出发,结合有限积分变换法及其对应的逆变换法推导出正交各向异性矩形中厚板弯曲问题的解析解.该解析解统一适用于计算各向同性/正交各向异性矩形薄板、中厚板和厚板的弯曲问题,并且通过具体算例验证了所得解析解的正确性.