复杂的群桩习题，尚存在一些问题，欢迎指点！

huahua2005

最近做一个群桩的练习题（自己结合以往工程做的练习），考虑了所有的桩（几十根）及其上建筑物、附近一部分填土区域，桩与土之间、建筑物与土之间、土层与土层之间均设置接触（感觉ADINA-m建模，加上用fortran写出建群体的命令流与有规律的群接触面，还是很爽的，比flac要好用很多），最终网格是划起来了（划分时基本采用了4面体，边界为advacing front的划分方式，划不出来的地方如body88采用边界为delaunay的划分方式)，也可以算，但一直不收敛（换了n种point-size，一直不收敛，这点上就比flac差些了），现将in文件发上来(由于网格划分时间较长，将划分前的命令发上来，划分后即可算），期待高手高见，也欢迎大家切磋～～

mingzzeng

请问，关于 fortran和adina之间开发的资料有哪些？
能否共享？？

huahua2005

不是开发资料，是自己编写的生成命令流的程序，很简单的～～
给个生成群桩的例子如下：
PROGRAM MAIN
         DIMENSION DIAN(6,7,2),POINT(3,10,2)
        OPEN(1,FILE='DATA.IN')
        XI=0.5
        YI=0.5
        XD=3.0
        YD=3.0
      IC=5
        IW=6
        RC=0.25
         DD=-3.8
         DG=-14
         CC=DD-DG
         ZC=(DD+DG)/2.0

         DD1=-5.1
         DG1=-14
         CC1=DD1-DG1
         ZC1=(DD1+DG1)/2.0

         KK=0
         DO 1,I=1,IC+1
         DO 1,J=1,IW+1
       DIAN(I,J,1)=XI+(I-1)*XD
         DIAN(I,J,2)=YI+(J-1)*YD
                   KK=KK+1
        WRITE(1,'(A,I2,2(A,F5.2),A,/3(A,F6.2),A)')'BODY CYLINDER NAME=',KK       
     $,' OPTION=CENTERED POSITION=VECTOR ORIENTAT=SYSTEM CX1=',DIAN(I,J,
     $1),' CX2=',DIAN(I,J,2),',',' CX3=',ZC,' SYSTEM=0 AXIS=ZL RADIUS=',
     $RC,' LENGTH=',CC,' SHEET=NO'
1     CONTINUE
      DO 2,I=1,7
          KK=KK+1
        AX=DIAN(6,7,1)+1.5
        BY=YI+(I-1)*YD
        WRITE(1,'(A,I2,2(A,F5.2),A,/3(A,F6.2),A)')'BODY CYLINDER NAME=',KK       
     $,' OPTION=CENTERED POSITION=VECTOR ORIENTAT=SYSTEM CX1=',AX,' CX2=
     $',BY,',',' CX3=',ZC1,' SYSTEM=0 AXIS=ZL RADIUS='
     $,RC,' LENGTH=',CC1,' SHEET=NO'
2    CONTINUE
         DO 3,I=1,3
         DO 3,J=1,10
       POINT(I,J,1)=AX+1.92+(I-1)*3.8
         POINT(I,J,2)=YI+(J-1)*2
                   KK=KK+1
        WRITE(1,'(A,I2,2(A,F5.2),A,/3(A,F6.2),A)')'BODY CYLINDER NAME=',KK       
     $,' OPTION=CENTERED POSITION=VECTOR ORIENTAT=SYSTEM CX1=',POINT(I,J
     $,1),' CX2=',POINT(I,J,2),',',' CX3=',ZC1,' SYSTEM=0 AXIS=ZL RADIUS
     $=',RC,' LENGTH=',CC1,' SHEET=NO'
3     CONTINUE
         END

武汉游侠

楼主的30万个单元的确很考验机子
我的机子承受不起
不过好像模型里面body80的face5没有约束
呵呵

huahua2005

后来加上了，也还是不行，哎，仍然谢谢楼上的，呵呵～