复合材料层合板自由振动分析的无网格自然邻接点Petrov-Galerkin法

基于一阶剪切变形理论,提出了复合材料层合板自由振动分析的无网格自然邻接点Petrov-Galerkin法。计算时在复合材料层合板中面上仅需要布置一系列的离散节点,并利用这些节点构建插值函数。在板中面上的局部多边形子域上,采用加权余量法建立复合材料层合板自由振动分析的离散化控制方程,并且这些子域可由Delaunay三角形方便创建。自然邻接点插值形函数具有Kronecker delta函数性质,因而无需经过特别处理就能准确地施加本质边界条件。对不同边界条件、不同跨厚比、不同材料参数和不同铺设角度的复合材料层合板,由本文提出的无网格自然邻接点Petrov-Galerkin法进行自由振动分析时均可得到满意的结果。数值算例结果表明,本文方法求解复合材料层合板的自由振动问题是行之有效的。     

基于逆复合二次径向基函数的对称复合材料层合板自由振动分析

利用逆复合二次径向基函数无网格配点法对Reddy的高阶剪切变形理论进行离散,预测了对称复合材料层合板的自由振动特性.将不同材料参数、几何尺寸和边界条件的层合板固有频率计算结果与相关文献中的结果进行对比,结果表明:逆复合二次径向基函数在对称复合材料层合板自由振动分析方面具有收敛性好及精度高等一系列优点.     

纤维曲线铺放的复合材料层合板的自由振动分析

假设曲线纤维的方向角沿板的长度方向按照线性规律变化,导出了纤维曲线铺设时的参考路径,将参考路径沿板的宽度方向平移可得一种曲线纤维增强复合材料单层板,当这种纤维曲线铺设的单层板对称铺放时即可得相应的曲线纤维增强复合材料层合板．基于弹性薄板的小挠度理论,建立了曲线纤维增强复合材料层合板自由振动问题的基本方程,采用微分求积法进行数值求解,得到了层合板的自振频率及相应的振型．与已有文献计算结果的比较,验证了本文计算结果的正确性．通过数值算例分析了微分求积法求解本问题时的收敛性,研究了纤维铺放路径和边界条件的不同对曲线纤维增强复合材料层合板频率及振型的影响,研究结果可为该种结构的设计提供一定的参考．     

反对称铺设复合材料层合板非线性后屈曲分析

基于层合板壳理论,考虑反对称铺设层合板的拉弯耦合效应和后屈曲过程中的非线性几何变形,推导了由应力函数和挠度表示的复合材料层合板的后屈曲控制方程。引入无量纲参数对控制方程和边界条件进行无量纲化,以消除材料参数及几何尺寸对分析结果的影响。采用摄动法将无量纲的非线性控制方程及边界条件展开成一系列非齐次线性摄动方程组,分析各阶摄动方程的通解与特解的构造,并逐次求解,建立了反对称铺设复合材料层合板受单向均布压力作用的临界屈曲荷载及后屈曲平衡路径的理论解。进而运用ABAQUS软件对复合材料层合板在面内压缩载荷作用下的屈曲和后屈曲进行有限元分析,结果表明理论解与ABAQUS结果十分接近,验证了理论解的正确性。在此基础上进一步讨论了铺设角度、铺设层数和拉弯耦合效应等对层合板后屈曲性能的影响。研究发现层合板的屈曲载荷受铺设角度与层数的影响较为显著,而拉弯耦合效应使板的屈后强度大大降低。     

基于自然单元法的功能梯度中厚板自由振动分析

自然单元法是一种基于自然邻近插值的无网格数值方法．相对于移动最小二乘近似而言,自然邻近插值不涉及到复杂的矩阵求逆运算,也不需要任何人为的参数．基于一阶剪切变形板理论,利用自然单元法对功能梯度中厚板的自由振动进行了数值分析．功能梯度板材料属性沿厚度方向呈梯度连续变化．由于自然邻近插值函数具有Kronecker delta函数性质,可以直接施加本质边界条件．通过本文给出的方法,对不同梯度指数和不同边界条件的功能梯度中厚板的振动频率进行了计算．通过与文献结果的对比验证了自然单元法求解的有效性．     

基于Kriging插值无网格法求解轴对称弹性力学问题

