CN115691712B 基于累积剪切应变损伤变量的微观疲劳裂纹扩展模拟方法 （上海理工大学）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN115691712B(45)授权公告日2025.07.04

(21)申请号202211279471.1

(22)申请日2022.10.19

(65)同一申请的已公布的文献号申请公布号CN115691712A

(43)申请公布日2023.02.03

(73)专利权人上海理工大学

地址200082上海市杨浦区军工路516号

专利权人华东理工大学

(72)发明人杨杰郭昊函张显程杨荟楠王润梓卢荣生

GO6F30/23(2020.01)

GO6F30/25(2020.01)GO6F111/10(2020.01)

GO6F119/02(2020.01)

GO6F119/14(2020.01)

(56)对比文件

CN114662356A,2022.06.24CN115062378A,2022.09.16审查员钟嘉欣

(74)专利代理机构北京众达德权知识产权代理

有限公司11570专利代理师邢升

(51)Int.CI.

G16C60/00(2019.01)权利要求书2页说明书6页附图2页

(54)发明名称

基于累积剪切应变损伤变量的微观疲劳裂纹扩展模拟方法

(57)摘要

该发明涉及疲劳裂纹扩展数值模拟技术领域，尤其涉及一种基于累积剪切应变损伤变量的微观疲劳裂纹扩展模拟方法，包括：生成Voronoi图；完成代表性体积单元的建立；编写UDMGINI损伤子程序并将其嵌入UMAT子程序；在有限元计算时嵌入UMAT子程序，得到晶体塑性结果参量；UDMGINI子程序对晶体塑性结果参量进行调用，计算各滑移系的累积剪切应变，选择累积剪切应变最大的滑移系为最活跃滑移系；当最活跃滑移系上累积剪切应变达到临界值时，判定裂纹扩展。该方法实现了多晶材料在微观尺度下的疲劳

CN

CN115691712B

CN115691712B权利要求书1/2页

2

1.基于累积剪切应变损伤变量的微观疲劳裂纹扩展模拟方法，其特征在于，包括以下步骤：

S1、通过Matlab语言生成多晶材料的Voronoi图和对应的随机晶粒取向；

S2、在ABAQUS中通过布尔运算对多晶材料建立代表性体积单元，通过Python语言读取S1中Voronoi图的拓扑信息和每个晶粒的取向信息，随同材料参数批量赋予到每个晶粒中；

S3、编写用户自定义损伤子程序UDMGINI,将其嵌入至用户自定义材料子程序UMAT中；

S4、将UMAT子程序嵌入ABAQUS软件，进行有限元计算，生成每一增量步对应的应力增量和塑性应变增量，并基于此对当前时刻的应力应变状态值进行更新；

S5、UMAT子程序将应力应变状态、增量信息和材料参数作为输入变量，生成并旋转得到全局弹性矩阵，进一步迭代求解得到晶体塑性结果参量，并将其储存至雅可比矩阵DDSDDE和可共享的状态变量组STATEV中，完成材料响应的生成；

S6、UDMGINI子程序对S5中包含结果参量的STATEV组进行调用，并对裂纹尖端前方单元的累积剪切应变及最活跃滑移系进行计算与判定，累积剪切应变为：

式中：γ为累积剪切应变，是塑性剪切应变率“基于一定时间间隔的积分；

S7、若最活跃滑移系a上的累积剪切应变γa与临界总剪切应变γcr的比值小于1,即,则返回S4继续进行有限元计算；若累积剪切应变Ya与临界总剪切应变γcr的比值等于1,即裂纹达到开裂标准，执行裂纹扩展命令，扩展方向平行于当前最活跃滑移系的滑动方向；

S8、UDMGINI子程序将S7中的扩展信息反馈至ABAQUS求解器中；

S9、ABAQUS求解器根据反馈信息对各个单元的刚度矩阵进行组装；

S10、对已组装的刚度矩阵进行有限元求解；

S11、对求解结果的收敛性进行判定：若结果收敛，则根据单元平均解判断并更新裂纹集合，成功实现裂纹扩展，并将结果储存以进行下一个增量步的迭代运算；若不收敛，则自动采用较小增量步重复以上步骤直至收敛或达到最小增量步长为止。

2.如权利要求1所述的基于累积剪切应变损伤变量的微观疲劳裂纹扩展模拟方法，其特征在于：

所述S3中，UMAT子程序框架基于晶体塑性理论，采用基于滑移的幂函数型晶体塑性本构方程，硬化