dtnoda_stifint.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	dtnoda_stifint (ms, stifn, dt_stifint)

Function/Subroutine Documentation

dtnoda_stifint()

subroutine dtnoda_stifint	(	dimension(numnod), intent(in)	ms,
		dimension(numnod), intent(in)	stifn,
		intent(inout)	dt_stifint )

Definition at line 30 of file dtnoda_stifint.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_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      "com04_c.inc"
C-----------------------------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real , INTENT(INOUT) ::  dt_stifint
      my_real , INTENT(IN) ::  ms(numnod) ,stifn(numnod)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER 
     .   N
      my_real 
     .   dtnod, dtfac, dtn
C----------------------------------------------------------
C    Estimation Time step needed for contact stiffness computation 
C----------------------------------------------------------
 
        dtfac = zep9
        dtnod = ep20
        DO n=1,numnod
           IF(stifn(n)>zero.AND.ms(n)>zero)THEN 
             dtn = dtfac*sqrt(two * ms(n) / stifn(n)) 
             dtnod = min(dtnod,dtn) 
           ENDIF                      
        ENDDO
C
        dt_stifint = dtnod
 
      RETURN