{"title":"Robust Inverse Material Design With Physical Guarantees Using the Voigt-Reuss Net","authors":"Sanath Keshav, Felix Fritzen","doi":"10.1002/nme.70296","DOIUrl":null,"url":null,"abstract":"<p>We apply the Voigt-Reuss net, a spectrally normalized neural surrogate introduced in [38], for forward and inverse mechanical homogenization with a key guarantee that all predicted effective stiffness tensors satisfy Voigt-Reuss bounds in the Löwner sense during training, inference, and gradient-driven optimization. The approach operates in a bounded spectral space by reparametrizing each effective tensor relative to its Voigt and Reuss bounds, ensuring that all outputs reside within a unit-cube domain and are mapped back via a deterministic inverse transform to physically admissible tensors. For 3D elasticity, a fully connected Voigt-Reuss net is trained on <span></span><math>\n <semantics>\n <mrow>\n <mo>∼</mo>\n <mn>1</mn>\n <mo>.</mo>\n <mn>18</mn>\n </mrow>\n <annotation>$$ \\sim 1.18 $$</annotation>\n </semantics></math> million high-fidelity homogenization labels of stochastic biphasic microstructures. The model ingests 236 image-derived morphological descriptors and phase parameters that encode bulk and shear moduli, enabling a single surrogate to represent material combinations spanning 4 orders of magnitude. Due to the rotational invariance of the descriptor set, the surrogate recovers the isotropic projection of the effective stiffness (<span></span><math>\n <semantics>\n <mrow>\n <msup>\n <mrow>\n <mi>R</mi>\n </mrow>\n <mrow>\n <mn>2</mn>\n </mrow>\n </msup>\n <mo>≥</mo>\n <mn>0</mn>\n <mo>.</mo>\n <mn>998</mn>\n </mrow>\n <annotation>$$ {R}^2\\ge 0.998 $$</annotation>\n </semantics></math> for isotropy-related components). However, anisotropy-revealing entries remain unlearnable from the available features. At the tensor level, the relative Frobenius error exhibits a median of approximately 1.8% (mean approximately 3.6%) and approaches the irreducible isotropic-projection floor, far outperforming all alternative surrogates considered. In 2D plane-strain elasticity, spectral normalization is integrated with a differentiable microstructure renderer and a convolutional regressor, yielding a surrogate that maps generator parameters to effective stiffness tensors for highly anisotropic and high-contrast composites. Voigt-Reuss net is compared against vanilla and Cholesky regressors trained with identical architectures, data, and training procedures. The unconstrained surrogates frequently violate Voigt/Reuss bounds and even positive definiteness, whereas Voigt-Reuss net produces no violations by design and also enables robust inverse design. Utilizing the surrogate with batched first-order optimization, the approach can match prescribed target tensors and optimize tensor functionals, recovering classical optima while avoiding nonphysical and spurious designs observed with unconstrained surrogates.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"127 7","pages":""},"PeriodicalIF":2.9000,"publicationDate":"2026-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.70296","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"International Journal for Numerical Methods in Engineering","FirstCategoryId":"5","ListUrlMain":"https://onlinelibrary.wiley.com/doi/10.1002/nme.70296","RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"ENGINEERING, MULTIDISCIPLINARY","Score":null,"Total":0}
引用次数: 0
Abstract
We apply the Voigt-Reuss net, a spectrally normalized neural surrogate introduced in [38], for forward and inverse mechanical homogenization with a key guarantee that all predicted effective stiffness tensors satisfy Voigt-Reuss bounds in the Löwner sense during training, inference, and gradient-driven optimization. The approach operates in a bounded spectral space by reparametrizing each effective tensor relative to its Voigt and Reuss bounds, ensuring that all outputs reside within a unit-cube domain and are mapped back via a deterministic inverse transform to physically admissible tensors. For 3D elasticity, a fully connected Voigt-Reuss net is trained on million high-fidelity homogenization labels of stochastic biphasic microstructures. The model ingests 236 image-derived morphological descriptors and phase parameters that encode bulk and shear moduli, enabling a single surrogate to represent material combinations spanning 4 orders of magnitude. Due to the rotational invariance of the descriptor set, the surrogate recovers the isotropic projection of the effective stiffness ( for isotropy-related components). However, anisotropy-revealing entries remain unlearnable from the available features. At the tensor level, the relative Frobenius error exhibits a median of approximately 1.8% (mean approximately 3.6%) and approaches the irreducible isotropic-projection floor, far outperforming all alternative surrogates considered. In 2D plane-strain elasticity, spectral normalization is integrated with a differentiable microstructure renderer and a convolutional regressor, yielding a surrogate that maps generator parameters to effective stiffness tensors for highly anisotropic and high-contrast composites. Voigt-Reuss net is compared against vanilla and Cholesky regressors trained with identical architectures, data, and training procedures. The unconstrained surrogates frequently violate Voigt/Reuss bounds and even positive definiteness, whereas Voigt-Reuss net produces no violations by design and also enables robust inverse design. Utilizing the surrogate with batched first-order optimization, the approach can match prescribed target tensors and optimize tensor functionals, recovering classical optima while avoiding nonphysical and spurious designs observed with unconstrained surrogates.
我们应用Voigt-Reuss网络([38]中引入的频谱归一化神经代理)进行正向和逆机械均匀化,并在训练、推理和梯度驱动优化期间保证所有预测的有效刚度张量在Löwner意义上满足Voigt-Reuss界。该方法通过相对于Voigt和Reuss边界重新参数化每个有效张量来在有界光谱空间中操作,确保所有输出都位于单位立方体域中,并通过确定性逆变换映射回物理上可接受的张量。对于3D弹性,在~ 1上训练一个完全连接的voight - reuss网络。18 $$ \sim 1.18 $$万个高保真均质随机双相微观结构标签。该模型摄取236个图像衍生的形态描述符和相位参数,编码体积和剪切模量,使单个代理能够表示跨越4个数量级的材料组合。由于描述子集的旋转不变性,代理恢复有效刚度(r2≥0)的各向同性投影。998 $$ {R}^2\ge 0.998 $$与各向同性相关的组件)。然而,揭示各向异性的条目仍然无法从可用的特征中学习。在张量水平上,相对Frobenius误差的中位数约为1.8% (mean approximately 3.6%) and approaches the irreducible isotropic-projection floor, far outperforming all alternative surrogates considered. In 2D plane-strain elasticity, spectral normalization is integrated with a differentiable microstructure renderer and a convolutional regressor, yielding a surrogate that maps generator parameters to effective stiffness tensors for highly anisotropic and high-contrast composites. Voigt-Reuss net is compared against vanilla and Cholesky regressors trained with identical architectures, data, and training procedures. The unconstrained surrogates frequently violate Voigt/Reuss bounds and even positive definiteness, whereas Voigt-Reuss net produces no violations by design and also enables robust inverse design. Utilizing the surrogate with batched first-order optimization, the approach can match prescribed target tensors and optimize tensor functionals, recovering classical optima while avoiding nonphysical and spurious designs observed with unconstrained surrogates.
期刊介绍:
The International Journal for Numerical Methods in Engineering publishes original papers describing significant, novel developments in numerical methods that are applicable to engineering problems.
The Journal is known for welcoming contributions in a wide range of areas in computational engineering, including computational issues in model reduction, uncertainty quantification, verification and validation, inverse analysis and stochastic methods, optimisation, element technology, solution techniques and parallel computing, damage and fracture, mechanics at micro and nano-scales, low-speed fluid dynamics, fluid-structure interaction, electromagnetics, coupled diffusion phenomena, and error estimation and mesh generation. It is emphasized that this is by no means an exhaustive list, and particularly papers on multi-scale, multi-physics or multi-disciplinary problems, and on new, emerging topics are welcome.