Esope source files

   1 : C RIGLI2    SOURCE    CHAT      05/01/13    03:00:02     5004           
   2 : C RIGLI2    SOURCE    DJER      94/12/19    21:15:07     654
   3 :       SUBROUTINE RIGLI2(REL,LRE,BPSS,WORK)
   4 : C====================================================================           
   5 : C
   6 : C     MODIFIE REL (LRE,LRE) MATRICE DE RIGIDITE EN INTRODUISANT DES
   7 : C     RAIDEUR LOCALES EQUIVALENTES : Klux Kluy Kluz Klrx Klry Klrz
   8 : C====================================================================           
   9 :       IMPLICIT INTEGER(I-N)
  10 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  11 :       DIMENSION REL(LRE,LRE),BPSS(3,3),WORK(9),C(6),RE(6,6),C1(3)
  12 : C
  12 :       CALL ZERO(REL,LRE,LRE)
  12 :       CALL ZERO(RE,6,6)
  13 :       XZER = 0.D0
  14 : C      
  15 :       DO I =1,6
  16 :         RE(I,I)=WORK(I)
  17 :       END DO
  18 : C
  19 :       DO 5 I =1,3
  20 :          DO 6 I1=1,3
  21 :             C1(I1)=BPSS(I1,I)
  22 :     6    CONTINUE
  23 :          DO 1 J=1,3
  24 :             S1=XZER
  25 :             S2=XZER
  26 :             DO 2 K=1,3
  27 :                S1=S1+C1(K)*RE(K,J)
  28 :                S2=S2+C1(K)*RE(K+3,J+3)
  29 :     2       CONTINUE
  30 :             C(J)  =S1
  31 :             C(J+3)=S2
  32 :     1    CONTINUE
  33 :          DO 3 J=1,3
  34 :             S1=XZER
  35 :             S2=XZER
  36 :             DO 4 K=1,3
  37 :                S1=S1+C(K)*BPSS(K,J)
  38 :                S2=S2+C(K+3)*BPSS(K,J)
  39 :     4       CONTINUE
  40 :             REL(I,J)    =S1
  41 :             REL(I+3,J+3)=S2
  42 :             REL(I+6,J+6)=S1
  43 :             REL(I+9,J+9)=S2
  44 :     3    CONTINUE
  45 :     5 CONTINUE
  46 : 
  47 : *-------------------------------------------------------------
  48 : *      WRITE(*,*) '**  REL APRES ( dans RIGLI2) ** '
  49 : *      DO 112 IA=1,LRE
  50 : *       WRITE(*,1991) (REL(IA,IC),IC=1,LRE)
  51 : *  112 CONTINUE
  52 : *1991  FORMAT(18(1X,E12.5))
  53 : *-------------------------------------------------------------
  54 :       RETURN                                                                    
  55 :       END                                                                       
  56 :  
  57 :