Heterogeneous Stress Field Research Articles

The effect of stress barriers on unconventional-singularity-driven frictional rupture

Whether or not energy dissipation is localized in the vicinity of the rupture tip, and whether any distal energy dissipation far from the crack tip has a significant influence on rupture dynamics are key questions in the description of frictional ruptures, in particular regarding the application of Linear Elastic Fracture Mechanics (LEFM) to earthquakes. These questions are investigated experimentally using a 40-cm-long experimental frictional interface. Three independent pistons apply a normal load with a fourth piston applying a shear load, enabling the application of a heterogeneous stress state and stress barriers. After loading the frictional interface to a near-critical state, subsequent unloading of one normal-load piston leads to dynamic ruptures which propagate into the heterogeneous stress fields. The ruptures in these experiments are found to be driven by unconventional singularities, characterized by an ever-increasing breakdown work with slip, and as a result do not conform to the assumptions of LEFM. As these experimental stress barriers inhibit slip, they therefore also reduce the breakdown work occurring outside of the cohesive zone. It is shown that this distal weakening, far from the crack tip, must be considered for the accurate prediction of rupture arrest length. These experiments are performed in the context of a proposed stimulation technique for Enhanced Geothermal Systems (EGSs). It has previously been suggested, through theoretical arguments, that stress barriers could be induced through the manipulation of pore pressure such that there is reduced seismic hazard during the shear stimulation of EGSs. This stimulation technique, known as preconditioning, is demonstrated here to reduce the mechanical energy flux to the crack tip, G, while also increasing the fracture energy, Gc. Preconditioning is shown to be capable of arresting seismic rupture and reducing co-seismic slip, slip velocity, and seismic moment at preconditioning stresses which are reasonably achievable in the field. Due to the fully-coupled nature of seismic rupture and fault slip, preconditioning also reduces distal weakening and its contribution to the propagation of induced seismic ruptures. In a similar vein, heterogeneous pore pressure fields associated with some seismic swarms can be used to explain changes in stress drop within the swarm without recourse to material or total-stress heterogeneity.

Source analysis of low frequency seismicity at Mt. Vesuvius by a hybrid moment tensor inversion

Seismicity at Mt. Vesuvius has been relatively weak in the last decades. While the occurrence of shallow volcano-tectonic (VT) events at Mt. Vesuvius is well known, the occurrence of deeper low frequency events (LF) was only recently recognized. Previous source studies only targeted VT events, which were found to have quite heterogeneous focal mechanisms. In this paper, we perform for the first time the source inversion of LF seismicity at Mt. Vesuvius, analysing 27 LF events recorded from 2012 to 2021 with the aim to investigate their source processes. Given the challenges of analysing weak LF earthquakes, we implement a specific moment tensor (MT) inversion approach that combines the fit of displacement seismograms in the time domain and amplitude spectra in the frequency domain. The inversion is simultaneously performed for the source depth and moment tensor components in the 2–7 and 2–5 Hz frequency band, assuming either a full or deviatoric MT representation. Source parameter uncertainties are estimated by using a Bayesian bootstrapping scheme. Our results confirm a larger depth of LF events compared to VTs and show a strong heterogeneity of the LF seismic sources, which present various rupture types, different orientations and heterogeneous, whilst poorly resolved, non-double-couple components. The MT variability is qualitatively confirmed by significant differences among the recorded waveforms. The heterogeneity of both VT and LF source processes is attributed to complex source processes in a highly fractured seismogenic volume submitted to a heterogeneous stress field.

Laboratory Earthquakes Simulations—Emergence, Structure, and Evolution of Fault Heterogeneity

