Strain Energy Research Articles

Comparison of Experimental and Numerical Investigation of Mono-Composite and Metal Leaf Spring

The automotive industry is increasingly focused on reducing vehicle weight, leading to the widespread adoption of composite materials with high strength-to-weight ratios in both aviation and automotive sectors. These materials are gradually replacing traditional options like steel. Leaf springs, one of the oldest and most common suspension components, continue to be widely used in vehicles. This study aims to replace conventional multi-leaf steel springs with mono-composite leaf springs while preserving the same load-carrying capacity and stiffness. Composite materials, such as glass fiber and epoxy resin, provide advantages including higher elastic strain energy storage, superior strength-to-weight ratios, and enhanced corrosion resistance compared to steel. Consequently, the weight of leaf springs can be reduced without sacrificing performance. The steel and mono-composite leaf springs were modeled using Catia software, and their performance was evaluated using ANSYS 15.0 software.

Differential growth and shape formation of a flower-shaped structure

Morphogenesis, which is a complex interplay of biological, chemical, and physical processes, facilitates differential growth in various biological systems, particularly in plant organs, such as petals and leaves. Although recent studies have increasingly delved into the mechanical aspects of these biological structures, there remains a lack of a comprehensive quantitative understanding of shape formation induced by differential growth in plant organs. Thus, this study addressed this gap by employing a multifaceted approach encompassing theoretical analyses and comprehensive finite element analysis of the effect of differential growth on the shape of a cylinder into a flower shape. Based on the derivation of the strain energy expressions for the axisymmetric and asymmetric configurations, the shape function of the growing flower-like structures were calculated through mathematical optimization. The findings of this study shed light on the influence of the growth function, geometric characteristics, and material properties. As the wave number increased, the final configuration tended to have smaller waves, whereas a longer cylinder buckled more easily and the thickness had a minimal effect. This study offers insights that can pave the way for innovative geometrical designs, thereby providing inspiration for applications on both micro-and macroscales, such as in the realms of self-assembly of soft robotics and flexible electronics.

Study on the effects of electrode volume fraction and electrode particle radius on dynamic electrochemical-thermal-aging coupling characteristics of lithium cell based on cylindrical cell design method

Electrode volume fraction and electrode particle radius are important electrode parameters for the lithium cell design and deeply influence the electrochemical-physical processes inside lithium cell, which is still insufficiently clear. In order to obtain the effects of electrode parameters on lithium cell, a combine method of the dynamic electrochemical-thermal-aging coupling model and cylindrical cell design was proposed. Simulation was conducted at charge rate 1C. Cell temperature, current density and Solid-Electrolyte-Interface at ambient temperature 25 °C, and total strain energy density and lithium plating at ambient temperature −5 °C were analyzed. As electrode volume fraction or electrode particle radius increases, maximum temperature and maximum temperature difference of lithium cell increase. When electrode volume fraction increases from 0.3 to 0.7, maximum temperature of cell increases by 2.6 °C at least. Current density peak increases with the increase of the positive electrode volume fraction or electrode particle radius. When electrode volume fraction increases from 0.3 to 0.7, current density peak increases by 44 %∼74 %. As electrode volume fraction increases, total strain energy density Ws decreases in the early stage of charge, while increases in the late stage of charge. As the negative (or positive) electrode particle radius increases, Ws of the negative (or positive) electrode increase, while Ws of the positive (or negative) electrode decreases. Solid-Electrolyte-Interface thickness increases with the increase of the electrode volume fraction or electrode particle radius. The effect range of negative electrode volume on lithium plating is 0.5 ∼ 0.7. The lithium metal thickness increase with the increase of the positive electrode volume fraction. As the negative or positive electrode particle radius increases, the lithium metal thickness increase or decrease, respectively.

Evolution of Atomic-Level Interfacial Fracture Mechanics in Magnesium-Zinc Compounds Used for Bioresorbable Vascular Stents.

