双椭球模型

/PREP7
!******************************
!定义参数
!******************************
LENGTH=0.11
WIDTH=0.05
HEIGHT=0.01
DIA=0.006
ah=DIA/2
bh=0.0025
chf=2*ah
chb=6*ah
rf=2*chf/(chf+chb)
rb=2*chb/(chf+chb)
U=38
I=180
EFF=0.85
W=EFF*U*I
STEP=0.002
COUNT=0.03/STEP
V=0.00235
!******************************
!定义单元类型
!******************************
ET,1,SOLID70
ET,2,Surf152
DOF,TEMP
KEYOPT,2,3,0
KEYOPT,2,4,1
KEYOPT,2,5,0
KEYOPT,2,6,0
KEYOPT,2,7,0
KEYOPT,2,8,3 !根据固体表面温度计算对流系数KEYOPT,2,9,0
!******************************
!定义材料参数
MPTEMP,,,,,,,,
MPTEMP,1,25
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400 MPTEMP,6,500 MPTEMP,7,600 MPTEMP,8,700 MPTEMP,9,800 MPTEMP,10,900 MPTEMP,11,1000 MPTEMP,12,1100 MPTEMP,13,1200 MPTEMP,14,1300 MPTEMP,15,1400 MPTEMP,16,1500 MPTEMP,17,2000 MPDATA,KXX,1,,24.5 MPDATA,KXX,1,,26.5 MPDATA,KXX,1,,27.0 MPDATA,KXX,1,,27.9 MPDATA,KXX,1,,28.1 MPDATA,KXX,1,,28.4 MPDATA,KXX,1,,28.3 MPDATA,KXX,1,,30.5 MPDATA,KXX,1,,31.2 MPDATA,KXX,1,,28.1 MPDATA,KXX,1,,28.0 MPDATA,KXX,1,,28.2 MPDATA,KXX,1,,28.3 MPDATA,KXX,1,,28.3 MPDATA,KXX,1,,29.5 MPDATA,KXX,1,,59.1 MPDATA,KXX,1,,59.1 MPTEMP,,,,,,,, MPTEMP,1,25 MPTEMP,2,100 MPTEMP,3,200 MPTEMP,4,300 MPTEMP,5,400 MPTEMP,6,500 MPTEMP,7,600 MPTEMP,8,700 MPTEMP,9,750 MPTEMP,10,800 MPTEMP,11,900 MPTEMP,12,1000 MPTEMP,13,1100
MPTEMP,15,1300 MPTEMP,16,1400 MPTEMP,17,1500 MPTEMP,18,2000 MPDATA,C,1,,431 MPDATA,C,1,,483 MPDATA,C,1,,521 MPDATA,C,1,,573 MPDATA,C,1,,626 MPDATA,C,1,,693 MPDATA,C,1,,795 MPDATA,C,1,,966 MPDATA,C,1,,1100 MPDATA,C,1,,952 MPDATA,C,1,,748 MPDATA,C,1,,676 MPDATA,C,1,,722 MPDATA,C,1,,738 MPDATA,C,1,,739 MPDATA,C,1,,741 MPDATA,C,1,,750 MPDATA,C,1,,750 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7780 MPTEMP,,,,,,,, MPTEMP,1,25 MPTEMP,2,100 MPTEMP,3,200 MPTEMP,4,300 MPTEMP,5,400 MPTEMP,6,500 MPTEMP,7,600 MPTEMP,8,700 MPTEMP,9,750 MPTEMP,10,800 MPTEMP,11,850 MPTEMP,12,900 MPTEMP,13,1000 MPTEMP,14,1100 MPTEMP,15,1200 MPTEMP,16,1300 MPTEMP,17,1400
MPTEMP,19,2000
MPDATA,ENTH,1,,0
MPDATA,ENTH,1,,2.6666E8 MPDATA,ENTH,1,,6.672E8 MPDATA,ENTH,1,,1.08278E9 MPDATA,ENTH,1,,1.54919E9 MPDATA,ENTH,1,,2.06228E9 MPDATA,ENTH,1,,2.64112E9 MPDATA,ENTH,1,,3.32621E9 MPDATA,ENTH,1,,3.72817E9 MPDATA,ENTH,1,,4.12757E9 MPDATA,ENTH,1,,4.47554E9 MPDATA,ENTH,1,,4.78373E9 MPDATA,ENTH,1,,5.33803E9 MPDATA,ENTH,1,,5.88209E9 MPDATA,ENTH,1,,6.45017E9 MPDATA,ENTH,1,,7.02509E9 MPDATA,ENTH,1,,7.60111E9 MPDATA,ENTH,1,,8.18076E9 MPDATA,ENTH,1,,1.10963E10
!综合换热系数MPTEMP,,,,,,,,
MPTEMP,1,25
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,700
MPTEMP,9,800
MPTEMP,10,900
MPTEMP,11,1000
MPTEMP,12,1100
MPTEMP,13,1200 MPTEMP,14,1300
MPTEMP,15,1400
MPTEMP,16,1500
MPTEMP,17,2000
MPDATA,HF,1,,1.67
MPDATA,HF,1,,6.68
MPDATA,HF,1,,13.36
MPDATA,HF,1,,20.04
MPDATA,HF,1,,26.72
MPDATA,HF,1,,33.4
MPDATA,HF,1,,56.5
MPDATA,HF,1,,79.6
MPDATA,HF,1,,102.7
MPDATA,HF,1,,125.8
MPDATA,HF,1,,148.9
MPDATA,HF,1,,172
MPDATA,HF,1,,195.1
MPDATA,HF,1,,218.2
MPDATA,HF,1,,241.3
MPDATA,HF,1,,264.4
MPDATA,HF,1,,379.9
SAVE
!******************************
!建立几何模型
!******************************
WPSTYL,,,,,,,,0 !设定工作平面风格BLOCK,0,-WIDTH,0,LENGTH+0.001,0,HEIGHT, !生成长方体焊材1 wpoff,-0.001,,
wpoff,0,0,0.009
VSBW, 1
wprot,0,0,120
VSBW, 2
wprot,0,0,-120
wpoff,0,0,-0.0015
VSBW, 4
VDELE, 2, , ,1
KWPAVE, 16
wprot,0,0,90
VSBW, 3
wprot,0,0,-90
wpoff,0.0015,0,0
KWPAVE, 20
VSBW, 1
save
wpoff,-0.008,0,0
wprot,0,0,90
VSBW, 2
KWPAVE, 26
wprot,0,0,30
VSBW, 3
VSBW, 6
wprot,0,0,-30
wpoff,0,0,-0.01
VSBW, 7
wprot,0,0,30
VSBW, 2
VSBW, 3
!实体模型切割完毕进行网格划分
save
FLST,5,2,6,ORDE,2
!5代表所选实体对应操作命令的第四个参数，2代表后续FITEM的行数，6代表体，orde为有序，2代表FITEM的个数
FITEM,5,4
FITEM,5,-5 !选择4,5体，负5和前面FITEM命令的4构成一个范围，即4到5；CM,_Y,VOLU
VSEL, , , ,P51X !第四个参数为P51X
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 1, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!焊缝单元
LESIZE,22, , ,4, , , , ,1
LESIZE,25, , ,4, , , , ,1
LESIZE,18, , ,4, , , , ,1
LESIZE,32, , ,4, , , , ,1
LESIZE,36, , ,4, , , , ,1
LESIZE,38, , ,6, , , , ,1
LESIZE,40, , ,6, , , , ,1
EXTOPT,ACLEAR,1
EXTOPT,VSWE,TETS,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,110,0
!*
VSWEEP,4,18,17
VSWEEP,5,25,24
!热影响区单元划分
FLST,5,3,6,ORDE,3
FITEM,5,1
FITEM,5,8
FITEM,5,-9
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 1, , 1, 0 CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
LESIZE,48, , ,4, , , , ,1 LESIZE,10, , ,6, , , , ,1 LESIZE,59, , ,16, , , , ,1 LESIZE,61, , ,16, , , , ,1 LESIZE,47, , ,16, , , , ,1 LESIZE,13, , ,16, , , , ,1 LESIZE,33, , ,16, , , , ,1 LESIZE,56, , ,16, , , , ,1 EXTOPT,ESIZE,110,0 VSWEEP,1,7,5
VSWEEP,9,40,19
VSWEEP,8,38,37
!母材网格划分
FLST,5,3,6,ORDE,3
FITEM,5,6
FITEM,5,-7
FITEM,5,12
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 1, , 1, 0 CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
LESIZE,66, , ,2, , , , ,1 LESIZE,3, , ,2, , , , ,1 LESIZE,43, , ,2, , , , ,1 LESIZE,21, , ,2, , , , ,1 LESIZE,74, , ,4, , , , ,1
LESIZE,7, , ,4, , , , ,1 LESIZE,51, , ,10, , , , ,1 LESIZE,73, , ,10, , , , ,1 LESIZE,64, , ,10
, , , , ,1 LESIZE,30, , ,10, , , , ,1 EXTOPT,ESIZE,40,0 VSWEEP,6,41,42
VSWEEP,12,53,10
VSWEEP,7,33,32
!过渡区网格划分
MSHKEY,0
MSHAPE,1,3d
VMESH,10
MSHKEY,0
MSHAPE,1,3d
VMESH,2
MSHKEY,0
MSHAPE,1,3d
VMESH,11
wprot,0,0,-30
KWPAVE, 18
FLST,3,11,6,ORDE,4
FITEM,3,1
FITEM,3,-2
FITEM,3,4
FITEM,3,-12
VSYMM,X,P51X, , , ,0,0 FLST,2,11,6,ORDE,4
FITEM,2,1
FITEM,2,-2
FITEM,2,4
FITEM,2,-12
VCLEAR,P51X
FLST,2,11,6,ORDE,4
FITEM,2,1
FITEM,2,-2
FITEM,2,4
FITEM,2,-12
VDELE,P51X, , ,1
ALLSEL,ALL
ASEL,S,EXT
FLST,5,2,5,ORDE,2
FITEM,5,66
FITEM,5,73
ASEL,U, , ,P51X
ALLSEL,BELOW,AREA
FLST,5,34,5,ORDE,17
FITEM,5,12
FITEM,5,56
FITEM,5,-57
FITEM,5,59
FITEM,5,62
FITEM,5,64
FITEM,5,67
FITEM,5,-70
FITEM,5,72
FITEM,5,76
FITEM,5,-80
FITEM,5,83
FITEM,5,-85
FITEM,5,87
FITEM,5,-88
FITEM,5,90
FITEM,5,-102
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 1, , 2, 0, CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
amesh,all
allsel,all
esel,s,type,,2 !
sf,all,conv,-1,25 ! allsel,all
nsel,all
finish
******************************
!循环加载、求解
!******************************
!施加载荷
*DO,I,0,COUNT !循环求解开始
DISY=I*STEP !定义参数disy
T=(DISY+STEP)/V !定义参数t
/SOLU !进入求解模块
ANTYPE,4 !设定分析类型
TRNOPT,FULL !设定为完全瞬态分析
TUNIF,200, !施加均匀初始温度
OUTRES,BASIC,LAST !输出最后一个子步的基本数据TIME,T !设定求解结束时间
NROPT,FULL,,ON !设定Newton-Raphson选项AUTOTS,1 !打开自动时间步
NSUBST,5,,,1 !设定子步数5
KBC,1 !阶跃加载方式
LNSRCH,1 !打开线性搜索
TSRES,ERASE !清除时间步选项
TIMINT,1,THERM !设定时间积分
!*******************************
!焊接过程开始
!*******************************
!双椭球热源的加载
!*******************************
!*******************************
!前半球节点的选取
!*******************************
LOCAL,11,2,0,DISY,0.0075, , , ,chf/ah,bh/ah, !生成局部坐标系11 NSEL,S,LOC,X,,ah !选择x=0~dia的节点
NSEL,R,LOC,Y,,90 !在当前选择集中再选出y=0~90的节点NSEL,R,LOC,Z,,90 !在当前选择集中再选出z=0～90的节点CM,NODEF,NODE !生成节点组件nodef
!*******************************
!前半球函数、表格定义
!**
*****************************
LOCAL,12,0,0,DISY,0.0075, , , ,1,1, !生成局部坐标系12
!删除以下参数变量
*DEL,TABLENAME
*DEL,TABLEMTID
*DEL,TABLE_ROW1
*DEL,TABLE_ROW2
*DEL,TABLE_ROW3
*DEL,TABLE_ROW4
*DEL,TABLE_ROW5
*DEL,TABLECSYS
*DEL,'TEMPF'
!定义表格各行
*SET,TABLENAME,'TEMPF'
*DIM,TABLE_ROW1,,1
*DIM,TABLE_ROW2,,1
*DIM,TABLE_ROW3,,1
*DIM,TABLE_ROW4,,1
*DIM,TABLE_ROW5,,1
!表格各行赋值
*SET,TABLE_ROW1(1),rf
*SET,TABLE_ROW2(1),w
*SET,TABLE_ROW3(1),chf
*SET,TABLE_ROW4(1),ah
*SET,TABLE_ROW5(1),bh
*SET,TABLECSYS,12
!表格各数据赋值
*DIM,%TABLENAME%,TABLE,7,38,1,,,,%TABLECSYS% !定义一个载荷表格*SET,%TABLENAME%(0,0,1), 0.0,-999 !表格赋值
*SET,%TABLENAME%(2,0,1), 0.0
*SET,%TABLENAME%(3,0,1), %TABLE_ROW1(1)%
*SET,%TABLENAME%(4,0,1), %TABLE_ROW2(1)%
*SET,%TABLENAME%(5,0,1), %TABLE_ROW3(1)%
*SET,%TABLENAME%(6,0,1), %TABLE_ROW4(1)%
*SET,%TABLENAME%(7,0,1), %TABLE_ROW5(1)%
*SET,%TABLENAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%TABLENAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%TABLENAME%(0,3,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%TABLENAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%TABLENAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%TABLENAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%TABLENAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1
*SET,%TABLENAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%TABLENAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%TABLENAME%(0,10,1), 0.0, -3, 0, 1, 2, 17, -1
*SET,%TABLENAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%TABLENAME%(0,12,1), 0.0, -5, 0, 1, 19, 17, -1
*SET,%TABLENAME%(0,13,1), 0.0, -1, 0, 1, -3, 4, -5
*SET,%TABLENAME%(0,14,1), 0.0, -3, 0, 2, 0, 0, 3
*SET,%TABLENAME%(0,15,1), 0.0, -5, 0, 1, 3, 17, -3
*SET,%TABLENAME%(0,16,1), 0.0, -3, 0, 2, 0, 0, 20
*SET,%TABLENAME%(0,17,1), 0.0, -6, 0, 1, 20, 17, -3
*SET,%TABLENAME%(0,18,1), 0.0, -3, 0, 1, -5, 4, -6
*SET,%TABLENAME%(0,19,1), 0.0, -5, 0, 1, -1, 1, -3
*SET,%TABLENAME%(0,20,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%TABLENAME%(0,21,1), 0.0, -3, 0, 1, 4, 17, -1
*SET,%TABLENAME%(0,22,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%TABLENAME%(0,23,1), 0.0, -6, 0, 1, 21, 17, -1
*SET,%TABLENAME%(0,24,1), 0.0, -1, 0, 1, -3, 4, -6
*SET,%TABLENAME%(0,25,1), 0.0, -3, 0, 1, -5, 1, -1
*SET,%TABLENAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%TABLENAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -3
*SET,%TABLENAME%(0,28,1), 0.0, -1, 0, 1, -4, 3, -5
*SET,%TABLENAME%(0,29,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%TABLENAME%(0,30,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%TABLENAME%(0,31,1), 0.0, -1, 0, 1, 20, 3, 21
*SET,%TABLENAME%(0,32,1), 0.0, -2, 0, 1, -1, 3, 19
*SET,%TABLENAME%(0,33,1), 0.0, -1, 0, 3.14, 0, 0, -2
*SET,%TABLENAME%(0,34,1), 0.0, -4, 0, 1, -2, 3, -1
*SET,%TABLENAME%(0,35,1), 0.0, -1, 16, 3.14, 0, 0, 0
*SET,%TABLENAME%(0,36,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%TABLE
NAME%(0,37,1), 0.0, -1, 0, 1, -3, 4, -2
*SET,%TABLENAME%(0,38,1), 0.0, 99, 0, 1, -1, 0, 0
!施加表格热源载荷
BF,NODEF,HGEN,%TEMPF% !在节点组件施加热生成载荷CMDELE,NODEF !删除节点组件nodef
ALLSEL,ALL !选择全部图元
!*******************************
!后半球节点的选取
!*******************************
LOCAL,13,2,0,DISY,0.0075, , , ,chb/ah,bh/ah, !生成局部坐标系13 NSEL,S,LOC,X,,ah !选择x=0~dia的节点
NSEL,R,LOC,Y,-90,0 !在当前选择集中再选出y=-90~0的节点NSEL,R,LOC,Z,,90 !在当前选择集中再选出z=0～90的节点CM,NODEB,NODE !生成节点组件nodeb
!*******************************
!后半球函数、表格定义
!*******************************
LOCAL,14,0,0,DISY,0, , , ,1,1, !生成局部坐标系14
!删除以下参数变量
*DEL,TABLENAME
*DEL,TABLEMTID
*DEL,TABLE_ROW1
*DEL,TABLE_ROW2
*DEL,TABLE_ROW3
*DEL,TABLE_ROW4
*DEL,TABLE_ROW5
*DEL,TABLECSYS
*DEL,'TEMPB'
!定义表格各行
*SET,TABLENAME,'TEMPB'
*DIM,TABLE_ROW1,,1
*DIM,TABLE_ROW2,,1
*DIM,TABLE_ROW3,,1
*DIM,TABLE_ROW4,,1
*DIM,TABLE_ROW5,,1
!各行赋初值
*SET,TABLE_ROW1(1),rb
*SET,TABLE_ROW2(1),w
*SET,TABLE_ROW3(1),chb
*SET,TABLE_ROW4(1),ah
*SET,TABLE_ROW5(1),bh
*SET,TABLECSYS,14
!表格各数据赋值
*DIM,%TABLENAME%,TABLE,7,38,1,,,,%TABLECSYS% !定义一个表格载荷*SET,%TABLENAME%(0,0,1), 0.0,-999 !表格赋值
*SET,%TABLENAME%(2,0,1), 0.0
*SET,%TABLENAME%(3,0,1), %TABLE_ROW1(1)%
*SET,%TABLENAME%(4,0,1), %TABLE_ROW2(1)%
*SET,%TABLENAME%(5,0,1), %TABLE_ROW3(1)%
*SET,%TABLENAME%(6,0,1), %TABLE_ROW4(1)%
*SET,%TABLENAME%(7,0,1), %TABLE_ROW5(1)%
*SET,%TABLENAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%TABLENAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%TABLENAME%(0,3,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%TABLENAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%TABLENAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%TABLENAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%TABLENAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1
*SET,%TABLENAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%TABLENAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%TABLENAME%(0,10,1), 0.0, -3, 0, 1, 2, 17, -1
*SET,%TABLENAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%TABLENAME%(0,12,1), 0.0, -5, 0, 1, 19, 17, -1
*SET,%TABLENAME%(0,13,1), 0.0, -1, 0, 1, -3, 4, -5
*SET,%TABLENAME%(0,14,1), 0.0, -3, 0, 2, 0, 0, 3
*SET,%TABLENAME%(0,15,1), 0.0, -5, 0, 1, 3, 17, -3
*SET,%TABLENAME%(0,16,1), 0.0, -3, 0, 2, 0, 0, 20
*SET,%TABLENAME%(0,17,1), 0.0, -6, 0, 1, 20, 17, -3
*SET,%TABLENAME%(0,18,1), 0.0, -3, 0, 1, -5, 4, -6
*SET,%TABLENAME%(0,19,1), 0.0, -5, 0, 1, -1, 1, -3
*SET,%TABLENAME%(0,20,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%TABLENAME%(0,21,1), 0.0, -3, 0,
1, 4, 17, -1
*SET,%TABLENAME%(0,22,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%TABLENAME%(0,23,1), 0.0, -6, 0, 1, 21, 17, -1
*SET,%TABLENAME%(0,24,1), 0.0, -1, 0, 1, -3, 4, -6
*SET,%TABLENAME%(0,25,1), 0.0, -3, 0, 1, -5, 1, -1
*SET,%TABLENAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%TABLENAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -3
*SET,%TABLENAME%(0,28,1), 0.0, -1, 0, 1, -4, 3, -5
*SET,%TABLENAME%(0,29,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%TABLENAME%(0,30,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%TABLENAME%(0,31,1), 0.0, -1, 0, 1, 20, 3, 21
*SET,%TABLENAME%(0,32,1), 0.0, -2, 0, 1, -1, 3, 19
*SET,%TABLENAME%(0,33,1), 0.0, -1, 0, 3.14, 0, 0, -2
*SET,%TABLENAME%(0,34,1), 0.0, -4, 0, 1, -2, 3, -1
*SET,%TABLENAME%(0,35,1), 0.0, -1, 16, 3.14, 0, 0, 0
*SET,%TABLENAME%(0,36,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%TABLENAME%(0,37,1), 0.0, -1, 0, 1, -3, 4, -2
*SET,%TABLENAME%(0,38,1), 0.0, 99, 0, 1, -1, 0, 0
!
BF,NODEB,HGEN,%TEMPB% !在节点组件nodeb上施加热生成载荷CMDELE,NODEB !删除节点组件nodeb
ALLSEL,ALL !选择全部图元
SOLVE !开始求解
BFDELE,ALL,HGEN !删除所有节点上的热生成载荷
*ENDDO !循环求解结束
SAVE !保存数据库。