基于Kriging插值无网格法,提出了实际应用中复杂轴对称弹性力学问题求解的一条新途径.Kriging插值无网格法是一种新型的无网格法,该方法构造的形函数满足Kronecker delta函数性质,可以直接施加本质边界条件.采用Kriging插值无网格法分析轴对称问题,得到了轴对称问题的无网格离散方程,并编制了相应的计算程序.通过厚壁圆筒的静力学和动力学分析,对所提方法进行了检验.数值算例结果表明,提出的方法对求解轴对称弹性力学问题是行之有效的.     

变角度纤维复合材料环扇形层合板的自由振动分析

与传统的直线纤维增强复合材料相比，变角度纤维复合材料具有更强的可设计性，为改善结构性能提供了更大的可能．鉴于此，本文将研究纤维的变角度铺设对复合材料环扇形层合板的自振频率及振动模态的影响．假设纤维的方向角沿环扇形板的径向线性变化，基于经典的层合板理论，采用微分求积法获得了环扇形层合板自由振动问题的数值解．通过与现有文献及ABAQUS有限元结果的比较验证了本文模型及方法的正确性和收敛性，并详细分析了纤维起始角和终止角的变化对层合板的自振频率及振动模态的影响．研究结果表明：与常刚度层合板相比，变角度纤维复合材料层合板的基频具有更大的调整空间，通过合理选择纤维起始角和终止角可有效提高层合板的基频．研究结果可为该种新型复合材料结构的优化设计提供一定的参考．     

具损伤压电智能层合板的动力分析

基于应变能等效原理、高阶剪切变形理论和Hamilton变分原理,考虑复合材料铺设层内的损伤效应,建立了具损伤压电智能层合板的运动控制方程,并运用Galerkin方法进行求解.数值算例中,讨论了损伤效应、厚跨比及压电层厚度与层合板总厚度之比对四边简支压电智能层合板自由振动频率的影响和外部控制电压对其动力响应的影响.     

基于各向异性修正偶应力理论的Reddy型层合板的自由振动

本文基于各向异性修正偶应力理论建立了只含一个尺度参数的Reddy型复合材料层合板的自由振动模型。同见诸于文献的细观尺度Kirchhoff薄板偶应力模型相比,本文提出的新模型能够更精确的预测细观尺度下的中、厚层合板的自振频率。基于Hamilton原理推导了细观尺度下Reddy型复合材料层合板的运动微分方程以及边界条件,并以正交铺设的四边简支复合材料层合方板为例进行了解析求解,分析了尺度参数对自振频率的影响并对比了Kirchhoff、Mindlin和Reddy等三种板模型计算结果的异同。算例结果表明本文所给出的模型能够捕捉到复合材料层合板自由振动问题的尺度效应。另外,在细观尺度下Kirchhoff板模型所预测的自振频率相对于Mindlin板模型和Reddy板模型总是过高,且越接近厚板三者的差别就越大,这与经典理论中三种板模型的对比情况是一致的。     

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

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

复合材料正交加筋板在横向载荷下挠度计算的半解析法

本文采用一种新的半解析法,即独特利用Heaviside函数建立与加筋板等效的变刚度模型来开展复合材料双向正交加筋板在横向载荷下的弯曲挠度分析．此模型可以准确地描述筋条在板面上的分布,以及由于筋条的存在而导致的板面刚度不均匀分布．使用Galerkin加权残值法求解该模型的控制方程,得到不同边界条件和载荷情况下的级数解．对于双向正交加筋板,将此半解析法的结果与传统均匀化方法和使用商业有限元软件ABAQUS建立的有限元模型所得到的弯曲挠度结果比较,验证了此方法的准确性和优越性．不同于传统均匀化方法,本双向正交加筋板的弯曲挠度半解析法可精确、有效地获取加筋间的局部弯曲挠度,可以促进复合材料结构的设计分析与优化的研究进展．     

Three-dimensional free vibration analysis of rotating laminated conical shells: layerwise differential quadrature (LW-DQ) method

This paper focuses on the free vibration analysis of thick, rotating laminated composite conical shells with different boundary conditions based on the three-dimensional theory, using the layerwise differential quadrature method (LW-DQM). The equations of motion are derived applying the Hamilton’s principle. In order to accurately account for the thickness effects, the layerwise theory is used to discretize the equations of motion and the related boundary conditions through the thickness of the shells. Then, the equations of motion as well as the boundary condition equations are transformed into a set of algebraic equation applying the DQM in the meridional direction. This study demonstrates the applicability, accuracy, stability and the fast rate of convergence of the present method, for free vibration analyses of rotating thick laminated conical shells. The presented results are compared with those of other shell theories obtained using conventional methods and a special case where the angle of the conical shell approaches zero, that is, a cylindrical shell and excellent agreements are achieved.     

