Re: Movement Heat Flux using DFLUX subroutine

http://imechanica.org/node/18499



Am 18.01.2019 um 18:05 schrieb Rami Ghorbel [email protected] [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]
>
>