intstamp_dt.F File Reference

#include "implicit_f.inc"
#include "intstamp_c.inc"
#include "param_c.inc"
#include "sms_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	intstamp_dt (intstamp, ipari, neltst, ityptst, dt2t, nodnx_sms, diag_sms, ms, v, stifn, stifr)

Function/Subroutine Documentation

intstamp_dt()

subroutine intstamp_dt	(	type(intstamp_data), dimension(*)	intstamp,
		integer, dimension(npari,*)	ipari,
		integer	neltst,
		integer	ityptst,
			dt2t,
		integer, dimension(*)	nodnx_sms,
			diag_sms,
			ms,
			v,
			stifn,
			stifr )

Definition at line 30 of file intstamp_dt.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE intstamp_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      "intstamp_c.inc"
#include      "param_c.inc"
#include      "sms_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IPARI(NPARI,*), NELTST, ITYPTST, NODNX_SMS(*)
C     REAL
      my_real
     .   dt2t, diag_sms(*), ms(*), v(3,*), stifn(*), stifr(*)
      TYPE(INTSTAMP_DATA) INTSTAMP(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER KDIR, NN, K, IROT, MAIN, INTDAMP
C     REAL
      my_real
     .   mass, stf, inm, str, alpha, vis, 
     .   c, cvis(nintstamp), cvisr(nintstamp)
C-----------------------------------------------
      IF(idtmins/=2.AND.idtmins_int==0)THEN
        DO nn=1,nintstamp
          cvis(nn) =zero
          cvisr(nn)=zero
        ENDDO
        DO nn=1,nintstamp
          intdamp=intstamp(nn)%INTDAMP
          mass =intstamp(nn)%MASS
          alpha=intstamp(nn)%DAMP
          stf  =intstamp(nn)%STF
          c    =alpha*sqrt(four*mass*stf)
          cvis(nn)=cvis(nn)+c
          IF(intdamp/=0)THEN
            cvis(intdamp)=cvis(intdamp)+c
          END IF
          irot=intstamp(nn)%IROT
          IF(irot/=0)THEN
            inm=min(intstamp(nn)%IN(1),
     .              intstamp(nn)%IN(2),
     .              intstamp(nn)%IN(3))
            str  =intstamp(nn)%STR
            alpha=intstamp(nn)%DAMPR
            c    =alpha*sqrt(four*inm*str)
            cvisr(nn)=cvisr(nn)+c
            IF(intdamp/=0)THEN
              cvisr(intdamp)=cvisr(intdamp)+c
            END IF
          END IF
        ENDDO
        DO nn=1,nintstamp
          intdamp=intstamp(nn)%INTDAMP
          mass =intstamp(nn)%MASS
          stf  =intstamp(nn)%STF
          main =intstamp(nn)%MSR
          IF(stf/=zero)THEN
            alpha=cvis(nn)/sqrt(four*mass*stf)
            stf  =stf*(sqrt(one+alpha*alpha)+alpha)**2
            stifn(main)=stifn(main)+stf
          END IF
          irot=intstamp(nn)%IROT
          IF(irot/=0)THEN
            inm=min(intstamp(nn)%IN(1),
     .              intstamp(nn)%IN(2),
     .              intstamp(nn)%IN(3))
            str =intstamp(nn)%STR
            IF(str/=zero)THEN
              alpha=cvisr(nn)/sqrt(four*inm*str)
              str  =str*(sqrt(one+alpha*alpha)+alpha)**2
              stifr(main)=stifr(main)+str
            END IF
          END IF
        ENDDO
      ELSE
        DO nn=1,nintstamp
          cvis(nn) =zero
          cvisr(nn)=zero
        ENDDO
        DO nn=1,nintstamp
          intdamp=intstamp(nn)%INTDAMP
          mass =intstamp(nn)%MASS
          alpha=intstamp(nn)%DAMP
          stf  =intstamp(nn)%STF
          c    =alpha*sqrt(four*mass*stf)
          cvis(nn)=cvis(nn)+c
          IF(intdamp/=0)THEN
            cvis(intdamp)=cvis(intdamp)+c
          END IF
          irot=intstamp(nn)%IROT
          IF(irot/=0)THEN
            inm=min(intstamp(nn)%IN(1),
     .              intstamp(nn)%IN(2),
     .              intstamp(nn)%IN(3))
            str  =intstamp(nn)%STR
            alpha=intstamp(nn)%DAMPR
            c    =alpha*sqrt(four*inm*str)
            cvisr(nn)=cvisr(nn)+c
            IF(intdamp/=0)THEN
              cvisr(intdamp)=cvisr(intdamp)+c
            END IF
          END IF
        ENDDO
        DO nn=1,nintstamp
          main=intstamp(nn)%MSR
          IF(nodnx_sms(main)==0)THEN
            mass=intstamp(nn)%MASS
            stf =intstamp(nn)%STF
            IF(stf/=zero)THEN
              alpha=cvis(nn)/sqrt(four*mass*stf)
              stf  =stf*(sqrt(one+alpha*alpha)+alpha)**2
              stifn(main)=stifn(main)+stf
            END IF
            irot=intstamp(nn)%IROT
            IF(irot/=0)THEN
              inm=min(intstamp(nn)%IN(1),
     .                intstamp(nn)%IN(2),
     .                intstamp(nn)%IN(3))
              str =intstamp(nn)%STR
              IF(str/=zero)THEN
          alpha=cvisr(nn)/sqrt(four*inm*str)
          str  =str*(sqrt(one+alpha*alpha)+alpha)**2
                stifr(main)=stifr(main)+str
              END IF
            END IF
          ELSE
C
            mass =intstamp(nn)%MASS
            stf  =intstamp(nn)%STF
            vis  =cvis(nn)
C
C stabilite des amortisseurs (masse augment e)::
C DM=HALF*C*DTMINS/DTFACS <=> 2*DM/C=DTMINS/DTFACS 
            diag_sms(main)=diag_sms(main)+half*vis*dtmins/dtfacs
C
            irot=intstamp(nn)%IROT
            IF(irot/=0)THEN
              inm=min(intstamp(nn)%IN(1),
     .                intstamp(nn)%IN(2),
     .                intstamp(nn)%IN(3))
              str =intstamp(nn)%STR
              IF(str/=zero)THEN
          alpha=cvisr(nn)/sqrt(four*inm*str)
          str  =str*(sqrt(one+alpha*alpha)+alpha)**2
                stifr(main)=stifr(main)+str
              END IF
            END IF
          END IF
        ENDDO
      END IF
C
      RETURN