Bioresorbable magnesium-metal vascular stents are gaining popularity due to their biodegradable nature and good biological and mechanical properties. They are also suitable candidate materials for biodegradable stents. Due to the rapid degradation rate of Mg metal vascular scaffolds, a Mg/Zn bilayer composite was formed by a number of means, such as magnetron sputtering and physical vapor deposition, thus delaying the degradation time of the Mg metal vascular scaffolds while providing good radial support for the stenotic vessels. However, the interlaminar compounds at the metal interface have an essential impact on the mechanical properties of the bi-material interface, especially the cracking and delamination of the Mg matrix Zn coating vascular stent in the radially expanded process layer. Intermetallic compounds (IMCs) are commonly found in dual-layer composites, such as Mg/Zn composites and multi-layer structures. They are frequently overlooked in simulations aiming to predict mechanical properties. This paper analyses the interfacial failure processes and evolutionary mechanisms of interfacial fracture mechanics of a Mg/Zn interface with an intermetallic compound layer between coated Zn and Mg matrix metallic vascular stents. The simulation results show that the fracture mode in the Mg/Zn interface with an intermetallic compound involves typical ductile fracture under static tensile conditions. The dislocation line defects mainly occur on the side of the Mg, which induces the Mg/Zn interfacial crack to expand along the interface into the pure Mg. The stress intensity factor and the critical strain energy release rate decrease as the intermetallic compound layer's thickness gradually increases, indicating that the intensity of stress and the force of the crack extending and expanding along the crack tip are weakened. The presence of intermetallic compounds at the interface can significantly strengthen the mechanical properties of the material interface and alleviate the crack propagation between the interfaces.

Phenyl, propyl polysiloxane modified epoxy resin II: The co-continuous phase, concurrently strengthening and toughening mechanism

The applications of epoxy resins are restricted in certain areas due to their brittleness, while improving the epoxy toughness always compromises their strength. We reported a phenyl, propyl polysiloxane modified epoxy resins with excellent tensile strength and toughness before. In this study, the mechanism of strengthening and toughening are studied and elaborated. Compared to cured EP resin (EP), cured P/E resin (P/E20) significantly improved the tensile strength, elongation at break, impact strength, and critical strain energy release rate (GIc) by 45 %, 248 %, 18 %, and 512 %, respectively. Fracture analysis of the P/E20 revealed the appearance of parabolic markings, typical and convincing evidence of ductile fracture, indicating good toughness of the resin. AFM height and Young's modulus plots reveal the complex details of the dimples and the co-continuous phases, which provide the strong evidence to explain the reinforcement and toughening mechanism.

Allylic Epoxides Increase the Strain Energy of Cyclic Olefin Monomers for Ring-Opening Metathesis Polymerization.

Ring-opening metathesis polymerization (ROMP) is an effective method for synthesizing functional polymers, but since the technique typically relies on high ring strain cyclic olefins, the most common monomers are norbornene derivatives. The reliance on one class of monomer limits the obtainable properties of ROMP polymers. In this work, we investigate new bicyclic monomers synthesized via epoxidation of commercial dienes. DFT estimates of these monomers' ring strains suggests a significant increase in strain for cyclic olefins containing allylic epoxides. We found that the eight-membered (3,4-COO) and five-membered (CPO) cyclic olefins were particularly effective for ROMP. CPO was of especially intriguing due to its excellent polymerizability when compared to the limited reactivity of other five-membered rings. Unlike polynorbornenes, the resulting polymers of both monomers displayed glass transition temperatures well below room temperature. Interestingly, poly(3,4-COO) showed both high stereo- and regioregularity while poly(CPO) showed little regularity. Both polymers could be readily modified via post-polymerization ring-opening of the reactive allylic epoxides. With a high epoxide density in poly(CPO), CPO is an exciting new ROMP monomer that is easily synthesized, can be polymerized to high conversion at room temperature, and may be facilely modified to yield a wide range of functional materials.

Allylic Epoxides Increase the Strain Energy of Cyclic Olefin Monomers for Ring‐Opening Metathesis Polymerization

Ring‐opening metathesis polymerization (ROMP) is an effective method for synthesizing functional polymers, but since the technique typically relies on high ring strain cyclic olefins, the most common monomers are norbornene derivatives. The reliance on one class of monomer limits the obtainable properties of ROMP polymers. In this work, we investigate new bicyclic monomers synthesized via epoxidation of commercial dienes. DFT estimates of these monomers’ ring strains suggests a significant increase in strain for cyclic olefins containing allylic epoxides. We found that the eight‐membered (3,4‐COO) and five‐membered (CPO) cyclic olefins were particularly effective for ROMP. CPO was of especially intriguing due to its excellent polymerizability when compared to the limited reactivity of other five‐membered rings. Unlike polynorbornenes, the resulting polymers of both monomers displayed glass transition temperatures well below room temperature. Interestingly, poly(3,4‐COO) showed both high stereo‐ and regioregularity while poly(CPO) showed little regularity. Both polymers could be readily modified via post‐polymerization ring‐opening of the reactive allylic epoxides. With a high epoxide density in poly(CPO), CPO is an exciting new ROMP monomer that is easily synthesized, can be polymerized to high conversion at room temperature, and may be facilely modified to yield a wide range of functional materials.

