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.  
  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. DO 1 I=1,LRE
  34. DO 2 J=1,LRE
  35. REL(I,J)=0.D0
  36. 2 CONTINUE
  37. 1 CONTINUE
  38. *
  39. IF(IFOUR.EQ.0.OR.(IFOUR.EQ.1.AND.NIFOUR.EQ.0)) THEN
  40. FAC=2*XPI*ES/R
  41. ELSE
  42. FAC=XPI*ES/R
  43. ENDIF
  44. IF(IFOUR.EQ.0) THEN
  45. REL(1,1)=FAC
  46. ELSE IF(IFOUR.EQ.1) THEN
  47. AN= DBLE(NIFOUR)
  48. REL(1,1)=FAC
  49. REL(1,3)=FAC*AN
  50. REL(3,1)=REL(1,3)
  51. REL(3,3)=FAC*AN*AN
  52. ELSE
  53. KERRE=2
  54. RETURN
  55. ENDIF
  56. RETURN
  57. END
  58.  
  59.  
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  

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