Stress Field Research Articles

Singular stress field induced two-dimensional topological polar structures in SrTiO3: Topological number manipulation via loading modes

Topological structures in condensed matter attract much attention of the scientific community due to their unique and significant properties. In contrast to the well-explored topological spin structures such as meron, and skyrmion in ferromagnetic materials, ferroelectrics face challenges in forming topological polar structures due to the absence of a spontaneous driving force. Nevertheless, the coupling effect between polarization and strain (gradient), e.g. piezoelectricity and flexoelectricity, provides a mechanical approach to induce topological polar structures that reflect the intricate features of the strain (gradient) field. This study demonstrates the creation of a topological polar structure through a mode II crack-induced singular stress field in SrTiO3 by using phase-field simulations. Furthermore, under mixed-mode loading of mode I and mode II, alterations in the mode ratio induce a variety of topological structures, including a novel structure (winding number nv=4) that is absent under the single-mode loading. Our study demonstrates that topological strain field can be a general and effective tool to generate and control topological polar structures.

Research on the Interaction Mechanism of Multi-Fracture Propagation in Hydraulic Fracturing

During the hydraulic-fracturing process, stress interference occurs among multiple wells and fractures, potentially affecting the trajectory of hydraulic fracture propagation. Previous studies have largely overlooked the influence of proppant support stresses on the trajectories of fracture propagation. This paper establishes a mathematical model, grounded in the boundary element method, designed to compute the propagation of multiple fractures, considering both proppant support on the fracture surface and dynamic perturbations within the local stress field. The findings of this research reveal that the stress field induced by hydraulic fracturing exhibits dynamic evolution characteristics, necessitating a comprehensive study of the fracture initiation and extension across the entire fracturing time domain. The effect of the residual fracture width under proppant action on the in situ stress field cannot be ignored. During simultaneous fracturing, hydraulic fractures are inclined to propagate in the direction of the maximum horizontal principal stress, particularly as the in situ differential stress escalates. Staggered fracturing between wells has been proven to be more effective than head-to-head fracturing. Simply increasing the well spacing cannot solve the problem of inter-well fracture interference. In zipper fracturing, adjusting the fracturing sequence can inhibit the fracture intersections between wells, thereby controlling the trajectory of fracture propagation. The aforementioned research has considerable significance in guiding the control of fracture morphology during hydraulic-fracturing processes.

4D Electrical Resistivity Imaging of Stress Perturbations Induced During High‐Pressure Shear Stimulation Tests

AbstractFluid flow through fractured media is typically governed by the distribution of fracture apertures, which are in turn governed by stress. Consequently, understanding subsurface stress is critical for understanding and predicting subsurface fluid flow. Although laboratory‐scale studies have established a sensitive relationship between effective stress and bulk electrical conductivity in crystalline rock, that relationship has not been extensively leveraged to monitor stress evolution at the field scale using electrical or electromagnetic geophysical monitoring approaches. In this paper we demonstrate the use time‐lapse 3‐dimensional (4D) electrical resistivity tomography to image perturbations in the stress field generated by pressurized borehole packers deployed during shear‐stimulation attempts in a 1.25 km deep metamorphic crystalline rock formation.

Investigation of the Strain–Stress Field in Nanoscale Multilayer Systems by the Phase Plane Method

This paper presents the results of the study of stress relaxation fields, deformation, and temperature of the system of nanostructured multilayer coatings. In the work, a nonlinear relationship between strain and stress was used to take into account nonlinear effects in the mechanism of nanostructure formation. The paper assumes that a friction surface is provided by the self-organization of shear components: both stress and strain on the one hand, and temperature on the other. The studied objects are described in the adiabatic approximation, taking into account the fact of the evolution of stresses and strains. With the help of phase portraits of the system, the dependence of the deformation processes on the stresses arising in the system without coating and with coating is shown. It is shown that the rate of change of deformation depends on the characteristics of the mechanical impact on the coating and on the amount of stress and deformation. A conclusion is drawn regarding the transition process in the presence of two regions (Hooke and plastic deformation) in the corresponding phase portrait of the strain–stress field of the system. The results of the work can be used to determine the effective parameters of a coating in the analysis of experimental time dependences of stresses.

R-Phase Transformation Evolution in NiTi SMA Wires Studied via the Internal Friction Technique

