新手请教如何调用ABAQUS内部函数？

whua8

以下是论坛上的邓肯张模型的一小部分，其中CALL SPRINC(STRESS,PS,1,NDI,NSHR)是调用主应力函数吧，请问STRESS,PS,1,NDI,NSHR变量的书写是不是一定的啊 ，那其中的“1”代表什么意思呢？小弟刚学请指教

        S1S3O=STATEV(1)
        S3O=STATEV(2)
        SSO=STATEV(3)
        CALL SPRINC(STRESS,PS,1,NDI,NSHR)
        write(*,*)"PS(1)=",PS(1)

zhangmeng

你可以在帮助中找到

whua8

下面这段UMAT是再论坛上下载的邓肯张模型UMAT，
其中的语句S1S3O=STATEV(1)  
        S3O=STATEV(2)
        SSO=STATEV(3)
C 程序中给出的注释为
    S1S3O为历时上最大主应力差，S3O为前期固结应力，SSO为历史上最大应力水平
楼上的兄弟请问怎么知道STATEV(1) ， STATEV(2)和STATEV(3)分别代表的是S1S3O，S3O，SSO呢。  谢谢了

  SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
     1 RPL,DDSDDT,DRPLDE,DRPLDT,
     2 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
     3 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
     4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
C
      INCLUDE 'ABA_PARAM.INC'
C
      CHARACTER*80 CMNAME
      DIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),
     1 DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),
     2 TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3),DROT(3,3),
     3 DFGRD0(3,3),DFGRD1(3,3)
C-----------------------------------------------------------------------------------
C    求解平面应变问题,存储4个应力,σ11、σ22、σ33、σ12；
C                         4个应变,ε11、ε22、ε33、γ12；
C    求解过程中在*.inp文件中使用*INITIAL CONDITIONS,TYPE=SOLUTION定义状态变量初值
C                           或采用子程序SDVINI来定义
C-----------------------------------------------------------------------------------
C    AK即k,AN即n,RF破坏比,c即凝聚力,FA即内摩擦角,PA大气压力,VKB即Kb,VNB即nb,AUR为卸荷Kur/k
      DIMENSION PS(NDI)
        AK=PROPS(1)
        AN=PROPS(2)
        RF=PROPS(3)
        C=PROPS(4)
        FA=PROPS(5)/180*3.1415926
        PA=PROPS(6)
        VKB=PROPS(7)
        VNB=PROPS(8)
        AUR=PROPS(9)
C     S1S3O为历时上最大主应力差，S3O为前期固结应力，SSO为历史上最大应力水平
        S1S3O=STATEV(1)
        S3O=STATEV(2)
        SSO=STATEV(3)
      write(110,*)'S1S30=',S1S30
        CALL SPRINC(STRESS,PS,1,NDI,NSHR)
        write(*,*)"PS(1)=",PS(1)

zhangmeng

邓肯张模型很多人都开发过

whua8

我这个邓肯张是再论坛上下载的，由问题啊，运行的话只计算几步就没了，楼上的兄弟帮忙看看我的问题如何？ 谢谢了

哪位大侠能不能上传一个能正常运行的邓肯张模型inp和for文件。跪谢啊

zhangmeng

呵呵，小Case。邓肯张很多人都做过，你多看几个版本就会理解更深了

whua8

楼上的兄弟：“你多看几个版本就会理解更深了”   请问哪里下载呢？

下面是运行邓肯张模型的时候出现负特征值，请教怎么解决呢？

***WARNING: THE SYSTEM MATRIX HAS 4 NEGATIVE EIGENVALUES.

***WARNING: SOLVER PROBLEM.  ZERO PIVOT WHEN PROCESSING NODE 3 INSTANCE
             PART-1-1 D.O.F. 1
  Typically a zero pivot is associated with an overconstraint in the model.
However, ABAQUS was not able to indentify any chains of constraints that might
lead to an overconstraint at this node. This suggests that the model might be
in fact insufficiently constrained in this particular degree of freedom.
Please check the constraints used in this model.


***WARNING: Solver problem. Zero pivot when processing D.O.F. 1 of 1 nodes.
             The nodes have been identified in node set
             WarnNodeSolvProbZeroPiv_1_1_3_1_1.

zhangmeng

你到维谱上搜就有

whua8

邓肯张子程序发表的文章是很多的，但我要的是inp和umat啊  ，新手想好好学习学习，兄弟能否提供一个谢谢了，ytdoc@163.com

zhangmeng

那些文章后面就附有程序呀，你太不仔细看了！

VibInfo

这个应该是fortran77的格式吧

其中的数字不是变量，应该是续行

文件应该是.for的

linqus

从上面的计算来看，
似乎不是.for文件的连接问题，
还请  zhangmeng  兄指点。
谢谢 :)

ayao

楼主您现在邓肯张模型做的怎么样了？小弟也想用这个，希望您不吝赐教，能否给一个inp和umat啊，谢谢，我的邮箱是ayao521@126.com

zhangmeng

自定义的材料模型很难考虑稳定性的问题
所以经常出现矩阵奇异
我去开ABAQUS的发布会听说6.7版本Explicit也可以定义自定义材料了
但是以单元二次开发的形式出现

vbgdvu

回复 14 # zhangmeng 的帖子

本人使用邓肯张就出现数值奇异，不收敛的问题，还列出了很多单元极度变形。