考虑表面改质的地下工程构件微动裂纹扩展速率预测方法与流程

本发明涉及地下矿井结构健康监测、断裂力学评估、表面改性力学分析，具体涉及一种考虑表面改质的地下工程构件微动裂纹扩展速率预测方法，更涉及一种结合表面弹性基本解、微动接触理论与疲劳断裂理论，对带有表面改质层的地下工程构件微动裂纹扩展速率及剩余寿命进行定量预测的方法。

背景技术：

1、隧道、综合管廊、矿山巷道、深基坑等地下工程中的装配式构件及连接节点，长期承受围岩蠕变、地层交变压力、施工爆破和交通振动等复杂荷载作用。构件接触面在法向集中力、界面剪切力与交变地应力耦合作用下，发生微米级往复滑移，形成典型的微动现象。微动过程中，接触边缘区域易萌生微裂纹；在地下水及侵蚀性离子的协同腐蚀作用下，裂纹持续扩展，最终可能导致构件开裂、结构渗漏甚至支护失稳等重大安全事故。为改善构件接触面的耐磨、抗滑移及抗疲劳性能，工程上普遍采用表面改质处理技术，如表面滚压强化、改性涂层喷涂、渗碳处理及喷砂硬化等，在构件表层形成力学性能异于基体的改质层。该改质层具备相对独立的表面弹性与表面刚度，能够显著改变集中力作用下的表面应力分布，消除传统弹性解中的应力奇异性，从而影响界面切应力、微动滑移幅值及微动裂纹的萌生与扩展行为。

2、然而，现有地下工程微动裂纹预测技术仍存在以下不足：其一，现有表面改质特征长度的计算模型仅引入改性层与基体剪切模量单一参数，未考虑滚压工艺塑性变形引入的表层残余压应力强化效应，同时忽略微动接触界面表面粗糙度带来的局部应力集中弱化作用；单一弹性参数无法表征滚压强化、粗糙微动接触、界面剪切刚度多因素耦合作用，采用该模型计算得到的改性特征长度与地下工程构件微动裂纹实际扩展行为偏差较大，裂纹扩展速率预测精度低。其二，传统flamant平面集中力基本解未考虑表面改质层的表面弹性效应，存在应力奇异性，导致计算得到的接触应力与切应力偏离构件实际受力状态，影响后续裂纹预测的准确性。其三，现有微动裂纹预测模型多仅基于基体材料力学参数，未建立表面改质参数-接触应力-微动滑移-裂纹扩展的定量耦合关系，无法适配带改性表层的工程构件。其四，现有方法多依赖经验取值与常规疲劳模型，预测精度较低，难以有效指导表面改质工艺优化及结构运维决策。

3、有鉴于此，亟需结合gurtin-murdoch表面弹性理论、平面表面改质集中力基本解、微动力学模型及疲劳裂纹扩展理论，构建出一套适用于带表面改质层地下构件的微动裂纹预测方法，具有重要的理论价值与工程应用前景。

技术实现思路

1、鉴于背景技术中存在的技术问题，本发明提供了一种考虑表面改质的地下工程构件微动裂纹扩展速率预测方法，引入平面表面改质集中力基本解修正接触应力场，建立表面改质参数、界面接触应力、微动滑移幅值、等效应力幅与裂纹扩展速率的全链条定量计算模型，实现带表面改质层地下构件微动裂纹扩展速率、裂纹发展趋势及剩余疲劳寿命的精准预测，同时可为构件表面改质工艺参数优化、微动损伤防控提供理论依据，适用于深部矿井锚杆、管片螺栓等关键承力构件的健康监测与防护。

2、第一方面，本发明实施例提供了一种考虑表面改质的地下工程构件微动裂纹扩展速率预测方法，其包括如下步骤：

3、s1、参数采集与标定：对经表面改质处理的地下工程微动易发构件进行现场工况参数实测、基体与表面改质层力学参数修正与标定、界面与疲劳参数标定；并基于所述力学参数计算得到表面改质特征长度、临界表面改质特征长度；

