Télécharger cerrig.eso

Retour à la liste

Numérotation des lignes :

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

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