Stress Field Research Articles

IGA Approximations of Elastic Interfaces and their Defects in an Elastic Medium with Couple Stress

AbstractThe objective of this work is to develop and implement a computational algorithm for calculating stress and couple-stress fields induced by bulk and interfacial line defects such as dislocations and generalized disclinations within phase/grain boundaries. The thermodynamic driving forces on line and planar defects, responsible for coupled plasticity and interface motion, fully coupled to stress, couple-stress and applied boundary conditions are also computed. A continuum approach for small deformations is considered following the formulation outlined in Acharya and Fressengeas (Continuum mechanics of the interaction of phase boundaries and dislocations in solids. In: Differential Geometry and Continuum Mechanics, pages 123–165, 2015), extended herein in the thermodynamics to accommodate physically necessary ingredients that arose in the modeling in Zhang and Acharya (J Mech Phys Solids 119:188–223, 2018), Zhang et al. (J Mech Phys Solids 114:258–302, 2018). Constitutive relations are derived from kinematics, balance laws and from the use of the second law of thermodynamics in global form. One of the challenges presented by this approach is the inclusion of couple stresses, Toupin (Arch Rational Mech Anal 17(2):85–112, 1964), and the consequent treatment of 4th order systems arising from the equations of balances of linear and angular momentum. In order to deal with these equations, the classical FEM approach is replaced by iso-geometric analysis (IGA), as proposed by Hughes et al. (Comput Methods Appl Mech Eng 194(39):4135–4195, 2005). Results on stress and couple stress fields of the various defects involved are computed. Further, the coupling of the dislocation density with the eigenwall allows for the capturing of the shear parallel to the grain boundary, which has been observed experimentally and through molecular dynamics.

Stress state and strength of welded joints with V-shaped elements

It has been established that in the vicinity of V-shaped elements with a sharp apex, the stress field has a singularity of type σ→∞, which requires the use of the concept of “stress intensity factor” in strength calculations. This need arose due to different dimensions of the SIF and fracture toughness. Special criteria for the destruction of structures with V-shaped stress concentrators have been developed, which include, as special cases, the criteria of classical mechanics and the criteria of fracture mechanics. An experimental test of the proposed criteria and calculation methods was carried out. The information obtained is useful for analyzing welded joints on main pipelines and other equipment.

Analytical study on plastic loosening zone of weak surrounding rocks tunnel in high cold areas

AbstractThe analysis of plastic loosening zone of tunnels with cold condition, high in situ stress, and weak stratum is the key scientific problem of tunnel construction in cold areas. Based on the unified strength theory and considering the coupling effect of high in situ stress and low‐temperature field, the analytical model of the plastic loosening zone of tunnel surrounding rocks was proposed, and the influence laws of physical parameters (initial ground stress, freezing temperature, excavation radius, and intermediate principal stress) was explored. The results illustrate that: (1) with the coupling of low temperature and stress fields, the in‐situ stress has a more significant effect on the radius of the plastic loosening zone of the tunnel surrounding rocks, showing a nonlinear correlation. When the in situ stress is small, the influence of the temperature field on the radius of the plastic loosening zone is more obvious. (2) The radius of the plastic loosening zone of the tunnel surrounding rocks increases synchronously with the excavation radius and gradually decreases with the increase of the intermediate principal stress. (3) The greater the temperature difference in the low‐temperature field, the tunnel stability has a more significant response to the intermediate principal stress and is positively correlated with the coefficient of the intermediate principal stress. (4) In the actual project, the proposed model can well describe the mechanical behavior and deformation characteristics of tunnel surrounding rocks in low‐temperature environment, showing applicability. The research results are expected to provide a theoretical basis for the study of tunnel stability in cold areas.

Glacial Hydro-Isostatic Adjustment at Mid-Ocean Ridges

