考虑横法向热应变的C~0型Reddy板理论和三角形板单元

考虑横法向热变形,建议了C0型Reddy理论,并用于分析复合材料层合/夹层板热膨胀问题。虽然考虑了横法向热应变,但不增加额外的位移变量。此理论位移场不含有横向位移一阶导数,构造有限元时仅需C0插值函数。基于这一模型,运用虚位移原理推导了复合材料板平衡方程以及构造了6节点三角形板单元,并分析了简支复合材料层合/夹层板的热膨胀问题。数值结果表明,建立的模型能准确分析复合材料层合/夹层板热膨胀问题,而忽略横法向热应变的理论分析热膨胀问题误差较大。     

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

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

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

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

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

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

一种预测层合板层间剪应力的新后处理方法

层合板是航空航天领域典型的承力构件,过大的层间应力是导致其分层失效的主要原因.准确的层间应力预测往往依赖于三维平衡方程后处理方法(TPM).然而,该方法需要计算面内应力的一阶导,使得基于C0型板理论构造的线性单元无法使用TPM计算横向剪应力.本文在三维平衡方程后处理方法的基础上,提出了一种新后处理方法(NPM).新后处理方法通过虚功等效法消除了三维平衡方程后处理方法中产生的位移参数的高阶导.基于提出的新后处理方法和C0型板理论,仅需使用线性单元就可以预测层合板的横向剪应力.为了验证所提方法的有效性,本文基于修正锯齿理论(RZT)和所提方法构造了一种C0连续的三节点三角形线性板单元.数值算例表明,所提方法和三维平衡方程后处理方法具有相同的计算精度,提出的板单元能够准确高效地预测层合板的横向剪应力.此外,所提方法便于结合现有的有限元商用软件使用,基于商用软件中板壳单元获得的节点位移,使用新后处理方法极易获得准确的层间剪应力.     

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

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

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

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

复合材料层合板层间应力分析的混合有限元法

本文根据整体-局部变分原理,构造了一种适用于复合材料层合板层间应力分析的混合有限单元,本文以层合厚板和典型Pagano问题为例证明了该单元的有效性和可靠性。     

薄壁弯曲结构有限元计算精度研究

基于h-p型有限元精度计算法,以薄壁弯曲结构为研究对象,系统地介绍了实体单元常见的分类方法及优缺点;通过理论公式推导了薄壁弯曲结构发生弹性和弹塑性变形时的位移和应力理论解;采用有限元法计算数值解,研究了影响有限元计算精度的因素和规律,并用算例证实了研究结果的合理性。研究结果表明:当单元类型、积分方式、阶次、长高比相同时,只有1层实体单元情况下得到的计算误差总是大于多层单元;只要严格控制单元长高比为1左右,单元层数不小于4层,采用一阶全积分六面体单元就可以控制位移及应力误差在5%以内;当采用一阶减缩积分六面体单元,只需2层单元就可以控制弹性位移误差在1%左右,但此时应力误差达30%以上,对于塑性变形,单元层数达6层时其位移误差仍达8%以上;对于二阶六面体及二阶四面体单元,只需2层单元,且不需严格控制单元长高比为1左右就可以使位移及应力计算误差在5%以内。     

新修正偶应力理论Reddy型层合板稳定分析

基于新修正偶应力理论建立了一个Reddy型复合材料层合板稳定性模型。该理论中曲率张量不对称,而偶应力矩张量对称。Reddy型层合板模型能够满足横向剪切应力为0的自由表面条件,而且横向剪切为二次函数,避免了常剪力一阶理论需要引入的剪力修正系数。为了便于工程应用,通过虚功原理推导了只含纤维材料尺度参数正交铺设的Reddy型层合板偶应力模型的稳定性方程,并以微尺度正交铺设四边简支层合方板为例,分析了不同铺设角和轴向载荷作用时临界载荷的细观尺度效应,并且与一阶剪切变形和Kirchhoff板理论结果对比。结果表明,本文建立的新修正偶应力Reddy型层合板模型更适合分析较厚的复合材料层合板稳定性的尺度效应。     

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.     

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.     

板结构计算模型的新发展

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

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.     

Bending analysis of functionally graded CNT reinforced doubly curved singly ruled truncated rhombic cone

The bending analysis of functionally graded carbon nanotube (CNT) reinforced doubly curved singly ruled truncated rhombic cone is investigated. In this study, a simple C⁰ isoparametric finite element formulation based on third order shear deformation theory is presented. To characterize the membrane-flexure behavior observed in a CNT reinforced truncated rhombic cone, a displacement field involving higher-order terms in in-plane fields is considered. The proposed kinematics field incorporates for transverse shear deformation and nonlinear variation of the in-plane displacement field through the thickness to predict the overall response of the CNT reinforced truncated rhombic cone in an accurate sense. The material properties of the CNT reinforced truncated rhombic cone are estimated according to the rule of mixture. The present model eliminates the need of shear correction factor and imposed zero-transverse shear strain at upper and lower surface of the truncated rhombic cone. The new feature in present model is simultaneous inclusion of twist curvature in strain field as well as curvature in displacement field that makes it suitable for moderately thick and deep truncated rhombic cone. The proposed new mathematical model is implemented in finite element code written in FORTRAN. The proposed model has been validated with analytical, experimental, and finite element results from the literature. This is first attempt to study bending response of CNT reinforced doubly curved singly ruled truncated rhombic cone. The effect of CNT distribution, boundary condition, loading pattern, and other geometric parameters are also examined.     

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.