当前位置：首页 > news >正文

GTN损伤模型在金属成型中的实战改造手记

news 2026/6/17 8:19:45

修正的考虑剪切和静水压力的各向异性和各向同性 GTN 模型 Abaqus_Vumat 子程序以及跑通的案例，包括三个子程序 for 文件（各向同性，自己推导修改以及各向异性）、cae 文件和 inp 文件等，还有一些论文资料。注：1. 各项异性模型为 hill48（三维应力状态） 2. abaqus 版本为 2020，更低版本可以导入 inp 文件查看及求解 3. 子程序是在一个各向同性子程序基础上改的，推导了相关公式，修改了感觉不太对的地方子程序都会提供，可以自己推导。 4.基质屈服应力考虑了温度和应变率，另外程序中提供了两种常用的本构。

搞金属成型仿真的老司机应该都踩过GTN模型的坑——这玩意儿在剪切和静水压力下的预测总是飘忽不定。最近为了项目需求，硬刚了半个月的VUMAT子程序，总算把各向同性和各向异性的修正版GTN模型都跑通了。这里分享些实战心得，附上三个版本的子程序文件（各向同性原版、各向同性修正版、Hill48各向异性版）和验证案例。

一、屈服函数的暴力改造

原版各向同性模型最坑的地方在于q2参数的处理。在静水压力主导的工况下，屈服函数中的指数项会直接崩掉。我们团队推导时加入了修正系数η：

! 原版屈服函数 Phi = (q/sigmabar)**2 + 2.0*q1*fstar*cosh(1.5*q2*p/sigmabar) - (1.0 + q3*fstar**2) ! 修正后加入压力修正项 eta = 1.0/(1.0 + alpha*(p/sigmabar)**beta) Phi = (q/sigmabar)**2 + 2.0*q1*fstar*cosh(1.5*q2*eta*p/sigmabar) - (1.0 + q3*fstar**2)

这个eta函数在高压状态下能有效抑制cosh项爆炸（参数α和β需要通过单轴压缩试验标定）。调试时发现当β取0.5时，304不锈钢的预测误差从23%降到了7%以内。

二、Hill48各向异性的嫁接手术

各向异性改造的核心是把等效应力替换成Hill48准则。这里有个巨坑：Hill的应力分量顺序和Abaqus内置顺序不一致！

! Hill48等效应力计算（注意应力分量顺序！） F = 0.5*(HillF*(S(2)-S(3))**2 + HillG*(S(3)-S(1))**2 + HillH*(S(1)-S(2))**2) F = F + HillL*S(4)**2 + HillM*S(5)**2 + HillN*S(6)**2 seq = sqrt(F)

这里S(1)~S(6)对应Abaqus的应力分量顺序是σ11,σ22,σ33,σ12,σ13,σ23。当初没注意这个，直接套用论文公式导致预测结果完全反常识，浪费了两天时间查错...

三、温度应变率的双杀处理

基质材料的屈服应力采用修正Johnson-Cook模型，支持两种模式切换：

! 模式1: 标准J-C模型 IF(MODEL_TYPE == 1) THEN sigmay = (A + B*epbar**n) * (1.0 + C*log(epdot/epdot0)) * (1.0 - (T-Troom)/(Tmelt-Troom)**m) ! 模式2: 考虑动态回复的修正版 ELSE sigmay = A*(1.0 + epbar/ep0)**B * exp(-C*(T-Troom)) * (1.0 + (epdot/epdot0)**D) ENDIF

在高温高速工况下（比如热冲压），模式2比传统模型精度提升约15%。但要注意D参数不能为负，否则会在低应变率时出现应力骤降的异常现象。