Télécharger la plaquette Cast3M


fr en

Télécharger calcyf01.eso

Retour à la liste

Numérotation des lignes : 

calcyf01

  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