Télécharger crit1.eso

Retour à la liste

Numérotation des lignes :

  1. C CRIT1 SOURCE CHAT 05/01/12 22:30:38 5004
  2. FUNCTION CRIT1(XM,XP,THETA,XM0,XP0)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. C Include contenant quelques constantes dont XPI :
  6. -INC CCREEL
  7. C=======================================================================
  8. C CALCULE LA VALEUR DE LA FONCTION DE CHARGE
  9. C
  10. C ENTREES:XM VALEUR DU MOMENT
  11. C XP VALEUR DE L EFFORT NORMAL
  12. C THETA ANGLE DE LA FISSURE
  13. C
  14. C SORTIE:CRIT1 VALEUR DE LA FONCTION
  15. C=======================================================================
  16. CRI=2.D0*XP0*(1.D0-THETA/XPI)*XM
  17. CRI=CRI+XM0*SIN(THETA)*XP
  18. IF(CRI.GE.XZERO) THEN
  19. GO TO 1
  20. ELSE
  21. GO TO 2
  22. ENDIF
  23. 1 CONTINUE
  24. CRIT1=SIN(THETA)/2.D0
  25. CRIT1=CRIT1-COS((THETA+XPI*XP/XP0)/2.D0)
  26. CRIT1=CRIT1+XM/XM0
  27. RETURN
  28. 2 CONTINUE
  29. CRIT1=SIN(THETA)/2.D0
  30. CRIT1=CRIT1-COS((THETA-XPI*XP/XP0)/2.D0)
  31. CRIT1=CRIT1-XM/XM0
  32. RETURN
  33. END
  34.  
  35.  
  36.  
  37.  
  38.  
  39.  

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