Procedures | Cast3M

   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 :