Time-dependent fiber bundles with local load sharing. II. General Weibull fibers

Phys Rev E Stat Nonlin Soft Matter Phys. 2009 Dec;80(6 Pt 2):066115. doi: 10.1103/PhysRevE.80.066115. Epub 2009 Dec 17.

Abstract

Fiber bundle models (FBMs) are useful tools in understanding failure processes in a variety of material systems. While the fibers and load sharing assumptions are easily described, FBM analysis is typically difficult. Monte Carlo methods are also hampered by the severe computational demands of large bundle sizes, which overwhelm just as behavior relevant to real materials starts to emerge. For large size scales, interest continues in idealized FBMs that assume either equal load sharing (ELS) or local load sharing (LLS) among fibers, rules that reflect features of real load redistribution in elastic lattices. The present work focuses on a one-dimensional bundle of N fibers under LLS where life consumption in a fiber follows a power law in its load, with exponent rho , and integrated over time. This life consumption function is further embodied in a functional form resulting in a Weibull distribution for lifetime under constant fiber stress and with Weibull exponent, beta. Thus the failure rate of a fiber depends on its past load history, except for beta=1 . We develop asymptotic results validated by Monte Carlo simulation using a computational algorithm developed in our previous work [Phys. Rev. E 63, 021507 (2001)] that greatly increases the size, N , of treatable bundles (e.g., 10(6) fibers in 10(3) realizations). In particular, our algorithm is O(N ln N) in contrast with former algorithms which were O(N2) making this investigation possible. Regimes are found for (beta,rho) pairs that yield contrasting behavior for large N. For rho>1 and large N, brittle weakest volume behavior emerges in terms of characteristic elements (groupings of fibers) derived from critical cluster formation, and the lifetime eventually goes to zero as N-->infinity , unlike ELS, which yields a finite limiting mean. For 1/2<or=rho<or=1 , however, LLS has remarkably similar behavior to ELS (appearing to be virtually identical for rho=1 ) with an asymptotic Gaussian lifetime distribution and a finite limiting mean for large N. The coefficient of variation follows a power law in increasing N but, except for rho=1, the value of the negative exponent is clearly less than 1/2 unlike in ELS bundles where the exponent remains 1/2 for 1/2<rho<or=1. For sufficiently small values 0<rho1, a transition occurs, depending on beta , whereby LLS bundle lifetimes become dominated by a few long-lived fibers. Thus the bundle lifetime appears to approximately follow an extreme-value distribution for the longest lived of a parallel group of independent elements, which applies exactly to rho=0. The lower the value of beta , the higher the transition value of rho , below which such extreme-value behavior occurs. No evidence was found for limiting Gaussian behavior for rho>1 but with 0<beta(rho+1)<1, as might be conjectured from quasistatic bundle models where beta(rho+1) mimics the Weibull exponent for fiber strength.

Publication types

  • Research Support, U.S. Gov't, Non-P.H.S.

MeSH terms

  • Algorithms
  • Cluster Analysis
  • Materials Testing
  • Models, Statistical
  • Monte Carlo Method
  • Normal Distribution
  • Physics / methods*