复合材料柔性梁剖面刚度特性分析的哈密尔顿求解体系.pdfVIP

复合材料柔性梁剖面刚度特性分析的哈密尔顿求解体系。 张呈林刘勇 (南京航空航天大学飞行器系南京210016) 摘要本文提出一种基于哈密尔顿体系的求解复台材料柔性粱剖面剐度特性的有限 元列式。充分考虑了复合材料柔性梁的变形特点，分析时采用了Giavotto等人“1的翘曲位移 模式，对于翘曲位移来加任何限制，直接引入三维翘曲t完整地保留截面的三维直力状态，从 而得到完整的应变能列式。井充分利用哈密尔顿体系的特性，得到了辛型的求解算法，算蜊 表明了方法的高精度。 关键调复台材料柔性粱剖面特性哈密尔顿系统有限元 1引言 现在，先进的直升机大多采用了复合材料桨叶．无轴承旋翼(无轴承尾桨)的应用 大大简化，优化了直升机结构。同时，也对分析方法提出了新的要求。复合材料由于其 各向异性的弹性性质将拉伸和各个剪切都关联起来了，故复合材料浆升柔性粱段的扭转、 弯曲及拉伸的总体运动形态是互相耦合作用的，翘曲(分面内和面外)及其耦台不容忽 视；再加之其大扭转的工作特点，就导致了复合材料柔性粱分析的复杂性。 目前常用的成功的做法“”1是：将桨叶模拟为梁．根据渐次变分原理可以将整个同题 分解为线性的二维剖面特性分析和一维的非线性运动分析。其中削面特性分析是整个筒 题的基础。～方面通过二维特性分析可以确定柔性梁的剖面刚度特性；另一方面可以作 为一种通用的分析复合材料梁剖面特性的工具，为后面的一维非线性分析提供本构关系． 复合材料柔性梁其剖面刚度特性是用 6×6的对称矩阵描述的，它一方面包含了 材料铺层的弹性耦台，另一方面也包含了 ；．．一埘滢渤驾。。。0 柔性梁剖面翘曲变形(及其耦合)所贡献 的弹性。其核心是对于翘曲的分析。二维 剖面线性分析十几年来一直沿袭Giavotto 等人“]发展的有限元技术，采用桨叶应变 能表示面内外翘曲函数，由二维有限元离 散导向二阶常系数微分方程，方程的特征 图1，坚标采 值给出截面翘曲。但是，钟万勰等m”最近的研究表明，Giavotto等人的模型是属于弹性 力学的哈密尔顿理论体系的，Giavotto等人未能自觉地利用哈密尔顿体系的很多优良的 ·车文得到航空酃“九五”、课题和航空基金资助 638 特性，如守恒等，以致使原本守恒的体系人为地化入了耗散体系，大大降低了问题的求 解精度。哈密尔顿体系的研究表明，将传统的拉格朗日体系转入哈密尔顿体系，在新的 pleetie)矩阵中自然地包含了系统全部的属性，这两类矩阵的特征分析也有相应的保结构 算法。 本文分析时采用了Giavotto等人的翘曲位移模式，完整地保留截面的三维应力状 态，从而得到完整的应变能列式，以分析翘曲对于剖面刚度特性，特别是对于剪切刚度 的影响。分析中采用了有限元等参列式。 2复合材料梁拉、弯、扭圣维南问题的哈密尔顿列式 等，应变、位移、应力向量为： ￡={‘勺y”ynt。￡。}7={矿扩)7；目={uvw)7 口；{以oo，o：o。口。}7={，p7)7 (1) 应变一位移关系为： 鑫。 0 箬 o 8 a o ayaT s—Eo口+El口，Eo一 ^一殴]^一鼬㈣ 0 0 0 o o 0 a一如a一如o 应力一应变关系为： ‘56‘s3] c51。52f54c55 f56] 。：；c。，。。：：。，，c：：：I!c66‰l，e一|‰‰％‰‰l