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}
Pub Date : 2026-01-16DOI: 10.1016/j.ijsolstr.2026.113854
Gregoire Markey , Etienne Barthel , Catalin R. Picu
Glass wool is used extensively for thermal and acoustic insulation in buildings, appliances, automotive and aerospace applications. In this work we study the response of fiberglass networks crosslinked using polymeric binders to small amplitude oscillatory excitation in the low acoustic frequency range. Two types of binders are used, with low and high damping, and a range of fiber densities and binder volume fractions are considered. We report the scaling of storage and loss moduli of the network with these parameters. Further, we evaluate the loss factor and find that it is independent of the fiber mass density but is power law related to the binder volume fraction. We develop a model that provides a mechanistic basis for the experimental results reported and may guide material design.
{"title":"Mechanical dissipation in random fiberglass networks under small strains","authors":"Gregoire Markey , Etienne Barthel , Catalin R. Picu","doi":"10.1016/j.ijsolstr.2026.113854","DOIUrl":"10.1016/j.ijsolstr.2026.113854","url":null,"abstract":"<div><div>Glass wool is used extensively for thermal and acoustic insulation in buildings, appliances, automotive and aerospace applications. In this work we study the response of fiberglass networks crosslinked using polymeric binders to small amplitude oscillatory excitation in the low acoustic frequency range. Two types of binders are used, with low and high damping, and a range of fiber densities and binder volume fractions are considered. We report the scaling of storage and loss moduli of the network with these parameters. Further, we evaluate the loss factor and find that it is independent of the fiber mass density but is power law related to the binder volume fraction. We develop a model that provides a mechanistic basis for the experimental results reported and may guide material design.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"329 ","pages":"Article 113854"},"PeriodicalIF":3.8,"publicationDate":"2026-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146036750","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-16DOI: 10.1016/j.ijsolstr.2026.113833
Thao Nguyen, Darby J. Luscher, Miles Buechler, Michael B. Prime, Saryu J. Fensin
The original Tepla (TEnsile PLAsticity) ductile damage model, based on the Gurson yield surface, has long been used to model damage evolution and material failure under dynamic loading. Unfortunately, Tepla suffered from mesh sensitivity, numerical instability, and limited predictive capability. Here, we theoretically reformulate Tepla to address these issues. We especially focus on the prediction of porosity, which is the key state variable used for modeling ductile damage, as compared to more easily measured surface velocities, which are at best an indirect measure of damage. Key model changes include separating the viscosity during volumetric void growth from underlying shear strength behavior and switching to an iterative bisection solver. The new Tepla is then calibrated on incipient spall experiments on half-hard copper and tantalum, which demonstrate its ability to simultaneously fit the model to recovered porosity distributions and measured surface velocities, a stringent test. Improved numerical behavior, such as greatly reduced mesh sensitivity, is also shown in those simulations. Finally, the new Tepla model is applied to several high-explosive loaded, sweeping wave experiments, showing the ability of the model to predict behavior on tests with significantly different loading conditions and histories than the calibration data.
{"title":"Reformulating a Gurson-based dynamic damage model and demonstrating improved predictive power and numerical robustness","authors":"Thao Nguyen, Darby J. Luscher, Miles Buechler, Michael B. Prime, Saryu J. Fensin","doi":"10.1016/j.ijsolstr.2026.113833","DOIUrl":"10.1016/j.ijsolstr.2026.113833","url":null,"abstract":"<div><div>The original Tepla (TEnsile PLAsticity) ductile damage model, based on the Gurson yield surface, has long been used to model damage evolution and material failure under dynamic loading. Unfortunately, Tepla suffered from mesh sensitivity, numerical instability, and limited predictive capability. Here, we theoretically reformulate Tepla to address these issues. We especially focus on the prediction of porosity, which is the key state variable used for modeling ductile damage, as compared to more easily measured surface velocities, which are at best an indirect measure of damage. Key model changes include separating the viscosity during volumetric void growth from underlying shear strength behavior and switching to an iterative bisection solver. The new Tepla is then calibrated on incipient spall experiments on half-hard copper and tantalum, which demonstrate its ability to <em>simultaneously</em> fit the model to recovered porosity distributions and measured surface velocities, a stringent test. Improved numerical behavior, such as greatly reduced mesh sensitivity, is also shown in those simulations. Finally, the new Tepla model is applied to several high-explosive loaded, sweeping wave experiments, showing the ability of the model to predict behavior on tests with significantly different loading conditions and histories than the calibration data.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"329 ","pages":"Article 113833"},"PeriodicalIF":3.8,"publicationDate":"2026-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146036726","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-14DOI: 10.1016/j.ijsolstr.2026.113852
Sichen Liu, Wenbin Yu
In practical engineering applications, slender structures often have complex geometries, such as aircraft wings and rotor blades. Slender structures with low aspect ratios were often analyzed using three-dimensional (3D) Finite Element Analysis (FEA) as beam models were widely considered inadequate. However, analyzing slender structures with complex geometries was usually impractical due to the large number of degrees of freedom (DOFs) needed. This underscored the need to evaluate whether beam models can preserve accuracy while substantially reducing computational cost. To perform this assessment, a semi-analytical solution of the Timoshenko beam model, based on Variational Asymptotic Beam Sectional Analysis (VABS) considering all possible couplings, was used to predict the natural frequencies of general composite structures in free vibration analysis. The natural frequencies obtained from the semi-analytical solution were compared against the results from 3D FEA to assess the accuracy of the VABS-based Timoshenko model. The findings indicated that, for general composite beams, the Timoshenko model provided accurate predictions of the first several modes, excluding non-beam modes, mitigating the need for resource-intensive 3D FEA.
{"title":"Semi-analytical solution of VABS-based Timoshenko Beam Model for free vibration of composite structures","authors":"Sichen Liu, Wenbin Yu","doi":"10.1016/j.ijsolstr.2026.113852","DOIUrl":"10.1016/j.ijsolstr.2026.113852","url":null,"abstract":"<div><div>In practical engineering applications, slender structures often have complex geometries, such as aircraft wings and rotor blades. Slender structures with low aspect ratios were often analyzed using three-dimensional (3D) Finite Element Analysis (FEA) as beam models were widely considered inadequate. However, analyzing slender structures with complex geometries was usually impractical due to the large number of degrees of freedom (DOFs) needed. This underscored the need to evaluate whether beam models can preserve accuracy while substantially reducing computational cost. To perform this assessment, a semi-analytical solution of the Timoshenko beam model, based on Variational Asymptotic Beam Sectional Analysis (VABS) considering all possible couplings, was used to predict the natural frequencies of general composite structures in free vibration analysis. The natural frequencies obtained from the semi-analytical solution were compared against the results from 3D FEA to assess the accuracy of the VABS-based Timoshenko model. The findings indicated that, for general composite beams, the Timoshenko model provided accurate predictions of the first several modes, excluding non-beam modes, mitigating the need for resource-intensive 3D FEA.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"328 ","pages":"Article 113852"},"PeriodicalIF":3.8,"publicationDate":"2026-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145978764","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-13DOI: 10.1016/j.ijsolstr.2026.113850
Xinyue Wu , Shengwang Hao
Strain rate effects represent a critical aspect of the constitutive relationship describing the plastic flow behaviour of a ductile material owing to the heterogeneous distributions of strain and strain rate induced by necking localisation. This study quantitatively evaluated hardening by measuring the evolutions of the local true stresses and strains in high-density polyethylene (HDPE) specimens subjected to tension at different strain loading rates. The local true stress–strain curves exhibited three typical behaviours in the plastic flow stage owing to the different distances of each section from the centre of necking. The strain and strain rate hardening coefficients were estimated using statistics describing the relationships between the local true stresses, strains, and strain rates measured at different cross-sections, where local true strains were at the same values. The strain rate hardening coefficient was constant when the local true strains were moderate and decreased logarithmically according to the applied strain loading rate. Based on these results, the strain hardening coefficient was calculated, and the constitutive equation of plastic flow behaviours was discussed. The relationship of local true stress and true strain in the necked zone was divided into four stages: the early viscoelastic stage, constant strain rate hardening coefficient stage (necking forming stage), transition stage, and necking propagation stage. Different mechanical response mechanisms host the deformation behaviour in each stage. The strain rate hardening coefficient in the second stage decreases as a logarithmic dependence on the applied (average) strain rate. The dependence of local stress on local true strain in the first and second stages follows different expressions. In the third (transition) stage, the strain rate hardening coefficient increases with the increase in local strain. Finally, these results are verified against the measured true stress and strain gradients as well as the relationship between them. The results of this study provide insights into the constitutive relationship for the whole deformation process of polymers and suggest an approach for describing mechanical behaviours through measurements of local mechanical properties.
{"title":"Revealing the strain and strain rate hardening of HDPE according to strain loading rate by measuring local stresses and strains","authors":"Xinyue Wu , Shengwang Hao","doi":"10.1016/j.ijsolstr.2026.113850","DOIUrl":"10.1016/j.ijsolstr.2026.113850","url":null,"abstract":"<div><div>Strain rate effects represent a critical aspect of the constitutive relationship describing the plastic flow behaviour of a ductile material owing to the heterogeneous distributions of strain and strain rate induced by necking localisation. This study quantitatively evaluated hardening by measuring the evolutions of the local true stresses and strains in high-density polyethylene (HDPE) specimens subjected to tension at different strain loading rates. The local true stress–strain curves exhibited three typical behaviours in the plastic flow stage owing to the different distances of each section from the centre of necking. The strain and strain rate hardening coefficients were estimated using statistics describing the relationships between the local true stresses, strains, and strain rates measured at different cross-sections, where local true strains were at the same values. The strain rate hardening coefficient was constant when the local true strains were moderate and decreased logarithmically according to the applied strain loading rate. Based on these results, the strain hardening coefficient was calculated, and the constitutive equation of plastic flow behaviours was discussed. The relationship of local true stress and true strain in the necked zone was divided into four stages: the early viscoelastic stage, constant strain rate hardening coefficient stage (necking forming stage), transition stage, and necking propagation stage. Different mechanical response mechanisms host the deformation behaviour in each stage. The strain rate hardening coefficient in the second stage decreases as a logarithmic dependence on the applied (average) strain rate. The dependence of local stress on local true strain in the first and second stages follows different expressions. In the third (transition) stage, the strain rate hardening coefficient increases with the increase in local strain. Finally, these results are verified against the measured true stress and strain gradients as well as the relationship between them. The results of this study provide insights into the constitutive relationship for the whole deformation process of polymers and suggest an approach for describing mechanical behaviours through measurements of local mechanical properties.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"328 ","pages":"Article 113850"},"PeriodicalIF":3.8,"publicationDate":"2026-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146035200","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-12DOI: 10.1016/j.ijsolstr.2026.113846
Sepideh Ebad Sichani, Xin Ning
Inspired by the continuous shells and graded porous interiors of natural bird bones, this study presents a framework to design, optimize, and additively manufacture bird-bone-like materials for a new class of aircraft wing designs without traditional components such as ribs, spars, and stiffeners. Additive manufacturing, including fused deposition modeling (FDM), enables the rapid fabrication of these complex bio-inspired geometries with minimal material waste but introduces significant anisotropy due to its layer-by-layer deposition process. We implemented a transversely isotropic material model with Hill’s yield criterion to capture the directional dependence of FDM-printed polylactic acid (PLA). Using the Covariance Matrix Adaptation Evolution Strategy (CMA-ES), the bird-bone-inspired materials were optimized to minimize wing mass while maximizing load-carrying capacity. This framework achieved substantial improvements in structural efficiency, with 48–54 % for wings with lattice-based internal structures and 23–37 % for foam-based internal structures compared to reference designs. Experimental validation through structural testing of 3D-printed wings showed strong agreement with numerical predictions, with differences in effective stiffness and load-carrying capacity within 1.4–3.3 % and 1.2–13.5 %, respectively, of simulated values. The results confirm the effectiveness of this integrated framework for designing lightweight, high-performance bird-bone-inspired materials for aerospace applications.
{"title":"Optimization, additive manufacturing, and testing of bird-bone-inspired materials for aircraft wing designs","authors":"Sepideh Ebad Sichani, Xin Ning","doi":"10.1016/j.ijsolstr.2026.113846","DOIUrl":"10.1016/j.ijsolstr.2026.113846","url":null,"abstract":"<div><div>Inspired by the continuous shells and graded porous interiors of natural bird bones, this study presents a framework to design, optimize, and additively manufacture bird-bone-like materials for a new class of aircraft wing designs without traditional components such as ribs, spars, and stiffeners. Additive manufacturing, including fused deposition modeling (FDM), enables the rapid fabrication of these complex bio-inspired geometries with minimal material waste but introduces significant anisotropy due to its layer-by-layer deposition process. We implemented a transversely isotropic material model with Hill’s yield criterion to capture the directional dependence of FDM-printed polylactic acid (PLA). Using the Covariance Matrix Adaptation Evolution Strategy (CMA-ES), the bird-bone-inspired materials were optimized to minimize wing mass while maximizing load-carrying capacity. This framework achieved substantial improvements in structural efficiency, with 48–54 % for wings with lattice-based internal structures and 23–37 % for foam-based internal structures compared to reference designs. Experimental validation through structural testing of 3D-printed wings showed strong agreement with numerical predictions, with differences in effective stiffness and load-carrying capacity within 1.4–3.3 % and 1.2–13.5 %, respectively, of simulated values. The results confirm the effectiveness of this integrated framework for designing lightweight, high-performance bird-bone-inspired materials for aerospace applications.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"328 ","pages":"Article 113846"},"PeriodicalIF":3.8,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145978765","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-12DOI: 10.1016/j.ijsolstr.2026.113849
Shah Wasif Sazzad , Sanjay Dharmavaram , Luigi E. Perotti
The ability to unfold three-dimensional curved surfaces to flat templates has many applications ranging from space exploration, to communication, to mapping and image processing. In this context, we propose a new algorithm to unfold curved surfaces to a planar template that is selected under chosen design criteria. The given surface is first tessellated and potential cutlines are identified by joining the topological defects in the tessellation. These cutlines isolate regions (or unit patches) of lower Gauss curvature, which can be flattened with smaller distortion. Vice versa, regions of high Gauss curvature are driven toward the boundaries of the template, where area can be more easily added or subtracted with the same goal of minimizing distortion. Based on the determined cutlines, a graph is constructed where the nodes correspond to the unit patches and the edges to the patches’ connectivity. The edge weights are assigned based on chosen design criteria so that the graph’s minimum spanning tree determines the connections between the unit patches in the unfolded template. In this work, we consider criteria to avoid overlapping and based on the area or shape of the unfolded template, or leading to compact refolding. Each unit patch is mapped to the flat template and linked to its adjacent subunits following the minimum spanning tree. An elastic energy minimization scheme is applied to reduce distortion. The unfolding procedure can be reversed using a separate path to achieve compact refolding, which may be advantageous for transportation and storage. The proposed strategy is demonstrated in the unfolding of icosahedral shells, geodesic domes, and a paraboloid according to different design criteria. In all these examples, the limited distortion of the unfolded template with respect to the original surface is presented.
{"title":"Using topological defects to unfold thin structures: A graph-based approach with energy-driven distortion minimization","authors":"Shah Wasif Sazzad , Sanjay Dharmavaram , Luigi E. Perotti","doi":"10.1016/j.ijsolstr.2026.113849","DOIUrl":"10.1016/j.ijsolstr.2026.113849","url":null,"abstract":"<div><div>The ability to unfold three-dimensional curved surfaces to flat templates has many applications ranging from space exploration, to communication, to mapping and image processing. In this context, we propose a new algorithm to unfold curved surfaces to a planar template that is selected under chosen design criteria. The given surface is first tessellated and potential cutlines are identified by joining the topological defects in the tessellation. These cutlines isolate regions (or unit patches) of lower Gauss curvature, which can be flattened with smaller distortion. Vice versa, regions of high Gauss curvature are driven toward the boundaries of the template, where area can be more easily added or subtracted with the same goal of minimizing distortion. Based on the determined cutlines, a graph is constructed where the nodes correspond to the unit patches and the edges to the patches’ connectivity. The edge weights are assigned based on chosen design criteria so that the graph’s minimum spanning tree determines the connections between the unit patches in the unfolded template. In this work, we consider criteria to avoid overlapping and based on the area or shape of the unfolded template, or leading to compact refolding. Each unit patch is mapped to the flat template and linked to its adjacent subunits following the minimum spanning tree. An elastic energy minimization scheme is applied to reduce distortion. The unfolding procedure can be reversed using a separate path to achieve compact refolding, which may be advantageous for transportation and storage. The proposed strategy is demonstrated in the unfolding of icosahedral shells, geodesic domes, and a paraboloid according to different design criteria. In all these examples, the limited distortion of the unfolded template with respect to the original surface is presented.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"328 ","pages":"Article 113849"},"PeriodicalIF":3.8,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145978840","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-12DOI: 10.1016/j.ijsolstr.2026.113848
Fernando Ramirez , Arturo Rodriguez-Herrera , Paul R. Heyliger
The natural vibration behavior of traction-free nanoparticles using an integral formulation of non-local elasticity theory is reported. Frequency spectra of sphere- and cubed-shaped particles composed of silicon, carbon, and germanium were computed using the finite element method. Various particle sizes, material internal lengths, and non-local weighting factors were considered. It was found that non-local frequencies are consistently lower than those obtained using local elasticity, indicating a material softening effect introduced by the non-local theory. Additionally, non-local frequencies approach those calculated using classical local elasticity as the local weighting factor increases and/or the material-internal-length to particle-size ratio decreases. Finally, the non-local frequency-radius product varies with particle size, indicating that the frequency scale invariance holding in classical elasticity is not valid. Instead, it was found that normalized non-local frequencies remain constant for a given material-internal-length to particle-size ratio, regardless of the particle size. This result introduces an alternative concept of scale invariance within the framework of non-local elasticity.
{"title":"Novel non-local frequency scale invariance for nanostructures","authors":"Fernando Ramirez , Arturo Rodriguez-Herrera , Paul R. Heyliger","doi":"10.1016/j.ijsolstr.2026.113848","DOIUrl":"10.1016/j.ijsolstr.2026.113848","url":null,"abstract":"<div><div>The natural vibration behavior of traction-free nanoparticles using an integral formulation of non-local elasticity theory is reported. Frequency spectra of sphere- and cubed-shaped particles composed of silicon, carbon, and germanium were computed using the finite element method. Various particle sizes, material internal lengths, and non-local weighting factors were considered. It was found that non-local frequencies are consistently lower than those obtained using local elasticity, indicating a material softening effect introduced by the non-local theory. Additionally, non-local frequencies approach those calculated using classical local elasticity as the local weighting factor increases and/or the material-internal-length to particle-size ratio decreases. Finally, the non-local frequency-radius product varies with particle size, indicating that the frequency scale invariance holding in classical elasticity is not valid. Instead, it was found that normalized non-local frequencies remain constant for a given material-internal-length to particle-size ratio, regardless of the particle size. This result introduces an alternative concept of scale invariance within the framework of non-local elasticity.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"328 ","pages":"Article 113848"},"PeriodicalIF":3.8,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145978762","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}