Comparing Predictive Accuracy and Computational Costs for Viscoelastic Modeling of Spinal Cord Tissues

2019 ◽  
Vol 141 (5) ◽  
Author(s):  
Nicole L. Ramo ◽  
Kevin L. Troyer ◽  
Christian M. Puttlitz

Abstract The constitutive equation used to characterize and model spinal tissues can significantly influence the conclusions from experimental and computational studies. Therefore, researchers must make critical judgments regarding the balance of computational efficiency and predictive accuracy necessary for their purposes. The objective of this study is to quantitatively compare the fitting and prediction accuracy of linear viscoelastic (LV), quasi-linear viscoelastic (QLV), and (fully) nonlinear viscoelastic (NLV) modeling of spinal-cord-pia-arachnoid-construct (SCPC), isolated cord parenchyma, and isolated pia-arachnoid-complex (PAC) mechanics in order to better inform these judgements. Experimental data collected during dynamic cyclic testing of each tissue condition were used to fit each viscoelastic formulation. These fitted models were then used to predict independent experimental data from stress-relaxation testing. Relative fitting accuracy was found not to directly reflect relative predictive accuracy, emphasizing the need for material model validation through predictions of independent data. For the SCPC and isolated cord, the NLV formulation best predicted the mechanical response to arbitrary loading conditions, but required significantly greater computational run time. The mechanical response of the PAC under arbitrary loading conditions was best predicted by the QLV formulation.

2022 ◽  
Vol 119 (3) ◽  
pp. e2117232119
Author(s):  
Giulio Franchini ◽  
Ivan D. Breslavsky ◽  
Francesco Giovanniello ◽  
Ali Kassab ◽  
Gerhard A. Holzapfel ◽  
...  

Experimental data and a suitable material model for human aortas with smooth muscle activation are not available in the literature despite the need for developing advanced grafts; the present study closes this gap. Mechanical characterization of human descending thoracic aortas was performed with and without vascular smooth muscle (VSM) activation. Specimens were taken from 13 heart-beating donors. The aortic segments were cooled in Belzer UW solution during transport and tested within a few hours after explantation. VSM activation was achieved through the use of potassium depolarization and noradrenaline as vasoactive agents. In addition to isometric activation experiments, the quasistatic passive and active stress–strain curves were obtained for circumferential and longitudinal strips of the aortic material. This characterization made it possible to create an original mechanical model of the active aortic material that accurately fits the experimental data. The dynamic mechanical characterization was executed using cyclic strain at different frequencies of physiological interest. An initial prestretch, which corresponded to the physiological conditions, was applied before cyclic loading. Dynamic tests made it possible to identify the differences in the viscoelastic behavior of the passive and active tissue. This work illustrates the importance of VSM activation for the static and dynamic mechanical response of human aortas. Most importantly, this study provides material data and a material model for the development of a future generation of active aortic grafts that mimic natural behavior and help regulate blood pressure.


Polymers ◽  
2019 ◽  
Vol 11 (4) ◽  
pp. 581 ◽  
Author(s):  
Imran Syed ◽  
Pascal Vouagner ◽  
Frank Fleck ◽  
Jorge Lacayo-Pineda

Nonlinear material response is analysed with the Fourier transform (FT) of the raw signal measured by a high-frequency dynamic mechanical analyzer (HF DMA). It is known from rheological behaviour of elastomers that reinforcing fillers additionally induce nonlinearity in an already inherently nonlinear system. This behaviour is often described in terms of a mechanical response of strain sweeps, essentially the transition from the linear viscoelastic (LVE) to the nonlinear viscoelastic (NVE) region. In the current investigation, the NVE region could be observed with respect to frequency under low-amplitude deformation. A foldover effect was observed, whereby the material exhibited a nonlinear dependency in relation to the increment of the filler amount above the percolation threshold. In addition, an apparent superharmonic resonance was observed within higher orders of vibrational modes which is further indication of nonlinearity. In this paper, the analytical approach is presented as a novel method to characterise the behaviour of the polymer–filler interaction by HF DMA.


