如何把单个单元的等效应力、等效应变输出到指定文件...

liyunliang_zi

各位高手：
      仿效书上的例子，想把编号为1904的单元状态输出出来：

      iext=nelmntid(i,0)
      if(iext.eq.1904) then
          open(1,file="out.txt")
          write(1,*) 等效应力代号，等效应变代号
      endif
可是，调试通不过。
在dyna970中，i代表单元内部编号，在subroutine umat41(cm,eps,sig,hisv,dt1,capa,etype,time, temp,i,ixs,x,k,j)传参中出现了i。

各位高手，你们怎么实现的？

liyunliang_zi

找不到好办法，我用以下方式调试：
hsv(3)=hsv(3)+1
if (hsv.gt.100) then
   write(*,*) scle, aj1, ak fac1, depi, hsv(1)
   hsv(3)=0
endif

但是，该方法的缺点会造成计算特别慢，而且，不知道输出的是哪个单元的信息。