AbstractSeismic faults are known to exhibit a high level of spatial and temporal complexity, and the causes and consequences of this complexity have been the topic of numerous research works in the past decade. In this paper, we investigate the origins and the structure of this complexity by considering a numerical model of laboratory earthquake experiment, where we introduce a fault with homogeneous mechanical properties but allow it to evolve spontaneously to its natural level of complexity. This is achieved by coupling the elastic deformability of the off‐fault medium (and therefore allowing for heterogeneous stress fields to develop) and the discrete degradation and gouge formation at the fault plane (and therefore allowing for structural heterogeneity to develop). Numerical results show the development of persistent stress, damage, and gouge thickness heterogeneities, with a much larger variability in space than in time. Strong positive correlations are found between these quantities, which suggest a positive feedback between local normal stress and damage rate, only mildly mitigated by the mobility of the granular gouge in the interface. For a wide range of confining stresses, after a sufficient number of seismic cycles, the fault reaches a state of established disorder with a constant roughness, a certain amount of periodicity at the millimetric scale, and a power law decay of the Power Spectral Density at smaller spatial scales. The typical height‐to‐wavelength ratio of geometrical asperities and the correlations between stress and damage profiles are in good agreement with previous field or lab estimates.

The evolution of stresses and shear resistance on rough faults at large slip

The deviation of natural faults from planarity results in geometric asperities and a heterogeneous stress field, affecting the shear resistance and slip behavior. This study numerically examines the effects of macro-scale roughness, local friction, and wear on the evolution of the fault stresses and shear resistance with fault slip. In contrast to previous analytical and numerical studies, our quasistatic simulations allow for slip at the scale of roughness wavelengths, large roughness levels, and fault opening. The simulations indicate that analytical solutions, which predict a linear increase in shear resistance with slip, significantly overestimate the additional shear resistance from roughness, i.e., the roughness drag, at slips larger than several percent of the minimum roughness wavelength. Starting with mated fault surfaces, the average shear resistance on the fault initially rises rapidly, but the rate of increase diminishes significantly as slip increases and the faults open. While initially showing a larger increase in shear resistance than self-similar faults, at large slip, self-affine faults exhibit significantly smaller shear resistance, and the background stresses are bounded without the need for off-fault inelastic deformation. Although smaller than predicted by analytical solutions, the simulation results show that the additional resistance from macro-scale geometrical irregularities is important, enabling rough faults to operate under shear resistance significantly larger than that from local friction only. The impact of macro-scale roughness increases with the roughness amplitude and the friction coefficient, while wear reduces fault stresses and average shear resistance. The combined effect of local friction and roughness is nonlinear, involving a direct contribution to shear resistance from roughness, as observed in prior studies, and an additional contribution from the increasing average normal stress with slip, which amplifies the resistance from friction and becomes more pronounced at large slip.

Predicting mechanical failure of polycrystalline dual-phase nickel-based alloys by numerical homogenization using a phase field damage model

Brazing of nickel-based alloys plays a major role in the assembly of turbine components, e.g., abradable sealing systems. In a brazed joint of nickel-based alloys a composition of brittle and ductile phases can be formed if the brazing conditions are not ideal. This heterogeneous microstructure is a crucial challenge for predicting the damage behavior of a brazed joint. The initiation and evolution of microdamage inside of the brittle phase of a virtual dual-phase microstructure representing the material in a brazed joint is studied by means of numerical simulations. A phase field approach for brittle damage is employed on the microscale. The simulation approach is capable of depicting phenomena of microcracking like kinking and branching due to heterogeneous stress and strain fields on the microscale. No information regarding the initiation sites and pathways of microcracks is needed a priori. The reliability of calculating the effective critical energy quantities as a microstructure-based criterion for macroscopic damage is assessed. The effective critical strain energy density and the effective critical energy release rate are evaluated for single-phase microstructures, and the approach is transferred to dual-phase microstructures. The local critical strain energy density turns out to be better suited as a model input parameter on the microscale as well as for a microstructure-based prediction of macroscopic damage compared to a model employing the energy release rate. Regarding the uncertainty of the model prediction, using the effective critical energy release rate leads to a standard deviation which is five times larger than the standard deviation in the predicted effective critical strain energy density.

Impact of soft minerals on crack propagation in crystalline rocks under uniaxial compression: A grain‐based numerical investigation

