复合材料层合圆柱壳的有限元分析

本文根据Reissner-Mindlin型的全局位移场(一阶和三阶)，应用有限元预测一修正法，数值计算和分析了机械载荷作用下复合材料层合圆柱壳的挠度和横向剪应力。首先按照一般的有限元分析过程(没有引入剪切修正系数)计算出层合圆柱壳的挠度预测值；然后利用Lagrange插值构造横向剪应力的一般形式，使得满足层间连续和表面上为零的条件，通过最小二乘法拟合三维应力平衡方程获得横向剪应力；最后在单元上计算和引入剪切修正系数，再经过有限元分析计算出层合圆柱壳的挠度修正值。数值计算结果与三维线弹性解的比较表明，挠度修正值和横向剪应力的精度是十分满意的。     

一种分析对称铺设层合板的新方法

本文对多层纤维对称铺设层合板的分析提出了一种新方法。此方法以假设板厚方向横向剪应力分布为前提,来建立层合板的微分方程组。其特点是反映了层合板间剪应变跳跃式变化这一实际情况,满足层间剪应力平衡条件。文章通过算例说明此法对分析层间剪切刚度较低的层合板简便、有效。     

复合材料层合壳有限元分析的预测-修正法

对于复合材料层合壳的有限元分析，本文根据Reissner-Mindlin型的全局位移场给出了一个预测一修正法。首先按照一般的有限元分析过程(没有引入剪切修正系数)计算出全局响应(如挠度，频率和屈曲载荷等)的预测值；然后利用Lagrange插值构造横向剪应力的一般形式，使得满足层间连续和表面上为零的条件，通过最小二乘法拟合三维应力平衡方程获得横向剪应力；最后在单元上计算和引入剪切修正系数，再经过有限元分析计算出全局响应的修正值。     

粘弹层合板的稳态振动和层间应力

应用混合分层理论和Ressiner混合变分原理，在板厚方向取二次位移插值函数和三次、四次横向应力插值函数推导出粘弹层合板的动力学方程，得出简支粘弹层合板稳态振动的解。不仅得出与三层弹性板精确的自振频率吻合良好的解，而且对于粘弹层合板，所计算的自振频率和结构损耗因子也与三维结果吻合较好。计算了自由阻尼层合板对应的低阶法向位移响应幅值和层问横向应力的幅值。结果表明，较高的层间横向正应力是低频稳态振动中引起粘弹层合板分层破坏的主要因素，采用适当模量和厚度的粘弹性材料将有效地降低粘弹层合板的层间横向正应力的幅值。     

一种可以计算声波的反射和透射的层合板模型

本文给出一种层合板的模型,可以用来同时计算声波在层合板上的反射与透射,并且是严格满足层间位移和横向剪应力连续条件,然后给出一系列数值例子,与相应的准确解比较。以及说明在不同情况下声反射和透射的一些特点。     

新型厚薄通用的压电层合板三角形CDST-S6E单元

基于压电复合材料层合板一阶剪切变形理论及叠层理论,构造了一种新型三角形三节点压电层合板单元,简记为CDST-S6E单元.该单元采用压电耦合的运动方程求解位移场及电势场,层合板主体结构用一阶剪切变形理论模拟,其剪应变场及单元转角场由结点包含有两个剪切自由度的DST-S6单元理论确定,电势作为附加自由度,应用叠层理论对压电层合板的电势场沿厚度方向进行线性插值.该CDST-S6E单元不需要借助减缩积分、假设应力或应变等辅助数学手段,也不会产生对稳定性带来影响的附加零能模式,可较好避免厚薄板单元的剪切闭锬问题且具有简洁的表达形式.数值算例表明,CDST-S6E单元具有较高的精度,可以较为精确地预测压电层合板的变形及电势场,是一种厚薄通用的优质压电层合板单元.     

???????????????????????????

主要基于三维弹性力学和状态空间法给出了四边简支各向同性矩形层合板自由振动和 强迫振动问题的精确解. 首先基于三维弹性力学建立了层合板的基本方程,利用状态空间法 解决了层合板的自由振动问题,然后根据Lagrange动力学方程求解了层合板受到横向冲击时 的强迫振动响应.     

具有自由边的复合材料层合板的解析解

