Flexible multibody approach in bone strain estimation during physical activity:

Bone strain plays a major role as the activation signal for the bone (re)modeling process, which is vital for keeping bones healthy. Maintaining high bone mineral density reduces the chances of fracture in the event of an accident. Numerous studies have shown that bones can be strengthened with physical exercise. Several hypotheses have asserted that a stronger osteogenic (bone producing) effect results from dynamic exercise than from static exercise. These previous studies are based on short-term empirical research, which provide the motivation for justifying the experimental results with a solid mathematical background. The computer simulation techniques utilized in this work allow for non-invasive bone strain estimation during physical activity at any bone site within the human skeleton. All models presented in the study are threedimensional and actuated by muscle models to replicate the real conditions accurately. The objective of this work is to determine and present loading-induced bone strain values resulting from physical activity. It includes a comparison of strain resulting from four different gym exercises (knee flexion, knee extension, leg press, and squat) and walking, with the results reported for walking and jogging obtained from in-vivo measurements described in the literature. The objective is realized primarily by carrying out flexible multibody dynamics computer simulations. The dissertation combines the knowledge of finite element analysis and multibody simulations with experimental data and information available from medical field literature. Measured subject-specific motion data was coupled with forward dynamics simulation to provide natural skeletal movement. Bone geometries were defined using a reverse engineering approach based on medical imaging techniques. Both computed tomography and magnetic resonance imaging were utilized to explore modeling differences. The predicted tibia bone strains during walking show good agreement with invivo studies found in the literature. Strain measurements were not available for gym exercises; therefore, the strain results could not be validated. However, the values seem reasonable when compared to available walking and running invivo strain measurements. The results can be used for exercise equipment design aimed at strengthening the bones as well as the muscles during workout. Clinical applications in post fracture recovery exercising programs could also be the target. In addition, the methodology introduced in this study, can be applied to investigate the effect of weightlessness on astronauts, who often suffer bone loss after long time spent in the outer space.

Comparison of beam theory and finite-element analysis with in vivo bone strain data from the alligator cranium

The mechanical behavior of the vertebrate skull is often modeled using free-body analysis of simple geometric structures and, more recently, finite-element (FE) analysis. In this study, we compare experimentally collected in vivo bone strain orientations and magnitudes from the cranium of the American alligator with those extrapolated from a beam model and extracted from an FE model. The strain magnitudes predicted from beam and FE skull models bear little similarity to relative and absolute strain magnitudes recorded during in vivo biting experiments. However, quantitative differences between principal strain orientations extracted from the FE skull model and recorded during the in vivo experiments were smaller, and both generally matched expectations from the beam model. The differences in strain magnitude between the data sets may be attributable to the level of resolution of the models, the material properties used in the FE model, and the loading conditions (i.e., external forces and constraints). This study indicates that FE models and modeling of skulls as simple engineering structures may give a preliminary idea of how these structures are loaded, but whenever possible, modeling results should be verified with either in vitro or preferably in vivo testing, especially if precise knowledge of strain magnitudes is desired. (c) 2005 Wiley-Liss, Inc.

Bone stress and strain after use of a miniplate for molar protraction and uprighting: a 3-dimensional finite element analysis

The aim of this study was to use the finite element method to evaluate the distribution of stresses and strains on the local bone tissue adjacent to the miniplate used for anchorage of orthodontic forces. Methods: A 3-dimensional model composed of a hemimandible and teeth was constructed using dental computed tomographic images, in which we assembled a miniplate with fixation screws. The uprighting and mesial movements of the mandibular second molar that was anchored with the miniplate were simulated. The miniplate was loaded with horizontal forces of 2, 5, and 15 N. A moment of 11.77 N.mm was also applied. The stress and strain distributions were analyzed, and their correlations with the bone remodeling criteria and miniplate stability were assessed. Results: When orthodontic loads were applied, peak bone strain remained within the range of bone homeostasis (100-1500 mu m strain) with a balance between bone formation and resorption. The maximum deformation was found to be 1035 mu m strain with a force of 5 N. At a force of 15 N, bone resorption was observed in the region of the screws. Conclusions: We observed more stress concentration around the screws than in the cancellous bone. The levels of stress and strain increased when the force was increased but remained within physiologic levels. The anchorage system of miniplate and screws could withstand the orthodontic forces, which did not affect the stability of the miniplate.

