ABAQUS焊接模拟-Python编写移动高斯热源子程序-不带生死单元

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

潜热     J/kg	固相线温度     oC	液相线温度     oC		对流换热系数J/m2/s/oC	辐射换热系数		绝对零度     oC	波尔兹曼常数
390000	615	655		80	0.85		-273.15	5.67E-08
热导率     W/m/oC	温度     oC		密度     kg/m3			弹性模量     Pa	泊松比	温度     oC
119	20		2700			6.67E+10	0.33	20
121	100					6.08E+10	0.33	100
126	200					5.68E+10	0.33	150
130	300					5.44E+10	0.33	200
138	400					5.1E+10	0.33	250
145	2000					4.31E+10	0.33	300
						3E+10	0.33	500
						1E+10	0.33	2000
热膨胀系数     m/m/oC	温度     oC		比热     J/kg/oC	温度     oC		屈服应力     Pa	塑性应变	温度
2.23E-05	20		900	20		2.5E+08	0	25
2.28E-05	100		921	100		2.4E+08	0	50
2.47E-05	200		1005	200		2.25E+08	0	100
2.55E-05	300		1047	300		1.9E+08	0	200
2.67E-05	400		1089	400		1.57E+08	0	250
2.70E-05	500		1129	2000		1.33E+08	0	300
2.70E-05	1000					54300000	0	350
						20800000	0	400
						9930000	0	450
						8600000	0	500
						5060000	0	550
						3300000	0	600
						1000	0	2000