考虑横法向热应变的C~0型整体-局部高阶层合梁理论

通过考虑横法向热变形,本文建立了预先满足层间应力连续的C0型整体-局部高阶层合梁理论,并用于分析复合材料层合梁热膨胀和热弯曲问题。虽然考虑了横法向应变,不增加额外的位移变量。此理论位移场不含有横向位移一阶导数,便于构造多节点高阶单元。基于虚功原理推导了复合材料层合梁平衡方程,并分析了简支多层复合材料梁热膨胀和热弯曲问题。数值结果表明,建立的模型能准确分析复合材料层合梁热膨胀和热弯曲问题,忽略横法向应变的理论分析热膨胀问题误差较大。     

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

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

C~0型高阶锯齿理论及复合材料厚梁热应力分析

基于已有锯齿理论构造单元时,需使用满足单元间C1连续的插值函数,难于构造多节点高阶单元,而且精度较低。针对已有锯齿理论存在的问题,本文首先发展了C0型锯齿理论。通过虚位移原理推导出在热载荷作用下复合材料梁的平衡方程,并给出了简支复合材料层合梁解析解。基于发展的锯齿理论分析了复合材料夹层梁和层合梁热膨胀问题,并与其他理论结果对比。数值结果表明,发展的C0型锯齿理论能克服已有锯齿理论的难题。     

复合材料厚梁精化锯齿理论及有限元分析

由于具有预先满足层间应力连续的优点,锯齿理论被广泛研究和应用。然而,至今锯齿理论仍然存在如下难题：基于锯齿理论构造单元时,需使用满足单元间C1连续的插值函数,难于构造多节点高阶单元,而且精度较低。如果这些问题不被重视和解决,应用此类理论分析复合材料力学问题可能得出不恰当的结论。通过发展高精度的考虑横法向应变的C0型锯齿理论,本文将克服已有锯齿理论遇到的上述难题。基于发展的锯齿理论,构造三节点梁单元验证发展理论模型的性能。     

正交复合材料面层夹层板非线性理论及应用

针对复合材料面层夹层板的构造和变形特点，考虑横向剪应力在面层和芯层粘结处的连续条件，应用Ｈａｍｉｌｔｏｎ原理建立了基于五个未知函数的正交铺设复合材料面层夹层板的非线性精化理论。对静力学问题，控制方程化简为由四个基本未知函数表述。文中还分析了简支正交铺设复合材料面层夹层板的非线性弯曲，给出了载荷—挠度特征关系和板中应力的分布状况。数值计算表明，夹层板面层和芯层粘结处的层间剪应力在工程设计中是十分重要的。     

用于多种夹层板等效的有限元分析法

针对夹层板力学性能解析法难于计算复杂结构的夹层板且通用性差的问题,本文采用有限元分析法研究了夹层板性能的等效方法。对夹层板的代表体单元模型施加位移约束,模拟弯曲变形时线性独立的应变分量和弯曲内力;根据夹层板内力与应变的本构关系,求出刚度矩阵;最后由刚度矩阵得出宏观等效弹性常数,从而把夹层板等效成连续材料的单层板单元。将该方法与解析法计算结果进行比较得到的夹层板单元四个主要弹性常数误差在0.2%以内,验证了该方法的有效性;另外采用该方法等效三种典型结构夹层板,比较实际模型和等效模型的弯曲响应,得到的误差均在1.4%以内,表明该方法在不考虑复杂多变的夹芯结构时具有通用性。     

基于Layerwise离散层理论的复合材料夹层板屈曲分析

复合材料夹层结构由于面板和芯层力学特性差异较大,屈曲分析时要分层考虑各层的剪切变形。基于Reddy的Layerwise离散层理论,假设每一层变形服从一阶剪切变形理论,在统一的位移场描述下,推导建立了一种用于复合材料夹层结构屈曲分析的四节点四边形板单元,并采用混合插值方法对单元的剪切锁定进行了修正。分别对三种典型的夹层板结构进行线性屈曲有限元分析,并将计算结果与文献中已有结果进行了对比。结果表明：本文的分析方法能离散考虑各层的力学特性,将结构离散为多层时,计算结果与三维弹性理论或高阶板理论吻合;将结构等效为单层时,计算结果与基于一阶剪切变形理论的文献结构吻合,验证了单元的有效性。     

考虑界面损伤层合梁板的层间应力分析

基于精确应力分析的广义六自由度板理论,应用变分原理和损伤力学中 的应变等效原理,考虑复合材料铺设层内和层间界面处的损伤效应,建立了具两种损伤模式 的复合材料层合板的三维非线性平衡微分方程,且运用有限差分法对考虑损伤简支层合梁板 的层间应力进行了求解.     

考虑横向拉伸影响的FGM夹层板动不稳定性分析

以受多种形式面内载荷作用下的功能梯度夹层矩形板为研究对象,夹层板的功能梯度层组分材料沿厚度方向按幂律分布,考虑组分材料物性参数的温度相关性.研究中所用的位移场在Reddy高阶剪切变形理论基础上考虑了横向拉伸的影响.利用能量原理和Galerkin法得到四边简支功能梯度材料夹层矩形板的动力学模型.运用Bolotin法求得了...     

基于高阶变形理论的夹层板的弯曲与振动

考虑面板和夹芯的面内刚度和横向剪切刚度以及抗弯刚度,考虑了高阶剪切变形,根据横向剪应变分布情况给出横向剪切转角的位移函数,基于哈密尔顿原理,推导了基于高阶变形理论、适用于软、硬夹芯情况夹层板的基本方程。作为算例,以四边简支条件下的夹层板的弯曲与振动,在不同的面板与芯层的弹性模量比和厚度比下进行了计算,并与Reissner理论、Hoff理论以及邓宗白基于Reissner理论的修正模型的计算结果进行了对比。与前述理论与方法相比,本文方法考虑因素更为全面,对夹层板的适用范围更为广泛,计算结果更为精确。针对Nastran软件计算夹层板的振动问题,对其适用范围作了简要分析。     

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 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.     

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.     

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.     

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.     

基于改进位移模式的一维C1有限元超收敛算法

提出了基于改进位移模式的一维C1有限元超收敛算法。利用单元内部需满足平衡方程的条件,推导了超收敛计算解析公式的显式,即将高阶有限元解的位移模式用常规有限元解的位移模式表示。用常规有限元解的位移模式与高阶有限元解的位移模式之和构造新的位移模式。采用积分形式推导了单元刚度矩阵。该算法在前处理阶段使用了超收敛计算公式,在常规试函数的基础上,增加了高阶试函数,使得单元内平衡方程的残差减少,从而达到提高精度的目标。对于Hermite单元,本文的结点和单元的位移、导数都达到了h4阶的超收敛精度。     

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.