Investigating Mechanical Characteristics of Rocks Under Freeze–Thaw Cycles Using Grain-Based Model

AbstractBased on the discrete element method (DEM), a water-contained grain-based model (GBM) is developed in this study to evaluate the effects of freeze–thaw cycles (FTCs) on the mechanical characteristics of the rock. A set of freeze–thaw and uniaxial compression tests is carried out to explore the impact of micro damage caused by FTCs on the mechanical prosperities of rock samples. By monitoring the development and distribution of micro-cracks during freeze–thaw test and uniaxial compression test, the damage mechanism of FTCs is revealed from a microscopic perspective, which shows that FTCs deteriorate the strength and brittleness parameters as exponential functions. The parametric analysis is carried out to explore the influence of porosity and mineral components on the mechanical behaviors of rock against freeze–thaw and uniaxial loadings. Based on the parametric analysis results, it is found that UCS, Young’s modulus, and total strain energy at peak stress decrease with the increase of porosity and clay content, which emphasizes the contributions of porosity and mineral components on the mechanical properties of rock samples. It is proved that the water-contained grain-based model developed in this study can capture the damage caused by FTCs on the mechanical performance of rock from a microscopic perspective, which provides novel perspectives on the phenomenon of rock degradation in response to fluctuations in temperature.

Fault slip and seismic source response characteristics under induced stress wave

As coal mining depth and structural complexity increase, the failure severity of mining-induced fault activation leading to rock bursts escalates. Fault fracture zones, consisting of various secondary structures, are essential for understanding tectonic evolution and seismic wave propagation. The interaction between mining-induced stress waves and in-situ stress complicates the dynamics and activation markers of fault fracture zones. This paper aims to clarify the patterns of stick–slip instability and stress evolution in these zones under coal mining-induced stress waves. Using theoretical models and continuous-discrete coupling simulations, the study explores the dynamic response of fault fracture zones in the meta-instability stage, the micro-scale fracture force chain structure, and the distribution and synergistic instability of fault plane seismic sources. The research shows that fault fracture zones changes stress wave propagation path and energy redistribution due to its heterogeneity leads to multiple times decomposition of wave field. Minor delamination between the hanging wall and footwall can trigger unstable slip. Tensional seismic sources with extensive failure drive fault slip. Dynamic stress waves disrupt force chains near seismic sources, increasing hazard levels. Static in-situ stress affects shear seismic sources by influencing porosity and crack angles. Seismic sources along the fault strike transition from divergent to concentrated, trending horizontally. In the meta-instability stage, the initial fracture area locks, with strong structures causing cohesive fractures in the middle of the fault zone, releasing strain energy and forming a sudden increase in slip intensity. These findings provide crucial insights into fault fracture-development-activation markers, fault region stress deformation evolution mechanisms, and risk control of fault-induced rock bursts.

Nonlinear dynamics of wing-like structures using a momentum subspace-based Koiter-Newton reduction

Cantilevers find a wide range of applications in the design of scientific equipment and large-scale engineering structures such as aircraft wings. Analysis techniques based on linearization approximations are unable to capture the large amplitude oscillation behaviour of such structures and thus, necessitates development of dedicated nonlinear methods. In this work, the recent developments in the Koiter-Newton model reduction method are utilized to obtain nonlinear reduced order models (ROMs) from full finite element structural models in order to simulate large amplitude dynamics of cantilevers. The method describes a nonlinear system of governing equations comprising quadratic and cubic terms which are obtained as higher order derivatives of the in-plane strain energy. To ensure that the large rotations in cantilevers and the resultant foreshortening effect is also accounted for, a ROM updating algorithm is adopted where the ROM parameters are varied with the structural deflections. Linear eigenmodes of the structure are utilized to form the reduction subspace. To validate the methodology, the ROM solution is compared against experimental results and a convergence study is conducted to identify the number of modes needed to replicate the nonlinear response. Finally, a composite wingbox structure is considered for which time domain simulations are conducted and frequency response curves, obtained through a frequency sweep, are presented.

