TIG焊接温度场模拟

应该是你的子程序写错了，建议再核实

老师，这个用你给的也是这样，上面是资料里面实例给的，都是热源中心在4个角。。。。

那就是你模型的坐标与你的程序没对应上

SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP, $ TEMP,PRESS,SNAME) INCLUDE 'ABA_PARAM.INC' DIMENSION COORDS(3),FLUX(2),TIME(2) CHARACTER*80 SNAME C wu,焊接电压 C wi,焊接电流 C effi,焊接效率系数 C q,电弧有效热功率W C v,焊接速度m/s C q,电弧有效热功率W wu=16 wi=80 effi=0.5 v=0.005 q=wu*wi*effi d=v*TIME(1) x=COORDS(1) y=COORDS(2) z=COORDS(3) C 从坐标x0,y0,z0开始，沿着z方向移动 x0=-0.05 y0=0 z0=0.003 C a1,a2,b,c为双椭球的形状参数 a1=0.002 a2=0.003 b=0.002 c=0.003 C f1为热源分配系数 f1=1.0 PI=3.1415926 heat1=6.0*sqrt(3.0)*q/(a1*b*c*PI*sqrt(PI))*f1 heat2=6.0*sqrt(3.0)*q/(a2*b*c*PI*sqrt(PI))*(2.0-f1) shape1=exp(-3.0*(x-x0-d)**2/(a1)**2-3.0*(y-y0)**2/b**2 $ -3.0*(z-z0)**2/c**2) shape2=exp(-3.0*(x-x0-d)**2/(a2)**2-3.0*(y-y0)**2/b**2 $ -3.0*(z-z0)**2/c**2) C JLTYP＝1，表示为体热源 JLTYP=1 IF(x .GE.(x0+d)) THEN FLUX(1)=heat1*shape1 ELSE FLUX(1)=heat2*shape2 ENDIF RETURN END

1，你需要把程序截图出来，我现在也记不清程序里的内容了，不然没法回答你。2，是的，你理解是对的。3，需要你把双椭球的区域用数学函数表示出来，进而编程，让其中心移动即可。4，你理解是对的，主要由于公式里是R2，因此不需要对其开方。