复合材料航天结构的应力分析.pdfVIP

4 力学与工程一杜庆华院士八十寿展庆贺史集北京：清华大学出版社．1999 复合材料航天结构的应力分析 薛明德 胡亍 曹 阳 程乐锦 (清华大学工程力学系， lo0084，北京) 摘要本支针对复合材料航天结构提出一些有效的算法．包括在三维实体wilson非协 调元基础上，构造能适用于多层，变厚度壳的三层32结点相对自由度壳单元：可是活反 映材料铺层方向变化的材料弹性矩阵与压力的算法：以压可实现持贝与旋转的子结构算 法等．实例分析说明所提出的算法建可靠和高技的，在微机上实施了对火箭仅嚣魁的应 力分析。 关键词复合材料；航天蛄构：有限单元法 一、引言 航天结构要承受起飞与着陆时的巨大推力、侧压与热载，在空间长期安全可靠地 工作，但又要尽量降低结构重量。复合材料具有高比强度、高比刚度、能经受各种复 杂恶劣环境条件和防热隔热的优良性质，被广泛用于航天结构。但其使用给结构分析 带来许多复杂的问题。复合材料具有很强的各向异性，材料主方向随结构内点位变 化；通常由许多薄层构成，层厚与结构尺寸比极小；薄层内的拉伸模量与强度非常 高，但横剪、横拉模量很低，强度更低，故需对那些通常认为较次要的应力分量给予 准确计算；各结构的连接部位(如法兰与壳体连接处等)与开孔处需补强，应考虑不同 材料界面上应力、应变的间断。凡此种种，使其结构分析远比金属材料复杂得多。 迄今为止，对复合材料结构各国学者通常采用二维板壳元进行分析，但单层理论 虽然可反映结构的总体行为，计算稳定性与动力学问题，却无法准确反映结构横向应 力。Reddy等发展了二维多层理论，但目前往往只限于解决典型问题，如等厚度矩形 板，圆柱壳等，用于复杂的实际结构仍有许多困难。特别是对于不同材料构件连接的 应力集中区，采用各种二维单元很难准确反映引起结构破坏的局部应力；而由于各层 厚度与结构尺寸相比极薄，传统的观点认为利用三维实体元几乎无法克服巨大的计算 规模障碍。为了对复台材料与金属混合的航天结构进行分析，本文完成了以下工作： 1．选用与发展了～系列高效、高精度、可计算横向应力、便于连接的单元。 2．给出能灵活反映复合材料纤维铺层主方向变化的有限元刚度阵与应力算法。 3．采用子结构技术。几何形状相同的子结构能拷贝与旋转，大大节省了计算机时 与存储量，使超大型结构计算得以在微机上实施。 力学与工程一杜庆华院士八十寿展庆贺支集 二、几种特殊的单元 本文12^三维实体元为基础，发展能适用于变厚度多层薄壳的单元。wilson“1提出 在单元插值函数中附加内部自由度的非协调项可提高单元精度，并证明了其收敛性； 例如∞1三维8结点wilson单元可达20结点单元同级精度，结点数仅为其2，5。他们81还 将三维实体元的节点位移用位移差与平均值代替，克服了将实体元用于壳体时刚度矩 阵病态的数值困难，称之为相对自由度壳元。如】6结点三维实体单元内位移插值为 ”别牲q翰 M“…“… 此处M，vj，Ⅵ以及M是结点的绝对位移和绝对插值函数；M忙1，2，…，16)的表达式同 常规的16结点等参单元，见文12驯；a，，屈，，，(j=1，2，…，5)是单元的内部自由度 Gl={(1一f2) G2=”(1一矿) G=1一f2 G．=勿(1一{2) G，={q(1一叩2) 引入相对位移珥，巧，可如下 = “ “=1，2，3，4，9，lO，l1，12) 一坼一峨 4 ●了Jl } )“ ”地 !芏％ E，羁的表达式与上式类似。单元内点的绝对位移改为节点相对位移插值，为 {习=善冠馁}+窑G悟} (4) 式中两为相对插值函数，可由一般三维等参16结点单元的插值函数M得到： {苎。Ⅳ