Flexible Multibody Simulation Approach in the Dynamic Analysis of Bone Strains Activity

The objective of this study is to show that bone strains due to dynamic mechanical loading during physical activity can be analysed using the flexible multibody simulation approach. Strains within the bone tissue play a major role in bone (re)modeling. Based on previous studies, it has been shown that dynamic loading seems to be more important for bone (re)modeling than static loading. The finite element method has been used previously to assess bone strains. However, the finite element method may be limited to static analysis of bone strains due to the expensive computation required for dynamic analysis, especially for a biomechanical system consisting of several bodies. Further, in vivo implementation of strain gauges on the surfaces of bone has been used previously in order to quantify the mechanical loading environment of the skeleton. However, in vivo strain measurement requires invasive methodology, which is challenging and limited to certain regions of superficial bones only, such as the anterior surface of the tibia. In this study, an alternative numerical approach to analyzing in vivo strains, based on the flexible multibody simulation approach, is proposed. In order to investigate the reliability of the proposed approach, three 3-dimensional musculoskeletal models where the right tibia is assumed to be flexible, are used as demonstration examples. The models are employed in a forward dynamics simulation in order to predict the tibial strains during walking on a level exercise. The flexible tibial model is developed using the actual geometry of the subject’s tibia, which is obtained from 3 dimensional reconstruction of Magnetic Resonance Images. Inverse dynamics simulation based on motion capture data obtained from walking at a constant velocity is used to calculate the desired contraction trajectory for each muscle. In the forward dynamics simulation, a proportional derivative servo controller is used to calculate each muscle force required to reproduce the motion, based on the desired muscle contraction trajectory obtained from the inverse dynamics simulation. Experimental measurements are used to verify the models and check the accuracy of the models in replicating the realistic mechanical loading environment measured from the walking test. The predicted strain results by the models show consistency with literature-based in vivo strain measurements. In conclusion, the non-invasive flexible multibody simulation approach may be used as a surrogate for experimental bone strain measurement, and thus be of use in detailed strain estimation of bones in different applications. Consequently, the information obtained from the present approach might be useful in clinical applications, including optimizing implant design and devising exercises to prevent bone fragility, accelerate fracture healing and reduce osteoporotic bone loss.

The Dynamic Analysis of a Human Skeleton Using the Flexile Multibody Simulation Approach

The human motion study, which relies on mathematical and computational models ingeneral, and multibody dynamic biomechanical models in particular, has become asubject of many recent researches. The human body model can be applied to different physical exercises and many important results such as muscle forces, which are difficult to be measured through practical experiments, can be obtained easily. In the work, human skeletal lower limb model consisting of three bodies in build using the flexible multibody dynamics simulation approach. The floating frame of reference formulation is used to account for the flexibility in the bones of the human lower limb model. The main reason of considering the flexibility inthe human bones is to measure the strains in the bone result from different physical exercises. It has been perceived the bone under strain will become stronger in order to cope with the exercise. On the other hand, the bone strength is considered and important factors in reducing the bone fractures. The simulation approach and model developed in this work are used to measure the bone strain results from applying raising the sole of the foot exercise. The simulation results are compared to the results available in literature. The comparison shows goof agreement. This study sheds the light on the importance of using the flexible multibody dynamic simulation approach to build human biomechanical models, which can be used in developing some exercises to achieve the optimalbone strength.

Time series analysis of jaw muscle contraction and tissue deformation during mastication in miniature pigs

Masticatory muscle contraction causes both jaw movement and tissue deformation during function. Natural chewing data from 25 adult miniature pigs were studied by means of time series analysis. The data set included simultaneous recordings of electromyography (EMG) from bilateral masseter (MA), zygomaticomandibularis (ZM) and lateral pterygoid muscles, bone surface strains from the left squamosal bone (SQ), condylar neck (CD) and mandibular corpus (MD), and linear deformation of the capsule of the jaw joint measured bilaterally using differential variable reluctance transducers. Pairwise comparisons were examined by calculating the cross-correlation functions. Jaw-adductor muscle activity of MA and ZM was found to be highly cross-correlated with CD and SQ strains and weakly with MD strain. No muscle’s activity was strongly linked to capsular deformation of the jaw joint, nor were bone strains and capsular deformation tightly linked. Homologous muscle pairs showed the greatest synchronization of signals, but the signals themselves were not significantly more correlated than those of non-homologous muscle pairs. These results suggested that bone strains and capsular deformation are driven by different mechanical regimes. Muscle contraction and ensuing reaction forces are probably responsible for bone strains, whereas capsular deformation is more likely a product of movement.