Summary Recent studies have suggested a link between ice age sea level fluctuations and variations in magma production and crustal faulting along mid-ocean ridges based on the detection of Milankovitch cycle frequencies in topography off several ridges. These fluctuations have also been connected to variability in hydrothermal metal fluxes near ridges. Ice age sea level calculations have shown that the sea level change across glacial cycles will be characterized by significant geographic variability, that is, departures from eustasy, due to the gravitational, deformation and rotational effects of the glacial isostatic adjustment (GIA) process. Using a state-of-the-art GIA simulation that incorporates 3-D variations in Earth viscoelastic structure, including plate boundaries, and updated constraints on the magnitude and geometry of ice mass fluctuations, we predict global sea level changes from Last Glacial Maximum (LGM, 26 ka) to present and from the Penultimate Glacial Maximum (PGM, 143 ka) to the Last Interglacial (LIG, 128 ka). We focus on the results along three ridges: the Mid-Atlantic Ridge, Juan de Fuca Ridge and East Pacific Rise, which are examples of slow, intermediate and fast spreading ridges, respectively. Sea level change across the Mid-Atlantic Ridge shows the greatest variability, ranging from a sea level fall greater than 200 m in Iceland to a maximum rise of ∼150 m in the South Atlantic, with significant non-monotonicity north of the Equator as the ridge weaves across the field of sea level changes. We also calculate changes in crustal normal stress from LGM to present-day across the Mid-Atlantic and Juan de Fuca Ridges and the East Pacific Rise. These results indicate that the contribution from ice mass changes to the crustal stress field can be significant well away from the location of ancient ice complexes. We conclude that any exploration of the hypothesized links to magma production and crustal faulting must consider both ocean and ice loading effects and, more generally, the profound geographic variability of the GIA process.

A two‐scale numerical analysis of intra‐yarn hybrid natural/synthetic woven composites

AbstractThis paper investigates how combining natural and synthetic fibers within yarns, altering the degree of hybridisation, and varying weave architectures affect the homogenized elastic properties and stress distributions in 2D woven composite laminae using a two‐scale homogenization scheme. The novelty lies in integrating a micro‐scale representative volume element model with a meso‐scale repeating unit cell model to capture intra‐yarn natural/synthetic fiber hybridisation. The study focuses on 2/2 twill and 5‐harness satin woven laminae with flax/E‐glass, hemp/E‐glass and basalt/E‐glass hybrid yarns, all with a total fiber volume fraction of 0.6. Fiber distributions are created through a random sequential expansion algorithm for hybrid RVEs, while periodic meso‐structures are used to define weave architectures for RUCs. Results show that yarn‐level fiber hybridisation significantly influences matrix stress distributions, with variations of up to 22%. In contrast, weave architecture has a minimal effect, with variations in homogenized properties below 2%. Varying fiber types, degrees of hybridisation and weave architectures allow for the tailoring of yarn and lamina properties and altering the stress distributions at micro‐ and meso‐scales and potentially influencing damage mechanisms. The modeling approach for analyzing the mechanical behavior of intra‐yarn hybrid natural/synthetic woven laminae could potentially be used to tailor their tolerance to damage.Highlights Two‐scale homogenization integrates RVE and RUC for 2D woven hybrid composites. Analyzed flax/E‐glass, hemp/E‐glass, basalt/E‐glass in twill and satin weaves. Intra‐yarn hybridisation significantly affects properties and stress fields. Hybridizing fibers offers tailored properties and may improve damage tolerance.

Experimental and Numerical Indentation Tests to Study the Crack Growth Properties Caused by Various Indenters Under High Temperature

ABSTRACTThis paper investigates the influence of indenter shape and rock texture on the crack growth properties and rock hardness. For this purpose, two different rock specimens such as basalt and marble with various textures were prepared and tested by Vickers indentation hardness device with rhombus indenter shape under two different temperatures of 35°C and 100°C. Concurrent with the experimental test, Vickers indentation simulations have been done on three calibrated rock models with nine different indenter shapes. The tensile strength of marble was 8 MPa, while basalt had a tensile strength of 10.8 MPa. Regarding compressive strength, marble exhibited 71 MPa, whereas basalt had a compressive strength of 153 MPa. Marble and basalt had elastic moduli of 45 and 95 GPa, respectively. The physical loading was applied vertically at a rate of 0.004 mm/min. The rock's texture and temperature significantly impact Vickers indentation hardness. Penetration by the Vickers indenter creates an elastic‐plastic stress field, leading to radial cracks due to exceeding the critical stress level. Ring cracks are caused by bulging of the test material and nested cracks directly under the indentation due to high shear and bending stresses in the region. The indentation shape affected the extent of the damage zone below it, leading to increased crack growth with smaller indentation diameters. The radial fracture number increased when the cross‐section changed from circular to quadrilateral shape. The crack growth and damage zone area increased with higher rock temperature due to increased rock brittleness.

