Two-Way Coupled Multiscale Model for Predicting Mechanical Behavior of Bone Subjected to Viscoelastic Deformation and Fracture Damage

Author(s):  
Taesun You ◽  
Yong-Rak Kim ◽  
Taehyo Park

This paper presents a two-way linked computational multiscale model and its application to predict the mechanical behavior of bone subjected to viscoelastic deformation and fracture damage. The model is based on continuum thermos-mechanics and is implemented through the finite element method (FEM). Two physical length scales (the global scale of bone and local scale of compact bone) were two-way coupled in the framework by linking a homogenized global object to heterogeneous local-scale representative volume elements (RVEs). Multiscaling accounts for microstructure heterogeneity, viscoelastic deformation, and rate-dependent fracture damage at the local scale in order to predict the overall behavior of bone by using a viscoelastic cohesive zone model incorporated with a rate-dependent damage evolution law. In particular, age-related changes in material properties and geometries in bone were considered to investigate the effect of aging, loading rate, and damage evolution characteristics on the mechanical behavior of bone. The model successfully demonstrated its capability to predict the viscoelastic response and fracture damage due to different levels of aging, loading conditions (such as rates), and microscale damage evolution characteristics with only material properties of each constituent in the RVEs.

2021 ◽  
Vol 2021 ◽  
pp. 1-18
Author(s):  
Hao Liu ◽  
Lulin Zheng ◽  
Yujun Zuo ◽  
Zhonghu Wu ◽  
Wenjibin Sun ◽  
...  

The different directions of joints in rock will lead to great differences in damage evolution characteristics. This study utilizes DIP (digital image processing) technology for characterizing the mesostructure of sandstone and combines DIP technology with RFPA2D. The mesoscale fracture mechanics behavior of 7 groups of jointed sandstones with various dip angles was numerically studied, and its reliability was verified through theoretical analysis. According to digital image storage principle and box dimension theory, the box dimension algorithm of rock mesoscale fracture is written in MATLAB, the calculation method of fractal dimension of mesoscale fracture was proposed, and the corresponding relationship between mesoscale fractal dimension and fracture damage degree was established. Studies have shown that compressive strength as well as elastic modulus of sandstone leads to a U-shaped change when joint dip increases. There are a total of six final failure modes of joint samples with different inclination angles. Failure mode and damage degree can be quantified by D (fractal dimension) and ω (mesoscale fracture damage degree), respectively. The larger the ω, the more serious the damage, and the greater the D, the more complex the failure mode. Accumulative AE energy increases exponentially with the increase of loading step, and the growth process can be divided into gentle period, acceleration period, and surge period. The mesoscale fracture damage calculation based on the fractal dimension can be utilized for quantitatively evaluating the spatial distribution characteristics of mesoscale fracture, which provides a new way to study the law of rock damage evolution.


2020 ◽  
Vol 62 (6) ◽  
pp. 573-583
Author(s):  
Andreas Lutz ◽  
Lukas Huber ◽  
Claus Emmelmann

Author(s):  
Joseph Park ◽  
Andrew Shin ◽  
Somaye Jafari ◽  
Joseph L. Demer

AbstractThe optic nerve (ON) is a recently recognized tractional load on the eye during larger horizontal eye rotations. In order to understand the mechanical behavior of the eye during adduction, it is necessary to characterize material properties of the sclera, ON, and in particular its sheath. We performed tensile loading of specimens taken from fresh postmortem human eyes to characterize the range of variation in their biomechanical properties and determine the effect of preconditioning. We fitted reduced polynomial hyperelastic models to represent the nonlinear tensile behavior of the anterior, equatorial, posterior, and peripapillary sclera, as well as the ON and its sheath. For comparison, we analyzed tangent moduli in low and high strain regions to represent stiffness. Scleral stiffness generally decreased from anterior to posterior ocular regions. The ON had the lowest tangent modulus, but was surrounded by a much stiffer sheath. The low-strain hyperelastic behaviors of adjacent anatomical regions of the ON, ON sheath, and posterior sclera were similar as appropriate to avoid discontinuities at their boundaries. Regional stiffnesses within individual eyes were moderately correlated, implying that mechanical properties in one region of an eye do not reliably reflect properties of another region of that eye, and that potentially pathological combinations could occur in an eye if regional properties are discrepant. Preconditioning modestly stiffened ocular tissues, except peripapillary sclera that softened. The nonlinear mechanical behavior of posterior ocular tissues permits their stresses to match closely at low strains, although progressively increasing strain causes particularly great stress in the peripapillary region.


Energies ◽  
2021 ◽  
Vol 14 (15) ◽  
pp. 4585
Author(s):  
Marian Bulla ◽  
Stefan Kolling ◽  
Elham Sahraei

