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

在7.1版的產品新功能說明中，描述如下：

- 內部MPC算法可以用于接觸單元

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

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

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

內部MPC 方法可以克服傳統接觸算法和ANSYS中的其它多點約束工具的缺點，例如：

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

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

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


終于測試成功，原來在使用SHSD前，使用了nummrg,kp命令，出現錯誤。計算最大位移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

此主題相關圖片如下：



全部使用實體單元，計算結果很接近，如下

/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#分頁標題#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

此主題相關圖片如下：