Télécharger vonmis.eso

Retour à la liste

Numérotation des lignes :

  1. C VONMIS SOURCE CHAT 05/01/13 04:09:37 5004
  2. FUNCTION VONMIS(SIGMA,ICLE,ALFAH,COMIS)
  3. C
  4. C CALCUL DE LA CONTRAINTE DE VON MISES
  5. C
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8 (A-H,O-Z)
  8. DIMENSION SIGMA(*)
  9. DIMENSION COMIS(*)
  10. C****
  11. C
  12. C ICLE=1 CAS DES ELEMENTS MASSIFS
  13. C ICLE=2 CAS DES COQUES
  14. C ICLE=3 CAS DES MEMBRANES
  15. C ICLE=4 CAS DES CABLES
  16. C ICLE=5 CAS QUELCONQUE
  17. C ICLE=6 CAS DES CONTRAINTES PLANES
  18. C ICLE=7 CAS DES COQUES A NU=0. OU CONTRAINTES PLANES
  19. C ICLE=8 CAS DES MEMBRANES A NU=0. OU CONTRAINTES PLANES
  20. C ICLE=9 CAS DES COQUES MASSIFS
  21. C ICLE=10 CAS DES JOINTS
  22. C ICLE=11 CAS DES POUTRES
  23. C ICLE=12 CAS DES TUYAUX
  24. C ICLE=13 CAS DES COQUES AVEC CISAILLEMENT
  25. C ICLE=14 CAS DES ELEMENTS MASSIFS EN 1D (UNID)
  26. C
  27. C****
  28. C
  29. C ALFAH COEFFICIENT MULTIPLICATEUR DE LA FLEXION DANS LES COQUES
  30. C COMIS COEFFICIENT DES CONTRAINTES DANS LE VON MISES
  31. C
  32. C****
  33. GO TO (1,2,3,4,5,6,7,8,9,10,11,12,13,14),ICLE
  34. 1 CONTINUE
  35. 1SQRT(ABS(SIGMA(1)*SIGMA(1)+SIGMA(2)*SIGMA(2)+SIGMA(3)*SIGMA(3)
  36. 1-SIGMA(1)*SIGMA(2)-SIGMA(2)*SIGMA(3)-SIGMA(3)*SIGMA(1)+3.D0*
  37. 2(SIGMA(4)*SIGMA(4)+SIGMA(5)*SIGMA(5)+SIGMA(6)*SIGMA(6)))
  38. 3 )
  39. RETURN
  40. 2 VONMIS=
  41. . SQRT((SIGMA(4)*SIGMA(4)+SIGMA(5)*SIGMA(5)-SIGMA(4)*SIGMA(5)+
  42. 1 3.*SIGMA(6)*SIGMA(6))*ALFAH+SIGMA(1)*SIGMA(1)+SIGMA(2)*SIGMA(2)
  43. 2 -SIGMA(1)*SIGMA(2)+3.*SIGMA(3)*SIGMA(3))
  44. RETURN
  45. 3 VONMIS=SQRT(SIGMA(1)*SIGMA(1)+SIGMA(2)*SIGMA(2)-SIGMA(1)*SIGMA(2)
  46. 1 +3.D0*SIGMA(3)*SIGMA(3))
  47. RETURN
  48. 4 VONMIS=ABS(SIGMA(3))
  49. RETURN
  50. 5 VONMIS=SQRT(0.5*((COMIS(2)+COMIS(3))*SIGMA(1)*SIGMA(1)
  51. 1+(COMIS(1)+COMIS(3))*SIGMA(2)*SIGMA(2)+(COMIS(1)+COMIS(2))*
  52. 1 SIGMA(3)*SIGMA(3))
  53. 1 -COMIS(1)*SIGMA(2)*SIGMA(3)-COMIS(2)*SIGMA(1)*SIGMA(3)-COMIS(3)*
  54. 1 SIGMA(1)*SIGMA(2)+0.5*COMIS(4)*SIGMA(4)*SIGMA(4)+COMIS(5)*0.5*
  55. 1 SIGMA(5)*SIGMA(5)+0.5*COMIS(6)*SIGMA(6)*SIGMA(6) )
  56. RETURN
  57. 6 VONMIS=SQRT(SIGMA(1)*SIGMA(1)+SIGMA(2)*SIGMA(2)-SIGMA(1)*SIGMA(2)
  58. 1+3.*SIGMA(4)*SIGMA(4))
  59. GO TO 999
  60. 7 VONMIS=SQRT(SIGMA(1)*SIGMA(1)+ALFAH*SIGMA(4)*SIGMA(4))
  61. GO TO 999
  62. 8 VONMIS=ABS(SIGMA(1))
  63. GO TO 999
  64. C
  65. C*** ELEMENT COQUE EPAISSE
  66. C
  67. 9 VONMIS=SQRT((SIGMA(2)*SIGMA(2)+SIGMA(3)*SIGMA(3)-SIGMA(2)*
  68. 1 SIGMA(3)+3.D0*(SIGMA(4)*SIGMA(4)+SIGMA(5)*SIGMA(5)+
  69. 2 SIGMA(6)*SIGMA(6)))*ALFAH+SIGMA(8)*SIGMA(8)+
  70. 3 SIGMA(9)*SIGMA(9)-SIGMA(8)*SIGMA(9)+3.D0*(SIGMA(10)*
  71. 4 SIGMA(10)+SIGMA(11)*SIGMA(11)+SIGMA(12)*SIGMA(12)))
  72. GO TO 999
  73. 10 VONMIS=SQRT(SIGMA(1)*SIGMA(1)+SIGMA(3)*SIGMA(3))
  74. * +COMIS(1)*SIGMA(2) + COMIS(2)
  75. GO TO 999
  76. C
  77. C POUTRES FORMULATION GLOBALE
  78. C
  79. 11 VONMIS=SQRT(SIGMA(1)*SIGMA(1)+
  80. . SIGMA(5)*SIGMA(5)+
  81. . SIGMA(6)*SIGMA(6)+
  82. . SIGMA(4)*SIGMA(4))
  83. GO TO 999
  84. C
  85. C TUYAUX ( MEME CONVENTION QUE LES POUTRES ET
  86. C PAS DE PRESSION POUR LE MOMENT )
  87. C
  88. 12 VONMIS=SQRT(SIGMA(1)*SIGMA(1)+
  89. . SIGMA(5)*SIGMA(5)+
  90. . SIGMA(6)*SIGMA(6)+
  91. . SIGMA(4)*SIGMA(4))
  92. GO TO 999
  93. C
  94. C COQUES AVEC CISAILLEMENT TRANSVERSE
  95. C
  96. 13 VONMIS=SQRT(SIGMA(1)*SIGMA(1)+SIGMA(2)*SIGMA(2)-SIGMA(1)*SIGMA(2)
  97. . +3.*SIGMA(4)*SIGMA(4)+3.*SIGMA(5)*SIGMA(5)+3.*SIGMA(6)*SIGMA(6))
  98. GO TO 999
  99. C
  100. C= Modes de calcul UNIDIMENSIONNELs (1D)
  101. 14 VONMIS=SQRT(ABS( SIGMA(1)*SIGMA(1)+SIGMA(2)*SIGMA(2)
  102. . +SIGMA(3)*SIGMA(3)-SIGMA(1)*SIGMA(2)
  103. . -SIGMA(2)*SIGMA(3)-SIGMA(3)*SIGMA(1) ))
  104. C GOTO 999
  105. 999 RETURN
  106. END
  107.  
  108.  
  109.  
  110.  
  111.  

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