[原创] 命令流共享－FLAC3D模型导入ANSYS中

haitang

在学习了版主dynamax 的《用tecplot对flac3d计算结果进行后处理》后，借鉴了用一个字符串进行文件输入输出的思想，编制了以下的命令流，由于flac3d中不好剖分网格，所以将FLAC3D模型导入ANSYS中进行加密处理，再利用我前面编得从ansys转到flac3d的程序，输入到flac3d中进行计算，会是一个比较好的办法，可以供大家学习。

;********************** get nodetotal & zone total *********************

;; Initialization
def initialization
        IO_READ  = 0
        IO_WRITE = 1
        IO_FISH  = 0
        IO_ASCII = 1
        SuferFile='FLAC-ANSYS_MODEL.dat'
        array message(1)
end
initialization

def writehead
    pnt1=zone_head
    pnt2=gp_head
    Z_total=0
    P_total=0
    loop while pnt1 # null
    Z_total=Z_total+1
    pnt1=z_next(pnt1)
    endloop

    loop while pnt2 # null
    P_total=P_total+1
    pnt2=gp_next(pnt2)
    endloop
    message(1)='!------------------------------------------------------------------ \n'
    message(1)='!This is a free soft for converting the model from Flac3d to Ansys! \n'
    message(1)='!                Code by haitang 2006-05-30 \n'
    message(1)='!------------------------------------------------------------------ \n'
    message(1)='!The Number of Nodes and Elements \n'
    message(1)=message(1)+'!'+string(P_total)+' , '+string(Z_total)+'\n'
    message(1)=message(1)+'/prep7 \n'
    message(1)=message(1)+'ET,1,SOLID45'
    status = write(message,1)
end

define GetElements
   Z1_p=zone_head

   loop n (1,Z_total)
     Z1_code=z_code(Z1_p)
     Z1_No=z_id(Z1_p)
     Z1_Group=z_group(Z1_p)
     if Z1_code=0 then
       message(1)='MAT, '+string(Z1_Group)+'\n'   
       message(1)=message(1)+'E, '+string(gp_id(z_gp(Z1_p,1)))
       message(1)=message(1)+','+ string(gp_id(z_gp(Z1_p,2)))
       message(1)=message(1)+','+ string(gp_id(z_gp(Z1_p,5)))
       message(1)=message(1)+','+ string(gp_id(z_gp(Z1_p,3)))
       message(1)=message(1)+','+ string(gp_id(z_gp(Z1_p,4)))
       message(1)=message(1)+','+ string(gp_id(z_gp(Z1_p,7)))
       message(1)=message(1)+','+ string(gp_id(z_gp(Z1_p,8)))
       message(1)=message(1)+','+ string(gp_id(z_gp(Z1_p,6)))
       status = write(message,1)
       Z1_p=z_next(Z1_p)
     endif
     if Z1_code=1 then
       message(1)='MAT, '+string(Z1_Group)+'\n'
       message(1)=message(1)+'E, '+string(gp_id(z_gp(Z1_p,1)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,2)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,4)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,4)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,3)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,5)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,6)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,6)))
       status = write(message,1)
       Z1_p=z_next(Z1_p)
     endif
     if Z1_code=4 then
       message(1)='MAT, '+string(Z1_Group)+'\n'
       message(1)=message(1)+'E, '+string(gp_id(z_gp(Z1_p,1)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,2)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,3)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,3)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,4)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,4)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,4)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,4)))
       status = write(message,1)
       Z1_p=z_next(Z1_p)
     endif               
     if Z1_code=2 then
       message(1)='MAT, '+string(Z1_Group)+'\n'
       message(1)=message(1)+'E, '+string(gp_id(z_gp(Z1_p,1)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,2)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,5)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,3)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,4)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,4)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,4)))
       message(1)=message(1)+','+string(gp_id(z_gp(Z1_p,4)))
       status = write(message,1)
       Z1_p=z_next(Z1_p)
     endif
   endloop
END


DEFINE GetNodes
message(1)=''
    P1_p=gp_head
    loop n (1,P_total)   
        node_id=gp_id(p1_p)
        message(1)='N, '+string(node_id)
           message(1)=message(1)+','+string(gp_xpos(p1_p))
           message(1)=message(1)+','+string(gp_ypos(p1_p))
           message(1)=message(1)+','+string(gp_zpos(p1_p))
           status = write(message,1)                  
        P1_p=gp_next(P1_p)
     endloop
END

; Main Function
def haitang       
  status = close
        status = open(SuferFile,IO_WRITE,IO_ASCII)
        if status = 0 then
            writehead
                message(1)='The Node coordination of "X" "Y" "Z" \n'
            GetNodes
            message(1)='!The Nodes id nums of Elements. \n'
            GetElements
            message(1)='/PNUM,MAT,1 \n'
            message(1)=message(1)+ '/REPLOT '
            status = write(message,1)
            status = close
            ii = out('Successfully Write Data Into File ' + SuferFile)
        else
            ii=out('Open File Error! Status = ' + string(status))                       
        endif                       
end
haitang

lakewater

加分鼓励！

lbt18872

值得加分鼓励，谢谢

wangweipeng

学海无涯

zdhshan

谢谢！

lc_rw

海棠楼主,不仅学术水平高,人品也很好,感谢您的共享!
真正的高手

feiyang52013

您好啊
我现在了你的ansys2flac的程序，可是老是出问题，那个程序是不是对flac的版本，单元的是数目和形状都有要求阿

sch

兄弟很服你是有心人！可以和你多交流！

ansys2007

楼主，你以后必然成为大家！！！
绝对的大家风范！！很感激你！

xuexi2007

好像历来都是水平越高的人越是胸襟宽广 ，，知否心宽路亦宽。额是初学者，正不知道如何下手 ，莽莽撞撞的碰到了这个论坛，今天就尽情表达一下新人额的惊喜，崇拜和感激吧。

zhongnan

多谢多谢
以后还要向你请教！

smallkin

楼主是在是太牛了，佩服死了

zty20041044

haitang大哥的必顶。偶像

hmli19770604

不错，好东西，谢谢

jiangjwhl

相当受启发，FISH用的太漂亮了

xiedaolei

高手，学习了

yzhanguo

高手就是高手
学习了

sunxinpo

还是不知道怎么用啊，期待详细说明

litengfa

谢谢楼主  认真学习下