Sources Esope | Cast3M

rigli2
C RIGLI2    SOURCE    CHAT      05/01/13    03:00:02     5004
C RIGLI2    SOURCE    DJER      94/12/19    21:15:07     654
      SUBROUTINE RIGLI2(REL,LRE,BPSS,WORK)
C====================================================================
C
C     MODIFIE REL (LRE,LRE) MATRICE DE RIGIDITE EN INTRODUISANT DES
C     RAIDEUR LOCALES EQUIVALENTES : Klux Kluy Kluz Klrx Klry Klrz
C====================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION REL(LRE,LRE),BPSS(3,3),WORK(9),C(6),RE(6,6),C1(3)
C
      CALL ZERO(REL,LRE,LRE)
      CALL ZERO(RE,6,6)
      XZER = 0.D0
C
      DO I =1,6
        RE(I,I)=WORK(I)
      END DO
C
      DO 5 I =1,3
         DO 6 I1=1,3
            C1(I1)=BPSS(I1,I)
    6    CONTINUE
         DO 1 J=1,3
            S1=XZER
            S2=XZER
            DO 2 K=1,3
               S1=S1+C1(K)*RE(K,J)
               S2=S2+C1(K)*RE(K+3,J+3)
    2       CONTINUE
            C(J)  =S1
            C(J+3)=S2
    1    CONTINUE
         DO 3 J=1,3
            S1=XZER
            S2=XZER
            DO 4 K=1,3
               S1=S1+C(K)*BPSS(K,J)
               S2=S2+C(K+3)*BPSS(K,J)
    4       CONTINUE
            REL(I,J)    =S1
            REL(I+3,J+3)=S2
            REL(I+6,J+6)=S1
            REL(I+9,J+9)=S2
    3    CONTINUE
    5 CONTINUE
 
*-------------------------------------------------------------
*      WRITE(*,*) '**  REL APRES ( dans RIGLI2) ** '
*      DO 112 IA=1,LRE
*       WRITE(*,1991) (REL(IA,IC),IC=1,LRE)
*  112 CONTINUE
*1991  FORMAT(18(1X,E12.5))
*-------------------------------------------------------------
      RETURN
      END