Sources Esope | Cast3M

calcymi
C CALCYMI   SOURCE    BR232186  12/10/04    21:15:17     7520
      SUBROUTINE CALCYMI (EPSIM,XGAMMT,XGAMMC,
     &                    XLAMM,XMUM,XDI,XEPS33,YMI)
C
C====&===1=========2=========3=========4=========5=========6=========7==
C Commentaires :   Subroutine permettant de calculer la partie du seuil
C                  lie à la partie membrane YM1
C                  Appelé par la routine lcgldm.eso
C
C Auteur       : B. Richard (Dr - Ing.) - CEA/DEN/DANS/DM2S/SEMT/EMSI
C====&===1=========2=========3=========4=========5=========6=========7==
C
C----DECLARATION GENERALES----------------------------------------------
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
      REAL*8 EPSIM(3,3),VAPM(3),VEPM(3,3)
C
C-----CALCUL DE LA PARTIE DU SEUIL YM1
C
      XTRA = EPSIM(1,1) + EPSIM(2,2)
 
      CALL JACOB4(EPSIM,2,VAPM,VEPM)
 
      CALL CALCGM(XTRA,XGAMMT,XGAMMC,XGM)
      CALL CALCGM(VAPM(1),XGAMMT,XGAMMC,XGM1)
      CALL CALCGM(VAPM(2),XGAMMT,XGAMMC,XGM2)
 
      YMI = (1.0D0/(1.0D0+XDI))**2.0D0
 
      YMI = YMI*((XLAMM/4.0D0*(XTRA+XEPS33)**2.0D0)*XGM+
     &      (XMUM/2.0D0)*(VAPM(1)**2.0D0)*XGM1+
     &      (XMUM/2.0D0)*(VAPM(2)**2.0D0)*XGM2)
 
      RETURN
      END