abaqus 多道焊焊接热源问题

随风191225

一条焊缝分为四段，设四个分析步，只有第一段有热源，剩下的三段都没有热源，试了各种修改的程序都不对，哪位大神能给指点一下，感激不尽。附上我的程序：
      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

      wu=380.0        !电流
      wi=2.0          !电压
      effi=0.8        !焊接热效率
      v=0.01      !焊接速度
     
      q=wu*wi*effi  
        
      d=v*TIME(2)     !位移
     
      x=COORDS(1)
      y=COORDS(2)
      z=COORDS(3)
      
      x0=0            !热源1原点坐标
      y0=0
      z0=0.002
      
      x1=0        !热源2原点坐标
      y1=-0.01
      z1=0.002
      
       x2=0       !热源2原点坐标
      y2=-0.02
      z2=0.002
      
       x3=0         !热源2原点坐标
      y3=-0.03
      z3=0.002
      
          a=0.002     !热源形状参数
          b=0.002
          c=0.002
       
       PI=3.1415926
      
          heat=6.0*sqrt(3.0)*q/(a*b*c*PI*sqrt(PI))        !热源公式        
          shape1=exp(-3.0*(x-x0)**2/a**2-3.0*(y-y0+d)**2/b**2
     $        -3.0*(z-z0)**2/c**2)
     
       Shape2=exp(-3.0*(x-x1)**2/a**2-3.0*(y-y1+d)**2/b**2
     $        -3.0*(z-z1)**2/c**2)
     
      Shape3=exp(-3.0*(x-x2)**2/a**2-3.0*(y-y2+d)**2/b**2
     $        -3.0*(z-z2)**2/c**2)
     
      Shape4=exp(-3.0*(x-x3)**2/a**2-3.0*(y-y3+d)**2/b**2
     $        -3.0*(z-z3)**2/c**2)
     
C     JLTYP＝1，表示为体热源
      JLTYP=1
              IF(JSTEP .EQ.3) THEN        !3：分析步编号
          FLUX(1)=heat*shape1
          ENDIF
       IF(JSTEP .EQ.4) THEN        !4：分析步编号
          FLUX(1)=heat*shape2+0.04
          ENDIF
          IF(JSTEP .EQ.5) THEN        !5：分析步编号
          FLUX(1)=heat*shape2
          END IF
      IF(JSTEP .EQ.6) THEN        !6：分析步编号
          FLUX(1)=heat*shape3+0.04
          ENDIF
          IF(JSTEP .EQ.7) THEN        !7：分析步编号
          FLUX(1)=heat*shape3
          END IF
      IF(JSTEP .EQ.8) THEN        !8：分析步编号
                FLUX(1)=heat*shape4+0.04
          ENDIF
          IF(JSTEP .EQ.9) THEN        !9：分析步编号
          FLUX(1)=heat*shape4
          END IF
      RETURN
       END

b9638032

樓主讓你參考看看我的可以成功的
在我看來你的程序應該是沒有問題
你要不要檢查一下你CAE設定

heat=6*sqrt(3.0)*q/(a*b*c*PI*sqrt(PI))
        shape1=exp(-3*(x-x0-d)**2/a**2-3*(y-y0)**2/b**2-3*(z-z0)**2/c**2)
      
      shape2=exp(-3*(x-x1)**2/a**2-3*(y-y1+d)**2/b**2-3*(z-z0)**2/c**2)
      
      shape3=exp(-3*(x-x2-d)**2/a**2-3*(y-y2)**2/b**2-3*(z-z0)**2/c**2)
      
      shape4=exp(-3*(x-x3)**2/a**2-3*(y-y3-d)**2/b**2-3*(z-z0)**2/c**2)
      
      
      JLTYP=1
        if (JSTEP. eq. 1) then
      FLUX(1)=heat*shape1
      endif
      
      if (JSTEP. eq. 2) then
      FLUX(1)=heat*shape2
      endif
      
      if (JSTEP. eq. 3) then
      FLUX(1)=heat*shape3
      endif
      
      if (JSTEP. eq. 4) then
      FLUX(1)=heat*shape4
      endif
       RETURN
      END
     

diaoerlangdang

两位你们好，我有个问题想问你们一下。就是if判断里面的那个JSTEP是什么意思啊？

随风191225

JSTEP是分析步的意思

随风191225

b9638032 发表于 2016-9-2 15:56
樓主讓你參考看看我的可以成功的
在我看來你的程序應該是沒有問題
你要不要檢查一下你CAE設定

这个问题已经解决了，请问环形焊缝的模拟有没有做过啊，做过的话指导一下吧

ABAQUS初学者

随风191225 发表于 2016-9-28 15:09
这个问题已经解决了，请问环形焊缝的模拟有没有做过啊，做过的话指导一下吧 ...

楼主不说下是怎么出错的？

小兰子mst

确实没有问题

hhkayhh

:P:P:P:P

imrky

感谢楼主分享

liuzhen.neu

请问楼主，问题解决了吗？

石榴点心

谢谢楼主，学习了

chenhy6418802

厉害了
学习了

科研小白

IF(JSTEP .EQ.8) THEN        !8：分析步编号
                FLUX(1)=heat*shape4+0.04
          ENDIF

你好，请问这一步是什么意思呀？

ab62017

感謝樓主分享