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

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

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

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

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

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

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

　　ansys中有nsel，esel這兩個(gè)命令很強(qiáng)大，幾乎能選擇你想要的節(jié)點(diǎn)集和單元集，選擇好集合之后導(dǎo)入到一個(gè)數(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個(gè)數(shù)據(jù)，因此通過以上辦法得到的文本文件須處理成abaqus要求的格式，可以用fortran或matlab實(shí)現(xiàn)(我是利用fortran的格式輸入輸出實(shí)現(xiàn)的)。

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