Télécharger la plaquette Cast3M


fr en

Télécharger probdef.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  probdef.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *==================================================================*
  5. *                                                                  *
  6. *       CAS TEST PROBDEF                                           *
  7. *       ----------------
  8. *                                                                  *
  9. *       CALCUL IDEALISE D UNE PROBABILITE DE DEFAILLANCE           *
  10. *                                                                  *
  11. *                                                                  *
  12. *                                                                  *
  13. * Probabilite qu une resistance    (R) soit inferieure a           *
  14. *                une sollicitation (S)                             *
  15. *                                                                  *
  16. *                                                                  *
  17. *  R et S sont des variables aléatoires log-normales               *
  18. *                                                                  *
  19. *                                                                  *
  20. *  La probabilité de défaillance et l'indice de fiabilité          *
  21. *  calculés sont comparés au résultat analytique                   *
  22. *                                                                  *
  23. *                                                                  *
  24. *                                                                  *
  25. *==================================================================*
  26. opti echo 0;
  27. *
  28. *
  29. *
  30. * Moyenne, ecart-type et coeff. de variation de R
  31. muR = 15.;
  32. cvR = 0.1;
  33. sigR = muR*cvR;
  34. *
  35. * Moyenne, ecart-type et coeff. de variation de S
  36. muS = 5.;
  37. cvS = 0.1;
  38. sigS = muS*cvS;
  39. *
  40. * Corrélation 
  41. * ATTENTION : le cas RHO différent de 0 n est pas traite !
  42. RHO = 0.;
  43. *
  44. *beta et proba theorique pour 2 lois log-normales
  45. *cvR = sigR / muR;
  46. *cvS = sigS / muS;
  47. lR = log(cvR**2+1);
  48. lS = log(cvS**2+1);
  49. lRS = log(RHO*cvR*cvS+1);
  50. *
  51. BETA_th = (log(muR/muS)*(((cvS**2+1)/(cvR**2+1))**0.5))/
  52. ((lR+lS-(2.*lRS))**0.5);
  53. *
  54. Pf_th = 1-(PROB 0. 1. 0. 3. BETA_th);
  55. *
  56. *
  57. * Calcul des 4 premiers moments statistiques de R et S
  58. *
  59. NpR = 10;
  60. NpS = 10;
  61. *
  62. TabR = QUADRATU 'LOGN' muR sigR NpR;
  63. TabS = QUADRATU 'LOGN' muS sigS NpS;
  64. *
  65. mr sr rr br = PARASTAT TabR;
  66. ms ss rs bs = PARASTAT TabS;
  67. *
  68. *
  69. MESS '%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%';
  70. MESS '%%% CALCUL DES PARAMETRES STATISTIQUES ';
  71. MESS '%%% MOYENNE       ' mr ms;
  72. MESS '%%% ECART_TYPE    ' sr ss;
  73. MESS '%%% SYMETRIE      ' rr rs;
  74. MESS '%%% APLATISSEMENT ' br bs;
  75. MESS '%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%'; 
  76. *
  77. *
  78. * Calcul de la proba de défaillance et de l indice de fiabilité
  79. Pf_RS = PROBABRS mr sr rr br ms ss rs bs  ;
  80. Beta_RS = INDIBETA -20. 20. Pf_RS;
  81. *
  82. * Calcul des erreurs
  83. *
  84. ErBeta = 100 *  (ABS ((Beta_RS - Beta_th) / Beta_th));
  85. *
  86. MESS '%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%';
  87. MESS '                       RESULTATS';
  88. MESS '%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%';
  89. MESS '%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%';
  90. MESS ' ';
  91. MESS 'Probabilite de defaillance par PROBABRS      ' Pf_RS;
  92. MESS 'Probabilite de defaillance theorique         ' Pf_th;
  93. MESS ' ';
  94. MESS 'Indice de fiabilite par PROBABRS       ' Beta_RS;
  95. MESS 'Indice de fiabilite theorique          ' Beta_th;
  96. MESS ' ';
  97. MESS '%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%';
  98. MESS ' ';
  99. MESS '  ERREUR par rapport a Beta Theorique  ' ErBeta '%';
  100. *==================================================================*            
  101. *                  code fonctionnement                             *            
  102. *==================================================================*            
  103. *                                                                               
  104. *     L'ecart maximum entre valeur theorique et calculee doit etre              
  105. * inferieure a 0.2 %.
  106. *
  107. SI (ErBeta <EG 0.2);
  108.     ERRE 0;                                                                     
  109. SINON;                                                                          
  110.     ERRE 5;                                                                     
  111. FINSI;                                                                          
  112. *
  113. *==================================================================*            
  114. *                 Temps de calcul et fin                           *            
  115. *==================================================================*            
  116. *                                                                               
  118. OPTI ECHO 1;
  119. *
  120. fin;
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  

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