Construction of rectangular displacement-based element for thick/thin plates

In this paper, a method of constructing displacement-based element for thick/thin plates is developed by using the technique of generalized compatibility, and a rectangular displacement based element with 12 degrees of freedom for thick/thin plates is presented. This method enjoys a good accuracy with simple formulation and is free of shear locking as the thickness of the plate approaches zero. The project supported by National Natural Science Foundation of China through Grant No. 59208075     

一种薄壳厚壳通用的轴对称曲壳单元

有限元法分析轴对称壳体时，常需区分薄壳和厚壳并选用不同的单元，给计算带来不便．为此，通过对现有的轴对称壳体单元的研究，基于加权残值法，将广义协调条件引入剪应变场，构造了一种挠度和转角各自独立的新型轴对称曲壳单元．算例结果表明，新单元具有很高的精度，既可用于厚壳，也可用于薄壳．该单元可用于轴对称壳体结构的计算分析．     

A NEW QUADRILATERAL THIN PLATE ELEMENT BASED ON THE MEMBRANE-PLATE SIMILARITY THEORY

ＩｎｔｒｏｄｕｃｔｉｏｎＢｅｃａｕｓｅｏｆｔｈｅｒｅｑｕｉｒｅｍｅｎｔｏｆｃ1ｃｏｎｔｉｎｕｉｔｙ ,ｉｔｉｓｖｅｒｙｄｉｆｆｉｃｕｌｔｔｏｃｏｎｓｔｒｕｃｔｃｏｎｆｏｒｍｉｎｇＫｉｒｃｈｈｏｆｆｐｌａｔｅｂｅｎｄｉｎｇｅｌｅｍｅｎｔｓ.Ｉｎｏｒｄｅｒｔｏｏｖｅｒｃｏｍｅｔｈｅｄｉｆｆｉｃｕｌｔｙ ,ｍａｎｙａｐｐｒｏａｃｈｅｓｈａｖｅｂｅｅｎｐｒｅｓｅｎｔｅｄ .Ｉｎｔｈｅｓｅａｐｐｒｏａｃｈｅｓ,ｔｈｅｒｅｑｕｉｒｅｍｅｎｔｏｆｃ1ｉｓｒｅｌｅａｓｅｄｏｎｃｏｎｄｉｔｉｏｎｏｆｔｈｅｃｏｎｖｅｒｇｅｎ…     

复合材料层合厚板锯齿理论及有限元分析

对于较厚复合材料弯曲问题,已有锯齿型厚板理论最大误差超过35%。为了合理地分析较厚复合材料弯曲问题,发展了准确高效的锯齿型厚板理论。此理论位移变量个数独立于层合板层数,其面内位移不含有横向位移一阶导数,构造有限元时仅需C0插值函数,故称此理论为C0型锯齿厚板理论。基于发展的锯齿理论,构造了六节点三角形单元并推导了复合材料层合/夹层板弯曲问题有限元列式。为验证C0型锯齿厚板理论性能,分析了复合材料层合/夹层厚板弯曲问题,并与已有C1型锯齿理论对比。结果表明,本文的C0型锯齿厚板理论最大误差15%,比已有锯齿型厚板理论准确高效。     

General non-linear finite element analysis of thick plates and shells

A non-linear finite element analysis is presented, for the elasto-plastic behavior of thick shells and plates including the effect of large rotations. The shell constitutive equations developed previously by the authors [Voyiadjis, G.Z., Woelke, P., 2004. A refined theory for thick spherical shells. Int. J. Solids Struct. 41, 3747–3769] are adopted here as a base for the formulation. A simple C⁰ quadrilateral, doubly curved shell element developed in the authors’ previous paper [Woelke, P., Voyiadjis, G.Z., submitted for publication. Shell element based on the refined theory for thick spherical shells] is extended here to account for geometric and material non-linearities. The small strain geometric non-linearities are taken into account by means of the updated Lagrangian method. In the treatment of material non-linearities the authors adopt: (i) a non-layered approach and a plastic node method [Ueda, Y., Yao, T., 1982. The plastic node method of plastic analysis. Comput. Methods Appl. Mech. Eng. 34, 1089–1104], (ii) an Iliushin’s yield function expressed in terms of stress resultants and stress couples [Iliushin, A.A., 1956. Plastichnost’. Gostekhizdat, Moscow], modified to investigate the development of plastic deformations across the thickness, as well as the influence of the transverse shear forces on plastic behaviour of plates and shells, (iii) isotropic and kinematic hardening rules with the latter derived on the basis of the Armstrong and Frederick evolution equation of backstress [Armstrong, P.J., Frederick, C.O., 1966. A mathematical representation of the multiaxial Bauschinger effect. (CEGB Report RD/B/N/731). Berkeley Laboratories. R&D Department, California.], and reproducing the Bauschinger effect. By means of a quasi-conforming technique, shear and membrane locking are prevented and the tangent stiffness matrix is given explicitly, i.e., no numerical integration is employed. This makes the current formulation not only mathematically consistent and accurate for a variety of applications, but also computationally extremely efficient and attractive.     

