北航固体力学报告.doc

固体力学课程报告 《加筋壁板结构热屈曲分析方法》 姓名： 学号： 2013年12月 加筋壁板结构热屈曲分析方法 1，绪论 高超声速技术综合了航空航天领域内诸多学科的新技术，为未来航空航天领域的重要研究发展方向，被认为是既隐身技术之后的又一重点技术。高超声飞行器以高速在大气内飞行时，周围的空气受到压缩并产生巨大的摩擦作用，使飞行器的动能大部分转化为热能。如在飞行中，机身前缘和机翼前缘等部位受到严重的气动加热，温度急剧升高，可达1400℃以上，大面积气动再入温度也在400℃～1100 复合材料加筋板是近代飞机和复合材料航天器结构中最重要的结构形式之一。该结构不仅保持了复合材料所具有的比强度高和比刚度高等优点，薄板与加强筋的组合方式还极大地提高了板壳结构的后屈曲承载能力，比单纯增加板壳厚度更为经济 当飞行器在太空飞行时，由于太阳辐射关系可以产生足够的热量，使层合壁板升温以至于发生屈曲，求知屈曲的临界温度对于热结构设计者来说是至关重要的。目前对于求解复合板的热屈曲问题主要有两种方法： 1，准确解法 当平板变形时，在x，y和z方向的位移分量一般可用幂级数来表示： （1） 以上幂级数所包括的项数，由问题准确性的要求决定。假定w与z无关，且令u与v为： （2） 式中的为中间面的位移；是根据切应力沿层合板厚度的分布给出。 对任一单层K来说，应力-应变关系为： （3） 层合板上合力和力矩可通过沿其厚度的应力积分获得： （4） 其中可求板上的热力和热力矩： （5） 欲求式（4）-（5）的合力与力矩，必须事先给定满足横向剪切应力要求的，对于层合板屈曲分析，还需要给出边界条件。 国内学者孙良新和徐泰然对正交对称矩形层合板受均匀温度分布的问题进行了计算，得到了板的屈曲温度为： （6） 解析解法虽然看起来简洁，然而计算起来并不简单，尤其是对不同的载荷和层板的尺寸和铺层角等考虑起来更加麻烦了，因此用有限元法计算便成为一个较为便捷的出路。 2，有限元分析法 孙良新教授利用最小位能原理，推导和建立了在高阶位移模式基础上的有限元公式，结构的总位能为 （7） 式中为单元总数，对于每个单元沿厚度方向积分，有 （8） 孙良新通过把有限元法与分析法计算结果和经典解进行了比较，说明了有限元法有足够的精度。本文利用有限单元法求解不同温度场环境下结构的屈曲载荷，也就相当于给定边界条件与载荷约束，可求对应的临界失稳温度。, 3，加筋板热屈曲有限元模型分析 3.1热应力分析基本原理 热应力研究在已知温度场的情况下确定模型的应力应变。设物体内部存在温差的分布，那么这个温差将引起热膨胀，其膨胀量为，为热膨胀系数（thermal expansion coefficient），则该物体的物理方程由于增加了热膨胀量将变为： （9） 其他平衡方程。边界条件和几何方程与普通弹性问题相同，弹性问题的虚功原理的一般表达式为，也就是 （10） 将上面的物理方程带入，可得 （11） 进一步可以写成 （12） 该式即为热应力问题的虚功原理。 3.2热屈曲分析基本原理 一般的特征屈曲问题是寻找载荷使得结构刚度矩阵奇异，即使得问题 有非显示解。是施加载荷后的期限刚度矩阵，即为非显示位移解，载荷可以为压力载荷、集中载荷、位移载荷以及热载荷。特征屈曲分析一般用于估计加筋结构的临界屈曲载荷。 屈曲载荷的计算依赖与结构的基本状态（the base state），若屈曲分析是结构分析的第一个分析，那基本状态就是结构的初始条件，否则基本状态是上一个分析步结束时的状态，由此，基本状态可以包含预载，预载为零时即为经典特征屈曲问题。 在特征屈曲分析里施加一个额外的载荷，其值多少并不重要，屈曲载荷通过与相乘得到，值从下式求解： （13） 式中是依赖于基本状态的刚度矩阵，包含了预载的影响； 是基于外载荷得到的起始应力和载荷刚度矩阵； 是特征值 是屈曲模态 和为模型自由度，为模态阶数； 则临界