AbstractVarying external conditions in the metallogenetic process of crystalline rocks contribute to the complex mineral and textural characteristics, rendering the mechanical properties highly heterogeneous at the mineral scale. This research focused on the influences of minerals with relatively low strength and stiffness (soft minerals) in crystalline rocks on their cracking behavior. A particle‐based discrete element method (DEM) was first employed to establish random grain‐based models (GBM) of crystalline rocks containing soft minerals with different contents, strengths, and stiffnesses. On this basis, responses of the mechanical properties and crack propagation to these parameters were systematically investigated and an optimized micro‐parameter calibration method for real CT (computed‐tomography) ‐based GBM was then proposed considering the characteristics of soft minerals. The results demonstrate that with the decrease of the strength of soft mica, the intragranular cracks are prone to initiate and propagate inside the soft minerals and lead to the final crack coalescence. The decrease in the stiffness of soft minerals enhances their controlling effects on the cracking propagation. Based on the CT‐based granite model, it was found that a heterogeneous stress field is produced due to the spatial distribution of the soft (mica) and hard (quartz and feldspar) minerals, and the mica minerals tend to terminate cracks or force cracks to deflect or bypass individual grains during crack propagation. This study sheds light on the damage and failure processes of crystalline rocks with contrast mineral components.

Effect of pore on the deformation behaviors of NiTi shape memory alloys: A crystal-plasticity-based phase field modeling

Porous NiTi shape memory alloys (SMA) have been developed into important biomedical and damping materials. However, the effect of pore on the martensitic microstructure evolution, functional properties, and plastic deformation has not been fully clear yet. In this work, to reveal the functional properties including superelasticity (SE), elastocaloric effect (eCE), one-way shape memory effect (OWSME), and stress-assisted two-way shape memory effect (SATWSME), as well as plastic deformation of porous NiTi SMAs, the polycrystalline SMA system is considered as a composite containing the grain-interior (GI) and grain-boundary (GB) phases, a thermo-mechanically coupled, grain-size-dependent and crystal-plasticity-based phase field model is proposed for the GI phase, and an elasto-viscoplastic constitutive model is established for the GB phase. The effects of porosity, pore diameter, pore aspect ratio, pore orientation, and pore distribution, as well as the interaction between pore and grain size, are comprehensively revealed and discussed. The simulations demonstrate that the pores cause a heterogeneous stress field, and the interaction among pores and that between pores and grain boundaries can enhance such heterogeneity, which can promote the initiation of martensitic transformation and reorientation. Moreover, the geometrical constraints caused by the pores can significantly affect the morphology and evolution of martensite microstructures. However, due to the diverse microstructure evolutions during SE, OWSME, and SATWSME, their interactions with the pores are quite different, resulting in various dependencies of each functional property and plastic deformation on the pore, and such dependencies are strongly influenced by porosity, pore diameter, pore aspect ratio, pore orientation, and pore distribution. The simulated results in this work are helpful to comprehensively and deeply understand the pore effect and its microscopic mechanism.

3D Dynamic Rupture Modeling of the 6 February 2023, Kahramanmaraş, Turkey Mw 7.8 and 7.7 Earthquake Doublet Using Early Observations

Abstract The 2023 Turkey earthquake sequence involved unexpected ruptures across numerous fault segments. We present 3D dynamic rupture simulations to illuminate the complex dynamics of the earthquake doublet. Our models are constrained by observations available within days of the sequence and deliver timely, mechanically consistent explanations of the unforeseen rupture paths, diverse rupture speeds, multiple slip episodes, heterogeneous fault offsets, locally strong shaking, and fault system interactions. Our simulations link both earthquakes, matching geodetic and seismic observations and reconciling regional seismotectonics, rupture dynamics, and ground motions of a fault system represented by 10 curved dipping segments and embedded in a heterogeneous stress field. The Mw 7.8 earthquake features delayed backward branching from a steeply branching splay fault, not requiring supershear speeds. The asymmetrical dynamics of the distinct, bilateral Mw 7.7 earthquake are explained by heterogeneous fault strength, prestress orientation, fracture energy, and static stress changes from the previous earthquake. Our models explain the northward deviation of its eastern rupture and the minimal slip observed on the Sürgü fault. 3D dynamic rupture scenarios can elucidate unexpected observations shortly after major earthquakes, providing timely insights for data-driven analysis and hazard assessment toward a comprehensive, physically consistent understanding of the mechanics of multifault systems.

