www.pudn.com > tempD.rar > tempD.for, change:2016-11-11,size:5974b


 	SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP, 
	1                  TEMP,PRESS,SNAME) 
C 
	INCLUDE 'ABA_PARAM.INC' 
C 
C 
	DIMENSION COORDS(3),FLUX(2),TIME(2) 
	CHARACTER*80 SNAME 
C 
      SELECT CASE(JSTEP) 
C       !A00-1 
      CASE(3) 
      wu=22.0 
	wi=150.0 
	x0=0 
	y0=1.9 
	aa=8 
	a=20 
	b=6 
	c=4 
	v=3 
C      !A01-1 
	CASE(6) 
      wu=32.0 
	wi=500.0 
	x0=-2.0 
	y0=5 
	aa=8 
	a=20 
	b=15 
	c=6 
	v=7.5 
C  !A01-2 
	CASE(9) 
      wu=32.0 
	wi=500.0 
	x0=2.0 
	y0=5 
	aa=8 
	a=20 
	b=15 
	c=6 
	v=7.5 
C     !A02-1 
      CASE(12) 
      wu=32.0 
	wi=600.0 
	x0=-3.5 
	y0=8.1 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	   !A02-2 
      CASE(15) 
      wu=32.0 
	wi=600.0 
	x0=3.5 
	y0=8.1 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	    !B00-1 
      CASE(18) 
      wu=22.0 
	wi=180.0 
	x0=0 
	y0=-1.9 
	aa=8 
	a=20 
	b=6 
	c=4 
	v=3 
C	   !B01-1 
      CASE(21) 
      wu=32.0 
	wi=500.0 
	x0=2.0 
	y0=-5 
	aa=8 
	a=20 
	b=15 
	c=6 
	v=7.5 
C	   !B01-2 
      CASE(24) 
      wu=32.0 
	wi=500.0 
	x0=-2.0 
	y0=-5 
	aa=8 
	a=20 
	b=15 
	c=6 
	v=7.5 
C	  !B02-1 
      CASE(27) 
      wu=32.0 
	wi=600.0 
	x0=3.5 
	y0=-8.1 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C    !B02-2 
      CASE(30) 
      wu=32.0 
	wi=600.0 
	x0=-3.5 
	y0=-8.1 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C     !B03-1 
	CASE(33) 
      wu=32.0 
	wi=600.0 
	x0=4.0 
	y0=-11.2 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C    !B03-2 
	CASE(36) 
      wu=32.0 
	wi=600.0 
	x0=-4.0 
	y0=-11.2 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C      !B04-1 
      CASE(39) 
      wu=32.0 
	wi=600.0 
	x0=4.2 
	y0=-14.3 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	    !B04-2 
      CASE(42) 
      wu=32.0 
	wi=600.0 
	x0=-4.2 
	y0=-14.3 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	    !A03-1 
      CASE(45) 
      wu=32.0 
	wi=600.0 
	x0=-4.0 
	y0=11.2 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	    !A03-2 
      CASE(48) 
      wu=32.0 
	wi=600.0 
	x0=4.0 
	y0=11.2 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	    !A04-1 
      CASE(51) 
      wu=32.0 
	wi=600.0 
	x0=-4.2 
	y0=14.3 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	    !A04-2 
      CASE(54) 
      wu=32.0 
	wi=600.0 
	x0=4.2 
	y0=14.3 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C     !A05-1 
      CASE(57) 
      wu=32.0 
	wi=600.0 
	x0=-4.4 
	y0=17.4 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C     !A05-2 
	CASE(60) 
      wu=32.0 
	wi=600.0 
	x0=4.4 
	y0=17.4 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C      !B05-1 
	CASE(63) 
      wu=32.0 
	wi=600.0 
	x0=4.4 
	y0=-17.4 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C      !B05-2 
      CASE(66) 
      wu=32.0 
	wi=600.0 
	x0=-4.4 
	y0=-17.4 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	     !B06-1 
      CASE(69) 
      wu=32.0 
	wi=600.0 
	x0=4.8 
	y0=-20.5 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	     !B06-2 
      CASE(72) 
      wu=32.0 
	wi=600.0 
	x0=-4.8 
	y0=-20.5 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	     !B07-1 
      CASE(75) 
      wu=32.0 
	wi=600.0 
	x0=5.1 
	y0=-23.6 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	      !B07-2 
      CASE(78) 
      wu=32.0 
	wi=600.0 
	x0=-5.1 
	y0=-23.6 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	     !B08-1 
      CASE(81) 
      wu=32.0 
	wi=600.0 
	x0=5.4 
	y0=-26.7 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C       !B08-2 
      CASE(84) 
      wu=32.0 
	wi=600.0 
	x0=-5.4 
	y0=-26.7 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C     !B09-1 
	CASE(87) 
      wu=32.0 
	wi=600.0 
	x0=5.7 
	y0=-29.8 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C     !B09-2 
	CASE(90) 
      wu=32.0 
	wi=600.0 
	x0=-5.7 
	y0=-29.8 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C      !B10-1 
      CASE(93) 
      wu=32.0 
	wi=650.0 
	x0=6.0 
	y0=-32.9 
	aa=8 
	a=20 
	b=18 
	c=8 
	v=7.5 
