Télécharger la plaquette Cast3M


fr en

Télécharger cerksg.eso

Retour à la liste

Numérotation des lignes : 

cerksg

  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.  
  21. -INC PPARAM
  22. -INC CCOPTIO
  23. -INC CCREEL
  24.       DIMENSION REL(LRE,*),COORD(3,*)
  25. *
  26.       KERRE=0
  27.       R = COORD(1,1)
  28.       IF(ABS(R).LE.XPETIT) THEN
  29.          KERRE=1
  30.          RETURN
  31.       ENDIF
  32. *
  33.       IF(IFOUR.EQ.0.OR.(IFOUR.EQ.1.AND.NIFOUR.EQ.0)) THEN
  34.          FAC=2*XPI*EFFORT/R
  35.       ELSE
  36.          FAC=XPI*EFFORT/R
  37.       ENDIF
  38.       IF(IFOUR.EQ.0) THEN
  39.          REL(1,1)=FAC
  40.       ELSE IF(IFOUR.EQ.1) THEN
  41.          AN= DBLE(NIFOUR)
  42.          AN2=AN*AN
  43.          REL(1,1)=FAC *(1.D0+AN2)
  44.          REL(3,3)=REL(1,1)
  45.          REL(1,3)=FAC*AN*2.D0
  46.          REL(3,1)=REL(1,3)
  47.       ELSE
  48.          KERRE=2
  49.          RETURN
  50.       ENDIF
  51.       RETURN
  52.       END
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales