一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法
【技术领域】
[0001] 本发明提供一种预测含初始缺陷复合材料奖叶剩余寿命寿命的渐进损伤方法，属 于复合材料结构疲劳可靠性技术领域。
【背景技术】
[0002] 复合材料W质量轻、强度高等优点被广泛地应用于航空航天领域，尤其是在直升 机结构的设计中，复合材料已占据了直升机的主要部分。而奖叶作为直升机的主要承力构 件，在一些恶劣工作环境下使用时，很容易受到损伤，且其疲劳性能关系到整个直升机的飞 行安全。现今对含初始缺陷复合材料结构疲劳寿命预测的方法主要是疲劳试验和有限元计 算，试验方法需试样和试验成本高，且不方便设计使用；而一些有限元计算方法，需要通过 大量的试验确定模型参数，导致预测成本高，计算复杂，误差大，计算精度和效率低；目前尚 缺乏预测含初始缺陷复合材料奖叶结构疲劳寿命的简便方法，为此，本发明提出一种预测 含初始缺陷复合材料奖叶结构疲劳寿命的方法，该方法简单实用、需要的模型参数少、计算 精度高，并能合理给出奖叶结构在疲劳载荷作用下的裂纹扩展路径等优点，具有重要的学 术意义和工程应用价值。
【发明内容】
[0003] 1、目的:本发明的目的是提供一种预测含初始缺陷复合材料奖叶剩余寿命寿命的 渐进损伤方法，该方法简单实用、需要的模型参数少、计算精度高，并能合理给出奖叶结构 在疲劳载荷作用下的裂纹扩展路径等优点。
[0004] 2、技术方案:本发明提供一种预测含初始缺陷复合材料奖叶剩余寿命寿命的渐进 损伤方法，该方法具体步骤如下：
[0005] 步骤一、建立含不同初始缺陷复合材料奖叶结构有限元模型
[0006] 复合材料奖叶结构主要由碳纤维复合材料蒙皮和两根复合材料单向带大梁组成， 具体结构如图1所示。为提高计算精度，有限元建模时，奖叶蒙皮和大梁分别采用solid46层 单元和solid45^维实体单元进行模拟。初始缺陷包括中屯、圆孔、两侧缺口 W及冲击损伤的 Ξ种不同缺陷类型。复合材料奖叶结构承受载荷W压强的形式施加到奖叶的上下翼面，复 合材料奖叶结构边界条件为奖叶端面固支，如图2所示。
[0007] 步骤二、含不同初始缺陷复合材料奖叶结构渐进损伤分析
[000引含不同初始缺陷复合材料奖叶结构渐进损伤分析主要包括：（1)疲劳载荷下结构 有限元应力分布计算(2)单元剩余强度失效判别准则(3)失效单元材料性能退化方式。具体 为：
[0009] (1)疲劳载荷下结构有限元应力分布计算。采用步骤一中的Ξ维有限元模型，可W 准确得到奖叶结构的Ξ维应力分布状态。由于奖叶属于柔性结构，使用过程中会有几何大 变形产生，因此应力分析时需考虑几何大变形对应力分布的影响。
[0010] (2)单元剩余强度失效判别准则。复合材料结构受疲劳载荷加载时，其剩余强度是 随着循环次数的增加而逐渐减小的。本文采用应力控制剩余强度的方程计算各单元的剩余 强度值，具体表达式为
[00川 n = C(s-S0r[R0-R(n)]b (1)
[0012] 式（1)中η为单元疲劳寿命，S为疲劳应力，So为拟合疲劳极限，Ro为初始剩余强度 值，R(n)为经过η次循环后材料的剩余强度，m和C为由疲劳试验拟合参数。应力控制剩余强 度的方程可W由复合材料层合板疲劳试验得到。进行渐进损伤分析时，若单元应力小于其 剩余强度，则认为单元失效。
[0013] (2)失效单元材料性能退化方式。若单元未发生失效，而疲劳载荷循环次数不断增 加，则按照式（1)对未失效单剩余强度进行退化;若但元发生失效，即单元各方向应力小于 其剩余强度，则用生死单元技术杀死失效单元，即单元刚度降为0。
[0014] 渐进损伤分析具体流程算法如图3所示:输入奖叶单元类型、单元属性、几何尺寸、 边界条件、疲劳应力等，建立奖叶有限元模型;然后计算所有复合材料单元各主方向应力分 布;使用剩余失效准则，对每个奖叶单元进行失效判定。如果奖叶单元无任何失效发生，贝U 增加循环增量Ni，用疲劳剩余寿命模型式（1)计算该单元各主方向剩余强度，并采用材料性 能逐渐退化方法进行材料剩余强度退化;如果奖叶单元有失效发生，则杀死失效单元，令其 刚度将为0，并判定奖叶结构是否发生彻底破坏，若奖叶结构没有彻底破坏，则在此循环次 数下进行单元应力分析，继续判定有无新的单元失效发生，直到该循环下奖叶结构无新的 失效单元出现为止;增加循环增量Ni并按照式（1)对复合材料单元进行剩余强度退化，进行 下一次循环;当判定奖叶结构发生彻底破坏，则将当前循环数作为奖叶结构疲劳寿命。
[0015] 步骤Ξ、含不同初始缺陷复合材料奖叶结构剩余寿命累计损伤计算
[0016] 根据步骤一和二，计算不同恒幅疲劳载荷水平下的疲劳寿命。奖叶结构的载荷-寿 命表达式如式(2)
[0017]
(2)
[0018] 式(2)中，F为疲劳载荷最大值，N为计算得到的奖叶疲劳寿命，mi和Cl为拟合参数。 式(2)中的拟合参数可W通过最小二乘法拟合得到。
[0019] 由Miner累积损伤理论
[0020]
巧）
[0021] 式(3)中m为疲劳载荷谱中第i级载荷的循环次数，Ni为奖叶结构在第i级疲劳载荷 下发生失效的循环次数，j为疲劳载荷谱中载荷级数，T为载荷谱块数。
[0022] 将式(1)代入式(3)可W得到复合材料奖叶结构谱载下的疲劳寿命为
(4)
[0023]
[0024] 3、优点及功效:本发明提供一种预测含初始缺陷复合材料奖叶剩余寿命的渐进损 伤方法，其特点是简单实用、操作方便、计算精度高。
[0025] 附图表说明
[00%]图1为复合材料奖叶结构示意图。
[0027]图2为复合材料奖叶结构载荷边界条件示意图。
[00%]图3为本发明所述流程框图。
[0029] 图中符号说明如下：
[0030] 图3中的Ni为奖叶有限元模型中第i个单元的循环次数，F为疲劳载荷，N为复合材 料奖叶结构不同恒幅疲劳载荷寿命。
【具体实施方式】
[0031] 图3为本发明所述方法的流程框图，本发明分Ξ大步骤实现，具体为：
[0032] 步骤一、建立含不同初始缺陷复合材料奖叶结构有限元模型
[0033] 复合材料奖叶结构主要由碳纤维复合材料蒙皮和两根复合材料单向带大梁组成， 具体结构如图1所示。为提高计算精度，有限元建模时，奖叶蒙皮和大梁分别采用solid46层 单元和solid45^维实体单元进行模拟。初始缺陷包括中屯、圆孔、两侧缺口 W及冲击损伤的 Ξ种不同缺陷类型。复合材料奖叶结构承受载荷W压强的形式施加到奖叶的上下翼面，复 合材料奖叶结构边界条件为奖叶端面固支，如图2所示。
[0034] 步骤二、含不同初始缺陷复合材料奖叶结构渐进损伤分析
[0035] 含不同初始缺陷复合材料奖叶结构渐进损伤分析主要包括：（1)疲劳载荷下结构 有限元应力分布计算(2)单元剩余强度失效判别准则(3)失效单元材料性能退化方式。具体 为：
[0036] (1)疲劳载荷下结构有限元应力分布计算。采用步骤一中的Ξ维有限元模型，可W 准确得到奖叶结构的Ξ维应力分布状态。由于奖叶属于柔性结构，使用过程中会有几何