CAE聯(lián)盟( http://forums.caenet.cn,) NaCl原創(chuàng),轉(zhuǎn)載請注明出處.
論壇現(xiàn)在傳不上圖片,過幾天再試.抱歉.不過應(yīng)該也可以看懂的.
關(guān)于MARC中子結(jié)構(gòu)的應(yīng)用,幫助文檔里有三個例子,比較典型的是齒輪分析和兩個圓柱體接觸分析的
例子,齒輪的例子中指出,完整模型求解為1002S,利用子結(jié)構(gòu)方法求解,總時間為264S,效率得高近
3/4,可見效果較好,本人最近想利用子結(jié)構(gòu)解決問題,苦于沒有參考書籍及文獻,文獻中多以ANSYS的
子結(jié)構(gòu)進行分析,而MARC的幫助文檔中沒有一個完整的操作設(shè)置過程,在網(wǎng)上搜到的也都是求助,著
實為難了一番,研究數(shù)日,小有收獲,與大家分享.費話不多說,上例子。
參照兩個圓柱接觸的例子，做了兩個長方體接觸的模型，見圖1，下端固定,上端X向固定,Y向施加向下位移,
其它完整模型分析不詳述，求解共計32S。
子結(jié)構(gòu)求解步驟：
1，將每個變形體分為兩部分，分別為子結(jié)構(gòu)部分和超單元部分，如圖2。
2，求解超單元部分的剛度矩陣。將完整模型文件另存一個,如sub1,只保留上部變形體的超單元部
分，其余單元刪除，在超單元上邊界和下邊界（與子結(jié)構(gòu)相連的節(jié)點）添加邊界條件DOF-SET
NODES，固定其X，Y兩個方向的位移。命令：boundary conditions-〉structural-〉more-〉dof-
set nodes-〉displacement x/y。
3，超單元設(shè)置。在loadcase中，將類型設(shè)置為superelement。命令：loadcase-〉type-〉
superelement/
在proerties中，選中dof-set nodes，并在matrix name中設(shè)置一個矩陣名稱（很重要，后面要用
到）。
4，設(shè)置DMIG輸出，在JOBS中的DMIG OUT菜單中，選中stiffness matrix，在彈出菜單中，依次選
中element stiffness/global stiffness 及frequency中的once，coord system中的global，
elements中的all elements。
5，保存模型進行求解，計算時間約5S，正常退出號為3022。
6,相同的方法,把完整模型另存sub2,將下部變形體的超單元部分剛度矩陣求出,注意matrix name中
設(shè)置一個矩陣名稱要換成不同于上部的,在文件夾中會分別生成兩個剛度矩陣的存儲文件
XXX_XXX_dmigst_0000
7,兩個超單元的剛度矩陣求出后,就開始計算關(guān)鍵部位了.再把完整模型另存為sub,將定義為超單元
的兩個部分刪掉,注意將帶有邊界條件的節(jié)點保留.
8.在JOBS菜單中,進入include file菜單,在下面的model section欄里選add file,分別將上面生成
的兩個剛度矩陣存儲文件導(dǎo)入.OK確認(rèn).
9.進入input file text菜單，依次輸入如下4行命令：
k2gg，XXX（第一個剛度矩陣名稱），1，
include，XXX_XXX_dmigst_0000（第一個剛度矩陣存儲文件名稱）
k2gg，XXX（第二個剛度矩陣名稱），1，
include，XXX_XXX_dmigst_0000（第二個剛度矩陣存儲文件名稱）
10.接下來就可以求解了,時間約為5S,這樣求解剛度矩陣各5S,主程序計算5S,共計約15S,節(jié)省了50%
以上的計算時間,模型越大,效果越名顯,可以看以一下結(jié)果對比圖.