从三维弹性力学基本方程出发,通过假设自由边的边界位移函数,建立了正交异性层合板的状态方程,给出了对边自由,对边简支矩形板的解析解.此解满足层合板的基本方程和层间连续条件.用本文的方法比较容易处理层合板的自由边.算例表明,数值结果具有较高的精度.     

基于Legendre多项式的双压电层合板三维压电谱元法模拟

谱单元作为一种高阶单元具有计算效率高和精度高的特点。本文在基于Legendre正交多项式的三维谱单元基础上提出了三维压电谱单元模型,用于压电层合板静力和动力性能模拟研究。压电层合板结构中的位移和电势自由度均离散到三维压电谱单元Gauss-Lobatto-Legendre(GLL)配置节点上,并且未对电势沿压电层厚度方向上的变化做任何假定。通过提高谱单元中沿压电层合板厚度方向上的形函数阶数的方法,来削弱三维谱单元在模拟薄板结构中出现的剪切闭锁现象。为验证单元的计算精度,取双压电层合板结构进行静力和动力行为模拟验证,并将计算结果与现有文献中的其它单元模型及有限元结果进行对比。结果表明,压电谱单元可有效模拟压电层合板的静力和动力行为,且提高谱单元形函数阶数可提高数值模拟精度。     

复合材料层合板粘弹性固化残余应力分析

基于Ｓｃｈａｐｅｒｙ积分型粘弹性本构关系，推导了考虑横向剪切效应的复合材料层合板线性热粘弹性有限元分析列式，对层合板的粘弹性响应和加工成型过程中的残余应力进行了分析，给出了一些有意义的结果。     

An advanced higher-order theory for laminated composite plates with general lamination angles

This paper proposes a higher-order shear deformation theory to predict the bending response of the laminated composite and sandwich plates with general lamination configurations.The proposed theory a priori satisfies the continuity conditions of transverse shear stresses at interfaces.Moreover,the number of unknown variables is independent of the number of layers.The first derivatives of transverse displacements have been taken out from the inplane displacement fields,so that the C 0 shape functions are only required during its finite element implementation.Due to C 0 continuity requirements,the proposed model can be conveniently extended for implementation in commercial finite element codes.To verify the proposed theory,the fournode C 0 quadrilateral element is employed for the interpolation of all the displacement parameters defined at each nodal point on the composite plate.Numerical results show that following the proposed theory,simple C 0 finite elements could accurately predict the interlaminar stresses of laminated composite and sandwich plates directly from a constitutive equation,which has caused difficulty for the other global higher order theories.     

A new zig-zag theory and C<Superscript>0</Superscript> plate bending element for composite and sandwich plates

A higher-order zig-zag theory for laminated composite and sandwich structures is proposed. The proposed theory satisfies the interlaminar continuity conditions and free surface conditions of transverse shear stresses. Moreover, the number of unknown variables involved in present model is independent of the number of layers. Compared to the zig-zag theory available in literature, the merit of present theory is that the first derivatives of transverse displacement have been taken out from the in-plane displacement fields, so that the C⁰ interpolation functions is only required during its finite element implementation. To obtain accurately transverse shear stresses by integrating three-dimensional equilibrium equations within one element, a six-node triangular element is employed to model the present zig-zag theory. Numerical results show that the present zig-zag theory can predict more accurate in-plane displacements and stresses in comparison with other zig-zag theories. Moreover, it is convenient to obtain transverse shear stresses by integration of equilibrium equations, as the C⁰ shape functions is only used when implemented in a finite element.     

A C<Superscript>0</Superscript>-type zig–zag theory and finite element for laminated composite and sandwich plates with general configurations

In order to conveniently develop C⁰ continuous element for the accurate analysis of laminated composite and sandwich plates with general configurations, this paper develops a C⁰-type zig–zag theory in which the interlaminar continuity of transverse shear stresses is a priori satisfied and the number of unknowns is independent of the number of layers. The present theory is applicable not only to the cross-ply but also to the angle-ply laminated composite and sandwich plates. On the premise of retaining the merit of previous zig–zag theories, the derivatives of transverse displacement have been taken out from the displacement fields. Therefore, based on the proposed zig–zag theory, it is very easy to construct the C⁰ continuous element. To assess the performance of the proposed model, the classical quadratic six-node triangular element with seven degrees of freedom at each node is presented for the static analysis of laminated composite and sandwich plates. The typical examples are taken into account to assess the performance of finite element based on the proposed zig–zag theory by comparing the present results with the three-dimensional elasticity solutions. Numerical results show that the present model can produce the more accurate deformations and stresses compared with the previous zig–zag theories.     

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

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

