dtnodamp.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "scr18_c.inc"
#include "com08_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	dtnodamp (itab, ms, in, stifn, stifr, dt2t, weight, igrnod, dampr, istop, i_damp_rdof_tab, icontact, ixc, x)

Function/Subroutine Documentation

dtnodamp()

subroutine dtnodamp	(	integer, dimension(*)	itab,
			ms,
			in,
			stifn,
			stifr,
			dt2t,
		integer, dimension(*)	weight,
		type (group_), dimension(ngrnod)	igrnod,
			dampr,
		integer	istop,
		integer, dimension(*)	i_damp_rdof_tab,
		integer, dimension(*)	icontact,
		integer, dimension(nixc,*)	ixc,
			x )

Definition at line 31 of file dtnodamp.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE groupdef_mod
      use element_mod , only : nixc
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "scr18_c.inc"
#include      "com08_c.inc"
#include      "param_c.inc"
C-----------------------------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ITAB(*),WEIGHT(*),ISTOP,
     .   I_DAMP_RDOF_TAB(*),ICONTACT(*),IXC(NIXC,*)
C     REAL
      my_real dt2t, stifn(*), stifr(*), 
     .   ms(*) ,in(*) ,dampr(nrdamp,*),x(3,*)
C-----------------------------------------------
      TYPE (GROUP_)  , DIMENSION(NGRNOD) :: IGRNOD
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER N,I,J,ISK,IGR,ND,K,L,LL
      my_real factb,dampbr,dampar,d_tstart,d_tstop,dampt,
     .   dt2n(numnod),bb,coeff,dt_cible,eta,beta,
     .   alpha,dampa0          
C----------------------------------------------------------
 
C--------> TAG of elements
        DO nd=1,ndamp
          IF (dampr(20,nd)>0) THEN
          DO i=1,numelc
            DO j=1,4
              k=ixc(j+1,i)
                IF (icontact(k)==1) THEN
                DO l=1,4
            ll = ixc(l+1,i)
            IF (ll/=k) THEN
              i_damp_rdof_tab(ll) = dampr(19,nd)
            ENDIF  
          END DO      
                ENDIF
            ENDDO
          ENDDO
            ENDIF
        ENDDO
 
        DO nd=1,ndamp
C-------->
          IF (ncycle==0) THEN
              alpha = dtmin1(11)/dtfac1(11)
              IF (alpha>0) THEN
                dampa0=max(dampr(9,nd),dampr(11,nd),dampr(13,nd))      
                dampr(9,nd)=min(dampr(9,nd),1/alpha)
                dampr(11,nd)=min(dampr(11,nd),1/alpha)
                dampr(13,nd)=min(dampr(13,nd),1/alpha)
                dampar=max(dampr(9,nd),dampr(11,nd),dampr(13,nd))
                IF (dampa0/=dampar) THEN          
                  print *,"WARNING ALPHA reduced to",1/alpha
                ENDIF            
              ENDIF
            ENDIF    
C-------->    
          igr   = nint(dampr(2,nd))
          isk   = nint(dampr(15,nd))
          factb = one
          d_tstart = dampr(17,nd)
          d_tstop  = dampr(18,nd)
          IF (tt>=d_tstart .AND. tt<=d_tstop) THEN
            dampar = max(dampr(9,nd),dampr(11,nd),dampr(13,nd))    
            dampbr = max(dampr(10,nd),dampr(12,nd),dampr(14,nd))      
C-------->      
            DO i=1,igrnod(igr)%NENTITY
              n=igrnod(igr)%ENTITY(i)
              IF(stifr(n)<=zero)THEN
                istop=-itab(n)
              ELSEIF(in(n)>zero)THEN
                dt2n(n)=dtfac1(11)*sqrt(two*in(n)/stifr(n))
              ENDIF                               
            ENDDO      
C-------->
            DO i=1,igrnod(igr)%NENTITY
              n=igrnod(igr)%ENTITY(i)
              dampt  = min(dt1,dt2n(n))*factb
                beta = min(dampbr,dampt)
                alpha = dampar                   
                bb = beta+half*half*alpha*dt2n(n)*dt2n(n)
                eta = sqrt(bb*bb + dt2n(n)*dt2n(n)) - bb        
                IF (idtmin(11)==3) THEN      
                  dt_cible = eta+0*max(eta,dtmin1(11))
              ELSE
                  dt_cible = eta
              ENDIF              
                coeff = (one+(2*beta/dt_cible))/((1-half*alpha*dt_cible))
              IF (dampr(19,nd)/=0) THEN                                               
                  IF ((i_damp_rdof_tab(n)+icontact(n))/=0) THEN
                    stifr(n) = stifr(n)*coeff                           
                  ENDIF
                ENDIF                            
            ENDDO
C-------->      
          ENDIF    
               
          END DO
 
C
      RETURN