我做的热结构耦合的一个简单算例,求版主加分。
!establish model,get corrodiation of the point,then draw the line&area.
/PREP7
K,1,0.010075,,, !~~~~~~~~~~~~~standard units,m,kg,n,s,oc
K,2,0.01221,,,
K,3,0.01375,0.002,,
K,4,0.01375,0.004,,
K,5,0.0148,0.004,,
K,6,0.0148,0.012,,
K,7,0.011075,0.012,,
K,8,0.011075,0.0105,,
K,9,0.010075,0.0105,,
K,10,0,0.0105,,
K,11,0,0.0045,,
K,12,0.010075,0.0045,,
K,13,0.010075,0.0033,,
K,14,0.010275,0.0033,,
K,15,0.010275,0.0013,,
K,16,0.010075,0.0013,,
K,17,0.010075,0.0095,,
K,18,0.012075,0.0095,,
K,19,0.012075,0.009,,
K,20,0.010075,0.009,,
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~lines establishment
LSTR, 1, 2
LSTR, 2, 3
LSTR, 3, 4
LSTR, 4, 5
LSTR, 5, 6
LSTR, 6, 7
LSTR, 7, 8
LSTR, 8, 9
LSTR, 9, 10
LSTR, 10, 11
LSTR, 11, 12
LSTR, 12, 13
LSTR, 13, 14
LSTR, 14, 15
LSTR, 15, 16
LSTR, 16, 1
LSTR, 9, 17
LSTR, 17, 18
LSTR, 18, 19
LSTR, 19, 20
LSTR, 20, 12
LSTR, 17, 20
FLST,2,18,4
FITEM,2,12
FITEM,2,13
FITEM,2,14
FITEM,2,15
FITEM,2,16
FITEM,2,1
FITEM,2,2
FITEM,2,3
FITEM,2,4
FITEM,2,5
FITEM,2,6
FITEM,2,7
FITEM,2,8
FITEM,2,17
FITEM,2,18
FITEM,2,19
FITEM,2,20
FITEM,2,21
AL,P51X
FLST,2,6,4
FITEM,2,9
FITEM,2,17
FITEM,2,22
FITEM,2,21
FITEM,2,11
FITEM,2,10
AL,P51X
/REPLOT,RESIZE
!define the parameter of streel 304 and Na-Ga glass
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,1,,16.2
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,1,,502
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,1.93e3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,2.52e3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,2,,800
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,2,,1.03
!~~~~~~~~~~~~~~~~~~~slect the FEA element PLANE55
ET,1,PLANE55
KEYOPT,1,1,0
KEYOPT,1,3,1
KEYOPT,1,4,0
KEYOPT,1,8,0
KEYOPT,1,9,0
KEYOPT,1,1,0
KEYOPT,1,3,1
KEYOPT,1,4,0
KEYOPT,1,8,0
KEYOPT,1,9,0
!~~~~~~~~~~~~~~~~~~~~Mesh
ESIZE,0.0005,0,
MSHAPE,0,2D
MSHKEY,0
TYPE, 1
MAT, 1
asel,s,,,1
AMESH,all
MAT, 2
asel,s,,,2
AMESH,all
finish
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
/solu !coming solution
ALLSEL
NSEL,ALL
IC,ALL,TEMP,30 !define the intial temperature
FLST,2,15,4,ORDE,4
FITEM,2,1
FITEM,2,-9
FITEM,2,11
FITEM,2,-16
DL,P51X, ,TEMP,150,0 !define the boundary condition with 150oc.
antype,trans !analysis type
solcontrol,on !激活非线性求解优化
!~~~~~~~~~~~~~~~~~~~瞬态非线性热分析求解选项
time,10 !定义载荷步终止时间,本步时间长度为 8 秒
deltim,0.1,0.1,2.0,on !设定初始时间及最小、最大时间间隔(即开始求解时间、最小步长、最大步长)
kbc,1 !指定为阶跃载荷
autots,on !打开自动时间步长开关
timint,on !打开时间积分效应
tintp,,,,1 !指定为 1 阶积分参数
lnsrch,1 !打开线性搜索
outres,nsol,1 !每隔1个子步输出节点解
allsel,all
save
solve
finish
/post1
set,last
/graphics,full
plnsol,temp,,0 !绘制温度云图
!plnsol,tg,sum,0 !绘制温度梯度云图
!finish
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~mechanical analysis
/prep7
ET,2,PLANE182
!结构材料定义
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,1.93e11
MPDATA,PRXY,1,,0.27
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,1,REFT,,,
MPDATA,ALPX,1,,17.3e-5
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,2,REFT,,,
MPDATA,ALPX,2,,5.10e-7
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,0.716e11
MPDATA,PRXY,2,,0.23
allsel,all
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~change the elemet typle from planel55 to plane182 which is mechanical element
emodif,all,type,2
finish
/solu
antype,0
/pbf,temp,,1
LDREAD,TEMP,last,,, ,'file','rth',' '
!原来的,增加1,不能运算结果 LDREAD,TEMP,last,, ,1,'file','rth',' '
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~fix the bottom face
FLST,2,1,4,ORDE,1
FITEM,2,1
DL,P51X, ,ALL,0
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~define the symmetry displement on the axi-line
DL, 10, ,SYMM
alls
solve
finish
/post1
set,last
/graphics,full
plnsol,u,sum,0,1.0 !绘制总位移
plnsol,s,eqv,0,1.0 !绘制等效应力
/CONT,1,9,AUTO
/REPLOT |