包含接触的热－力直接耦合例子

两个同心圆管之间有一个小间隙，内管中突然流入一种热流体，求经过3分钟后外表面的温度分布和应力分布，已知条件：
管材弹性模量：2E11N/(m*m)
热膨胀系数：5E-5 /F
泊松比：0.3
导热系数：10W/mC
密度：7880Kg/(m*m)
比热：500 J/Kg*C
外管外半径：0.131m
外管内半径：0.121m
内管外半径：0.12m
内管内半径：0.11m
流体温度：300 C
流体与内管内壁对流系数：300W/m*m*C
内、外管接触导热：0.1W/C （其中C表示摄氏度）
命令流如下：
/prep7
et,1,13,4,,1 !直接选用耦合单元plane13，自由度为ux、uy、temp，轴对称
et,2,48 !定义结构接触单元
keyopt,2,1,1 !设定接触的相应选项
keyopt,2,2,1
keyopt,2,7,1
r,2,2e11,0,0.0001,,,0.1 !定义接触单元的实常数
mp,ex,1,2e11 !定义管材的属性
mp,alpx,1,5e-5
mp,prxy,1,0.3
mp,kxx,1,10
mp,dens,1,7880
mp,c,1,500
rect,0.11,0.12,0,0.02 !建模
rect,0.121,0.131,0,0.02
amesh,all
nsel,s,loc,x,0.11
cp,1,ux,all
cp,2,temp,all !耦合内管内壁的X位移和温度自由度
nsel,s,loc,x,0.12
cp,3,ux,all
cp,4,temp,all !耦合内管外壁的X位移和温度自由度
nsel,s,loc,x,0.121
cp,5,ux,all
cp,6,temp,all !耦合外管内壁的X位移和温度自由度
nsel,s,loc,x,0.131
cp,7,ux,all
cp,8,temp,all !耦合外管外壁的X位移和温度自由度
nsel,s,loc,y,0.02
nsel,r,loc,x,0,0.12
cp,9,uy,all !耦合内管顶部的Y位移
nsel,s,loc,y,0.02
nsel,r,loc,x,0.121,0.131
cp,10,uy,all !耦合外管顶部的Y位移
nsel,s,loc,x,0.12
cm,cont,node
nsel,s,loc,x,0.121
cm,targ,node
type,2
real,2
gcgen,cont,targ,3 !创建接触单元
/solu
antype,trans !进行瞬态分析
tunif,20 !定义初始平均温度
tref,20 !定义参考温度
sfl,4,conv,300,,300 !内管内壁对流边界条件
sfl,6,conv,10,,20 !外管外壁对流边界条件
nsel,s,loc,y,0 !约束底边的Y向位移
d,all,uy,0
time,180 !定义分析时间3分钟
deltime,10,5,15 !时间步长
outres,all,all
kbc,1
autots,on !打开自动时间步长
allsel
solve
/post1
plnsol,temp !查看温度分布
plnsol,s,eqv !查看等效应力分布


48用于点面接触，无需预先知道确切的接触位置，接触面之间也不需保持一致的网格。
171和169用于面面接触，与点－面接触单元比，具有很多优点，例如：需要较少的接触单元，因而只需较少的磁盘空间和CPU时间，有利于可视化...
但是如果通过一组节点来定义接触面，生成多个单元，那么就可以通过点－面接触单元来模拟面－面的接触问题