Atomistic study of the impact response of bicontinuous nanoporous gold as a protection medium: Effect of porous-nonporous interface on failure evolution

Compared to their non-porous counterparts, metallic foams are known to exhibit improved functionality (e.g., damping capacity) when subjected to high-speed impact loading. Here, we report the results of a molecular dynamics study of bicontinuous nanoporous gold (NPG) subjected to impact loading. Additionally, we investigate the effects of a heterogeneous impact zone on the fate of a witness specimen initially protected by an NPG target. A cube-shaped flyer object (L0= 408.6 Å) composed of full-density f.c.c. gold (FDG) strikes with impact speed Uflyer= 1.0 km s−1 an initially stationary equal-sized NPG target, which subsequently transmits a dispersed shock wave into a protected FDG witness specimen located at the downstream end of the target. The NPG target is stochastic, with porosity ϕ= 0.5 and mean ligament diameter D¯L=64±6 Å. A corresponding simulation for an FDG target serves as a baseline case. As anticipated, the sharp, planar shock imparted by the FDG flyer rapidly becomes highly curved and broadened along the shock direction in the NPG target. Intense plastic and convective flows (ejecta and jetting) lead to spatially heterogeneous mass flow and energy localization, which persist across the entire length of the NPG target studied. Whereas the shock transmitted by the FDG target leaves the witness specimen intact and essentially undamaged, the heterogeneous stress and flow fields imparted by the NPG target results in destruction of the witness. Independent simulations for the same cube-shaped NPG target shocked on the three statistically equivalent {100} target faces reveal modest run-to-run variability in the target and witness responses. The difference in the effects of the NPG and FDG targets on the witness response motivates additional study of the failure evolution inside the witness. It appears from the preliminary results that the porous-nonporous interface might induce a failure pattern that, in some aspects, resembles failure waves observed in brittle solids.

Strain glass transition in Nb nanowire toughened NiTiHf shape memory alloy composite wires

In this work, we successfully prepared Nb nanowire toughened NiTiHf composite wires with different diameters. The strain glass transition was found in the as-drawn 0.1 mm diameter wire characterized by the Vogel-Fulcher type frequency-dependent modulus anomaly. However, the strain glass transition was not detected in 0.3 mm diameter wires. Microstructural characterization reveals that densely distributed Nb nanowires with small diameter and discontinuous orientation impose both physical confinement and heterogeneous stress fields on the matrix, which interrupt the long-range ordered martensitic transformation, leading to a continuous, high-order like strain glass transition. This work provides a new design strategy for searching new strain glass systems.

Simulated Rupture Dynamics and Radiated Energy on Heterogeneously Damaged Faults

AbstractRupture dynamics along a heterogeneous fault is studied in the framework of the recently developed damage‐breakage rheological model. The model utilizes a fault structure accommodating shear and wear in a heterogeneous weak fault‐core a few centimeters thick, separating two blocks with damage level gradually decreasing toward the host rock. We first demonstrate the similarity of the static stress field around a fault, represented by a three‐body tribosystem, to the one developed around a rough frictional interface. We show that both models predict heterogeneous stress field pattern. We then apply simulations of rupture on heterogeneously damaged faults. We show that increasing initial heterogeneity amplitudes is associated with smaller events with lower slip rates. The simulations further allow us to quantify the amount of accumulated damage correlative with wearing. During the propagating rupture, the strength of the fault‐core evolves, leading to higher wear generation along relatively strong zones or barriers. The total wear production in a given event is strongly dominated by the initial damage heterogeneity. Processing of synthetic seismograms shows excess energy radiation of high‐frequency seismic waves comparing to the expected radiation from planar faults. This radiation is enhanced with the increase in the variability of fault strength heterogeneity. Further calculations of the scaled energy show good fit with previous seismological and laboratory observations and demonstrate the impact of fault heterogeneity on radiated energy during dynamic rupture. Therefore, initial fault heterogeneity, manifested by fault‐core strength or by geometrical irregularity controls many aspects of earthquake rupture, including slip displacement and velocity.

