含冲击损伤复合材料加筋板的剩余强度分析.docVIP

 PAGE 7 含冲击损伤复合材料加筋板的剩余强度分析 林智育，许希武 南京航空航天大学航空宇航学院结构强度研究所，南京，210016 Email:  HYPERLINK mailto:linzhiyu@nuaa.edu.cn linzhiyu@nuaa.edu.cn Fax: 025Phone: 025摘要：采用含刚度折减椭圆形弹性核复合材料层合板低速冲击损伤模型，利用含任意椭圆核各向异性板杂交应力有限元来模拟含损伤区域，杆单元来模拟筋条，钉单元模拟铆钉或胶层，建立起含低速冲击损伤复合材料加筋层板的力学分析方法。利用基于特征曲线概念的点应力判据、最大正应力判据、最大剪应力判据分别预计蒙皮的破坏、筋条的破坏和胶层或铆钉的破坏，获得了与试验结果相吻合的结果。 关键词：复合材料、加筋板、冲击损伤、剩余强度 引言 复合材料由于其具有的优良特性在现代飞机结构中得到了广泛的应用。复合材料加筋层板结构是航空航天等工程结构中常用的主承力部件。加筋层板经常会在制造过程中和使用过程中由于外来物的低能量冲击损伤，这些冲击损伤的存在将严重削弱加筋层板的承载能力，特别是压缩承载能力，给飞机结构的安全使用构成严重威胁。含冲击损伤复合材料加筋层板在复杂载荷作用下的剩余强度一直是人们关注的研究课题之一。 本文采用椭圆形弹性核复合材料低速冲击损伤模型，将冲击损伤简化为一椭圆形弹性核。采用含任意椭圆核各向异性板杂交应力有限元（模拟含损伤区域），结合杆单元（模拟筋条）、钉单元（模拟胶层或铆钉）来分析含冲击损伤复合材料加筋壁板的应力分布，采用基于特征曲线概念的点应力判据、最大正应力判据、最大剪应力判据分别预计蒙皮的破坏、筋条的破坏和胶层或铆钉的破坏，获得了与试验结果相吻合的结果。 含弹性核有限大复合材料加筋板的应力分析 采用含刚度折减椭圆形弹性核复合材料层合板低速冲击损伤模型，将含冲击损伤复合材料加筋壁板的应力分析问题化归为含椭圆形弹性核复合材料加筋壁板的应力分析问题进行求解。在建模分析过程中，采用含任意椭圆核各向异性板杂交应力有限元[1]（模拟受损伤板区域）结合以杆单元（模拟筋条）和钉单元（模拟铆钉或胶层）建立起应力分析模型。 含任意椭圆核各向异性板杂交应力有限元 单元内假设应力场和位移场 单元内的应力场和位移场采用满足平衡方程、几何方程和物理方程的复变函数级数解，该级数解在椭圆核边界处满足位移协调条件和应力连续条件。 含椭圆孔有限大各向异性板，当孔边应力主矢为零时，文献[2]给出了其复势函数表达式： （1） 在椭圆孔边界上有，其复势函数在孔边的边界值为： （2） 考虑椭圆形各向异性板，文献[3]中给出了其复势函数： （3） 将边界值代入复势函数表达式，可以得到椭圆形各向异性板复势函数的边界值： （4） 设在椭圆形孔边作用有应力或位移，则其边界条件可以写为： （5） 位移边界条件： ， 应力边界条件： ， 其中：外边界取正值，内边界取负值。 有限元单元列式 单元依照板的能量泛函[4]导出，其表达式为： （6） 式中是板所占的区域，是板的边界，是已知的应力边界，位移向量是定义在边界处，应力分量在边界处已知，位移向量和应力向量定义在内，是柔度矩阵。 由于假设的应力场满足平衡方程，式(6)可以简化为： （7） 单元列式时，单元内应力场、位移场以应力参数进行插值，单元边界处位移场以节点位移进行插值 （8） 式中和分别为应力参数和位移参数矩阵，是由在单元边界处的形状函数矩阵，本文中采用常规有限元的形状函数矩阵。 将上式代入泛函式（7），有 （9） 为等效节点力。在单元上缩聚掉参数 ， （10） 即得到单元刚度矩阵 （11） 钉单元列式 钉单元指结构连接件中的螺栓和铆钉，主要用于受剪状态的模拟。本文中将钉元推广应用于加筋壁板中的铆接和胶接的模拟。在结构钉元局部坐标系下的钉元刚度矩阵[5]为： （12） 其中：k为钉元的刚度系数。 其受力状态方程可写为： （13） 式中：S为钉元两端的剪应力列阵；K为结构坐标系下钉元的刚度矩阵；δ为钉元节点的位移阵。 含冲击损伤复合材料加筋壁板的剩余强度分析 分别采用基于特征曲线概念的点应力判据、最大正应力判据和最大剪应力判据来分别预计蒙皮的破坏、筋条的破坏和铆钉或胶层的破坏。 3.1 蒙皮的剩余强度失效判据 根据应力分析结果，将特征长度与点应力判据相结合，采用Yamada-Sun准则，来预测含损伤层板的剩余强度。 Yamada-Sun椭圆准则的失效判据为： （14） σ1和τ12分别为每一层的正轴正应力和剪应力。X为单层板纵向强度，当σ1＞0时，为纵向拉伸强度XT；当σ1＜0时，为纵向压缩强度XC