Ansys中，SHSD用于殼-實(shí)體裝配，我試了好長(zhǎng)時(shí)間，也沒有成功。不知哪位兄弟有算例或好的方法？

在7.1版的產(chǎn)品新功能說明中，描述如下：

- 內(nèi)部MPC算法可以用于接觸單元

一個(gè)新的內(nèi)部多點(diǎn)約束(MPC) 算法可以用于下列接觸單元：CONTA171, CONTA172, CONTA173, CONTA174, 和 CONTA175。利用該方法，程序可以基于接觸運(yùn)動(dòng)關(guān)系創(chuàng)建內(nèi)部MPC方程。

你可以使用MPC方法(KEYOPT(2) = 2)結(jié)合綁定接觸定義(KEYOPT(12) = 5或6)來定義下列類型的接觸裝配和運(yùn)動(dòng)約束關(guān)系：

%26#8226;實(shí)體-實(shí)體裝配 – 接觸面與目標(biāo)面定義在實(shí)體單元表面
%26#8226;殼-殼裝配 – 接觸面與目標(biāo)面定義在殼單元面
%26#8226;殼-實(shí)體裝配 – 接觸面定義在殼單元面上，目標(biāo)面定義在實(shí)體單元面上
%26#8226;剛性約束面 –接觸節(jié)點(diǎn)被約束為由控制節(jié)點(diǎn)定義的剛體位移 (類似于CERIG命令)
%26#8226;分配力表面 –施加在控制節(jié)點(diǎn)上的力或位移以一種基于形函數(shù)意義上的平均分配到接觸節(jié)點(diǎn)上 (類似于RBE3命令)
%26#8226;梁-實(shí)體裝配 – 一根梁的端部節(jié)點(diǎn)作為控制節(jié)點(diǎn)并與實(shí)體或殼面相連(使用剛性約束面或者分配力表面類型的MPC)

內(nèi)部MPC 方法可以克服傳統(tǒng)接觸算法和ANSYS中的其它多點(diǎn)約束工具的缺點(diǎn)，例如：

%26#8226;接觸面節(jié)點(diǎn)的自由度被消去，可以減小系統(tǒng)方程求解的波前大小
%26#8226;不需要輸入接觸剛度。對(duì)于小變形問題，它代表真實(shí)的線性接觸行為；求解系統(tǒng)方程時(shí)不需要平衡迭代。對(duì)于大變形問題，MPC方程在每個(gè)平衡迭代中不斷進(jìn)行更新，克服了傳統(tǒng)約束方程只適用于小應(yīng)變的限制條件
%26#8226;平動(dòng)和轉(zhuǎn)動(dòng)自由度都可以約束
%26#8226;內(nèi)部MPC的生成非常簡(jiǎn)單，因?yàn)椴捎玫氖墙佑|對(duì)定義方法
%26#8226;自動(dòng)考慮形函數(shù)，對(duì)于分配力多點(diǎn)約束，如果采用了高階單元或軸對(duì)稱單元，不需要加權(quán)因子(類似于RBE3命令)。除了力以外，這種類型的MPC還可以在控制節(jié)點(diǎn)上施加位移。

一個(gè)新的網(wǎng)格工具， SHSD命令，可幫助您生成用于殼－實(shí)體裝配的殼－實(shí)體界面。

詳細(xì)信息參見ANSYS Structural Analysis Guide中的Using the Internal MPC Approach for Assemblies and Kinematic Constraints 。


終于測(cè)試成功，原來在使用SHSD前，使用了nummrg,kp命令，出現(xiàn)錯(cuò)誤。計(jì)算最大位移0.002，命令流如下：

/prep7
et,1,45
et,2,63
R,2,1,1,1,1,0,0,
mp,ex,1,1e7
mp,prxy,1,.2
block,0,5,-0.5,0.5,-0.5,0.5
wprot,,-90
rect,5,10,-.5,.5 ! area at center of two blocks
esize,.25
vmesh,all
type,2
real,2
mshape,2
esize,.5
aslv,u
amesh,all
et,3,TARGE170
keyopt,3,5,2
et,4,Conta175
KEYOPT,4,2,2
KEYOPT,4,12,5
keyopt,4,5,0
KEYOPT,4,11,1
R,3
nslv,s,1
nsel,r,loc,x,5
TYPE,3
real,3
ESURF
ESEL,S,TYPE,,2
nsle,s
nsel,r,loc,x,5
TYPE,4
real,3
ESURF
save

alls
/PSYMB,ESYS,1
/view,1,1,1,1
eplot

shsd,3,CREATE

/PSF,PRES,NORM,2,0,1
nsel,s,loc,x,10
esln
nsle
sf,all,pres,-10
alls
nsel,s,loc,x
d,all,all
nsel,all
save,model,db
fini
/solu
solve
fini
/post1
/contour,,12
plnsol,u,sum

此主題相關(guān)圖片如下：



全部使用實(shí)體單元，計(jì)算結(jié)果很接近，如下

/prep7
et,1,45
mp,ex,1,1e7
mp,prxy,1,.2
block,,10,-0.5,.5,-.5,.5
esize,.25
vmesh,all
nummrg,all
numcmp,all
nsel,s,loc,x,10-0.5,10
nsel,r,loc,y,0.5 #p#分頁(yè)標(biāo)題#e#
/PSF,PRES,NORM,2,0,1
sf,all,pres,10
alls
nsel,s,loc,x
d,all,all
nsel,all
save,model,db
fini
/solu
solve
fini
save

此主題相關(guān)圖片如下：