Dynamic analysis of composite plate with multiple delaminations based on higher-order zigzag theory

In a recent paper, Cho and Kim [Journal of Applied Mechanics] proposed a higher-order cubic zigzag theory of laminated composites with multiple delaminations. The proposed theory is not only accurate but also efficient because it work with a minimal number of degrees of freedom with the application of interface continuity conditions as well as bounding surface conditions of transverse shear stresses including delaminated interfaces. In this work, we investigate the dynamic behavior of laminated composite plates with multiple delaminations. A four-node finite element based on the efficient higher-order zigzag plate theory of laminated composite plates with multiple delaminations is developed to refine the prediction of frequencies, mode shape, and time response. Through the dynamic version of the variational approach, the dynamic equilibrium equations and variationally consistent boundary conditions are obtained. Natural frequency prediction and time response analysis of a composite plate with multiple delaminations demonstrate the accuracy and efficiency of the present finite element method. To prevent penetration violation at the delamination interfaces, unilateral contact constraints by Lagrange multiplier method are applied in the time response analysis. The present finite element is suitable for the prediction of dynamic response of thick composite plates with multiple and arbitrary shaped delaminations.     

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.     

NONLINEAR THREE-DIMENSIONAL ANALYSIS OF COMPOSITELAMINATED PLATES

ＮＯＮＬＩＮＥＡＲＴＨＲＥＥ－ＤＩＭＥＮＳＩＯＮＡＬＡＮＡＬＹＳＩＳＯＦＣＯＭＰＯＳＩＴＥＬＡＭＩＮＡＴＥＤＰＬＡＴＥＳ￥(江晓禹，张相周，陈百屏)ＪｉａｎｇＸｉａｏｙｕ；（ＳｏｕｔｈｗｅｓｔｅｒｎＪｉａｏｔｏｎｇＵｎｉｖｅｒｓｉｔｙ，Ｃｈｅｎｇｄｕ６...     

FREE VIBRATION OF ANISOTROPIC RECTANGULAR PLATES BY GENERAL ANALYTICAL METHOD

According to the differential equation for transverse displacement function of anisotropic rectangular thin plates in free vibration, a general analytical solution is established. This general solution, composed of the composite solutions of trigonometric function and hyperbolic function, can satisfy the problem of arbitrary boundary conditions along four edges. The algebraic polynomial with double sine series solutions can also satisfy the problem of boundary conditions at four corners. Consequently, this general solution can be used to solve the vibration problem of anisotropic rectangular plates with arbitrary boundaries accurately. The integral constants can be determined by boundary conditions of four edges and four corners. Each natural frequency and vibration mode can be solved by the determinate of coefficient matrix from the homogeneous linear algebraic equations equal to zero. For example, a composite symmetric angle ply laminated plate with four edges clamped has been calculated and discussed.     

Sensitivity analysis of composite laminated plates with bonding imperfection in Hamilton system

Sensitivity analysis of composite laminated plates with bonding imperfection is carried out based on the radial point interpolation method （RPIM） in a Hamilton system. A set of hybrid governing equations of response and sensitivity quantities is reduced using the spring-layer model and the modified Hellinger-Reissner （H-R） variational principle. The analytical method （AM）, the semi-analytical method （SAM）, and the finite difference method （FDM） are used for sensitivity analysis based on the reduced set of hybrid governing equations. A major advantage of the hybrid governing equations is that the convolution algorithm is avoided in sensitivity analysis. In addition, sensitivity analysis using this set of hybrid governing equations can obtain response values and sensitivity coefficients simultaneously, and accounts for bonding imperfection of composite laminated plates.     

高阶剪切变形理论下三边夹紧一边铰支复合材料层板的几何非线性分析

首先用虚位移原理推导出以位移形式表达的Reddy型高阶剪变形理论复合材料层板的非线性控制方程及相应的边界条件。选定的五个位移函数均满足三边夹紧一边铰支边界条件，用Galerkin方法把无量纲化之后的控制方程转化为一组非线性代数方程组，用线性化的方法和可调节参数的修正迭代法求解这组方程。最后求出了不同复合材料的挠度和弯矩值。