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

whua8 · 发表于 2006-8-3 21:05

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

以下是论坛上的邓肯张模型的一小部分，其中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 · 发表于 2006-8-4 11:55

你可以在帮助中找到

whua8 · 发表于 2006-8-4 22:24

下面这段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 · 发表于 2006-8-5 13:36

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

whua8 · 发表于 2006-8-6 11:13

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

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

zhangmeng · 发表于 2006-8-7 10:40

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

whua8 · 发表于 2006-8-7 21:42

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

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

***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 · 发表于 2006-8-8 19:06

你到维谱上搜就有

whua8 · 发表于 2006-8-9 10:27

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

zhangmeng · 发表于 2006-8-9 18:44

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

VibInfo · 发表于 2006-8-10 07:23

这个应该是fortran77的格式吧

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

文件应该是.for的

linqus · 发表于 2006-8-10 09:07

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

ayao · 发表于 2007-5-6 09:02

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

zhangmeng · 发表于 2007-5-14 12:08

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

vbgdvu · 发表于 2012-6-21 12:28

回复 14 # zhangmeng 的帖子

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

账号		自动登录	找回密码
密码			我要加入