Defining Temperature as a moving boundary condition in ABAQUS.
Defining Temperature as a moving boundary condition in ABAQUS.
(OP)
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
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





RE: Defining Temperature as a moving boundary condition in ABAQUS.
Where are x0, y0 and rlas defined?
RE: Defining Temperature as a moving boundary condition in ABAQUS.
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
RE: Defining Temperature as a moving boundary condition in ABAQUS.