Télécharger calcyf01.eso

Retour à la liste

Numérotation des lignes :

  1. C CALCYF01 SOURCE BR232186 12/10/04 21:15:14 7520
  2. SUBROUTINE CALCYF01 (EPSIF,XGAMF,XLAMF,XMUF,YF01)
  3. C
  4. C====&===1=========2=========3=========4=========5=========6=========7==
  5. C Commentaires : Subroutine permettant de calculer la partie du seuil
  6. C lie à la partie flexion YF01
  7. C Appelé par la routine lcgldm.eso
  8. C
  9. C Auteur : B. Richard (Dr - Ing.) - CEA/DEN/DANS/DM2S/SEMT/EMSI
  10. C====&===1=========2=========3=========4=========5=========6=========7==
  11. C
  12. C----DECLARATION GENERALES----------------------------------------------
  13. C
  14. IMPLICIT INTEGER(I-N)
  15. IMPLICIT REAL*8(A-H,O-Z)
  16.  
  17. REAL*8 EPSIF(3,3),VAPF(3),VEPF(3,3)
  18. C
  19. C-----CALCUL DE LA PARTIE DU SEUIL YF01
  20. C
  21.  
  22. C
  23. C Calcul des courbures principales (partie flexion)
  24. C
  25. CALL JACOB4(EPSIF,2,VAPF,VEPF)
  26. C
  27. C Calcul de XTRAK
  28. C
  29. XTRAK = EPSIF(1,1)+EPSIF(2,2)
  30. C
  31. C Calcul de YF01
  32. C
  33. IF (XTRAK.GE.0.0D0) THEN
  34. IF ((VAPF(1).GE.0.0D0).AND.(VAPF(2).GE.0.0D0)) THEN
  35. YF01 = 0.0D0
  36. ENDIF
  37.  
  38. IF ((VAPF(1).LT.0.0D0).AND.(VAPF(2).LT.0.0D0)) THEN
  39. YF01 = XMUF*(VAPF(1)**2.0D0)+XMUF*(VAPF(2)**2.0D0)
  40. ENDIF
  41.  
  42. IF ((VAPF(1).GE.0.0D0).AND.(VAPF(2).LT.0.0D0)) THEN
  43. YF01 = XMUF*(VAPF(2)**2.0D0)
  44. ENDIF
  45.  
  46. IF ((VAPF(1).LT.0.0D0).AND.(VAPF(2).GE.0.0D0)) THEN
  47. YF01 = XMUF*(VAPF(1)**2.0D0)
  48. ENDIF
  49.  
  50. YF01 = (1.0D0-XGAMF)*YF01
  51.  
  52. ELSE
  53.  
  54. IF ((VAPF(1).GE.0.0D0).AND.(VAPF(2).GE.0.0D0)) THEN
  55. YF01 = 0.0D0
  56. ENDIF
  57.  
  58. IF ((VAPF(1).LT.0.0D0).AND.(VAPF(2).LT.0.0D0)) THEN
  59. YF01 = XMUF*(VAPF(1)**2.0D0)+XMUF*(VAPF(2)**2.0D0)
  60. ENDIF
  61.  
  62. IF ((VAPF(1).GE.0.0D0).AND.(VAPF(2).LT.0.0D0)) THEN
  63. YF01 = XMUF*(VAPF(2)**2.0D0)
  64. ENDIF
  65.  
  66. IF ((VAPF(1).LT.0.0D0).AND.(VAPF(2).GE.0.0D0)) THEN
  67. YF01 = XMUF*(VAPF(1)**2.0D0)
  68. ENDIF
  69.  
  70. YF01 = YF01 + (XLAMF/2.0D0)*(XTRAK**2.0D0)
  71.  
  72. YF01 = (1.0D0-XGAMF)*YF01
  73. ENDIF
  74.  
  75.  
  76. RETURN
  77. END
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  

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