4、s2、接触应力场求解：基于表面改质层力学模型，采用平面表面改质集中力修正基本解，计算构件接触面法向应力与极坐标径向应力，获得改质层作用下的真实接触应力分布；

5、s3、界面微动参数计算：由所述接触面法向应力结合摩擦模型计算界面摩擦力，求解界面平均切应力，并依据滑移理论计算界面微动滑移幅值；

6、s4、微动裂纹萌生模式判定：将表面改质特征长度与临界表面改质特征长度进行比较，判定微动裂纹的萌生位置及开裂模式，选择对应的微动疲劳常数；

7、s5、微动耦合力学参数计算：将远场交变地应力与界面切应力进行耦合，得到综合等效应力幅，结合几何形状修正系数与瞬时裂纹长度，求解裂纹尖端应力强度因子幅；

8、s6、裂纹扩展速率求解：将应力强度因子幅与对应模式的微动疲劳常数代入微动工况适配的paris公式，计算当前载荷循环下的裂纹扩展速率；

9、s7、剩余疲劳寿命计算：以初始裂纹长度为积分下限、临界失效裂纹长度为积分上限，对所述裂纹扩展速率进行积分，得到构件的剩余疲劳寿命；

10、s8、分级预警与输出：根据所述裂纹扩展速率和/或剩余疲劳寿命进行分级预警，制定差异化处置方案并输出评估报告。

11、作为本发明的进一步改进，步骤s4中，当表面改质特征长度≤临界表面改质特征长度时，判定最大切应力位于构件表面，裂纹以表面开口型模式萌生，选用第一组微动疲劳参数；当表面改质特征长度>临界表面改质特征长度时，判定最大切应力迁移至构件次表面，裂纹以亚表面剪切型模式萌生，选用第二组微动疲劳参数。

12、作为本发明的进一步改进，所述表面改质特征长度的计算过程为：基于基体剪切模量与表面改质层剪切模量的比值，确定理想表面改质特征长度；引入表层残余压应力的强化修正，得到残余应力修正后的表面改质特征长度；引入表面粗糙度的弱化修正，得到最终的表面改质特征长度；计算公式如下：

13、；

14、其中，

β为表面改质特征长度，

β

0为理想表面改质特征长度，

σ

res为表层残余压应力，

r

q为表面均方根粗糙度，

k

1为残余应力拟合系数，

k

2为表面粗糙度拟合系数，

σ

s为应力量纲参数，

l

0为长度量纲参数。

15、作为本发明的进一步改进，所述临界表面改质特征长度通过标准试样微动疲劳试验标定得到，标定规则包括：制备多组仅改变表面改质强度、其余参数一致的标准试样；对每组试样进行微动疲劳加载，观测裂纹萌生位置；拟合裂纹萌生位置发生表面/次表面转变时对应的表面改质特征长度，作为所述临界表面改质特征长度。

16、作为本发明的进一步改进，还包括：

17、s9、表面改质工艺反向调优：预设构件目标剩余疲劳寿命，将步骤s7的计算过程逆向迭代，反演求解满足目标寿命所需的表面改质特征长度最优值，并据此输出表面改质层剪切模量或改质厚度的调控范围，作为表面改质工艺优化的定量依据。

18、作为本发明的进一步改进，步骤s1中，现场工况参数包括：构件的接触面法向线集中力、远场交变地应力幅、接触特征长度、实际接触面积及初始裂纹长度，并确定临界失效裂纹长度 ；

19、界面与疲劳参数包括：界面摩擦系数、微动耦合修正系数、几何形状修正系数、微动疲劳常数；

20、基体与表面改质层力学参数包括：构件基体剪切模量、表面改质层剪切模量、表层残余应力、表面均方根粗糙度、残余应力影响系数及表面粗糙度影响系数。

21、作为本发明的进一步改进，步骤s2中，所述表面改质层力学模型为gurtin–murdoch表面弹性模型，所述修正基本解为平面半无限体表面集中力修正基本解；