An Analytical Approach for the Near‐Tip Field Around V‐Notch in Orthotropic Materials

ABSTRACTThis study introduces a novel theoretical framework for determining the notch stress intensity factors (NSIFs) in composite materials. The approach involves utilizing Irwin's complex stress functions for Mode I and Mode II loading conditions. By investigating displacement and singular stress components around V‐notch tips, the research aims to provide a comprehensive understanding of NSIFs in composite materials. Three different composite materials were considered in the study to explore the relationship between stress and material properties. The stress singularity order, denoted as , is determined through eigenequations. Moreover, explicit solutions for displacement and singular stress fields are derived for a more thorough understanding of the behavior around notches in composite materials. To validate the proposed theoretical framework, a rigorous numerical study is conducted using FEM on finite‐size notches. The results obtained from the numerical simulations are compared with theoretical predictions to assess the accuracy and reliability of the developed approach.

Constraint on the background stress in the source region of the 2016 Kumamoto earthquake sequence based on temporal changes in elastic strain energies and coseismic stress rotation

Summary We investigated the deviatoric stress magnitude of the background stress fields before the 2016 Kumamoto earthquake sequence in Japan based on temporal changes in elastic strain energies caused by the mainshock and coseismic stress rotation. We modelled the six components of background stress fields from stress orientations together with the parameter of effective friction coefficients (μ’ = 0.3, 0.15, 0.1, 0.05, and 0.03) using the 3-D Mohr diagram. We computed the absolute stress fields immediately following the primary events (the largest foreshock and mainshock) of the earthquake sequence by combining coseismic stress change fields with background stress fields. The total amount of elastic strain energy released by the mainshock increased with the effective friction coefficient. Considering the energy balance in which some part of the released energy must be consumed as radiated energy, we understood that the model with μ’ = 0.03 was unrealistic. We also examined the dependence of coseismic stress rotation on the effective friction coefficient. Furthermore, we applied a stress inversion method to moment tensor data of earthquakes to directly estimate coseismic stress rotation. Comparing the theoretical and estimated coseismic stress rotations, we concluded that the models with μ’ = 0.3 and 0.15 were more consistent with the observations than those with μ’ < 0.1. In the reasonable models, the deviatoric stress magnitudes were 37–65 and 39−70 MPa at a depth of 10 km on the northern and southern source faults, respectively.

Using Resistivity Structure to Study the Seismogenic Mechanism of the 2021 Luxian Ms6.0 Earthquakes

Over the past few years, there has been a noticeable change in the occurrence of seismic disasters in Sichuan, China. The focus has shifted from Western Sichuan to the previously more stable Southeastern Sichuan. The recent Ms6.0 earthquake in Luxian, Southeastern Sichuan, on 16 September 2021, has once again captured the interest of scholars, who are closely examining the seismogenic environment and potential seismic hazards in the region. We conducted a magnetotelluric (MT) array survey in the Luxian earthquake area to explore the deep seismogenic environment of the 2021 Luxian Ms6.0 earthquake zone and understand the potential effects of industrial extraction on seismic activities. Here are the insights we obtained: Underneath the anticline in the Luxian Ms6.0 earthquake area, there is a structure that mainly exhibits high resistance. On the other hand, beneath the syncline, a structure with medium to low resistance is observed. The epicenter of the mainshock was identified near the intersection of high- and low-resistance media within the Fuji syncline area. Smaller aftershocks that followed the mainshock were mainly concentrated in the low-resistance layers at depths of 3–5 km in the Fuji syncline area. MT survey results have confirmed the existence of a detachment zone in the shallow crust near the epicenter of the Luxian Ms6.0 earthquake. It is believed that this detachment layer played a significant role in the seismogenic process of the Luxian Ms6.0 earthquake. During different stress conditions, this layer became active and caused the compression and faulting of a hidden fault below, resulting in the Luxian Ms6.0 earthquake. After the main earthquake, a series of smaller aftershocks with varying focal mechanisms occurred as the stress fields continued to release. It is important to note that the Luxian Ms6.0 earthquake highlights the ongoing high stress levels in the southern region of the Sichuan Basin. This emphasizes the need for continued monitoring and consideration of potential seismic hazards in the southern Sichuan area.