The specific damping capacity variation of heat-treated NiTi was observed during a pseudoelasticity test. The detailed B2 → R-phase transformation process in cold-drawn NiTi wires undergoing middle-temperature aging was studied via X-ray diffraction, transmission electron microscope, and internal friction technique. Results show that, as aging time increased at 450 °C, the dynamic phase transition splitting from B2 → R to B2 → R1 and B2 → R2 became evident. However, such a splitting process was not observed for the sample after aging at 400 °C. The reason for R-phase generation is attributed to non-uniformly distributed stress fields. The splitting of the internal friction peak, in conjunction with high-resolution transmission electron microscope and mechanic results, suggests a substantial occurrence of short-range segregation of Ni atoms in the B2-NiTi matrix. Furthermore, the specific damping capacity (SDC) exhibits a gradual increase with prolonged annealing time. Specifically, the sample with significant dynamic phase transition splitting reaches an SDC value of 0.60.

Application of laboratory micro X-ray fluorescence devices for X-ray topography

It is demonstrated that high-resolution energy-dispersive X-ray fluorescence mapping devices based on a micro-focused beam are not restricted to high-speed analyses of element distributions or to the detection of different grains, twins and subgrains in crystalline materials but can also be used for the detection of dislocations in high-quality single crystals. Si single crystals with low dislocation densities were selected as model materials to visualize the position of dislocations by the spatially resolved measurement of Bragg-peak intensity fluctuations. These originate from the most distorted planes caused by the stress fields of dislocations. The results obtained by this approach are compared with laboratory-based Lang X-ray topographs. The presented methodology yields comparable results and it is of particular interest in the field of crystal growth, where fast chemical and microstructural characterization feedback loops are indispensable for short and efficient development times. The beam divergence was reduced via an aperture management system to facilitate the visualization of dislocations for virtually as-grown, non-polished and non-planar samples with a very pronounced surface profile.

Induced Stress Fields in Isotropic Elastic Discs Due to Contact Shear Stresses and Their Implications for Brazilian Test Strength Characterization

Induced Stress Fields in Isotropic Elastic Discs Due to Contact Shear Stresses and Their Implications for Brazilian Test Strength Characterization

Building a carbon dioxide storage portfolio for the Barrow-Dampier sub-basins through regional screening – an integrated geoscience approach

Safe geological carbon dioxide (CO2) storage (GCS) requires rocks with suitable injectivity, capacity and sealing properties to ensure secure long-term containment of injected CO2. A regional understanding of the subsurface is essential to determine the potential for GCS of a basin and to select target sites. This is best addressed by integrating the basin’s tectono-stratigraphic evolution, its gross depositional environment, and its hydrodynamic, thermal and stress regimes. A basin-scale GCS assessment for the Barrow-Dampier sub-basins was conducted by Geoscience Australia and SLB. The objective of the study is to high-grade geological intervals and sites for potential GCS and to understand potential storage capacity and key risk factors. Stratigraphic and structural mapping of key storage intervals was performed using the reprocessed seismic volume and well database associated with the project. Analysis of critically stressed faults was used to estimate the likelihood of reactivation based on the far-field regional stress field and fault mechanical properties. Pressure, temperature, porosity, permeability, and water geochemistry data have been screened for >500 wells for assessing the storage unit intervals and predicting the hydrodynamic regime. Calibrated 2D basin models provide information on the regional pressure-temperature regime, porosity/permeability distribution, and sealing effectiveness. Potentiometric surface maps for the aquifer systems inform the distribution of potential CO2 plume migration. Results of this integrated regional basin study are used to quantify the risk of identified storage containers and to map the chance of success for GCS at a regional scale. The project results are publicly available via NOPIMS.

Fatigue Life Prediction and Verification of Railway Fastener Clip Based on Critical Plane Method

