| Deformation within the Earth's lithosphere is largely controlled by the rheology of the rock. Ductile behavior in rocks is often associated with plasticity due to dislocation motion or diffusion under high pressures and temperatures. However, ductile behavior can also occur in brittle materials. An example would be cataclastic flow associated with folding at shallow crustal levels, steep subduction zones, and large-scale deformation at plate boundaries. Engineers utilize damage mechanics to model the continuum deformation of brittle materials. We utilize a modified form of damage mechanics where damage represents a reduction in frictional strength and includes a yield stress. We use this empirical approach to simulate the bending of the lithosphere.; We use numerical simulations to obtain elastostatic solutions for plate bending and where the stress exceeds a yield stress, we apply damage to reduce the elastic moduli. Damage is calculated at each time step by a power-law relationship of the ratio of the yield stress to stress and the yield strain to the strain. To test our method, we apply our damage rheology to a plate deforming under applied shear, a constant bending moment, and a constant load. We simulate a wide range of behaviors from slow relaxation to instantaneous failure, over timescales that span six orders of magnitude. Stress relaxation produces elastic-perfectly plastic behavior in cases where failure does not occur. For cases of failure, we observe a rapid increase in damage leading to failure. The changes in the rate of damage accumulation in failure cases are similar to the changes in b-values of acoustic emissions observed in triaxial compression tests of fractured rock and b-value changes prior to some large earthquakes. Thus continuum damage mechanics can simulate ductile behavior due to brittle mechanisms as well as observations of laboratory experiments and seismicity. |
