Nonlinear Material Parameter Estimation Using Inverse Finite Element Analysis

Author(s):  
Jeffrey E. Bischoff

Constitutive parameters for biological materials are ideally regressed against data from well-designed experiments in which the loading and boundary conditions give rise to a homogeneous region of deformation. Such conditions may exist for healthy tissue within the context of in vitro tests, but rarely are met when attempting to measure parameters physiologically or noninvasively, due to complex boundary conditions or heterogeneous material structure and properties. The ability to estimate parameters in these situations is essential in many clinically relevant studies, including determination of tendon/ligament parameters in whole knee studies, non-destructive evaluation of evolving material parameters in laboratory studies, and estimation of heterogeneous parameters due to local normal or pathologic disruptions in tissue microstructure. In such cases, computational algorithms must be used to regress material parameters for a given constitutive model against the available data, in which the experimental conditions are modeled as accurately as possible without significant regard to complexity. The work presented here is focused on development of an iterative, inverse finite element (FE) algorithm for estimation of material parameters from experimental data obtained from tests with nonlinear complexities from contact, large deformations, and constitutive models.

2001 ◽  
Vol 123 (4) ◽  
pp. 333-340 ◽  
Author(s):  
W. L. Dunbar, ◽  
K. U¨n ◽  
P. S. Donzelli ◽  
R. L. Spilker

We have developed an approximate method for simulating the three-dimensional contact of soft biphasic tissues in diarthrodial joints under physiological loading. Input to the method includes: (i) kinematic information describing an in vitro joint articulation, measured while the cartilage is deformed under physiological loads, (ii) geometric properties for the relaxed (undeformed) cartilage layers, obtained for the analyses in this study via stereophotogrammetry, and (iii) material parameters for the biphasic constitutive relations used to represent cartilage. Solid models of the relaxed tissue layers are assembled in physiological positions, resulting in a mathematical overlap of the cartilage layers. The overlap distribution is quantified and converted via the biphasic governing equations into applied traction boundary conditions for both the solid and fluid phases for each of the contacting layers. Linear, biphasic, three-dimensional, finite element analysis is performed using the contact boundary conditions derived for each of the contacting layers. The method is found to produce results consistent with the continuity requirements of biphasic contact. Comparison with results from independent, biphasic contact analyses of axisymmetric problems shows that the method slightly underestimates the contact area, leading to an overestimation of the total traction, but yields a good approximation to elastic stress and solid phase displacement.


2008 ◽  
Vol 36 (1) ◽  
pp. 63-79 ◽  
Author(s):  
L. Nasdala ◽  
Y. Wei ◽  
H. Rothert ◽  
M. Kaliske

Abstract It is a challenging task in the design of automobile tires to predict lifetime and performance on the basis of numerical simulations. Several factors have to be taken into account to correctly estimate the aging behavior. This paper focuses on oxygen reaction processes which, apart from mechanical and thermal aspects, effect the tire durability. The material parameters needed to describe the temperature-dependent oxygen diffusion and reaction processes are derived by means of the time–temperature–superposition principle from modulus profiling tests. These experiments are designed to examine the diffusion-limited oxidation (DLO) effect which occurs when accelerated aging tests are performed. For the cord-reinforced rubber composites, homogenization techniques are adopted to obtain effective material parameters (diffusivities and reaction constants). The selection and arrangement of rubber components influence the temperature distribution and the oxygen penetration depth which impact tire durability. The goal of this paper is to establish a finite element analysis based criterion to predict lifetime with respect to oxidative aging. The finite element analysis is carried out in three stages. First the heat generation rate distribution is calculated using a viscoelastic material model. Then the temperature distribution can be determined. In the third step we evaluate the oxygen distribution or rather the oxygen consumption rate, which is a measure for the tire lifetime. Thus, the aging behavior of different kinds of tires can be compared. Numerical examples show how diffusivities, reaction coefficients, and temperature influence the durability of different tire parts. It is found that due to the DLO effect, some interior parts may age slower even if the temperature is increased.


Sensors ◽  
2021 ◽  
Vol 21 (4) ◽  
pp. 1388
Author(s):  
Daniele Oboe ◽  
Luca Colombo ◽  
Claudio Sbarufatti ◽  
Marco Giglio

