law104_upd.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	law104_upd (ifailg, nuparam, nuparf, uparam, uparf, nloc_dmg, imat, mlaw_tag, ipm, matparam)

Function/Subroutine Documentation

law104_upd()

subroutine law104_upd	(	integer	ifailg,
		integer	nuparam,
		integer	nuparf,
		intent(inout)	uparam,
		intent(in)	uparf,
		type (nlocal_str_)	nloc_dmg,
		integer	imat,
		type(mlaw_tag_), intent(inout)	mlaw_tag,
		integer, dimension(npropmi,*)	ipm,
		type(matparam_struct_), intent(inout)	matparam )

Definition at line 32 of file law104_upd.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE nlocal_reg_mod
      USE elbuftag_mod 
      USE matparam_def_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IFAILG,NUPARAM,NUPARF,IMAT,IPM(NPROPMI,*)
      my_real, DIMENSION(NUPARF) , INTENT(IN)  :: uparf
      my_real, DIMENSION(NUPARAM), INTENT(INOUT) :: uparam
      TYPE (NLOCAL_STR_) :: NLOC_DMG 
      TYPE(MLAW_TAG_),INTENT(INOUT)    :: MLAW_TAG
      TYPE(MATPARAM_STRUCT_) ,INTENT(INOUT)    :: MATPARAM
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER ILOC,IGURSON
      my_real 
     .  q1,q2,q3,epn,as,kw,f0,fc,fr,rlen,hkhi
C=======================================================================
       igurson = 0
      IF (ifailg == 1) THEN
        iloc  = matparam%NLOC
        !IGURSON = 0  no damage (default) 
        IF (iloc == 0) THEN
          igurson = 1 ! local damage
        ELSEIF (iloc == 1) THEN
          igurson = 2 ! non local forest(micromorphic)
        ELSEIF (iloc == 2) THEN
          igurson = 3 ! Non local (Peerling)
          mlaw_tag%NUVAR = mlaw_tag%NUVAR + 1
        ENDIF
        ipm(8,imat) = mlaw_tag%NUVAR
        ! Gurson yield criterion parameters
        q1    = uparf(2) 
        q2    = uparf(3) 
        q3    = uparf(4) 
        ! Trigger plastic strain for damage nucleation
        epn   = uparf(5) 
        ! Nucleation rate
        as    = uparf(6) 
        ! Nahshon-Hutchinson shear parameter
        kw    = uparf(7) 
        ! Void volume fraction at fracture
        fr    = uparf(8)
        ! Critical void volume fraction 
        fc    = uparf(9)
        ! Initial void volume fraction
        f0    = uparf(10)
        ! Non-local internal length
        rlen  = uparf(11)
        ! Micromorphic penalty parameter
        hkhi  = uparf(12)
c
        ! Storage of the non-local internal length
        IF (iloc>0) THEN
          nloc_dmg%LEN(imat) = max(nloc_dmg%LEN(imat), rlen)
          CALL get_lemax(nloc_dmg%LE_MAX(imat),nloc_dmg%LEN(imat))
          mlaw_tag%G_PLANL  = 1
          mlaw_tag%L_PLANL  = 1
          mlaw_tag%G_EPSDNL = 1
          mlaw_tag%L_EPSDNL = 1
        ENDIF
c
        ! Tag for damage output
        ! -> Number of output modes (stored in DMG(NEL,I), I>1)
        matparam%NMOD  = 5
        ! Total number of damage outputs
        ! -> DMG(NEL,1) = Global damage output 
        ! -> DMG(NEL,2:NMOD+1) = Damage modes output
        mlaw_tag%G_DMG = 1 + matparam%NMOD
        mlaw_tag%L_DMG = 1 + matparam%NMOD
        ! -> Modes allocation and definition
        ALLOCATE(matparam%MODE(matparam%NMOD))
        matparam%MODE(1) = "Void growth volume fraction fg"
        matparam%MODE(2) = "Nucleation volume fraction fn"
        matparam%MODE(3) = "Shear growth volume fraction fsh"
        matparam%MODE(4) = "Total void volume fraction ft"
        matparam%MODE(5) = "Effective void volume fraction f*"
c
        ! Storage of damage parameters
        uparam(30) = igurson
        uparam(31) = q1
        uparam(32) = q2
        uparam(33) = q3
        uparam(34) = epn
        uparam(35) = as
        uparam(36) = kw
        uparam(37) = fr
        uparam(38) = fc
        uparam(39) = f0
        uparam(40) = hkhi 
      ENDIF
c-----------
      RETURN