Télécharger vnmisd.eso

Retour à la liste

Numérotation des lignes :

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

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