Heterogeneity Versus Anisotropy and the State of Stress in Stable Cratons: Observations From a Deep Borehole in Northeastern Alberta, Canada

AbstractGeophysical logs collected from a deep borehole drilled into the North American Craton in Northeastern Alberta shed valuable light on the state of stress in stable cratons. Observed breakout (BO) azimuths rotate between three depth intervals, ranging from 100° at 1,650–2,000 m to 173° at 2,000–2,210 m, and finally to 145° at the bottom. No apparent fractures that might have disturbed the stress field were found. Two potential causes of these rotating BOs can be either a heterogeneous stress field or elastic and strength anisotropy of the formation. The latter interpretation is favored since observed BO azimuths are strongly controlled by rock metamorphic textures, as validated by their close correlations with dip directions of foliation planes and polarization directions from dipole sonic logs. Monte Carlo realizations further demonstrate that anisotropic metamorphic rocks subjected to a uniform horizontal stress direction could result in the observed azimuth‐rotating BOs. Stress magnitudes inferred from this analysis, which incorporates rock anisotropy and weak foliation failure planes, suggest a normal faulting regime and a maximum horizontal compressive direction consistent with that in the overlying Western Canadian Sedimentary Basin (NE‐SW) and the motion of the North American plate. The inferred stress magnitudes are low and Mohr‐Coulomb analyses demonstrate that the formation is not near the critical condition for slip on weak planes. However, more detailed investigations should be conducted since Monte Carlo calculations indicate that analyses from BO widths, particularly when a conventional Kirsch‐based formula is employed, are highly nonunique, allowing for large variations in potential stress states.

Tectonic stress changes related to plate spreading prior to the 2021 Fagradalsfjall eruption in SW Iceland

The Reykjanes Peninsula in SW Iceland is a part of the Mid-Atlantic plate boundary. It forms its transtensional segment with several volcanic and faulting systems. We focus on the 2017 seismicity that occurred in the central part of Reykjanes at the place of the Fagradalsfjall volcano prior to its eruption on March 19, 2021. We invert well-determined focal mechanisms of the 2017 seismicity and provide mapping of tectonic stress in space and time. Our results disclose heterogeneous stress field manifested by mix of shear, tensile and compressive fracturing. The prominent stress direction was in the azimuth of 120° ± 8°, which represents the overall extension related to rifting in the Reykjanes Peninsula. The activity initiated on the transform fault segment with predominantly shear strike-slip events. The non-shear fractures occurred later being associated with normal dip-slips and corresponding to the opening of volcanic fissures trending in the azimuth of 30–35°, perpendicular to the extension. The dip-slips were mainly located above an aseismic dike detected in the centre of the 2017 swarm. This dike represents a zone of crustal weakening during a preparatory phase of the future 2021 Fagradalsfjall volcanic eruption located at the same place. Moreover, we detected local variation of stress when the stress axes abruptly interchanged their directions in the individual stress domains. These stress changes are interpreted in a consequence of plate spreading and upcoming fluid flow during a preparatory phase of a rifting episode.

Three-dimensional printing of biomimetic variable stiffness composites with controlled orientations and volume fraction of fibers

In this paper, biomimetic continuous fiber reinforced composite structures were modeled and printed on the basis of variable distributions of fiber orientation and the fiber volume fraction, which are similar to the wood structure in the vicinity of a knot. The variable distributions of curved fibers in the variable stiffness composites (VSCs) without gaps and overlaps were realized during the 3D printing process by the control of positioning a nozzle and the adaptive feed of plastic. To demonstrate efficiency of the modeling method and the 3D printing process, a VSC plate with a hole under tensile loading was designed, printed and tested. The results of tensile tests showed that the main failure of the composite plates occurs far from the hole. Furthermore, numerical and experimental strain fields obtained by the finite element method and digital image correlation, respectively, are similar. The efficiency in fiber savings for the composite plates was improved by the change in reinforcement from unidirectional fibers to curved fibers adapted to heterogeneous stress fields.

Prince Rupert’s Drops: An analysis of fragmentation by thermal stresses and quench granulation of glass and bubbly glass

