Télécharger probabrs.procedur

Retour à la liste

Numérotation des lignes :

  1. * PROBABRS PROCEDUR AM 99/07/13 21:15:06 3631
  2. 'DEBPROC' PROBABRS MuR*'FLOTTANT' SigR*'FLOTTANT'
  3. rR*'FLOTTANT' bR*'FLOTTANT'
  4. MuS*'FLOTTANT' SigS*'FLOTTANT'
  5. rS*'FLOTTANT' bS*'FLOTTANT';
  6. *|=====================================================================
  7. *|
  8. *| << OBJET >> :
  9. *|
  10. *| Procedure permettant le calcul de la probabilite de defaillance
  11. *| idealisee :
  12. *| Prob(R<S)
  13. *| Probabilite que la resistance (R) soit inferieure a
  14. *| la sollicitation (S)
  15. *|
  16. *| (la resistance et la sollicitation sont definies par leur
  17. *| quatre premiers moments statistiques)
  18. *|
  19. *|
  20. *| << ENTREE >> :
  21. *| -- ARGUMENTS OBLIGATOIRES DANS TOUS LES CAS --
  22. *|
  23. *| MuR, MuS : objets de type REEL designant la valeur moyenne
  24. *| de R et S.
  25. *|
  26. *| SigR, SigS : objets de type REEL designant l'ecart-type
  27. *| de R et S.
  28. *|
  29. *| rR, rS : objets de type REEL designant la valeur du coefficient
  30. *| de symétrie de R et S.
  31. *|
  32. *| bR, bS : objets de type REEL designant la valeur du coefficient
  33. *| d'aplatissement de R et S.
  34. *|
  35. *| << SORTIE >> :
  36. *|
  37. *| Pf : objet de type REEL designant la probabilite de défaillance
  38. *|
  39. *|=====================================================================
  40. *
  41. * ______________________
  42. *
  43. * DECLARATIONS GENERALES
  44. * ______________________
  45. *
  46. SI = 'MOT' 'SI' ; SINON = 'MOT' 'SINON' ; FINSI = 'MOT' 'FINSI' ;
  47. REPETER = 'MOT' 'REPETER'; QUITTER = 'MOT' 'QUITTER'; FIN = 'MOT' 'FIN';
  48. EXISTE = 'MOT' 'EXISTE' ; EGA = 'MOT' 'EGA' ;
  49. ERREUR = 'MOT' 'ERREUR' ;
  50. ENTIER = 'MOT' 'ENTIER ' ;
  51. FLOTTANT = 'MOT' 'FLOTTANT' ;
  52. MOT = 'MOT' 'MOT ' ;
  53. * _____
  54. *
  55. * TESTS
  56. * _____
  57. *
  58. * _________________
  59. *
  60. * DEBUT DES CALCULS
  61. * _________________
  62. *
  63. *
  64. Pf =0.;
  65. oldS = 0.;
  66. FR = 0.;
  67. NPAS = 600;
  68. tabmin = PROG (MuS-(6*SigS)) (MuR-(6*SigR));
  69. tabmax = PROG (MuS+(6*SigS)) (MuR+(6*SigR));
  70. binf = mini tabmin;
  71. bsup = maxi tabmax;
  72. h = (bsup-binf)/NPAS;
  73. REPE BOUC NPAS;
  74. absx = binf + (&BOUC*h);
  75. valR = PROB MuR SigR rR bR absx;
  76. valS = PROB MuS SigS rS bS absx;
  77. fS = (valS-oldS)/h;
  78. FR = valR;
  79. Pf = Pf + (h*fS*FR);
  80. oldS = valS;
  81. FIN BOUC;
  82. 'FINPROC' Pf;
  83.  

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