Damage rate is a predictor of fatigue life and creep strain rate in fatigue of tensile human cortical bone samples

We present results on the growth of damage in 29 fatigue tests of human femoral cortical bone from four individuals, aged 53–79. In these tests we examine the interdependency of stress, cycles to failure, rate of creep strain, and rate of modulus loss. The behavior of creep rates has been reported recently for the same donors as an effect of stress and cycles (Cotton, J. R., Zioupos, P., Winwood, K., and Taylor, M., 2003, "Analysis of Creep Strain During Tensile Fatigue of Cortical Bone," J. Biomech. 36, pp. 943–949). In the present paper we first examine how the evolution of damage (drop in modulus per cycle) is associated with the stress level or the "normalized stress" level (stress divided by specimen modulus), and results show the rate of modulus loss fits better as a function of normalized stress. However, we find here that even better correlations can be established between either the cycles to failure or creep rates versus rates of damage than any of these three measures versus normalized stress. The data indicate that damage rates can be excellent predictors of fatigue life and creep strain rates in tensile fatigue of human cortical bone for use in practical problems and computer simulations.

Analysis of creep strain during tensile fatigue of cortical bone

During fatigue tests of cortical bone specimens, at the unload portion of the cycle (zero stress) non-zero strains occur and progressively accumulate as the test progresses. This non-zero strain is hypothesised to be mostly, if not entirely, describable as creep. This work examines the rate of accumulation of this strain and quantifies its stress dependency. A published relationship determined from creep tests of cortical bone (Journal of Biomechanics 21 (1988) 623) is combined with knowledge of the stress history during fatigue testing to derive an expression for the amount of creep strain in fatigue tests. Fatigue tests on 31 bone samples from four individuals showed strong correlations between creep strain rate and both stress and “normalised stress” (σ/E) during tensile fatigue testing (0–T). Combined results were good (r2=0.78) and differences between the various individuals, in particular, vanished when effects were examined against normalised stress values. Constants of the regression showed equivalence to constants derived in creep tests. The universality of the results, with respect to four different individuals of both sexes, shows great promise for use in computational models of fatigue in bone structures.

Numerical simulation of strain-adaptive bone remodelling in the ankle joint

Background: The use of artificial endoprostheses has become a routine procedure for knee and hip joints while ankle arthritis has traditionally been treated by means of arthrodesis. Due to its advantages, the implantation of endoprostheses is constantly increasing. While finite element analyses (FEA) of strain-adaptive bone remodelling have been carried out for the hip joint in previous studies, to our knowledge there are no investigations that have considered remodelling processes of the ankle joint. In order to evaluate and optimise new generation implants of the ankle joint, as well as to gain additional knowledge regarding the biomechanics, strain-adaptive bone remodelling has been calculated separately for the tibia and the talus after providing them with an implant. Methods: FE models of the bone-implant assembly for both the tibia and the talus have been developed. Bone characteristics such as the density distribution have been applied corresponding to CT scans. A force of 5,200 N, which corresponds to the compression force during normal walking of a person with a weight of 100 kg according to Stauffer et al., has been used in the simulation. The bone adaptation law, previously developed by our research team, has been used for the calculation of the remodelling processes. Results: A total bone mass loss of 2% in the tibia and 13% in the talus was calculated. The greater decline of density in the talus is due to its smaller size compared to the relatively large implant dimensions causing remodelling processes in the whole bone tissue. In the tibia, bone remodelling processes are only calculated in areas adjacent to the implant. Thus, a smaller bone mass loss than in the talus can be expected. There is a high agreement between the simulation results in the distal tibia and the literature regarding. Conclusions: In this study, strain-adaptive bone remodelling processes are simulated using the FE method. The results contribute to a better understanding of the biomechanical behaviour of the ankle joint and hence are useful for the optimisation of the implant geometry in the future.

