按键盘上方向键 ← 或 → 可快速上下翻页，按键盘上的 Enter 键可回到本书目录页，按键盘上方向键 ↑ 可回到本页顶部！
————未阅读完？加入书签已便下次继续阅读！
alls
outres；all；all
solve
save；'1。11thermal_analysis_initial_temp'；'db'
！查看初始温度场，应该全为25摄氏度
！***************************************************************
！按顺序模拟温度场，以第一道为例
！***************************************************************
！**************************************************第1道焊缝填充
esel；s；type；；1 ！在第一类单元中加热流载荷
esel；r；live
eplot
*do；im；0；lnum；1
！****************************选择有效半径内节点为当前节点
vsel；s；；；v1 ！选择焊缝1上表面单元
eslv；r
nsle
asel；s；；；a1
nsla；r；1
esln；r；0
！热源加载； 将各段后点（0…lnum）为热源中心；先加载后删除热源
！即先消除上段所加高斯热源；并将上段的温度值作为下段的初始值
！*******************************************包括起始加热点
c=V*tinc*im ！热源位置
b=thic/cos（arf/2）*sqrt（2）/2
tm=dt1＋im*tinc
time；tm
antype；4；rest
timint；on
autots；on
kbc；1
deltim；0。01；0。01；0。2 ！时间步长
！***********************************以下为施加热流密度载荷
*do；i；nemin；nemax；1
*if；esel（i）；eq；1；then ！如果该单元在上面的集合里，则
xsy=centrx（i） ！读取该单元的中心坐标
ysy=centry（i）
zsy=centrz（i）
rr=abs（sqrt（xsy*xsy＋（ysy…b）*（ysy…b）＋（zsy…c）*（zsy…c））） ！该单元中心距离热源中心
*if；rr；le；radius；then ！该单元中心在加热半径范围内
qr=qmax*exp（…k*rr*rr） ！该单元中心处的热流大小
sfe；i；nmface（i）；hflux；；qr ！在每个单元所指定的面上施加热源载荷
*endif
*endif
*enddo
allsel ！全选求解
outres；all；all
solve
esel；s；type；；1
esel；r；live
vsel；s；；；v1 ！选择焊缝1上表面单元
eslv；r
nsle
asel；s；；；a1
nsla；r；1
esln；r；0
！！！！！！！！！！！！！！输入先暂时停止一下
！插入〃/post1plnsol；temp/psf；hflux；；1；1；on/replot〃的内容进行查看
！问题：热流载荷的分布和温度场的分布不一致（该完全一致才对），有点偏移？
！若要继续，读入〃finish/solu〃
*do；i；1；6
sfedele；all；i；hflux ！删除每个单元六个面上的热载荷
*enddo
*enddo
alls
esel；s；live
eplot ！查看温度场
save；'1。11thermal_analysis_alive1'；'db'
！=
！生成节点温度分布云图动画文件的程序段
！=
/post1
/seg；dele ！允许图形数据存储到局部终端存储器
/dscale；1；1。0 ！显示位移时，设置放大系数
avprin；0；0 ！规定主项量和总向量如何被计算
avres；1 ！规定当powergraphies可用时，结果数据如何被平均
/seg；multi；1。11temp…10；0。2 ！存动画；分节存储后来的显示；动画名字；每个画面间延迟系数
esel；s；live
！*************************************************读取初始稳态温度场结果
set；；；1；；dt0 ！从结果文件读；载荷步；子步；比例因子实部；比例因子虚部；时间
plnsol；temp
！*****************************************读取第1道焊缝过程中的温度场结果
*do；im；0；lnum；1
tm=dt1＋im*tinc
set；；；1；；tm
plnsol；temp
*enddo
/seg；off；1。11temp…10；0。2
anim；1；1 ！*********************************生成动画完毕
！****************************************存储第1道焊缝过程中的温度场结果图片
esel；s；live
*do；im；0；lnum；1
*if；im；eq；0；then
tm=dt1
set；；；1；；tm
plnsol；temp
/image；save；1。11temp0；jpeg
*endif
tm=dt1＋im*tinc
set；；；1；；tm
plnsol；temp
*if；tm；eq；dt1＋lnum/2*tinc；then
set；；；1；；tm
plnsol；temp6。5；jpeg
*endif
*enddo
/image；save；1。11temp12。5；jpeg
！温度场完毕以第一道为例