兩個(gè)同心圓管之間有一個(gè)小間隙，內(nèi)管中突然流入一種熱流體，求經(jīng)過(guò)3分鐘后外表面的溫度分布和應(yīng)力分布，已知條件：
管材彈性模量：2E11N/(m*m)
熱膨脹系數(shù)：5E-5 /F
泊松比：0.3
導(dǎo)熱系數(shù)：10W/mC
密度：7880Kg/(m*m)
比熱：500 J/Kg*C
外管外半徑：0.131m
外管內(nèi)半徑：0.121m
內(nèi)管外半徑：0.12m
內(nèi)管內(nèi)半徑：0.11m
流體溫度：300 C
流體與內(nèi)管內(nèi)壁對(duì)流系數(shù)：300W/m*m*C
內(nèi)、外管接觸導(dǎo)熱：0.1W/C （其中C表示攝氏度）
命令流如下：
/prep7
et,1,13,4,,1 !直接選用耦合單元plane13，自由度為ux、uy、temp，軸對(duì)稱(chēng)
et,2,48 !定義結(jié)構(gòu)接觸單元
keyopt,2,1,1 !設(shè)定接觸的相應(yīng)選項(xiàng)
keyopt,2,2,1
keyopt,2,7,1
r,2,2e11,0,0.0001,,,0.1 !定義接觸單元的實(shí)常數(shù)
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 !耦合內(nèi)管內(nèi)壁的X位移和溫度自由度
nsel,s,loc,x,0.12
cp,3,ux,all
cp,4,temp,all !耦合內(nèi)管外壁的X位移和溫度自由度
nsel,s,loc,x,0.121
cp,5,ux,all
cp,6,temp,all !耦合外管內(nèi)壁的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 !耦合內(nèi)管頂部的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 !創(chuàng)建接觸單元
/solu
antype,trans !進(jìn)行瞬態(tài)分析
tunif,20 !定義初始平均溫度
tref,20 !定義參考溫度
sfl,4,conv,300,,300 !內(nèi)管內(nèi)壁對(duì)流邊界條件
sfl,6,conv,10,,20 !外管外壁對(duì)流邊界條件
nsel,s,loc,y,0 !約束底邊的Y向位移
d,all,uy,0
time,180 !定義分析時(shí)間3分鐘
deltime,10,5,15 !時(shí)間步長(zhǎng)
outres,all,all
kbc,1
autots,on !打開(kāi)自動(dòng)時(shí)間步長(zhǎng)
allsel
solve
/post1
plnsol,temp !查看溫度分布
plnsol,s,eqv !查看等效應(yīng)力分布


48用于點(diǎn)面接觸，無(wú)需預(yù)先知道確切的接觸位置，接觸面之間也不需保持一致的網(wǎng)格。
171和169用于面面接觸，與點(diǎn)－面接觸單元比，具有很多優(yōu)點(diǎn)，例如：需要較少的接觸單元，因而只需較少的磁盤(pán)空間和CPU時(shí)間，有利于可視化...
但是如果通過(guò)一組節(jié)點(diǎn)來(lái)定義接觸面，生成多個(gè)單元，那么就可以通過(guò)點(diǎn)－面接觸單元來(lái)模擬面－面的接觸問(wèn)題