Pub Date : 2026-02-01Epub Date: 2025-11-20DOI: 10.1016/j.ijsolstr.2025.113780
Guangfei Zhu, Rumeng Liu, Lifeng Wang
Strain engineering is a crucial approach for controlling the physical properties of van der Waals (vdW) materials. However, understanding the mechanical response of structures under strain remains a significant challenge. This study investigates the in-plane shear-induced wrinkling process of a bilayer vdW structure by molecular dynamics method. The mechanical mechanism is further validated with the assistance of a continuum-discrete model. The results demonstrated that out-of-plane deformation induces interlayer sliding, causing sliding instability. Interlayer sliding induces the formation of SP barrier rings composed of domain walls. The threefold rotational symmetry of the interlayer sliding breaks the symmetry of the wrinkling pattern. The shear direction also introduces anisotropic effects on wrinkle formation. The sliding instability of the twisted structure is correspondingly weakened. In addition, the onset of wrinkling and the transition timing during dynamic shear processes are significantly influenced by temperature and shear strain rate. These results provide new insights into the out-of-plane deformation mechanisms of layered vdW materials and offer novel approaches for strain engineering.
{"title":"Shear-induced anisotropic wrinkling and sliding instability in van der Waals Bilayers","authors":"Guangfei Zhu, Rumeng Liu, Lifeng Wang","doi":"10.1016/j.ijsolstr.2025.113780","DOIUrl":"10.1016/j.ijsolstr.2025.113780","url":null,"abstract":"<div><div>Strain engineering is a crucial approach for controlling the physical properties of van der Waals (vdW) materials. However, understanding the mechanical response of structures under strain remains a significant challenge. This study investigates the in-plane shear-induced wrinkling process of a bilayer vdW structure by molecular dynamics method. The mechanical mechanism is further validated with the assistance of a continuum-discrete model. The results demonstrated that out-of-plane deformation induces interlayer sliding, causing sliding instability. Interlayer sliding induces the formation of SP barrier rings composed of domain walls. The threefold rotational symmetry of the interlayer sliding breaks the symmetry of the wrinkling pattern. The shear direction also introduces anisotropic effects on wrinkle formation. The sliding instability of the twisted structure is correspondingly weakened. In addition, the onset of wrinkling and the transition timing during dynamic shear processes are significantly influenced by temperature and shear strain rate. These results provide new insights into the out-of-plane deformation mechanisms of layered vdW materials and offer novel approaches for strain engineering.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"326 ","pages":"Article 113780"},"PeriodicalIF":3.8,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145615475","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-02-01Epub Date: 2025-11-03DOI: 10.1016/j.ijsolstr.2025.113744
Duc Chung Vu
Despite the crucial role of breakable particles in numerous natural processes and industrial applications, accurately simulating particle breakage and its distinct variants remains a significant challenge for realistic discrete element method (DEM) simulations. In this work, we employ the bonded cell method (BCM), in which the particle is modeled as an aggregate of polyhedral cells whose common surfaces obey the Griffith fracture criterion, to simulate the fracture behavior of a single particle impacting a rigid plane. We focus on the influence of the restitution coefficient of cohesive bonds between cells on the particle fracture regimes and dissipated energy. We find that the crossover values of the damage potential , separating the three fracture regimes, follow a power-law relationship with the intercell restitution coefficient. Interestingly, at low values of corresponding to the first regime, where the particle undergoes the elastic rebound without crack formation, the effective restitution coefficient is independent of the intercell restitution coefficient. Our simulation data also reveal that the evolution of normalized energy dissipated by contact inelasticity and friction is well captured by a power-law function of . In the fragmented state, the power-law exponent is approximately 1 and remains independent of the intercell restitution coefficient, implying a linear dependence between the dissipated energy and the supplied kinetic energy. We show that the fraction of supplied energy lost to inelastic and frictional dissipation increases from about 50% to nearly 95% as the intercell restitution coefficient decreases. In contrast, in the damaged state, the power-law exponent is greater than 1 and decreases with decreasing intercell restitution coefficient. Finally, the dependence of several physical variables such as particle damage and fracture efficiency on the intercell restitution coefficient and impact velocity is also investigated.
{"title":"Impact-induced breakage of a single particle: Effect of the intercell restitution coefficient","authors":"Duc Chung Vu","doi":"10.1016/j.ijsolstr.2025.113744","DOIUrl":"10.1016/j.ijsolstr.2025.113744","url":null,"abstract":"<div><div>Despite the crucial role of breakable particles in numerous natural processes and industrial applications, accurately simulating particle breakage and its distinct variants remains a significant challenge for realistic discrete element method (DEM) simulations. In this work, we employ the bonded cell method (BCM), in which the particle is modeled as an aggregate of polyhedral cells whose common surfaces obey the Griffith fracture criterion, to simulate the fracture behavior of a single particle impacting a rigid plane. We focus on the influence of the restitution coefficient of cohesive bonds between cells on the particle fracture regimes and dissipated energy. We find that the crossover values of the damage potential <span><math><mi>ω</mi></math></span>, separating the three fracture regimes, follow a power-law relationship with the intercell restitution coefficient. Interestingly, at low values of <span><math><mi>ω</mi></math></span> corresponding to the first regime, where the particle undergoes the elastic rebound without crack formation, the effective restitution coefficient is independent of the intercell restitution coefficient. Our simulation data also reveal that the evolution of normalized energy dissipated by contact inelasticity and friction is well captured by a power-law function of <span><math><mi>ω</mi></math></span>. In the fragmented state, the power-law exponent is approximately 1 and remains independent of the intercell restitution coefficient, implying a linear dependence between the dissipated energy and the supplied kinetic energy. We show that the fraction of supplied energy lost to inelastic and frictional dissipation increases from about 50% to nearly 95% as the intercell restitution coefficient decreases. In contrast, in the damaged state, the power-law exponent is greater than 1 and decreases with decreasing intercell restitution coefficient. Finally, the dependence of several physical variables such as particle damage and fracture efficiency on the intercell restitution coefficient and impact velocity is also investigated.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"326 ","pages":"Article 113744"},"PeriodicalIF":3.8,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145428920","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This paper addresses the numerical implementation algorithm for an advanced anisotropic plasticity and damage continuum model. The framework of the proposed theory is based on the introduction of effective undamaged configurations, where no damage occurs, and the damaged configurations that account for elastic–plastic deformation and damage. The anisotropic plastic behavior is characterized by the Hoffman yield condition. The onset of damage is defined by a combination of the first and second deviatoric stress invariants related to the growth and coalescence of micro-defects (micro-voids and micro-shear-cracks). A stress-state-dependent damage strain rate tensor is introduced to capture the damage evolution caused by tension- and shear-induced mechanisms. The constitutive rate equations are numerically integrated using an explicit inelastic (plastic or plastic-damage) predictor-elastic corrector method. The consistent tangent modulus is derived and used to ensure quadratic convergence in the global finite element method. Moreover, numerical calculations for various biaxial loading conditions, including shear- and tension-induced damage mechanisms, demonstrate the accuracy and efficiency of the numerical algorithm. Numerical results are compared with experimental data at both the global load–displacement curve and the local strain fields, measured using the digital image correlation (DIC) technique. Scanning electron microscopy (SEM) is employed to compare the numerically predicted damage mechanism by examining fracture surfaces.
{"title":"Numerical analysis of anisotropic plasticity and damage based on the inelastic predictor-elastic corrector method","authors":"Zhichao Wei , Sanjeev Koirala , Steffen Gerke , Michael Brünig","doi":"10.1016/j.ijsolstr.2025.113770","DOIUrl":"10.1016/j.ijsolstr.2025.113770","url":null,"abstract":"<div><div>This paper addresses the numerical implementation algorithm for an advanced anisotropic plasticity and damage continuum model. The framework of the proposed theory is based on the introduction of effective undamaged configurations, where no damage occurs, and the damaged configurations that account for elastic–plastic deformation and damage. The anisotropic plastic behavior is characterized by the Hoffman yield condition. The onset of damage is defined by a combination of the first and second deviatoric stress invariants related to the growth and coalescence of micro-defects (micro-voids and micro-shear-cracks). A stress-state-dependent damage strain rate tensor is introduced to capture the damage evolution caused by tension- and shear-induced mechanisms. The constitutive rate equations are numerically integrated using an explicit inelastic (plastic or plastic-damage) predictor-elastic corrector method. The consistent tangent modulus is derived and used to ensure quadratic convergence in the global finite element method. Moreover, numerical calculations for various biaxial loading conditions, including shear- and tension-induced damage mechanisms, demonstrate the accuracy and efficiency of the numerical algorithm. Numerical results are compared with experimental data at both the global load–displacement curve and the local strain fields, measured using the digital image correlation (DIC) technique. Scanning electron microscopy (SEM) is employed to compare the numerically predicted damage mechanism by examining fracture surfaces.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"326 ","pages":"Article 113770"},"PeriodicalIF":3.8,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145569849","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-02-01Epub Date: 2025-11-07DOI: 10.1016/j.ijsolstr.2025.113748
X. Kuci, M.G.D. Geers, V.G. Kouznetsova
This work proposes a novel framework that combines dynamic computational homogenization with an extended transfer matrix method (TMM) to efficiently model wave propagation in locally resonant metamaterials (LRMs) with arbitrary microstructures. Unlike other methods in the literature, which assume specific symmetries and normal incidences, the presented approach addresses general multi-layered LRM setups for 2D and 3D wave propagation, including oblique incidences. First, the dynamic computational homogenization is applied to an LRM to extract the effective homogenized inertial and mechanical characteristics, yielding a macro-scale homogenized enriched continuum description. The enriched continuum description provides frequency-dependent properties, such as the effective dynamic impedance tensor, revealing wave attenuation behaviors near resonance frequencies. Wave propagation is then analyzed in multi-layered LRM setups with acoustic and/or elastic incoming media. A constrained dispersion equation is solved numerically to accurately model interface interactions without relying on analytical simplifications. The framework is validated against direct numerical simulations (DNS) through several representative case studies, demonstrating its versatility and significant computational efficiency. This novel approach paves the way for efficient wave impedance control and transmission analyses, providing new insights into the design and functionality of LRMs for advanced acoustic devices, such as acoustic filters and waveguides.
{"title":"Efficient wave analysis in multi-layered locally resonant metamaterials: A semi-analytical approach integrating dynamic homogenization","authors":"X. Kuci, M.G.D. Geers, V.G. Kouznetsova","doi":"10.1016/j.ijsolstr.2025.113748","DOIUrl":"10.1016/j.ijsolstr.2025.113748","url":null,"abstract":"<div><div>This work proposes a novel framework that combines dynamic computational homogenization with an extended transfer matrix method (TMM) to efficiently model wave propagation in locally resonant metamaterials (LRMs) with arbitrary microstructures. Unlike other methods in the literature, which assume specific symmetries and normal incidences, the presented approach addresses general multi-layered LRM setups for 2D and 3D wave propagation, including oblique incidences. First, the dynamic computational homogenization is applied to an LRM to extract the effective homogenized inertial and mechanical characteristics, yielding a macro-scale homogenized enriched continuum description. The enriched continuum description provides frequency-dependent properties, such as the effective dynamic impedance tensor, revealing wave attenuation behaviors near resonance frequencies. Wave propagation is then analyzed in multi-layered LRM setups with acoustic and/or elastic incoming media. A constrained dispersion equation is solved numerically to accurately model interface interactions without relying on analytical simplifications. The framework is validated against direct numerical simulations (DNS) through several representative case studies, demonstrating its versatility and significant computational efficiency. This novel approach paves the way for efficient wave impedance control and transmission analyses, providing new insights into the design and functionality of LRMs for advanced acoustic devices, such as acoustic filters and waveguides.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"326 ","pages":"Article 113748"},"PeriodicalIF":3.8,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145518000","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-02-01Epub Date: 2025-11-10DOI: 10.1016/j.ijsolstr.2025.113753
Amit Bhowmick , Xiang Gao , Wenquan Lu , Jeevanjyoti Chakraborty , Jun Xu
The growing interest in sodium-ion batteries (SIBs) is fueled by their abundant resources and environmentally friendly nature, with amorphous silicon (a-Si) emerging as a promising anode material for enhancing capacity. However, the key challenge lies in sustaining reversible capacity during cycling. In this work, we developed a multiscale electrochemical model incorporating an a-Si anode to elucidate the performance parameters of SIBs. Additionally, we integrated an electro-chemo-mechanical model at the particle level to capture stress generation, an essential factor in the degradation of high-capacity electrodes. Unlike existing models, our approach accounts for large-deformation chemo-mechanics at the particle scale and includes simulations under varying charge rates to explore multiscale behavior. The results reveal that coupled sodiation significantly prolongs complete cycling times and influences discharge dynamics, indicating that neglecting this coupling leads to an underestimation of actual capacity. Furthermore, we observed pronounced polarization effects at higher charge rates, resulting in heterogeneous stress distributions across the anode. With the identification of critical failure parameters for both active particles and binder materials, offering novel insights for mitigating degradation in high-capacity electrode systems.
{"title":"Multiscale electro-chemo-mechanical model of high-capacity amorphous silicon anode material in sodium-ion batteries","authors":"Amit Bhowmick , Xiang Gao , Wenquan Lu , Jeevanjyoti Chakraborty , Jun Xu","doi":"10.1016/j.ijsolstr.2025.113753","DOIUrl":"10.1016/j.ijsolstr.2025.113753","url":null,"abstract":"<div><div>The growing interest in sodium-ion batteries (SIBs) is fueled by their abundant resources and environmentally friendly nature, with amorphous silicon (a-Si) emerging as a promising anode material for enhancing capacity. However, the key challenge lies in sustaining reversible capacity during cycling. In this work, we developed a multiscale electrochemical model incorporating an a-Si anode to elucidate the performance parameters of SIBs. Additionally, we integrated an electro-chemo-mechanical model at the particle level to capture stress generation, an essential factor in the degradation of high-capacity electrodes. Unlike existing models, our approach accounts for large-deformation chemo-mechanics at the particle scale and includes simulations under varying charge rates to explore multiscale behavior. The results reveal that coupled sodiation significantly prolongs complete cycling times and influences discharge dynamics, indicating that neglecting this coupling leads to an underestimation of actual capacity. Furthermore, we observed pronounced polarization effects at higher charge rates, resulting in heterogeneous stress distributions across the anode. With the identification of critical failure parameters for both active particles and binder materials, offering novel insights for mitigating degradation in high-capacity electrode systems.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"326 ","pages":"Article 113753"},"PeriodicalIF":3.8,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145517999","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-02-01Epub Date: 2025-11-15DOI: 10.1016/j.ijsolstr.2025.113765
Yiming Chen, Dongqi An, Jinbao Li, Guangping Gong, Rui Li
Investigating the buckling behaviors of cracked plates carries substantial significance, since crack existence induces remarkable modifications to plate mechanical properties, potentially leading to significant degradation of structural load-carrying capability. This study develops a novel analytic solution framework that integrates the finite integral transform (FIT) method with an elementary domain decomposition strategy for solving buckling problems of single-edge-cracked rectangular thin plates. The through-thickness edge crack is modeled as an internal free boundary. The proposed framework exhibits universal applicability to plates with arbitrary combinations of simply supported, clamped, and free edges, and requires no assumptions regarding the form of the solutions throughout the derivation. The framework briefly comprises four key steps: decomposition of a single-edge-cracked plate into four elementary domains, followed by the application of a double cosine FIT to the governing equation of each domain; enforcement of all boundary and continuity conditions pertaining to Kirchhoff shear forces and rotations to eliminate a subset of the unknowns; substitution of inverse transforms into unapplied bending moment and deflection conditions to formulate the complete system of linear algebraic equations; determination of analytic solutions by solving the equations. Comprehensive buckling load/mode solutions of representative single-edge-cracked plates are presented as new benchmarks. A comparison of the solutions with other methods is conducted to validate the effectiveness of the FIT-based new solution framework. Utilizing the derived analytic solutions, a parametric study is conducted to quantitatively investigate the influences of boundary conditions, crack length ratio, crack location, and aspect ratio on the buckling behaviors.
{"title":"A unified analytic solution framework for buckling analysis of single-edge-cracked rectangular plates","authors":"Yiming Chen, Dongqi An, Jinbao Li, Guangping Gong, Rui Li","doi":"10.1016/j.ijsolstr.2025.113765","DOIUrl":"10.1016/j.ijsolstr.2025.113765","url":null,"abstract":"<div><div>Investigating the buckling behaviors of cracked plates carries substantial significance, since crack existence induces remarkable modifications to plate mechanical properties, potentially leading to significant degradation of structural load-carrying capability. This study develops a novel analytic solution framework that integrates the finite integral transform (FIT) method with an elementary domain decomposition strategy for solving buckling problems of single-edge-cracked rectangular thin plates. The through-thickness edge crack is modeled as an internal free boundary. The proposed framework exhibits universal applicability to plates with arbitrary combinations of simply supported, clamped, and free edges, and requires no assumptions regarding the form of the solutions throughout the derivation. The framework briefly comprises four key steps: decomposition of a single-edge-cracked plate into four elementary domains, followed by the application of a double cosine FIT to the governing equation of each domain; enforcement of all boundary and continuity conditions pertaining to Kirchhoff shear forces and rotations to eliminate a subset of the unknowns; substitution of inverse transforms into unapplied bending moment and deflection conditions to formulate the complete system of linear algebraic equations; determination of analytic solutions by solving the equations. Comprehensive buckling load/mode solutions of representative single-edge-cracked plates are presented as new benchmarks. A comparison of the solutions with other methods is conducted to validate the effectiveness of the FIT-based new solution framework. Utilizing the derived analytic solutions, a parametric study is conducted to quantitatively investigate the influences of boundary conditions, crack length ratio, crack location, and aspect ratio on the buckling behaviors.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"326 ","pages":"Article 113765"},"PeriodicalIF":3.8,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145569850","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-02-01Epub Date: 2025-11-04DOI: 10.1016/j.ijsolstr.2025.113742
Yuhang Zhang , Xiuming Liu , Yiqun Hu , Suhang Ding
Metallic materials with nanoscale voids have emerged as a novel class of structural and functional materials due to their unique properties. This study employs molecular dynamics simulations to investigate the strain-rate-dependent mechanical behaviors of Cu50Zr50 metallic glass stochastic network nanostructure (MGSNN) under uniaxial tension and compression over the strain rates from 5 × 106 to 5 × 109 s–1. It is found that the MGSNN exhibits strain-rate-independent elasticity and prominent strain-rate-sensitive plasticity. The Young’s modulus remains nearly constant, whereas the yield strength and ultimate tensile strength (UTS) significantly increase with increasing strain rate. This phenomenon originates from intrinsic deformation mechanisms: elastic response is governed by bond stretching that is inherently strain-rate-independent, whereas plastic deformation involves shear transformation zone activation and extension, which requires strain energy release and is strongly strain-rate-dependent. At higher strain rates, the release of energy is more difficult, and thus, the activation and appreciation of plastic events are restricted. As a result, the yield strength, yield strain, and UTS are enhanced at higher strain rates. The sequential yielding, necking, and breakage of individual nanowires are hindered at higher strain rates, resulting in the delayed global fracture of the MGSNN. A modified Gibson-Ashby relation comprising the strain rate effect precisely predicts the yield strength. The findings provide fundamental insights into the deformation mechanisms of amorphous porous nanostructures and establish guidelines for designing metallic glass nanofoams with tailored mechanical properties for structural and functional applications.
{"title":"Strain-rate-independent elasticity and strain-rate-sensitive plasticity in metallic glass stochastic network nanostructure","authors":"Yuhang Zhang , Xiuming Liu , Yiqun Hu , Suhang Ding","doi":"10.1016/j.ijsolstr.2025.113742","DOIUrl":"10.1016/j.ijsolstr.2025.113742","url":null,"abstract":"<div><div>Metallic materials with nanoscale voids have emerged as a novel class of structural and functional materials due to their unique properties. This study employs molecular dynamics simulations to investigate the strain-rate-dependent mechanical behaviors of Cu<sub>50</sub>Zr<sub>50</sub> metallic glass stochastic network nanostructure (MGSNN) under uniaxial tension and compression over the strain rates from 5 × 10<sup>6</sup> to 5 × 10<sup>9</sup> s<sup>–1</sup>. It is found that the MGSNN exhibits strain-rate-independent elasticity and prominent strain-rate-sensitive plasticity. The Young’s modulus remains nearly constant, whereas the yield strength and ultimate tensile strength (UTS) significantly increase with increasing strain rate. This phenomenon originates from intrinsic deformation mechanisms: elastic response is governed by bond stretching that is inherently strain-rate-independent, whereas plastic deformation involves shear transformation zone activation and extension, which requires strain energy release and is strongly strain-rate-dependent. At higher strain rates, the release of energy is more difficult, and thus, the activation and appreciation of plastic events are restricted. As a result, the yield strength, yield strain, and UTS are enhanced at higher strain rates. The sequential yielding, necking, and breakage of individual nanowires are hindered at higher strain rates, resulting in the delayed global fracture of the MGSNN. A modified Gibson-Ashby relation comprising the strain rate effect precisely predicts the yield strength. The findings provide fundamental insights into the deformation mechanisms of amorphous porous nanostructures and establish guidelines for designing metallic glass nanofoams with tailored mechanical properties for structural and functional applications.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"326 ","pages":"Article 113742"},"PeriodicalIF":3.8,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145464478","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-02-01Epub Date: 2025-11-23DOI: 10.1016/j.ijsolstr.2025.113785
Qi Yao , Stephan Rudykh
This paper demonstrates the deformation-activated negative group velocity (NGV) state in soft fiber-reinforced composites (soft FCs). We exhibit the mechanism of inducing the NGV state in the vicinity of buckling. Based on extensive simulations, we construct a map of NGV intervals versus deformation level and fiber spacing ratio, revealing clear trends and tunable behavior. In the transition region of buckling direction, NGV states can be selectively activated for shear waves with different polarizations. The presence and order of NGV polarizations are found to correlate with the direction of post-buckling fiber morphology. NGV states in both polarizations indicate non-principal, helical buckling patterns and complex evolutions. The relative width of NGV intervals of two polarizations mirrors the magnitude of helical deformation in each direction. These findings suggest NGV behavior as a predictive indicator for post-buckling configurations. This work offers a new pathway for understanding NGV states in soft FCs, with potential applications in reconfigurable materials.
{"title":"Deformation-activated polarization-selective negative group velocity in rectangularly-arranged fiber-reinforced soft composite","authors":"Qi Yao , Stephan Rudykh","doi":"10.1016/j.ijsolstr.2025.113785","DOIUrl":"10.1016/j.ijsolstr.2025.113785","url":null,"abstract":"<div><div>This paper demonstrates the deformation-activated negative group velocity (NGV) state in soft fiber-reinforced composites (soft FCs). We exhibit the mechanism of inducing the NGV state in the vicinity of buckling. Based on extensive simulations, we construct a map of NGV intervals versus deformation level and fiber spacing ratio, revealing clear trends and tunable behavior. In the transition region of buckling direction, NGV states can be selectively activated for shear waves with different polarizations. The presence and order of NGV polarizations are found to correlate with the direction of post-buckling fiber morphology. NGV states in both polarizations indicate non-principal, helical buckling patterns and complex evolutions. The relative width of NGV intervals of two polarizations mirrors the magnitude of helical deformation in each direction. These findings suggest NGV behavior as a predictive indicator for post-buckling configurations. This work offers a new pathway for understanding NGV states in soft FCs, with potential applications in reconfigurable materials.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"326 ","pages":"Article 113785"},"PeriodicalIF":3.8,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145615477","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-02-01Epub Date: 2025-11-23DOI: 10.1016/j.ijsolstr.2025.113782
Yizhou Huang , Zhiyang Yao , Cheng Qiu , Yile Hu , Jinglei Yang
Understanding failure mechanisms of microcapsules is crucial for self-healing efficiency. However, complex failure modes of micrpcapsules makes numerical simulations challenging. Moreover, experimental methods struggle to accurately explore how mechanical properties influence microcapsule failure mechanisms. In this study, the enhanced peridynamic (PD) method with correcting cross-material interface bonds under non-uniform discretization was developed for microcapsule-modified concrete. Tensile tests are conducted on models containing various material components, and the PD calculation results agree well with finite element method results, validating the accuracy of the proposed method. Subsequently, a parametric analysis is performed to investigate the effects of interface strength matching between the microcapsule shell and mortar on the failure mode. PD results show that stronger interface ensures microcapsule rupture rather than debonding. Additionally, the study of characteristic core-to-shell ratios (ranging from 20:1 to 1:1.5) demonstrated that the shell thickness has little impact on whether a crack can penetrate the microcapsule. Finally, composite models reveal that increasing volume fraction of microcapsules leads to a decrease in stiffness and strength due to preferential crack propagation pathways provided by microcapsules. Furthermore, the method proposed in this study successfully captures all the failure modes observed in the CT scans of the compression test, further validating the effectiveness of the proposed approach. The method proposed in this paper is a new and effective tool for analyzing microcapsule-modified concretes.
{"title":"Failure mechanism analysis of microcapsule-modified concrete via state-based peridynamics","authors":"Yizhou Huang , Zhiyang Yao , Cheng Qiu , Yile Hu , Jinglei Yang","doi":"10.1016/j.ijsolstr.2025.113782","DOIUrl":"10.1016/j.ijsolstr.2025.113782","url":null,"abstract":"<div><div>Understanding failure mechanisms of microcapsules is crucial for self-healing efficiency. However, complex failure modes of micrpcapsules makes numerical simulations challenging. Moreover, experimental methods struggle to accurately explore how mechanical properties influence microcapsule failure mechanisms. In this study, the enhanced peridynamic (PD) method with correcting cross-material interface bonds under non-uniform discretization was developed for microcapsule-modified concrete. Tensile tests are conducted on models containing various material components, and the PD calculation results agree well with finite element method results, validating the accuracy of the proposed method. Subsequently, a parametric analysis is performed to investigate the effects of interface strength matching between the microcapsule shell and mortar on the failure mode. PD results show that stronger interface ensures microcapsule rupture rather than debonding. Additionally, the study of characteristic core-to-shell ratios (ranging from 20:1 to 1:1.5) demonstrated that the shell thickness has little impact on whether a crack can penetrate the microcapsule. Finally, composite models reveal that increasing volume fraction of microcapsules leads to a decrease in stiffness and strength due to preferential crack propagation pathways provided by microcapsules. Furthermore, the method proposed in this study successfully captures all the failure modes observed in the CT scans of the compression test, further validating the effectiveness of the proposed approach. The method proposed in this paper is a new and effective tool for analyzing microcapsule-modified concretes.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"326 ","pages":"Article 113782"},"PeriodicalIF":3.8,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145615473","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Concrete damage estimation is significantly influenced by the inherent stochastic flaws within its microstructure, yet conventional models often oversimplify these complexities. Aggregates’ location and geometric attributes play a crucial role in defining fracture locality preference. This study introduces a refined three-phase mesoscale model that integrates geometric irregularities and spatial distributions of aggregate morphology, enabling more accurate simulations of damage evolution under static and impact loading. Uncertainty quantification is incorporated to assess the probabilistic variations in flaw distributions across the cement matrix, interfacial transition zones (ITZ), and aggregates to provide deeper insights into failure mechanisms with an analysis of computational complexity. Numerical analysis reveals the critical role of Weibull-idealized population flaw sizes in determining structural failure reliability, highlighting the interplay between aggregate characteristics and relative impedances of mortar and ITZ. It provided better resolutions in fracture processing zones, crack initiation locations, and crack patterns at various strain levels. The findings contribute to enhancing predictive capabilities, offering a more robust framework for optimizing concrete design and resilience under diverse loading conditions.
{"title":"Uncertainty-driven damage estimation of concrete: Integrating realistic aggregate morphology and stochastic flaw behaviour","authors":"Abhinov Bharadwaj , Mohit Kr Sharma , Sukanta Chakraborty , Debojit Biswas","doi":"10.1016/j.ijsolstr.2025.113775","DOIUrl":"10.1016/j.ijsolstr.2025.113775","url":null,"abstract":"<div><div>Concrete damage estimation is significantly influenced by the inherent stochastic flaws within its microstructure, yet conventional models often oversimplify these complexities. Aggregates’ location and geometric attributes play a crucial role in defining fracture locality preference. This study introduces a refined three-phase mesoscale model that integrates geometric irregularities and spatial distributions of aggregate morphology, enabling more accurate simulations of damage evolution under static and impact loading. Uncertainty quantification is incorporated to assess the probabilistic variations in flaw distributions across the cement matrix, interfacial transition zones (ITZ), and aggregates to provide deeper insights into failure mechanisms with an analysis of computational complexity. Numerical analysis reveals the critical role of Weibull-idealized population flaw sizes in determining structural failure reliability, highlighting the interplay between aggregate characteristics and relative impedances of mortar and ITZ. It provided better resolutions in fracture processing zones, crack initiation locations, and crack patterns at various strain levels. The findings contribute to enhancing predictive capabilities, offering a more robust framework for optimizing concrete design and resilience under diverse loading conditions.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"326 ","pages":"Article 113775"},"PeriodicalIF":3.8,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145615467","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}