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

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

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

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

含脱层损伤复合材料层合梁的冲击动力屈曲

针对含初始缺陷和脱层损伤的复合材料层合梁的轴向冲击动力屈曲问题进行了分析。基于Hamilton原理导出了考虑初始缺陷、轴向和横向惯性、横向剪切变形以及转动惯性影响时含脱层损伤复合材料梁的非线性动力屈曲控制方程;基于B-R准则,采用有限差分方法求解了受轴向冲击载荷作用下含脱层损伤复合材料梁的动力屈曲问题;讨论了冲击速度、初始几何缺陷、铺层角度以及脱层长度等因素对复合材料层合梁动力屈曲的影响。     

轴向运动SMA层合梁的主共振分析

研究了均布横向载荷作用下轴向运动SMA(形状记忆合金)层合梁的横向非线性振动。考虑轴向运动效应、轴力等因素的综合影响,利用力平衡条件、变形协调方程及SMA多项式函数的本构关系,建立了SMA层合梁在均布横向载荷作用下的动力学方程。针对两端简支边界条件,通过伽辽金积分得到了轴向运动SMA层合梁横向振动微分方程。应用平均法得到了横向载荷作用下系统主共振幅频响应方程,对理论结果进行了数值验证;分析了轴向运动速度、温度、激励参数对系统稳态响应的影响。结果表明:轴向速度、轴向载荷的变化只对系统共振频率产生影响;在外激励幅值较大时,温度增加和SMA层增厚对系统产生了相同的减振效果。     

受分布载荷复合材料层合梁应力分析的一般理论

为了克服层合梁经典理论的缺点，提高层间应力的计算精度，提出了受分布载荷层合梁应力分析的一般理论。首先根据叠加原理将原始受力状态分解成对称与反对称受力状态。然后用正交完备的三角级数和勒让德级数构造这两种受力状态中每一铺层与层间胶层的位移场，并应用广义势能原理确定位移场中的待定系数，从而确定层合梁的位移场和应力场。同时，单层梁与单层梁之间的胶层被视为各向同性材料并且与其它材料层具有相似的力学特性，即具有有限厚度、有限弹性常数。计算结果显示,这种解法的收敛性非常好，根据物理方程与根据平衡方程得到的横向剪应力和正应力分布非常一致。     

压电智能复合材料层合梁的力电耦合模型及层间应力分析

由于非凡的物理性能,石墨烯纳米片(GPL)被认为是最有吸引力的复合材料增强材料之一.GPL增强材料可以明显提高聚偏氟乙烯(PVDF)压电性能和力学性能.在力电载荷作用下,对含均匀石墨烯薄片增强(GSR)智能压电复合材料层合梁层间应力预测至关重要.若对受到力电耦合作用且层与层之间材料性能突变的压电层合梁层间剪切变形预测有误,则其层间应力过大可能导致层间失效.因此,论文提出一种适于分析此类问题且满足层与层之间相容性条件的有效力电耦合模型,用于含GSR致动器的复合材料层合梁层间应力分析.应用Reissner混合变分原理(RMVT),可以提高考虑力电耦合效应的横向剪应力预测精度.三维(3D)弹性理论和所选模型计算结果将用于评估所提梁模型性能.此外,还从力电载荷、压电层厚度、石墨烯体积分数和长厚比等方面对含GSR致动器复合材料层合梁力学响应特性进行了系统的研究.     

复合材料层合梁在横向线载荷下的脱层扩展

基于可动边界变分原理对层合梁脱层扩展进行了分析;考虑了脱层间的接触效应,建立了层合梁在横向线载荷作用下的非线性控制微分方程及相应的定解条件;应用Griffith准则导出了脱层前缘各点处的能量释放率表达式;通过算例讨论了脱层长度、脱层深度、几何尺寸、材料性质等因素对脱层扩展的影响.研究表明:脱层越长、越深、横向载荷越大,脱层越容易扩展;梁的长高比L/h及材料的E_(11)/E_(22)越大,脱层越不易发生扩展.     

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

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

