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

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

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

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

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

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

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

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

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

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


終于測試成功，原來在使用SHSD前，使用了nummrg,kp命令，出現(xiàn)錯誤。計算最大位移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)圖片如下：



全部使用實體單元，計算結(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#分頁標(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)圖片如下：