Bone Remodeling Algorithm Incorporating Various Quantities as Mechanical Stimulus and Assuming Initial Microcrack in Bone

2017 ◽  
Vol 754 ◽  
pp. 189-193
Author(s):  
Libor Borák ◽  
Petr Marcián

It is widely accepted that bones have the ability to adapt to new biomechanical environment by changing their material properties, geometry and inner architecture. Bones have also an exceptional ability to self-repair, to remove microcracks and to prevent the bone damage caused by the fatigue failure. These abilities are enabled through coupled processes of bone resorption and bone formation, the processes collectively referred to as bone remodeling. Numerous studies have shown that bone remodeling is governed by combination of mechanical stimulus (strains) and its frequency, both sensed by sensor cells (osteocytes). Through mechanotransduction, the stimulus is transmitted to actor cells (osteoclasts, osteoblasts) that actually do the bone resorption or formation. Several theories have been proposed to predict bone remodeling and several finite-element-based algorithms have been introduced. The vast majority of them uses strain energy density as the mechanical stimulus. The purpose of this paper is to investigate and discuss the applicability of also other strain-based representations of the mechanical stimulus in simulations of remodeling of bone with an initial microcrack. The need for developing more reliable models is essential for both clinicians and engineers who are interested, for instance, in prediction of bone performance when various implants are involved.

1996 ◽  
Vol 436 ◽  
Author(s):  
R. P. Vinci ◽  
J. C. Bravman

AbstractWe have modeled the effects of grain aspect ratio on strain energy density in (100)-oriented grains in a (111)-textured Cu film on a Si substrate. Minimization of surface energy, interface energy, and strain energy density (SED) drives preferential growth of grains of certain crystallographic orientations in thin films. Under conditions in which the SED driving force exceeds the surface- and interface-energy driving forces, Cu films develop abnormally large (100) oriented grains during annealing. In the elastic regime the SED differences between the (100) grains and the film average arise from elastic anisotropy. Previous analyses indicate that several factors (e.g. elimination of grain boundaries during grain growth) may alter the magnitude of the SED driving force. We demonstrate, using finite element modeling of a single columnar (100) grain in a (111) film, that changes in grain aspect ratio can significantly affect the SED driving force. A minimum SED driving force is found for (100) Cu grains with diameters on the order of the film thickness. In the absence of other stagnation mechanisms, such behavior could cause small grains to grow abnormally and then stagnate while large grains continue to grow. This would lead to a bimodal grain size distribution in the (100) grains preferred by the SED minimization.


1971 ◽  
Vol 93 (2) ◽  
pp. 138-145 ◽  
Author(s):  
B. R. Simon ◽  
A. S. Kobayashi ◽  
D. E. Strandness ◽  
C. A. Wiederhielm

Possible relations between arterial wall stresses and deformations and mechanisms contributing to atherosclerosis are discussed. Necessary material properties are determined experimentally and from available data in the literature by assuming the arterial response to be a static finite deformation of a thick-walled cylinder constrained in a state of plane strain and composed of an incompressible, nonlinear elastic, transversely isotropic material. Experimental justification from the literature and supporting theoretical considerations are presented for each assumption. The partial derivative of the strain energy density function δW1/δI , necessary for in-plane stress calculation, is determined to be of exponential form using in situ biaxial test results from the canine abdominal aorta. An axisymmetric numerical integration solution is developed and used as a check for finite element results. The large deformation finite element theory of Oden is modified to include aortic material nonlinearity and directional properties and is used for a structural analysis of the aortic cross section. Results of this investigation are: (a) Fung’s exponential form for the strain energy density function of soft tissues is found to be valid for the aorta in the biaxial states considered; (b) finite deformation analyses by the finite element method and numerical integration solution reveal that significant tangential stress gradients are present in arteries commonly assumed to be “thin-walled” tubes using linear theory.


