首页 > 最新文献
paperpass

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

英文 中文
Unveiling the evolving fracture toughness: The R-curve phenomenon in double network hydrogels 揭示断裂韧性的演化:双网状水凝胶中的r曲线现象
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-04-01 Epub Date: 2026-01-27 DOI: 10.1016/j.jmps.2026.106531
Wenjing Lu , Shuai Xu , Chong Wang , Zishun Liu
Double network (DN) hydrogels exhibit remarkable mechanical properties, most notably their high fracture toughness, which expands their potential applications across various fields. To guarantee the reliability of DN gels in practical uses, an in-depth comprehension of their fracture characteristics is crucial. Previous studies on the fracture behavior of DN gels mainly concentrated on the crack initiation process, often treating fracture toughness as a constant characteristic of the material. In contrast, our experimental research unveils the evolving fracture toughness as a fracture resistance curve (R-curve) in DN gels, indicating that fracture toughness increases with crack propagation until reaching a plateau of steady-state, rather than remaining constant. This phenomenon is attributed to mechanisms of crack tip softening and stress de-concentration resulting from the fracture of the brittle PAMPS network in DN gels. We have identified three critical parameters, i.e., the initiation fracture toughness (Ginit), steady state toughness (Gss), steady state crack extension length (Lss), to quantify this behavior. Additionally, we explore the effect of pre-damage of PAMPS network on the R-curve and develop a theoretical model linking the degree of pre-damage (h) to Gss. Our model effectively predicts the R-curve of DN gels under various pre-damage conditions. Furthermore, by considering different PAMPS network formulations, we establish a scaling law linking (Gss−Ginit) with Lss, thus creating a unified framework for understanding the internal physical properties of diverse DN gels. These insights also differentiate the toughening mechanisms in DN gels: energy dissipation from the PAMPS network affects (Gss−Ginit), while the Ginit reflects the intrinsic fracture properties governed by the elastic PAAm network. We anticipate that our experimental results and theoretical models will enhance the understanding of fracture in soft tough materials.
双网(DN)水凝胶具有优异的力学性能,特别是其高断裂韧性,这扩大了其在各个领域的潜在应用。为了保证DN凝胶在实际应用中的可靠性,深入了解其断裂特性至关重要。以往对DN凝胶断裂行为的研究主要集中在裂纹萌生过程上,往往将断裂韧性作为材料的恒定特征。相比之下,我们的实验研究揭示了断裂韧性在DN凝胶中的演变为断裂阻力曲线(r曲线),表明断裂韧性随着裂纹的扩展而增加,直到达到稳态状态的平台,而不是保持不变。这一现象归因于DN凝胶中脆性PAMPS网络断裂导致的裂纹尖端软化和应力去集中机制。我们确定了三个关键参数,即起裂韧性(Ginit),稳态韧性(Gss),稳态裂纹扩展长度(Lss),以量化这种行为。此外,我们还探讨了PAMPS网络预损伤对r曲线的影响,并建立了预损伤程度(h)与Gss之间的理论模型。我们的模型有效地预测了DN凝胶在各种预损伤条件下的r曲线。此外,通过考虑不同的PAMPS网络配方,我们建立了连接(Gss−Ginit)与Lss的标度定律,从而为理解不同DN凝胶的内部物理性质创建了统一的框架。这些见解也区分了DN凝胶中的增韧机制:PAMPS网络的能量耗散影响(Gss - Ginit),而Ginit反映了弹性PAAm网络控制的固有断裂特性。我们期望我们的实验结果和理论模型能增进对软韧材料断裂的认识。
{"title":"Unveiling the evolving fracture toughness: The R-curve phenomenon in double network hydrogels","authors":"Wenjing Lu ,&nbsp;Shuai Xu ,&nbsp;Chong Wang ,&nbsp;Zishun Liu","doi":"10.1016/j.jmps.2026.106531","DOIUrl":"10.1016/j.jmps.2026.106531","url":null,"abstract":"<div><div>Double network (DN) hydrogels exhibit remarkable mechanical properties, most notably their high fracture toughness, which expands their potential applications across various fields. To guarantee the reliability of DN gels in practical uses, an in-depth comprehension of their fracture characteristics is crucial. Previous studies on the fracture behavior of DN gels mainly concentrated on the crack initiation process, often treating fracture toughness as a constant characteristic of the material. In contrast, our experimental research unveils the evolving fracture toughness as a fracture resistance curve (R-curve) in DN gels, indicating that fracture toughness increases with crack propagation until reaching a plateau of steady-state, rather than remaining constant. This phenomenon is attributed to mechanisms of crack tip softening and stress de-concentration resulting from the fracture of the brittle PAMPS network in DN gels. We have identified three critical parameters, i.e., the initiation fracture toughness (<span><math><msub><mi>G</mi><mtext>init</mtext></msub></math></span>), steady state toughness (<span><math><msub><mi>G</mi><mtext>ss</mtext></msub></math></span>), steady state crack extension length (<span><math><msub><mi>L</mi><mtext>ss</mtext></msub></math></span>), to quantify this behavior. Additionally, we explore the effect of pre-damage of PAMPS network on the R-curve and develop a theoretical model linking the degree of pre-damage (h) to <span><math><msub><mi>G</mi><mtext>ss</mtext></msub></math></span>. Our model effectively predicts the R-curve of DN gels under various pre-damage conditions. Furthermore, by considering different PAMPS network formulations, we establish a scaling law linking <span><math><mrow><mo>(</mo><mrow><msub><mi>G</mi><mtext>ss</mtext></msub><mo>−</mo><msub><mi>G</mi><mtext>init</mtext></msub></mrow><mo>)</mo></mrow></math></span> with <span><math><msub><mi>L</mi><mtext>ss</mtext></msub></math></span>, thus creating a unified framework for understanding the internal physical properties of diverse DN gels. These insights also differentiate the toughening mechanisms in DN gels: energy dissipation from the PAMPS network affects <span><math><mrow><mo>(</mo><mrow><msub><mi>G</mi><mtext>ss</mtext></msub><mo>−</mo><msub><mi>G</mi><mtext>init</mtext></msub></mrow><mo>)</mo></mrow></math></span>, while the <span><math><msub><mi>G</mi><mtext>init</mtext></msub></math></span> reflects the intrinsic fracture properties governed by the elastic PAAm network. We anticipate that our experimental results and theoretical models will enhance the understanding of fracture in soft tough materials.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"210 ","pages":"Article 106531"},"PeriodicalIF":6.0,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146056278","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
Multiscale modeling of coupled thermo-hydro-mechanical- chemical behavior in hydrate-bearing sediment 含水沉积物热-水-力-化学耦合行为的多尺度模拟
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-04-01 Epub Date: 2026-01-09 DOI: 10.1016/j.jmps.2026.106512
Jidu Yu , Jidong Zhao , Weijian Liang
This study presents a hybrid continuum-discrete multiscale computational framework that integrates the material point method (MPM) and the discrete element method (DEM) to model fully coupled thermo-hydro-mechanical-chemical (THMC) behavior of hydrate-bearing sediments (HBS). Key innovation of the framework lies in its direct use of DEM to model microscale mechanisms, such as hydrate bond degradation, particle rearrangement, and pore evolution, thereby bypassing the need for conventional elastoplastic constitutive models to define effective stress. We show that a simple hydrate saturation-dependent contact model within the DEM can effectively reproduce characteristic shear and volumetric responses of HBS under various hydrate saturation and confining stresses. By embedding a DEM-based representative volume element (RVE) at each material point in the MPM grid, microscale mechanical behaviors are seamlessly homogenized to inform large-deformation macroscale multiphysics processes. Numerical simulations of biaxial compression and indenter penetration demonstrate the framework’s capability to capture critical phenomena, including shear band formation, shear-induced dilation, and the generation of negative excess pore pressure that drives localized hydrate dissociation. The results further reveal that while higher hydrate saturation enhances shear strength, it also promotes brittle failure and intensified dissociation. Conversely, increased confining stress suppresses volumetric dilation and stabilizes the sediment by mitigating the development of negative pore pressure. This multiscale approach provides a powerful new tool for elucidating complex THMC interactions in HBS, with important implications for assessing hydrate-related geohazards and optimizing gas extraction strategies.
本文提出了一种结合物质点法(MPM)和离散元法(DEM)的连续-离散多尺度混合计算框架,对含水沉积物(HBS)的热-水-机械-化学(THMC)完全耦合行为进行建模。该框架的关键创新在于直接使用DEM来模拟微观尺度机制,如水合物键降解、颗粒重排和孔隙演化,从而绕过了传统弹塑性本构模型来定义有效应力的需要。研究表明,在不同水合物饱和度和围应力条件下,简单的水合物饱和度依赖接触模型可以有效地再现HBS的剪切和体积特征响应。通过在MPM网格中的每个材料点嵌入基于dem的代表性体积单元(RVE),可以无缝地均匀化微尺度力学行为,从而为大变形的宏观多物理场过程提供信息。双轴压缩和压头穿透的数值模拟表明,该框架能够捕获关键现象,包括剪切带形成、剪切引起的膨胀以及产生负超孔隙压力,从而驱动局部水合物解离。结果进一步表明,较高的水合物饱和度在提高抗剪强度的同时,也促进了脆性破坏,加剧了解离。相反,围应力的增加抑制了体积膨胀,并通过减缓负孔隙压力的发展来稳定沉积物。这种多尺度方法为阐明HBS中复杂的THMC相互作用提供了强大的新工具,对评估水合物相关地质灾害和优化天然气开采策略具有重要意义。
{"title":"Multiscale modeling of coupled thermo-hydro-mechanical- chemical behavior in hydrate-bearing sediment","authors":"Jidu Yu ,&nbsp;Jidong Zhao ,&nbsp;Weijian Liang","doi":"10.1016/j.jmps.2026.106512","DOIUrl":"10.1016/j.jmps.2026.106512","url":null,"abstract":"<div><div>This study presents a hybrid continuum-discrete multiscale computational framework that integrates the material point method (MPM) and the discrete element method (DEM) to model fully coupled thermo-hydro-mechanical-chemical (THMC) behavior of hydrate-bearing sediments (HBS). Key innovation of the framework lies in its direct use of DEM to model microscale mechanisms, such as hydrate bond degradation, particle rearrangement, and pore evolution, thereby bypassing the need for conventional elastoplastic constitutive models to define effective stress. We show that a simple hydrate saturation-dependent contact model within the DEM can effectively reproduce characteristic shear and volumetric responses of HBS under various hydrate saturation and confining stresses. By embedding a DEM-based representative volume element (RVE) at each material point in the MPM grid, microscale mechanical behaviors are seamlessly homogenized to inform large-deformation macroscale multiphysics processes. Numerical simulations of biaxial compression and indenter penetration demonstrate the framework’s capability to capture critical phenomena, including shear band formation, shear-induced dilation, and the generation of negative excess pore pressure that drives localized hydrate dissociation. The results further reveal that while higher hydrate saturation enhances shear strength, it also promotes brittle failure and intensified dissociation. Conversely, increased confining stress suppresses volumetric dilation and stabilizes the sediment by mitigating the development of negative pore pressure. This multiscale approach provides a powerful new tool for elucidating complex THMC interactions in HBS, with important implications for assessing hydrate-related geohazards and optimizing gas extraction strategies.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"210 ","pages":"Article 106512"},"PeriodicalIF":6.0,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956831","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
Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime 金属晶格超材料的动态破碎:激波模式图和向拓扑无关压实状态的过渡
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-04-01 Epub Date: 2026-01-23 DOI: 10.1016/j.jmps.2026.106515
Brandon K. Zimmerman, Eric B. Herbold, Mukul Kumar, Jonathan Lind
Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. To investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities 0.1≤ρ¯≤0.9 and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.
增材制造的晶格超材料在强度和能量吸收方面提供了设计上的通用性,并通过选择晶格拓扑提供了额外的自由度。在准静态加载下,单胞结构可以强烈地影响刚度、屈服和后屈服行为,但晶格拓扑的影响是否以及在多大程度上持续到动态加载场景,直到压实冲击状态,还没有确定。利用LLNL的ALE3D氢代码对冲击速度为0.25 ~ 2.25 km/s的多种晶格类型(包括回转体、八面体、Schwarz D和菱形十二面体)的动载荷进行了计算研究。我们生成并比较了每个晶格拓扑的震荡Hugoniots,这表明在临界速度以上,结构之间的区别可能不会持久,而压缩晶格的行为相似。为了研究临界速度以上拓扑相关准静态压缩和拓扑无关压缩之间的转变,建立了晶格材料的一维弹性线性硬化塑性致密固体(E-LHP-DS)冲击模型,该模型依赖于受限压缩将准静态力学和冲击力学联系起来。与类似的工作不同,该模型不假设在屈服之前的刚性行为或致密化时的锁定行为,允许对晶格力学进行更丰富的探索。该分析模型只有6个参数,同时适合相对密度0.1≤ρ¯≤0.9的准静态受限压缩模拟,并预测了动态压缩行为,以穿越几种不同的冲击模式,每种冲击模式都由临界冲击速度(相当于临界应力)定义。将数值结果与一维E-LHP-DS冲击模型预测结果进行比较表明,强冲击下的拓扑独立性与致密化的发生有关,致密化可以基于准静态受限压缩结果进行预测。
{"title":"Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime","authors":"Brandon K. Zimmerman,&nbsp;Eric B. Herbold,&nbsp;Mukul Kumar,&nbsp;Jonathan Lind","doi":"10.1016/j.jmps.2026.106515","DOIUrl":"10.1016/j.jmps.2026.106515","url":null,"abstract":"<div><div>Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. To investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities <span><math><mrow><mn>0.1</mn><mo>≤</mo><mover><mi>ρ</mi><mo>¯</mo></mover><mo>≤</mo><mn>0.9</mn></mrow></math></span> and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"210 ","pages":"Article 106515"},"PeriodicalIF":6.0,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146033264","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
Contact in strain gradient elasticity: The rigid flat punch problem 应变梯度弹性中的接触:刚性平冲问题
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-04-01 Epub Date: 2026-01-23 DOI: 10.1016/j.jmps.2026.106527
P.A. Gourgiotis , Th. Zisis , A.E. Giannakopoulos , H.G. Georgiadis
This work revisits the classical flat-punch indentation problem within the framework of Mindlin’s form-II strain-gradient elasticity, uncovering new phenomena driven by the interplay between contact geometry and material length scale. The analysis is carried out under plane strain conditions, and the corresponding mixed boundary value problem is solved using integral equation techniques. Two distinct contact regimes are examined. In the first, assuming full contact beneath the rigid indenter, the pressure distribution exhibits hypersingular behavior with tensile (adhesive-like) tractions near the contact edges. These arise purely from kinematic constraints, without invoking any cohesive law or surface energy. The second regime emerges by relaxing the flatness assumption, allowing for partial separation beneath the punch. In this case, contact is sustained only within a central region, flanked by separation gaps near the edges and balanced by concentrated edge reactions. The resulting pressure is entirely positive and exhibits a classical square-root singularity. Both the contact width and edge forces are shown to depend sensitively on Poisson’s ratio and the material length scale. Beyond a critical length, the contact region collapses, and the problem reduces to the superposition of two Flamant-type concentrated contact solutions. These findings reveal a rich class of indentation responses naturally captured by strain gradient elasticity-phenomena inaccessible to classical continuum models. They may have important implications for nano/micro-indentation experiments on materials with pronounced internal length scales, such as polymers, ceramics, composites, cellular solids, masonry, and biological tissues.
这项工作在Mindlin的ii型应变梯度弹性框架内重新审视了经典的平冲压痕问题,揭示了接触几何形状和材料长度尺度之间相互作用驱动的新现象。在平面应变条件下进行了分析,采用积分方程技术求解了相应的混合边值问题。研究了两种不同的接触机制。在第一种情况下,假设刚性压头下方完全接触,压力分布表现出超奇异行为,在接触边缘附近具有拉伸(类似粘合剂)牵引力。这些纯粹是由运动学约束产生的,而不调用任何内聚定律或表面能。第二种机制是通过放松平面假设,允许冲孔下方部分分离而出现的。在这种情况下,接触仅在中心区域内维持,两侧是靠近边缘的分离间隙,并通过集中的边缘反应来平衡。由此产生的压力完全是正的,并表现出经典的平方根奇点。两者的接触宽度和边缘力显示敏感地依赖于泊松比和材料长度尺度。超过临界长度后,接触区域坍塌,问题就变成了两种flamant型浓缩接触溶液的叠加。这些发现揭示了一种丰富的压痕响应类别,自然地被应变梯度弹性捕获-经典连续体模型无法获得的现象。它们可能对具有明显内部长度尺度的材料(如聚合物、陶瓷、复合材料、细胞固体、砌体和生物组织)的纳米/微压痕实验具有重要意义。
{"title":"Contact in strain gradient elasticity: The rigid flat punch problem","authors":"P.A. Gourgiotis ,&nbsp;Th. Zisis ,&nbsp;A.E. Giannakopoulos ,&nbsp;H.G. Georgiadis","doi":"10.1016/j.jmps.2026.106527","DOIUrl":"10.1016/j.jmps.2026.106527","url":null,"abstract":"<div><div>This work revisits the classical flat-punch indentation problem within the framework of Mindlin’s form-II strain-gradient elasticity, uncovering new phenomena driven by the interplay between contact geometry and material length scale. The analysis is carried out under plane strain conditions, and the corresponding mixed boundary value problem is solved using integral equation techniques. Two distinct contact regimes are examined. In the first, assuming full contact beneath the rigid indenter, the pressure distribution exhibits hypersingular behavior with tensile (adhesive-like) tractions near the contact edges. These arise purely from kinematic constraints, without invoking any cohesive law or surface energy. The second regime emerges by relaxing the flatness assumption, allowing for partial separation beneath the punch. In this case, contact is sustained only within a central region, flanked by separation gaps near the edges and balanced by concentrated edge reactions. The resulting pressure is entirely positive and exhibits a classical square-root singularity. Both the contact width and edge forces are shown to depend sensitively on Poisson’s ratio and the material length scale. Beyond a critical length, the contact region collapses, and the problem reduces to the superposition of two Flamant-type concentrated contact solutions. These findings reveal a rich class of indentation responses naturally captured by strain gradient elasticity-phenomena inaccessible to classical continuum models. They may have important implications for nano/micro-indentation experiments on materials with pronounced internal length scales, such as polymers, ceramics, composites, cellular solids, masonry, and biological tissues.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"210 ","pages":"Article 106527"},"PeriodicalIF":6.0,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146033263","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
Generalized configurational force for thermoelasticity with application to isogeometric analysis of thermoelastic crack propagation 热弹性广义构形力及其在热弹性裂纹扩展等几何分析中的应用
IF 6 2区 工程技术 Q2 MATERIALS SCIENCE, MULTIDISCIPLINARY Pub Date : 2026-04-01 Epub Date: 2026-02-03 DOI: 10.1016/j.jmps.2026.106541
Chung-Shuo Lee, Pei-En Chou, Ganesh Subbarayan
Accurately predicting the forces driving crack propagation is critical for engineering applications where fracture is an important failure criterion. An easy to apply configurational force theory for estimating crack driving force in thermoelastic solids is derived in this study. The derivation relies on the notion of a configurational optimization problem for heterogeneities inserted into a homogeneous domain. The configurational optimization allows arbitrary objectives subject to constraints of variational principles governing the solid behavior. The optimality criteria obtained by solving an adjoint problem naturally yields a generalized Eshelby Energy-Momentum tensor, which leads to the classical J−, L− and M− integrals of fracture mechanics when Helmholtz free energy is chosen as the objective and the problem is constrained to be isothermal. We calculate the configurational force using an isogeometric computational framework that relies on a recently developed enrichment procedure for multi-material corner singularities, of which cracks are a special case. Enriched displacement and temperature functions applied at crack tips are illustrated and demonstrated through several numerical examples on thermoelastic solids. Examples include (a) crack driving force estimation under adiabatic/isothermal conditions, (b) edge crack propagation simulation and (c) bimaterial wedge under varying mechanical and thermal loadings.
在断裂是重要破坏准则的工程应用中,准确预测驱动裂纹扩展的力至关重要。本文推导了一种易于应用的热弹性固体裂纹驱动力计算的构形力理论。推导依赖于插入到齐次域中的异构的构型优化问题的概念。结构优化允许任意目标服从于控制固体行为的变分原则的约束。通过求解伴随问题得到的最优性准则自然产生广义的Eshelby能量-动量张量,当选择Helmholtz自由能作为目标并将问题限制为等温时,可得到经典的断裂力学J−,L−和M−积分。我们使用一个等几何计算框架来计算构型力,该框架依赖于最近开发的多材料角点奇点的浓缩程序,其中裂纹是一个特殊情况。通过几个热弹性固体的数值例子,说明了裂纹尖端的丰富位移和温度函数。例子包括(a)绝热/等温条件下的裂纹驱动力估计,(b)边缘裂纹扩展模拟和(c)不同机械和热载荷下的双材料楔。
{"title":"Generalized configurational force for thermoelasticity with application to isogeometric analysis of thermoelastic crack propagation","authors":"Chung-Shuo Lee,&nbsp;Pei-En Chou,&nbsp;Ganesh Subbarayan","doi":"10.1016/j.jmps.2026.106541","DOIUrl":"10.1016/j.jmps.2026.106541","url":null,"abstract":"<div><div>Accurately predicting the forces driving crack propagation is critical for engineering applications where fracture is an important failure criterion. An easy to apply configurational force theory for estimating crack driving force in thermoelastic solids is derived in this study. The derivation relies on the notion of a configurational optimization problem for heterogeneities inserted into a homogeneous domain. The configurational optimization allows arbitrary objectives subject to constraints of variational principles governing the solid behavior. The optimality criteria obtained by solving an adjoint problem naturally yields a generalized Eshelby Energy-Momentum tensor, which leads to the classical <span><math><mrow><mi>J</mi><mo>−</mo></mrow></math></span>, <span><math><mrow><mi>L</mi><mo>−</mo></mrow></math></span> and <span><math><mrow><mi>M</mi><mo>−</mo></mrow></math></span> integrals of fracture mechanics when Helmholtz free energy is chosen as the objective and the problem is constrained to be isothermal. We calculate the configurational force using an isogeometric computational framework that relies on a recently developed enrichment procedure for multi-material corner singularities, of which cracks are a special case. Enriched displacement and temperature functions applied at crack tips are illustrated and demonstrated through several numerical examples on thermoelastic solids. Examples include (a) crack driving force estimation under adiabatic/isothermal conditions, (b) edge crack propagation simulation and (c) bimaterial wedge under varying mechanical and thermal loadings.</div></div>","PeriodicalId":17331,"journal":{"name":"Journal of The Mechanics and Physics of Solids","volume":"210 ","pages":"Article 106541"},"PeriodicalIF":6.0,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146110149","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
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
期刊
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