- 积分
- 0
- 注册时间
- 2024-6-2
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 sfl769 于 2024-6-3 21:48 编辑
各位好,请问如何在热力耦合过程中,应力场运算中使用usdfld子程序实现材料温度大于一定值时单元删除,我编写了一份子程序,使用2个非独立变量,控制单元删除变量为第二个,虽然能输出结果云图,但并没有显示单元删除,请问是子程序的问题还是别的设置的问题,万分感谢,子程序如下:
- subroutine usdfld(field,statev,pnewdt,direct,t,celent,time,dtime,
- 1 cmname,orname,nfield,nstatv,noel,npt,layer,kspt,kstep,kinc,
- 2 ndi,nshr,coord,jmac,jmtyp,matlayo,laccflg)
- c
- include 'aba_param.inc'
- c
- character*80 cmname,orname
- character*3 flgray(15)
- dimension field(nfield),statev(nstatv),direct(3,3),t(3,3),time(2),
- * coord(*),jmac(*),jmtyp(*)
- dimension array(15),jarray(15)
- c
- c Get temperatures from previous increment
- call getvrm('TEMP',array,jarray,flgray,jrcd,
- $ jmac, jmtyp, matlayo, laccflg)
-
- field(1)=array(1)
- if (field(1).GE.1530) then
- statev(2)=0.0
- else
- statev(2)=1.0
- end if
-
- return
- end
复制代码 |
|