Fatigue life predictions for welded boiler water walls

Boiler water walls experience in-phase thermo-mechanical loading during start-ups and shut-downs, leading to low cycle fatigue (LCF) failure. This study aims at establishing an FEA-based failure prediction method for estimating the fatigue performance and service life of the welded water walls. The developed model is validated for predicting failure in the defect-free uniaxial fatigue specimens. Stress-mechanical strain hysteresis loops and accumulated inelastic strain energy density per cycle parameters are extracted from fatigue tests at 0.4%, 0.6%, and 0.7% strains. A combination of cyclic plasticity and continuous damage mechanics (CDM) theory is utilized to predict fatigue crack initiation sites and estimate the specimen fatigue life. Accumulated damage has been calculated for the life cycle of each specimen. FEA model predicted failure and service life agrees well with the experimental results. The established failure analysis parameters are then transferred from the specimen level to the water wall component level, thereby estimating the service life of defect-free water walls at 750 cycles.

Molecular Chain Rearrangement-Induced In Situ Formation of Nanofibers for Improving the Strength and Toughness of Hydrogels.

Although poly(vinyl alcohol) (PVA) hydrogel has high elasticity and is suitable for cartilage tissue engineering, it is difficult to have both high strength and toughness. In this study, a simple and universal strategy is proposed to prepare strong and tough PVA hydrogels by in situ forming nanofibers on the original network structure induced by a molecular chain rearrangement. Quenching-tempering alteratively in ethanol and water several times is carried out to strengthen PVA hydrogels (PVA-Etn hydrogels) due to the advantages of noncovalent bonds in adjustability and reversibility. The results show that, after three quenching-tempering cycles, PVA-Et3 hydrogel with water content up to 79 wt % shows comprehensive improved mechanical properties. The compression modulus, tensile modulus, fracture strength, tensile strain, and tear energy of the PVA-Et3 hydrogel are 270, 250, 260, 130, and 180% of the initial PVA hydrogel, respectively. The improved mechanical properties of the PVA-Et3 hydrogel are attributed to the strong cross-linked PVA chains and hydrogen bond-reinforced nanofibers. This study not only provides a simple and efficient solution for the preparation of strong and tough polymer scaffolds in tissue engineering but also provides new insights for understanding the mechanism of improving the mechanical properties of polymer hydrogels by adjusting the molecular structure.

Damping optimization of a fully viscoelastic structure with enforced mass reduction

Current damping methods employed in industry, such as Constrained Layer Damping (CLD), all require the introduction of additional mass and a secondary material to the structure. In transportation industries, this weight addition is often correlated with an increase in operational costs. Additionally, in applications with large temperature changes, such as in space, the secondary material introduces an additional source of stress during thermal expansion and contraction. This work presents a methodology to significantly improve the damping performance of a fully viscoelastic, single material structure while also enforcing a mass reduction. The work expands upon the modal strain energy (MSE) method to calculate the damping performance of a fully viscoelastic structure and develops a novel weighting scheme to make the optimization more robust. The methodology is then implemented on a technology demonstration thermoplastic lunar rover. This implementation demonstrates the ability to optimize a fully viscoelastic structure to significantly improve the damping performance while also achieving a large decrease in mass. The implementation also demonstrates the robustness of the weighting scheme as it allows the natural frequencies to shift between iterations. The skin thickness of the rover’s base panel was optimized. After optimization, the base panel of the rover achieved a 20% damping improvement for Mode 1 with a 9% mass reduction. The skin thickness of the full rover was optimized and achieved a 15% damping improvement for Mode 2 with an 11% mass reduction.

Bound-constrained optimization using Lagrange multiplier for a length scale insensitive phase field fracture model

