Local Stress State Research Articles

Modified Centroid of Root Projection Method for Determining the Center of Resistance of a Tooth

Abstract Center of resistance (CR) has been widely accepted in dentistry as a reference point for controlling tooth moment, which depends on the direction of loading and the morphology of the periodontal ligament (PDL). In clinical practice, dentists estimate the location of CR based on the morphology of the root of teeth, which may lead to a misestimation of orthodontic treatment. A quick method was proposed to efficiently determine the CR by identifying the centroid of the root projection (CRP), according to the orthodontic force. However, the original CRP method was limited to single-rooted teeth and it did not provide a strategy for handling the overlapping roots projection of multi-rooted teeth. To address this issue, we expanded the CRP method to accommodate multi-rooted teeth by applying the weighted average method. We further validated the modified CRP method using finite element simulation (FEA) for both single-rooted and multi-rooted teeth considering mesial-distal and buccal-lingual force directions. The evaluation of displacement distribution along the projection direction allowed us to assess translation and rotation movements, which confirmed that the centroid of root projection can accurately serve as the CR for the multi-rooted teeth. Additionally, we observed heterogeneous stress distributions in the multi-rooted teeth. Considering the well-acknowledged bone remodeling effect in response to local stress states, this indicated that the comprehensive indexes beyond the CR are desired for evaluating or controlling the tooth movement.

The role of irradiation-enhanced interstitial diffusion in over-pressurizing fission gas bubbles in UO2

Fission gas bubbles in UO2 nuclear fuel have been observed to exhibit pressures in excess of the equilibrium bubble pressure; however, the cause of bubble over-pressurization has not yet been demonstrated. The mechanical interaction between a bubble and the surrounding matrix or grain boundary depends on the internal pressure of the bubble and local stress state, such that over-pressurized bubbles are thought to be responsible for fragmentation and pulverization, when exposed to a temperature ramp. Here, we investigate the role of U interstitials, produced through irradiation, in over-pressurizing bubbles by using a combined molecular dynamics (MD) and cluster dynamics approach. Firstly, the energies for the capture of interstitials and vacancies by bubbles have been determined from MD as a function of the ratio of gas atoms to vacancies that make up the bubble. Secondly, these reaction energies have been implemented in the cluster dynamics code Centipede to predict bubble over-pressurization as a function of temperature for typical fission rates. It was found that there is a transition from low pressure bubbles (at high temperatures) to high pressure bubbles (at lower temperatures). The cause of this behavior was shown to be the creation of irradiation-induced interstitials that are highly mobile relative to vacancies at low temperature; whereas, vacancies are sufficiently mobile at high temperatures to limit bubble pressures. This result supports the hypothesis that over-pressurized bubbles form during steady-state operation and that this behavior is highly sensitive to the local pellet temperature.

Optimized Designation of Mesoscopic Configuration of Zr/Ti Layered Metal Composites for Strength-ductility Synergy Improvement

Zirconium (Zr) and its alloys are highly promising for extreme service environments due to their excellent radiation and corrosion resistance. However, a major challenge in Zr-based materials lies in the strength-ductility trade-off, which limits their broader structural applications. This study aims to address this challenge by fabricating Zr/Ti layered metal composites (LMCs) with different layer thickness ratios (LTRs) and investigating the effects of LTR on the mechanical properties of Zr/Ti LMCs. Zr/Ti LMCs with varying LTRs were fabricated to explore how mesoscopic configurations influence their deformation mechanisms. The results reveal that an optimized layer thickness (~50 μm for Zr and ~40 μm for Ti) maximizes the hetero-deformation-induced (HDI) strengthening effect, resulting in a synergy improvement in strength and ductility. Stable and constrained cracks observed in the Zr layers contributed to local stress relief and improved ductility. The anisotropic deformation between the Zr and Ti layers transformed the local stress states into triaxial states. In addition, the localized stress state and inter-layer strain transmission also promote prismatic <a> slip near the interfaces. This study provides key insights into the design of Zr-based composites with improved mechanical performance, offering potential solutions for their application in demanding environments.

Simulation of Localized Stress Impact on Solidification Pattern during Plasma Cladding of WC Particles in Nickel-Based Alloys by Phase-Field Method