When volcanic eruptions involve interaction with external water (hydrovolcanism), the result is an ash-rich and energetic volcanic plume, as illustrated dramatically by the January 2022 Tonga eruption. The origin of the high explosive energy of these events remains an important question. We investigate this question by studying Prince Rupert's Drops (PRDs)-tadpole-shaped glass beads formed by dripping molten glass into water-which have long fascinated materials scientists because the great strength of the head contrasts with the explosivity of the metastable interior when the tail is broken. We show that the fragment size distribution (FSD) produced by explosive fragmentation changes systematically with PRD fragmentation in air, water, and syrup. Most FSDs are fractal over much of the size range, scaling that can be explained by the repeated fracture bifurcation observed in three-dimensional images from microcomputed tomography. The shapes of constituent fragments are determined by their position within the original PRD, with platey fragments formed from the outer (compressive) shell and blocky fragments formed by fractures perpendicular to interior voids. When molten drops fail to form PRDs, the glass disintegrates by quench granulation, a process that produces fractal FSDs but with a larger median size than explosively generated fragments. Critically, adding bubbles to the molten glass prevents PRD formation and promotes quench granulation, suggesting that granulation is modulated by heterogeneous stress fields formed around the bubbles during sudden cooling and contraction. Together, these observations provide insight into glass fragmentation and potentially, processes operating during hydrovolcanism.

Effect of natural defects on the fracture behaviors and failure mechanism of basalt through mesotesting and FDEM modeling

As an important geological medium for engineering construction in Southwest China, basalt typically contains natural defects, such as hidden joints and amygdales, which significantly affect its mechanical properties. In this work, the fracture behaviors of basalt with different defects were studied in terms of its strength, deformation and failure characteristics, and crack evolution laws using a comprehensive mesotesting scheme and finite-discrete element method (FDEM) modeling, revealing the failure mechanism of basalt. The results showed that: (1) Intact specimens exhibiting high strength and smooth stress–strain curves burst instantly when loaded to the peak strength, and the corresponding acoustic emission events are abnormally concentrated, indicating a fragmentation failure mode. (2) The strength of fractured samples decreases significantly, the stress–strain curve is bimodal or multimodal, and the acoustic emission events are very active during the entire failure process, which is closely related to the local damage generated during each significant stress drop, with tension dominant, shear dominant, or tension–shear composite failure modes. (3) For amygdaloidal basalt samples with a serrated stress–strain curve, each small-scale stress drop is accompanied by the initiation, propagation, and coalescence of microcracks, resulting in relatively scattered acoustic emission events, indicating a splitting failure mode. (4) Because of the small particle size, dense arrangement, good uniformity, and high strength at the grain scale, intact basalt exhibits a relatively uniform stress field easily induced under loading, leading to an extremely high strength, fragmentation failure mode, and significant brittleness. However, a heterogeneous stress field near the original defects, such as hidden joints or amygdales, is prominent, and it significantly affects the fracture paths, which is the main reason for the strong randomness of the mechanical properties, evident strength reduction, and significant change in the failure modes. The research results can provide reference for revealing the occurrence mechanism, guiding warning systems, and formulating control methods for high-stress disasters in hard brittle rocks.

Aftershock distributions, moment tensors and stress evolution of the 2016 Iniskin and 2018 Anchorage Mw 7.1 Alaskan intraslab earthquakes

SUMMARY We present a detailed study of two Mw 7.1 intraslab earthquakes that occurred in southern Alaska: the Iniskin earthquake of 24 January 2016, and the Anchorage earthquake of 30 November 2018. We have relocated and recovered moment tensors for hundreds of aftershocks following both events, and inverted for stress histories. The aftershock distribution of the Iniskin earthquake suggests that the rupture propagated updip along a fault dipping steeply into the Pacific Plate and terminated at a stratigraphic horizon, inferred to be either the interface or Moho of the subducting slab. In addition, four earthquakes ruptured the main fault in the preceding two years and had similar moment tensors to the mainshock. This evidence suggests that the mainshock likely reactivated a pre-existing, outer-rise fault. The Anchorage earthquake sequence is complex due to its location near the boundary of the subducting Yakutat and Pacific plates, as evidenced by the aftershock distribution. Aftershock hypocentres form two main clusters that appear to correspond to orthogonal, conjugate faults, consistent with the two nodal planes of the dominant focal mechanisms. Both geographic groups display many focal mechanisms similar to the mainshock, which could indicate simultaneous rupture on conjugate planes. The time dependence in stress ratio for the Iniskin sequence can be interpreted in terms of pore-pressure evolution within the mainshock fault zone. In particular, our observations are consistent with a dehydration-assisted transfer mechanism where fluids are produced during rupture through antigorite dehydration and raised to high pore pressures through matrix collapse and/or thermal pressurization. The Anchorage sequence exhibits a more complex stress ratio evolution that may be associated with stress adjustments within a distributed fault network, or reflect a strongly heterogeneous stress field.

