A multiscale fatigue analysis model is developed for brittle composite materials. The mathematical homogenization theory is generalized to account for multiscale damage effects in heterogeneous media and a closed form expression relating nonlocal microphase fields to the overall strains and damage is derived. The evolution of fatigue damage is approximated by the first order initial value problem with respect to the number of load cycles. An efficient integrator is developed for the numerical integration of the continuum damage based fatigue cumulative law. The accuracy and computational efficiency of the proposed model for both low-cycle and high-cycle fatigue are investigated by numerical experimentation.
IntroductionIn recent years several fatigue models have been developed within the framework of continuum damage mechanics (CDM) [10][12] [29][35] [36]. Within this framework, internal state variables are introduced to model the fatigue damage. The degradation of material response under cyclic loading is simulated using constitutive equations which couple damage cumulation and mechanical responses. The microcrack initiation and growth are lumped together in the form of the evolution of damage variables from zero to some critical value. Most of the existing CDM based fatigue damage models are based on the classical (local) continuum damage theory even though it is well known that the accumulation of damage leads to strain softening and loss of ellipticity in quasi-static problems and hyperbolicity in dynamic problems (see, for example, [45] revealed the intrinsic links between the nonlocal CDM theory and fracture mechanics providing a possibility for building a unified framework to simulate crack initiation, propagation and overall structural failure under cyclic loading.When applying the CDM based fatigue model to life prediction of engineering systems, the coupling between mechanical response and damage cumulation poses a major computational challenge. This is because the number of cycles to failure, especially for high-cycle fatigue, is usually as high as tens of millions or more, and therefore, it is practically not feasible to carry out a direct cycle-by-cycle simulation for the fully coupled models even with today's powerful computers. An efficient integrator, or the so-called