ZWT本构模型二次开发，计算出现负体积

kb2018

有没有做过的前辈可以解答下，我最近在做ZWT本构模型的二次开发，只是编了
[img]，前面σ= C0 ε+αε2 +βε3 一小部分，一计算就出现负体积，单位制检查过没有问题，希望做过的前辈解答下，不甚感激。 p1~p4:
                                         1173.0(Mpa)        –131120(Mpa)    6.329×106(Mpa)      0.28

1.       subroutine umat41 (cm,eps,sig,epsp,hsv,dt1,capa,
   
2.      . etype,time,temp,failel,crv,cma)
   
3.       dimension  cm(*),eps(*),sig(*),hsv(*)
   
4.       character*(*) etype
   
5.       hsv(1)=hsv(1)+eps(1)
   
6.       hsv(2)=hsv(2)+eps(2)
   
7.       hsv(3)=hsv(3)+eps(3)
   
8.       hsv(4)=hsv(4)+eps(4)
   
9.       hsv(5)=hsv(5)+eps(5)
   
10.       hsv(6)=hsv(6)+eps(6)
    
11.       ga=1./((1.+cm(4))*(1.-2.*cm(4)))
    
12.       hsv(7)=ga*((1.-cm(4))*hsv(1)+cm(4)*(hsv(2)+hsv(3)))
    
13.       hsv(8)=ga*((1.-cm(4))*hsv(2)+cm(4)*(hsv(1)+hsv(3)))
    
14.       hsv(9)=ga*((1.-cm(4))*hsv(3)+cm(4)*(hsv(1)+hsv(2)))
    
15.       hsv(10)=ga*((1.-2.*cm(4))/2.*hsv(4))
    
16.       hsv(11)=ga*((1.-2.*cm(4))/2.*hsv(5))
    
17.       hsv(12)=ga*((1.-2.*cm(4))/2.*hsv(6))
    
18.       sig(1)=sig(1)+(cm(1)+2.*cm(2)*hsv(7)+
    
19.      . 3.*cm(3)*(hsv(7))**2)
    
20.      . *ga*((1.-cm(4))*eps(1)+cm(4)*(eps(2)+eps(3)))           
    
21.       sig(2)=sig(2)+(cm(1)+2.*cm(2)*hsv(8)
    
22.      . +3.*cm(3)*(hsv(8))**2)
    
23.      . *ga*((1.-cm(4))*eps(2)+cm(4)*(eps(1)+eps(3)))
    
24.       sig(3)=sig(3)+(cm(1)+2.*cm(2)*hsv(9)
    
25.      . +3.*cm(3)*(hsv(9))**2)
    
26.      . *ga*((1.-cm(4))*eps(3)+cm(4)*(eps(1)+eps(2)))
    
27.       sig(4)=sig(4)+(cm(1)+2.*cm(2)*hsv(10)
    
28.      . +3.*cm(3)*(hsv(10))**2)
    
29.      . *ga*((1.-2.*cm(4))/2*eps(4))
    
30.       sig(5)=sig(5)+(cm(1)+2.*cm(2)*hsv(11)
    
31.      . +3.*cm(3)*(hsv(11))**2)
    
32.      . *ga*((1.-2.*cm(4))/2*eps(5))
    
33.       sig(6)=sig(6)+(cm(1)+2.*cm(2)*hsv(12)
    
34.      . +3.*cm(3)*(hsv(12))**2)
    
35.      . *ga*((1.-2.*cm(4))/2*eps(6))            
    
36.       return
    
37.       end

复制代码

nicky_wang22

请问一下您是怎么通过license认证的呀