As materials science continues to advance, the correlation between microstructure and macroscopic properties has garnered growing interest for optimizing and predicting material performance under various operating conditions. The phase-field method has emerged as a crucial tool for investigating the interplay between microstructural characteristics and internal material properties. In this study, we propose a phase-field approach to couple two-phase growth with stress–strain elastic energy at the mesoscale, enabling the simulation of local stress effects on the solidified structure during the plasma cladding of WC particles and nickel-based alloys. This model offers a more precise prediction of microstructural evolution influenced by stress. Initially, the phase field of WC-Ni binary alloys was modeled, followed by simulations of actual local stress conditions and their impacts on WC particles and nickel-based alloys with ProCAST and finite element analysis software. The results indicate that increased stress reduces grain boundary migration, decelerates WC particle dissolution and diffusion, and diminishes the formation of reaction layers and Ostwald ripening. Furthermore, experimental validation corroborated that the model’s predictions were consistent with the observed microstructural evolution of WC particles and nickel-based alloy composites.

Uncertainty Analysis of Tensile Strength of Scarf Adhesive Joints Using Numerical Method

This study uses three-dimensional finite element (FE) analysis to investigate stress distribution and fatigue strength in adhesive scarf joints, which are complex due to their oblique geometry and varying material properties. ANSYS was used to model stress components within the adhesive and at the adhesive-adherent interface. The FE models, created with eight-node solid elements, utilized fine meshing to ensure accurate stress representation. A mesh size of 0.025 mm was optimal. The models, assuming linearly elastic materials, were validated by comparing simulated deformations with experimental results. Findings enhance understanding of local stress states, improving adhesive joint design reliability. Key Words: Adhesive scarf joints, finite element analysis, stress distribution, fatigue strength, ANSYS,

Accounting for defect position in ultrasonic fatigue test specimen with heterogeneous stress distribution

To account for or neglect the defect position, i.e. the fatigue initiating defect location, both radially and axially, is evaluated for hourglass-shaped ultrasonic fatigue specimen. The commonly used analytical equations to calculate the stress is compared against a finite element (FE) based approach, which is able to fully considering the stress state at the defect position. Notably, the effects on several common fatigue analyses are evaluated: the fatigue strength distribution, the stress–life and the stress–defect relationships. Fracture mechanical assessment is also performed, for a comprehensive VHCF characterization of the EN-GJS-500-7 ductile cast iron used in the study. The VHCF properties are characterized up to 3⋅108 cycles, using the Step-Stress fatigue testing method under fully reversed loading. The FE-model and Weibull distribution as the choice of fatigue strength distribution, enables size effect evaluation by Weakest-link effective volume with the highly stressed volume method as benchmark. The work shows that it is imperative to use the local stress state at the defect position, as the distribution of failures can diverge largely from the center of the specimen, and that neglecting this causes systematic error and flawed potentially results.

Design optimization of advanced tow-steered composites with manufacturing constraints

Tow steering technologies, such as automated fiber placement, enable the fabrication of composite laminates with curvilinear fiber, tow, or tape paths. Designers may therefore tailor tow orientations locally according to the expected local stress state within a structure, such that strong and stiff orientations of the tow are (for example) optimized to provide maximal mechanical benefit. Tow path optimization can be an effective tool in automating this design process, yet has a tendency to create complex designs that may be challenging to manufacture. In the context of tow steering, these complexities can manifest in defects such as tow wrinkling, gaps, overlaps. In this work, we implement manufacturing constraints within the tow path optimization formulation to restrict the minimum tow turning radius and the maximum density of gaps between and overlaps of tows. This is achieved by bounding the local value of the curl and divergence of the vector field associated with the tow orientations. The resulting local constraints are effectively enforced in the optimization framework through the Augmented Lagrangian method. The resulting optimization methodology is demonstrated by designing 2D and 3D structures with optimized tow orientation paths that maximize stiffness (minimize compliance) considering various levels of manufacturing restrictions. The optimized tow paths are shown to be structurally efficient and to respect imposed manufacturing constraints. As expected, the more geometrical complexity that can be achieved by the feedstock tow and placement technology, the higher the stiffness of the resulting optimized design.

Future change of humid heat extremes and population exposure in Turkey

