這篇文章的目的是利用ansys強大的apdl命令編程建模，然后利用abaqus豐富的材料屬性，本構(gòu)關(guān)系。這樣就充分利用了兩個軟件的特長，而workbench在這里面很好的起到了一個橋梁作用，建議大家直接看方法二，簡單有效

　　有興趣的也可以看下方法一，over!

　　這段時間在做一個關(guān)于材料隨機分配的問題，發(fā)現(xiàn)ansys的apdl能實現(xiàn)材料屬性隨機分配，但是要用到abaqus的非線性求解器，但是ansys和abaqus之間不知道怎么轉(zhuǎn)換，所以想了個笨辦法：

　　abaqus中inp文件主要有兩部分內(nèi)容，一部分是節(jié)點坐標(biāo)和單元所包含的節(jié)點編號，以及單元編號集節(jié)點編號集等，另外一部分就是利用第一部分的單元節(jié)點集來設(shè)置一些加載邊界條件，及設(shè)置材料屬性，求解及輸出控制等。

　　inp文件的第二部分用手寫就可以實現(xiàn)，而第一部分比較難實現(xiàn)(模型復(fù)雜就更不用說了)。

　　首先是把ansys模型中所有的節(jié)點坐標(biāo)讀到數(shù)組里，再把每個單元中的節(jié)點編號輸?shù)揭粋€二維數(shù)組中(NELEM命令)，最后通過*vwrite命令輸出到文本文件中。

　　節(jié)點集和單元集的實現(xiàn)：

　　ansys中有nsel，esel這兩個命令很強大，幾乎能選擇你想要的節(jié)點集和單元集，選擇好集合之后導(dǎo)入到一個數(shù)組中，通過*vwrite輸出到文本文件中。

　　例如下面是選擇材料材料屬性為1的單元集，最后把這些單元都輸出到數(shù)組emat1(namt1)中了。

　　esel,s,mat,,1

　　*GET,nmat1,elem,0,COUNT

　　*DIM,emat1,ARRAY,nmat1

　　*GET,emat_min,elem,0,NUM,MIN

　　*DO,i,1,nmat1,1

　　emat1(i)=emat_min

　　emat_min=ELNEXT(emat_min)

　　*ENDDO

　　其中要注意的是abaqus中集合每行不能超過16個數(shù)據(jù)，因此通過以上辦法得到的文本文件須處理成abaqus要求的格式，可以用fortran或matlab實現(xiàn)(我是利用fortran的格式輸入輸出實現(xiàn)的)。

　　另外一個辦法是用Workbench FE model可以直接輸出abaqus inp文件，但是這種辦法只能輸出模型的節(jié)點坐標(biāo)和單元所包含的節(jié)點，節(jié)點集和單元集都需要另想辦法實現(xiàn)。