Robustness of specimen design criteria for identification of anisotropic mechanical behaviour from heterogeneous mechanical fields

To reduce the number of tests required to characterize anisotropic elastoplastic constitutive models, an approach is to design specimen geometries to diversify the stress states generated in a single test. The experiments are then processed using an inverse identification method based on full-field measurements to achieve the full potential of that specific test. Recent optimization methods were able to design complex specimens in which highly heterogeneous stress fields were generated. However, the specimen design is only assessed based on numerical simulations and does not consider the effect of the biases introduced by the full-field measurement method. The goal of this work is therefore to take into account some of the most frequently observed measurement biases in the specimen selection process. The proposed approach uses synthetic test images generated with numerical simulations. Four specimen geometries have been ranked based on two selection criteria. The first one is an indicator of the heterogeneity of the stress fields obtained by finite element simulations (unbiased data). The second one quantifies error for the identification procedure due to measurement biases. The two criteria provide different rankings for the set of specimens. It is concluded that the design with the most heterogeneous stress fields (first criterion) is not necessarily the more robust design in terms of measurement noise (second criterion), so the optimized geometry should be selected based on a compromise between these two criteria.

Approach to failure through record breaking avalanches in a heterogeneous stress field

We study how the competition of the disordered local strength and of the evolving inhomogeneous stress field affects the evolution of the series of breaking avalanches accompanying the fracture of heterogeneous materials. To generate fracture processes, we use a fiber bundle model of localized load sharing where the degree of strength disorder is controlled by varying two parameters of the distribution of the breaking threshold of fibers. Analyzing the record statistics of avalanches of breaking fibers, we demonstrate that both for low and high disorders the series of crackling events remains stationary until global failure making the collapse of the system unpredictable. Based on computer simulations, we determine a region of the parameter plane of strength disorder where global failure is preceded by an accelerating breaking activity. We show that the record avalanche with the longest lifetime can be used to identify the onset of acceleration of the fracture process towards the catastrophic failure. Comparison of the results to their equal load sharing counterparts reveals that the accelerating regime is shorter than in case of a homogeneous stress field due to the higher degree of brittleness of the system caused by stress localization.

Preliminary study of stress changes evolution on central part of Sumatran Fault influenced by large interplate earthquakes and tectonic stress rates

The inland seismic activity in Great Sumatran Fault (GSF) has significantly increased over the past several decades after the occurrence of historical large interplate earthquakes along the plate boundary. This condition led to some occurrences of historical intraplate earthquakes along Sumatran fault. To quantitatively examine the physical mechanisms between intraplate earthquakes and interplate earthquakes, we estimated the static coseismic stress changes of Coulomb failure function (ΔCFF) using receiver fault approach from large historical-recorded interplate earthquakes and the increase in tectonic stress rates. We examined this research in the central part of GSF since this zone is assumed to have the most heterogeneous stress field and thus became our focus study area. The cumulative ΔCFF models showed almost all segments in the central part of GSF suffered negative changes (<-0.1 MPa) which assumed to be unlikely to rupture in short time. However, the preliminary analysis of the increase in tectonic stress rate indicated that large intraplate earthquakes occurred on Angkola and Siulak segments were dominantly influenced by the increase in interseismic stress rate just after the series of large subduction earthquake occurrences, apart from the decreased stress changes from those major interplate earthquakes.