|
|||||||||||
|
Introduction |
|
|
Bone failure projects |
|
![]() |
Bone remodeling projects |
![]() |
Methods |
| Group Members | |
| Recent Publications |
Our group’s research is dedicated to simulate aspects of bone in its multiple scales. The temporal scales range from milliseconds (10-3s) for fracture whole bone fracture to several years (107s) for bone architecture changes, as in the case of osteoporosis. The spatial scale on the other hand ranges from a few micrometers (10-5m), which isthe size of a single trabecula to several centimeters (10-2m), for an entire vertebra. Our research covers and combines biomechanics in the range of all these scales. We investigate multiscale biological models with computer simulations and develop experimental methods for their validation. ↑Top↑
Bone failure projects lay in the milliseconds temporal scale. Customarily, we further differentiate these projects by their spatial scale.
Today, a patient's risk for osteoporotic fractures is estimated using DXA. In research, high-resolution quantitative computed tomography, such as HR-pQCT, is used to investigate bone quality. The aim of this project is to establish alternative assessments based on large-scale FE simulations and HR-pQCT. Ideally, this project will allow better understanding of the underlying changes in the osteopenic bone.
In the course of this project, a parallel linear FE solver ParFE is extended to account for large deformations and nonlinear material properties. It remains to be evaluated if and to what extent nonlinear simulations will allow a more accurate prediction of patient specific fracture risk.
On a micro-architectural level we are investigating networks of trabecular bone. Linear elastic FE solvers are the tool we use to predict apparent stiffness of a trabecular network. Some of the apparent limitations are the the ability to define fracture relevant parameters, such as the ultimate strength or the local failure of a trabecular network. Within this project, we explore different methods to enhance the predictability of failures in trabecular networks.
Nonlinear post-elastic mechanical properties of trabecular bone tissue play a key role in understanding fracture onset. This project aims to overcome the small-sample limitations in mechanical testing by using a novel micro-bending setup adapted to a single trabecula. The combination of fluorescence microscopy and FE simulation permits the measurement of characteristic load-displacement curves and the indirect modeling of material tissue properties in the plastic region. After a first validation study, this approach will allow a better insight into the material fracture mechanism, as well as into the role of collagen and calcium density in fracture.
↑Top↑
A larger time scale is used to simulate architectural changes of bone during a lifetime. Bone remodeling projects simulate different biological and biomechanical effects on the microstructure of whole bones. Depending on the species, the size of the investigated whole bone samples varies from a few millimeters to several centimeters.
In the recent years, bone remodeling algorithms based on theoretical bone adaptation models have been developed to simulate the process in bone and its microstructure. However, restricted computational possibilities have limited such simulations in size, resolution and timescale. Taking advantage of the high computational power of recently developed parallel computer tools we hope to overcome these limitations. With the simulation of bone adaptation of whole human vertebrae at high resolutions we cross a border and model the long term microstructural changes of the whole bone due to osteoporosis.
Furthermore, bone adaptation simulations are expanded to simulate long term effects of clinical interventions. Vertebroplasty is a minimally invasive technique to stabilize fractured vertebra during which the vertebral body is being filled with bone cement. Efforts are made to simulate the bone behavior under the changed mechanical conditions. Finally, the anabolic or catabolic effect of signal agents is modeled in the bone adaptation simulations.
Although a promising method, in silico modeling of the bone adaptation currently lacks experimental validation and is mostly dismissed as a supplement for the investigations of the existing biological theories. Further development of the simulation algorithm, and most importantly its validation against experimental data, can help establish in silico modeling as a self-standing technique, thus reducing the number of animal studies, as well as opening a new window of possible applications in the clinics. With this project we attempt to validate in silico simulations of the effects of aging, disease (osteoporosis) and associated treatments with in vivo experimental studies conducted in mice. We thus hope to establish this technique as a reliable computational model of bone adaptation.
Our group focuses on the development of computational and experimental methods. We have established a physics-based in silicoframework for the simulation of bone remodeling, as well as a nonlinear FE code for predicting failure onset. For the model validation we have constructed a mechanical setup, which allows the examination of tissue properties and study of bone failure behavior. Finally, standard post-processing of the experimental data has been further improved with an image-based algorithm, which allows determination of the 3D strains based on micro-CT images.
The aim of this work is to develop a nonlinear finite element (FE) solver, which is based on the existing parallel FE package ParFE, introduced earlier at ETH Zurich (http://sourceforge.net/projects/parfe). ParFE is using the Trilinos framework (http://trilinos.sandia.gov/) and the CG solver AztecOO in combination with a smoothed aggregation based algebraic multigrid preconditioner. Specifically, changes to the simulation code will involve the implementation of a Newton-Raphson (NR) nonlinear solving scheme using the routines of ParFE for the individual linear solves. This resembles the approach of modern commercial FE solvers used for example in structural mechanics. However, a special-purpose FE solver had to be developed in order to handle the large number of elements (order of billions) and for efficient use of very large number of processors. Previous findings suggest that nonlinear geometrical effects in combination with complex tissue properties play an important role in the mechanical failure of the bone microstructure. Although the best material model for bone still has to be determined, we strongly believe that micro-FE simulations could be significantly improved in their accuracy by incorporating nonlinear effects.
Using mechanical testing in combination with high-resolution 3D imaging allows the observation of potential failure mechanisms of the bone microstructure. Additionally, strain mapping with deformable image registration allows determination of local displacements and strains, thus enabling a quantitative assessment of the relevant failure modes.
The current validation study suggests that the precision of the deformable image registration is roughly four-fold better than the voxel size (7 µm and 30 µm, respectively). A validation of the computed strains is underway.
Bone microstructure is often studied with the help of high-resolution three-dimensional computer tomography scans. The biophysical models are implemented into flexible simulation tools that allow easy adaptation and exchange of the underlying theoretical principles. The simulation tools are based on advection equations and micro FE analysis (ParFE). In order to handle large input data sets and extensive computation steps the simulations are performed at the Swiss National Supercomputing Center,which allows parallel computation.
We have selected a bending test method to measure the mechanical properties of single trabecula. To facilitate such experiments a prototype micro bending machine was designed. The micromechanical measurements are carried out under direct visual observation using a light microscope.
Networks of trabecular bone are investigated with a combination of a micro-compression device and synchrotron radiation computer tomography (Swiss Light Source, SLS). This procedure, known as Image Guided Failure Assessment (IGFA), is planned for the investigation of 3D failure mechanisms in human vertebral bone specimens. Apart from obvious mechanical and biological relevance, we hope that this data will allow validation of our mechanical simulation models.
↑Top↑
Müller, Ralph
Lorenzetti, Silvio
Badilatti, Sandro
Carretta, Roberto
Christen, David
Levchuk, Alina
Li, Zihui
Ruffoni, Davide
Rickenbacher Dominik
Schulte Friederike
Zwahlen, Alexander
Back to Search | Style:
Simple ·
Quote | Export:
EndNote XML ·
BibTex
View more details
Link to Website / File
20 records found
[1] 2 Next 10 items »
| Title | Author(s) | Year | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Strain-adaptive in silico modeling of bone adaptation - A computer simulation validated by in vivo micro-computed tomography data ![]() |
Schulte, F. A.; Zwahlen, A.; Lambers, F. M.; Kuhn, G.; Ruffoni, D.; Betts, D.; Webster, D. J.; Müller, R. | 2013 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Export: EndNote XML · BibTex |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
In vivo validation of predictive models for bone remodeling and mechanobiology ![]() |
Levchuk, A.; Müller, R. | 2013 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Export: EndNote XML · BibTex |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Strain-adaptive in silico modeling of bone adaptation - A computer simulation validated by in vivo micro-computed tomography data ![]() |
Schulte, F. A.; Zwahlen, A.; Lambers, F. M.; Kuhn, G.; Ruffoni, D.; Betts, D.; Webster, D. J.; Müller, R. | 2012 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Export: EndNote XML · BibTex |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Imaging of cellular spread on a three-dimensional scaffold by means of a novel cell-labeling technique for high-resolution computed tomography ![]() |
Thimm, B. W.; Hofmann, S.; Schneider, P.; Carretta, R.; Müller, R. | 2012 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Export: EndNote XML · BibTex |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Increased marrow adiposity in premenopausal women with idiopathic osteoporosis ![]() |
Cohen, A.; Dempster, D. W.; Stein, E. M.; Nickolas, T. L.; Zhou, H.; McMahon, D. J.; Müller, R.; Kohler, T.; Zwahlen, A.; Lappe, J. M.; Young, P.; Recker, R. R.; Shane, E. | 2012 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Export: EndNote XML · BibTex |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Deformable image registration and 3D strain mapping for the quantitative assessment of cortical bone microdamage ![]() |
Christen, D.; Levchuk, A.; Schori, S.; Schneider, P.; Boyd, S. K.; Müller, R. | 2012 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Export: EndNote XML · BibTex |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Towards patient‐specific material modeling of trabecular bone post‐yield behavior ![]() |
Carretta, R.; Lorenzetti, S.; Müller, R. | 2012 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Export: EndNote XML · BibTex |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Abnormal bone microarchitecture and evidence of osteoblast dysfunction in premenopausal women with idiopathic osteoporosis ![]() |
Cohen, A.; Dempster, D. W.; Recker, R. R.; Stein, E. M.; Lappe, J. M.; Zhou, H.; Wirth, A. J.; van Lenthe, G. H.; Kohler, T.; Zwahlen, A.; Müller, R.; Rosen, C. J.; Cremers, S.; Nickolas, T. L.; McMahon, D. J.; Rogers, H.; Staron, R. B.; Lemaster, J.; Shane, E. | 2011 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Export: EndNote XML · BibTex |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Computational finite element bone mechanics accurately predicts mechanical competence in the human radius of an elderly population ![]() |
Mueller, T. L.; Christen, D.; Sandercott, S.; Boyd, S. K.; van Rietbergen, B.; Eckstein, F.; Lochmuller, E. M.; Müller, R.; van Lenthe, G. H. | 2011 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Export: EndNote XML · BibTex |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Towards validation of computational analyses of peri-implant displacements by means of experimentally obtained displacement maps ![]() |
Basler, S. E.; Mueller, T. L.; Christen, D.; Wirth, A. J.; Müller, R.; van Lenthe, G. H. | 2011 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Export: EndNote XML · BibTex |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
[1] 2 Next 10 items »
Back to Search | Style:
Simple ·
Quote | Export:
EndNote XML ·
BibTex
Wichtiger Hinweis:
Diese Website wird in älteren Versionen von Netscape ohne
graphische Elemente dargestellt. Die Funktionalität der
Website ist aber trotzdem gewährleistet. Wenn Sie diese
Website regelmässig benutzen, empfehlen wir Ihnen, auf
Ihrem Computer einen aktuellen Browser zu installieren. Weitere
Informationen finden Sie auf
folgender
Seite.
Important Note:
The content in this site is accessible to any browser or
Internet device, however, some graphics will display correctly
only in the newer versions of Netscape. To get the most out of
our site we suggest you upgrade to a newer browser.
More
information