A thermodynamic based constitutive model considering the mutual influence of multiple physical fields

In multiple physical fields, the mutual influence among these fields can significantly impact material elastoplasticity. This paper proposes a thermodynamic-based constitutive model that incorporates the mutual influence of multiple physical fields. Rather than treating physical field characteristics as adjustable “parameters” affecting material coefficients, the proposed model employs a thermodynamic dissipation potential derived from the Onsager reciprocity relations, accounting for thermodynamic forces coupling. This dissipation potential ensures that the thermodynamic flow in the stress field is influenced by both stress field and other physical fields thermodynamic forces, which describes the plastic flow under multiple physical fields, while preserving thermodynamic duality. The paper begins with the formulation of a generalized thermodynamic model applicable to diverse materials and types of coupled fields, which is then degraded to a specific model for AA5182-O AlMg alloy under the influence of temperature and strain rate fields coupling. Given the universal applicability of the generalized model, such degradation provides a structured approach framework for developing thermodynamics-based constitutive models. For different materials encountered in practical engineering, new thermodynamic forces can be introduced to describe their unique mechanical properties while preserving the overarching thermodynamics-based model framework, thereby facilitating model scalability. The paper concludes with a validation example, showing that within the Portevin-Le Chatelie (PLC) regime, the plastic flow stress of AA5182-O AlMg alloy decreases with increasing strain rate at low temperatures but increases at high temperatures. The accurate simulation of these distinct strain rate effects crucially relies on integrating the mutual influence of temperature field and strain rate field.

An enhanced numerical model for considering coupled strain-softening and seepage effects on rock masses surrounding a submarine tunnel

The seepage of groundwater and the strain-softening of rock mass in a submarine tunnel expand the plastic region of rock, thereby affecting its overall stability. It is therefore essential to study the stress and strain fields in the rocks surrounding the submarine tunnel by considering the coupled effect of strain-softening and seepage. However, the evolution equation for the hydro-mechanical parameters in the existing fully coupled solution is a uniform equation that is unable to reproduce the characteristics of rock mass in practice. In this study, an updated numerical procedure for the submarine tunnel is derived by coupling strain-softening and seepage effect based on the experimental results. According to the hydro-mechanical coupling theory, the hydro-mechanical parameters such as elastic modulus, Poisson's ratio, Biot's coefficient and permeability coefficient of rocks are characterized by the fitting equations derived from the experimental data. Then, the updated numerical procedure is deduced with the governing equations, boundary conditions, seepage equations and fitting equations. The updated numerical procedure is verified accurately compared with the previous analytical solution. By utilizing the updated numerical procedure, the characteristics of stress field and the influences of initial pore water pressure, Biot's coefficient, and permeability coefficient on the stress, displacement and water-inflow of the surrounding rocks are discussed. Regardless of the variations in hydro-mechanical parameters, the stress distribution has a similar trend. The initial permeability coefficient exerts the most significant influence on the stress field. With the increases in initial pore water pressure and Biot's coefficient, the plastic region expands, and the water-inflow and displacement increase accordingly. Given the fact that the stability of the tunnel is more sensitive to the seepage force controlled by the hydraulic parameters, it is suggested to dewater the ground above the submarine tunnel to control the initial pore water pressure.

Characteristics of rock strength failure and identification of hazardous areas in the roof of deep mining stope

The roof of deep mining stopes is notoriously unstable and challenging to manage due to high in-situ stress and deteriorating rock conditions. To address this issue, a comprehensive approach was employed, incorporating indoor physical experiments, numerical simulations, and field tests to investigate the stope roof rock composition characteristics and strength failure modes in the Laoyingshan deep coal mine. The study revealed the distribution characteristics of the high in-situ stress field in this mining area and delineated different levels of hazards in the roof. The results indicate that, under identical lithological conditions, the mineral composition of deep rock formations is complex. Changes in some mineral components and properties were observed, leading to increased expansiveness and water absorption in clay minerals such as sodium montmorillonite and kaolinite within the basic roof layer. In simulated tests of surrounding rock strength in deep mining stopes, as stress levels approached 2 MPa, the damage factor increased, showing a trend of localized concentration. Internal damage began to exhibit noticeable spatial evolutionary expansion patterns and gradual nucleation. When stress exceeded 8 MPa, rocks transitioned from macroscopic fracture surfaces to microscopic fragmentation zones, with a sudden decrease in stiffness degradation values, internal instantaneous collapse, and phenomena akin to “rock burst” under high stress disturbance. In-situ stress testing determined that the stress field in the deep mining zone of Laoyingshan is in a thrust fault stress state, with horizontal stress predominance. The maximum stress direction is predominantly southwest-northeast, with stress magnitude reaching 25.49 MPa. By arranging post-mining and pre-mining borehole observation stations and utilizing image grayscale fracture recognition technology in MATLAB, the roof was classified into four danger levels: “extremely broken,” “broken,” “significant fractures,” and “minor fractures.” Main control measures for coupling grouting in deep-shallow zonation rigidity were proposed based on these classifications.

