板结构计算模型的新发展

现代复合材料层合板具有高强和轻型的突出优点,从而在军工和民用等诸多领域发挥着重要作用。这种板结构的特点是随着纤维走向的不同,层间材料的物理-力学特性发生剧烈变化。沿板厚方向变形的梯度比较陡峭,并在层间结合面处发生强不连续,呈现zig-zag(锯齿状)现象。这导致横向剪应变在板的静态和动态响应中发生重要作用,不计横向变形的经典组合板计算模型CLPT难以适应现代多层板计算分析的需要。考虑横向剪切变形影响的板的计算模型得到重视和发展。需要指出,现有各种考虑剪切变形影响的计算模型虽然有了很大的发展,但在全面和准确性上仍然存在一定的不足,难以适应现代多层组合板横向力和物理性能多变的情况。模型预测的沿板厚方向位移和应力的变化规律难以通过严格的检验。本文提出的以比例边界有限元为基础的正交各向异性板的数值计算模型,同时可适用于各种薄板与厚板的分析,对现代复合材料层合板的分析具有特殊的优越性。所得到的板的位移、正应力和剪应力沿板厚方向的变化,与三维弹性理论的标准解高度吻合。数值算例进一步表明,随着层间纤维走向的变化,板内位移场和应力场沿板厚方向剧烈变化所呈现的锯齿现象均可以精准地进行模拟。据此,本文建议方法对现代板分析的广泛适应性和高度准确性得到了充分论证。     

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

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

计入横向剪切变形时层合板位移差分解法

在考虑横向剪切变形对层合板弹性解的影响时，本文提出一种数值计算方法。由边界条件给出边界结点位移的表达式，根据薄板的经典理论和一阶横向剪切变形理论导出位移增量所满足的平衡微分方程，引用经典理论计算的横向剪力修正了荷载列阵。致使在较粗的网格划分时、宽广的层合板长厚比范围内，仍能得到与解析解颇为一致的数值解。     

基于新修正偶应力理论的Mindlin层合板稳定性分析

基于新的各向异性修正偶应力理论提出一个Mindlin复合材料层合板稳定性模型。该理论包含纤维和基体两个不同的材料长度尺度参数。不同于忽略横向剪切应力的修正偶应力Kirchhoff薄板理论,Mindlin层合板考虑横向剪切变形引入两个转角变量。进一步建立了只含一个材料细观参数的偶应力Mindlin层合板工程理论的稳定性模型。计算了正交铺设简支方板Mindlin层合板的临界载荷。计算结果表明该模型可以用于分析细观尺度层合板稳定性的尺寸效应。     

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

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

不同边界条件正交各向异性中厚板的振动与稳定分析

本文采用高阶剪切变形理论对正交各向异性中厚矩形板进行振动与稳定分析,数值计算采用样条有限点法,得出了六种不同边界条件矩形板的自振频率和屈曲载荷,并与相应的经典板理论的结果进行比较.结果说明横向剪切变形对复合材料层合板的影响与板的各向异性程度、板的宽厚比(b/h)、层合板的层数和板的支承条件有关,它随着层合板各向异性程度的增加而增加,随着层合板宽厚比的增加而逐渐消失.     

变角度纤维复合材料层合斜板的颤振分析

假设纤维方向角沿层合板的长度方向线性变化，研究了变角度纤维复合材料层合斜板的颤振．通过坐标变换将斜板变换为正方形板，采用层合板表面连续变化的速度环量来模拟空气对其的作用，速度环量分布利用Cauchy积分公式计算．建立了系统的Lagrange方程并采用Ritz法得到了层合板的自振频率和颤振/不稳定性分离临界速度．通过数值算例验证了本文模型和方法的正确性和收敛性，分析了各个铺层内纤维方向角的变化对自振频率和颤振/不稳定性分离临界速度的影响．研究结果表明，通过纤维的变角度铺设，可有效地提高层合板的基频和颤振/不稳定性分离临界速度．经合理设计的变角度复合材料层合板具有抑制颤振的作用．     

基于连续丝束剪切技术的变角度复合材料层合板的热屈曲分析

连续丝束剪切(Continuous Tow Shearing, CTS)铺放技术是一种新的变角度层合板制作技术，这种新技术能显著减少丝束铺放过程中产生的丝束重叠和间隙等缺陷，然而，采用CTS技术铺设时，层合板的厚度将随着纤维角度的变化而变化．本文基于一阶剪切变形理论并应用Chebyshev-Ritz法对这种变厚度的变角度复合材料层合板的热屈曲问题进行了研究．假设纤维方向角沿板的长度方向按照线性变化，获得了均匀温度变化时变厚度层合板的临界热屈曲荷载．通过与现有文献的比较验证了本文方法的正确性，并进一步讨论了纤维铺设技术、纤维方向角的变化以及边界条件的不同对变角度复合材料层合板的临界屈曲温度的影响．研究结果表明，在体积相同的情况下，采用CTS铺设较传统的自动丝束铺放(AFP)可以进一步提升变角度层合板的临界屈曲温度．本文的研究结果可为变角度复合材料的设计提供一定的参考．     

