Hi everyone I am performing a simulation of casting where I am considering the Heat transfer coefficient as a function of time.

For which I am using the GAPCON subroutine.

Please can anyone go through my subroutine and tell if there is any issue in it ?

subroutine gapcon(ak,d,flowm,temp,predef,time,ciname,
1 slname,msname,coords,noel,node,npred,kstep,kinc)
c
include 'aba_param.inc'
c
character*80 ciname,slname,msname
dimension flowm(2),temp(2),predef(2,*),
1 ak(5),d(2),coords(3),time(2)
c
c define gap conductivity dependent on temperature.
 
real*8 :: GC_slope, GC_intercept, GAP_slope, GAP_intercept
 
if (temp(1) <= 550.00) then
ak(1) = 2700.00
elseif (temp(1) >= 550.00 .and. temp(1) <= 610.00) then
ak(1) = 2.73 * temp(1) + 1034.7
elseif(temp(1) > 610.00 .and. temp(1) <= 720.00) then
ak(1) = 3000.00
end if
 

What is happening is my simulation is not matching with my exact process. If I am increasing the step time also then my acting is not going below the solidification temperature. Please can anyone help me with this?

Thanks
 
ak(2) =0.0
ak(3) = 0.0
ak(4) = 2.73
ak(5) = 0.0
return
end

Abaqus ​​​​​​​