模拟T型接头的温度场,算完结果后,感觉计算每一步,都会讲前一步的温度场删除,这是为什么?命令流如下:
/CLEAR
/FILNAME,T,1
/TITLE,WELD ANALYSIS BY "ELEMENT BIRTH AND DEATH"
/PREP7
/UNIT,SI
!*************定义单元类型及材料属性(Q345C)***********
ET,1,SOLID70,
MPTEMP,1,20,100,300,450,550,600 !温度
MPTEMP,7,720,800,1000,1200,1500,5000
MPDATA,KXX,1,1,51.9,51.1,46.1,41,37.5,35.6 !热导率
MPDATA,KXX,1,7,30.6,26,29.4,29.7,29.7,42
MPDATA,C,1,1,450,499,565,630,705,773 !比热
MPDATA,C,1,7,1080,931,437,400,400,400
MPDATA,PRXY,1,1,0.28,0.31,0.33,0.34,0.36,0.37 !泊松比
MPDATA,PRXY,1,7,0.37,0.42,0.47,0.49,0.49,0.49
MPDATA,ALPX,1,1,1.1E-5,1.1E-5,1.2E-5,1.3E-5,1.4E-5,1.4E-5 !膨胀系数
MPDATA,ALPX,1,7,1.4E-5,1.5E-5,1.5E-5,1.5E-5,1.5E-5,1.5E-5
MPDATA,EX,1,1,2.07E11,2.03E11,2E11,1.62E11,0.98E11,0.6E11 !弹性模量
MPDATA,EX,1,7,0.4E11,0.2E11,0.1E11,0.1E11,0.1E11,0.1E11,
!*************定义腹板和面板尺寸尺寸***********
l1=0.3 !面板长度
b1=0.2 !面板宽度
h1=0.016 !面板高度
l2=0.2 !腹板长度
b2=0.02 !腹板宽度
h2=0.2 !腹板高度
h3=0.01 !焊趾高度
arf=45 !焊趾角度
a=2*h3 !加密网格尺寸
a1=h3/2 !过度网格尺寸
!*************建立模型和网格划分***********
K,1,-l1/2,0,0
K,2,-l1/2,h1,0
K,3,-b2/2-a-a1,0,0
K,4,-b2/2-a-a1,h1,0
K,5,-b2/2-a,0,0
K,6,-b2/2-a,h1,0
K,7,b2/2+a,0,0
K,8,b2/2+a,h1,0
K,9,-b2/2,h1+h3,0
K,10,-b2/2-h3,h1,0
K,11,-b2/2,h1,0
K,12,b2/2,h1,0
K,13,-b2/2,h1+a,0
K,14,b2/2,h1+a,0
K,15,-b2/2,h1+a+a1,0
K,16,b2/2,h1+a+a1,0
K,17,-b2/2,h1+h2,0
K,18,b2/2,h1+h2,0
K,19,-l1/2,0,b1
/PNUM,KP,1
A,1,2,4,3
A,3,4,6,5
A,5,6,8,7
A,11,12,14,13
A,13,14,16,15
A,15,16,18,17
A,11,9,10
APLOT
/PNUM,KP,0
/PNUM,AREA,1
ASEL,S,AREA,,1,2,1
ASEL,A,AREA,,7
ARSYM,X,ALL, , , ,0,0
AGLUE,ALL
L,1,19
ALLSEL,ALL
VDRAG,ALL,,,,,,34
VGLUE,ALL
NUMMRG,ALL
ASEL,S,AREA,,64,68,4
ASEL,A,AREA,,7,10,3
LSLA,S,ALL
LESIZE,ALL,h3/5
ALLSEL,ALL
/PNUM,LINE,1
/PNUM,AREA,0
LSEL,S,LINE,,1,3,2
LSEL,A,LINE,,16,19,3
LSEL,A,LINE,,24,26,2
LESIZE,ALL,,,4
LSEL,S,LINE,,5,7,2
LSEL,A,LINE,,28,30,2
LSEL,A,LINE,,15,17,2
LESIZE,ALL,,,1
LSEL,S,LINE,,2,4,2
LSEL,A,LINE,,25,27,3
LSEL,A,LINE,,18,20,2
LESIZE,ALL,,,10,
ESIZE,b1/20
VSWEEP,ALL
!*************定义组件***********
/PNUM,VOLU,1
/PNUM,LINE,O
/TRLCY,VOLU,1,ALL
VPLOT
VSEL,S,,,7 !右侧焊缝
ESLV,S,
NSLV,S,1
CM,e1,ELEM
CM,n1,NODE
/PNUM,VOLU,0
/PNUM,AREA,1
ASEL,S,,,38
ESLA,S !右侧焊缝表面节点
NSLA,S,1
CM,a1,ELEM
CM,a1node,NODE
VSEL,S,,,11
ESLV,S,
NSLV,S,1
CM,e2,ELEM
CM,n2,NODE
ASEL,S,,,51
ESLA,S
NSLA,S,1
CM,a2,ELEM
CM,a2node,NODE
ALLSEL,ALL
ASEL,S,EXT !对流组件
ESLA,S
CM,aconv,ELEM
ALLSEL,ALL
!*************求解温度场***********
/SOLU
TUNIF,20 !初始温度20
LNSRCH,1 !打开线性搜索
TIMINT,1 !时间积分瞬态效应
AUTOTS,1 !自动时间步打开
ANTYPE,4, !瞬态分析
CMSEL,S,aconv !选择面组件aconv,需要施加对流边界的面
SFA,ALL,1,CONV,8,20 !给选择的面上施加对流换热系数
ALLSEL,ALL
CMSEL,S,e1 !选择组件e1,即代表第一道焊缝单元
CMSEL,A,e2 !选择组件e2,即代表第二道焊缝单元
EKILL,ALL
ESEL,S,LIVE
EPLOT
ALLSEL,ALL
j=0 !定义时间变量
*DO,i,0,0.18,0.02
j=j+6 !加热时间为6s,即热源每次移动时间停留6秒
CMSEL,S,n1 !选择第一道焊缝的所有节点
NSEL,R,LOC,z,i,i+0.02
ESLN,S,1 !选择与节点相关的单元
EALIVE,ALL !
BF,ALL,HGEN,1.41E9 !选定节点施加相同的热生成率,即采用均匀体热源的方式施加
ALLSEL,ALL
TIME,j
DELTIM,0.3
OUTRES,NSOL,ALL
KBC,1
SOLVE
BFDELE,ALL,HGEN
*ENDDO
|