Esope source files

   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
  24 :       CALL JACOB4(EPSIF,2,VAPF,VEPF)
  25 : C
  26 : C     Calcul de XTRAK
  27 : C
  28 :       XTRAK = EPSIF(1,1)+EPSIF(2,2)  
  29 : C
  30 : C     Calcul de YF01
  31 : C
  32 :       IF (XTRAK.GE.0.0D0) THEN
  33 :          IF ((VAPF(1).GE.0.0D0).AND.(VAPF(2).GE.0.0D0)) THEN
  34 :             YF01 = 0.0D0
  35 :          ENDIF
  36 :          
  37 :          IF ((VAPF(1).LT.0.0D0).AND.(VAPF(2).LT.0.0D0)) THEN
  38 :             YF01 = XMUF*(VAPF(1)**2.0D0)+XMUF*(VAPF(2)**2.0D0)
  39 :          ENDIF
  40 :          
  41 :          IF ((VAPF(1).GE.0.0D0).AND.(VAPF(2).LT.0.0D0)) THEN
  42 :             YF01 = XMUF*(VAPF(2)**2.0D0)            
  43 :          ENDIF
  44 :          
  45 :          IF ((VAPF(1).LT.0.0D0).AND.(VAPF(2).GE.0.0D0)) THEN
  46 :             YF01 = XMUF*(VAPF(1)**2.0D0)
  47 :          ENDIF
  48 :          
  49 :          YF01 = (1.0D0-XGAMF)*YF01
  50 :          
  51 :       ELSE
  52 :       
  53 :          IF ((VAPF(1).GE.0.0D0).AND.(VAPF(2).GE.0.0D0)) THEN
  54 :             YF01 = 0.0D0
  55 :          ENDIF
  56 :          
  57 :          IF ((VAPF(1).LT.0.0D0).AND.(VAPF(2).LT.0.0D0)) THEN
  58 :             YF01 = XMUF*(VAPF(1)**2.0D0)+XMUF*(VAPF(2)**2.0D0)
  59 :          ENDIF
  60 :          
  61 :          IF ((VAPF(1).GE.0.0D0).AND.(VAPF(2).LT.0.0D0)) THEN
  62 :             YF01 = XMUF*(VAPF(2)**2.0D0)            
  63 :          ENDIF
  64 :          
  65 :          IF ((VAPF(1).LT.0.0D0).AND.(VAPF(2).GE.0.0D0)) THEN
  66 :             YF01 = XMUF*(VAPF(1)**2.0D0)
  67 :          ENDIF
  68 :          
  69 :          YF01 = YF01 + (XLAMF/2.0D0)*(XTRAK**2.0D0)
  70 :          
  71 :          YF01 = (1.0D0-XGAMF)*YF01
  72 :       ENDIF
  73 : 
  74 : 
  75 :       RETURN
  76 :       END
  77 :       
  78 :       
  79 :       
  80 :       
  81 :