Esope source files

   1 : C CERRIG    SOURCE    CHAT      05/01/12    21:52:29     5004           
   2 :       SUBROUTINE CERRIG(REL,LRE,ES,COORD,KERRE)
   3 : C-----------------------------------------------------------------------
   4 : C       ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE
   5 : C                      DE L ELEMENT CERCE
   6 : C-----------------------------------------------------------------------
   7 : C  ENTREE
   8 : C          LRE     DIMENSION DU TABLEAU REL
   9 : C          ES      YOUN * SECT
  10 : C          COORD   COORDONNEES DU NOEUD
  11 : C     
  12 : C  SORTIE
  13 : C          REL     MATRICE DE RIGIDITE ELEMENTAIRE
  14 : C          KERRE   0 SI PAS DE PB
  15 : C                  1 SI RAYON NUL
  16 : C                  2 CAS IMPOSSIBLE
  17 : C-----------------------------------------------------------------------
  18 :       IMPLICIT INTEGER(I-N)
  19 :       IMPLICIT REAL*8(A-H,O-Z)
  20 : -INC CCOPTIO
  21 : -INC CCREEL
  22 :       DIMENSION REL(LRE,*),COORD(3,*)
  23 : *     
  24 :       KERRE=0
  25 :       R = COORD(1,1)
  26 :       IF(ABS(R).LE.XPETIT) THEN
  27 :          KERRE=1
  28 :          RETURN
  29 :       ENDIF
  30 : *     
  31 :       DO 1 I=1,LRE
  32 :          DO 2 J=1,LRE
  33 :             REL(I,J)=0.D0
  34 :  2       CONTINUE
  35 :  1    CONTINUE
  36 : *     
  37 :       IF(IFOUR.EQ.0.OR.(IFOUR.EQ.1.AND.NIFOUR.EQ.0)) THEN
  38 :          FAC=2*XPI*ES/R
  39 :       ELSE
  40 :          FAC=XPI*ES/R
  41 :       ENDIF
  42 :       IF(IFOUR.EQ.0) THEN
  43 :          REL(1,1)=FAC
  44 :       ELSE IF(IFOUR.EQ.1) THEN
  45 :          AN= DBLE(NIFOUR)
  46 :          REL(1,1)=FAC
  47 :          REL(1,3)=FAC*AN
  48 :          REL(3,1)=REL(1,3)
  49 :          REL(3,3)=FAC*AN*AN
  50 :       ELSE
  51 :          KERRE=2
  52 :          RETURN
  53 :       ENDIF
  54 :       RETURN
  55 :       END
  56 : 
  57 : 
  58 : 
  59 : 
  60 :  
  61 :  
  62 :