Search
Defining temperatur...

Clear all

# Defining temperature as a moving boundary condition in abaqus?

2 Posts
1 Users
0 Reactions
902 Views
Posts: 2
Topic starter
(@golbed121)
New Member
Joined: 8 years ago

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

Topic Tags
Posts: 2
Topic starter
(@golbed121)
New Member
Joined: 8 years ago

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

Topic Tags