ABAQUS平板對接-不帶生死單元。Python編寫移動高斯熱源子程序（包括高斯面熱源、雙橢球熱源）

模型作如下假設：材料為各向同性材料，不考慮熔池流動及相變影響。

考慮到過來學習的大多都是和我一樣的學生黨，因此設置了一個大家都能接受的價格。

如果視頻中有什么錯誤或沒講清的大家可以留言！！

捕獲.png

Welding-Temp_20170311141538.gif

Welding-Stress_20170311141621.gif

詳細操作視頻講解請查看：http://www.yqgqt.org.cn/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