Accurately predicting the fatigue life of fastener clips is crucial for improving material processing technology, enhancing the anti-fatigue design level, and guiding the maintenance and repair of track structures. Conventional methods that solely evaluate the fatigue life of fastener clips based on the uniaxial stress index under displacement loads may lead to significant errors. In this paper, the stress field of a type-III clip under displacement loading conditions was numerically simulated based on fatigue test standards, and the fatigue life of the clip was analyzed using the Fatemi–Socie (FS) multiaxial fatigue criterion based on the critical plane method. A comparison with standard fatigue test results revealed that, under non-resonance conditions, the predicted position of the fatigue critical plane of the fastener clip coincided with the fracture surface observed at the middle of the measured small arc, with a life prediction error of 7.3%. To further investigate the predictive capability of the FS multiaxial fatigue criterion for the fatigue life of the fastener clip under resonance conditions, the stress levels of the clip under non-resonance and resonance conditions were compared through on-site testing, and the effects of inertial loads caused by vertical and lateral vibration acceleration on the stress field of the clip were analyzed in numerical simulation according to the results of clip acceleration tests; the fastener clip’s resonance fracture position and fatigue life were also predicted based on the aforementioned multiaxial fatigue criterion. To verify the accuracy of the predicted results, a testing method was proposed that equates the high-frequency resonant inertial load of the fastener clip to a low-frequency additional preload under the premise of consistent stress fields. A comparison with the numerical simulation results shows that considering only vertical inertial loads would result in a discrepancy between the measured fracture location and the actual one. Considering both vertical and lateral inertial loads, the fracture location (at the heel of the small arc) under resonance conditions could be accurately determined, with a life prediction error of 13.8%. Compared to the non-resonant displacement loading condition, the inertial loads caused by acceleration under resonance conditions led to a reduction in fatigue life of approximately 77.8%.

Image analysis of acoustic data and interpretation of rock stress orientations for geothermal exploration in Gothenburg borehole GE-1, SW Sweden

Abstract This study contributes to geothermal exploration in 1660–1520 Ma old, reworked bedrock in Sweden. Our primary objectives are to constrain the orientation of horizontal stresses, and to discuss implications for geothermal exploration. High-resolution acoustic televiewer image data reveal the downhole distribution of stress indicators (borehole breakouts, drilling-induced fractures and petal centreline fractures) and pre-existing structures (natural fractures, foliation). About 135 m of stress indicators are measured from 0.2–1.0 km. The results suggest a uniform NNW–SSE mean maximum horizontal stress orientation. A total of 1525 pre-existing structures (natural fractures, foliation) are mapped in borehole GE-1. The prevailing stress regime controls whether natural fractures and foliation are well-oriented for stimulation. For strike-slip and normal faulting stress regimes, well-oriented fractures steeply dip towards the WSW. For a reverse faulting stress regime, shallow dipping fractures are well-oriented for stimulation. The downhole distribution of stress indicators and other stress measurements in the region and other parts of Fennoscandia tentatively suggest a strike-slip stress regime, but additional studies are needed to constrain the complete stress field at study depth and towards engineered geothermal systems reservoir target depths. Our secondary objective is to highlight that interpretation of high-resolution acoustic data, particularly in metamorphic crystalline rocks, is subjective and that more guidelines for data interpretation are needed. The interactive interpretation of the images is based on visual analyses of complex pre-existing structures and stress indicators with highly variable shapes. The application of three methods for data analyses in the GE-1 borehole proposes that drilling-induced fractures are little influenced by the method applied. Interpretations on individual borehole breakout azimuths may, however, result in over 10° differences in orientation.

Fabricating high-performance biomedical PLLA/PVDF blend micro bone screws through in situ structuring of oriented PVDF submicron fibers in microinjection molding

Poly(l-lactide) (PLLA) is regarded as a polymer with excellent biocompatibility, but its inherent brittleness property greatly restricts its application in the biomedical engineering field. Blending PLLA with other polymers is one of the sufficiently viable methods of property improvement. In this paper, the Poly(l-lactide) (PLLA)/Polyvinylidene fluoride (PVDF) blend was first prepared by melt-compounding, and then microinjection molded into the blend micro bone screw with high strength and toughness. The experimental results show that the PVDF dispersed phases could in situ form the highly oriented fibers with shish-kebab structures under the effect of strong shear stress field of microinjection molding along flow direction. These parallelly aligning PVDF fibers could simultaneously realize the significant enhancement in both toughness and strength of the micro bone screw. Specifically, in the bending test the blend micro bone screw with oriented PVDF fibers can withstand the compression displacement of nearly 2 mm, while the neat PLLA one would break within only 0.5 mm displacement. In addition, a bending recovery phenomenon was observed in the load-displacement curve of the blend, demonstrating the mechanical transition from brittle fracture to ductile fracture. Moreover, with addition of 30 % PVDF, the elongation at break increases from 7.8 % to 57.8 %, while the tensile strength increases from 60.9 MPa to 74.3 MPa. The PLLA/PVDF micro screws possess better toughness and excellent biocompatibility compared to pure PLLA products. They also exhibit great potential for stimulating the proliferation of cells through piezoelectric output, opening up new possibilities for the development of next-generation fracture fixation materials.