Numerical investigations on the strain-adaptive bone remodelling in the periprosthetic femur: Influence of the boundary conditions

Background: There are several numerical investigations on bone remodelling after total hip arthroplasty (THA) on the basis of the finite element analysis (FEA). For such computations certain boundary conditions have to be defined. The authors chose a maximum of three static load situations, usually taken from the gait cycle because this is the most frequent dynamic activity of a patient after THA. Materials and methods: The numerical study presented here investigates whether it is useful to consider only one static load situation of the gait cycle in the FE calculation of the bone remodelling. For this purpose, 5 different loading cases were examined in order to determine their influence on the change in the physiological load distribution within the femur and on the resulting strain-adaptive bone remodelling. First, four different static loading cases at 25%, 45%, 65% and 85% of the gait cycle, respectively, and then the whole gait cycle in a loading regime were examined in order to regard all the different loadings of the cycle in the simulation. Results: The computed evolution of the apparent bone density (ABD) and the calculated mass losses in the periprosthetic femur show that the simulation results are highly dependent on the chosen boundary conditions. Conclusion: These numerical investigations prove that a static load situation is insufficient for representing the whole gait cycle. This causes severe deviations in the FE calculation of the bone remodelling. However, accompanying clinical examinations are necessary to calibrate the bone adaptation law and thus to validate the FE calculations.

Modification of the in vivo four-point loading model for studying mechanically induced bone adaptation

We modified the noninvasive, in vivo technique for strain application in the tibiae of rats (Turner et al,, Bone 12:73-79, 1991), The original model applies four-point bending to right tibiae via an open-loop, stepper-motor-driven spring linkage, Depending on the magnitude of applied load, the model produces new bone formation at periosteal (Ps) or endocortical surfaces (Ec.S). Due to the spring linkage, however, the range of frequencies at which loads can be applied is limited. The modified system replaces this design with an electromagnetic vibrator. A load transducer in series with the loading points allows calibration, the loaders' position to be adjusted, and cyclic loading completed under load central as a closed servo-loop. Two experiments were conducted to validate the modified system: (1) a strain gauge was applied to the lateral surface of the right tibia of 5 adult female rats and strains measured at applied loads from 10 to 60 N; and (2) the bone formation response was determined in 28 adult female Sprague-Dawley rats. Loading was applied as a haversine wave with a frequency of 2 Hz for 18 sec, every second day for 10 days. Peak bending loads mere applied at 33, 40, 52, and 64 N, and a sham-loading group tr as included at 64 N, Strains in the tibiae were linear between 10 and 60 N, and the average peak strain at the Ps.S at 60 N was 2664 +/- 250 microstrain, consistent with the results of Turner's group. Lamellar bone formation was stimulated at the Ec.S by applied bending, but not by sham loading. Bending strains above a loading threshold of 40 N increased Ec Lamellar hone formation rate, bone forming surface, and mineral apposition rate with a dose response similar to that reported by Turner et al, (J Bone Miner Res 9:87-97, 1994). We conclude that the modified loading system offers precision for applied loads of between 0 and 70 N, versatility in the selection of loading rates up to 20 Hz, and a reproducible bone formation response in the rat tibia, Adjustment of the loader also enables study of mechanical usage in murine tibia, an advantage with respect to the increasing variety of transgenic strains available in bone and mineral research. (Bone 23:307-310; 1998) (C) 1998 by Elsevier Science Inc. All rights reserved.

Increased tissue resistance in the nude mouse against Candida albicans without altering strain-dependent differences in susceptibility

Strain differences in tissue responses to infection with Candida albicans were examined in nude mice having susceptible (CBA/CaH) and resistant (BALB/c) parentage. Homozygous (nu/nu) mice of both strains were more resistant to systemic infection with C. albicans than heterozygous (nu/+) littermates as indicated by a reduction in both the severity of tissue damage and colony counts in the brain and kidney. However, the tissue lesions in nu/nu CBA/CaH mice were markedly more severe than those in nu/nu mice with the BALB/c background. This pattern was reflected in the greater fungal burden in the CBA/CaH strain. Analysis of cDNA from infected tissues using a competitive polymerase chain reaction excluded interferon-gamma (IFN-gamma), tumour necrosis factor-alpha (TNF-alpha), and interleukin 6 (IL-6) as mediators of the enhanced resistance of the nude mice. The results confirm that the different patterns of lesion severity in BALB/c and CBA/CaH mice do not involve T lymphocyte-mediated pathology, and are consistent with the hypothesis that strain-dependent tissue damage is not dependent on the effector function of macrophages or their precursors.

