分析例子为<挢旗扦渌;ansys13.0热力学有限元分析从入门到精通>15.2的例子：两厚璧筒热应力分析
需要注意的是单元的匹配，比如用77单元采用8节点网格，那么结构单元也须采用8节点单元183。否则会出现错误或警告。
apdl运行通过，单元网格尺寸改为0.01
apdl通过运行，结果正确无警告。
!*Thermal-StructuralCouplingAnalysis
fini
/cle
/filname,exce,1
/title,lovz
/prep7
!physicsthermal
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
!physicsstruct
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
沿径向的径向应力SX，周向应力SZ的分布曲线plpath,sx,sz
沿径向的径向应力SX图形显示plpagm,sx,,node
等效应力分布
等效应力分布3d扩展
