Télécharger la plaquette Cast3M


fr en

Télécharger voneps.eso

Retour à la liste

Numérotation des lignes : 

voneps

  1. C VONEPS    SOURCE    PV        22/04/15    13:20:15     11344          
  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(*)
  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
  38.      . SQRT(EPSIL(1)*EPSIL(1)+EPSIL(2)*EPSIL(2)+EPSIL(3)*EPSIL(3) +
  39.      . 0.5D0*(EPSIL(4)*EPSIL(4)+EPSIL(5)*EPSIL(5)+EPSIL(6)*EPSIL(6)))
  40.       GO TO 999
  41.     2 VONEPS= SQRT((EPSIL(4)*EPSIL(4)+EPSIL(5)*EPSIL(5) +
  42.      .  0.5D0*EPSIL(6)*EPSIL(6))*UNALF+EPSIL(1)*EPSIL(1)+
  43.      .  EPSIL(2)*EPSIL(2) + 0.5D0*EPSIL(3)*EPSIL(3))
  44.       GO TO 999
  45.     3 VONEPS=SQRT(EPSIL(1)*EPSIL(1)+EPSIL(2)*EPSIL(2)
  46.      . +0.5D0*EPSIL(3)*EPSIL(3))
  47.       GO TO 999
  48.     4 VONEPS=ABS(EPSIL(3))
  49.       GO TO 999
  50.     5 VONEPS=SQRT(0.5*((COMIS(2)+COMIS(3))*EPSIL(1)*EPSIL(1)
  51.      .+(COMIS(1)+COMIS(3))*EPSIL(2)*EPSIL(2)+(COMIS(1)+COMIS(2))*
  52.      . EPSIL(3)*EPSIL(3))
  53.      . -COMIS(1)*EPSIL(2)*EPSIL(3)-COMIS(2)*EPSIL(1)*EPSIL(3)-COMIS(3)*
  54.      . EPSIL(1)*EPSIL(2)+0.5*COMIS(4)*EPSIL(4)*EPSIL(4)+COMIS(5)*0.5*
  55.      . EPSIL(5)*EPSIL(5)+0.5*COMIS(6)*EPSIL(6)*EPSIL(6)   )
  56.       GO TO 999
  57.     6 VONEPS=SQRT(EPSIL(1)*EPSIL(1)+EPSIL(2)*EPSIL(2)+EPSIL(3)*EPSIL(3)
  58.      . +0.5D0*EPSIL(4)*EPSIL(4))
  59.       GO TO 999
  60.     7 VONEPS=SQRT(EPSIL(1)*EPSIL(1)+UNALF*EPSIL(4)*EPSIL(4))
  61.       GO TO 999
  62.     8 VONEPS=ABS(EPSIL(1))
  63.       GO TO 999
  64. C
  65. C***  ELEMENT COQUE EPAISSE
  66. C
  67.     9 VONEPS=SQRT((EPSIL(2)*EPSIL(2)+EPSIL(3)*EPSIL(3)
  68.      .             +0.5D0*(EPSIL(4)*EPSIL(4)+EPSIL(5)*EPSIL(5)+
  69.      .             EPSIL(6)*EPSIL(6)))*UNALF+EPSIL(8)*EPSIL(8)+
  70.      .             EPSIL(9)*EPSIL(9) +0.5D0*(EPSIL(10)*
  71.      .             EPSIL(10)+EPSIL(11)*EPSIL(11)+EPSIL(12)*EPSIL(12)))
  72.       GO TO 999
  73.    10 VONEPS=SQRT(EPSIL(1)*EPSIL(1)+EPSIL(3)*EPSIL(3))
  74.      .      +COMIS(1)*EPSIL(2) + COMIS(2)
  75.       GO TO 999
  76. C
  77. C  POUTRES  FORMULATION GLOBALE
  78. C  ATTENTION  VERIFIER LE COEFEP
  79. C
  80.    11 VONEPS=SQRT(EPSIL(1)*EPSIL(1)+
  81.      .             EPSIL(5)*EPSIL(5)+
  82.      .             EPSIL(6)*EPSIL(6)+
  83.      .             EPSIL(4)+EPSIL(4))
  84.       GO TO 999
  85. C
  86. C  TUYAUX
  87. C
  88.    12 VONEPS=SQRT(EPSIL(1)*EPSIL(1)+
  89.      .             EPSIL(5)*EPSIL(5)+
  90.      .             EPSIL(6)*EPSIL(6)+
  91.      .             EPSIL(4)+EPSIL(4))
  92.   999 VONEPS=VONEPS*COEFEP
  93.       RETURN
  94.       END
  95.  
  96.  
  97.  

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