Sources Esope | Cast3M

rigbex
C RIGBEX    SOURCE    CHAT      05/01/13    02:56:35     5004
C RIGBEX    SOURCE    CHAT      89/05/26    22:52:13
      SUBROUTINE RIGBEX(REL,LRE,ES,AL,TGL1)
C-----------------------------------------------------------------------
C       ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE
C                      DE L ELEMENT DE BARRE
C-----------------------------------------------------------------------
C  ENTREE
C          LRE     DIMENSION DU TABLEAU REL
C          ES      YOUN * SECT
C
C  SORTIE
C          REL     MATRICE DE RIGIDITE ELEMENTAIRE
C          KERRE   0 SI PAS DE PB
C                  1 SI LONGUEUR NULLE
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      PARAMETER(XZER=0.D0)
C
      DIMENSION REL(LRE,LRE),TGL1(2,LRE),C(2),RE(2,2),C1(2)
C
      FAC2=ES/AL
      RE(1,1)=FAC2
      RE(1,2)=-FAC2
      RE(2,1)=-FAC2
      RE(2,2)=FAC2
C
      DO 5 I =1,LRE
C
      DO 6 I1=1,2
        C1(I1)=TGL1(I1,I)
    6 CONTINUE
 
      DO 1 J=1,2
         S=XZER
      DO 2 K=1,2
        S=S+C1(K)*RE(K,J)
    2 CONTINUE
        C(J)=S
    1 CONTINUE
C
      DO 3 J=1,LRE
        S=XZER
      DO 4 K=1,2
        S=S+C(K)*TGL1(K,J)
    4 CONTINUE
        REL(I,J)=S
    3 CONTINUE
C
    5 CONTINUE
*      WRITE(*,*) ' matrice  K  elem ='
*      WRITE(*,1190) ((REL(IA,IC),IA=1,LRE),IC=1,LRE)
*1190   FORMAT(12(12(1X,E8.2),/))
      RETURN
      END