AbstractGlobal climate projections show that humid heat extremes will expand toward the higher latitudes, making the midlatitudes hotspots for these extremes. Therefore, a thorough explanation of their regional characteristics becomes crucial, given that the changes in these extremes can potentially render a large proportion of the global population at risk. Here, we perform the first analysis of historical and projected changes in the intensity and frequency of humid heat extremes and quantify the population exposure to these extremes in Turkey, using long‐term simulations from the non‐hydrostatic mesoscale model of Consortium for Small‐scale Modeling (COSMO‐CLM) under the RCP8.5 emission scenario. We portray not only the nationwide changes in the humid heat extremes and population exposure but also their regional aspects by exploiting the K‐means clustering algorithm. Our results suggest significant future increases in the intensity and frequency of these extremes over a wide geographical area, which includes the surroundings of Adana, Antalya, Izmir, Sakarya, Ordu and Diyarbakir, most of which are coastal locations. Over most of these regions, severe humid heat stress is expected to last nearly a month every year, with almost 56% of the land area is projected to experience local historical upper tail heat stress conditions for at least an additional 10 consecutive hours. Further, we explicate a significant rise in the number of people exposed to severe humid heat stress, concentrated along most coastal regions, by as much as 1.6 million person‐days. More than 20% of Turkey's population may confront severe humid heat stress for at least 1 h, with that percentage falling to 4.15% for at least five consecutive hours, which indicates that people will not only endure more intense humid heat stress but also be exposed to these conditions consecutively over a period of many hours.

Acoustic parameters and wave velocity field evolution characteristics in the high-stress region of granite

The risks of underground engineering originate from disturbances caused by the alterations to the stress field within rock masses. To understand the variations in the stress field within rock masses under local stress conditions, a set of double-hole pullout tests was used to design. The MS/AE location and tomography were used to investigate the evolution characteristics of the internal stress concentration area in granite. Results show that acoustic emission events transition from scattered and nucleated to regionally aggregated, then developing into arc-shaped patterns, ultimately exhibiting a trend of semi-spherical distribution. The increase in the fractal dimension of the spatial distribution of acoustic emission events indicates increased stress concentration degree, whereas a decrease signifies dispersion or transfer. Tomography results demonstrate the dynamic development process of local stress zones, with initial velocity anomalies as the core, continuously aggregating and expanding, showing an increasing difference trend compared to the velocity outside the region.

Impact of stress regime change on the permeability of a naturally fractured carbonate buildup (Latemar, the Dolomites, northern Italy)

Abstract. Changing stress regimes control fracture network geometry and influence porosity and permeability in carbonate reservoirs. Using outcrop data analysis and a displacement-based linear elastic finite-element method, we investigate the impact of stress regime change on fracture network permeability. The model is based on fracture networks, specifically fracture substructures. The Latemar, predominantly affected by subsidence deformation and Alpine compression, is taken as an outcrop analogue for an isolated (Mesozoic) carbonate buildup with fracture-dominated permeability. We apply a novel strategy involving two compressive boundary loading conditions constrained by the study area's NW–SE and N–S stress directions. Stress-dependent heterogeneous apertures and effective permeability were computed in the 2D domain by (i) using the local stress state within the fracture substructure and (ii) running a single-phase flow analysis considering the fracture apertures in each fracture substructure. Our results show that the impact of the modelled far-field stresses at (i) subsidence deformation from the NW–SE and (ii) Alpine deformation from N–S increased the overall fracture aperture and permeability. In each case, increasing permeability is associated with open fractures parallel to the orientation of the loading stages and with fracture densities. The anisotropy of permeability is increased by the density and connectedness of the fracture network and affected by shear dilation. The two far-field stresses simultaneously acting within the selected fracture substructure at a different magnitude and orientation do not necessarily cancel each other out in the mechanical deformation modelling. These stresses affect the overall aperture and permeability distributions and the flow patterns. These effects – potentially ignored in simpler stress-dependent permeability – can result in significant inaccuracies in permeability estimation.

Accurate determination of active slip systems using a novel lattice rotation analysis: Quasi-in-situ EBSD/ECCI study on the homogeneous/heterogeneous deformation of polycrystalline zirconium