A NEW HYBRID QUADRILATERAL FINITE ELEMENT FOR MINDLIN PLATE

ＡＮＥＷＨＹＢＲＩＤＱＵＡＤＲＩＬＡＴＥＲＡＬＦＩＮＩＴＥＥＬＥＭＥＮＴＦＯＲＭＩＮＤＬＩＮＰＬＡＴＥＣｈｉｎＹｉ（秦奕）（ＴｉａｎｊｉｎＡｒｃｈｉｔｅｃｔｕｒａｌＤｅｓｉｇｎＩｎｓｔｉｔｕｔｅ，Ｔｉａｎｊｉｎ）ＺｈａｎｇＪｉｎｇ－ｙｕ（张敬宇）（Ｉ...     

Boundary element method for orthotropic thick plates

The fundamental solutions of the orthotropic thick plates taking into account the transverse shear deformation are derived by means of Hörmander's operator method and a plane-wave decomposition of the Dirac δ-function in this papey. The boundary integral equations of the thick plates have been formulated which are adapted to arbitrary boundary conditions and plane forms. The numerical calculation of the fundamental solutions is discussed in detail. Some numerical examples are analyzed with BEM.     

各向异性弹塑性中厚度板壳的有限元分析

本文在文献［２，３］的基础上，提出了一个解各向异性弹塑性中厚度板壳问题的有限元方法。考虑材料各向异性的特点，采用了Ｈｉｌｌ推广的Ｈｕｂｅｒ－Ｍｉｓｅｓ屈服准则；借用Ｏｗｅｎ的剪切修正系数，正确计及了叠层复合材料壳体的横向剪切效应；为了避免“自锁”现象，文中采用了９节点的Ｈｅｔｅｒｏｓｉｓ二次壳单元；特别是本文利用插值外推的思想，提出了一个带预测的弧长增量控制法，显著提高了确定变形路径的计算效率。几个数值算例表明本文给出的有限元方法对于各向异性中厚度板壳的弹塑性分析有较好的精度，尤其是对具有复杂变形路径的结构计算，收敛速度提高更快。     

A hybrid/mixed model finite element analysis for eigenvalue problem of moderately thick plates

The buckling and free vibration problems of moderately thick plate are considered in this paper by using the hybrid/mixed finite element model. A modified Reissner principle which only requires C⁰ continuity is derived. No lockling phenomenon is observed. Linear interpolation is used for all independent unknown function. Finally a displacement generalized eigenvalue equation is obtained, in which the stiffness matrix is symmetric and positively definite. The calculated results show that the method proposed is simple, reliable and satisfactory.     

中厚板统计分析的随机边界元法

本文建立了分析含随机材料参数并具厚度不均匀性的中厚板问题的随机边界元法,基于Taylor级数展开技术,分析和到广义位移的均值和一阶偏差的积分方程,其中将材料参数的随机性和厚度的不均匀性作为等效荷载处理,从而得到广义边界位移或面力的均值和协方差,并进一步求出部点广义位移和内力的均值和协方差,最后用本文方法计算了两个数例,并对所得结果进行了分析,探讨。     

一种适用于简支多边形薄板问题的有限单元法

简支多边形薄板双调和定解问题可分解为二个互不耦合的Ｐｏｉｓｓｏｎ方程的定解问题，通过比拟，可用二维平面应力（应变）单元进行有限元分析。数值算例表明，本文采用的方法在较粗的网格下即可得到较高的精度，而且具有速度快、占用计算机内存少等优点。     

薄板几何非线性中的精化元方法及膜闭锁问题

基于放松单元间协调条件的大变形变分原理和全局拉格朗日方法，推导了几何非线性精化三角形薄板单元。对几何刚度矩阵，通过引入特殊的单元位移函数，有效地消除了薄板弯曲问题中伴生的膜闭锁现象。数值结果表明该单元在几何非线性分析中既能消除膜闭锁又具有较高精度。     

求解不连续中厚板自由振动的微分容积单元法

基于区域叠加原理和微分容积法,发展了一种新型的数值方法——微分容积单元法,用以分析具有不连续几何特征的中厚板的自由振动。根据板的不连续情况将其划分为若干单元,在每个单元内用微分容积法将控制微分方程离散成为一组线性代数方程．在相邻的单元连接处应用位移连续条件和平衡条件,引入边界约束条件后得到一套关于各配点位移的齐次线性代数方程,由此可导出求解系统固有频率的特征方程。本文用子空间迭代法求解特征方程,并以开孔板、混合边界条件板和突变厚度板为例研究了方法的收敛性和计算精度。     

Optimizing formulation of hybrid model for thin and moderately thick plates

A 20 — DOF hybrid stress element based upon Mindlin plate theory is developed using the optimization design method for thin and moderately thick plates. Numerical tests consist of the convergency and performance to the plates with arbitrary thickness and shape and of the ultimate thin plate problems.Projects Supported by the National Natural Science Foundation of China.     

Static and stability analyses of arbitrary straight-sided quadrilateral thin plates by DQM

