Télécharger la plaquette Cast3M


fr en

Télécharger weibull.procedur

Retour à la liste

Numérotation des lignes : 

  1. * WEIBULL   PROCEDUR  CHAT      96/06/07    21:15:07     2167           
  2. * WEIBULL                      
  3. *--------------------------------------------------------------                 
  4. *                                                                               
  5. *  PROBABILITY OF FAILURE FOR BRITTLE MATERIALS
  6. *  WITH THE FOLLOWING EXPRESSION                
  7. *                                                                               
  8. *  SEGUENTE :  PR= 1 - EXP( -V0 * < ((SIG - SIGU)/SIG0)**N > )                  
  9. *                     
  10. *  Weibull Weakest Link Theory
  11. *  Principle of Indipendent Action (PIA) 
  12. *--------------------------------------------------------------               
  13. *    ENTREE  DANS L'ORDRE                                                       
  14. *                                                                               
  15. *      ZSIG     CHAMP PAR ELEMENT DE CONTRAINTES                      
  16. *      ZOBAFF   OBJET MODELE ASSOCIE A ZSIG                 
  17. *      ZV0      FATTORE CHE TIENE CONTO DEL  
  18. *               VOLUME D'INTEGRAZIONE EFFETTIVO   
  19. *      ZSIGU    TENSIONE LIMITE
  20. *      ZSIG0    TENSIONE DI NORMALIZZAZIONE 
  21. *      ZXN      Modulo di WEIBULL      
  22. *                                                                               
  23. *    SORTIE DANS L'ORDRE                                                        
  24. *                                                                               
  25. *      ZXP     PROBABILITE                                                      
  26. *      ZSIGF   CHAMP PAR ELEMENT  QU ON INTEGRE                                 
  27. *                                                                               
  28. *--------------------------------------------------------------                 
  29. DEBPROC WEIBULL  
  30.   ZSIG*MCHAML     ZOBMO*MMODEL                                      
  31.   ZV0*FLOTTANT    ZSIGU*FLOTTANT                                                
  32.   ZSIG0*FLOTTANT  ZXN*FLOTTANT ;                                                
  33. *--------------------------------------------------------------                 
  34. imode1 = ega (vale mode) 'AXIS' ;
  35. *
  36. * Principal stresses
  37. tpr = prin zsig zobmo ;
  38. t1  = exco si11 tpr   ;
  39. t2  = exco si22 tpr   ;
  40. t3  = exco si33 tpr   ;
  41. *
  42. * Tensile stresses only
  43. COEF11 = T1 MASQUE SUPERIEUR ZSIGU ;
  44. COEF22 = T2 MASQUE SUPERIEUR ZSIGU ;
  45. COEF33 = T3 MASQUE SUPERIEUR ZSIGU ;
  46. *
  47. SI1   = COEF11 * T1 * (1. / ZSIG0) ;
  48. SII1  = COEF22 * T2 * (1. / ZSIG0) ;
  49. SIII1 = COEF33 * T3 * (1. / ZSIG0) ;
  50. *
  51. SITR   = REDUIRE SI1   ZOBMO ;
  52. SIITR  = REDUIRE SII1  ZOBMO ;
  53. SIIITR = REDUIRE SIII1 ZOBMO ;
  54. *
  55. * Calculation of the integral
  56. SIG2 = ( SITR ** ZXN ) + (SIITR ** ZXN ) + ( SIIITR ** ZXN ) ;
  57. SI IMODE1 ;
  58.   zr   = coor 1 zsig ;
  59.   sig2 = sig2 * zr   ;
  60. FINSI;
  61. ZSIGF= SIG2 * ZV0 ;
  62.  
  63. XPP  = INTG ZOBMO ZSIGF si11;
  64. xpp= xpp + (INTG ZOBMO ZSIGF si22); 
  65. xpp = xpp + (INTG ZOBMO ZSIGF si33); 
  66. *                                                         
  67. * On volume integral value, probability of survival and probability
  68. * of failure.
  69. XPM = -1. * XPP ;                                                              
  70. XPQ = EXP  XPM  ; 
  71. ZXP = 1. - XPQ  ;                                                 
  72. *            
  73. FINPROC  ZXP ZSIGF ;
  74. * 
  75.  
  76.  

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