Workflow Helps Predict Casing Deformation During Hydraulic Fracturing in Shale Gas

_ This article, written by JPT Technology Editor Chris Carpenter, contains highlights of paper URTeC 4043054, “Casing-Deformation-Risk Prediction With Numerical Simulation During Hydraulic Fracturing in Deep Shale Gas Reservoirs,” by Jianfa Wu, Xuewen Shi, and Qiyong Gou, PetroChina, et al. The paper has not been peer reviewed. _ In this study, a numerical simulation workflow was established based on comprehensive modeling of a time-lapse stress field, hydraulic fracturing simulation, natural fractures, and other weak interfaces to analyze casing-deformation mechanisms of shale gas horizontal wells during the multistage hydraulic fracturing process. The authors write that they aim to guide the prediction of, and provide mitigating solutions for, casing-deformation risks while improving stimulation efficiency. Introduction The burial depth of deep shale gas reservoirs in the southern Sichuan Basin is approximately 3500 to 4500 m. These have complex geological conditions with developed faults and fractures and high in-situ horizontal stresses. During multistage hydraulic fracturing operations in horizontal wells, establishing an effective prediction method for casing-deformation risks in these reservoirs is a challenge. An integrated workflow is proposed to evaluate and predict casing-deformation risks, effectively support hydraulic fracturing operations, and reduce and mitigate casing failure in these reservoirs. Modeling Methodology It is essential to consider the effects of multiscale natural fractures when constructing a 3D stress model. Natural-fracture stability can be evaluated and used to identify unstable fractures that may cause casing failures. The fundamental workflow of predicting and preventing the casing failure is as follows: 1. Build a 3D geological model characterizing the spatial distribution and geometry of faults and natural fractures 2. Construct a 3D stress model to characterize the in-situ stress state 3. Perform numerical simulations of hydraulic fracturing and 4D coupled stress simulation to obtain the stress field during and after hydraulic fracturing 4. Quantitatively analyze the shear slip risks of natural fractures by calculating slip tolerance as an indicator to determine risk levels of casing deformation 5. Optimize and adjust engineering parameters for hydraulic fracturing to reduce and mitigate casing deformation Natural-Fracture Distribution and Model Construction. In this study, natural fractures parameters such as density, strike, inclination, length of extension, and height were obtained. The natural fractures were calibrated and verified with cores and wellbore images. Natural-fracture distribution was then validated with the expected tectonic history and present-day stress regime. In Fig. 1, the left-hand plot depicts the ant-tracking results by seismic attribute, the central plot portrays natural fractures validated by microseismic events, and the right-hand plot illustrates the 3D natural-fracture model. Results showed two major fracture strikes, one in the northeast and the other in the northwest, forming a mesh network distribution that will be used in the stress model, fracturing simulation, and slip analysis.

Southeastern extension of Singhbhum Shear Zone, Eastern Indian Shield: A detailed appraisal

