Télécharger cerksg.eso

Retour à la liste

Numérotation des lignes :

  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.  
  58.  

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