Sources Esope | Cast3M

hbmzw
C HBMZW     SOURCE    OF166741  26/05/11    21:15:17     12538          
 
*=======================================================================
* calcule la derivee Z,w
*      Z,w = [ Z0,w  0     0  ... ]
*            [      Z1,w   0  ... ]
*            [            ...     ]
*            [               Zh,w ]
*      Zj,w = [ -2j²wM     jC   ]
*             [ -jC      -2j²wM ]
*=======================================================================
 
      SUBROUTINE HBMZW(NT,NHBM,NDDL,KTQ,KTKAM,KdZw)
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
-INC TMDYNC
 
      INTEGER NT,NHBM,NDDL
 
      INTEGER I,J
c      REAL*8 AA,BB,Mi,Ci
      REAL*8 AA,BB
 
* Segment des variables
      MTQ=KTQ
* Segment des proprietes mecaniques
      MTKAM=KTKAM
* Segment local : matrice dZw
      MdZw = KdZw
 
*      Initialisation
      DO J = 1,NT
        DO I = 1,NT
          MdZw.MATRC(I,J)=0.D0
        ENDDO
      ENDDO
 
* Remplissage
      DO J=2,2*NHBM,2
        r_za = -2.D0 * (OMEG*(J/2)**2)
        r_zb = (J/2)
        j_z2 = NDDL*(1+(J-2))
        j_z1 = NDDL*(1+(J-1))
 
        DO I=1,NDDL
c          Mi = XM(I,1)
c          Ci = XASM(I,1)
c          AA = -2.*(OMEG*(J/2)**2)*Mi
c          BB = (J/2)*Ci
          AA = r_za * XM(I,1)
          BB = r_zb * XASM(I,1)
          i_z2 = j_z2 + I
          i_z1 = j_z1 + I
 
          MdZw.MATRC(i_z2,i_z2) = AA
          MdZw.MATRC(i_z2,i_z1) = BB
          MdZw.MATRC(i_z1,i_z2) = -BB
          MdZw.MATRC(i_z1,i_z1) = AA
        ENDDO
      ENDDO
 
      RETURN
      END