usdfld子程序实现条件单元删除

sfl769

各位好，请问如何在热力耦合过程中，应力场运算中使用usdfld子程序实现材料温度大于一定值时单元删除，我编写了一份子程序，使用2个非独立变量，控制单元删除变量为第二个，虽然能输出结果云图，但并没有显示单元删除，请问是子程序的问题还是别的设置的问题，万分感谢，子程序如下：

1.       subroutine usdfld(field,statev,pnewdt,direct,t,celent,time,dtime,
   
2.      1 cmname,orname,nfield,nstatv,noel,npt,layer,kspt,kstep,kinc,
   
3.      2 ndi,nshr,coord,jmac,jmtyp,matlayo,laccflg)
   
4. c
   
5.       include 'aba_param.inc'
   
6. c
   
7.       character*80 cmname,orname
   
8.       character*3  flgray(15)
   
9.       dimension field(nfield),statev(nstatv),direct(3,3),t(3,3),time(2),
   
10.      * coord(*),jmac(*),jmtyp(*)
    
11.       dimension array(15),jarray(15)
    
12. c
    
13. c Get temperatures from previous increment
    
14.       call getvrm('TEMP',array,jarray,flgray,jrcd,
    
15.      $     jmac, jmtyp, matlayo, laccflg)
    
16.          
    
17.           field(1)=array(1)
    
18.           if (field(1).GE.1530) then
    
19.             statev(2)=0.0
    
20.           else
    
21.             statev(2)=1.0
    
22.           end if
    
23.          
    
24.       return
    
25.       end

复制代码