一种基于模态应力场的制动盘疲劳裂纹快速预测方法

本发明属于机械结构疲劳寿命仿真预测，特别是涉及一种基于模态应力场的制动盘疲劳裂纹快速预测方法。

背景技术：

1、在机械结构疲劳寿命仿真预测领域，制动盘在摩擦制动过程中承受剧烈的热机耦合载荷，由此产生的复杂残余应力场是引发其表面热疲劳裂纹萌生与扩展的主要原因。为预测其疲劳寿命，当前先进方法通常采用一个物理意义明确的两阶段分析框架：首先进行完全热机耦合分析以获取制动盘的残余应力场；随后将该应力场作为载荷条件，驱动第二阶段的裂纹扩展分析。然而，该框架在工程实践中面临一个显著瓶颈：第一阶段计算所得的残余应力场是一个包含数十万乃至数百万节点应力分量的庞大空间数据集。将此海量数据完整、精确地传递并施加至第二阶段的裂纹扩展模型，不仅过程繁琐，更会致使裂纹扩展模型因初始数据规模巨大而计算效率低下，这在需要进行多参数优化或快速设计迭代的场景中，往往导致难以接受的时间成本。

2、为提升计算效率，现有技术中存在一些简化处理方式。例如，采用均匀压力、简化温度梯度等理想载荷来替代真实的残余应力场，但这类方法会严重牺牲预测结果的物理精度。另一种思路是采用子模型技术，通过截取局部区域进行分析以缩小模型规模，但在子区域选取与边界条件传递方面仍较大程度依赖工程经验，且未能从根本上解决核心分析区域（即裂纹潜在萌生与扩展区域）的残余应力数据集规模过大的问题。因此，亟待发展一种新型的残余应力场处理方法，该方法需在充分保留其影响疲劳寿命的关键力学特征的前提下，实现对原始数据规模的极大压缩，从而为后续裂纹扩展分析构建高效、精确的载荷输入，以突破当前工程仿真的效率瓶颈。

技术实现思路

1、本发明的目的在于提供一种基于模态应力场的制动盘疲劳裂纹快速预测方法，解决现有技术中因残余应力场数据规模庞大导致后续裂纹扩展分析效率低下、且现有简化方法难以兼顾精度与效率的问题，实现基于物理降阶的高效、精度可控的疲劳裂纹快速仿真预测。

2、本发明所采用的技术方案是，一种基于模态应力场的制动盘疲劳裂纹快速预测方法，包括以下步骤：

3、s1：建立制动盘与摩擦块的三维装配体模型，定义材料参数与接触属性，施加制动载荷与自然冷却边界条件，通过完全热机耦合瞬态分析模拟制动全过程，提取制动盘冷却至室温后的单元高斯积分点上的原始残余应力场；

4、s2：从所述三维装配体模型中提取制动盘单体结构，并对该单体结构进行自由状态下的模态分析，提取其前阶模态振型及对应的模态应力场，构建模态应力场基；

5、s3：基于步骤s1与步骤s2中制动盘网格的一致性，将所述原始残余应力场投影到所述模态应力场基上，基于高斯数值积分方法进行数值积分以组装线性方程组，通过求解该线性方程组得到模态参与系数向量；

6、s4：基于所述模态参与系数向量与所述模态应力场基，通过线性叠加重构降阶残余应力场，并在制动盘单体结构的单元高斯积分点上执行逐点加权求和计算，以生成重构残余应力场；

7、s5：将所述重构残余应力场施加于仅包含制动盘结构的有限元模型，执行静力分析以验证应力场平衡；若应力场平衡不满足预设条件，则调整模态阶数并返回步骤s2；待应力场平衡满足预设条件后，预设初始裂纹，采用扩展有限元法模拟裂纹在残余应力驱动下的扩展，基于裂纹扩展模型积分得到制动盘的预测疲劳寿命。

8、进一步地，s2中，