2002 ◽  
Vol 124 (4) ◽  
pp. 403-410 ◽  
Author(s):  
J. Lau ◽  
Z. Mei ◽  
S. Pang ◽  
C. Amsden ◽  
J. Rayner ◽  
...  

Thermal reliability of the solder sealing ring of Agilent Technologies’ bubble-actuated photonic cross-connect switches has been investigated in this paper. Emphasis is placed on the determination of the thermal-fatigue life of the solder sealing ring under shipping/storing/handling conditions. The solder ring is assumed to obey the Garofalo-Arrhenius creep constitutive law. The nonlinear responses such as the deflections, stresses, creep strains, and creep strain energy density of the 3-D photonic package have been determined with a commercial finite element code. In addition, isothermal fatigue tests have been performed to obtain the relationship between the number of cycle-to-failure and the strain energy density. Thus, by combining the finite element results and the isothermal fatigue test results, the average thermal-fatigue life of the solder sealing ring is readily determined and is found to be more than adequate for shipping/storing/handling the photonic switches.


2008 ◽  
Vol 131 (1) ◽  
Author(s):  
In Gwun Jang ◽  
Il Yong Kim ◽  
Byung Man Kwak

In bone-remodeling studies, it is believed that the morphology of bone is affected by its internal mechanical loads. From the 1970s, high computing power enabled quantitative studies in the simulation of bone remodeling or bone adaptation. Among them, Huiskes et al. (1987, “Adaptive Bone Remodeling Theory Applied to Prosthetic Design Analysis,” J. Biomech. Eng., 20, pp. 1135–1150) proposed a strain energy density based approach to bone remodeling and used the apparent density for the characterization of internal bone morphology. The fundamental idea was that bone density would increase when strain (or strain energy density) is higher than a certain value and bone resorption would occur when the strain (or strain energy density) quantities are lower than the threshold. Several advanced algorithms were developed based on these studies in an attempt to more accurately simulate physiological bone-remodeling processes. As another approach, topology optimization originally devised in structural optimization has been also used in the computational simulation of the bone-remodeling process. The topology optimization method systematically and iteratively distributes material in a design domain, determining an optimal structure that minimizes an objective function. In this paper, we compared two seemingly different approaches in different fields—the strain energy density based bone-remodeling algorithm (biomechanical approach) and the compliance based structural topology optimization method (mechanical approach)—in terms of mathematical formulations, numerical difficulties, and behavior of their numerical solutions. Two numerical case studies were conducted to demonstrate their similarity and difference, and then the solution convergences were discussed quantitatively.


Author(s):  
Michał Nowak ◽  
Jan Sokołowski ◽  
Antoni Żochowski

The trabecular bone can adapt its form to mechanical loads and form structures that are both lightweight and very stiff. In this sense, it is a problem similar to structural optimization, especially topology optimization. The natural phenomenon leading to mechanical optimization of the bone structures is called trabecular bone remodeling. The main assumption and the benchmark for the numerical models of the phenomenon is the observation that the strain energy density on the structural surface is constant. This constant value corresponds to the homeostatic strain energy density, the state of bone tissue with a perfect balance of the loss, and gain of the bone mass. We assumed that the trabecular bone can form an optimal structure. The idea behind the investigation is to carry out studies on the ground of mechanics and to interpret clinical observations in the context of the results obtained from the optimization studies. In this way, clinical observations have been confirmed by strict arguments based on mechanics, leading to the unequivocal conclusion that equalization of the strain energy density on the trabecular bone surface allows minimizing the strain energy in the whole structure of the bone. This proves the veracity of the assumption that the remodeling process leads to the formation of the structure with the highest stiffness. In addition, this article elaborates on two new aspects of the remodeling phenomenon resulting directly from the considerations in the field of shape optimization important for numerical simulation. The first one concerns the influence of surface curvature on the remodeling process. The second one concerns the role of the bone surface where different loads are analyzed. Both aspects show the need of actual trabecular bone geometry model for the simulation of the trabecular bone remodeling phenomenon.


Sign in / Sign up

Export Citation Format

Share Document