A differential quadrature (DQ) methodology is employed for the static and stability analysis of irregular quadrilateral straight-sided thin plates. A four-noded super element is used to map the irregular physical domain into a square computational domain. Second order transformation schemes with relative ease and low computational effort are employed to transform the fourth order governing equations of thin plates between the domains. Within the domain, the displacements are the only degrees of freedom whereas, along the boundaries, the displacements as well as the second order derivatives of the displacements with respect to the associated normal coordinate variables in the computational domain are the two sets of degrees of freedom. The implementation procedures for different boundary conditions including free-edge boundaries are formulated. To demonstrate the accuracy, convergency and stability of the methodology, detailed studies of skewed and trapezoidal plates for different geometries under different boundary and loading conditions are made. Good agreement is achieved between the results of the present methodology and those of other DQ methodologies or other comparable numerical algorithms.     

A quadrilateral element based on refined global-local higher-order theory for coupling bending and extension thermo-elastic multilayered plates

In this paper a refined higher-order global-local theory is presented to analyze the laminated plates coupled bending and extension under thermo-mechanical loading. The in-plane displacement fields are composed of a third-order polynomial of global coordinate z in the thickness direction and 1,2–3 order power series of local coordinate ζ_k in the thickness direction of each layer, which is identical to the 1,2–3 global-local higher-order theory by Li and Liu [Li, X.Y., Liu, D., 1997. Generalized laminate theories based on double superposition hypothesis. Int. J. Numer. Methods Eng. 40, 1197–1212] Moreover, a second-order polynomial of global coordinate z in the thickness direction is chosen as transverse displacement field. The transverse shear stresses can satisfy continuity at interfaces, and the number of unknowns does not depend on the layer numbers of the laminate.Based on this theory, a quadrilateral laminated plate element satisfying the requirement of C¹ continuity is presented. By solving both bending and thermal expansion problems of laminates, it can be found that the present refined theory is very accurate and obviously superior to the existing 1,2–3 global-local higher-order theory. The most attractive feature of this theory is that the transverse shear stresses can be accurately predicted from direct use of constitutive equations without any post-processing method. It is also shown that the present quadrilateral element possesses higher accuracy.     

Variational principles of perforated thin plates and finite element method for buckling and post-buckling analysis

On the basis of the general theory of perforated thin plates under large deflections^{[1, 2]}, variational principles with deflectionw and stress functionF as variables are stated in detail. Based on these principles, finite element method is established for analysing the buckling and post-buckling of perforated thin plates. It is found that the property of element is very complicated, owing to the multiple connexity of the region. Project supported by National Natural Science Foundation of China.     

A new approach to an analysis of large deflections of thin elastic plates

For thin plates undergoing large deflections, a modified energy expression has been suggested and a new set of differential equations has been obtained in a decoupled form. Accuracy of the equations has been tested for a circular and a square plate with immovable as well as movable edge conditions under a uniform static load. Results obtained are in excellent agreement with other known results. These new equations are more advantageous than Berger's decoupled equations which fail to give meaningful results for movable edge conditions.     

An exact element method for bending of nonhomogeneous thin plates

In this paper, based on the step reduction method, a new method, the exact element method for constructing finite element, is presented. Since the new method doesn't need the variational principle, it can be applied to solve non-positive and positive definite partial differential equations with arbitrary variable coefficient. By this method, a triangle noncompatible element with 6 degrees of freedom is derived to solve the bending of nonhomogeneous plate. The convergence of displacements and stress resultants which have satisfactory numerical precision is proved. Numerical examples are given at the end of this paper, which indicate satisfactory results of stress resultants and displacements can be obtained by the present method.     

On the analysis of thick rectangular plates

Summary Thick rectangular plates are investigated using the method of initial functions proposed by Vlasov. The governing equations are derived from the three-dimensional elasticity equations using a MacLaurin series approach. As the governing equations can be obtained in the form of series, approximate theories of any desired order can be constructed easily by proper truncation. An exact solution is obtained for an allround simply supported thick plate using a Navier type solution. A Levy type solution for higher order theories is illustrated for the case of a thick plate with two opposite edges simply supported and other two edges clamped. Numerical results obtained are compared with those of classical, Reissner and Srinivas et al. solutions.

---

Übersicht Mit Hilfe der Methode der Initial-Funktionen von Vlasov werden rechteckige Platten untersucht. Die zugehörigen Gleichungen werden aus den Gleichungen für das dreidimensionale Problem durch eine Entwicklung in MacLaurin-Reihen gewonnen. Durch Abbrechen dieser Reihen können Näherungen beliebiger Ordnung erhalten werden. Für den Fall einer allseitig einfach gelagerten dicken Platte wird eine exakte Lösung erhalten, bei der eine Lösung vom Navier-Typ verwendet wird. Eine Lösung vom Levy-Typ höherer Ordnung wird am Beispiel einer dicken Platte abgeleitet, von der zwei gegenüberliegende Ecken einfach gelagert, die anderen fest eingespannt sind. Die numerischen Ergebnisse werden mit den klassischen, von Reissner, Srinivas u. a. erhaltenen Resultaten verglichen.