一种复合材料胶接修补结构的剩余疲劳寿命的预测方法技术

本发明专利技术提供了一种复合材料胶接修补结构的剩余疲劳寿命的预测方法，步骤为：1）确定未修复结构的材料常数；2）对修复结构进行建模；模型建立后，进行实体装配、划分网格、创建约束；3）求解应力强度因子幅值与裂纹长度之间的关系式；4）根据金属材料的断裂韧性，基于软件输出的应力强度因子确定修复结构的最终失效裂纹长度；5）把材料常数值、最终失效裂纹长度、应力强度因子变化幅值和裂纹长度的关系式代入Paris公式，积分求得修复结构的剩余疲劳寿命。本发明专利技术简化了计算步骤，显著提高了复合材料胶接修复结构的剩余疲劳寿命预测效率，降低计算成本，促进复合材料胶接修复技术的发展。

Method for predicting residual fatigue life of composite adhesive repair structure

The invention provides a prediction method, the residual fatigue life of a composite patch structure comprises the following steps: 1) to determine the material constants of unrepaired structure; 2) modeling of the repair structure; model, entity assembly, grid, create constraints; 3) to solve the relationship between stress intensity the amplitude factor and crack length; 4) according to the fracture toughness of metallic materials, the software output of the stress intensity factor of the crack length to determine the final failure repair based on the structure; 5) the material constant value, the final failure relation into Paris formula of crack length, stress intensity factor amplitude and crack length are obtained by integrating. The residual fatigue life of the repaired structure. The invention simplifies the calculation steps, and significantly improves the prediction efficiency of the residual fatigue life of the composite adhesive repair structure, reduces the calculation cost, and promotes the development of the adhesive repair technology of the composite material.

【技术实现步骤摘要】