The accurate determination of slip systems is crucial for understanding deformation mechanisms of metals in particular, how slip/twining transfer or blocking occurs. This paper proposes a novel crystallographic lattice rotation analysis, which combines quasi-in-situ Electron Back Scatter Diffraction (EBSD) and Electron channeling contrast imaging (ECCI) techniques to precisely predict the active slip system with specific slip directions in polycrystalline Zr. Large data sets for hundreds of grains were quantitatively analyzed by this method to statistically gain deep insight on the deformation mechanism of Zr. Our results show that prismatic 〈a〉 slip is the primary slip system in most grains, and pyramidal 〈a〉 slip plays a secondary role in coordinating the local stress state of plastic deformation. The predicted slip direction is the same as the direction in which dislocation density increases. This method can serve as a complementary tool to the intragranular misorientation axis (IGMA) method, bridging the gap between macro-mechanical response and microstructural deformation mechanisms.

Evidence of Dislocation Mixed Climb in Quartz From the Main Central and Moine Thrusts: An Electron Tomography Study

AbstractIn this study we apply electron tomography to characterize 3D dislocation microstructures in two quartz mylonite specimens from the Moine and Main Central Thrusts, both of which accommodated displacements by dislocation creep in the presence of water. Both specimens show dislocation activity with dislocation densities of the order of 3–4 × 1012 m−2 and evidence of recovery from the presence of subgrain boundaries. 〈a〉 slip occurs predominantly on pyramidal and prismatic planes (〈a〉 basal glide is not active). [c] Glide is not significant. On the other hand, we observe a very high level of activation of 〈c + a〉 glide on the , , (n = 1,2) and even planes. Approximately 60% of all dislocations show evidence of climb with a predominance of mixed climb, a deformation mechanism characterized by dislocations moving in a plane intermediate between the glide and the climb planes. This atypical mode of deformation demonstrates comparable glide and climb efficiency under natural deformation conditions. It promotes dislocation glide in planes not expected for the quartz structure, probably by inhibiting lattice friction. Our quantitative characterization of the microstructure enables us to assess the strain that dislocations can generate. We show that glide systems indicated by the observed dislocations are sufficient to accommodate any arbitrary 3D strain by themselves. Although historically dislocation glide has been regarded as being primarily responsible for producing strain, activation of climb can also directly contribute to the finite strain. On the basis of this characterization, we propose a numerical modeling approach for attempting to characterize the local stress state that gave rise to the observed microstructure.

Dynamic mechanical response and constitutive model of (Ti37.31Zr22.75Be26.39Al4.55Cu9)94Co6 high-entropy bulk metallic glass

In this work, the mechanical response and fracture characteristics of (Ti37.31Zr22.75Be26.39Al4.55Cu9)94Co6 high-entropy bulk metallic glass (HE-BMG) were investigated in detail over a wide range of strain rates (10−4–105 s−1). The HE-BMG exhibited a negative strain rate sensitivity under uniaxial compression, with the strength showing more significant rate dependence under dynamic conditions. The shear band behavior translated from the dominance of multiple shear bands propagations under quasi-static compression to the rapid propagation of a single shear band to form cracks under dynamic compression. Under dynamic loading, the shearing velocity increased along an arc-shaped displacement path, and the local stress state on the shear fracture surface shifted from compressive-shear to tensile-shear, accompanied by changes in fracture morphologies. The spall strength of HE-BMG decreased as flyer impact speed increased, while the long-term dependence of spall strength on strain rate may be positive. With the increase of impact speed, the main microstructural features of the spalling surface translated from flat regions accompanied by dimples to cup-cone structures. Furthermore, the complete parameters of the Johnson–Holmquist II (JH-2) model for HE-BMG were obtained based on experimental data. Numerical simulations of planar impact, penetration, and hypervelocity impact are in good agreement with experimental results, demonstrating the validity of the JH-2 model parameters. The current work has important guiding value for the application of HE-BMG in space debris protection.

Quantifying the expansion rates of aftershock zones for magnitude-7 class earthquakes around the Japanese archipelago

