×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!
  • Students Click Here

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Jobs

DFLUX subroutine

DFLUX subroutine

DFLUX subroutine

(OP)
I am using DFLUX for concrete heat generation. Flux(1) is depending on eqvage. My code is below. This subroutine is working perfectly when there is only one concrete part. If I add some insulation, It will not work right. The eqvage is not adding up correctly for some reason if there is more than one part? I cant figure out what is the difference. Heat source is never changed.


subroutine dflux(flux,sol,kstep,kinc,time,noel,npt,coords,
1 jltyp,temp,press,sname)
c
include 'aba_param.inc'
c
dimension flux(2), time(2), coords(3)
character*80 sname
c
common xys(100000)
xys(time(1)) = exp(4070.3D0*(1/296D0-1/(273D0+(sol-32D0)/1.8D0)))
eqvage = 0
c
c
do 100 n = 1,time(1),1
100 eqvage = eqvage+xys(n)
c
if (eqvage<=7.61) then
flux(1) = 70D0 + 15.11D0*eqvage
else if (eqvage<=17) then
flux(1) = 185D0 - 11.82D0*(eqvage-7.61)
else if (eqvage<=100) then
flux(1) = 74D0 - 0.651D0*(eqvage - 17)
else if (eqvage<=300) then
flux(1) = 20D0 - 0.1D0*(eqvage - 100)
else
flux(1) = 0
end if

c
write(7,*) '!!!',time(1),xys(time(1)), eqvage,sol
return
end

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members!


Resources