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