JSRUN 用代码说话
精选fortran
全部fortran
面热源子程序 SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,JLTYP, 1 TEMP,PRESS,SNAME) C INCLUDE 'ABA_PARAM.INC' C DIMENSION COORDS(3),FLUX(2),TIME(2) CHARACTER*80 SNAME C U=16.5 AI=60. v=250./60 yita=0.75 R0=5. qm=yita*U*AI*1000./3.14/R0/R0 dx=v*time(1) rr=((COORDS(1)-dx)**2+COORDS(2)**2 FLUX(1)=3*qm*EXP(-(3*rr)/R0**2) RETURN END