Esope source files

   1 : C CERKSG    SOURCE    CHAT      05/01/12    21:52:23     5004           
   2 :       SUBROUTINE CERKSG(REL,LRE,EFFORT,COORD,KERRE)
   3 : C-----------------------------------------------------------------------
   4 : C       ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE GEOMETRIQUE
   5 : C                      DE L ELEMENT CERCE
   6 : C-----------------------------------------------------------------------
   7 : C  ENTREE
   8 : C          LRE     DIMENSION DU TABLEAU REL
   9 : C          EFFORT  EFFORT DANS LA CERCE ( SUPPOSE TOUKOURS AXI )
  10 : C          COORD   COORDONNEES DU NOEUD
  11 : C     
  12 : C  SORTIE
  13 : C          REL     MATRICE DE RIGIDITE GEOMETRIQUE 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 :       IF(IFOUR.EQ.0.OR.(IFOUR.EQ.1.AND.NIFOUR.EQ.0)) THEN
  32 :          FAC=2*XPI*EFFORT/R
  33 :       ELSE
  34 :          FAC=XPI*EFFORT/R
  35 :       ENDIF
  36 :       IF(IFOUR.EQ.0) THEN
  37 :          REL(1,1)=FAC
  38 :       ELSE IF(IFOUR.EQ.1) THEN
  39 :          AN= DBLE(NIFOUR)
  40 :          AN2=AN*AN
  41 :          REL(1,1)=FAC *(1.D0+AN2)
  42 :          REL(3,3)=REL(1,1)
  43 :          REL(1,3)=FAC*AN*2.D0
  44 :          REL(3,1)=REL(1,3)
  45 :       ELSE
  46 :          KERRE=2
  47 :          RETURN
  48 :       ENDIF
  49 :       RETURN
  50 :       END
  51 : 
  52 : 
  53 : 
  54 : 
  55 :  
  56 :  
  57 :