Workbench的焊接模擬過(guò)程（高斯移動(dòng)熱源）

2016-08-26 by:CAE仿真在線來(lái)源:互聯(lián)網(wǎng)

一直以來(lái),都是用ANSYS經(jīng)典做焊接模擬的,前兩天用ANSYS Workbench做了一個(gè)焊接模擬過(guò)程,我查閱了很多文獻(xiàn)和網(wǎng)絡(luò)資料,沒(méi)有人用Workbeach做焊接模擬,所以也不確定我做的是否完全正確,望高手斧正。

問(wèn)題闡述:三維平板堆焊焊接殘余應(yīng)力分析,先分析焊接溫度場(chǎng),然后轉(zhuǎn)化為結(jié)構(gòu)分析焊接殘余應(yīng)力,板材模型為0.2×0.2×0.006m,焊接速度10mm/s?？紤]其對(duì)稱性,建立了板材的一半模型進(jìn)行分析。
模型如圖所示。

一半模型

靠近焊縫的部位網(wǎng)格要細(xì)化。

網(wǎng)格劃分

對(duì)加熱的表面named selection為A1、A2,因?yàn)楹竺娴母咚篃嵩醇虞d采用的是APDL,為后面加高斯熱源作準(zhǔn)備。

加載熱源面

設(shè)置好瞬態(tài)熱分析的步長(zhǎng)和對(duì)流條件,插入Command定義高斯熱源。
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNCCSYS
*SET,_FNCNAME,'GAOSI'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*SET,_FNC_C1(1),2000
*SET,_FNC_C2(1),0.01
*SET,_FNC_C3(1),0.007
*SET,_FNCCSYS,0
! /INPUT,HANJIE.func,,,1
*DIM,%_FNCNAME%,TABLE,6,19,1,,,,%_FNCCSYS%
! Begin of equation: Qm*exp(-3*({X}^2 ({Y}-V*{TIME})^2)/R^2)
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 1, 18, 3, 1
*SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, 3, 2, -1
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,11,1), 0.0, -5, 0, 1, -4, 17, -1
*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, -3, 1, -5
*SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 1, 19, 17, -1
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 1, -3, 4, -2
*SET,%_FNCNAME%(0,17,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,18,1), 0.0, -2, 0, 1, 17, 3, -1
*SET,%_FNCNAME%(0,19,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: Qm*exp(-3*({X}^2 ({Y}-V*{TIME})^2)/R^2)
D,A1,TEMP,%GAOSI%
D,A2,TEMP,%GAOSI%瞬態(tài)熱分析設(shè)置完畢,然后就是求解,大概需要10分鐘。
不同時(shí)刻時(shí)的溫度場(chǎng)結(jié)果:

溫度分布

焊接溫度場(chǎng)的動(dòng)畫:

溫度場(chǎng)動(dòng)畫

然后把瞬態(tài)熱分析的結(jié)果轉(zhuǎn)化為瞬態(tài)結(jié)構(gòu)分析,把溫度場(chǎng)作為焊接殘余應(yīng)力分析的熱載荷。

在瞬態(tài)結(jié)構(gòu)分析模塊里,設(shè)置好焊接約束條件,設(shè)置好時(shí)間和步長(zhǎng),讀取每一步溫度載荷作為熱載荷進(jìn)行結(jié)構(gòu)分析,要注意設(shè)置熱載荷讀取的時(shí)間和結(jié)構(gòu)分析時(shí)間一致(默認(rèn)條件下結(jié)構(gòu)分析讀取熱分析的最后一步作為熱載荷)。
求解時(shí)間大概半個(gè)小時(shí),不同時(shí)刻的應(yīng)力場(chǎng)分布:

應(yīng)力場(chǎng)