Trans-crustal, northerly dipping, arcuate Singhbhum Shear Zone (SSZ) constituting the interface of Proterozoic North Singhbhum Mobile Belt (NSMB) and Singhbhum Craton (SC) is established over a strike length of about ∼200 km from Baharagorato Chakradharpur in the west beyond which it bifurcates around the wedge shaped Chakradharpur granite. Based on geological mapping and meso to micro scale structural analysis, present work examines the SSZ’s southeastern extension beyond Baharagora, indicating its southward continuation for about 50 km to Dubukidihi within the NSMB lithopackage. Bangriposi Shear Zone (BSZ) which defines the terrane boundary between the SC and NSMB in this sector, splays towards north as Bangriposi West Shear Zone (BWSZ) and both transects the Dhanjori Basin longitudinally. BSZ and BWSZ are interpreted to be cogenetic with SSZ. NSMB in this area comprises of pelitic and psamopelitic rocks intruded by mafic–ultramafic magmatism and granitoids. Considering a north to south compression regime during accretion of NSMB with SC, the study area, which forms the interface between the craton and mobile belt, is at an acute to obtuse angle with the direction of overall stress field. It is marked by several high strain zones characterized by N-S shear band cleavages with moderate easterly dips, and northerly pitching stretching/mineral lineation. Consistent westerly vergent folds on mylonitic plane and orientation of stretching lineation indicates top to WSW movement in the shear zones. Tectonically interleaved mafic–ultramafic schist within the Romapahari granite are sites of copper mineralization with enhanced nickel concentration pointing towards future exploration sites for epigenetic copper, uranium mineralization which is so characteristic of SSZ.

Lithological control and structural inheritance on faults growth in multilayer foreland sequences

Foreland sectors and foredeep-forebulge systems are affected, as the orogenic wedge migrates, by successive stages of stress states and tectonic deformation, resulting in the development of complex fault networks, even if characterized by limited deformation. The role played by structural inheritance and changes in stress field through time, in influencing the successive re-activations of fault segments, is still a topic to be thoroughly investigated. In this work, thanks to an extensive database made available by courtesy of Energean, we were able to investigate a foreland sector at the margin of the southern Apennines. By means of thickness analysis of the Neogene foredeep sequence and of displacement analysis along the fault network, we documented a shift from forebulge-related extension, in Early Pliocene, to a new tectonic phase, since Middle Pliocene, related to a strike slip stress field, possibly related to the activity of the Tremiti Fault Zone. We also characterized the geometry and connectivity of the cover-restricted faults, developing above propagating normal faults and observed a clear correlation between fault propagation tendency and lithological/mechanical layering within the cover units.

Tectonomagmatic evolution of Pune – Nasik Deccan Dykes: Insights from structure and dimension scaling

The Deccan Continental Flood Basalt of the Indian Peninsula is characterized by extensive basaltic eruptions ornamented with three spectacular distinct dyke swarms: the Pune – Nasik, Narmada – Tapi, and Western Coastal dyke swarms. Our study area is the Pune – Nasik dyke swarm, which has ~465 mappable dykes. These dykes exhibit different orientations with a predominant trend of N101° and vary in length from less than 1 km to ~64 km. These dykes are massively jointed and occasionally contain vesicles filled with secondary minerals like quartz and calcite. The host rock is weathered basalt of various older Deccan flows. In this study, we have calculated magmatic overpressures and magma chamber depths using the aspect ratios (length/thickness) of the dykes. The average estimated source depth is ~13 km, based on an average Young's modulus for the host rock basalt (Eavg, 7.5 GPa). Additionally, we compared the inferred magma source depths of the Pune – Nasik, and Narmada-Tapi dyke swarms which include the Nandurbar – Dhule, and Pachmarhi dykes of the Deccan Flood Basalt Province. Our findings indicate that the magma chamber source depth is greater in the Pune – Nasik dyke swarm compared with other dyke swarms. The variation in strike distribution of the Pune-Nasik dyke swarm may be attributed to several factors, including a larger magma chamber, local stress fields generated by shallow magma chamber, or the superimposition of tectonic stress fields (N-S and E-W extension) during the emplacement of dykes. This contrasts with the commonly held belief that the dykes are solely connected to a central edifice of the Reunion hotspot.

Critical length prediction for fiber breaking in injection molded of long-fiber reinforced thin-walled products

Fiber length is one of the vital factors in determining the properties of long fiber reinforced thermoplastics (LFRT). However, many factors may generate fiber breakage in the injection molding process of LFRT, such as friction, melt shear, injection rate, etc. Especially for thin-walled products, owing to the relatively narrow mold cavity and frozen layer near the mold wall, it’s usually difficult for the fibers to move freely with the melt, leading to fiber breakage under the shear stress of the flow field eventually. In such a situation, the existing buckling breakage model based on the assumption fibers immersed in the flow field being in an unconstrained state is no longer applicable. In this paper, a new model for calculating the critical length of constrained fiber in thin-walled injection molded products is developed. A typical long glass fiber reinforced polypropylene thin-walled product was employed to verify the reliability of the model. The high-precision fiber length analyzer was used to calculate the fiber length distribution at different distances from the injection port. At the same time, the length distribution at the corresponding position was predicted via the model. The average prediction accuracy in the three sampling sections is 87.8%, 91.4% and 83.6%, respectively. The results show that the proposed model can effectively predict the fiber length distribution of injection molded thin-walled products. This conclusion provides important support for further accurate calculation of the mechanical properties of products.

