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

時間:2008-09-15 11:30:04 來源:

這篇文章的目的是利用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é)點坐標和單元所包含的節(jié)點編號，以及單元編號集節(jié)點編號集等，另外一部分就是利用第一部分的單元節(jié)點集來設(shè)置一些加載邊界條件，及設(shè)置材料屬性，求解及輸出控制等。
     inp文件的第二部分用手寫就可以實現(xiàn)，而第一部分比較難實現(xiàn)（模型復(fù)雜就更不用說了）。
     首先是把ansys模型中所有的節(jié)點坐標讀到數(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é)點坐標和單元所包含的節(jié)點，節(jié)點集和單元集都需要另想辦法實現(xiàn)。#p#分頁標題#e#
       我總算找到了從ansys前處理到abaqus求解器的一條小路（作用有限）^_^ ^_^

上一篇： ANSYS中用表面效應(yīng)單元加任意方向的荷載

下一篇：任意位置繪制橢圓兩法ZZ

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

相關(guān)文章