本专利技术属于复合结构疲劳寿命预测
，具体涉及一种基于材料常数不变假设条件的复合材料胶接修补结构剩余疲劳寿命预测方法。
技术介绍
断裂力学中常用Paris公式描述裂纹扩展速率与裂纹尖端应力强度因子变化幅值关系，其表达式为其中为疲劳载荷作用一个循环裂纹扩展的长度，即裂纹扩展速率；ΔK为裂纹尖端的应力强度因子；C和m为待定常数，与材料，环境等等有关。采用Paris公式积分求解含裂纹金属板复合材料胶接修补结构的剩余疲劳寿命过程中，如何确定材料常数C和m，是该方法中的一个关键和难点。现有技术的常用做法是采用试验测试得到修补结构的裂纹扩展速率，结合有限元方法求得修复结构的应力强度因子幅值ΔK，将Paris公式等式两边取自然对数，求得试验条件下的材料常数C和m。但是该方法要求在求解剩余疲劳寿命之前首先通过实验得到修复结构的疲劳寿命试验数据，在很多情况下针对复杂的复合材料胶接修补结构该类实验是无法实施的，而且即便可以实施也需要反复计算修复结构中的应力强度因子幅值，步骤繁琐，计算量大。
技术实现思路
本专利技术的目的是解决剩余疲劳寿命预测的技术问题，应用Abaqus有限元软件平台，基于Paris公式，假定修复前后金属板的材料常数C和m不变，结合有限元得到的修复结构中的应力强度因子值，提供一种切实可行的复合材料胶接修补结构的剩余疲劳寿命预测方法，步骤为：1)在等幅交变载荷作用条件下，对未修补的含裂纹复合材料的疲劳行为进行试验，得到裂纹扩展速率与裂纹长度的关系，然后采用多项式拟合，得到Paris公式中的材料常数；2)对含裂纹的复合材料进行胶接修复，然后在Abaqus有限元软件平台上对上述修复结构进行有限元建模；模型建立后，进行实体装配、划分网格、创建约束；3)求解应力强度因子幅值与裂纹长度之间的关系式；4)根据金属材料的断裂韧性，基于软件输出的应力强度因子确定修复结构的最终失效裂纹长度，当修复结构的断裂韧性等于金属板的断裂韧性时对应的裂纹长度，为修复结构的最终失效裂纹长度；5)把材料常数值、最终失效裂纹长度、应力强度因子变化幅值和裂纹长度的关系式代入Paris公式，积分求得修复结构的剩余疲劳寿命。进一步的，所述步骤2)中建模是根据修复结构中各材料的特点，从Abaqus单元库中选择适合各自结构特点的单元类型，进行建模。进一步的，所述步骤2)在划分网格的过程中，裂纹处进行网格密度细化。进一步的，所述步骤2)在创建约束过程中不同组件的界面节点处采用绑定约束；边界条件采用金属板一端固定，另一端施加载荷的约束形式。进一步的，所述步骤3)的具体步骤是：在预置裂纹处的两个界面处，根据结构特点，设置可在载荷作用下分离的代表裂纹的缝合面Seam，定义为初始裂纹，通过Abaqus软件的断裂力学模块，输出裂纹尖端的应力强度因子；改变裂纹长度，输出不同裂纹长度下的裂纹尖端的应力强度因子；通过输出修复结构在最大载荷和最小载荷作用下应力强度因子的大小，求得裂纹尖端的应力强度因子变化幅值；通过拟合应力强度因子幅值与裂纹长度之间的关系，得到两者的表达式。与现有技术相比，本专利技术的优点和积极效果是：本专利技术以未修复板的材料常数作为修复结构的材料常数，提出了一套切实可行的剩余疲劳寿命预测方案。只需通过试验方法获得结构相对简单的未修补板的材料常数，结合有限元软件输出得到的应力强度因子幅值，并通过多项式拟合应力强度因子幅值ΔK与裂纹长度a的关系，即可实现对修复结构的剩余寿命的预测，使得步骤大为简化，可显著提高复合材料胶接修复结构的剩余疲劳寿命预测精度并且降低计算成本，为该修复技术的推广提供理论支撑。附图说明图1.本专利技术的流程示意图；图2.实施例1的网格划分效果图；图3.实施例2的网格划分效果图。具体实施方式下面结合附图及具体实施方式对本专利技术的技术方案作进一步详细的说明。本实施例的技术方案是：1)确定未修复结构的材料常数：在交变载荷作用条件下，对未修补的含裂纹复合材料进行疲劳加载，测得损伤未修补裂纹板的裂纹扩展速率并通过查阅应力强度因子手册得到对应的含裂纹损伤板的ΔK，在双对数坐标系下拟合所得到的数据点得到Paris公式中在本试验条件下的材料常数C和m。2)对修复结构进行建模：对含裂纹的复合材料进行胶接修复，然后在Abaqus有限元平台上对上述修复结构进行建模。根据修复结构中各材料的特点，从Abaqus单元库中选择适合各自结构特点的单元类型，进行建模。当各组件的模型建立后，进行实体装配、划分网格、创建约束。需要强调的是，在划分网格的过程中，裂纹处要进行网格密度细化，经过反复尝试，建议裂纹附近的网格划分尺寸为0.08-0.12mm，可以满足计算精度要求又不至于带来过大的计算量。考虑复合材料胶接修补结构为不同组件的组合体，在不同组件的界面节点处采用绑定约束，以模拟各组件之间协同变形的状态。边界条件为裂纹板一端施加固定约束，另一端施加载荷。3)求解复合材料胶接修复结构中应力强度因子幅值ΔK与裂纹长度a之间的关系式ΔK～a。其具体步骤是：在预置裂纹处的两个界面处，根据结构特点，设置可在载荷作用下分离的代表裂纹的缝合面Seam，将其定义为初始裂纹。通过Abaqus软件的断裂分析模块，直接输出裂纹尖端的应力强度因子；改变裂纹长度，即可以得到不同裂纹长度下的裂纹尖端的应力强度因子。通过输出修复结构在最大载荷和最小载荷作用下应力强度因子Kmax和Kmin，即可求得裂纹尖端的应力强度因子变化幅值ΔK＝Kmax-Kmin。通过拟合，得到应力强度因子幅值ΔK与裂纹长度a之间的关系表达式ΔK～a。4)确定修复结构最终失效时的裂纹扩展临界值：通过实验或者查阅相关材料力学性能手册可以获得所修复的金属材料对应的断裂韧性KIC，根据第3)步得到的ΔK～a，当ΔK＝KIC时所对应的裂纹长度a即为修复结构失效时所对应的裂纹扩展临界值acr。5)预测修复结构的剩余疲劳寿命：把步骤1)～4)得到的C和m、ΔK～a关系式代入式(1)，以考察的当前裂纹长度a0为积分下限，以裂纹扩展临界值acr为积分上限，对Paris公式进行积分，求得具有长度为a0裂纹的复合材料胶接修补结构的剩余疲劳寿命。实施例一采用T300/E51碳纤维复合材料补片对含中心裂纹LY12CZ铝合金金属损伤板进行单面胶接修补，预测修复结构的疲劳寿命。首先，测试在交变载荷作用下，损伤裂纹板裂纹扩展速率与应力强度因子幅值的关系，在双对数坐标系下拟合数据点得到在本试验条件下的Paris公式材料常数C和m，计算结果为C＝6.3496×10-12,m＝3.72972。其次，对含裂纹复合材料单面胶接修补结构进行有限元建模。根据修补结构中各材料的特点，从Abaqus软件库中选出适合金属板、胶层和复合材料补片的三维实体单元C3D8R、连续壳单元SC8R和三维粘性单元COH3D8进行实体建模，根据实际结构中各结构的位置，对实体进行装配，对模型进行网格划分，布局全局种子，并细化裂纹处的网格密度，将裂纹附近的网格尺寸划分为0.1mm，网格划分效果如图2所示。图2中，金属板1与胶层、胶层与复合材料2补片的接触面采用绑定约束。损伤金属板一端固定，另一端施加载荷。接着，求解应力强度因子幅值与裂纹长度之间的关系式。其具体步骤是：在模型中心垂直载荷方向设置一个长度本文档来自技高网...

