Subroutine for coefficient of thermal expansion

[anoop_1989]

Dear all,
I am trying to run a coefficient of thermal expansion subroutine using the UEXP subroutine. For this, like a first step I tried a simple thermal stress analysis on a 0.2*0.1steel bar giving a body heat flux. The calculated stress values are matching with the theoretical analysis. Now instead of giving the thermal expansion as a direct input, I would like to enter the same as a subroutine. The attached subroutine is given below. But the problem aborted due to compilation error. Can anyone help to solve the issue? For the next stage, I have to incorporate the change of volumetric strain and the thermal strain to analyze the phase transformation in welding.

subroutine uexpan(expan,dexpandt,temp,time,dtime,predef,dpred,
statev,cmname,nstatv,noel)
c
include 'aba_param.inc'
c
character*80 cmname
c
dimension expan(*),dexpandt(*),temp(2),time(2),predef(*),
dpred(*),statev(nstatv)
c
alpha = 1.1d-05
c
expan(1) = alpha*temp(2)
c
return
end

 

[anoop_1989]

cae file

 

[anoop_1989]

jnl file

 

[FEA way]

Do you have proper version of Fortran compiler installed ? Try the following command to see if subroutines work correctly on your machine: abaqus verify -user_std

You can also use exemplary UEXPAN subroutine available in the Verification Guide (User subroutine verification —> UEXPAN).

 

[Mustaine3]

You didn't attach the subroutine and copying the unformatted code from your post won't help to answer your question.

The first example in the Verification Manual is basically what you are trying to do. Why don't you take that file and start from there?

 

[anoop_1989]

Dear FEA way and Mustaine3

THANK YOU FOR THE COMMENTS

Fortran is working fine. I have done welding analysis previously using the DFLUX subroutine.
I have used the same subroutine from the Abaqus manual to verify the code. It is working well with the Abaqus inp file. But showing an error with my inp file. So I changed a few formats in the file from the Abaqus manual and it is working. The corrected file is attached

Now I have to do include the phase transformation change due to volumetric change along with the thermal expansion(epsilon theral+epsilon volumetric). Does anyone have a standard file for the same?