C	     !B10-2 
      CASE(96) 
      wu=32.0 
	wi=650.0 
	x0=-6.0 
	y0=-32.9 
	aa=8 
	a=20 
	b=18 
	c=8 
	v=7.5 
C	    !B11-1 
      CASE(99) 
      wu=32.0 
	wi=700.0 
	x0=5.8 
	y0=-36 
	aa=8 
	a=20 
	b=25 
	c=5 
      v=7.5 
C	    !B11-2 
      CASE(102) 
      wu=32.0 
	wi=700.0 
	x0=-5.8 
	y0=-36 
	aa=8 
	a=20 
	b=25 
	c=5 
	v=7.5 
C	     !A06-1 
      CASE(105) 
      wu=32.0 
	wi=600.0 
	x0=-4.8 
	y0=20.5 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	    !A06-2 
      CASE(108) 
      wu=32.0 
	wi=600.0 
	x0=4.8 
	y0=20.5 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
 
C	     !A07-1 
      CASE(111) 
      wu=32.0 
	wi=600.0 
	x0=-5.1 
	y0=23.6 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	    !A07-2 
      CASE(114) 
      wu=32.0 
	wi=600.0 
	x0=5.1 
	y0=23.6 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C	   !A08-1 
      CASE(117) 
      wu=32.0 
	wi=600.0 
	x0=-5.4 
	y0=26.7 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C     !A08-2 
      CASE(120) 
      wu=32.0 
	wi=600.0 
	x0=5.4 
	y0=26.7 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C     !A09-1 
	CASE(123) 
      wu=32.0 
	wi=600.0 
	x0=-5.7 
	y0=29.8 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C    !A09-2 
	CASE(126) 
      wu=32.0 
	wi=600.0 
	x0=5.7 
	y0=29.8 
	aa=8 
	a=20 
	b=15 
	c=8 
	v=7.5 
C     !A10-1 
      CASE(129) 
      wu=32.0 
	wi=650.0 
	x0=-6.0 
	y0=32.9 
	aa=8 
	a=20 
	b=18 
	c=8 
	v=7.5 
C	    !A10-2 
      CASE(132) 
      wu=32.0 
	wi=650.0 
	x0=6.0 
	y0=32.9 
	aa=8 
	a=20 
	b=18 
	c=8 
	v=7.5 
C	    !A11-1 
      CASE(135) 
      wu=32.0 
	wi=700.0 
	x0=-5.8 
	y0=36 
	aa=8 
	a=20 
	b=25 
	c=5 
	v=7.5 
 
C	    !A11-2 
      CASE(138) 
      wu=32.0 
	wi=700.0 
	x0=5.8 
	y0=36 
	aa=8 
	a=20 
	b=25 
	c=5 
	v=7.5 
c	 
	case default  
	wi=0 
	wu=0 
	x0=0 
	y0=0 
	aa=8 
	a=20 
	b=6.5 
	c=6.2 
	v=0 
C 
	END SELECT 
 
C 
      x=COORDS(1) 
      y=COORDS(2) 
      z=aa-v*time(1) 
      eff1=0.7 
C 
 
C 
      q1=wu*wi*eff1*1000 
C 
      f1=0.6 
C 
      PI=3.1415926 
C 
      heat1=6.0*sqrt(3.0)*q1*f1/(aa*b*c*PI*sqrt(PI)) 
      heat2=6.0*sqrt(3.0)*q1*(2.0-f1)/(a*b*c*PI*sqrt(PI)) 
C 
 
C 
        shape1=exp(-3.0*(x-x0)**2/b**2-3.0*(y-y0)**2/c**2 
     $	-3.0*z**2/aa**2) 
        shape2=exp(-3.0*(x-x0)**2/b**2-3.0*(y-y0)**2/c**2 
     $	-3.0*z**2/a**2) 
C 
C     JLTYP 
      JLTYP=1 
 
	  IF(z .GE. 0) THEN 
	  FLUX(1)=heat1*shape1 
	  ELSE 
	  FLUX(1)=heat2*shape2 
	  ENDIF 
c 
   	  RETURN 
      END