- 积分
- 0
- 注册时间
- 2014-2-19
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 dongdong2011 于 2015-8-10 14:52 编辑
转载【转载】命令流共享-FLAC3D模型导入ANSYS中
base on Code by haitang 2006-05-30
使用方法:先使用flac读取网格文件FlacGrid.Flac3D,然后使用ANSYS文件导入功能 read input from 读取文件FLAC_to_ANSYS_MODEL.MAC
代码下载请回复
个人修改了一点点的代码
仅为共享
- <span style="background-color: rgb(255, 255, 255);">代码下载请回复</span>
复制代码 部分源代码
- ;********************** get nodetotal & zone total *********************
- ;使用方法:先使用flac读取网格文件FlacGrid.Flac3D,
- ;然后使用ANSYS文件导入功能 read input from 读取文件FLAC_ANSYS_MODEL.MAC
- ;; Initialization
- def initialization
- IO_READ = 0
- IO_WRITE = 1
- IO_FISH = 0
- IO_ASCII = 1
- SuferFile='FLAC-to-ANSYS_MODEL.MAC'
- 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)='! edit by dongdong2011 2015-08-10 \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
-
复制代码
haitang的源代码
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|