資源描述:
《ansys高斯熱源施加》由會員上傳分享,免費在線閱讀,更多相關內容在行業(yè)資料-天天文庫。
1、!平板對接開v型坡口多道焊成型(溫度場)!采用guass表面熱源finish/clear,nostart/filname,1.11thermal_analysis/title,3D_weld_thermal_analysis/units,si!定義國際制單位/prep7!**************************************************!定義焊接尺寸及焊接參數(shù)!**************************************************!********************
2、**********************構件尺寸len=100e-3!長0.1mwid=50e-3!寬0.05mthic=6e-3!厚0.006m!b=0.002!根部預留(m)!h=0!鈍邊(m)wweld=15e-3!焊接熱影響區(qū)寬度lnum=10!沿焊縫長度方向劃分份數(shù)!******************************************焊接參數(shù)I=200!焊接電流(A)U=30!焊接電壓(V)v=0.5/60!焊接速度(m/s)arf=60!坡口角度(°)pi=3.1415926!pi為圓周率eta=0.7
3、!焊接熱效率,手工電弧焊0.7,埋弧焊0.8K=1.2e4!熱流集中程度系數(shù)(/m^2)?them0=25!環(huán)境溫度及母材初始溫度,取室溫25℃!**********************************************定義高斯表面移動熱源參數(shù)!GAUSS熱源模型qr=qmax*exp(-k*r*r)!(單位W/m2)!假設熱源中心坐標為(a,b,c),a=0,b=每個焊縫的最高點的坐標y,c=v*dt!每個焊縫上表面任一點的坐標為(x,y,z),x>0,y>0,z>0!r^2=x^2+(y-b)^2+(z-c)^2
4、,每個焊縫上表面任一點距熱源中心距離qmax=0.24*k*eta*I*u/pi!加熱斑點中心熱流(W/m2)percent=0.99radius=abs(sqrt(1/k*log(1/(1-percent))))!有效熱半徑,單位為m!**********************************************************************!定義單元類型及材料屬性!********************************************************************
5、**et,1,solid70!計算單元類型!nothing!定義實常數(shù)mptemp,1,0,100,200,300,400,500!定義材料屬性溫度范圍mptemp,,1200,1250,1350,1400,1450,1500mptemp,,1505,1705,1905,2105,2305,2500mpdata,dens,1,1,7820,7800,7800,7800,7800,7800!定義材料密度(kg/m3)mpdata,dens,1,,7800,7800,7800,7800,7800,7800mpdata,dens,1,,78
6、00,7800,7800,7800,7800,7800mpdata,kxx,1,1,52,50.7,48.6,46.1,42.3,38.9!熱傳導系數(shù)(W/(m*K))mpdata,kxx,1,,30,30,30,30,30,30mpdata,kxx,1,,30,30,30,30,30,30mpdata,c,1,1,450,469,481,508.5,536,569!考慮相變潛熱的比熱容(J/kg*K)mpdata,c,1,,700,2172,5116,6589,8061,9533mpdata,c,1,,9533,7757,5982,
7、4206,2431,700!****************************************************************!建立有限元計算模型!*****************************************************************afun,deg!如果涉及到角度,將使用“度”(系統(tǒng)默認是弧度)K,1,0,0,0K,2,wweld,0,0K,3,wid,0,0K,4,wid,thic,0K,5,thic*tan(arf/2)+wweld,thic,0K,6,t
8、hic*tan(arf/2),thic,0K,7,0,0,lenA,1,3,4,6A,1,2,5,6cyl4,0,0,thic/cos(arf/2),arf,0,90cyl4,0,0,thic/cos(arf/2)*sq