ansys模型寫abaqus的inp文件的辦法

時(shí)間:2010-11-14 12:29:45 來(lái)源:未知

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

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

這段時(shí)間在做一個(gè)關(guān)于材料隨機(jī)分配的問(wèn)題，發(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)編號(hào)，以及單元編號(hào)集節(jié)點(diǎn)編號(hào)集等，另外一部分就是利用第一部分的單元節(jié)點(diǎn)集來(lái)設(shè)置一些加載邊界條件，及設(shè)置材料屬性，求解及輸出控制等。
     inp文件的第二部分用手寫就可以實(shí)現(xiàn)，而第一部分比較難實(shí)現(xiàn)（模型復(fù)雜就更不用說(shuō)了）。
     首先是把a(bǔ)nsys模型中所有的節(jié)點(diǎn)坐標(biāo)讀到數(shù)組里，再把每個(gè)單元中的節(jié)點(diǎn)編號(hào)輸?shù)揭粋€(gè)二維數(shù)組中（NELEM命令），最后通過(guò)*vwrite命令輸出到文本文件中。
    節(jié)點(diǎn)集和單元集的實(shí)現(xiàn)：
    ansys中有nsel，esel這兩個(gè)命令很強(qiáng)大，幾乎能選擇你想要的節(jié)點(diǎn)集和單元集，選擇好集合之后導(dǎo)入到一個(gè)數(shù)組中，通過(guò)*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中集合每行不能超過(guò)16個(gè)數(shù)據(jù)，因此通過(guò)以上辦法得到的文本文件須處理成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)。
       我總算找到了從ansys前處理到abaqus求解器的一條小路（作用有限）^_^ ^_^

上一篇： ANSYS 11.0安裝三法探討

下一篇：AutoCAD模型輸入Ansys

ansys模型寫abaqus的inp文件的辦法

相關(guān)文章