marecare02
Mechanical
- Mar 19, 2010
- 6
Hi!
I want to know if it's possible to write subroutine of 3 different time dependent surface film conditions at same time, with differnet h1, h2 and h3? they aren't temperature dependent. I wrote subroutine for one, but i am no sure if i can write for 3 sfilms. (i want to simulate 3 different intakes of energy, simultaneous, independent).
can you help me please. here is my code for one...
SUBROUTINE FILM(H,SINK,TEMP,KSTEP,KINC,TIME,NOEL,NPT,
1 COORDS,JLTYP,FIELD,NFIELD,SNAME,NODE,AREA)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION H(2),TIME(2),COORDS(3), FIELD(NFIELD)
CHARACTER*80 SNAME
C
SINK=1000
C
x=coords(1)
y=coords(2)
C
t2=time(2)
C
C
b=4
a=3
C
H(1)=100*exp(-(((y-t2)/a)**4+((x)/b)**4))/4*a*b
C
return
end
I want to know if it's possible to write subroutine of 3 different time dependent surface film conditions at same time, with differnet h1, h2 and h3? they aren't temperature dependent. I wrote subroutine for one, but i am no sure if i can write for 3 sfilms. (i want to simulate 3 different intakes of energy, simultaneous, independent).
can you help me please. here is my code for one...
SUBROUTINE FILM(H,SINK,TEMP,KSTEP,KINC,TIME,NOEL,NPT,
1 COORDS,JLTYP,FIELD,NFIELD,SNAME,NODE,AREA)
C
INCLUDE 'ABA_PARAM.INC'
C
DIMENSION H(2),TIME(2),COORDS(3), FIELD(NFIELD)
CHARACTER*80 SNAME
C
SINK=1000
C
x=coords(1)
y=coords(2)
C
t2=time(2)
C
C
b=4
a=3
C
H(1)=100*exp(-(((y-t2)/a)**4+((x)/b)**4))/4*a*b
C
return
end