Defining Temperature as a moving boundary condition in ABAQUS.

[Haider1040]

Dear all,
I am trying to define temperature as a moving boundary condition. I want to apply a temperature to a small circular area on a surface and then move that around on the surface. I was told that DISP subroutine. I have put together this subroutine but it doesn't work.

SUBROUTINE DISP(U,KSTEP,KINC,TIME,NODE,NOEL,JDOF,COORDS)
c
include 'ABA_PARAM.INC'
c
DIMENSION U(3),TIME(2),COORDS(3)
c Defining radius of the circular region for temperature application
x=COORDS(1)
y=COORDS(2)
r=sqrt((x-x0)**2+(y-y0)**2)
c
c Test of node position and temperature assignment
if (JDOF.eq.3) then
If (r.le.rlas) then
U(1)=1800
Else
U(1)=0
end if
U(2)=0
U(3)=0
End if
c
return
end

Can somebody please help me if I am thinking correctly or what am I doing wrong here.

Thank you
Haider

 

[Mustaine3]

Why UDISP instead of UTEMP?

Where are x0, y0 and rlas defined?

 

[Haider1040]

Dear Mustaine3,
Thank you for replying. Sorry I had forgotten to include the lines of x0,y0 and rlas in my post. Below is the complete subroutine. As for UTEMP or DISP i have no idea which one to use. Just started with DISP. Can you please help me. Cheers

SUBROUTINE DISP(U,KSTEP,KINC,TIME,NODE,NOEL,JDOF,COORDS)
c
include 'ABA_PARAM.INC'
c
DIMENSION U(3),TIME(2),COORDS(3)
real*8 x0,y0,x,y,z,rlas,r
x0=0
y0=0
rlas=0.05
c Defining radius of the circular region for temperature application
x=COORDS(1)
y=COORDS(2)
r=sqrt((x-x0)**2+(y-y0)**2)
c
c Test of node position and temperature assignment
if (JDOF.eq.3) then
If (r.le.rlas) then
U(1)=1800
Else
U(1)=0
end if
U(2)=0
U(3)=0
End if
c
return
end

 

[Mustaine3]

Do you have a model with translational DOF? Because that is what you apply with a DISP.