找回密码
 注册
Simdroid-非首页
查看: 378|回复: 6

[子程序] 请教,如何利用USDFLD实现渗透系数随应变改变

[复制链接]
发表于 2012-3-3 21:18:30 | 显示全部楼层 |阅读模式 来自 江苏南京
本帖最后由 peixingkai 于 2012-3-4 14:08 编辑

在《ABAQUS在隧道及地下工程中的应用》的6.6.1利用子程序USDFLD实现了渗透系数随应变的改变(子程序内容如下),可惜小弟菜鸟一只,有两个问题望高手指点:对于状态变量进行STATEV(2) =(0.1+EV)/(1+EV)及STATEV(3) =(8.26*(1+10*EV)**3)/(1+EV)/8.1赋值:
其原因是什么?这两个公式的意义?能否将其去掉?

SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,
     1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,
     2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO,
     3 LACCFLA)
C
      INCLUDE 'ABA_PARAM.INC'
C
      CHARACTER*80 CMNAME,ORNAME
      CHARACTER*3  FLGRAY(15)
      DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),
     1 T(3,3),TIME(2)
      DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*),
     1 COORD(*)
C
C Absolute value of current strain:
      CALL GETVRM('E',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,
     1 MATLAYO,LACCFLA)
      
      EPS1=ARRAY(1)
      EPS2=ARRAY(2)
      EPS3=ARRAY(3)
         
      EV=EPS1+EPS2+EPS3
C field variable
      FIELD(1)=EV

c        print*,eps1,eps2,ev

C state variable
      STATEV(1) =EV
      STATEV(2) =(0.1+EV)/(1+EV)
      STATEV(3) =(8.26*(1+10*EV)**3)/(1+EV)/8.1


C If error, write comment to .DAT file:
      IF(JRCD.NE.0)THEN
       WRITE(6,*) 'REQUEST ERROR IN USDFLD FOR ELEMENT NUMBER ',
     1     NOEL,'INTEGRATION POINT NUMBER ',NPT
      ENDIF
C
      RETURN
      END


 楼主| 发表于 2012-3-4 14:10:58 | 显示全部楼层 来自 江苏南京
Simdroid开发平台
求高手出手相助,已经被困几天了
回复 不支持

使用道具 举报

发表于 2012-12-21 14:33:37 | 显示全部楼层 来自 陕西西安
这是隧洞中那个例题吧。整体看看inp文件你就知道了
回复 不支持

使用道具 举报

发表于 2013-9-30 16:46:30 | 显示全部楼层 来自 广西南宁
peixingkai 发表于 2012-3-4 14:10
求高手出手相助,已经被困几天了

这个问题你解决了吗,如果解决了麻烦告诉我一下,谢谢了 qq357229580
回复 不支持

使用道具 举报

发表于 2014-1-18 13:31:59 | 显示全部楼层 来自 天津
楼主方便留下QQ吗,谢谢
回复 不支持

使用道具 举报

发表于 2019-3-18 16:15:27 | 显示全部楼层 来自 山东济南
:@:@:@:@:@:@
回复 不支持

使用道具 举报

发表于 2019-5-11 22:38:14 | 显示全部楼层 来自 黑龙江
qq353504387一起学习探讨一下
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 16:49 , Processed in 0.031768 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表