层合板六参量几何非线性高阶剪切理论

提出了层合板六参量的高阶剪切变形理论的位移场假定,以考虑在大变形条件下层合板法向变形和厚度的变化。同时对ｖｏｎ Ｋａｒｍａｎ应变位移简化假设进行补充修改,考虑某些有限变形条件下被忽略小量的影响,建立了对应于该文六参量模型和更加适凳大变形分析的层合板几何非线性关系,平衡方程和边界条件。利用该文模型分析了橡胶复合材料简支板的大变形弯曲行为,并对比Ｒｅｄｄｙ五参量几何非线性简单高阶剪切变形层合理论解和弹     

基于整体局部1,3位移模式的层合板层间应力研究

为提高层合板层间应力计算的准确性,对Reddy型高阶剪切理论的基本位移模式进行改进,提出整体-局部1,3高阶位移模式．在满足层间位移连续,层间剪切应力连续,以及上下表面自由的条件下,与前人提出的整体-局部1,2-3位移模式相比,层合板板结构每个节点的独立变量由13缩减到11,并且不随层数的增加而变化．将整体-局部1,3高阶位移模式位移和应力的数值解与解析解进行对比,验证了整体-局部1,3位移模式的准确性,可应用于复合材料层合板的位移和应力分析．     

Higher order zig-zag theory for fully coupled thermo-electric–mechanical smart composite plates

A higher order zig-zag plate theory is developed to refine the prediction of the mechanical, thermal, and electric behaviors fully coupled. Both in-plane displacement and temperature fields through the thickness are constructed by superimposing linear zig-zag field to the smooth globally cubic varying field. Smooth parabolic distribution through the thickness is assumed in the out-of-plane displacement field in order to consider transverse normal deformation. Linear zig-zag form is adopted in the electric potential. The layer-dependent degrees of freedom of displacement and temperature fields are expressed in terms of reference primary degrees of freedom by applying interface continuity conditions as well as bounding surface conditions of transverse shear stresses and transverse heat fluxes. Thus the proposed theory is not only accurate but also efficient. Through the numerical examples of coupled and uncoupled analysis, the accuracy and efficiency of the present theory are demonstrated. The present theory is suitable in the predictions of fully coupled behaviors of thick smart composite plate under mechanical, thermal, and electric loads combined.     

Corotational nonlinear analyses of laminated shell structures using a 4-node quadrilateral flat shell element with drilling stiffness

A new 4-node quadrilateral flat shell element is developed for geometrically nonlinear analyses of thin and moderately thick laminated shell structures. The fiat shell element is constructed by combining a quadrilateral area co- ordinate method （QAC） based membrane element AGQ6- II, and a Timoshenko beam function （TBF） method based shear deformable plate bending element ARS-Q12. In order to model folded plates and connect with beam elements, the drilling stiffness is added to the element stiffness matrix based on the mixed variational principle. The transverse shear rigidity matrix, based on the first-order shear deformation theory （FSDT）, for the laminated composite plate is evaluated using the transverse equilibrium conditions, while the shear correction factors are not needed. The conventional TBF methods are also modified to efficiently calculate the element stiffness for laminate. The new shell element is extended to large deflection and post-buckling analyses of isotropic and laminated composite shells based on the element independent corotational formulation. Numerical re- sults show that the present shell element has an excellent numerical performance for the test examples, and is applicable to stiffened plates.     

A new trigonometric shear deformation theory for isotropic,laminated composite and sandwich plates

A new trigonometric shear deformation theory for isotropic and composite laminated and sandwich plates, is developed. The new displacement field depends on a parameter “m”, whose value is determined so as to give results closest to the 3D elasticity bending solutions. The theory accounts for adequate distribution of the transverse shear strains through the plate thickness and tangential stress-free boundary conditions on the plate boundary surface, thus a shear correction factor is not required. Plate governing equations and boundary conditions are derived by employing the principle of virtual work. The Navier-type exact solutions for static bending analysis are presented for sinusoidally and uniformly distributed loads. The accuracy of the present theory is ascertained by comparing it with various available results in the literature. The results show that the present model performs as good as the Reddy’s and Touratier’s shear deformation theories for analyzing the static behavior of isotropic and composite laminated and sandwich plates.     

