Movement Heat Flux using DFLUX subroutine

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

Movement Heat Flux using DFLUX subroutine

Rami Ghorbel
Dear All,

I am currently working on Heat movement using DFLUX subroutine.
I wrote a subroutine using double ellipsoidal Gaussian distribution.
Unfortunately, results are false.

Attached my model and my subroutine.
Anyone can tell me what is the possible error?

Many Thanks;
Rami


--
*Cordialement;*

*Rami GHORBEL*

*Material Engineer*

*Phd Student*

*Phone :* (00216) 22 667 771

  ----------

      SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP,
     1                 TEMP,PRESS,SNAME)
C
      INCLUDE 'ABA_PARAM.INC'

      DIMENSION COORDS(3),FLUX(2),TIME(2)
      CHARACTER*80 SNAME
      Real t, x, y, z, x1, y1, z1, Lweld, v
      Real Q, a, b, c, pi, x0, y0, z0, heat, shape1
C time
      t=TIME(2)
      a=0.3
      b=0.5
      c=0.2
      cf=0.003
        pi=3.141592
       
      FLUX(1)=0
      FLUX(2)=0
      Lweld=60
      v=5
        x=COORDS(1)
        y=COORDS(2)
        z=COORDS(3)
c         try
      if (JSTEP. eq. 1) then
     
        x0=0
      y0=0
      z0=2
      Q=0.5*10*900
      x1=x-x0
        y1=y-v*t
        z1=z-z0

      heat=6*sqrt(3.0)*Q/(a*b*c*pi*sqrt(pi))
      shape1=exp(-3*(x1)**2/a**2-3*(y1)**2/b**2-3*(z1)**2/c**2)
     
      FLUX(1)=heat*shape1
     
      end if
     
      RETURN
      END


[Non-text portions of this message have been removed]

Reply | Threaded
Open this post in threaded view
|

Re: Movement Heat Flux using DFLUX subroutine

Frank Richter-2

http://imechanica.org/node/18499



Am 18.01.2019 um 18:05 schrieb Rami Ghorbel [hidden email] [Abaqus]:

>
> Dear All,
>
> I am currently working on Heat movement using DFLUX subroutine.
> I wrote a subroutine using double ellipsoidal Gaussian distribution.
> Unfortunately, results are false.
>
> Attached my model and my subroutine.
> Anyone can tell me what is the possible error?
>
> Many Thanks;
> Rami
>
> --
> *Cordialement;*
>
> *Rami GHORBEL*
>
> *Material Engineer*
>
> *Phd Student*
>
> *Phone :* (00216) 22 667 771
>
> ----------
>
> SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP,
> 1 TEMP,PRESS,SNAME)
> C
> INCLUDE 'ABA_PARAM.INC'
>
> DIMENSION COORDS(3),FLUX(2),TIME(2)
> CHARACTER*80 SNAME
> Real t, x, y, z, x1, y1, z1, Lweld, v
> Real Q, a, b, c, pi, x0, y0, z0, heat, shape1
> C time
> t=TIME(2)
> a=0.3
> b=0.5
> c=0.2
> cf=0.003
> pi=3.141592
>
> FLUX(1)=0
> FLUX(2)=0
> Lweld=60
> v=5
> x=COORDS(1)
> y=COORDS(2)
> z=COORDS(3)
> c try
> if (JSTEP. eq. 1) then
>
> x0=0
> y0=0
> z0=2
> Q=0.5*10*900
> x1=x-x0
> y1=y-v*t
> z1=z-z0
>
> heat=6*sqrt(3.0)*Q/(a*b*c*pi*sqrt(pi))
> shape1=exp(-3*(x1)**2/a**2-3*(y1)**2/b**2-3*(z1)**2/c**2)
>
> FLUX(1)=heat*shape1
>
> end if
>
> RETURN
> END
>
> [Non-text portions of this message have been removed]
>
>