ABAQUS平板对接-不带生死单元。Python编写移动高斯热源子程序（包括高斯面热源、双椭球热源）

模型作如下假设：材料为各向同性材料，不考虑熔池流动及相变影响。

考虑到过来学习的大多都是和我一样的学生党，因此设置了一个大家都能接受的价格。

如果视频中有什么错误或没讲清的大家可以留言！！

捕获.png

Welding-Temp_20170311141538.gif

Welding-Stress_20170311141621.gif

详细操作视频讲解请查看：http://www.jishulink.com/college/video/c10604

Python热源子程序：

1.双椭球热源

      SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,JLTYP,

     1                 TEMP,PRESS,SNAME)

C

      INCLUDE 'ABA_PARAM.INC'

      parameter(one=1.d0)

      DIMENSION COORDS(3),FLUX(2),TIME(2)

      CHARACTER*80 SNAME

C     Q,电弧有效热功率W

C     v,焊接速度m/s

C   d,当前时刻焊接斑点中心跟焊接初始位置的距离

    Q=800

    v=0.005

    d=v*TIME(2)

C

    x=COORDS(1)

    y=COORDS(2)

    z=COORDS(3)

C     焊接板厚度为0.005m,从坐标0,0,0.005开始，沿着x方向移动

    x0=0

    y0=0

    z0=0.005

C     a,b,c为椭球的半轴  

    a=0.0025

    b=0.0015

    c=0.002

    PI=3.1415

C

    heat=6*sqrt(3.0)*Q/(a*b*c*PI*sqrt(PI))

    shape=exp(-3*(x-x0-d)**2/a**2-3*(y-y0)**2/b**2-3*(z-z0)**2/c**2)

C     JLTYP＝1，表示为体热源

    JLTYP=1

    IF (KSTEP. eq. one) THEN

        FLUX(1)=heat*shape

      ENDIF

      RETURN

      END

2.高斯面热源

      SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,JLTYP,

     1                 TEMP,PRESS,SNAME)

C

      INCLUDE 'ABA_PARAM.INC'

      parameter(one=1.d0)

      DIMENSION COORDS(3),FLUX(2),TIME(2)

      CHARACTER*80 SNAME

C     Q,电弧有效热功率W

C     v,焊接速度m/s

C     Rh,加热斑点半径，95%的热量落在以Rh为半径的面积内

C   d,当前时刻焊接斑点中心跟焊接初始位置的距离

    Q=400

    v=0.005

    Rh=0.002

    d=v*TIME(2)

    x=COORDS(1)

    y=COORDS(2)

    z=COORDS(3)

C     焊接板厚度为0.005m,从坐标0,0开始，沿着x方向移动

    x0=0

    y0=0

C     

    PI=3.1415

C

    R=sqrt((x-x0-d)**2+(y-y0)**2)

C     JLTYP＝0，表示为面热源

    JLTYP=0

    IF (KSTEP. eq. one) THEN

        FLUX(1)=3*Q/(PI*Rh**2)*exp(-3*R**2/Rh**2)

      ENDIF

      RETURN

      END