[转帖]包含接触的热－力直接耦合例子

linda

搞到这个例子，感觉还不错！如果重复了，请斑竹删除哦！题目如下：
两个同心圆管之间有一个小间隙，内管中突然流入一种热流体，求经过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 !查看等效应力分布

[ 本帖最后由 rodge 于 2007-9-23 15:35 编辑 ]

yagn

gcgen,cont,targ,3 !创建接触单元？
这句是不是有点问题啊

cuanxinjia

好贴。我下下来学习了。呵呵呵。