I am doing welding analysis with phase transformation effects in Abaqus. For the developed UEXPAN subroutine, I have to make the system understand whether heating/cooling is taking place. To flag the elements, I used the USDFLD and GETVRM. But I have to calculate the temp difference (TEMP(1)-Temp_old) so that if the difference is negative, I can mark it as cooling, and if positive, I can mark it as Heating. It will help to give the EXPAN(1) at the right time. How can we get the old temperature for the same?
https://abaqus-docs.mit.edu/2017/English/SIMACAECAERefMap/simacae-t-plgexamplevolavetemp.htm
Himadri Nath Saha
Thank you for the response .
But I need to use the temperature values to mark the elements as Heating or cooling. This information has to be passed to the UEXPAN subroutine so that it can use the stain increment at the right time.
I tried to do a heat transfer analysis alone and checked the temperature(assigned as state variable 1) at the output .but its showing zero value.
SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,
1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,
2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO,LACCFLA)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME,ORNAME
CHARACTER*3 FLGRAY(15)
DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),
1 T(3,3),TIME(2)
DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*),COORD(*)
CALL GETVRM('TEMP',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,
1 MATLAYO,LACCFLA)
TEMP=ARRAY(1)
STATEV(1) = TEMP
RETURN
END
You have used Body heat flux, can u please provide me the subroutine .for file
- Badges
- Science topic