基于微分求积单元法的层合梁热弹性分析

基于微分求积单元法,开展了非均匀温度场中层合梁的热弹性分析．首先基于Fourier导热定律,分析一般热边界条件下层合梁的二维稳态温度场;然后基于二维热弹性力学理论,分析层合梁的热应力和变形．为求解热传导和热应力问题,将层合梁沿各层界面划分为若干空间子域,采用微分求积法对每一子域的控制方程和边界条件进行离散并求解．数值算例验证了本方法的收敛性,与已有文献结果的对比验证了本方法的正确性．最后,算例分析了非均匀热边界条件对夹层梁温度分布的影响,以及端部支承条件和长厚比对梁内位移和应力分布的影响．     

压电复合材料层合梁的分岔、混沌动力学与控制

研究了简支压电复合材料层合梁在轴向、横向载荷共同作用下的非线性动力学、分岔和混沌动力学响应. 基于vonKarman理论和Reddy高阶剪切变形理论,推导出了压电复合层合梁的动力学方程. 利用Galerkin法离散偏微分方程,得到两个自由度非线性控制方程,并且利用多尺度法得到了平均方程. 基于平均方程,研究了压电层合梁系统的动态分岔,分析了系统各种参数对倍周期分岔的影响及变化规律. 结果表明,压电复合材料层合梁周期运动的稳定性和混沌运动对外激励的变化非常敏感,通过控制压电激励,可以控制压电复合材料层合梁的振动,保持系统的稳定性,即控制系统产生倍周期分岔解,从而阻止系统通过倍周期分岔进入混沌运动,并给出了控制分岔图.      

一个改进的平面梁单元

根据有限单元法基本原理 ,提出了一个变截面平面梁单元 ,推导了其单元钢度矩阵。这一改进的梁单元用于分析梁高呈线性变化及二次抛物线变化的矩形截面梁 ,将得到准确解。文中给出了一个变截面悬臂梁算例 ,计算表明 ,这一改进的梁单元使变截面梁的分析大大简化     

风机塔筒结构横向振动特性的快速计算方法

本文提出了一种风机塔筒结构横向振动特性的快速计算方法．将机舱和叶片整体、连接法兰盘分别简化为集中质量,塔筒简化为非均匀悬臂梁,建立风机塔筒结构横向振动方程．给出了用假设模态法计算塔筒结构固有频率和模态函数的过程．通过与文献及有限元数值结果比较验证了方法的有效性．本文方法仅需给出结构的基本参数,如截面半径变化规律、法兰盘位置和质量、机舱及叶片质量,便可快速求解其频率和模态,无需建立其复杂的力学模型．     

圆形截面约束阻尼梁整体有限元建模与振动分析

为克服约束阻尼结构在工程应用中存在的由三维有限元建模造成的单元数目巨大的问题,本文提出了圆形实心截面梁附加约束阻尼层横向振动的整体有限元建模方法,并基于此模型,研究了材料的物理属性和几何因子对这类结构的固有频率的影响。通过与三维有限元解法相比较,证明了该方法的可行性与正确性,并给出了该方法的适用范围。     

加筋板壳稳定性分析中一种简单的有限元模式

采用考虑耦合和剪切效三结点十五参数罚函数层板壳单元和与之相适应的两结点非协调层合梁单元，适用于任意和筋复合材料层合板壳的稳定性分析。数值计算表明：该单元不仅对边界形状的适应性强，而且在工程精度范围内具有简单，经济，高效的优点。     

An aeroelastic model for composite rotor blades with straight and swept tips. Part I: Aeroelastic stability in hover

An analytical model for predicting the aeroelastic behavior of composite rotor blades with straight and swept tips is presented. The blade is modeled by beam type finite elements along the elastic axis. A single finite element is used to model the swept tip. The non-linear equations of motion for the finite element model are derived using Hamilton's principle and based on a moderate deflection theory and accounts for: arbitrary cross-sectional shape, pretwist, generally anisotropic material behavior, transverse shears and out-of-plane warping. Numerical results illustrating the effects of tip sweep, anhedral and composite ply orientation on blade aeroelastic behavior are presented. It is shown that composite ply orientation has a substantial effect on blade stability. At low thrust conditions, certain ply orientations can cause instability in the lag mode. The flap-torsion coupling associated with tip sweep can also induce aeroelastic instability in the blade. This instability can be removed by appropriate ply orientation in the composite construction.     