Earthquakes (mainshocks) trigger sequences of aftershocks, the frequency of which diminishes following a power-law decay, while the spatial domain of these aftershocks extends logarithmically over time. The delineation of the aftershock zone can be modulated by variables beyond the magnitude of the mainshock, encompassing the location of the fault (whether the fault is at a plate boundary), the depth at which the event occurs, and the prevailing local stress conditions. Here, we evaluate the expansion rate of aftershock zones by analyzing earthquakes of magnitude-7 class in the vicinity of the Japanese archipelago. Prior studies have offered approximate assessments of expansion rates; however, our approach involves the utilization of a straightforward algorithm for the automated estimation of this metric, facilitating the compilation of a catalog. Across the dataset, no pronounced correlations were discerned between the expansion rate and other examined parameters. Yet, an inverse relationship is identified between the expansion rate and the b value of aftershocks for mainshocks occurring at plate boundaries. This observation suggests that the expansion rate of aftershock zones predominantly mirrors the stress field following the mainshock. Such a pattern is not detected in mainshocks occurring within the plate's interior. While the expansion rate of aftershock zones is likely influenced by various factors, aftershock zones may expand more rapidly with higher differential stress in areas surrounding hypocenters of major interplate earthquakes of magnitude 8 or 9.

On Defect Evolution in EBM Additively Manufactured Ti-6Al-4V via In Situ Investigations.

This study concerned the in situ investigation of the defect evolution and fracture mechanism of additively manufactured (AM) Ti-6Al-4V under uniaxial tensile tests. In order to achieve this, microstructure characterization was initially carried out in order to identify the defects within the matrix of the candidate material. In situ testing was then performed, focusing on the spherical defect to observe its evolution under tensile loading. It was found that, before the fracture stage, the geometric evolution of the spherical defect towards an ellipse shape was dominated by the load in the tensile direction. In addition, the slip band density was found to be aggravated near the spherical defect due to the geometric discontinuity-induced stress concentration. During the fracture process, the defect geometry evolved as an irregular shape, which was mainly attributed to the micro-void-induced localized multi-axial stress state. The fracture analysis indicated that defects play a key role in crack initiation, leading to the fracture of LPBF materials.

Stress-constrained optimization of multiscale structures with parameterized microarchitectures using machine learning

A multiscale topology optimization framework for stress-constrained design is presented. Spatially varying microstructures are distributed in the macroscale where their material properties are estimated using a neural network surrogate model for homogenized constitutive relations. Meanwhile, the local stress state of each microstructure is evaluated with another neural network trained to emulate second-order homogenization. This combination of two surrogate models — one for effective properties, one for local stress evaluation — is shown to accurately and efficiently predict relevant stress values in structures with spatially varying microstructures. An augmented lagrangian approach to stress-constrained optimization is then implemented to minimize the volume of multiscale structures subjected to stress constraints in each microstructure. Several examples show that the approach can produce designs with varied microarchitectures that respect local stress constraints. As expected, the distributed microstructures cannot surpass density-based topology optimization designs in canonical volume minimization problems. Despite this, the stress-constrained design of hierarchical structures remains an important component in the development of multiphysics and multifunctional design. This work presents an effective approach to multiscale optimization where a machine learning approach to local analysis has increased the information exchange between micro- and macroscales.

Constitutive modelling of glassy polymers considering shear plasticity and craze yielding

The inelastic behaviour of thermoplastic polymers below the glass transition temperature can involve shear plasticity or crazing, depending on the strain rate and temperature. Shear plasticity is driven by local shear stresses and is essentially volume preserving. In contrast, crazing is a failure phenomenon occurring under tension and causes significant volume change. In some cases, crazing can also result in large inelastic deformations at macroscopic scale, referred to as craze yielding. The aim of this study is to propose a thermodynamically-consistent constitutive framework that accounts for both shear plasticity and craze yielding in glassy polymers. The main assumption is that shear plasticity and craze yielding occur exclusively from each other, depending on the local stress state. Both are modelled as thermally-activated processes with different activation stresses and flow rules. The theory is validated against experimental data for polylactic acid (PLA). Our model is capable of reproducing the stress–strain response including yielding, softening, and drawing under both shear plasticity and craze yielding, and also accurately predicts the volumetric deformation under tension. In particular, our simulations well capture the interesting phenomenon that craze yielding stabilises localised deformations and prevents necking. Our model can also simulate complex scenarios where both mechanisms occur simultaneously at different locations of the specimen.

