ansys-flac

ansys-flac 利用输入数据Node.dat及Element.dat产生Flac3D所需要的文件outforflac3d.dat

!ANSYS转FLAC3D程序，ANSYS中为SOLID45单元
!能转四面体，锲形块，和八面体
!group是以体来分的。
!by sun

nsel,all
esel,all
*get,NodeNum,node,,NUM,MAX
*get,EleNum,elem,,NUM,MAX
*dim,NodeData,array,NodeNum,3
*dim,EleData,table,EleNum,8 
 

*do,i,1,NodeNum,1 
  *get,NodeData(i,1),node,i,LOC,x 
  *get,NodeData(i,2),node,i,LOC,y 
  *get,NodeData(i,3),node,i,LOC,z 
*enddo
!写数据到文件 
*cfopen,Flac,TXT 
*vwrite, 
(';The file from ANSYS') 

!写入节点信息 

*vwrite,sequ,NodeData(1,1),NodeData(1,3),NodeData(1,2) 
g %G %G %G %G  

*do,i,1,EleNum,1 
  *get,E1,elem,i,NODE,1 
  *get,E2,elem,i,NODE,2 
  *get,E3,elem,i,NODE,3

  *get,E4,elem,i,NODE,4 
  *get,E5,elem,i,NODE,5 
  *get,E6,elem,i,NODE,6 
  *get,E7,elem,i,NODE,7 
  *get,E8,elem,i,NODE,8
!判断节点3，4的是否相同
  *if, E3,EQ,E4,then
   		*if, E5,EQ,E6, then
   		  	!写入单元信息 四面体
*vwrite,i,E1,E2,E3,E5
z t4 %G %G %G %G %G  size 1 1 1 rat 1 1 1 
  	   	*else
  	    		!写入单元信息 五面体
*vwrite,i,E1,E2,E5,E3,E6,E7
z w6 %G %G %G %G %G %G %G  size 1 1 1 rat 1 1 1
		*endif
	*else
     		!写入单元信息 八面体
*vwrite,i,E1,E2,E4,E5,E3,E8,E6,E7 
z b8 %G %G %G %G %G %G %G %G %G size 1 1 1 rat 1 1 1 
    *endif
*enddo


!分群，以体为基准，一个体分为一个群
!选体
*get,VolumeNum,volu,,NUM,MAX
*do,i,1,VolumeNum,1
  VSEL,S,,,i
  ESLV,R
  *get,EleNum_min,elem,,NUM,MIN
  *get,EleNum_max,elem,,NUM,MAX
  *vwrite
  
*GROUPS
  *vwrite,i
ZGROUP %G
	*do,j,EleNum_min,EleNum_max,10
	  *if,j+9,eq,EleNum_max,then
        *vwrite,j,j+1,j+2,j+3,j+4,j+5,j+6,j+7,j+8,j+9
%G %G %G %G %G %G %G %G %G %G %G
	  		*elseif,j+8,eq,EleNum_max,then
	  	  *vwrite,j,j+1,j+2,j+3,j+4,j+5,j+6,j+7,j+8
%G %G %G %G %G %G %G %G %G %G
	  		*elseif,j+7,eq,EleNum_max,then
	  		*vwrite,j,j+1,j+2,j+3,j+4,j+5,j+6,j+7
%G %G %G %G %G %G %G %G %G
	  		*elseif,j+6,eq,EleNum_max,then
	  		*vwrite,j,j+1,j+2,j+3,j+4,j+5,j+6
%G %G %G %G %G %G %G %G
	  		*elseif,j+5,eq,EleNum_max,then
	  		*vwrite,j,j+1,j+2,j+3,j+4,j+5
%G %G %G %G %G %G %G
	  		*elseif,j+4,eq,EleNum_max,then
	  		*vwrite,j,j+1,j+2,j+3,j+4
%G %G %G %G %G %G
	  		*elseif,j+3,eq,EleNum_max,then
	  		*vwrite,j,j+1,j+2,j+3
%G %G %G %G
	  		*elseif,j+2,eq,EleNum_max,then
	  		*vwrite,j,j+1,j+2
%G %G %G 
	  		*elseif,j+1,eq,EleNum_max,then
	  		*vwrite,j,j+1
%G %G 
	  	*else
				*vwrite,j,j+1,j+2,j+3,j+4,j+5,j+6,j+7,j+8,j+9
%G %G %G %G %G %G %G %G %G %G %G
			*endif
   *enddo
   allsel
*ENDDO
     
*cfclos