命令流共享-FLAC3D模型导入ANSYS中.doc

;* get nodetotal & zone total *; Initializationdef initializationIO_READ = 0 IO_WRITE = 1IO_FISH = 0IO_ASCII = 1SuferFile=FLAC-ANSYS_MODEL.datarray message(1)endinitializationdef 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)enddefine 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 GetNodesmessage(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) endloopEND; Main Functiondef haitang status = closestatus = 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)=mess