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

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

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

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

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

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

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

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

复合材料脱层梁的屈曲分析

本文研究了具有任意位置透型脱层的复合材料梁的屈曲问题。基于弹性理论建立了复合材料脱层梁的基本方程式。对脱层梁进行了分区处理,利用B样条函数作为位移型函数的基函数,方便地描述了脱层长度、脱层位置。考虑边界条件、区间位移连续性条件和弯矩剪力的平衡条件以及纵向内力的附加条件,对基本方程式进行了求解。得出了脱层位置不同,脱层长度不同的屈曲荷载的变化规律,并与轴对称脱层时的屈曲荷载进行了比较,认为层合梁考虑脱层对屈曲的影响是非常必要的。     

压电层合梁静动力学分析的高精度梁单元

给出了一个对复合材料压电层合梁进行数值分析的高精度压电层合梁单元。基于Shi三阶剪切变形板理论的位移场和Layer-wise理论的电势场,用力-电耦合的变分原理及Hamilton原理推导了压电层合梁单元列式。采用拟协调元方法推导了一个可显式给出单元刚度矩阵的两节点压电层合梁单元,并应用于压电层合梁的力-电耦合弯曲和自由振动分析。计算结果表明,该梁单元给出的梁挠度和固有频率与解析解吻合良好,并优于其它梁单元的计算结果,说明了本文所给压电层合梁单元的可靠性和准确性。研究结果可为力-电耦合作用下压电层合梁的力学分析提供一个简单、精确且高效的压电层合梁单元。     

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

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

形状记忆合金纤维复合材料薄壁梁的主动变形模型

提出具有变形主动驱动作用的SMA纤维混杂复合材料单闭室薄壁截面梁的力-位移本构关系模型。基于变分渐进法导出具有SMA主动纤维的复合材料薄壁空心梁的二维截面刚度系数以及截面内力（矩）与位移（转角）关系方程,含SMA纤维层合板材料性能由混合率进行预测。基于Tanaka的SMA应力应变关系以及Lin的线性相变动力模型,导出了SMA诱发的轴力、扭矩与弯矩的数学表达式。由本文建立的具有拉伸-扭转-弯曲静变形耦合的一般公式出发,讨论周向均匀刚度配置以及周向反对称刚度配置特殊情形,并给出了简化的本构方程。在不考虑SMA纤维含量和温度变化的情况下,本文的模型可以退化为普通纤维复合材料单闭室薄壁截面梁的已有结果。通过数值计算揭示了SMA对弯曲-扭转静变形特性的作用规律,分析了SMA纤维含量与初始应变、驱动温度和复合材料铺层角的影响。     

一种复合材料旋翼桨叶结构大变形分析方法的改进及验证

提出了一种适用于直升机旋翼复合材料桨叶大变形分析的改进方法。将旋翼桨叶变形分析分解为一维非线性分析和二维剖面特性分析,并考虑横向剪切、翘曲对剖面刚度及弹性耦合的影响;为使方法适用于旋翼气动弹性分析,将应变能中的广义应变用参考轴线处的弹性运动表示,保留所有非线性项,推导出计算复合材料桨叶大变形的公式;采用有限元法处理方程,对梁结构进行了分析,并将大变形状态下的位移计算结果与Princeton梁实验值、Minguet复合材料梁实验值以及中等变形梁理论计算结果进行了比较,验证了大变形状态下本文计算方法的正确性;此外与中等变形梁模型计算结果的对比,验证了本文方法在计算精度上的提高。     

形状记忆合金混杂复合材料薄壁梁主动变形模型

提出具有变形主动驱动作用的SMA纤维混杂复合材料单闭室薄壁截面梁的力-位移本构关系模型.基于变分渐近法导出具有SMA主动纤维的复合材料薄壁空心梁的二维截面刚度系数以及截面内力(矩)与位移(转角)关系方程,含SMA纤维层合板材料性能由混合率进行预测.基于Tanaka的SMA应力应变关系以及Lin的线性相变动力模型,导出了SMA诱发的轴力、扭矩与弯矩的数学表达式.由该文建立的具有拉伸-扭转-弯曲静变形耦合的一般公式出发,讨论周向均匀刚度配置以及周向反对称刚度配置特殊情形,并给出了简化的本构方程.在不考虑SMA纤维含量和温度变化的情况下,本文的模型可以退化为普通纤维复合材料单闭室薄壁截面梁的已有结果.通过数值计算揭示了SMA对弯曲-扭转静变形特性的作用规律,分析了SMA纤维含量、驱动温度和复合材料铺层角的影响.     

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.     

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

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

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.     

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.     

Measurement of thermal expansion coefficients of composites using strain gages

The measurement of the coefficients of thermal expansion (CTEs) of composite materials using electrical resistance strain gages is addressed. Analytical expressions for the CTEs of an orthotropic lamina are derived, accounting for the effects of transverse sensitivity and possible misalignment of the gages. Experiments are performed for the characterization of the thermal expansion behavior of a fiber-glass-reinforced epoxy unidirectional lamina using an invar specimen as reference material. Preliminary training cycles are performed for the determination of an optimal heating rate for the measurements, which ensures thermal equilibrium conditions. Three measurement cycles yield the principal CTEs of the lamina α₁, α₂ and α₁₂ with repeatability within ±0.34×10⁻⁶, ±0.85×10⁻⁶ and ±2.8×10⁻⁶/°C, respectively. It is noted that inhomogeneity of the specimen and variation in thermomechanical properties of the gages can cause a noticeable spead in the measurements.     

Thermal bending analysis of laminated orthotropic plates by the generalized differential quadrature method

The interlaminar stresses in a thin laminated rectangular orthotropic plate with four sides simply supported edges under bending was determined by using the generalized differential quadrature (GDQ) method involving the effects of thermal expansion strain and transverse load. The approximate stress and displacement solutions are obtained under the effects of thermal expansion force and uniform pressure load for eight-layer unidirectional laminates, symmetric cross-ply laminates. Numerical results on the dominant interlaminar stresses and displacement of bending analysis are compared to the Navier solution. The thermal induced forces have significant effect on the bending of plates.     

无网格伽辽金法求解平面偶应力问题

提出采用无网格伽辽金法(EFGM)求解偶应力问题,以避免有限元求解中因C¹连续要求可能引起的不便。推导了基于二次基和移动最小二乘技术的EFGM计算公式,通过计算受轴向均匀拉伸的带中心小孔无限大板和细长杆的偶应力问题,对所提方法进行了数值验证,与解析解相比结果令人满意,此外还讨论了节点密度和权函数对计算结果的影响。数值结果表明,所提算法可有效地求解平面偶应力问题。