16自由度的中厚板弯曲矩形单元

基于有限条带思想,引入结点扭率自由度,利用深梁单元的位移模式建立了一个4结点16自由度中厚板弯曲高阶单元,此单元是薄板单元BFS-16的推广形式,其特点是单元的横向位移、转角位移、剪应变位移模式直接构造,在边界上位移模式与深梁单元一致,方便与梁单元叠加,适应于带加劲肋的板弯曲问题分析,用于薄壁结构时可考虑翘曲。实例计算显示,此单元精度高,计算稳定,收敛快,无剪切闭锁现象,能较好地反映中厚板的边界效应。     

Elastodynamic analysis of low tension cables using a new curved beam element

In this paper, we address and overcome the difficulties associated with the use of the classic cable theory to treat low tension cables by developing a new three-noded locking-free nonlinear curved beam element. Based upon nonlinear generalized curved beam theory, large deformations and rotations in the new element are formulated in terms of Updated Lagrangian framework. Consistently coupled polynomial displacement fields are used to satisfy the membrane locking-free condition and the requirement of being able to recover the inextensible bending modes. Quintic transverse displacement interpolation functions are used to represent the bending deformation of the beam, while the axial and torsional displacement fields are derived by integration of the presumably linear membrane and torsional shear strain fields, which are coupled with the transverse displacement fields. Numerical results are presented to demonstrate the superior accuracy and the high convergence rate of the newly developed curved beam element. The stability and accuracy of the new element are further validated by experiments of an instrumented free-swinging steel cable experiencing slack and low tension. Good agreements in cable position and tension are observed between the experimental results and the finite element predictions.     

A new super convergent thin walled composite beam element for analysis of box beam structures

In this paper, a new composite thin wall beam element of arbitrary cross-section with open or closed contour is developed. The formulation incorporates the effect of elastic coupling, restrained warping, transverse shear deformation associated with thin walled composite structures. A first order shear deformation theory is considered with the beam deformation expressed in terms of axial, spanwise and chordwise bending, corresponding shears and twist. The formulated locking free element uses higher order interpolating polynomial obtained by solving static part of the coupled governing differential equations. The formulated element has super convergent properties as it gives the exact elemental stiffness matrix. Static and free vibration analyses are performed for various beam configuration and compared with experimental and numerical results available in current literature. Good correlation is observed in all cases with extremely small system size. The formulated element is used to study the wave propagation behavior in box beams subjected to high frequency loading such as impact. Simultaneous existence of various propagating modes are graphically captured. Here the effect of transverse shear on wave propagation characteristics in axial and transverse directions are investigated for different ply layup sequences.     

Analysis of laminated glass beams for photovoltaic applications

Laminated glass beams and plates are widely used in glazing and photovoltaic applications. One feature of these structures is a relatively thin and compliant polymeric layer for embedding solar cells. Proper design of photovoltaic glass modules requires an analysis of transverse shear strain distribution in polymeric encapsulant. In this paper a three layered beam with glass skins and a polymeric core is applied as a model structure to evaluate the mechanical properties. Robust relationships between the maximum deflection, the transverse shear strain of the core layer and the applied force in a three-point-bending test of laminated glass beam samples are derived. The first order shear deformation beam theory and a layer-wise type beam theory are applied. An expression for the transverse shear stiffness of the laminated glass beam is presented. The results for the maximum deflection are compared with the results discussed in the literature. Furthermore, a three-dimensional finite element analysis is performed to verify the applied beam theories. Three-point-bending tests for laminated glass beams with core layers from different polymeric materials are performed. The experimental data for the maximum deflection are compared with the derived expressions.     

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.