The spatial and temporal characterization of type IV creep fracture in notched specimen of modified 9Cr–1Mo steel

In order to observe type IV fracture occurring in the HAZ fine-grained zone of a modified 9Cr–1Mo steel at high resolution, the creep damage process of a single notched specimen was continuously observed in three dimensions. The initiation, growth and coalescence behaviour of individual creep voids were clearly observed. This was combined with an analysis of the local stress state using the finite element method. Furthermore, all creep voids that could be observed in the final loading stage were traced backward in time to determine the loading stages at which the creep voids were initiated together with the growth behaviour of the individual creep voids. Significant creep void formation was observed a little further in from the bottom of the notch, where the maximum principal and hydrostatic stresses were greatest, while the stress triaxiality was greatest a little closer to the centre of the specimen than the maximum principal and hydrostatic stresses, where significant creep void growth was observed. It was found that the creep voids nucleated at an early loading stage dominated the creep failure of the specimen. The growth behaviour could be well approximated by the model that combined the diffusion law and power creep law with the constraint effect from neighbouring grains. Engineering parameters that describe the progress of type IV damage were also investigated.

The mechanistic origins of heterogeneous void growth during ductile failure

In 1969, Rice and Tracey proposed that the rate of void growth during ductile failure is a strong function of the hydrostatic stress state. Numerous in-situ x-ray computed tomography (XCT) studies demonstrated that the average rate of void growth is well-predicted by modified versions of the Rice-Tracey equation. However, recent in-situ XCT studies of void growth demonstrated that individual voids grow in highly heterogeneous manners in a way that is not predicted by Rice-Tracey or similar models. Model-based studies using crystal plasticity finite element (CP-FE) suggest that local effects of grain orientation play a strong role during void growth, but we lack experimental data to test this hypothesis. The present study leverages recent advances in laboratory-based diffraction contrast tomography (Lab-DCT) and in-situ XCT to systematically examine the effects of grain orientation on void growth experimentally in an Al-2219 alloy. CP-FE modeling was used to assess the local stress states associated with voids and their influence on void growth rates. These data indicate that void growth is not simply controlled by the local hydrostatic stress or stress triaxiality. Additionally, no clear relationship between void growth rates and grain orientation were observed. Instead, void growth rates were highly stochastic in ways that could not be directly linked to the local stress or strain states. The combination of experimental and modeling data suggests that our current assumptions regarding the mechanisms of void growth are flawed and that additional factors, which may include the local dislocation density, affect the rate of void growth.

Effect of crustal stress state on magmatic stalling and ascent: case study from Puyehue-Cordón Caulle, Chile

The Southern Volcanic Zone (SVZ) in Chile is an active continental arc with a complex history of volcanism, where a range of magmatic compositions have been erupted in a variety of styles. In the Central SVZ, both monogenetic and polygenetic volcanoes exist, in close proximity to the Liquiñe-Ofqui Fault System (LOFS), but with variable local stress states. Previous studies have inferred varying crustal storage timescales, controlled by the orientation of volcanic centres relative to the N-S striking LOFS and σHMax in this region. To assess the relationship between volcanism and crustal stress states affected by large-scale tectonic structures and edifice controls, we present whole rock geochemical data, to ensure consistency in source dynamics and crustal processing, mineral-specific compositional data, thermobarometry, and Fe–Mg diffusion modelling in olivine crystals from mafic lavas, to assess ascent timescales, from the stratovolcanic edifice of Puyehue-Cordón Caulle and proximal small eruptive centres. Textural observations highlight differences in crystal maturation timescales between centres in inferred compression, transpression, and extension, yet source melting dynamics remain constant. Only samples from the stratovolcanic edifice (in regional compression) preserve extensive zonation in olivine macrocrysts; these textures are generally absent from proximal small eruptive centres in transtension or extension. The zonation in olivines from stratovolcanic lavas yields timescales on the order of a few days to a few weeks, suggesting that even in environments which inhibit ascent, timescales between unrest and eruption of mafic magmas may be short. Significantly, high-resolution compositional profiles from olivine grains in the studied samples record evidence for post-eruptive growth and diffusion, highlighting the importance of careful interpretation of diffusion timescales from zoned minerals in more slowly cooled lavas when compared with tephra samples.