An extended graphical solution for undrained cylindrical cavity expansion in K0‐consolidated Mohr–Coulomb soil

AbstractThis paper develops a general and complete solution for the undrained cylindrical cavity expansion problem in nonassociated Mohr‐Coulomb soil under nonhydrostatic initial stress field (i.e., arbitrary values of the earth pressure coefficient), by expanding a unique and efficient graphical solution procedure recently proposed by Chen and Wang in 2022 for the special in situ stress case with . It is interesting to find that the cavity expansion deviatoric stress path is always composed of a series of piecewise straight lines, for all different case scenarios of K0 being involved. When the cavity is sufficiently expanded, the stress path will eventually end, exclusively, in a major sextant with Lode angle θ in between and or on the specific line of . The salient advantage/feature of the present general graphical approach lies in that it can deduce the cavity expansion responses in full closed form, nevertheless being free of the limitation of the intermediacy assumption for the vertical stress and of the difficulty existing in the traditional zoning method that involves cumbersome, sequential determination of distinct Mohr–Coulomb plastic regions. Some typical results for the desired cavity expansion curves and the limit cavity pressure are presented, to investigate the impacts of soil plasticity parameters and the earth pressure coefficient on the cavity responses. The proposed graphical method/solution will be of great value for the interpretation of pressuremeter tests in cohesive‐frictional soils.

Characterization in the rat of the individual tendency to rely on alcohol to cope with distress and the ensuing vulnerability to drink compulsively

Abstract Only some vulnerable individuals who recreationally drink alcohol eventually develop the compulsive drinking pattern that characterizes alcohol use disorder. A new frontier in biomedical research lies in understanding the neurobehavioral mechanisms of this individual vulnerability, a necessary step towards developing novel effective therapeutic strategies. Translational research has been hindered by the lack of valid, reliable and robust approaches that enable the study of the influence of the reliance on alcohol to cope with stress or self-medicate negative emotional states on the subsequent transition to alcohol use disorder. We have therefore developed a behavioral task in the rat that enables the investigation of the neural and cellular basis of the exacerbation of the vulnerability to develop compulsive alcohol drinking by the use of alcohol to develop an adjunctive, anxiolytic, polydipsic drinking behavior in a schedule-induced polydipsia procedure. Hence, in our task, alcohol is introduced in the schedule-induced polydipsia context after several weeks of training with water so that rats are exposed to alcohol for the first time in a distressing context and learn to drink alcohol as a coping strategy. Capitalizing on this protocol, we have consistently been able to identify a subpopulation of rats that were unable to learn to cope with negative states by drinking water and relied on alcohol to do so. This maladaptive reliance on alcohol-drinking to cope with distress has been shown to be associated with an exacerbation of the subsequent transition to compulsive drinking. Furthermore, these vulnerable rats reached blood alcohol levels comparable to that of intoxication in humans, thereby developing two key features of alcohol use disorder, namely excessive alcohol intake and compulsive drinking. Altogether, this behavioral task provides a novel and unique tool for the investigation of the neurobehavioral mechanisms underlying the exacerbation of the individual vulnerability to developing compulsive alcohol drinking by the use of alcohol as a strategy to cope with distress, and for the evaluation of the efficacy of potential therapeutic strategies in a personalized medicine approach. This procedure, which focuses on an understudied but key factor of the development of alcohol use disorder, may become widely used as it benefits the fields of alcohol, emotion regulation and stress, the interest in which has substantially increased since the evidence of a profound exacerbation of alcohol use and alcohol-related negative consequences by the distress and social isolation engendered by the various measures implemented worldwide in response to the COVID-19 pandemic.

A new method for determining the ogden parameters of soft materials using indentation experiments