Author(s):  
Kaveh Laksari ◽  
Mehdi Shafieian ◽  
Kurosh Darvish ◽  
Keyanoush Sadeghipour

This study investigates the propagation of shock waves and self-preserving waves in soft tissues such as brain as a mechanism of injury in high rate loading conditions as seen in blast-induced neurotrauma (BINT). The derived mathematical models indicate that whereas linear viscoelastic models predict only decaying waves, instances of such phenomena as shock can be achieved in nonlinear media. In this study, a nonlinear viscoelastic material model for brain tissue was developed in compression. Furthermore, nonlinear viscoelastic wave propagation in brain tissue was studied and a criterion for the development of shock waves was formulated. It was shown that discontinuities in the acceleration that happen in blast loading conditions may evolve to shock waves, resulting in large discontinuities in strain and stress at the wave front leading to tissue injuries.


2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Amin Komeili ◽  
Akbar Rasoulian ◽  
Fatemeh Moghaddam ◽  
Marwan El-Rich ◽  
Le Ping Li

Abstract Background Linear elastic, hyperelastic, and multiphasic material constitutive models are frequently used for spinal intervertebral disc simulations. While the characteristics of each model are known, their effect on spine mechanical response requires a careful investigation. The use of advanced material models may not be applicable when material constants are not available, model convergence is unlikely, and computational time is a concern. On the other hand, poor estimations of tissue’s mechanical response are likely if the spine model is oversimplified. In this study, discrepancies in load response introduced by material models will be investigated. Methods Three fiber-reinforced C2-C3 disc models were developed with linear elastic, hyperelastic, and biphasic behaviors. Three different loading modes were investigated: compression, flexion and extension in quasi-static and dynamic conditions. The deformed disc height, disc fluid pressure, range of motion, and stresses were compared. Results Results indicated that the intervertebral disc material model has a strong effect on load-sharing and disc height change when compression and flexion were applied. The predicted mechanical response of three models under extension had less discrepancy than its counterparts under flexion and compression. The fluid-solid interaction showed more relevance in dynamic than quasi-static loading conditions. The fiber-reinforced linear elastic and hyperelastic material models underestimated the load-sharing of the intervertebral disc annular collagen fibers. Conclusion This study confirmed the central role of the disc fluid pressure in spinal load-sharing and highlighted loading conditions where linear elastic and hyperelastic models predicted energy distribution different than that of the biphasic model.


2011 ◽  
Vol 138-139 ◽  
pp. 466-470
Author(s):  
Long Yi ◽  
Yun Peng ◽  
Hou Quan Hong ◽  
Yu Liang Li

Based on the subroutine VUMAT, user-defined material model in the nonlinear finite element software ABAQUS/EXPLICIT, a nonlinear viscoelastic constitutive model is developed. The validify of the subroutine has been proven through the standard uniaxial tensile model. The shortage of finite element softwares which only have linear viscoelastic constitutive model is remedied. This paper presents the process of developing a material constitutive model and some useful technology. It can be referred for extending the material constitutive model in finite element softwares.


2021 ◽  
Vol 8 ◽  
Author(s):  
Shaopu Yang ◽  
Peng Wang ◽  
Yongqiang Liu ◽  
Xufeng Dong ◽  
Yu Tong ◽  
...  

To accurately characterize the mechanical behavior of magnetorheological elastomer (MRE) under a wide range of strain amplitude, excitation frequency, and magnetic field, the viscoelastic fractional derivative was introduced, and a modified Bouc-Wen model based on fractional derivative for MRE in a nonlinear viscoelastic region was established. The Bouc-Wen model can accurately describe the hysteretic characteristics of the MRE nonlinear viscoelastic region, but it cannot accurately simulate magneto-viscoelasticity and frequency dependence. The fractional derivative can express this characteristic with fewer parameters and higher accuracy. The model’s validity was verified by fitting the experimental data of stress and strain measured in shear mode. By analyzing the coupling relationship between the model parameters and strain amplitude, frequency, and magnetic flux densities, a method of parameter identification under multi-loading conditions was proposed, and the modified model parameters were identified. The results reveal that the modified Bouc-Wen model can accurately characterize the mechanical properties of the nonlinear viscoelastic region of MRE, and the fitting accuracy is significantly improved compared with the Bouc-Wen model. The expression of the model parameters obtained from the method of parameter identification under multi-loading conditions is accurate in a wide range of strain amplitude, frequency, and magnetic flux density. The fitness values of simulation data and experimental data under identified and non-identified conditions exceed 90%, showcasing the effectiveness of the modified Bouc-Wen model and the feasibility of the parameter identification method under multi-loading conditions.