The present study is focused on the development of a material model where the orthotropic-visco-elastic and orthotropic-visco-plastic mechanical behavior of a polymeric material is considered. The increasing need to reduce the climate-damaging exhaust gases in the automotive industry leads to an increasing usage of electric powered drive systems using Lithium-ion (Li-ion) batteries. For the safety and crashworthiness investigations, a deeper understanding of the mechanical behavior under high and dynamic loads is needed. In order to prevent internal short circuits and thermal runaways within a Li-ion battery, the separator plays a crucial role. Based on results of material tests, a novel material model for finite element analysis (FEA) is developed using the explicit solver Altair Radioss. Based on this model, the visco-elastic-orthotropic, as well as the visco-plastic-orthotropic, behavior until failure can be modeled. Finally, a FE simulation model of the separator material is performed, using the results of different tensile tests conducted at three different velocities, 0.1 mm·s−1, 1.0 mm·s−1 and 10.0 mm·s−1 and different orientations of the specimen. The purpose is to predict the anisotropic, rate-dependent stiffness behavior of separator materials in order to improve FE simulations of the mechanical behavior of batteries and therefore reduce the development time of electrically powered vehicles and consumer goods. The present novel material model in combination with a well-suited failure criterion, which considers the different states of stress and anisotropic-visco-dependent failure limits, can be applied for crashworthiness FE analysis. The model succeeded in predicting anisotropic, visco-elastic orthotropic and visco-plastic orthotropic stiffness behavior up to failure.


2005 ◽  
Vol 127 (7) ◽  
pp. 1158-1167 ◽  
Author(s):  
Adam H. Hsieh ◽  
Diane R. Wagner ◽  
Louis Y. Cheng ◽  
Jeffrey C. Lotz

In vivo rodent tail models are becoming more widely used for exploring the role of mechanical loading on the initiation and progression of intervertebral disc degeneration. Historically, finite element models (FEMs) have been useful for predicting disc mechanics in humans. However, differences in geometry and tissue properties may limit the predictive utility of these models for rodent discs. Clearly, models that are specific for rodent tail discs and accurately simulate the disc’s transient mechanical behavior would serve as important tools for clarifying disc mechanics in these animal models. An FEM was developed based on the structure, geometry, and scale of the mouse tail disc. Importantly, two sources of time-dependent mechanical behavior were incorporated: viscoelasticity of the matrix, and fluid permeation. In addition, a novel strain-dependent swelling pressure was implemented through the introduction of a dilatational stress in nuclear elements. The model was then validated against data from quasi-static tension-compression and compressive creep experiments performed previously using mouse tail discs. Finally, sensitivity analyses were performed in which material parameters of each disc subregion were individually varied. During disc compression, matrix consolidation was observed to occur preferentially at the periphery of the nucleus pulposus. Sensitivity analyses revealed that disc mechanics was greatly influenced by changes in nucleus pulposus material properties, but rather insensitive to variations in any of the endplate properties. Moreover, three key features of the model—nuclear swelling pressure, lamellar collagen viscoelasticity, and interstitial fluid permeation—were found to be critical for accurate simulation of disc mechanics. In particular, collagen viscoelasticity dominated the transient behavior of the disc during the initial 2200s of creep loading, while fluid permeation governed disc deformation thereafter. The FEM developed in this study exhibited excellent agreement with transient creep behavior of intact mouse tail motion segments. Notably, the model was able to produce spatial variations in nucleus pulposus matrix consolidation that are consistent with previous observations in nuclear cell morphology made in mouse discs using confocal microscopy. Results of this study emphasize the need for including nucleus swelling pressure, collagen viscoelasticity, and fluid permeation when simulating transient changes in matrix and fluid stress/strain. Sensitivity analyses suggest that further characterization of nucleus pulposus material properties should be pursued, due to its significance in steady-state and transient disc mechanical response.


2018 ◽  
Vol 35 (4) ◽  
pp. 441-454 ◽  
Author(s):  
M. Shishesaz ◽  
M. Hosseini

ABSTRACTIn this paper, the mechanical behavior of a functionally graded nano-cylinder under a radial pressure is investigated. Strain gradient theory is used to include the small scale effects in this analysis. The variations in material properties along the thickness direction are included based on three different models. Due to slight variations in engineering materials, the Poisson’s ratio is assumed to be constant. The governing equation and its corresponding boundary conditions are obtained using Hamilton’s principle. Due to the complexity of the governed system of differential equations, numerical methods are employed to achieve a solution. The analysis is general and can be reduced to classical elasticity if the material length scale parameters are taken to be zero. The effect of material indexn, variations in material properties and the applied internal and external pressures on the total and high-order stresses, are well examined. For the cases in which the applied external pressure at the inside (or outside) radius is zero, due to small effects in nano-cylinder, some components of the high-order radial stresses do not vanish at the boundaries. Based on the results, the material inhomogeneity indexn, as well as the selected model through which the mechanical properties may vary along the thickness, have significant effects on the radial and circumferential stresses.


Sign in / Sign up

Export Citation Format

Share Document