求邓肯－张模型的源程序

风零2

那编译好的你下载了吗？
我只看到过ADINA的，没看到ansys的
现在没办法，只有试着用宏做，还老是出错
要有。dll直接调用就好了

mushengyuan

就是啊，有源代码吗？供大家参考一下啊！

1360

我也在进行邓肯-张模型的二次开发，但不成功。下面是二维问题的邓肯-张源程序，主要根据应力计算土体的弹性模量和伯松比。


C
C     M A T E R I A L   P R O P E R T I E S
C
C            土的参数
C
C
              SSK=CTI(1)               ! K
              SSN=CTI(2)               ! n
              SRF=CTI(3)               ! Rf
              SSG=CTI(4)               ! G
              SSF=CTI(5)               ! F
              SSD=CTI(6)               ! D
              SKUR=CTI(7)              ! Kur
              SSC=CTI(8)               ! C
              SFAI=CTI(9)              ! 摩擦角(°)
              SFAI=SFAI*ATAN(1.0)/45.0
              PA=CTI(11)               ! 大气压力
              IF(PA.EQ.0.0) PA=9.801
C
C         应力分量和主应力
C
              STSX=STRESS(1)               !Sigma-x
              STSY=STRESS(2)               !Sigma-y
              STSS=STRESS(3)               !Tao-xy
              STS0=(STSX+STSY)/2.0
              SDIF=(STSX-STSY)/2.0
              STSR=SQRT(SDIF*SDIF+STSS*STSS)
              STS1=STS0-STSR               !Sigma-1
              STS3=STS0+STSR               !Sigma-3
C
C      计算切线模量和泊松比
C              
           IF(STS1.EQ.0.0) THEN            !第一次计算时，由于没有应力值，用切线值代替。
              YM=SSK
              PR=SSG
           ELSE
              A1=1.0-SIN(SFAI)
              BD=STS3-STS1                 !本来应该为Sigma1-Sigma3,但由于ADINA的应力符号规定与土力学不一致
              AC=-STS3
              IF(AC.LT.0.1) AC=0.1
              UU=A1*BD
              DD=2*SSC*COS(SFAI)+2.0*AC*SIN(SFAI)
              SRATIO=UU/DD                            ! 应力水平
              IF(SRATIO.GE.1.0) SRATIO=1.0
              B1=SSK*(ABS(AC/PA))**SSN
              YM=B1*(1.0-SRATIO*SRF)**2
              IF(TIME.EQ.0.0) YM=SSK
              C1=SSG-SSF*LOG10(ABS(AC/PA))           ! 初始波松比
              IF(C1.GT.0.48) C=0.48
              IF(C1.LT.0.01) C=0.01
              PR=C1+(0.48-C1)*SRATIO
              IF(TIME.EQ.0.0) PR=0.3
           END IF


希望开发成功以后能够共享。