2008 ◽  
Vol 36 (3) ◽  
pp. 211-226 ◽  
Author(s):  
F. Liu ◽  
M. P. F. Sutcliffe ◽  
W. R. Graham

Abstract In an effort to understand the dynamic hub forces on road vehicles, an advanced free-rolling tire-model is being developed in which the tread blocks and tire belt are modeled separately. This paper presents the interim results for the tread block modeling. The finite element code ABAQUS/Explicit is used to predict the contact forces on the tread blocks based on a linear viscoelastic material model. Special attention is paid to investigating the forces on the tread blocks during the impact and release motions. A pressure and slip-rate-dependent frictional law is applied in the analysis. A simplified numerical model is also proposed where the tread blocks are discretized into linear viscoelastic spring elements. The results from both models are validated via experiments in a high-speed rolling test rig and found to be in good agreement.


2018 ◽  
Vol 84 (10) ◽  
pp. 23-28
Author(s):  
D. A. Golentsov ◽  
A. G. Gulin ◽  
Vladimir A. Likhter ◽  
K. E. Ulybyshev

Destruction of bodies is accompanied by formation of both large and microscopic fragments. Numerous experiments on the rupture of different samples show that those fragments carry a positive electric charge. his phenomenon is of interest from the viewpoint of its potential application to contactless diagnostics of the early stage of destruction of the elements in various technical devices. However, the lack of understanding the nature of this phenomenon restricts the possibility of its practical applications. Experimental studies were carried out using an apparatus that allowed direct measurements of the total charge of the microparticles formed upon sample rupture and determination of their size and quantity. The results of rupture tests of duralumin and electrical steel showed that the size of microparticles is several tens of microns, the particle charge per particle is on the order of 10–14 C, and their amount can be estimated as the ratio of the cross-sectional area of the sample at the point of discontinuity to the square of the microparticle size. A model of charge formation on the microparticles is developed proceeding from the experimental data and current concept of the electron gas in metals. The model makes it possible to determine the charge of the microparticle using data on the particle size and mechanical and electrical properties of the material. Model estimates of the total charge of particles show order-of-magnitude agreement with the experimental data.


2021 ◽  
Vol 11 (12) ◽  
pp. 5509
Author(s):  
Hongjin Choi ◽  
Seonghwan Choi ◽  
Soo-Chang Kang ◽  
Myoung-Gyu Lee

A fully implicit stress integration algorithm is developed for the distortional hardening model, namely the e−HAH model, capable of simulating cross−hardening/softening under orthogonal loading path changes. The implicit algorithm solves a complete set of residuals as nonlinear functions of stress, a microstructure deviator, and plastic state variables of the constitutive model, and provides a consistent tangent modulus. The number of residuals is set to be 20 or 14 for the continuum or shell elements, respectively. Comprehensive comparison programs are presented regarding the predictive accuracy and stability with different numerical algorithms, strain increments, material properties, and loading conditions. The flow stress and r−value evolutions under reverse/cross−loading conditions prove that the algorithm is robust and accurate, even with large strain increments. By contrast, the cutting−plane method and partially implicit Euler backward method, which are characterized by a reduced number of residuals, result in unstable responses under abrupt loading path changes. Finally, the algorithm is implemented into the finite element modeling of large−size, S−rail forming and the springback for two automotive steel sheets, which is often solved by a hybrid dynamic explicit–implicit scheme. The fully implicit algorithm performs well for the whole simulation with the solely static implicit scheme.


Sign in / Sign up

Export Citation Format

Share Document