我也在用ansys建模,小程序:
Ansys to Flac3d.for内容:
real node(3,5000000)
integer element(8,5000000)
open(5, file='node.dat',status='old')
open(3, file='element.dat',status='old')
open(2, file='outforflac3d.dat')
write(2,*) 'title'
write(2,*) ' ansys to flac3d '
write(2,*) 'plo crea qqq '
write(2,*) 'plo add surf white'
write(2,*) 'plo add axes green '
write(2,*) 'plo set rot 0 0 35 '
write(2,*) 'plo set cen 0 10 0 '
write(2,*) 'plo set mag 0.64 '
write(2,*) 'plo set dis 26.59 '
write(2,*) 'plot show '
write(*,*) 'please input number of nodes '
read(*,*) numnode
write(*,*) 'please input number of elements '
read(*,*) numelement
! 读入单元结点信息
do 10 i=1,numnode
read(5,*) n,node(1,i),node(2,i),node(3,i)
10 continue
do 20 i=1,numelement
read(3,*) element(1,i),element(2,i),element(3,i),element(4,i),element(5,i),&
element(6,i),element(7,i),element(8,i),x1,x2,x3,x4,x5,i
20 continue
! 输出为Flac3D所要求的单元结点信息
do 30 i=1,numnode
write(2,994) i,node(1,i),node(3,i),node(2,i)
994 format(1x,'gen point id ',i5,' (',f8.4,',',f8.4,',',f8.4,')')
30 continue
do 40 i=1,numelement
write(2,996) element(4,i),element(1,i),element(8,i),element(3,i),element(5,i),&
element(7,i),element(2,i),element(6,i),i
996 format(1x,'gen zone brick p0 point',i5,' p1 point',i5,' p2 point',i5,' p3 point',i5,&
&' p4 point',i5,' p5 point',i5,' p6 point',i5,' p7 point',i5,' size 1 1 1 rat 1 1 1 &
& group ', i5)
40 continue
write(2,*) '; gen zon reflect dip 0 dd 0 origin 0 0 0 '
end |