首页 > 最新文献
paperpass

Journal of The Mechanics and Physics of Solids最新文献

英文 中文
Coupled grain-pore fabric evolutions in sheared granular materials: Anisotropy lagging and geometric emergence 剪切颗粒材料中耦合颗粒-孔隙结构演化:各向异性滞后和几何涌现
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-04-01 Epub Date: 2026-01-29 DOI: 10.1016/j.jmps.2026.106530
Amiya Prakash Das , Jidong Zhao , Thomas Sweijen
The morphological evolution of pore spaces is a critical yet poorly quantified microstructural determinant of the macroscopic mechanical and hydraulic behavior of granular materials. While the anisotropy of the grain contact network (Fc) is known to dictate material response, the concurrent evolution of pore space anisotropy (Fp) and its coupling with Fc remains inadequately understood. This study employs Minkowski moment tensor analysis within a Discrete Element Method (DEM) framework to bridge this gap. We systematically investigate dense and loose, monodisperse and polydisperse assemblies under cyclic triaxial loading to quantify the dynamic coupling between Fc and Fp. We demonstrate a moderate to strong correlation between Fc and Fp, with a systematic lag in the response of Fp attributed to hierarchical geometric emergence across scales. This lag is constrained by particle-scale free-volume reorganization and its kinematic compatibility with particle motion. Additionally, key pore-scale metrics, including inverse Voronoi cell fractions (ϕv−1), pore-scale porosity (ϕp), and pore shape anisotropy β^, are well described by gamma distributions across all packing densities and strain levels. Notably, the scaled ϕv−1 follows a k-gamma distribution, providing a statistically consistent descriptor for volume fluctuations. A strong correlation is also observed between the average pore shape factor (|β|avg) and global porosity, suggesting that |β|avg serves as a geometry-based descriptor linking collective pore deformation to packing density. These findings underscore the utility of the Minkowski tensor approach in capturing 3D fabric evolution and explicitly linking pore- and grain-scale interactions. The quantitative relationships and statistical descriptors presented here provide a new foundation for enhancing constitutive models in geotechnics and powder technology, offering insights relevant to future investigations into permeability evolution and shear band formation.
孔隙空间的形态演化是颗粒材料宏观力学和水力行为的一个关键但量化不足的微观结构决定因素。虽然已知颗粒接触网络(Fc)的各向异性决定了材料的响应,但孔隙空间各向异性(Fp)的同步演化及其与Fc的耦合仍然没有得到充分的了解。本研究采用离散元法(DEM)框架内的闵可夫斯基矩张量分析来弥补这一差距。我们系统地研究了密集和松散、单分散和多分散组件在循环三轴载荷下的动态耦合。我们证明了Fc和Fp之间存在中度到强烈的相关性,而Fp的响应存在系统性滞后,这归因于不同尺度的分层几何涌现。这种滞后受到粒子尺度自由体积重组及其与粒子运动的运动学相容性的限制。此外,关键的孔隙尺度指标,包括逆Voronoi细胞分数(ϕv−1),孔隙尺度孔隙率(ϕp)和孔隙形状各向异性β^,都可以通过所有充填密度和应变水平的伽马分布来很好地描述。值得注意的是,缩放后的ϕv−1遵循k-gamma分布,为体积波动提供了统计上一致的描述符。平均孔隙形状因子(|β|avg)与整体孔隙度之间也存在很强的相关性,表明|β|avg可以作为一个基于几何的描述符,将孔隙集体变形与充填密度联系起来。这些发现强调了闵可夫斯基张量方法在捕捉三维织物演化和明确连接孔隙和颗粒尺度相互作用方面的实用性。本文提出的定量关系和统计描述符为增强岩土工程和粉末技术中的本构模型提供了新的基础,为未来研究渗透率演化和剪切带形成提供了相关见解。
{"title":"Coupled grain-pore fabric evolutions in sheared granular materials: Anisotropy lagging and geometric emergence","authors":"Amiya Prakash Das ,&nbsp;Jidong Zhao ,&nbsp;Thomas Sweijen","doi":"10.1016/j.jmps.2026.106530","DOIUrl":"10.1016/j.jmps.2026.106530","url":null,"abstract":"<div><div>The morphological evolution of pore spaces is a critical yet poorly quantified microstructural determinant of the macroscopic mechanical and hydraulic behavior of granular materials. While the anisotropy of the grain contact network (<em>F<sub>c</sub></em>) is known to dictate material response, the concurrent evolution of pore space anisotropy (<em>F<sub>p</sub></em>) and its coupling with <em>F<sub>c</sub></em> remains inadequately understood. This study employs Minkowski moment tensor analysis within a Discrete Element Method (DEM) framework to bridge this gap. We systematically investigate dense and loose, monodisperse and polydisperse assemblies under cyclic triaxial loading to quantify the dynamic coupling between <em>F<sub>c</sub></em> and <em>F<sub>p</sub></em>. We demonstrate a moderate to strong correlation between <em>F<sub>c</sub></em> and <em>F<sub>p</sub></em>, with a systematic lag in the response of <em>F<sub>p</sub></em> attributed to hierarchical geometric emergence across scales. This lag is constrained by particle-scale free-volume reorganization and its kinematic compatibility with particle motion. Additionally, key pore-scale metrics, including inverse Voronoi cell fractions <span><math><mrow><mo>(</mo><msubsup><mi>ϕ</mi><mi>v</mi><mrow><mo>−</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow></math></span>, pore-scale porosity (<em>ϕ<sub>p</sub></em>), and pore shape anisotropy <span><math><mover><mi>β</mi><mo>^</mo></mover></math></span>, are well described by gamma distributions across all packing densities and strain levels. Notably, the scaled <span><math><msubsup><mi>ϕ</mi><mi>v</mi><mrow><mo>−</mo><mn>1</mn></mrow></msubsup></math></span> follows a <em>k</em>-gamma distribution, providing a statistically consistent descriptor for volume fluctuations. A strong correlation is also observed between the average pore shape factor (|<em>β</em>|<sub><em>avg</em></sub>) and global porosity, suggesting that |<em>β</em>|<sub><em>avg</em></sub> serves as a geometry-based descriptor linking collective pore deformation to packing density. These findings underscore the utility of the Minkowski tensor approach in capturing 3D fabric evolution and explicitly linking pore- and grain-scale interactions. The quantitative relationships and statistical descriptors presented here provide a new foundation for enhancing constitutive models in geotechnics and powder technology, offering insights relevant to future investigations into permeability evolution and shear band formation.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"210 ","pages":"Article 106530"},"PeriodicalIF":6.0,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146071614","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Toggling energy landscape to enable commutative and non-commutative responses in multistable structures 在多稳定结构中切换能量格局以启用交换和非交换响应
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-04-01 Epub Date: 2026-02-03 DOI: 10.1016/j.jmps.2026.106544
Lei Wu, Damiano Pasini
The deformation of a structure may be either insensitive or sensitive to the sequence of the applied forces: the former corresponds to a commutative response, whereas the latter indicates a non-commutative response. Among diverse physical systems, an elastic multistable structure tessellated from bistable units or bits can display commutative or non-commutative responses, each dictated by the manner of interaction between the units. In this work, we elucidate the underlying physics of commutative and non-commutative responses through the lens of the energy landscape of a multistable structure, and demonstrate that initially non-commutative units can be toggled to become commutative by altering the number and distribution of local energy minima in their state space, thus enabling reprogrammable sensitivity to the loading sequence. The concept is realized and demonstrated as a multistable switch in an electrical circuit capable of offering in situ adjustable levels of information protection through reprogrammable sequence sensitivity.
结构的变形可能对施加力的顺序不敏感或敏感:前者对应于交换响应,而后者则表示非交换响应。在不同的物理系统中,由双稳态单元或位镶嵌而成的弹性多稳态结构可以显示交换或非交换响应,每个响应都由单元之间的相互作用方式决定。在这项工作中,我们通过多稳定结构的能量景观镜头阐明了交换和非交换响应的基本物理,并证明了最初的非交换单元可以通过改变其状态空间中局部能量最小值的数量和分布来切换为交换单元,从而实现对加载序列的可重新编程灵敏度。该概念被实现并演示为电路中的多稳态开关,能够通过可重新编程的序列灵敏度提供现场可调的信息保护水平。
{"title":"Toggling energy landscape to enable commutative and non-commutative responses in multistable structures","authors":"Lei Wu,&nbsp;Damiano Pasini","doi":"10.1016/j.jmps.2026.106544","DOIUrl":"10.1016/j.jmps.2026.106544","url":null,"abstract":"<div><div>The deformation of a structure may be either insensitive or sensitive to the sequence of the applied forces: the former corresponds to a commutative response, whereas the latter indicates a non-commutative response. Among diverse physical systems, an elastic multistable structure tessellated from bistable units or bits can display commutative or non-commutative responses, each dictated by the manner of interaction between the units. In this work, we elucidate the underlying physics of commutative and non-commutative responses through the lens of the energy landscape of a multistable structure, and demonstrate that initially non-commutative units can be toggled to become commutative by altering the number and distribution of local energy minima in their state space, thus enabling reprogrammable sensitivity to the loading sequence. The concept is realized and demonstrated as a multistable switch in an electrical circuit capable of offering <em>in situ</em> adjustable levels of information protection through reprogrammable sequence sensitivity.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"210 ","pages":"Article 106544"},"PeriodicalIF":6.0,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146110148","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
3D phase-field modeling of H2→H3 phase transition induced intragranular cracking in single-crystal Ni-rich layered cathodes 富镍单晶层状阴极H2→H3相变致晶内裂纹的三维相场模拟
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-04-01 Epub Date: 2026-01-23 DOI: 10.1016/j.jmps.2026.106528
Zhengtao Liao , Xiangbiao Liao , Yuyang Lu , Chengcheng Cao , Linghui He , Yong Ni
Understanding the intragranular cracking mechanisms of single-crystal Ni-rich layered cathodes (SCNCM) during electrochemical cycling is essential for next-generation high-energy-density and long-life Li-ion batteries. However, the complex interplay among the factors driving crack initiation and propagation remains unclear. Herein, we present a fully coupled three-dimensional phase-field model that integrates anisotropic lithium diffusion, the H2-H3 structural phase transition, stress evolution, and fracture mechanics to elucidate intragranular fracture in SCNCM during deep delithiation. Simulations on representative particle geometries reveal that anisotropic diffusion and fracture energy alone can not initiate layer-parallel cracks, instead, layered delithiation pathways govern shape-dependent phase transition dynamics, while lattice mismatch from spatially heterogeneous phase transitions nucleates cracks preferentially at particle surfaces and drives their propagation along layered planes, in good agreement with experimental observations. Cubic particles with uniform layers form cracks prematurely, whereas spherical and octahedral particles—with shorter surface layers—delay crack initiation by 5.3% state of charge (SoC) due to retarded phase transition in central layers, thereby expanding the safe-charging window. Elevated charging rates accelerate central-layer phase transition, amplifying misfit and triggering earlier cracking, while reduced phase-transition eigenstrain or flatter aspect ratios mitigate stress concentrations and suppress damage. These results establish a predictive link between phase transition dynamics and intragranular fracture, providing design strategies for mechanically robust high-energy-density cathodes.
了解电化学循环过程中单晶富镍层状阴极(SCNCM)的晶内开裂机制对下一代高能量密度和长寿命锂离子电池至关重要。然而,驱动裂纹萌生和扩展的因素之间复杂的相互作用仍不清楚。在此,我们提出了一个全耦合的三维相场模型,该模型集成了各向异性锂扩散、H2-H3结构相变、应力演化和断裂力学,以阐明SCNCM在深度剥蚀过程中的晶内断裂。对典型颗粒几何形状的模拟表明,各向异性扩散和断裂能本身不能引发层平行裂纹,相反,层状衰减路径控制着形状相关的相变动力学,而空间非均质相变导致的晶格失配优先在颗粒表面形成裂纹,并驱动裂纹沿层状面扩展,这与实验观察结果吻合良好。具有均匀层的立方粒子过早形成裂纹,而具有较短表面层的球形和八面体粒子由于中心层的相变延迟,延迟了5.3%的荷电状态(SoC),从而扩大了安全充电窗口。升高的充注速率加速了中间层的相变,放大了错配,引发了早期的开裂,而降低的相变特征应变或平坦的长径比则缓解了应力集中,抑制了损伤。这些结果建立了相变动力学与晶内断裂之间的预测联系,为机械坚固的高能量密度阴极的设计提供了策略。
{"title":"3D phase-field modeling of H2→H3 phase transition induced intragranular cracking in single-crystal Ni-rich layered cathodes","authors":"Zhengtao Liao ,&nbsp;Xiangbiao Liao ,&nbsp;Yuyang Lu ,&nbsp;Chengcheng Cao ,&nbsp;Linghui He ,&nbsp;Yong Ni","doi":"10.1016/j.jmps.2026.106528","DOIUrl":"10.1016/j.jmps.2026.106528","url":null,"abstract":"<div><div>Understanding the intragranular cracking mechanisms of single-crystal Ni-rich layered cathodes (SC<img>NCM) during electrochemical cycling is essential for next-generation high-energy-density and long-life Li-ion batteries. However, the complex interplay among the factors driving crack initiation and propagation remains unclear. Herein, we present a fully coupled three-dimensional phase-field model that integrates anisotropic lithium diffusion, the H2-H3 structural phase transition, stress evolution, and fracture mechanics to elucidate intragranular fracture in SC<img>NCM during deep delithiation. Simulations on representative particle geometries reveal that anisotropic diffusion and fracture energy alone can not initiate layer-parallel cracks, instead, layered delithiation pathways govern shape-dependent phase transition dynamics, while lattice mismatch from spatially heterogeneous phase transitions nucleates cracks preferentially at particle surfaces and drives their propagation along layered planes, in good agreement with experimental observations. Cubic particles with uniform layers form cracks prematurely, whereas spherical and octahedral particles—with shorter surface layers—delay crack initiation by 5.3% state of charge (SoC) due to retarded phase transition in central layers, thereby expanding the safe-charging window. Elevated charging rates accelerate central-layer phase transition, amplifying misfit and triggering earlier cracking, while reduced phase-transition eigenstrain or flatter aspect ratios mitigate stress concentrations and suppress damage. These results establish a predictive link between phase transition dynamics and intragranular fracture, providing design strategies for mechanically robust high-energy-density cathodes.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"210 ","pages":"Article 106528"},"PeriodicalIF":6.0,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146033262","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A generalized higher-order phase-field model for brittle fracture in anisotropic rocks 各向异性岩石脆性断裂的广义高阶相场模型
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-04-01 Epub Date: 2026-01-27 DOI: 10.1016/j.jmps.2026.106526
Sijia Liu , Yunteng Wang , Xueyu Geng , Wei Wu
In this work, we formulate a generalized higher-order phase-field model within the rotated anisotropic framework for simulating brittle phenomena in anisotropic rocks. Our phase-field model accounts for both the anisotropic critical fracture energy release rate and the anisotropic degradation in stiffness. The innovative aspects of this model include (i) a fourth-order structural tensor enabling simulations of strongly anisotropic fractures with arbitrary, non-orthogonal symmetry axes for capturing the complexity of natural geological media; (ii) a volumetric–deviatoric coupling energy density for transitions from anisotropic responses in the undamaged state to isotropic responses in the damaged state; (iii) a patch-based Hessian recovery algorithm ensuring stable solutions of the higher-order PDEs to reduce the computational cost; and (iv) stochastic perturbations integrated into the anisotropic crack surface density function to capture microstructural heterogeneity. Several numerical benchmark examples are provided. The numerical results are compared with some laboratory experiments on brittle fracture in anisotropic rocks.
在这项工作中,我们建立了一个在旋转各向异性框架内的广义高阶相场模型来模拟各向异性岩石中的脆性现象。我们的相场模型考虑了各向异性临界断裂能量释放率和刚度的各向异性退化。该模型的创新之处包括(i)一个四阶结构张量,可以模拟具有任意非正交对称轴的强各向异性裂缝,以捕捉自然地质介质的复杂性;(ii)从未损伤状态的各向异性响应到损伤状态的各向同性响应的体积-偏差耦合能量密度;(iii)基于patch的Hessian恢复算法,保证高阶偏微分方程解稳定,降低计算成本;(iv)将随机扰动集成到各向异性裂纹表面密度函数中,以捕获微观结构的非均质性。给出了几个数值基准算例。数值结果与一些各向异性岩石脆性断裂的室内实验结果进行了比较。
{"title":"A generalized higher-order phase-field model for brittle fracture in anisotropic rocks","authors":"Sijia Liu ,&nbsp;Yunteng Wang ,&nbsp;Xueyu Geng ,&nbsp;Wei Wu","doi":"10.1016/j.jmps.2026.106526","DOIUrl":"10.1016/j.jmps.2026.106526","url":null,"abstract":"<div><div>In this work, we formulate a generalized higher-order phase-field model within the rotated anisotropic framework for simulating brittle phenomena in anisotropic rocks. Our phase-field model accounts for both the anisotropic critical fracture energy release rate and the anisotropic degradation in stiffness. The innovative aspects of this model include (i) a fourth-order structural tensor enabling simulations of strongly anisotropic fractures with arbitrary, non-orthogonal symmetry axes for capturing the complexity of natural geological media; (ii) a volumetric–deviatoric coupling energy density for transitions from anisotropic responses in the undamaged state to isotropic responses in the damaged state; (iii) a patch-based Hessian recovery algorithm ensuring stable solutions of the higher-order PDEs to reduce the computational cost; and (iv) stochastic perturbations integrated into the anisotropic crack surface density function to capture microstructural heterogeneity. Several numerical benchmark examples are provided. The numerical results are compared with some laboratory experiments on brittle fracture in anisotropic rocks.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"210 ","pages":"Article 106526"},"PeriodicalIF":6.0,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146056280","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
In search of constitutive conditions in isotropic hyperelasticity: polyconvexity versus true-stress-true-strain monotonicity 寻找各向同性超弹性的本构条件:多凸性与真应力-真应变单调性
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-03-01 Epub Date: 2025-12-04 DOI: 10.1016/j.jmps.2025.106465
Maximilian P. Wollner , Gerhard A. Holzapfel , Patrizio Neff
The polyconvexity of a strain-energy function is nowadays increasingly presented as the ultimate material stability condition for an idealized elastic response. While the mathematical merits of polyconvexity are clearly understood, its mechanical consequences have received less attention. In this contribution we contrast polyconvexity with the recently rediscovered true-stress-true-strain monotonicity (TSTS-M++) condition. By way of explicit examples, we show that neither condition by itself is strong enough to guarantee physically reasonable behavior for ideal isotropic elasticity. In particular, polyconvexity does not imply a monotone trajectory of the Cauchy stress in unconstrained uniaxial extension which TSTS-M++ ensures. On the other hand, TSTS-M++ does not impose a monotone Cauchy shear stress response in simple shear which is enforced by Legendre-Hadamard ellipticity and in turn polyconvexity. Both scenarios are proven through the construction of appropriate strain-energy functions. Consequently, a combination of polyconvexity, ensuring Legendre-Hadamard ellipticity, and TSTS-M++ seems to be a viable solution to Truesdell’s Hauptproblem. However, so far no isotropic strain-energy function has been identified that satisfies both constraints globally at the same time. Although we are unable to deliver a valid solution here, we provide several results that could prove helpful in the construction of such an exceptional strain-energy function. In memory of Miroslav Šilhavý (1949–2025)
应变-能量函数的多凸性现在越来越多地被认为是理想弹性响应的最终材料稳定性条件。虽然多凸性的数学优点已被清楚地理解,但它的力学后果却很少受到关注。在这篇贡献中,我们将多凸性与最近重新发现的真应力-真应变单调性(tsts - m++)条件进行了比较。通过明确的例子,我们表明,这两个条件本身都不足以保证理想各向同性弹性的物理合理行为。特别是,多凸性并不意味着柯西应力在无约束单轴拉伸中的单调轨迹,而tsts - m++保证了这一点。另一方面,tsts - m++在简单剪切中不施加单调的柯西剪应力响应,而单调的柯西剪应力响应是由legende - hadamard椭圆性和多凸性所施加的。通过构造适当的应变-能函数来证明这两种情况。因此,确保legende - hadamard椭圆性的多凸性与tsts - m++的组合似乎是解决Truesdell Hauptproblem的可行方案。然而,到目前为止,还没有发现同时满足这两个约束的各向同性应变能函数。虽然我们无法在这里提供一个有效的解决方案,但我们提供了几个结果,可以证明有助于构建这样一个特殊的应变-能量函数。纪念米罗斯拉夫Šilhavý (1949-2025)
{"title":"In search of constitutive conditions in isotropic hyperelasticity: polyconvexity versus true-stress-true-strain monotonicity","authors":"Maximilian P. Wollner ,&nbsp;Gerhard A. Holzapfel ,&nbsp;Patrizio Neff","doi":"10.1016/j.jmps.2025.106465","DOIUrl":"10.1016/j.jmps.2025.106465","url":null,"abstract":"<div><div>The polyconvexity of a strain-energy function is nowadays increasingly presented as the ultimate material stability condition for an idealized elastic response. While the mathematical merits of polyconvexity are clearly understood, its mechanical consequences have received less attention. In this contribution we contrast polyconvexity with the recently rediscovered true-stress-true-strain monotonicity (TSTS-M<sup>++</sup>) condition. By way of explicit examples, we show that neither condition by itself is strong enough to guarantee physically reasonable behavior for ideal isotropic elasticity. In particular, polyconvexity does not imply a monotone trajectory of the Cauchy stress in unconstrained uniaxial extension which TSTS-M<sup>++</sup> ensures. On the other hand, TSTS-M<sup>++</sup> does not impose a monotone Cauchy shear stress response in simple shear which is enforced by Legendre-Hadamard ellipticity and in turn polyconvexity. Both scenarios are proven through the construction of appropriate strain-energy functions. Consequently, a combination of polyconvexity, ensuring Legendre-Hadamard ellipticity, and TSTS-M<sup>++</sup> seems to be a viable solution to Truesdell’s Hauptproblem. However, so far no isotropic strain-energy function has been identified that satisfies both constraints globally at the same time. Although we are unable to deliver a valid solution here, we provide several results that could prove helpful in the construction of such an exceptional strain-energy function. <em>In memory of</em> Miroslav Šilhavý (1949–2025)</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"209 ","pages":"Article 106465"},"PeriodicalIF":6.0,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145689708","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Homogenization framework for rigid and non-rigid foldable origami metamaterials 刚性和非刚性可折叠折纸超材料的均匀化框架
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-03-01 Epub Date: 2026-01-16 DOI: 10.1016/j.jmps.2026.106519
Xuwen Li, Amin Jamalimehr, Mathias Legrand, Damiano Pasini
Origami metamaterials typically consist of folded sheets with periodic patterns, conferring them with remarkable mechanical properties. In the context of Continuum Mechanics, the majority of existing predictive methods are mechanism analogs which favor rigid folding and panel bending. While effective in predicting primary deformation modes, existing methods fall short in capturing the full spectrum of deformation of non-rigid foldable origami, such as the emergence of curvature along straight creases, local strain at vertices and warpage in panels. To fully capture the entire deformation spectrum and enhance the accuracy of existing methods, this paper introduces a homogenization framework for origami metamaterials where the faces are modeled as plate elements. Both asymptotic and energy-based homogenization methods are formulated and implemented. As a representative crease pattern, we examine the Miura origami sheet homogenized as an equivalent Kirchhoff-Love plate. The results reveal that certain effective elastic properties are nonlinearly related to both the initial fold angle and the crease stiffness. When benchmarked with results from fully resolved simulations, our framework yields errors up to 12.9%, while existing models, including the bar-and-hinge model and the rigid-panel model, show up to 161% error. The differences in errors are associated with the complex modes of crease and panel deformation in non-rigid origami, unexplored by the existing models. This work demonstrates a precise and efficient continuum framework for origami metamaterials as an effective strategy for predicting their elastic properties, understanding their mechanics, and designing their functionalities.
折纸超材料通常由具有周期性图案的折叠片组成,赋予它们非凡的机械性能。在连续介质力学的背景下,大多数现有的预测方法是机制类比,有利于刚性折叠和板弯曲。虽然现有的方法可以有效地预测原始变形模式,但在捕捉非刚性可折叠折纸的全变形谱方面存在不足,例如沿直折痕出现的曲率、顶点处的局部应变和面板中的翘曲。为了充分捕获整个变形谱,提高现有方法的精度,本文引入了一种折纸超材料的均匀化框架,其中面建模为板单元。提出并实现了渐近均匀化方法和基于能量的均匀化方法。作为一个代表性的折痕图案,我们研究了三浦折纸片均质为等效的基希霍夫-洛夫板。结果表明,某些有效弹性性能与初始折叠角和折痕刚度均呈非线性关系。当与完全解析的模拟结果进行基准测试时,我们的框架产生的误差高达12.9%,而现有模型,包括杆铰模型和刚性面板模型,显示的误差高达161%。误差的差异与非刚性折纸中折痕和面板变形的复杂模式有关,这是现有模型未探索的。这项工作为折纸超材料提供了一个精确而高效的连续体框架,作为预测其弹性特性、理解其力学和设计其功能的有效策略。
{"title":"Homogenization framework for rigid and non-rigid foldable origami metamaterials","authors":"Xuwen Li,&nbsp;Amin Jamalimehr,&nbsp;Mathias Legrand,&nbsp;Damiano Pasini","doi":"10.1016/j.jmps.2026.106519","DOIUrl":"10.1016/j.jmps.2026.106519","url":null,"abstract":"<div><div>Origami metamaterials typically consist of folded sheets with periodic patterns, conferring them with remarkable mechanical properties. In the context of Continuum Mechanics, the majority of existing predictive methods are mechanism analogs which favor rigid folding and panel bending. While effective in predicting primary deformation modes, existing methods fall short in capturing the full spectrum of deformation of non-rigid foldable origami, such as the emergence of curvature along straight creases, local strain at vertices and warpage in panels. To fully capture the entire deformation spectrum and enhance the accuracy of existing methods, this paper introduces a homogenization framework for origami metamaterials where the faces are modeled as plate elements. Both asymptotic and energy-based homogenization methods are formulated and implemented. As a representative crease pattern, we examine the Miura origami sheet homogenized as an equivalent Kirchhoff-Love plate. The results reveal that certain effective elastic properties are nonlinearly related to both the initial fold angle and the crease stiffness. When benchmarked with results from fully resolved simulations, our framework yields errors up to 12.9%, while existing models, including the bar-and-hinge model and the rigid-panel model, show up to 161% error. The differences in errors are associated with the complex modes of crease and panel deformation in non-rigid origami, unexplored by the existing models. This work demonstrates a precise and efficient continuum framework for origami metamaterials as an effective strategy for predicting their elastic properties, understanding their mechanics, and designing their functionalities.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"209 ","pages":"Article 106519"},"PeriodicalIF":6.0,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145995466","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
The effect of fiber plasticity on domain formation in soft biological composites—Part I: A bifurcation analysis 纤维塑性对软生物复合材料结构域形成的影响——第一部分:分岔分析
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-03-01 Epub Date: 2025-12-25 DOI: 10.1016/j.jmps.2025.106482
Michalis Agoras , Fernanda F. Fontenele , Nikolaos Bouklas
The main objective of this work is to shed light on the effect of fiber plasticity on the macroscopic response and domain formation in soft biological composites. This goal is pursued by analyzing the plane-strain response of two-phase laminates. In the context of this problem, the effect of fiber plasticity is accounted for by allowing the elastically stiffer layers (“fiber” phase) to also yield plastically and by taking the soft layers (“matrix” phase) to be purely elastic solids. The analysis is carried out at finite elastic and plastic strains, but it is restricted to unidirectional, non-monotonic loading paths, applied by initially increasing the macroscopic stretch along the direction of the layers up to a prescribed maximum value and then decreasing the same stretch down to a minimum value. A simple expression is derived for the critical conditions at which the homogenized behavior of the laminate loses strong ellipticity for the first time along the loading path. The relevance of this result stems from the fact that the loss of macroscopic ellipticity of these composites is known to coincide with the onset of bifurcations of the long-wavelength type. It follows from this result that, just like hyperelastic laminates, elastoplastic laminates may lose macroscopic ellipticity whenever their incremental strength in shear perpendicular to the layers vanishes for the first time. It is shown by means of specific numerical calculations that fiber plasticity has a softening effect on the critical shear strength of the laminate under monotonically increasing loading and a hardening effect under monotonically decreasing loading, whereas local elasticity has the opposite effects. Thus, the effects of local elasticity and fiber plasticity compete with each other at every stage of the deformation where both mechanisms are active and whether or not a specific laminate will lose its macroscopic ellipticity along a given loading path depends crucially on the relative strengths of these effects. In this regard, the influence of the relevant loading and material parameters is investigated in detail. For situations in which loss of macroscopic ellipticity does take place, a corresponding post-bifurcation solution for the homogenized behavior of the laminate is computed. The deformed state of the material described by this solution is characterized by twin lamellar domains that are formed at a length scale much larger than the width of the original, microscopic layers, but still much smaller than the overall dimensions of the macroscopic specimen under consideration. The macroscopic response of these composites is found to be much softer along the bifurcated equilibrium path than along the corresponding principal path. Nevertheless, no macroscopic unloading has been observed in the post-bifurcation regime, indicating accordingly that the macroscopic, post-bifurcation behavior of these materials is stable.
本研究的主要目的是揭示纤维塑性对软质生物复合材料宏观响应和结构域形成的影响。通过分析两相层合板的平面应变响应来实现这一目标。在这个问题的背景下,纤维塑性的影响是通过允许弹性较硬的层(“纤维”相)也具有塑性屈服,并通过将软层(“基体”相)作为纯弹性固体来解释的。该分析是在有限的弹塑性应变下进行的,但它仅限于单向、非单调的加载路径,通过首先沿层向增加宏观拉伸到规定的最大值,然后减少相同的拉伸到最小值来施加。导出了层合板沿加载路径均匀性首次失去强椭圆性临界条件的简单表达式。这一结果的相关性源于这样一个事实,即已知这些复合材料的宏观椭圆性的损失与长波长型分岔的开始相吻合。由此可知,弹塑性层合板与超弹性层合板一样,在垂直于层合板的剪切强度增量首次消失时,也会失去宏观椭圆性。通过具体的数值计算表明,在单调递增加载条件下,纤维塑性对层合板的临界抗剪强度有软化作用,在单调递减加载条件下,纤维塑性对层合板的临界抗剪强度有硬化作用,而局部弹性则相反。因此,局部弹性和纤维塑性的影响在变形的每个阶段都是相互竞争的,在这两个机制都是活跃的,一个特定的层合板是否会沿着给定的加载路径失去宏观椭圆性,关键取决于这些影响的相对强度。为此,详细研究了相关载荷和材料参数的影响。对于确实发生宏观椭圆度损失的情况,计算了相应的层合板均匀化行为的分岔后解。该溶液描述的材料变形状态的特征是双片层畴,其形成的长度尺度远大于原始微观层的宽度,但仍远小于所考虑的宏观试样的总体尺寸。这些复合材料的宏观响应沿分岔平衡路径比沿相应的主路径要软得多。然而,在分岔后没有观察到宏观卸载,这表明这些材料的宏观分岔后行为是稳定的。
{"title":"The effect of fiber plasticity on domain formation in soft biological composites—Part I: A bifurcation analysis","authors":"Michalis Agoras ,&nbsp;Fernanda F. Fontenele ,&nbsp;Nikolaos Bouklas","doi":"10.1016/j.jmps.2025.106482","DOIUrl":"10.1016/j.jmps.2025.106482","url":null,"abstract":"<div><div>The main objective of this work is to shed light on the effect of fiber plasticity on the macroscopic response and domain formation in soft biological composites. This goal is pursued by analyzing the plane-strain response of two-phase laminates. In the context of this problem, the effect of fiber plasticity is accounted for by allowing the elastically stiffer layers (“fiber” phase) to also yield plastically and by taking the soft layers (“matrix” phase) to be purely elastic solids. The analysis is carried out at finite elastic and plastic strains, but it is restricted to unidirectional, non-monotonic loading paths, applied by initially increasing the macroscopic stretch along the direction of the layers up to a prescribed maximum value and then decreasing the same stretch down to a minimum value. A simple expression is derived for the critical conditions at which the homogenized behavior of the laminate loses strong ellipticity for the first time along the loading path. The relevance of this result stems from the fact that the loss of macroscopic ellipticity of these composites is known to coincide with the onset of bifurcations of the long-wavelength type. It follows from this result that, just like hyperelastic laminates, elastoplastic laminates may lose macroscopic ellipticity whenever their incremental strength in shear perpendicular to the layers vanishes for the first time. It is shown by means of specific numerical calculations that fiber plasticity has a softening effect on the critical shear strength of the laminate under monotonically increasing loading and a hardening effect under monotonically decreasing loading, whereas local elasticity has the opposite effects. Thus, the effects of local elasticity and fiber plasticity compete with each other at every stage of the deformation where both mechanisms are active and whether or not a specific laminate will lose its macroscopic ellipticity along a given loading path depends crucially on the relative strengths of these effects. In this regard, the influence of the relevant loading and material parameters is investigated in detail. For situations in which loss of macroscopic ellipticity does take place, a corresponding post-bifurcation solution for the homogenized behavior of the laminate is computed. The deformed state of the material described by this solution is characterized by twin lamellar domains that are formed at a length scale much larger than the width of the original, microscopic layers, but still much smaller than the overall dimensions of the macroscopic specimen under consideration. The macroscopic response of these composites is found to be much softer along the bifurcated equilibrium path than along the corresponding principal path. Nevertheless, no macroscopic unloading has been observed in the post-bifurcation regime, indicating accordingly that the macroscopic, post-bifurcation behavior of these materials is stable.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"209 ","pages":"Article 106482"},"PeriodicalIF":6.0,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145845102","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Shear banding in high-entropy alloy 高熵合金中的剪切带
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-03-01 Epub Date: 2026-01-12 DOI: 10.1016/j.jmps.2026.106517
De-Ru Wang , Tong Li , Ming-Yao Su , Yan Chen , Hai-Ying Wang , Lan-Hong Dai
Shear banding is widely regarded as a precursor to dynamic failure of solids, yet what drives its emergence remains elusive. Here, we precisely capture the onset of shear bands in the most typical Cantor CoCrFeMnNi high-entropy alloy via the unique single-pulse Hopkinson torsion bar-based deformation-freezing hybridizing synchronized temperature measurement technique. It is surprisingly found that shear bands emerge well before the stress-strain curve peak and exhibit a distinctive nanoscale morphology, in sharp contrast to traditional consensus. Precise microstructural characterization reveals that these bands originate from discrete stacking faults, where atomic-scale dilatation at their tips reduces shear deformation resistance and triggers configurational softening. To disentangle and quantify the competing effect between thermal softening and configurational softening, we develop a non-equilibrium two-temperature continuum plastic flow model. Stability analysis establishes the nanoscale shear band instability criterion and demonstrates that configurational softening, rather than widely accepted adiabatic thermal softening, governs their emergence. These findings shed new insights into understanding shear banding and dynamic failure behaviors in metals.
剪切带被广泛认为是固体动态破坏的前兆,但导致其出现的原因仍然难以捉摸。在这里,我们通过独特的单脉冲霍普金森扭转杆变形-冻结杂交同步温度测量技术,精确捕捉了最典型的Cantor CoCrFeMnNi高熵合金中剪切带的开始。令人惊讶的是,剪切带出现在应力-应变曲线峰值之前,并表现出独特的纳米级形态,与传统的共识形成鲜明对比。精确的微观结构表征表明,这些条带起源于离散的堆积断层,其尖端的原子尺度膨胀降低了剪切变形阻力并触发了构型软化。为了解开和量化热软化和构型软化之间的竞争效应,我们建立了一个非平衡双温度连续体塑性流动模型。稳定性分析建立了纳米尺度剪切带失稳判据,并证明了它们的出现是由构型软化而不是广泛接受的绝热软化控制的。这些发现为理解金属的剪切带和动态破坏行为提供了新的见解。
{"title":"Shear banding in high-entropy alloy","authors":"De-Ru Wang ,&nbsp;Tong Li ,&nbsp;Ming-Yao Su ,&nbsp;Yan Chen ,&nbsp;Hai-Ying Wang ,&nbsp;Lan-Hong Dai","doi":"10.1016/j.jmps.2026.106517","DOIUrl":"10.1016/j.jmps.2026.106517","url":null,"abstract":"<div><div>Shear banding is widely regarded as a precursor to dynamic failure of solids, yet what drives its emergence remains elusive. Here, we precisely capture the onset of shear bands in the most typical Cantor CoCrFeMnNi high-entropy alloy via the unique single-pulse Hopkinson torsion bar-based deformation-freezing hybridizing synchronized temperature measurement technique. It is surprisingly found that shear bands emerge well before the stress-strain curve peak and exhibit a distinctive nanoscale morphology, in sharp contrast to traditional consensus. Precise microstructural characterization reveals that these bands originate from discrete stacking faults, where atomic-scale dilatation at their tips reduces shear deformation resistance and triggers configurational softening. To disentangle and quantify the competing effect between thermal softening and configurational softening, we develop a non-equilibrium two-temperature continuum plastic flow model. Stability analysis establishes the nanoscale shear band instability criterion and demonstrates that configurational softening, rather than widely accepted adiabatic thermal softening, governs their emergence. These findings shed new insights into understanding shear banding and dynamic failure behaviors in metals.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"209 ","pages":"Article 106517"},"PeriodicalIF":6.0,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956824","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A unified multi-perspective quadratic manifold for mitigating the Kolmogorov barrier in multiphysics damage 用于减轻多物理场损伤中Kolmogorov势垒的统一多视角二次流形
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-03-01 Epub Date: 2026-01-01 DOI: 10.1016/j.jmps.2025.106499
Qinghua Zhang , Stephan Ritzert , Jian Zhang , Jannick Kehls , Stefanie Reese , Tim Brepols
In multiphysics damage problems, material degradation is often modeled using local and/or global damage variables, whose evolution introduces strong nonlinearities and significant computational costs. Linear projection-based reduced-order models (ROMs) are widely used to accelerate these simulations but often fail to capture complex nonlinear damage evolution effectively. This limitation arises from the slow decay of the Kolmogorov n-width, which leads to a phenomenon known as the Kolmogorov barrier in linear approximation. To overcome this challenge, this study proposes a novel unified multi-perspective (multi-field and multi-state) quadratic manifold-based ROM framework for thermo-mechanically coupled damage-plasticity problems. A key feature lies in a multi-field and multi-state decomposition strategy that is grounded in the material’s physical response to guide the selection of mode numbers for each coupled field. Moreover, the framework decouples both material states and physical fields, providing clearer insights into the contributions and interactions of each field within the overall multiphysics simulation. Benchmark tests demonstrate that the proposed approach mitigates the Kolmogorov barrier of linear projection-based ROMs by ensuring a smooth and monotonic decrease in error as the number of modes increases. The proposed multi-perspective quadratic manifold framework offers a robust and flexible approach for efficiently reducing complex damage-involved multiphysics problems and shows the potential for industrial applications.
在多物理场损伤问题中,材料退化通常使用局部和/或全局损伤变量来建模,这些变量的演化引入了强非线性和显著的计算成本。基于线性投影的降阶模型(ROMs)被广泛用于加速这些模拟,但往往不能有效地捕捉复杂的非线性损伤演化。这种限制源于柯尔莫哥洛夫n-宽度的缓慢衰减,这导致线性近似中的柯尔莫哥洛夫势垒现象。为了克服这一挑战,本研究提出了一种新的统一的多视角(多场和多态)二次流形ROM框架,用于热-机械耦合损伤-塑性问题。其关键特征在于基于材料物理响应的多场多态分解策略,以指导每个耦合场的模态数选择。此外,该框架解耦了物质状态和物理场,为整个多物理场模拟中每个场的贡献和相互作用提供了更清晰的见解。基准测试表明,该方法通过确保误差随模态数量的增加而平滑单调地减小,从而减轻了基于线性投影的rom的Kolmogorov障碍。所提出的多视角二次流形框架为有效减少涉及复杂损伤的多物理场问题提供了一种鲁棒和灵活的方法,并显示出工业应用的潜力。
{"title":"A unified multi-perspective quadratic manifold for mitigating the Kolmogorov barrier in multiphysics damage","authors":"Qinghua Zhang ,&nbsp;Stephan Ritzert ,&nbsp;Jian Zhang ,&nbsp;Jannick Kehls ,&nbsp;Stefanie Reese ,&nbsp;Tim Brepols","doi":"10.1016/j.jmps.2025.106499","DOIUrl":"10.1016/j.jmps.2025.106499","url":null,"abstract":"<div><div>In multiphysics damage problems, material degradation is often modeled using local and/or global damage variables, whose evolution introduces strong nonlinearities and significant computational costs. Linear projection-based reduced-order models (ROMs) are widely used to accelerate these simulations but often fail to capture complex nonlinear damage evolution effectively. This limitation arises from the slow decay of the Kolmogorov <em>n</em>-width, which leads to a phenomenon known as the Kolmogorov barrier in linear approximation. To overcome this challenge, this study proposes a novel unified multi-perspective (multi-field and multi-state) quadratic manifold-based ROM framework for thermo-mechanically coupled damage-plasticity problems. A key feature lies in a multi-field and multi-state decomposition strategy that is grounded in the material’s physical response to guide the selection of mode numbers for each coupled field. Moreover, the framework decouples both material states and physical fields, providing clearer insights into the contributions and interactions of each field within the overall multiphysics simulation. Benchmark tests demonstrate that the proposed approach mitigates the Kolmogorov barrier of linear projection-based ROMs by ensuring a smooth and monotonic decrease in error as the number of modes increases. The proposed multi-perspective quadratic manifold framework offers a robust and flexible approach for efficiently reducing complex damage-involved multiphysics problems and shows the potential for industrial applications.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"209 ","pages":"Article 106499"},"PeriodicalIF":6.0,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145894534","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Electro-mechanical wrinkling of soft dielectric films bonded to hyperelastic substrates 超弹性衬底上软介质薄膜的机电起皱
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-03-01 Epub Date: 2025-12-30 DOI: 10.1016/j.jmps.2025.106490
Bin Wu , Linghao Kong , Weiqiu Chen , Davide Riccobelli , Michel Destrade
Active control of wrinkling in soft film-substrate composites using electric fields is a critical challenge in tunable material systems. Here, we investigate the electro-mechanical instability of a soft dielectric film bonded to a hyperelastic substrate, revealing the fundamental mechanisms that enable on-demand surface patterning. For the linearized stability analysis, we use the Stroh formalism and the surface impedance method to obtain exact and sixth-order approximate bifurcation equations that signal the onset of wrinkles. We derive the explicit bifurcation equations giving the critical stretch and critical voltage for wrinkling, as well as the corresponding critical wavenumber. We look at scenarios where the voltage is kept constant and the stretch changes, and vice versa. We provide the thresholds of the shear modulus ratio rc0 or pre-stretch λc0 below which the film-substrate system wrinkles mechanically, prior to the application of a voltage. These predictions offer theoretical guidance for practical structural design, as the shear modulus ratio r and/or the pre-stretch λ can be chosen to be slightly greater than rc0 and/or λc0, so that the film-substrate system wrinkles with a small applied voltage. Finally, we simulate the full nonlinear behavior using the Finite Element method (FEniCS) to validate our formulas and conduct a post-buckling analysis. This work advances the fundamental understanding of electro-mechanical wrinkling instabilities in soft material systems. By enabling active control of surface morphologies via applied electric fields, our findings open new avenues for adaptive technologies in soft robotics, flexible electronics, smart surfaces, and bioinspired systems.
利用电场对软膜基复合材料的起皱进行主动控制是可调材料系统中的一个关键挑战。在这里,我们研究了与超弹性衬底结合的软介质薄膜的机电不稳定性,揭示了实现按需表面图像化的基本机制。对于线性化的稳定性分析,我们使用Stroh形式和表面阻抗法获得精确的六阶近似分岔方程,这些分岔方程表示皱纹的开始。我们导出了显式分岔方程,给出了起皱的临界拉伸和临界电压,以及相应的临界波数。我们观察电压保持不变而拉伸变化的情况,反之亦然。我们提供了剪切模量比r0或预拉伸λc0的阈值,低于该阈值,薄膜-衬底系统在施加电压之前机械起皱。这些预测为实际结构设计提供了理论指导,因为剪切模量比r和/或预拉伸λ可以选择略大于r0和/或λc0,从而使薄膜-衬底系统在小的施加电压下起皱。最后,我们使用有限元方法(FEniCS)模拟了完整的非线性行为,以验证我们的公式并进行了屈曲后分析。这项工作促进了对软材料系统中机电起皱不稳定性的基本理解。通过应用电场实现表面形态的主动控制,我们的发现为软机器人、柔性电子、智能表面和仿生系统的自适应技术开辟了新的途径。
{"title":"Electro-mechanical wrinkling of soft dielectric films bonded to hyperelastic substrates","authors":"Bin Wu ,&nbsp;Linghao Kong ,&nbsp;Weiqiu Chen ,&nbsp;Davide Riccobelli ,&nbsp;Michel Destrade","doi":"10.1016/j.jmps.2025.106490","DOIUrl":"10.1016/j.jmps.2025.106490","url":null,"abstract":"<div><div>Active control of wrinkling in soft film-substrate composites using electric fields is a critical challenge in tunable material systems. Here, we investigate the electro-mechanical instability of a soft dielectric film bonded to a hyperelastic substrate, revealing the fundamental mechanisms that enable on-demand surface patterning. For the linearized stability analysis, we use the Stroh formalism and the surface impedance method to obtain exact and sixth-order approximate bifurcation equations that signal the onset of wrinkles. We derive the explicit bifurcation equations giving the critical stretch and critical voltage for wrinkling, as well as the corresponding critical wavenumber. We look at scenarios where the voltage is kept constant and the stretch changes, and vice versa. We provide the thresholds of the shear modulus ratio <span><math><msubsup><mi>r</mi><mrow><mrow><mi>c</mi></mrow></mrow><mn>0</mn></msubsup></math></span> or pre-stretch <span><math><msubsup><mi>λ</mi><mrow><mrow><mi>c</mi></mrow></mrow><mn>0</mn></msubsup></math></span> below which the film-substrate system wrinkles mechanically, prior to the application of a voltage. These predictions offer theoretical guidance for practical structural design, as the shear modulus ratio <em>r</em> and/or the pre-stretch <em>λ</em> can be chosen to be slightly greater than <span><math><msubsup><mi>r</mi><mrow><mrow><mi>c</mi></mrow></mrow><mn>0</mn></msubsup></math></span> and/or <span><math><msubsup><mi>λ</mi><mrow><mrow><mi>c</mi></mrow></mrow><mn>0</mn></msubsup></math></span>, so that the film-substrate system wrinkles with a small applied voltage. Finally, we simulate the full nonlinear behavior using the Finite Element method (<span>FEniCS</span>) to validate our formulas and conduct a post-buckling analysis. This work advances the fundamental understanding of electro-mechanical wrinkling instabilities in soft material systems. By enabling active control of surface morphologies via applied electric fields, our findings open new avenues for adaptive technologies in soft robotics, flexible electronics, smart surfaces, and bioinspired systems.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"209 ","pages":"Article 106490"},"PeriodicalIF":6.0,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145894540","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
期刊
Journal of The Mechanics and Physics of Solids
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
﹀
mjeditor
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1