Télécharger voneps.eso

Retour à la liste

Numérotation des lignes :

  1. C VONEPS SOURCE CHAT 05/01/13 04:09:30 5004
  2. FUNCTION VONEPS(EPSIL,ICLE,ALFAH,COMIS)
  3. C
  4. C CALCUL DE LA DEFORMATION DE VON MISES
  5. C
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8 (A-H,O-Z)
  8. DIMENSION EPSIL(*)
  9. DIMENSION COMIS(6)
  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
  25. C****
  26. C
  27. C ALFAH COEFFICIENT MULTIPLICATEUR DE LA FLEXION DANS LES COQUES
  28. C COMIS COEFFICIENT DES CONTRAINTES DANS LE VON MISES
  29. C
  30. CCCCCCCCCCCCCCCCCCCC
  31. COEFEP=SQRT(0.666666666666666666666667D0)
  32. UNALF=0.D0
  33. IF(ALFAH.GE.1.D-12) UNALF=1.D0/ALFAH
  34. C****
  35. GO TO (1,2,3,4,5,6,7,8,9,10,11,12),ICLE
  36. 1 CONTINUE
  37. . SQRT(EPSIL(1)*EPSIL(1)+EPSIL(2)*EPSIL(2)+EPSIL(3)*EPSIL(3) +
  38. . 0.5D0*(EPSIL(4)*EPSIL(4)+EPSIL(5)*EPSIL(5)+EPSIL(6)*EPSIL(6)))
  39. GO TO 999
  40. 2 VONEPS= SQRT((EPSIL(4)*EPSIL(4)+EPSIL(5)*EPSIL(5) +
  41. . 0.5D0*EPSIL(6)*EPSIL(6))*UNALF+EPSIL(1)*EPSIL(1)+
  42. . EPSIL(2)*EPSIL(2) + 0.5D0*EPSIL(3)*EPSIL(3))
  43. GO TO 999
  44. 3 VONEPS=SQRT(EPSIL(1)*EPSIL(1)+EPSIL(2)*EPSIL(2)
  45. . +0.5D0*EPSIL(3)*EPSIL(3))
  46. GO TO 999
  47. 4 VONEPS=ABS(EPSIL(3))
  48. GO TO 999
  49. 5 VONEPS=SQRT(0.5*((COMIS(2)+COMIS(3))*EPSIL(1)*EPSIL(1)
  50. .+(COMIS(1)+COMIS(3))*EPSIL(2)*EPSIL(2)+(COMIS(1)+COMIS(2))*
  51. . EPSIL(3)*EPSIL(3))
  52. . -COMIS(1)*EPSIL(2)*EPSIL(3)-COMIS(2)*EPSIL(1)*EPSIL(3)-COMIS(3)*
  53. . EPSIL(1)*EPSIL(2)+0.5*COMIS(4)*EPSIL(4)*EPSIL(4)+COMIS(5)*0.5*
  54. . EPSIL(5)*EPSIL(5)+0.5*COMIS(6)*EPSIL(6)*EPSIL(6) )
  55. GO TO 999
  56. 6 VONEPS=SQRT(EPSIL(1)*EPSIL(1)+EPSIL(2)*EPSIL(2)+EPSIL(3)*EPSIL(3)
  57. . +0.5D0*EPSIL(4)*EPSIL(4))
  58. GO TO 999
  59. 7 VONEPS=SQRT(EPSIL(1)*EPSIL(1)+UNALF*EPSIL(4)*EPSIL(4))
  60. GO TO 999
  61. 8 VONEPS=ABS(EPSIL(1))
  62. GO TO 999
  63. C
  64. C*** ELEMENT COQUE EPAISSE
  65. C
  66. 9 VONEPS=SQRT((EPSIL(2)*EPSIL(2)+EPSIL(3)*EPSIL(3)
  67. . +0.5D0*(EPSIL(4)*EPSIL(4)+EPSIL(5)*EPSIL(5)+
  68. . EPSIL(6)*EPSIL(6)))*UNALF+EPSIL(8)*EPSIL(8)+
  69. . EPSIL(9)*EPSIL(9) +0.5D0*(EPSIL(10)*
  70. . EPSIL(10)+EPSIL(11)*EPSIL(11)+EPSIL(12)*EPSIL(12)))
  71. GO TO 999
  72. 10 VONEPS=SQRT(EPSIL(1)*EPSIL(1)+EPSIL(3)*EPSIL(3))
  73. . +COMIS(1)*EPSIL(2) + COMIS(2)
  74. GO TO 999
  75. C
  76. C POUTRES FORMULATION GLOBALE
  77. C ATTENTION VERIFIER LE COEFEP
  78. C
  79. 11 VONEPS=SQRT(EPSIL(1)*EPSIL(1)+
  80. . EPSIL(5)*EPSIL(5)+
  81. . EPSIL(6)*EPSIL(6)+
  82. . EPSIL(4)+EPSIL(4))
  83. GO TO 999
  84. C
  85. C TUYAUX
  86. C
  87. 12 VONEPS=SQRT(EPSIL(1)*EPSIL(1)+
  88. . EPSIL(5)*EPSIL(5)+
  89. . EPSIL(6)*EPSIL(6)+
  90. . EPSIL(4)+EPSIL(4))
  91. 999 VONEPS=VONEPS*COEFEP
  92. RETURN
  93. END
  94.  
  95.  

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