Esope source files

   1 : C CALCYMI   SOURCE    BR232186  12/10/04    21:15:17     7520           
   2 :       SUBROUTINE CALCYMI (EPSIM,XGAMMT,XGAMMC,
   3 :      &                    XLAMM,XMUM,XDI,XEPS33,YMI)
   4 : C
   5 : C====&===1=========2=========3=========4=========5=========6=========7==
   6 : C Commentaires :   Subroutine permettant de calculer la partie du seuil
   7 : C                  lie Ã  la partie membrane YM1
   8 : C                  AppelÃ© par la routine lcgldm.eso
   9 : C
  10 : C Auteur       : B. Richard (Dr - Ing.) - CEA/DEN/DANS/DM2S/SEMT/EMSI
  11 : C====&===1=========2=========3=========4=========5=========6=========7==
  12 : C
  13 : C----DECLARATION GENERALES----------------------------------------------
  14 : C
  15 :       IMPLICIT INTEGER(I-N)
  16 :       IMPLICIT REAL*8(A-H,O-Z)
  17 :       
  18 :       REAL*8 EPSIM(3,3),VAPM(3),VEPM(3,3)
  19 : C
  20 : C-----CALCUL DE LA PARTIE DU SEUIL YM1
  21 : C
  22 :       XTRA = EPSIM(1,1) + EPSIM(2,2)
  23 : 
  23 :       CALL JACOB4(EPSIM,2,VAPM,VEPM)
  24 :       
  24 :       CALL CALCGM(XTRA,XGAMMT,XGAMMC,XGM)
  24 :       CALL CALCGM(VAPM(1),XGAMMT,XGAMMC,XGM1)
  24 :       CALL CALCGM(VAPM(2),XGAMMT,XGAMMC,XGM2)
  25 : 
  26 :       YMI = (1.0D0/(1.0D0+XDI))**2.0D0
  27 : 
  28 :       YMI = YMI*((XLAMM/4.0D0*(XTRA+XEPS33)**2.0D0)*XGM+
  29 :      &      (XMUM/2.0D0)*(VAPM(1)**2.0D0)*XGM1+
  30 :      &      (XMUM/2.0D0)*(VAPM(2)**2.0D0)*XGM2)
  31 : 
  32 :       RETURN
  33 :       END
  34 :       
  35 :       
  36 :       
  37 :       
  38 : 
  39 :