The inverse Finite Element Method (iFEM) is receiving more attention for shape sensing due to its independence from the material properties and the external load. However, a proper definition of the model geometry with its boundary conditions is required, together with the acquisition of the structure’s strain field with optimized sensor networks. The iFEM model definition is not trivial in the case of complex structures, in particular, if sensors are not applied on the whole structure allowing just a partial definition of the input strain field. To overcome this issue, this research proposes a simplified iFEM model in which the geometrical complexity is reduced and boundary conditions are tuned with the superimposition of the effects to behave as the real structure. The procedure is assessed for a complex aeronautical structure, where the reference displacement field is first computed in a numerical framework with input strains coming from a direct finite element analysis, confirming the effectiveness of the iFEM based on a simplified geometry. Finally, the model is fed with experimentally acquired strain measurements and the performance of the method is assessed in presence of a high level of uncertainty.


Author(s):  
Giuseppe Catania ◽  
Silvio Sorrentino

In the Rayleigh-Ritz condensation method the solution of the equation of motion is approximated by a linear combination of shape-functions selected among appropriate sets. Extensive literature dealing with the choice of appropriate basis of shape functions exists, the selection depending on the particular boundary conditions of the structure considered. This paper is aimed at investigating the possibility of adopting a set of eigenfunctions evaluated from a simple stucture as a general basis for the analysis of arbitrary-shaped plates. The results are compared to those available in the literature and using standard finite element analysis.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Charles Savoldelli ◽  
Elodie Ehrmann ◽  
Yannick Tillier

AbstractWith modern-day technical advances, high sagittal oblique osteotomy (HSOO) of the mandible was recently described as an alternative to bilateral sagittal split osteotomy for the correction of mandibular skeletal deformities. However, neither in vitro nor numerical biomechanical assessments have evaluated the performance of fixation methods in HSOO. The aim of this study was to compare the biomechanical characteristics and stress distribution in bone and osteosynthesis fixations when using different designs and placing configurations, in order to determine a favourable plating method. We established two finite element models of HSOO with advancement (T1) and set-back (T2) movements of the mandible. Six different configurations of fixation of the ramus, progressively loaded by a constant force, were assessed for each model. The von Mises stress distribution in fixations and in bone, and bony segment displacement, were analysed. The lowest mechanical stresses and minimal gradient of displacement between the proximal and distal bony segments were detected in the combined one-third anterior- and posterior-positioned double mini-plate T1 and T2 models. This suggests that the appropriate method to correct mandibular deformities in HSOO surgery is with use of double mini-plates positioned in the anterior one-third and posterior one-third between the bony segments of the ramus.


Materials ◽  
2020 ◽  
Vol 13 (15) ◽  
pp. 3280
Author(s):  
Raphaël Richert ◽  
Jean-Christophe Farges ◽  
Faleh Tamimi ◽  
Naim Naouar ◽  
Philippe Boisse ◽  
...  

Finite element (FE) models are widely used to investigate the biomechanics of reconstructed premolars. However, parameter identification is a complex step because experimental validation cannot always be conducted. The aim of this study was to collect the experimentally validated FE models of premolars, extract their parameters, and discuss trends. A systematic review was performed following Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. Records were identified in three electronic databases (MEDLINE [PubMed], Scopus, The Cochrane Library) by two independent reviewers. Twenty-seven parameters dealing with failure criteria, model construction, material laws, boundary conditions, and model validation were extracted from the included articles. From 1306 records, 214 were selected for eligibility and entirely read. Among them, 19 studies were included. A heterogeneity was observed for several parameters associated with failure criteria and model construction. Elasticity, linearity, and isotropy were more often chosen for dental and periodontal tissues with a Young’s modulus mostly set at 18–18.6 GPa for dentine. Loading was mainly simulated by an axial force, and FE models were mostly validated by in vitro tests evaluating tooth strains, but different conditions about experiment type, sample size, and tooth status (intact or restored) were reported. In conclusion, material laws identified herein could be applied to future premolar FE models. However, further investigations such as sensitivity analysis are required for several parameters to clarify their indication.


Sign in / Sign up

Export Citation Format

Share Document