Down-regulation of expression of osteoblast and osteocyte markers in periodontal tissues associated with the spontaneous alveolar bone loss of interleukin-10 knockout mice

The aim of this study was to unravel the mechanisms by which interleukin (IL)-10, a potent pleiotropic cytokine, modulates alveolar bone homeostasis in C57BL/6 wild-type (WT) and IL-10 knockout (IL-10 KO) mice, evaluated at 8, 24, and 48 wk of age. Interleukin-10 KO mice presented significant alveolar bone loss when compared with WT mice, and this was not associated with changes in leukocyte counts or bacterial load. The levels of expression of messenger RNA (mRNA) for tumor necrosis factor-alpha (TNF-alpha), IL-1 beta, IL-6, transforming growth factor-beta (TGF-beta), receptor activator of nuclear factor kappa B ligand (RANKL), osteoprotegerin (OPG), and matrix metalloproteinase 13 (MMP13) were similar between both strains, whereas a significant decrease of tissue inhibitor of metalloproteinase 1 (TIMP1) mRNA expression was found at 48 wk in IL-10 KO mice. The osteoblast markers core binding factor alpha1 (CBFA1) and type I collagen (COL-I) were expressed at similar levels in both strains, whereas the levels of alkaline phosphatase (ALP) and osteocalcin (OCN), and those of the osteocyte markers phosphate-regulating gene endopeptidases (PHEX) and dentin matrix protein 1 (DMP1) were significantly lower in IL-10 KO mice. Our results demonstrate that the alveolar bone loss in the absence of IL-10 was associated with a reduced expression of osteoblast and osteocyte markers, an effect independent of microbial, inflammatory or bone-resorptive pathways.

The influence of cone-beam computed tomography and periapical radiographic evaluation on the assessment of periapical bone destruction in dog`s teeth

Objective. The aim of this study was to determine the influence of periapical radiographs, cone beam computed tomography (CBCT) sections, and cone beam volumetric data on the determination of periapical bone destruction in endodontically treated distal root canals of premolar canine teeth. Nontreated mesial roots were used as controls. Study design. Enterococcus faecalis strain (ATCC 29212) was inoculated into 30 root canals of 2 mongrel dogs to induce apical periodontitis. After 60 days, the root canals of the distal roots of the 11 mandibular and 4 maxillary premolars were endodontically treated (n = 15). The mesial root canals were used as controls (no treatment). The bone destruction was evaluated after 6 months by 5 evaluators using periapical radiographs and by CBCT (coronal and sagittal sections). After the experimental period, the area of the lesions in periapical radiographs and CBCT sections were measured in mm(2) using the ImageTool software. A single evaluator measured the volumetric data using the OsiriX software. The comparison between the diagnosis methods in treated root canals and controls was performed using parametric and nonparametric criteria. The Pearson correlation coefficient was computed between radiographic values and CBCT volumetric data in treated root canals and controls. Results. The results showed the presence of chronic apical periodontitis in every inoculated tooth. After 6 months, periapical radiographs, coronal CBCT sections, and volumetric data showed lower bone destruction in endodontically treated teeth in comparison with the control group (P < .05). The 5 evaluators found no differences between the apical periodontitis area of treated teeth and controls when CBCT sagittal sections were used (P > .05). No correlation was found between x-ray and CBCT volumetric values in treated root canals. Conclusions. Although selected CBCT sagittal sections showed similar values of bone destruction in endodontically and nontreated root canals, volumetric CBCT data showed that periapical lesions of endodontically treated root canals had half of the volume of periapical lesions in nontreated root canals. No relationship could be found between the periapical values of bone destruction and volumetric data found in CBCT of treated rood canals. (Oral Surg Oral Med Oral Pathol Oral Radiol Endod 2011; 112: 272-279)