Pub Date : 2026-02-10DOI: 10.1016/j.ijsolstr.2026.113895
Lu Meng , Ji Qiu , Pei Li , Heng Zhang , Zhe Liu
The coupled thermo-electro-structural effect is a key factor affecting the performance and reliability of electrical devices. As the number of coupled physical fields increases, the computational complexity grows nonlinearly, rendering existing methods inadequate for achieving concurrent multiscale modeling of tri-field coupling, which severely constrains the optimal design for multiscale electrical devices under complex operating conditions. To this end, a coupled thermal-electro-structural direct finite element square (D-FE2) method is developed to perform full coupled multiphysics concurrent multiscale computation in heterogeneous materials and structures. Numerical results show that the proposed method not only accurately captures the Joule heating-induced thermal expansion phenomena in three-dimensional heterogeneous materials and structures under various temperature fields (i.e., steady-state and transient) and electrical excitations (i.e., direct current and alternating current), but also achieves over 91% reduction in intrinsic computational workload and above 95% savings in practical resource consumption compared to conventional direct numerical simulation (DNS) methods. This study provides a breakthrough solution for concurrent multiscale modeling and optimization of complex coupled multiphysics systems in large-scale electrical devices.
{"title":"Multiscale modeling and analysis of coupled thermo-electro-structural behavior in heterogeneous materials and structures with direct FE2 method","authors":"Lu Meng , Ji Qiu , Pei Li , Heng Zhang , Zhe Liu","doi":"10.1016/j.ijsolstr.2026.113895","DOIUrl":"10.1016/j.ijsolstr.2026.113895","url":null,"abstract":"<div><div>The coupled thermo-electro-structural effect is a key factor affecting the performance and reliability of electrical devices. As the number of coupled physical fields increases, the computational complexity grows nonlinearly, rendering existing methods inadequate for achieving concurrent multiscale modeling of tri-field coupling, which severely constrains the optimal design for multiscale electrical devices under complex operating conditions. To this end, a coupled thermal-electro-structural direct finite element square (D-FE<sup>2</sup>) method is developed to perform full coupled multiphysics concurrent multiscale computation in heterogeneous materials and structures. Numerical results show that the proposed method not only accurately captures the Joule heating-induced thermal expansion phenomena in three-dimensional heterogeneous materials and structures under various temperature fields (i.e., steady-state and transient) and electrical excitations (i.e., direct current and alternating current), but also achieves over 91% reduction in intrinsic computational workload and above 95% savings in practical resource consumption compared to conventional direct numerical simulation (DNS) methods. This study provides a breakthrough solution for concurrent multiscale modeling and optimization of complex coupled multiphysics systems in large-scale electrical devices.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"331 ","pages":"Article 113895"},"PeriodicalIF":3.8,"publicationDate":"2026-02-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146154252","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}
Chiral honeycombs are demonstrated to exhibit exceptional out–of–plane mechanical properties owing to their distinctive chirality. However, systematic investigations, particularly theoretical studies, remain limited. In this work, a comprehensive study is conducted on the out‑of‑plane crashworthiness of a prototypical chiral topology, namely the wavy tri‑chiral honeycomb (WTCH). The configuration of the WTCH can be interpreted as a transformation of the conventional hexagonal honeycomb (CHH), achieved by substituting its straight cell walls with curved segments. Finite element (FE) simulations are first performed to elucidate the out-of-plane crashworthiness and underlying deformation mechanism of the WTCH. Under quasi-static loading condition, three distinct collapse modes, i.e., progressive folding, concertina-like buckling, and global bending, are identified, governed by a trade-off between wave slope and cell wall slenderness ratio. The latter two modes are characterized by significantly elevated peak stresses, rendering them unfavorable for energy-absorbing applications. In contrast, under dynamic loading conditions, all WTCH configurations exhibit a progressive folding collapse mode. Based on the identified deformation mechanism, theoretical models for predicting the plateau stresses of the WTCH undergoing progressive folding collapse are developed using the Simplified Super Folding Element (SSFE) theory and validated against numerical simulations, showing strong agreement. Finally, a quantitative comparison with the CHH demonstrates that the WTCH achieves at least a 47 % improvement in specific energy absorption under identical cell size and relative density, along with consistently higher crush force efficiency. These findings provide valuable insights for the structural optimization and engineering applications of wavy chiral honeycombs in crash protection and safety contexts.
{"title":"Theoretical modeling of out-of-plane crashworthiness of a wavy tri-chiral honeycomb","authors":"Yilin Zhu , Miao Zhou , Xue Rui , Chao Yu , Yifeng Zhong , Chuanzeng Zhang","doi":"10.1016/j.ijsolstr.2026.113873","DOIUrl":"10.1016/j.ijsolstr.2026.113873","url":null,"abstract":"<div><div>Chiral honeycombs are demonstrated to exhibit exceptional out–of–plane mechanical properties owing to their distinctive chirality. However, systematic investigations, particularly theoretical studies, remain limited. In this work, a comprehensive study is conducted on the out‑of‑plane crashworthiness of a prototypical chiral topology, namely the wavy tri‑chiral honeycomb (WTCH). The configuration of the WTCH can be interpreted as a transformation of the conventional hexagonal honeycomb (CHH), achieved by substituting its straight cell walls with curved segments. Finite element (FE) simulations are first performed to elucidate the out-of-plane crashworthiness and underlying deformation mechanism of the WTCH. Under quasi-static loading condition, three distinct collapse modes, i.e., progressive folding, concertina-like buckling, and global bending, are identified, governed by a trade-off between wave slope and cell wall slenderness ratio. The latter two modes are characterized by significantly elevated peak stresses, rendering them unfavorable for energy-absorbing applications. In contrast, under dynamic loading conditions, all WTCH configurations exhibit a progressive folding collapse mode. Based on the identified deformation mechanism, theoretical models for predicting the plateau stresses of the WTCH undergoing progressive folding collapse are developed using the Simplified Super Folding Element (SSFE) theory and validated against numerical simulations, showing strong agreement. Finally, a quantitative comparison with the CHH demonstrates that the WTCH achieves at least a 47 % improvement in specific energy absorption under identical cell size and relative density, along with consistently higher crush force efficiency. These findings provide valuable insights for the structural optimization and engineering applications of wavy chiral honeycombs in crash protection and safety contexts.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"329 ","pages":"Article 113873"},"PeriodicalIF":3.8,"publicationDate":"2026-01-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146090419","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-01-23DOI: 10.1016/j.ijsolstr.2026.113866
Jean-Emmanuel Leroy, Valentin L. Popov
The paper builds upon a method previously developed and introduced by the authors for the normal contact of a rigid smooth indenter with a poroelastic half-space extending it to a single layer. Fundamental solutions in the Fourier and Laplace space for two types of the layer support, bonded and unbonded, are presented. It is shown how using these solutions discrete kernels can be constructed, enabling FFT-accelerated computation of the layer response to applied pressure distributions in both space and time, and therefore efficient solution of contact problems with indenters of arbitrary shape.
{"title":"Boundary element method for normal contacts of poroelastic layers","authors":"Jean-Emmanuel Leroy, Valentin L. Popov","doi":"10.1016/j.ijsolstr.2026.113866","DOIUrl":"10.1016/j.ijsolstr.2026.113866","url":null,"abstract":"<div><div>The paper builds upon a method previously developed and introduced by the authors for the normal contact of a rigid smooth indenter with a poroelastic half-space extending it to a single layer. Fundamental solutions in the Fourier and Laplace space for two types of the layer support, bonded and unbonded, are presented. It is shown how using these solutions discrete kernels can be constructed, enabling FFT-accelerated computation of the layer response to applied pressure distributions in both space and time, and therefore efficient solution of contact problems with indenters of arbitrary shape.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"329 ","pages":"Article 113866"},"PeriodicalIF":3.8,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146090418","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-01-23DOI: 10.1016/j.ijsolstr.2026.113865
Xiaojun Tan , Jian Ma , Shaohua Liu , Bo Cao , Xueyan Chen , Bing Wang , Muamer Kadic
High-performance and reusable energy-absorbing materials have tremendous potential in industrial applications. Achieving both high performance and reusability has long been a challenge due to their apparent incompatibility. To address this, we proposed a solid–liquid dual-state mechanical metamaterial. This metamaterial exhibits robust mechanical properties when the liquid metal is solid and achieves high energy absorption through its plastic deformation. Upon heating-induced solid–liquid state transition, its deformed state fully recovers its initial state, ensuring reusability. The metamaterials can be fabricated by injecting liquid metal into an hollow elastic lattice structure manufactured through additive manufacturing processes. The mechanical properties of solid–liquid dual-state mechanical metamaterials prepared from different liquid metal, such as gallium, Field’s metal, and Wood’s metal, are analyzed in this paper through a combined approach of experiments, theoretical analysis, and numerical simulations. The results reveal that the proposed metamaterial significantly outperforms all previously reported reusable energy-absorbing materials in specific energy absorption (SEA). This breakthrough driven by the solid–liquid state transition redefines the limits of reusable energy absorption and opens the path to develop a complete family of robust, reusable materials.
{"title":"Design, fabrication, and characterization of solid–liquid dual-state mechanical metamaterials","authors":"Xiaojun Tan , Jian Ma , Shaohua Liu , Bo Cao , Xueyan Chen , Bing Wang , Muamer Kadic","doi":"10.1016/j.ijsolstr.2026.113865","DOIUrl":"10.1016/j.ijsolstr.2026.113865","url":null,"abstract":"<div><div>High-performance and reusable energy-absorbing materials have tremendous potential in industrial applications. Achieving both high performance and reusability has long been a challenge due to their apparent incompatibility. To address this, we proposed a solid–liquid dual-state mechanical metamaterial. This metamaterial exhibits robust mechanical properties when the liquid metal is solid and achieves high energy absorption through its plastic deformation. Upon heating-induced solid–liquid state transition, its deformed state fully recovers its initial state, ensuring reusability. The metamaterials can be fabricated by injecting liquid metal into an hollow elastic lattice structure manufactured through additive manufacturing processes. The mechanical properties of solid–liquid dual-state mechanical metamaterials prepared from different liquid metal, such as gallium, Field’s metal, and Wood’s metal, are analyzed in this paper through a combined approach of experiments, theoretical analysis, and numerical simulations. The results reveal that the proposed metamaterial significantly outperforms all previously reported reusable energy-absorbing materials in specific energy absorption (SEA). This breakthrough driven by the solid–liquid state transition redefines the limits of reusable energy absorption and opens the path to develop a complete family of robust, reusable materials.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"329 ","pages":"Article 113865"},"PeriodicalIF":3.8,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146036727","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-01-22DOI: 10.1016/j.ijsolstr.2026.113871
Wabi Demeke , Sangryun Lee , Wonju Jeon , Seunghwa Ryu
Metamaterials composed of gradient-index (GRIN) phononic crystals (PnCs), which contain unit cells whose sizes vary perpendicular to the direction of wave propagation, serve as a means of focusing elastic waves for energy harvesting. Owing to the finite wavelength of the propagating wave, GRIN PnCs localize the wave within a finite-sized region, which experiences multimodal strains instead of focusing on a single focal point with a single mode of strain. Consequently, the energy harvesting of elastic waves across a localized, finite region has recently gained research interest. This is due to the advantage of harvesting more energy through a properly designed piezoelectric energy harvester (PEH) that is larger than the wavelength of the elastic wave. However, the design of unit cells in GRIN PnC has been predominantly limited to simple shapes. This study introduces random hole shapes in GRIN PnC to enhance the intensity of elastic energy localization across targeted finite-sized region, utilizing a data-efficient surrogate model through Bayesian optimization (BO). Additionally, the developed BO method identifies a unit cell design that offers wave focusing intensity comparable to that of benchmark deep neural network (DNN)-based optimization, while requiring only 5.9% of the dataset. This advancement in wave localization significantly enhances the wave localization intensity in the target region by 36% and improves power generation by up to 1.5 times compared to GRIN PnC design with simple circular hole.
{"title":"Enhancing elastic energy focusing in multimode strain regions via Bayesian optimization of gradient-index phononic crystals for energy harvesting","authors":"Wabi Demeke , Sangryun Lee , Wonju Jeon , Seunghwa Ryu","doi":"10.1016/j.ijsolstr.2026.113871","DOIUrl":"10.1016/j.ijsolstr.2026.113871","url":null,"abstract":"<div><div>Metamaterials composed of gradient-index (GRIN) phononic crystals (PnCs), which contain unit cells whose sizes vary perpendicular to the direction of wave propagation, serve as a means of focusing elastic waves for energy harvesting. Owing to the finite wavelength of the propagating wave, GRIN PnCs localize the wave within a finite-sized region, which experiences multimodal strains instead of focusing on a single focal point with a single mode of strain. Consequently, the energy harvesting of elastic waves across a localized, finite region has recently gained research interest. This is due to the advantage of harvesting more energy through a properly designed piezoelectric energy harvester (PEH) that is larger than the wavelength of the elastic wave. However, the design of unit cells in GRIN PnC has been predominantly limited to simple shapes. This study introduces random hole shapes in GRIN PnC to enhance the intensity of elastic energy localization across targeted finite-sized region, utilizing a data-efficient surrogate model through Bayesian optimization (BO). Additionally, the developed BO method identifies a unit cell design that offers wave focusing intensity comparable to that of benchmark deep neural network (DNN)-based optimization, while requiring only 5.9% of the dataset. This advancement in wave localization significantly enhances the wave localization intensity in the target region by 36% and improves power generation by up to 1.5 times compared to GRIN PnC design with simple circular hole.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"330 ","pages":"Article 113871"},"PeriodicalIF":3.8,"publicationDate":"2026-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146070913","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-01-20DOI: 10.1016/j.ijsolstr.2026.113856
Zhiyuan Li, Yuang Zhang, Wang Zhong, Yiqun Zhang, Zihan Sun, Naigang Hu
Deployable membrane antennas are extensively utilized in space exploration due to their high packing efficiency and lightweight properties. Prior to orbital deployment, these antennas must undergo repeated ground-based deployment tests to verify their reliability and performance during deployment. However, repeated deployment aggravates damage at the membrane creases, compromising structural integrity and deployment accuracy. To address this issue, this paper establishes an endurance degradation model for membrane creases under repeated folding based on experimental data, and derives a quantitative relationship between the crease characteristic dimension (i.e., crease angle) and the number of folding cycles. A deployment analysis model for the membrane structure is then developed based on the Flasher folding method, incorporating the effects of both repeated folding and initial bending stress at the creases on deployment behavior. Through the stress-bending moment relationship, it is demonstrated that an increasing number of folding cycles leads to higher initial bending stress at the creases, resulting in a significant reduction in deployment accuracy. Concurrently, the energy required for the deployment process, i.e., the driving force, increases substantially. Finally, repeated folding-deployment experiments were conducted on a prototype planar membrane antenna. The experimental results validate the crease endurance degradation model and confirm the mechanistic influence of repetitive folding on the unfolding process.
{"title":"Research on the repeated folding mechanism of membrane antennas based on crease endurance degradation","authors":"Zhiyuan Li, Yuang Zhang, Wang Zhong, Yiqun Zhang, Zihan Sun, Naigang Hu","doi":"10.1016/j.ijsolstr.2026.113856","DOIUrl":"10.1016/j.ijsolstr.2026.113856","url":null,"abstract":"<div><div>Deployable membrane antennas are extensively utilized in space exploration due to their high packing efficiency and lightweight properties. Prior to orbital deployment, these antennas must undergo repeated ground-based deployment tests to verify their reliability and performance during deployment. However, repeated deployment aggravates damage at the membrane creases, compromising structural integrity and deployment accuracy. To address this issue, this paper establishes an endurance degradation model for membrane creases under repeated folding based on experimental data, and derives a quantitative relationship between the crease characteristic dimension (i.e., crease angle) and the number of folding cycles. A deployment analysis model for the membrane structure is then developed based on the Flasher folding method, incorporating the effects of both repeated folding and initial bending stress at the creases on deployment behavior. Through the stress-bending moment relationship, it is demonstrated that an increasing number of folding cycles leads to higher initial bending stress at the creases, resulting in a significant reduction in deployment accuracy. Concurrently, the energy required for the deployment process, i.e., the driving force, increases substantially. Finally, repeated folding-deployment experiments were conducted on a prototype planar membrane antenna. The experimental results validate the crease endurance degradation model and confirm the mechanistic influence of repetitive folding on the unfolding process.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"329 ","pages":"Article 113856"},"PeriodicalIF":3.8,"publicationDate":"2026-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146036724","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}
The present work investigates the dynamic behaviour of viscoelastic composite shafts by developing a detailed mathematical model for woven fabric composites. Frequency-dependent properties of the epoxy resin are obtained using a Dynamic Mechanical Analyzer (DMA) and subsequently curve-fitted to derive a viscoelastic operator through a Genetic Algorithm–based optimization scheme. The key novelty of this study lies in proposing an operator-based composite modulus formulation, wherein elastic fibers are integrated into the viscoelastic matrix at the constitutive level, enabling a frequency-dependent representation of woven composites. This integration yields a theoretical operator modulus for the composite based on the rule of mixture concept that naturally introduces higher-order terms into the governing equations of motion. The proposed theoretical operator modulus is further experimentally validated by fabricating and testing a woven fabric composite in DMA. The validated model is then employed to analyze the dynamic behaviour of a viscoelastic shaft–rotor system. The influence of fiber orientation on system stability is examined, revealing a symmetric stability pattern, with minimum stability at 45° and maximum stability at 0° and 90° orientations relative to the spin axis. This behaviour is attributed to the alignment of fibers along principal stiffness directions. The findings underscore the critical importance of fiber orientation and provide a robust, experimentally supported modelling framework for designing dynamically stable composite shafts.
{"title":"Mathematical modelling and dynamic analysis of visco-elastic woven fabric composite shafts","authors":"Aditya Sharma , Krishanu Ganguly , Rajan Prasad , Jayanta Kumar Dutt","doi":"10.1016/j.ijsolstr.2026.113855","DOIUrl":"10.1016/j.ijsolstr.2026.113855","url":null,"abstract":"<div><div>The present work investigates the dynamic behaviour of viscoelastic composite shafts by developing a detailed mathematical model for woven fabric composites. Frequency-dependent properties of the epoxy resin are obtained using a Dynamic Mechanical Analyzer (DMA) and subsequently curve-fitted to derive a viscoelastic operator through a Genetic Algorithm–based optimization scheme. The key novelty of this study lies in proposing an operator-based composite modulus formulation, wherein elastic fibers are integrated into the viscoelastic matrix at the constitutive level, enabling a frequency-dependent representation of woven composites. This integration yields a theoretical operator modulus for the composite based on the rule of mixture concept that naturally introduces higher-order terms into the governing equations of motion. The proposed theoretical operator modulus is further experimentally validated by fabricating and testing a woven fabric composite in DMA. The validated model is then employed to analyze the dynamic behaviour of a viscoelastic shaft–rotor system. The influence of fiber orientation on system stability is examined, revealing a symmetric stability pattern, with minimum stability at 45° and maximum stability at 0° and 90° orientations relative to the spin axis. This behaviour is attributed to the alignment of fibers along principal stiffness directions. The findings underscore the critical importance of fiber orientation and provide a robust, experimentally supported modelling framework for designing dynamically stable composite shafts.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"329 ","pages":"Article 113855"},"PeriodicalIF":3.8,"publicationDate":"2026-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146090333","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}
Nanoindentation has been developed to map the local mechanical properties of heterogeneous solids. Although the local properties at the indentation site are required for mapping, the map resolution worsens because of the deformation of other components surrounding the indentation, such as the grains encompassing an indentation in heterogeneous polycrystals. To quantify the transitional change in a P-h curve with an increase in h, where P is the indentation load and h is the indentation penetration depth, this study proposes weight functions. These functions aim to quantify point-sharp nanoindentations on a circular column embedded in a matrix, which has an interface parallel to the indentation direction. The mechanical properties of the column can be accurately evaluated by extrapolating the weight functions to an indentation depth of zero, which improves the map resolution. When comparing the effects of the parallel and vertical interfaces, the impact of the vertical interface on the P-h curve is more significant than that of the parallel interface. This suggests that the correlation between the heterogeneity of local mechanical properties and the microstructure observed at the surface is controversial.
{"title":"Computational study of nanoindentation on an elastoplastic solid with an interface parallel to the indentation direction","authors":"Takashi Akatsu , Wataru Kubota , Yutaka Shinoda , Fumihiro Wakai","doi":"10.1016/j.ijsolstr.2026.113857","DOIUrl":"10.1016/j.ijsolstr.2026.113857","url":null,"abstract":"<div><div>Nanoindentation has been developed to map the local mechanical properties of heterogeneous solids. Although the local properties at the indentation site are required for mapping, the map resolution worsens because of the deformation of other components surrounding the indentation, such as the grains encompassing an indentation in heterogeneous polycrystals. To quantify the transitional change in a <em>P</em>-<em>h</em> curve with an increase in <em>h</em>, where <em>P</em> is the indentation load and <em>h</em> is the indentation penetration depth, this study proposes weight functions. These functions aim to quantify point-sharp nanoindentations on a circular column embedded in a matrix, which has an interface parallel to the indentation direction. The mechanical properties of the column can be accurately evaluated by extrapolating the weight functions to an indentation depth of zero, which improves the map resolution. When comparing the effects of the parallel and vertical interfaces, the impact of the vertical interface on the <em>P</em>-<em>h</em> curve is more significant than that of the parallel interface. This suggests that the correlation between the heterogeneity of local mechanical properties and the microstructure observed at the surface is controversial.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"329 ","pages":"Article 113857"},"PeriodicalIF":3.8,"publicationDate":"2026-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146036728","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-01-18DOI: 10.1016/j.ijsolstr.2026.113853
Hongzhao Li , Raj Kiran , Yu Tan , Haidong Fan , Qingyuan Wang , Peidong Li
Quasicrystals (QCs) have demonstrated excellent hydrogen storage capabilities, making them promising candidates for hydrogen-related applications. However, exposure to hydrogen can significantly affect their mechanical integrity. In this study, a phase-field fracture model (PFM) is proposed to predict fracture behaviors of polycrystalline QCs in a hydrogen-containing environment. The damage-field and interface-field variables are introduced to regularize the cracks and interfaces (or grain boundaries) in QCs, respectively. An equivalent critical energy release rate is introduced to characterize the influence of the interfaces (or grain boundaries) on hydrogen-assisted crack propagation. The present model is numerically implemented in Comsol Multiphysics based on the Weak Form PDE module. Several numerical examples are conducted to validate the ability of the proposed model to predict the fracture behaviors of QCs and to analyze the influence of hydrogen concentration on the fracture mechanism of QCs. Numerical results indicate that the hydrogen concentration and interfaces (or grain boundaries) significantly influence the crack propagation paths and the peak force and failure displacement in the fracture test of QCs. The existence of hydrogen decreases the critical fracture load and promotes the crack propagation in QCs significantly. The developed phase-field framework and numerical implementation approach provide a convenient tool for predicting fracture failure and assessing the service reliability of polycrystalline QCs in complex hydrogen-containing environments.
{"title":"Phase-field fracture modeling of polycrystalline quasicrystals subjected to hydrogen embrittlement","authors":"Hongzhao Li , Raj Kiran , Yu Tan , Haidong Fan , Qingyuan Wang , Peidong Li","doi":"10.1016/j.ijsolstr.2026.113853","DOIUrl":"10.1016/j.ijsolstr.2026.113853","url":null,"abstract":"<div><div>Quasicrystals (QCs) have demonstrated excellent hydrogen storage capabilities, making them promising candidates for hydrogen-related applications. However, exposure to hydrogen can significantly affect their mechanical integrity. In this study, a phase-field fracture model (PFM) is proposed to predict fracture behaviors of polycrystalline QCs in a hydrogen-containing environment. The damage-field and interface-field variables are introduced to regularize the cracks and interfaces (or grain boundaries) in QCs, respectively. An equivalent critical energy release rate is introduced to characterize the influence of the interfaces (or grain boundaries) on hydrogen-assisted crack propagation. The present model is numerically implemented in <span>Comsol</span> Multiphysics based on the W<span>eak</span> F<span>orm</span> PDE module. Several numerical examples are conducted to validate the ability of the proposed model to predict the fracture behaviors of QCs and to analyze the influence of hydrogen concentration on the fracture mechanism of QCs. Numerical results indicate that the hydrogen concentration and interfaces (or grain boundaries) significantly influence the crack propagation paths and the peak force and failure displacement in the fracture test of QCs. The existence of hydrogen decreases the critical fracture load and promotes the crack propagation in QCs significantly. The developed phase-field framework and numerical implementation approach provide a convenient tool for predicting fracture failure and assessing the service reliability of polycrystalline QCs in complex hydrogen-containing environments.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"329 ","pages":"Article 113853"},"PeriodicalIF":3.8,"publicationDate":"2026-01-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146090407","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-01-17DOI: 10.1016/j.ijsolstr.2026.113844
Rayehe Karimi Mahabadi, Hossein Salahshoor
While multimaterial additive manufacturing enables finely programmed heterogeneity, there remains no robust and objective-driven framework to assign materials across complex architectures under practical constraints. We introduce Data-Driven-Design () as a robust computational framework for multi-material lattice design, optimized with respect to a prescribed performance objective. The framework relies on representing physical constraints, material data, and design objectives as sets in a phase space and formulating the material selection problem as a distance minimization problem among the encoded sets. We showcase the approach in multi-material design of viscoelastic lattices provided with measurements of complex moduli as a function of frequency with the design objective of maximizing dissipation. For our numerical experiments, we import dynamic viscoelasticity measurement for twenty five different materials from literature, and show that multi-material designs can match or outperform the dissipation obtained from homogeneous designs made of the most dissipative material among the data registry. In a finite lattice example, we show that design yields a mechanical dissipation with 300% increase compared to best homogeneous design from a limited collection of materials. Beyond viscoelastic lattices, the framework generalizes naturally to multi-physics and multi-objective metastructure design, offering a unified, data-driven approach to optimal material selection under complex constraints.
{"title":"Data-Driven-Design (D3) of multi-material systems: A novel framework and its application to viscoelastic metamaterials","authors":"Rayehe Karimi Mahabadi, Hossein Salahshoor","doi":"10.1016/j.ijsolstr.2026.113844","DOIUrl":"10.1016/j.ijsolstr.2026.113844","url":null,"abstract":"<div><div>While multimaterial additive manufacturing enables finely programmed heterogeneity, there remains no robust and objective-driven framework to assign materials across complex architectures under practical constraints. We introduce Data-Driven-Design (<span><math><msup><mrow><mi>D</mi></mrow><mrow><mn>3</mn></mrow></msup></math></span>) as a robust computational framework for multi-material lattice design, optimized with respect to a prescribed performance objective. The framework relies on representing physical constraints, material data, and design objectives as sets in a phase space and formulating the material selection problem as a distance minimization problem among the encoded sets. We showcase the approach in multi-material design of viscoelastic lattices provided with measurements of complex moduli as a function of frequency with the design objective of maximizing dissipation. For our numerical experiments, we import dynamic viscoelasticity measurement for twenty five different materials from literature, and show that multi-material designs can match or outperform the dissipation obtained from homogeneous designs made of the <em>most dissipative material</em> among the data registry. In a finite lattice example, we show that <span><math><msup><mrow><mi>D</mi></mrow><mrow><mn>3</mn></mrow></msup></math></span> design yields a mechanical dissipation with 300% increase compared to best homogeneous design from a limited collection of materials. Beyond viscoelastic lattices, the <span><math><msup><mrow><mi>D</mi></mrow><mrow><mn>3</mn></mrow></msup></math></span> framework generalizes naturally to multi-physics and multi-objective metastructure design, offering a unified, data-driven approach to optimal material selection under complex constraints.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"329 ","pages":"Article 113844"},"PeriodicalIF":3.8,"publicationDate":"2026-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146090417","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}