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 doesnt 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
Sorry forgot the variables initialization lines in my post.
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