【技术保护点】
一种复合材料胶接修补结构的剩余疲劳寿命的预测方法，其特征在于，步骤为：1）在等幅交变载荷作用条件下，对未修补的含裂纹复合材料的疲劳行为进行试验，得到裂纹扩展速率与裂纹长度的关系，然后采用多项式拟合，得到Paris公式中的材料常数；2）对含裂纹的复合材料进行胶接修复，然后在Abaqus有限元软件平台上对上述修复结构进行有限元建模；模型建立后，进行实体装配、划分网格、创建约束；3）求解应力强度因子幅值与裂纹长度之间的关系式；4）根据金属材料的断裂韧性，基于软件输出的应力强度因子确定修复结构的最终失效裂纹长度，当修复结构的断裂韧性等于金属板的断裂韧性时对应的裂纹长度，为修复结构的最终失效裂纹长度；5）把材料常数值、最终失效裂纹长度、应力强度因子变化幅值和裂纹长度的关系式代入Paris公式，积分求得修复结构的剩余疲劳寿命。

【技术特征摘要】
1.一种复合材料胶接修补结构的剩余疲劳寿命的预测方法，其特征在于，步骤为：1）在等幅交变载荷作用条件下，对未修补的含裂纹复合材料的疲劳行为进行试验，得到裂纹扩展速率与裂纹长度的关系，然后采用多项式拟合，得到Paris公式中的材料常数；2）对含裂纹的复合材料进行胶接修复，然后在Abaqus有限元软件平台上对上述修复结构进行有限元建模；模型建立后，进行实体装配、划分网格、创建约束；3）求解应力强度因子幅值与裂纹长度之间的关系式；4）根据金属材料的断裂韧性，基于软件输出的应力强度因子确定修复结构的最终失效裂纹长度，当修复结构的断裂韧性等于金属板的断裂韧性时对应的裂纹长度，为修复结构的最终失效裂纹长度；5）把材料常数值、最终失效裂纹长度、应力强度因子变化幅值和裂纹长度的关系式代入Paris公式，积分求得修复结构的剩余疲劳寿命。2.根据权利要求1所述的一种复合材料胶接修补结构的剩余疲劳寿命的预测方法，其特征在于，所述步骤2）中建模是根据修复结构中各材料的特点，从Abaqus单元库中选...

【专利技术属性】
技术研发人员：穆志韬，金平，李旭东，刘治国，孔光明，马鸿儒，
申请(专利权)人：中国人民解放军海军航空工程学院青岛校区，
类型：发明
国别省市：山东;37