Personalized computational simulations have emerged as a vital tool to understand the biomechanical factors of a disease, predict disease progression, and design personalized intervention. Material modeling is critical for realistic biomedical simulations, and poor model selection can have life-threatening consequences for the patient. However, selecting the best model requires a profound domain knowledge and is limited to a few highly specialized experts in the field. Here we explore the feasibility of eliminating user involvement and automate the process of material modeling in finite element analyses. We leverage recent developments in constitutive neural networks, machine learning, and artificial intelligence to discover the best constitutive model from thousands of possible combinations of a few functional building blocks. We integrate all discoverable models into the finite element workflow by creating a universal material subroutine that contains more than 60,000 models, made up of 16 individual terms. We prototype this workflow using biaxial extension tests from healthy human arteries as input and stress and stretch profiles across the human aortic arch as output. Our results suggest that constitutive neural networks can robustly discover various flavors of arterial models from data, feed these models directly into a finite element simulation, and predict stress and strain profiles that compare favorably to the classical Holzapfel model. Replacing dozens of individual material subroutines by a single universal material subroutine—populated directly via automated model discovery—will make finite element simulations more user-friendly, more robust, and less vulnerable to human error. Democratizing finite element simulation by automating model selection could induce a paradigm shift in physics-based modeling, broaden access to simulation technologies, and empower individuals with varying levels of expertise and diverse backgrounds to actively participate in scientific discovery and push the boundaries of biomedical simulation.
Soft materials play an integral part in many aspects of modern life including autonomy, sustainability, and human health, and their accurate modeling is critical to understand their unique properties and functions. Today’s finite element analysis packages come with a set of pre-programmed material models, which may exhibit restricted validity in capturing the intricate mechanical behavior of these materials. Regrettably, incorporating a modified or novel material model in a finite element analysis package requires non-trivial in-depth knowledge of tensor algebra, continuum mechanics, and computer programming, making it a complex task that is prone to human error. Here we design a universal material subroutine, which automates the integration of novel constitutive models of varying complexity in non-linear finite element packages, with no additional analytical derivations and algorithmic implementations. We demonstrate the versatility of our approach to seamlessly integrate innovative constitutive models from the material point to the structural level through a variety of soft matter case studies: a frontal impact to the brain; reconstructive surgery of the scalp; diastolic loading of arteries and the human heart; and the dynamic closing of the tricuspid valve. Our universal material subroutine empowers all users, not solely experts, to conduct reliable engineering analysis of soft matter systems. We envision that this framework will become an indispensable instrument for continued innovation and discovery within the soft matter community at large.
more » « less- PAR ID:
- 10543181
- Publisher / Repository:
- Springer Science + Business Media
- Date Published:
- Journal Name:
- Engineering with Computers
- ISSN:
- 0177-0667
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
Nano-indentation is a promising method to identify the constitutive parameters of soft materials, including soft tissues. Especially when materials are very small and heterogeneous, nano-indentation allows mechanical interrogation where traditional methods may fail. However, because nano-indentation does not yield a homogeneous deformation field, interpreting the resulting load–displacement curves is non-trivial and most investigators resort to simplified approaches based on the Hertzian solution. Unfortunately, for small samples and large indentation depths, these solutions are inaccurate. We set out to use machine learning to provide an alternative strategy. We first used the finite element method to create a large synthetic data set. We then used these data to train neural networks to inversely identify material parameters from load–displacement curves. To this end, we took two different approaches. First, we learned the indentation forward problem, which we then applied within an iterative framework to identify material parameters. Second, we learned the inverse problem of directly identifying material parameters. We show that both approaches are effective at identifying the parameters of the neo-Hookean and Gent models. Specifically, when applied to synthetic data, our approaches are accurate even for small sample sizes and at deep indentation. Additionally, our approaches are fast, especially compared to the inverse finite element approach. Finally, our approaches worked on unseen experimental data from thin mouse brain samples. Here, our approaches proved robust to experimental noise across over 1000 samples. By providing open access to our data and code, we hope to support others that conduct nano-indentation on soft materials.more » « less
-
more » « less
Active origami-inspired designs, which incorporate active materials such as electroactive polymers and magnetoactive elastomers into self-folding structures, have shown good promise in engineering applications. In this article, finite element analysis models are developed for several bending and folding configurations that incorporate a combination of active and passive material layers, such as electroactive polymer–actuated unimorph benders based on single and multilayer electroactive polymers, electroactive polymer–actuated notched unimorph folding configurations, and a multi-field actuated bimorph configuration. Constitutive relations are developed for both electrostrictive and magnetoactive materials to model the coupled behaviors explicitly. Shell elements are adopted for their capacity of modeling thin films, relatively low computational cost, and ability to model the intrinsic coupled behaviors in the active materials under consideration. The electrostrictive coefficients are measured and then used as input in the constitutive modeling of the coupled behavior. The magnetization of the magnetoactive elastomer is measured and then used to calculate the magnetic torque as a function of the special orientation, which leads to spatial deformation of the magnetoactive elastomers. The objective of the study is to validate the constitutive models implemented through the finite element analysis method to simulate multi-field coupled behaviors of the active origami configurations. Through quantitative comparisons, simulation results show good agreement with experimental data. By investigating the impact of material selection, location, and geometric parameters, this finite element analysis approach can be used in design of self-folding structures, reducing trial-and-error iterations in experiments.
-
Abstract Entropy dynamics is a Bayesian inference methodology that can be used to quantify time-dependent posterior probability densities that guide the development of complex material models using information theory. Here, we expand its application to non-Gaussian processes to evaluate how fractal structure can influence fractional hyperelasticity and viscoelasticity in elastomers. We investigate how kinematic constraints on fractal polymer network deformation influences the form of hyperelastic constitutive behavior and viscoelasticity in soft materials such as dielectric elastomers, which have applications in the development of adaptive structures. The modeling framework is validated on two dielectric elastomers, VHB 4910 and 4949, over a broad range of stretch rates. It is shown that local fractal time derivatives are equally effective at predicting viscoelasticity in these materials in comparison to nonlocal fractional time derivatives under constant stretch rates. We describe the origin of this accuracy that has implications for simulating large-scale problems such as finite element analysis given the differences in computational efficiency of nonlocal fractional derivatives versus local fractal derivatives.more » « less
-
Constitutive models are important to biomechanics for two key reasons. First, constitutive modelling is an essential component of characterizing tissues’ mechanical properties for informing theoretical and computational models of biomechanical systems. Second, constitutive models can be used as a theoretical framework for extracting and comparing key quantities of interest from material characterization experiments. Over the past five decades, the Ogden model has emerged as a popular constitutive model in soft tissue biomechanics with relevance to both informing theoretical and computational models and to comparing material characterization experiments. The goal of this short review is threefold. First, we will discuss the broad relevance of the Ogden model to soft tissue biomechanics and the general characteristics of soft tissues that are suitable for approximating with the Ogden model. Second, we will highlight exemplary uses of the Ogden model in brain tissue, blood clot and other tissues. Finally, we offer a tutorial on fitting the one-term Ogden model to pure shear experimental data via both an analytical approximation of homogeneous deformation and a finite-element model of the tissue domain. Overall, we anticipate that this short review will serve as a practical introduction to the use of the Ogden model in biomechanics. This article is part of the theme issue ‘The Ogden model of rubber mechanics: Fifty years of impact on nonlinear elasticity’.more » « less