A full understanding of the material properties of skin tissue is crucial for exploring its tribo-mechanical behaviour. It has been widely accepted that the mechanical behaviour of skin tissue for both small and large deformations can be accurately described using a hyperelastic model, such as the one developed by Ogden. However, obtaining these Ogden parameters for in-vivo skin by in-vivo experiments no matter the indentation or suction tests is a significant challenge. The mathematical model used to describe the material behaviour during the test should consider not only the material nonlinearity but also the geometrical confinement of the tissue, the large deformations induced, and the fact that the specimens are relatively thin.A range of contact models is available to describe the contact behaviour during the indentation test. However, none of them can be used for hyperelastic materials with small thickness under large deformations. Simultaneously explaining material nonlinearity and geometric nonlinearity, either through theoretical equations or numerical calculations, poses a significant challenge.In this research, we propose a pragmatic method to obtain Ogden parameters for in-vivo skin tissue by combining experimental indentation results and numerical simulations. The indentation tests were used to obtain the force-indentation depth curves, while the numerical simulations were used to obtain the strain fields. The method assumes the material behaviour of specimens can be linearized in each small deformation increment, and the contact model developed by Hayes can be applied to accommodate each increment. Then, the linear elastic behaviour in each increment can be described by the elastic modulus E which were obtained using Hayes model, and the principal stresses in each increment were subsequently obtained using Hooke's law. By combining all stress fields, overall stress-strain curves can be constructed, from which the hyperelastic Ogden parameters can be obtained. A second numerical simulation of the hyperelastic indentation was then performed using the obtained Ogden parameters, allowing a comparison of the experimental and simulated relationships between force and indentation.

Deep learning as a powerful tool in digital photoelasticity: Developments, challenges, and implementation

Stress field evaluation, through fringe order maps, has always been of great importance in various engineering domains, providing essential insights into the mechanical response of loaded structures and enabling the prevention of potential failures, thus ensuring the safety and reliability of critical components. In recent years, the advent of new algorithms has ushered in a transformative era, particularly through computational techniques such as machine learning and, in particular, deep learning. These advanced approaches have proven instrumental in optimizing fringe demodulation processes, providing more accurate predictions, and uncovering complex patterns within structural behavior. This paper reports a careful review of papers and articles discussing the application of deep learning to photoelasticity studies. It explores the impact and applicability of deep learning in the field of digital photoelasticity, highlighting its benefits and advancements. In addition, the challenges associated with the implementation of this technique are examined to provide a critical insight. The paper also discusses experimental results and possible suggestions for future work in this new paradigm.

Real‐time topology optimization via learnable mappings

AbstractIn traditional topology optimization, the computing time required to iteratively update the material distribution within a design domain strongly depends on the complexity or size of the problem, limiting its application in real engineering contexts. This work proposes a multi‐stage machine learning strategy that aims to predict an optimal topology and the related stress fields of interest, either in 2D or 3D, without resorting to any iterative analysis and design process. The overall topology optimization is treated as regression task in a low‐dimensional latent space, that encodes the variability of the target designs. First, a fully‐connected model is employed to surrogate the functional link between the parametric input space characterizing the design problem and the latent space representation of the corresponding optimal topology. The decoder branch of an autoencoder is then exploited to reconstruct the desired optimal topology from its latent representation. The deep learning models are trained on a dataset generated through a standard method of topology optimization implementing the solid isotropic material with penalization, for varying boundary and loading conditions. The underlying hypothesis behind the proposed strategy is that optimal topologies share enough common patterns to be compressed into small latent space representations without significant information loss. Results relevant to a 2D Messerschmitt‐Bölkow‐Blohm beam and a 3D bridge case demonstrate the capabilities of the proposed framework to provide accurate optimal topology predictions in a fraction of a second.

A phase-field model for the brittle fracture of Euler–Bernoulli beams coupling stretching and bending

Damage gradient models seek to simulate fracture mechanics by modulating the material stiffness. Within this framework, a singular scalar field representing damage is commonly utilized to globally decrease elastic energy. However, when considering structural models like beams and plates, this approach often fails to adequately capture important aspects due to the interaction between stretching and bending contributions. In this work, we propose a model for planar Euler–Bernoulli beams that incorporates the following features: firstly, we utilize two phase-field damage parameters to describe material damaging, specifically addressing the ‘erosion’ occurring above and below the original, undamaged beam; secondly, we assume a simple linear dependence of the axial stress field on the thickness coordinate, along with linear dependence on the axial force and bending moment. By appropriately identifying the constitutive response, our model effectively considers the coupling between stretching and bending induced during through-the-thickness damage, demonstrating good agreement with 2D observations.

