1、分析例子为
2、需要注意的是单元的匹配，比如用77单元采用8节点网格，那么结构单元也须采用8节点单元183。否则会出现错误或警告。
3、apdl运行通过，单元网格尺寸改为0.01
4、apdl通过运行，结果正确无警告。
!*Thermal-Structural Coupling Analysis
fini
/cle
/filname,exce,1
/title,lovz
/prep7
!physics thermal
et,1,plane77,,,1
mp,kxx,1,2.2
mp,kxx,2,10.8
rect,0.1875,0.4,0,0.05
rect,0.4,0.6,0,0.05
aglue,all
numcmp,all
asel,,,,1
aatt,1,1,1
asel,,,,2
aatt,2,1,1
allsel
esize,0.01
mshk,2
amesh,all
!load
nsel,s,loc,x,0.1875
d,all,temp,200
nsel,s,loc,x,0.6
d,all,temp,70
nsel,all
allsel
physics,write,thermal
physics,clear
ddel,all,all
!physics struct
et,1,183,,,1
mp,ex,1,30e6
mp,alpx,1,0.65e-5
mp,nuxy,1,0.3
mp,ex,2,10.6e6
mp,alpx,2,1.35e-5
mp,nuxy,2,0.33
!load
nsel,s,loc,y,0.05
cp,1,uy,all !耦合顶部UY，使其具有共同的UY，用以模拟整个筒体的一小部分
nsel,s,loc,x,0.1875
cp,2,ux,all !耦合理由同上
nsel,s,loc,y
d,all,uy !约束目的为对称约束，可换为DSYM,SYMM,Y
nsel,all
tref,70
allsel
physics,write,struct
save
fini
/sol
physics,read,thermal
solve
fini
/post1
path,radial,2,30,100
ppath,1,,0.1875
ppath,2,,0.6
pasave,radial,filea
fini
/sol
physics,read,struct
ldread,temp,,,,,,rth
solve
fini
/post1
paresu,radial,filea
pmap,,mat
pdef,sx,s,x
pdef,sz,s,z
plpath,sx,sz
!plpagm,sx,,node
!plnsol,s,eqv,0,1
!/expand,27,axis,,,10
!plnsol,s,eqv,0,1
5、沿径向的径向应力SX，周向应力SZ的分布曲线 plpath,sx,sz
6、沿径向的径向应力SX图形显示 plpagm,sx,,node
7、等效应力分布
8、等效应力分布3d扩展
