dtcdk_reg.F File Reference

#include "implicit_f.inc"
#include "scr02_c.inc"
#include "scr18_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	dtcdk_reg (nloc_dmg, thk, nel, off, le, imat, nddl, dt2t)

Function/Subroutine Documentation

dtcdk_reg()

subroutine dtcdk_reg	(	type(nlocal_str_), target	nloc_dmg,
		intent(in)	thk,
		integer	nel,
		intent(in)	off,
		intent(in)	le,
		integer	imat,
		integer	nddl,
		intent(inout)	dt2t )

Definition at line 30 of file dtcdk_reg.F.

C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE nlocal_reg_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "scr02_c.inc"
#include      "scr18_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER                                    :: NEL,IMAT,NDDL
      my_real, INTENT(INOUT)                     ::
     .  dt2t
      my_real, DIMENSION(NEL), INTENT(IN)        :: 
     .  le,off,thk
      TYPE(NLOCAL_STR_), TARGET                  :: NLOC_DMG 
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real 
     .  l2,zeta,dtnl_th,dtnl,dampfac,le_max
      ! Coefficient for non-local stability to take into account damping
      my_real, PARAMETER :: cdamp = 0.7d0
C=======================================================================
      ! Recovering non-local parameters
      l2     = nloc_dmg%LEN(imat)**2 ! Non-local internal length ** 2
      zeta   = nloc_dmg%DENS(imat)   ! Non-local density
      le_max = nloc_dmg%LE_MAX(imat) ! Maximal length of convergence
c      
      !-----------------------------------------------------------------------
      ! Computing non-local timestep
      !-----------------------------------------------------------------------
      IF (nodadt == 0) THEN
        DO i = 1,nel
          ! If the element is not broken, normal computation
          IF (off(i)/=zero) THEN
            ! Non-local critical time-step in the plane
            dtnl = (two*(min(le(i),le_max))*sqrt(three*zeta))/
     .                  sqrt(twelve*l2 + (min(le(i),le_max))**2)
            ! Non-local critical time-step in the thickness
            IF (nddl>1) THEN
              IF (nddl > 2) THEN
                dtnl_th = (two*(min(thk(i)/nddl,le_max))*sqrt(three*zeta))/
     .                       sqrt(twelve*l2 + (min(thk(i)/nddl,le_max))**2)
              ELSE
                dtnl_th = (two*(min(thk(i),le_max))*sqrt(three*zeta))/
     .                       sqrt(twelve*l2 + (min(thk(i),le_max))**2)
              ENDIF
            ELSE
              dtnl_th = ep20
            ENDIF
            ! Retaining the minimal value
            dt2t = min(dt2t,dtfac1(1)*cdamp*dtnl_th,dtfac1(1)*cdamp*dtnl)
          ENDIF
        ENDDO
      ENDIF
c