A layerwise C0-type higher order shear deformation theory for laminated composite and sandwich plates

A novel layerwise C⁰-type higher order shear deformation theory (layerwise C⁰-type HSDT) for the analysis of laminated composite and sandwich plates is proposed. A C⁰-type HSDT is used in each lamina layer and the continuity of in-plane displacements and transverse shear stresses at inner-laminar layer is consolidated. The present layerwise theory retains only seven variables without increasing the number of variables when the number of lamina layers are intensified. The shear stresses through the plate thickness derived from the constitutive equation of the present theory have the same shape as those calculated from the equilibrium equation. In addition, the artificial constraints are added in the principle of virtual displacements (PVD) and are certainly fulfilled through a penalty approach. In this paper, two C⁰-continuity numerical methods, such as the Finite Element Method (FEM) and Bézier isogeometric element (BIEM) are utilized to solve a discrete system of equations derived from the PVD. Several numerical examples with various geometries, aspect ratios, stiffness ratios, and boundary conditions are investigated and compared with the 3D elasticity solution, the analytical, as well as, numerical solutions based on various plate theories.     

Accurate transverse stress evaluation in composite/sandwich thick laminates using a C0 HSDT and a novel post-processing technique

An efficient method for accurate evaluation of through-the-thickness distribution of transverse stresses in thick composite and sandwich laminates, using a displacement-based C⁰ higher-order shear deformation theory (HSDT), is presented. The technique involves a least square of error (LSE) method applied to the 3D equilibrium equations at the post-processing phase, after a primary finite element analysis is performed using the HSDT. This is distinctly different from the conventional method of integrating the 3D equilibrium equations, for transverse stress recovery in composite laminates during post-processing. Competence of the technique is demonstrated in the numerical examples through comparison with results from first-order shear deformation theory (FSDT), another HSDT and those from analytical and 3D elasticity solutions available in literature.     

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.     

Buckling of soft-core sandwich plates with angle-ply face sheets by means of a C0 finite element formulation

In order to avoid using C¹ interpolation functions in finite element implementation of the previous zig–zag theories, artificial constraints, in which the first derivatives of transverse displacement will be replaced by the assumed variables, are usually employed. However, such assumption will violate continuity conditions of transverse shear stresses at interfaces. Differing from previous work, this paper will propose a C⁰-type zig–zag theory for buckling analysis of laminated composite and sandwich plates with general configurations. The first derivatives of transverse displacement have been taken out from a displacement field of the proposed zig–zag theory. Thus, the C⁰ interpolation functions are only required in finite element implementations of the proposed model. Without use of any artificial constraints, an eight-node quadrilateral element based on the proposed model is presented by incorporating the terms associated with the geometric stiffness matrix. In order to verify performance of the proposed model, several buckling problems of sandwich plates with soft core have been analyzed. Numerical results show that the proposed model is able to predict accurately buckling loads of the soft-core sandwich plates with varying fiber orientations of face sheets.     

Free vibration of functionally graded sandwich plates using four-variable refined plate theory

This paper uses the four-variable refined plate theory (RPT) for the free vibration analysis of functionally graded material (FGM) sandwich rectangular plates.Unlike other theories, there are only four unknown functions involved, as compared to five in other shear deformation theories. The theory presented is variationally consistent and strongly similar to the classical plate theory in many aspects. It does not require the shear correction factor, and gives rise to the transverse shear stress variation so that the transverse shear stresses vary parabolically across the thickness to satisfy free surface conditions for the shear stress. Two common types of FGM sandwich plates are considered, namely, the sandwich with the FGM facesheet and the homogeneous core and the sandwich with the homogeneous facesheet and the FGM core. The equation of motion for the FGM sandwich plates is obtained based on Hamilton's principle. The closed form solutions are obtained by using the Navier technique. The fundamental frequencies are found by solving the eigenvalue problems. The validity of the theory is shown by comparing the present results with those of the classical, the first-order, and the other higher-ordex theories. The proposed theory is accurate and simple in solving the free vibration behavior of the FGM sandwich plates.