Higher order zig-zag theory for smart composite shells under mechanical-thermo-electric loading

A higher order zig-zag shell theory based on general tensor formulation is developed to refine the predictions of the mechanical, thermal, and electric behaviors. All the complicated curvatures of surface including twisting curvatures can be described in a geometrically exact manner in the present shell theory because the present theory is based on the geometrically exact surface representation. The in-surface displacement fields are constructed by superimposing the linear zig-zag field to the smooth globally cubic varying field through the thickness. Smooth parabolic distribution through the thickness is assumed in the out-of-plane displacement in order to consider transverse normal deformation and stress. The layer-dependent degrees of freedom of displacement fields are expressed in terms of reference primary degrees of freedom by applying interface continuity conditions as well as bounding surface free conditions of transverse shear stresses. Thus the proposed theory has only seven primary displacement unknowns and they do not depend upon the number of layers. To assess the validity of present theory, the developed theory is evaluated under the thermal and electric load as well as under the mechanical load of composite cylindrical shells. Through the numerical examples, it is demonstrated that the proposed smart shell theory is efficient because it has the minimal degrees of freedom. The present theory is suitable in the predictions of deformation and stresses of thick smart composite shells under the mechanical, thermal, and electric loads combined.     

A hybrid stress approach for laminated composite plates within the First-order Shear Deformation Theory

This paper presents a hybrid stress approach for the analysis of laminated composite plates. The plate mechanical model is based on the so called First-order Shear Deformation Theory, rationally deduced from the parent three-dimensional theory. Within this framework, a new quadrilateral four-node finite element is developed from a hybrid stress formulation involving, as primary variables, compatible displacements and elementwise equilibrated stress resultants. The element is designed to be simple, stable and locking-free. The displacement interpolation is enhanced by linking the transverse displacement to the nodal rotations and a suitable approximation for stress resultants is selected, ruled by the minimum number of parameters. The transverse stresses through the laminate thickness are reconstructed a posteriori by simply using three-dimensional equilibrium. To improve the results, the stress resultants entering the reconstruction process are first recovered using a superconvergent patch-based procedure called Recovery by Compatibility in Patches, that is properly extended here for laminated plates. This preliminary recovery is very efficient from the computational point of view and generally useful either to accurately evaluate the stress resultants or to estimate the discretization error. Indeed, in the present context, it plays also a key role in effectively predicting the shear stress profiles, since it guarantees the global convergence of the whole reconstruction strategy, that does not need any correction to accommodate equilibrium defects. Actually, this strategy can be adopted together with any plate finite element. Numerical testing demonstrates the excellent performance of both the finite element and the reconstruction strategy.     

Radial basis functions collocation and a unified formulation for bending,vibration and buckling analysis of laminated plates,according to a variation of Murakami’s zig-zag theory

In this paper, we propose a variation of the use of Murakami’s zig-zag theory for the analysis of laminated plates. The new theory accounts for through-the-thickness deformation, by considering a quadratic evolution of the transverse displacement with the thickness coordinate. The equations of motion and the boundary conditions are obtained by the Carrera’s Unified Formulation, and further interpolated by collocation with radial basis functions. This paper considers the analysis of static deformations, free vibrations and buckling loads on laminated composite plates.     

Generalized layerwise mechanics for the static and modal response of delaminated composite beams with active piezoelectric sensors

A coupled linear layerwise laminate theory and a beam FE are formulated for analyzing delaminated composite beams with piezoactuators and sensors. The model assumes zig-zag fields for the axial displacements and the electric potential and it treats the discontinuities in the displacement fields due to the delaminations as additional degrees of freedom. The formulation naturally includes the excitation of piezoelectric actuators, their interactions with the composite laminate, and the effect of delamination on the predicted sensory voltage. The quasistatic and modal response of laminated composite Gr/Epoxy beams with active or sensory layers having various delamination sizes is predicted. The numerical results illustrate the strong effect of delamination on the sensor voltage, on through the thickness displacement and on the stress fields. Finally, the effect of delamination on modal frequencies and shapes are predicted and compared with previously obtained experimental results.     

Buckling and post-buckling of laminated composite plates

This paper addresses the buckling and post-buckling of laminated composite plates using higher order shear deformation theory associated with Green–Lagrange non-linear strain–displacement relationships. All higher order terms arising from nonlinear strain–displacement relations are included in the formulation. The present plate theory satisfies zero transverse shear strain conditions at the top and bottom surfaces of the plate in von Karman sense. A C⁰ isoparametric finite element is developed for the present nonlinear model.     

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.     

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

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