The classical phase field model using the second-order geometric function α(φ)=φ2 (i.e., AT2 model), where φ∈0,1 is an auxiliary phase field variable representing material damage state, has wide applications in static and dynamic scenarios for brittle materials, but nonlinearity and inelasticity are found in its stress–strain curve. The phase field model using the linear geometric function α(φ)=φ (i.e., AT1 model), can avoid this, and a linear elastic threshold is available in its stress–strain curve. However, both AT2 and AT1 models are length scale sensitive phase field models, which could have difficulty in adjusting fracture strength and crack band simultaneously through a single parameter (the length scale). In this paper, a generalized quadratic geometric function (linear combination of AT1 and AT2 models) is used in the phase field model, where the extra parameter in this geometric function makes it a length scale insensitive phase field model. Similar to the AT1 model, negative phases can happen in the proposed generalized quadratic geometric function model. To solve this problem, a bound-constrained optimization using the Lagrange multiplier is derived, and the Karush–Kuhn–Tucker (KKT) conditions change from strain energy and maximum history strain energy (an indirect method acting on phase) to phase and Lagrange multiplier (a direct method acting on phase). Several simulations successfully validated the proposed model. A single element analysis and a bar under cyclic loading show the different stress–strain curves obtained from different models. A simulation of Mode I Brazilian test is compared with the experiment conducted by the authors, and two more simulations of Mode II shear test and mixed mode PMMA tensile test are compared with results from the literature.

Predicting fatigue life of additively manufactured lattice structures using the image-based Finite Cell Method and average strain energy density

A well-known Achilles' heel of laser powder bed fusion (L-PBF) additive manufactured lattice structures is the difficulty in predicting fatigue properties. The presence of manufacturing-induced defects significantly affects the fatigue resistance of the porous component and must be accurately captured by predictive models. To tackle this challenge, the as-built geometry of the lattice needs to be modeled, which introduces another challenge on the computational front. For the first time, a model based on the computer tomography (μ-CT) reconstruction of the as-built lattice geometry is simulated with the efficient finite cell method and combined with the average strain energy density (ASED) to obtain accurate fatigue predictions. This work presents a workflow for determining the fatigue resistance of lattice metamaterial, followed by a case study for method validation. The validation shows a good agreement between the predicted fatigue life and the experimental results.

A comparative study of low cyclic loading effects on plastic strain, nonlinear damping, and strength softening in fiber-reinforced recycled aggregate concrete

The incorporation of fibers significantly enhances the properties and structural performance of recycled aggregate concrete. This study introduces Fiber-Reinforced Recycled Aggregate Concrete (FRAC) by integrating steel fibers (SF) and polypropylene fibers (PPF) into the recycled aggregate concrete (RAC) matrix. A comprehensive series of cyclic loading tests, including reload and unload sequences, were conducted to meticulously analyze the deterioration of strength and energy dissipation capabilities of SF/PPF-reinforced RAC. The research examines the effects of fiber characteristics, volume ratios, lengths, and diameters on strength degradation after the peak stress along the stress-strain curve. Degradation models based on plastic strain were proposed for both SF-reinforced RAC and PPF-reinforced RAC. The study further investigates the fluctuation of strain energy in FRAC under low cyclic loading, revealing the correlation between plastic strain energy and equivalent viscous damping ratio. A nonlinear viscous damping model that accounts for the influence of the reinforcing index (RI) is proposed. This model effectively predicts the degradation of strength, the evolution of plastic strain, and the dynamics of nonlinear viscous damping in composite materials with varying fiber content. Ultimately, this research provides essential technical insights for the practical engineering application of RAC.

Experimental and numerical study on the effect of friction between crack faces on fracture parameters of hot mix asphalt concrete

