Esope source files

   1 : C RIGBEX    SOURCE    CHAT      05/01/13    02:56:35     5004           
   2 : C RIGBEX    SOURCE    CHAT      89/05/26    22:52:13
   3 :       SUBROUTINE RIGBEX(REL,LRE,ES,AL,TGL1)
   4 : C-----------------------------------------------------------------------        
   5 : C       ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE                 
   6 : C                      DE L ELEMENT DE BARRE                                    
   7 : C-----------------------------------------------------------------------        
   8 : C  ENTREE                                                                       
   9 : C          LRE     DIMENSION DU TABLEAU REL                                     
  10 : C          ES      YOUN * SECT                                                  
  11 : C                                                                               
  12 : C  SORTIE                                                                       
  13 : C          REL     MATRICE DE RIGIDITE ELEMENTAIRE                              
  14 : C          KERRE   0 SI PAS DE PB
  15 : C                  1 SI LONGUEUR NULLE                                          
  16 : C-----------------------------------------------------------------------        
  17 :       IMPLICIT INTEGER(I-N)
  18 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  19 :       PARAMETER(XZER=0.D0)
  20 : C
  21 :       DIMENSION REL(LRE,LRE),TGL1(2,LRE),C(2),RE(2,2),C1(2)
  22 : C                                                                               
  23 :       FAC2=ES/AL
  24 :       RE(1,1)=FAC2
  25 :       RE(1,2)=-FAC2
  26 :       RE(2,1)=-FAC2
  27 :       RE(2,2)=FAC2
  28 : C
  29 :       DO 5 I =1,LRE
  30 : C
  31 :       DO 6 I1=1,2
  32 :         C1(I1)=TGL1(I1,I)
  33 :     6 CONTINUE
  34 : 
  35 :       DO 1 J=1,2
  36 :          S=XZER
  37 :       DO 2 K=1,2
  38 :         S=S+C1(K)*RE(K,J)
  39 :     2 CONTINUE
  40 :         C(J)=S
  41 :     1 CONTINUE
  42 : C
  43 :       DO 3 J=1,LRE
  44 :         S=XZER
  45 :       DO 4 K=1,2
  46 :         S=S+C(K)*TGL1(K,J)
  47 :     4 CONTINUE
  48 :         REL(I,J)=S
  49 :     3 CONTINUE
  50 : C
  51 :     5 CONTINUE
  52 : *      WRITE(*,*) ' matrice  K  elem ='
  53 : *      WRITE(*,1190) ((REL(IA,IC),IA=1,LRE),IC=1,LRE)
  54 : *1190   FORMAT(12(12(1X,E8.2),/))
  55 :       RETURN                                                                    
  56 :       END                                                                       
  57 :  
  58 :