Geostatistical Analysis of Lineament Domains: The Study Case of the Apennine Seismic Province of Italy

Regional-scale swarms of subparallel linear topographic features, known as lineament domains, are a common feature of planetary surfaces. Lineament domains are superficial manifestations of the crustal stress field trajectory. Notably, one of the effects of active tectonics is seismicity. Italy is one of the most seismically active regions in the Mediterranean, with many destructive earthquakes that have occurred in past centuries. Here, we assess the seismic meaning of the main lineament domain in the tectonically active region of Central Italy. We describe the use of an automated analysis of satellite imagery coupled with spatial grid analysis to identify three lineament domains of the Central Apennines. Spatial and azimuthal comparisons of the main lineament domain (i.e., the Apennine Domain), with the known locations of earthquakes (moment magnitude of Mw > 5.5) that occurred during the past century, revealed the most seismically active tectonic areas and their spatial distributions. Further, we present a conceptual seismo-geodynamic model for the Central Apennines, which is characterized by regional arching and explains the presence of an extensional tectonic regime in the upper crustal layer of the active Apennines fold-and-thrust belt.

Study on the Interaction Propagation Mechanism of Inter-Cluster Fractures under Different Fracturing Sequences

Horizontal-well multi-cluster fracturing is one of the most important techniques for increasing the recovery rate in unconventional oil and gas reservoir development. However, under the influence of complex induced stress fields, the mechanism of interaction and propagation of fractures within each segment remains unclear. In this study, based on rock fracture criteria, combined with the boundary element displacement discontinuity method, a two-dimensional numerical simulation model of hydraulic fracturing crack propagation in a planar plane was established. Using this model, the interaction and propagation process of inter-cluster fractures under different fracturing sequences within horizontal well segments and the mechanism of induced stress field effects were analyzed. The influence mechanism of cluster spacing, fracture design length, and fracture internal pressure on the propagation morphology of inter-cluster fractures was also investigated. The research results indicate that, when using the alternating fracturing method, it is advisable to appropriately increase the cluster spacing to weaken the inhibitory effect of induced stress around the fractures created by prior fracturing on subsequent fracturing. Compared to the alternating fracturing method, the propagation morphology of fractures under the symmetrical fracturing method is more complex. At smaller cluster spacing, fractures created by prior fracturing are more susceptible to being captured by fractures from subsequent fracturing. The findings of this study provide reliable theoretical support for the optimization design of fracturing sequences and fracturing processes in horizontal well segments.

Analysis of Microwave-Induced Damage in Granite Aggregates Influenced by Mineral Texture

The use of microwave energy to recycle high-quality coarse aggregates from waste concrete or assist hard rock breakage in underground building engineering is promising. Controlling or promoting the damage of coarse aggregates, i.e., hard rocks, under microwave irradiation is a crucial issue faced by these techniques. Understanding the damage mechanisms of hard rocks exposed to microwaves is thus urgent. Fracture toughness is a significant mechanical parameter of rocks that reflects their ability to resist crack propagation and damage evolution. In this study, the fracture toughness degradation of microwave-heated granite was investigated by combining experimental investigations and numerical simulations. A three-point-bending (TPB) experiment was conducted on granite specimens after microwave irradiation. A coupled electromagnetic–thermal–mechanical model considering the actual mineral texture of the granite specimen was established. The evolution of the temperature gradient and stress field near the initial notch tip were investigated. The results suggest that the microwave-induced maximum temperature gradient and stress in granite are at the plagioclase–quartz (Pl–Qtz) interfaces or inside the Pl near the boundary. The region of cracking initiation was defined as the damage zone, which could be obtained by comparing the microwave-induced thermal stress with the critical value. The fracture toughness degradation, which corresponds to the evolution of the damage zones, can be divided into two stages. A relatively rapid decrease in fracture toughness in the first stage is primarily caused by the spread of the scattered damage zones along the Pl–Qtz interfaces; subsequently, a gentler fracture toughness degradation results mainly from the extension of the previous damage zones.