Residual stress release and corresponding microstructural changes in high-energy impact-modified GH4169 after aging at 425 °C and 650 °C

The nickel-based superalloy is a critical material in aero-engines, where enhanced durability and resistance to stress relaxation are essential. This study evaluates the effectiveness of a high-energy composite modification technique, applied through dual surface friction methods, to induce residual stress and examines the stability of this stress under medium-temperature aging conditions. GH4169 was selected and a foundational study was conducted to quantify residual stress relaxation under medium-temperature aging treatment. The surface strengthening method used was a high-energy composite modification known for its significant post-treatment effects. The specific temperatures for medium-temperature aging treatment were 425 °C and 650 °C. The test results of the residual stress field confirmed the phenomenon of residual stress relaxation under the influence of temperature. Correspondingly, the different degrees of stress relaxation caused by varying aging treatment conditions are closely related to the microstructural changes within the material. Areas with a higher number of annealing twins exhibited a more pronounced degree of residual stress relaxation. By elucidating the link between stress relaxation and microstructural changes, this research offers a theoretical foundation for optimizing surface-strengthening processes and setting optimal operational temperatures for advanced structural materials.

Advances in defect characterization techniques using polarized light observation in SiC wafers for power devices

This short review provides an overview of advancements in the characterization of defects in silicon carbide (SiC) wafers using polarized light observation. SiC, which has a wide bandgap and excellent thermal and electrical properties, is widely used in power devices. However, various defects such as threading screw dislocations (TSD), threading edge dislocations (TED), and basal plane dislocations (BPD) can significantly affect device performance and reliability. Polarized light observation offers a nondestructive method for visualizing these defects and analyzing the stress fields induced within the SiC crystal structure. This paper summarizes recent developments in this technique, including the application of analyzer rotation to enhance the contrast in defect visualization. Furthermore, the development of automated systems for rapid wafer evaluation is discussed, highlighting the role of polarized light observation in improving quality control and production efficiency in SiC power device manufacturing.

Effects of the 2022 Mw 7.0 Chihshang Earthquake on the Shallow Creeping of the Central Longitudinal Valley Fault in Eastern Taiwan: Insights from Precise Leveling Measurements

Abstract This study investigated the effects of the 2022 Mw 7.0 Chihshang earthquake in eastern Taiwan on the shallow creeping of the central Longitudinal Valley fault (LVF). Precise leveling surveys were conducted along the Ruisui-Dewu, Yuli, and Dongzhu routes over the LVF. Before the earthquake, previous leveling surveys had detected deformations caused by stable fault creeping of the LVF. Interseismic linear trends and coseismic deformations were removed from the leveling data. Two shallow faults at 1 km were estimated from the corrected-coseismic deformation for the Yuli route. A shallow fault at 1 km and a deep fault at 5 km were estimated from the interseismic deformation. A positive delta Coulomb failure function change was calculated for shallow faults due to the creeping movement on the deep fault before the earthquake, and the creeping movement on the deep fault before the earthquake may have promoted the rupture of the shallow fault. The strong Chihshang earthquake most likely triggered the rupture of the shallow fault, on which strain had been concentrated by the interseismic creeping of the deeper fault. In contrast, the accumulated total strain at the deeper part might be quite low due to a strain release process for a steady creeping movement on deep fault in the interseismic period. The significant slip was not induced in the coseismic period. Although significant slip occurred in the shallower part, the slip might not be induced in the deeper part adjacent to the asperity of the LVF. In the Dewu route, the observed large west-side-uplift is attributed to the coseismic rupture of the Chihshang earthquake. In the Dongzhu route, an interseismic east-side-uplift by a creeping fault and a coseismic west-side-uplift were observed at the same location, indicating a local tensile stress field by the Chihshang earthquake.