22、步骤s3中，所述摩擦模型为库仑摩擦模型，所述滑移理论为弹性滑移理论。

23、作为本发明的进一步改进，步骤s5中，所述综合等效应力幅由远场交变地应力幅与经微动耦合修正系数修正后的界面平均切应力叠加得到，计算公式如下：

24、；

25、其中，为综合等效应力幅，α为微动耦合修正系数，δσ为远场交变地应力幅，τ为界面平均切应力。

26、作为本发明的进一步改进，步骤s8中，所述分级预警包括：根据裂纹扩展速率和/或剩余疲劳寿命将构件工作状态划分为安全区、警戒区和危险区；所述差异化处置方案包括：安全区维持常规监测，警戒区加密检测频次并评估表面改质层完整性，危险区立即开展加固、补缝或界面再处理。

27、第二方面，本发明实施例提供了一种考虑表面改质的地下工程构件微动裂纹扩展速率预测系统，其用于执行上述考虑表面改质的地下工程构件微动裂纹扩展速率预测方法，包括：

28、参数采集与标定模块，用于获取经表面改质处理的深部矿井微动易发构件的几何参数、载荷参数、材料力学参数、表面改质参数及微动疲劳参数，并计算表面改质特征长度及标定临界表面改质特征长度；

29、接触应力场求解模块，用于基于表面改质集中力修正基本解计算构件接触面法向应力与极坐标径向应力；

30、界面微动参数计算模块，用于根据接触面法向应力计算界面摩擦力、界面平均切应力及微动滑移幅值；

31、裂纹萌生模式判定模块，用于将表面改质特征长度与临界表面改质特征长度进行比较，判定微动裂纹萌生位置并选择对应的微动疲劳常数；

32、耦合计算模块，用于将远场交变地应力与微动界面切应力耦合，得到综合等效应力幅，并计算裂纹尖端应力强度因子幅；

33、裂纹扩展与寿命计算模块，用于基于paris公式计算裂纹扩展速率，并通过积分获得剩余疲劳寿命；

34、分级预警与输出模块，用于根据裂纹扩展速率和剩余疲劳寿命进行分级预警，输出评估报告；

35、工艺反向调优模块，用于预设目标剩余疲劳寿命，逆向迭代反演求解满足目标寿命所需的表面改质特征长度最优值，并输出表面改质层剪切模量或改质厚度的调控范围。

36、有益效果：

37、本发明相比于现有技术，具备如下技术功效：

38、1、耦合滚压残余压应力、表面粗糙度两类核心工艺/界面参数：分别构建残余应力正向强化修正项、粗糙度负向弱化修正项，建立多物理量耦合的表面改质特征长度计算模型。模型可完全退化为仅考虑剪切模量的传统模型，兼容性强；完整表征界面剪切刚度、滚压残余压应力强化、粗糙微动应力集中共同作用下改性层有效影响范围，显著提升地下构件微动裂纹扩展速率的预测精度。

39、2、修正应力奇异性，计算精度大幅提升：引入gurtin-murdoch表面弹性模型与平面表面改质集中力基本解，解决了传统flamant解应力奇异的问题，精准表征表面改质层的应力均化作用，应力场计算结果与工程实际高度吻合。

40、3、模型耦合性完整，适配改性构件：建立表面改质参数-接触应力-微动滑移-裂纹扩展全链条定量模型，专门针对滚压、涂层、硬化等表面改性地下构件设计，填补了现有技术空白。

41、4、理论与工程结合，应用范围广：可应用于井下预制管片、连接螺栓、支护锚杆、钢拱架、基坑钢支撑、综合管廊拼接构件等各类带表面改质层的地下工程结构，通用性强。

42、5、双重指导价值：不仅可实现微动裂纹扩展速率与剩余寿命的精准预测，还可反向优化表面改质层剪切模量、改质厚度等工艺参数，从源头抑制微动裂纹萌生。

43、6、拓展性强：可直接对接现场传感监测系统，实现参数自动采集、模型在线运算，适配智能化地下工程健康监测平台。

44、上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的具体实施方式。