Jeremy A. McCulloch, Skyler R. St. Pierre, Kevin Linka, Ellen Kuhl
{"title":"On sparse regression, Lp-regularization, and automated model discovery","authors":"Jeremy A. McCulloch, Skyler R. St. Pierre, Kevin Linka, Ellen Kuhl","doi":"10.1002/nme.7481","DOIUrl":"10.1002/nme.7481","url":null,"abstract":"<p>Sparse regression and feature extraction are the cornerstones of knowledge discovery from massive data. Their goal is to discover interpretable and predictive models that provide simple relationships among scientific variables. While the statistical tools for model discovery are well established in the context of linear regression, their generalization to nonlinear regression in material modeling is highly problem-specific and insufficiently understood. Here we explore the potential of neural networks for automatic model discovery and induce sparsity by a hybrid approach that combines two strategies: regularization and physical constraints. We integrate the concept of <i>L</i><sub><i>p</i></sub> regularization for subset selection with constitutive neural networks that leverage our domain knowledge in kinematics and thermodynamics. We train our networks with both, synthetic and real data, and perform several thousand discovery runs to infer common guidelines and trends: <i>L</i><sub>2</sub> regularization or ridge regression is unsuitable for model discovery; <i>L</i><sub>1</sub> regularization or lasso promotes sparsity, but induces strong bias that may aggressively change the results; only <i>L</i><sub>0</sub> regularization allows us to transparently fine-tune the trade-off between interpretability and predictability, simplicity and accuracy, and bias and variance. With these insights, we demonstrate that <i>L</i><sub><i>p</i></sub> regularized constitutive neural networks can simultaneously discover both, interpretable models and physically meaningful parameters. We anticipate that our findings will generalize to alternative discovery techniques such as sparse and symbolic regression, and to other domains such as biology, chemistry, or medicine. Our ability to automatically discover material models from data could have tremendous applications in generative material design and open new opportunities to manipulate matter, alter properties of existing materials, and discover new materials with user-defined properties.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"125 14","pages":""},"PeriodicalIF":2.9,"publicationDate":"2024-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7481","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140588322","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Extension of the a posteriori finite element method (APFEM) to geometrical alterations and application to stochastic homogenisation","authors":"Yanis Ammouche, Antoine Jérusalem","doi":"10.1002/nme.7482","DOIUrl":"10.1002/nme.7482","url":null,"abstract":"<p>We recently proposed an efficient method facilitating the parametric study of a finite element mechanical simulation as a postprocessing step, that is, without the need to run multiple simulations: the a posteriori finite element method (APFEM). APFEM only requires the knowledge of the vertices of the parameter space and is able to predict accurately how the degrees of freedom of a simulation, i.e., nodal displacements, and other outputs of interests, for example, element stress tensors, evolve when simulation parameters vary within their predefined ranges. In our previous work, these parameters were restricted to material properties and loading conditions. Here, we extend the APFEM to additionally account for changes in the original geometry. This is achieved by defining an intermediary reference frame whose mapping is defined stochastically in the weak form. Subsequent deformation is then reached by correcting for this stochastic variation in the reference frame through multiplicative decomposition of the deformation gradient tensor. The resulting framework is shown here to provide accurate mechanical predictions for relevant applications of increasing complexity: (i) quantifying the stress concentration factor of a plate under uniaxial loading with one and two elliptical holes of varying eccentricities, and (ii) performing the stochastic homogenisation of a composite plate with uncertain mechanical properties and geometry inclusion. This extension of APFEM completes our original approach to account parametrically for geometrical alterations, in addition to boundary conditions and material properties. The advantages of this approach in our original work in terms of stochastic prediction, uncertainty quantification, structural and material optimisation and Bayesian inferences are all naturally conserved.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"125 14","pages":""},"PeriodicalIF":2.9,"publicationDate":"2024-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7482","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140588656","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Higher-order generalized-α methods for parabolic problems","authors":"Pouria Behnoudfar, Quanling Deng, Victor M. Calo","doi":"10.1002/nme.7485","DOIUrl":"10.1002/nme.7485","url":null,"abstract":"<p>We propose a new class of high-order time-marching schemes with dissipation control and unconditional stability for parabolic equations. High-order time integrators can deliver the optimal performance of highly accurate and robust spatial discretizations such as isogeometric analysis. The generalized-<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>α</mi>\u0000 </mrow>\u0000 <annotation>$$ alpha $$</annotation>\u0000 </semantics></math> method delivers unconditional stability and second-order accuracy in time and controls the numerical dissipation in the discrete spectrum's high-frequency region. We extend the generalized-<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>α</mi>\u0000 </mrow>\u0000 <annotation>$$ alpha $$</annotation>\u0000 </semantics></math> methodology to obtain high-order time marching methods with high accuracy and dissipation control in the discrete high-frequency range. Furthermore, we maintain the original stability region of the second-order generalized-<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>α</mi>\u0000 </mrow>\u0000 <annotation>$$ alpha $$</annotation>\u0000 </semantics></math> method in the new higher-order methods; we increase the accuracy of the generalized-<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>α</mi>\u0000 </mrow>\u0000 <annotation>$$ alpha $$</annotation>\u0000 </semantics></math> method while keeping the unconditional stability and user-control features on the high-frequency numerical dissipation. The methodology solves <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>k</mi>\u0000 <mo>></mo>\u0000 <mn>1</mn>\u0000 <mo>,</mo>\u0000 <mi>k</mi>\u0000 <mo>∈</mo>\u0000 <mi>ℕ</mi>\u0000 </mrow>\u0000 <annotation>$$ k>1,kin mathbb{N} $$</annotation>\u0000 </semantics></math> matrix problems and updates the system unknowns, which correspond to higher-order terms in Taylor expansions to obtain <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>(</mo>\u0000 <mn>3</mn>\u0000 <mo>/</mo>\u0000 <mn>2</mn>\u0000 <mi>k</mi>\u0000 <mo>)</mo>\u0000 <mtext>th</mtext>\u0000 </mrow>\u0000 <annotation>$$ left(3/2kright)mathrm{th} $$</annotation>\u0000 </semantics></math>-order method for even <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>k</mi>\u0000 </mrow>\u0000 <annotation>$$ k $","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"125 13","pages":""},"PeriodicalIF":2.9,"publicationDate":"2024-04-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7485","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140588590","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Richard Schussnig, Malte Rolf-Pissarczyk, Kathrin Bäumler, Thomas-Peter Fries, Gerhard A. Holzapfel, Martin Kronbichler
{"title":"On the role of tissue mechanics in fluid–structure interaction simulations of patient-specific aortic dissection","authors":"Richard Schussnig, Malte Rolf-Pissarczyk, Kathrin Bäumler, Thomas-Peter Fries, Gerhard A. Holzapfel, Martin Kronbichler","doi":"10.1002/nme.7478","DOIUrl":"10.1002/nme.7478","url":null,"abstract":"<p>Modeling an aortic dissection represents a particular challenge from a numerical perspective, especially when it comes to the interaction between solid (aortic wall) and liquid (blood flow). The complexity of patient-specific simulations requires a variety of parameters, modeling assumptions and simplifications that currently hinder their routine use in clinical settings. We present a numerical framework that captures, among other things, the layer-specific anisotropic properties of the aortic wall, the non-Newtonian behavior of blood, patient-specific geometry, and patient-specific flow conditions. We compare hemodynamic indicators and stress measurements in simulations with increasingly complex material models for the vessel tissue ranging from rigid walls to anisotropic hyperelastic materials. We find that for the present geometry and boundary conditions, rigid wall simulations produce different results than fluid–structure interaction simulations. Considering anisotropic fiber contributions in the tissue model, stress measurements in the aortic wall differ, but shear stress-based biomarkers are less affected. In summary, the increasing complexity of the tissue model enables capturing more details. However, an extensive parameter set is also required. Since the simulation results depend on these modeling choices, variations can lead to different recommendations in clinical applications.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"125 14","pages":""},"PeriodicalIF":2.9,"publicationDate":"2024-04-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7478","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140588323","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
D. Yago, G. Sal-Anglada, D. Roca, J. Cante, J. Oliver
{"title":"Machine learning in solid mechanics: Application to acoustic metamaterial design","authors":"D. Yago, G. Sal-Anglada, D. Roca, J. Cante, J. Oliver","doi":"10.1002/nme.7476","DOIUrl":"10.1002/nme.7476","url":null,"abstract":"<p>Machine learning (ML) and Deep learning (DL) are increasingly pivotal in the design of advanced metamaterials, seamlessly integrated with material or topology optimization. Their intrinsic capability to predict and interconnect material properties across vast design spaces, often computationally prohibitive for conventional methods, has led to groundbreaking possibilities. This paper introduces an innovative machine learning approach for the optimization of acoustic metamaterials, focusing on Multiresonant Layered Acoustic Metamaterial (MLAM), designed for targeted noise attenuation at low frequencies (below 1000 Hz). This method leverages ML to create a continuous model of the Representative Volume Element (RVE) effective properties essential for evaluating sound transmission loss (STL), and subsequently used to optimize the overall topology configuration for maximum sound attenuation using a Genetic Algorithm (GA). The significance of this methodology lies in its ability to deliver rapid results without compromising accuracy, significantly reducing the computational overhead of complete topology optimization by several orders of magnitude. To demonstrate the versatility and scalability of this approach, it is extended to a more intricate RVE model, characterized by a higher number of parameters, and is optimized using the same strategy. In addition, to underscore the potential of ML techniques in synergy with traditional topology optimization, a comparative analysis is conducted, comparing the outcomes of the proposed method with those obtained through direct numerical simulation (DNS) of the corresponding full 3D MLAM model. This comparative analysis highlights the transformative potential of this combination, particularly when addressing complex topological challenges with significant computational demands, ushering in a new era of metamaterial and component design.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"125 14","pages":""},"PeriodicalIF":2.9,"publicationDate":"2024-04-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7476","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140588950","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Prajwal Kammardi Arunachala, Sina Abrari Vajari, Matthias Neuner, Jay Sejin Sim, Renee Zhao, Christian Linder
{"title":"A multiscale anisotropic polymer network model coupled with phase field fracture","authors":"Prajwal Kammardi Arunachala, Sina Abrari Vajari, Matthias Neuner, Jay Sejin Sim, Renee Zhao, Christian Linder","doi":"10.1002/nme.7488","DOIUrl":"10.1002/nme.7488","url":null,"abstract":"<p>The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. The results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"125 13","pages":""},"PeriodicalIF":2.9,"publicationDate":"2024-04-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140588592","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Asymptotic homogenization of phase-field fracture model: An efficient multiscale finite element framework for anisotropic fracture","authors":"Pu-Song Ma, Xing-Cheng Liu, Xue-Ling Luo, Shaofan Li, Lu-Wen Zhang","doi":"10.1002/nme.7489","DOIUrl":"10.1002/nme.7489","url":null,"abstract":"<p>The intractable multiscale constitutives and the high computational cost in direct numerical simulations are the bottlenecks in fracture analysis of heterogeneous materials. In an attempt to achieve a balance between accuracy and efficiency, we propose a mathematically rigorous phase-field model for multiscale fracture. Leveraging the phase-field theory, the difficulty of discrete-continuous coupling in conventional cross-scale crack propagation analysis is resolved by constructing a continuum description of the crack. Based on the asymptotic expansion, an equivalent two-field coupled boundary-value problem is well-defined, from which we rigorously derive the macroscopic equivalent parameters, including the equivalent elasticity tensor and the equivalent fracture toughness tensor. In our approach, both the displacement field and the phase-field are simultaneously expanded, allowing us to obtain a fracture toughness tensor with diagonal elements of the corresponding matrix controlling anisotropic fracture behavior and non-diagonal elements governing crack deflection. This enables multiscale finite element homogenization procedure to accurately reproduce microstructural information, and capture the crack deflection angle in anisotropic materials without any a priori knowledge. From the numerical results, the proposed multiscale phase-field method demonstrates a significant reduction in computation time with respect to full-field simulations. Moreover, the method accurately reproduces physical consistent anisotropic fracture of non-centrosymmetric porous media, and the experimentally consistent damage response of fiber-reinforced composites. This work fuses well-established mathematical homogenization theory with the cutting-edge fracture phase-field method, sparking a fresh perspective for the fracture of heterogeneous media.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"125 13","pages":""},"PeriodicalIF":2.9,"publicationDate":"2024-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140588413","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Featured Cover","authors":"Qian Yang, Xiaoqin Shen, Jikun Zhao, Yumin Cheng","doi":"10.1002/nme.7484","DOIUrl":"https://doi.org/10.1002/nme.7484","url":null,"abstract":"<p>The cover image is based on the Research Article <i>An enriched virtual element method for 2D-3C generalized membrane shell model on surface</i> by Qian Yang et al., https://doi.org/10.1002/nme.7432.\u0000\u0000 <figure>\u0000 <div><picture>\u0000 <source></source></picture><p></p>\u0000 </div>\u0000 </figure></p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"125 8","pages":""},"PeriodicalIF":2.9,"publicationDate":"2024-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7484","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140333081","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Anderson Soares da Costa Azevêdo, Hao Li, Naouyuki Ishida, Lucas Oliveira Siqueira, Rômulo Luz Cortez, Emílio Carlos Nelli Silva, Shinji Nishiwaki, Renato Picelli
{"title":"Body-fitted topology optimization via integer linear programming using surface capturing techniques","authors":"Anderson Soares da Costa Azevêdo, Hao Li, Naouyuki Ishida, Lucas Oliveira Siqueira, Rômulo Luz Cortez, Emílio Carlos Nelli Silva, Shinji Nishiwaki, Renato Picelli","doi":"10.1002/nme.7480","DOIUrl":"10.1002/nme.7480","url":null,"abstract":"<p>Recent advancements in computational tools and additive manufacturing have expanded design possibilities on fluid devices and structures to also include aesthetics. However, traditional discrete density-based topology optimization methods usually use square and cubic regular meshes, resulting in jagged contour patterns that require mesh refinement and post-processing to numerical solution steps of complex physics problems. In this article, we propose a new isosurface boundary capture strategy for topology optimization in structural and fluid flow problems. The capture of smoothed boundaries is done via a simple strategy through element splitting and analysis domain remeshing. The smoothing procedure novelty lies on the optimizer regular mesh decomposition into smaller triangular or tetrahedral elements with a pseudo-density nodal function that produces implicit geometry boundaries. We employ the TOBS-GT (topology optimization of binary structures with geometry trimming) method to solve optimization problems for mean compliance and fluid flow energy dissipation, subject to a volume fraction constraint. Since, we perform discrete body-fitted optimization, the material interpolation models are expressed in linear form without penalty factors. Our method produces a lower computational cost topology evolution with high resolution and mitigated sharp details, achieved via smooth edges and surfaces, as demonstrated through benchmark numerical examples in two- and three-dimensional space. In addition, the proposed strategy facilitates the optimization of benchmark fluid flow examples for moderate Reynolds numbers.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"125 13","pages":""},"PeriodicalIF":2.9,"publicationDate":"2024-03-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140315799","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jakob Platen, Bennett Pauls, Atul Anantheswar, Thea Lautenschläger, Christoph Neinhuis, Michael Kaliske
{"title":"A nonlinear finite viscoelastic formulation relative to the viscous intermediate configuration applied to plants","authors":"Jakob Platen, Bennett Pauls, Atul Anantheswar, Thea Lautenschläger, Christoph Neinhuis, Michael Kaliske","doi":"10.1002/nme.7483","DOIUrl":"10.1002/nme.7483","url":null,"abstract":"<p>In the contribution at hand, a new formulation for finite strain viscosity relative to the viscous intermediate configuration is presented. The evolution of the viscous deformations is based upon a new numerical approach, which allows for a consistent consideration of anisotropic finite strain viscoelasticity, according to the authors knowledge. A standard Maxwell model is used to describe viscous behaviour at finite deformations. Furthermore, the orthotropic Yeoh material model is extended to include a distinction between behaviour under tensile and compression loading. The proposed formulation is validated, and parameters of the model are identified by material tests on <i>Sorghum bicolor</i> plants. Subsequently, numerical examples are shown to demonstrate the capabilities of the model. In general, the proposed Yeoh material formulation is shown to accurately represent the inability of fibres to carry compression loading. Furthermore, the viscoelastic approach, developed relative to the viscous intermediate configuration, is demonstrated to be capable of producing plausible results. Additionally, the mechanical behaviour of <i>Sorghum bicolor</i> plants is simulated using the introduced formulation. The results show that the contribution at hand describes a novel methodology to simulate the viscoelastic behaviour of plant materials reliably.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"125 13","pages":""},"PeriodicalIF":2.9,"publicationDate":"2024-03-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7483","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140201096","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}