24小时热门版块排行榜

返回列表

lj521775

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 2197.3
散金: 466
红花: 5
帖子: 436
在线: 194.6小时
虫号: 781653
注册: 2009-05-28
性别: GG
专业: 零件成形制造

[求助] DEFORM软件中韧性断裂子程序

本人采用DEFORM软件预测挤压成形过程中的表面裂纹，损伤模型为修正的Cockroft-Latham断裂准则，欲将临界损伤值修正为关于温度和应变速率的函数，故需要采用Fortran语言编写用户子程序。由于缺乏这方面的基础，望大家帮忙编写Cockroft-Latham断裂准则的用户子程序，谢谢！

   SUBROUTINE USRDM1(STS,EFSTS,EFEPS,DAMAG,STRLMT,DTIME)

C********************************************************************
C
C    User damage model
C
C    This routine is used to demonstrate the implementation of
C    damage routine.  All the real variables should be double
C    precision. the definition of arguments are described as follows
C
C    INPUT :
C
C    NRT    = DAMAGE MODEL NUMBER
C    STS    = STRESS
C    EFSTS    = EFFECTIVE STRESS
C    EFEPS    = EFFECTIVE STRAIN RATE
C    DAMAG    = PREVIOUS ACCUMULATED DAMAGE
C    STRLMT = STRAIN LIMIT
C    DTIME    = TIME INCREMENT
C
C    OUTPUT
C
C    DAMAG    = NEW VALUE OF ACCUMULATED DAMAGE
C
C********************************************************************

   IMPLICIT INTEGER*4 (I,J,K,L,M,N), REAL*8 (A-H,O-Z)
C  Additional element data available from the common block ELMCOM
C    COMMON /ELMCOM/
C
C       RZE    : Four corner coordinates
C       URZE : Velocity
C       STSE : Stress
C       EPSE : Strain rate
C       EFEPSE  : effective strain rate
C       EFSTSE  : Effective stress
C       TEPSE : Total effective strain
C       RDTYE : Density
C       TEMPE : Temperature
C       DTMPE : Temperature rate
C       DAMAGE  : Damage value
C       USRE1 : Element user state variable 1
C       USRE2 : Element user state variable 2
C       USRNE : Nodal user state variables 1,2 at 4 nodes
C       NODEE : Connectivity
C       KELE : Global element number
C       KELEL : Local element number
C       KGROUP  : Material group number
C
   COMMON /ELMCOM/ RZE(2,4),URZE(2,4),STSE(4),EPSE(4),EFEPSE,EFSTSE,
   +             TEPSE,RDTYE,TEMPE(4),DTMPE(4),DAMAGE,
   +             USRE1(1500),USRE2(1500),
   +             USRNE(1500,4),NODEE(4),KELE,KELEL,KGROUP
   DIMENSION STS(4)
C
C
User codes
C
   RETURN
   END

[ Last edited by lj521775 on 2012-7-26 at 17:38 ]

回复此楼