This article investigates the effect of friction between crack faces in ASCB (asymmetric semi-circular bend) specimen on fracture parameters of hot mix asphalt concrete. To achieve this, fracture tests were conducted on ASCB specimens containing cracks with different angles (0, 10, 20, 30, and 40°). In addition to the fracture tests, friction tests were also performed by designing a simple fixture. Subsequently, finite element analyses were carried out, and friction coefficients between crack faces were calculated based on the results of fracture and friction tests. Furthermore, two fracture criteria, ASED (average minimum strain energy density) and MTS (maximum tangential stress), were examined under two different scenarios: assuming zero-friction and non-zero friction between crack faces to predict fracture load and fracture strength. The fracture tests showed that with an increase in the crack angle, the fracture load initially decreased up to the crack angle of 20°, and then increased. Furthermore, the friction tests revealed that an increase in the normal load led to an increase in the friction coefficient. As per finite element results, the crack face contact occurred at a crack angle of 13.1°. Hence, crack faces in the ASCB specimens with the crack angles of 20, 30, and 40° had friction, which significantly affected the fracture parameters. The results also indicated that the average error values for the MTS and ASED criteria compared to the experimental results were 22.8 % and 18.8 %, respectively, when the friction` between crack faces is ignored, and 12.4 % and 4.0 % for the scenario that the friction effect is taken into account. The ASED criterion demonstrated higher accuracy than MTS criterion. Additionally, the average error value of fracture strength predicted by the ASED criterion with assuming non-zero friction was 4 %, indicating good predictive accuracy for fracture strength.

Analysis of high temperature and strain amplitude effects on low cycle fatigue behavior of pitting corroded killed E350 BR structural steel

High-tension structural steels are prone to accelerated fatigue damage from pitting corrosion and high-temperature. Despite adverse effects, research on their low cycle fatigue (LCF) behavior is limited. Specifically, studies analyzing temperature-dependent pit sensitivity effects, considering pit-related material’s susceptibility to surface topographic features variation and stress concentration are lacking. This study conducts LCF tests on pitting corroded killed E350 BR structural steel at multiple strain amplitudes and high temperatures. It develops temperature-dependent parameters, such as the cyclic softening pit sensitivity factor, suitable for integration into existing approaches like total cyclic plastic strain energy density (CPSED), power-law, average strain energy density (SED), Coffin-Manson, and pit stress intensity factor (pit-SIF). Further, it proposes multiple linear regression-based prediction models relating total CPSED and average SED with strain amplitude and temperature. Corroded specimens show higher plastic deformation and reduced peak stress, fatigue life, and total CPSED compared to uncorroded ones. The developed parameters, integrated with average SED approach, predicts LCF life within an error band of ±1.5, while power-law relationship reduces it to ±1.2. Moreover, pit-SIF approach estimates fatigue life within an error band of ±1.5. The findings provide critical knowledge for enhanced component design, leading to structural safety, performance, and fire resilience.

Vibration of composite open shell of hydrogen-electric fuselage with rectangular cutout in hygrothermal circumstances: theoretical and experimental research

An effective formulation is performed to approximate calculate the vibration of the typical composite fuselage battery compartment structures of the hydrogen-electric aircraft, which are highly sensitive to hygrothermal circumstances. The strain energy, kinetic energy and hygrothermal potential energy of CLOCSRC are deduced in frame of Donnell's shell assumption and a series of hypothetical elastic springs operating on the edges of divided components are introduced to simulate coupling connection and the elastic supports of the open shell, which contribute to total system energy in the form of elastic spring potential energy. The equations of motion for CLOCSRC subjected to distinct hygrothermal circumstances and arbitrary boundaries are eventually obtained in frame of the Rayleigh-Ritz method by means of the modified orthogonal polynomial as the displacement admissible function. After confirming the convergence and accuracy of the presented formulation by experiment and finite element method, several numerical examples are presented to investigative the effects of geometric parameters, elastic boundaries and hygrothermal circumstances on the vibration of CLOCSRC. Differ from traditional investigation, the present paper offers an effective approach to insight the hygrothermal mechanism of open shell with rectangular cutouts creatively. The relevant investigation demonstrated that the current strategy is beneficial for further research on the vibration characteristics of CLOCSRC.

Enhanced mechanical properties and atomic-scale mechanisms of ferroelastic domain switching for GdNbO4-La2Zr2O7 materials

Optimization of mechanical properties in La2Zr2O7 (LZO) ceramics, composites and coatings is an on-going requirement for their practical application. Herein, the contribution of monoclinic (La, Gd)NbO4 (m-LGNO) enhancement of fracture toughness by ∼56% reveals its capability to be a prominent toughening agent. Due to the ferroelastic nature of LGNO, ferroelastic switching takes place within the stress concentrated regions, giving rise to significant strain energy relaxation. Atomic-scale evidence reveals that ferroelastic 94°/86° domain switching can occur, yielding merged 94° domains and newly formed 86° domains. The relevant strains induced by ferroelastic domain switching are quantified up to 8.06% and 6.20% in shear and normal strain, respectively. Such domain switching strains highlight their contribution to accommodate external mechanical loading for the 50 mol% GNO-LZO composite. The results indicate that the unique ferroelastic nature and 94°/86° ferroelastic domain switching in m-LGNO cooperatively provide a significant toughening effect.