9、所述制动盘单体结构直接沿用步骤s1中制动盘部分的有限元网格；

10、通过收敛性判断策略确定模态阶数：从预设初始阶数开始，逐步增加模态阶数，利用提取的模态应力场基对所述原始残余应力场进行重构，当重构应力场与原始残余应力场在关键区域的均方根误差连续两次满足预设收敛条件时，将当前的模态阶数确定为最终阶数。

11、进一步地，所述s3具体包括：

12、将原始残余应力场与各阶模态应力场转换至同一套单元高斯积分点集；基于高斯数值积分方法，遍历制动盘结构的所有单元及其内部高斯积分点，在每个积分点处计算应力张量的双点积，并乘以该积分点的权重与雅可比行列式，对所有单元及积分点的计算结果进行求和，得到系数矩阵与右端项向量，构建线性方程组；求解该线性方程组，得到模态参与系数向量。

13、进一步地，所述s5包括：

14、基于步骤s2的制动盘单体结构建立新的制动盘有限元模型，将所述重构残余应力场作为预定义初始应力场施加至该模型的单元高斯积分点上；

15、所述静力分析在仅施加所述重构残余应力场且不施加任何其他外载荷与约束的条件下，检验制动盘模型的节点残余力范数是否小于容许阈值；若节点残余力范数大于容许阈值，则返回步骤s2增加模态阶数。

16、进一步地，所述s5还包括：

17、在重构残余应力场的最大主应力区域预设初始裂纹，采用扩展有限元法模拟裂纹在残余应力驱动下的准静态扩展；

18、在裂纹扩展过程中，采用相互作用积分法计算裂纹尖端应力强度因子，并基于paris疲劳裂纹扩展公式计算裂纹扩展速率；通过对扩展全过程进行数值积分，累加得到裂纹从初始裂纹长度扩展至临界长度的循环次数，该循环次数作为制动盘的预测疲劳寿命。

19、本发明的有益效果是：

20、1.本发明通过将需要传递的残余应力场从庞大的节点应力数据，压缩为少量模态系数与一组固定的模态基，大幅降低了裂纹扩展分析的初始数据规模、内存占用及单次求解耗时，从而显著缩短了制动盘疲劳裂纹的总体预测时间，计算效率可提升一个数量级以上，能够有效支撑多参数优化与快速设计迭代。

21、2.本发明基于结构动力学原理进行物理降阶，利用低阶模态表征结构主要的受力与变形模式；通过选取适当的模态阶数，能够在计算效率与预测精度之间实现优化平衡，确保重构应力场在最大应力值、应力梯度等关键力学特征上具有足够保真度，克服了简化载荷方法精度损失大的缺陷。

22、3.本发明所涉及的模态提取、应力场向模态空间投影、以及裂纹分析中的应力场重构步骤，均可通过脚本实现自动化操作，能够与主流商业有限元软件有效衔接，易于集成到自动化的仿真分析与设计流程中，提升了方法的工程实用性与可重复性。

23、4.本发明适用于作为一套高效、可靠的仿真预测方法进行实施与应用；其有效性与优越性可通过与传统全数据传递方法在同等条件下对比裂纹扩展路径、疲劳寿命及计算耗时等核心结果予以验证，该验证过程完全在数值仿真层面完成，不依赖于物理试验。

1.一种基于模态应力场的制动盘疲劳裂纹快速预测方法，其特征在于，包括以下步骤：

2.根据权利要求1所述的一种基于模态应力场的制动盘疲劳裂纹快速预测方法，其特征在于，s2中，

3.根据权利要求1所述的一种基于模态应力场的制动盘疲劳裂纹快速预测方法，其特征在于，所述s3具体包括：

4.根据权利要求1所述的一种基于模态应力场的制动盘疲劳裂纹快速预测方法，其特征在于，所述s5包括：

5.根据权利要求4所述的一种基于模态应力场的制动盘疲劳裂纹快速预测方法，其特征在于，所述s5还包括：