Télécharger calcymi.eso

Retour à la liste

Numérotation des lignes :

  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.  
  24. CALL JACOB4(EPSIM,2,VAPM,VEPM)
  25.  
  26. CALL CALCGM(XTRA,XGAMMT,XGAMMC,XGM)
  27. CALL CALCGM(VAPM(1),XGAMMT,XGAMMC,XGM1)
  28. CALL CALCGM(VAPM(2),XGAMMT,XGAMMC,XGM2)
  29.  
  30. YMI = (1.0D0/(1.0D0+XDI))**2.0D0
  31.  
  32. YMI = YMI*((XLAMM/4.0D0*(XTRA+XEPS33)**2.0D0)*XGM+
  33. & (XMUM/2.0D0)*(VAPM(1)**2.0D0)*XGM1+
  34. & (XMUM/2.0D0)*(VAPM(2)**2.0D0)*XGM2)
  35.  
  36. RETURN
  37. END
  38.  
  39.  
  40.  
  41.  
  42.  
  43.  
  44.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales