Next Article in Journal
Computational Analysis of Self-Expanding and Balloon-Expandable Transcatheter Heart Valves
Previous Article in Journal
Single-Subject Analyses Reveal Altered Performance and Muscle Activation during Vertical Jumping
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulating Metaphyseal Fracture Healing in the Distal Radius

1
Scientific Computing Centre Ulm (UZWR), Ulm University, 89081 Ulm, Germany
2
Centre for Trauma Research Ulm, Institute of Orthopaedic Research and Biomechanics, 89081 Ulm, Germany
3
Institute for Biomechanics, ETH Zurich, 8093 Zurich, Switzerland
4
Institute for Information Systems, FHNW University of Applied Sciences and Arts Northwestern Switzerland, 4600 Olten, Switzerland
5
Department of Trauma Surgery, Medical University Innsbruck, 6020 Innsbruck, Austria
*
Authors to whom correspondence should be addressed.
Submission received: 28 December 2020 / Revised: 18 February 2021 / Accepted: 22 February 2021 / Published: 25 February 2021

Abstract

:
Simulating diaphyseal fracture healing via numerical models has been investigated for a long time. It is apparent from in vivo studies that metaphyseal fracture healing should follow similar biomechanical rules although the speed and healing pattern might differ. To investigate this hypothesis, a pre-existing, well-established diaphyseal fracture healing model was extended to study metaphyseal bone healing. Clinical data of distal radius fractures were compared to corresponding geometrically patient-specific fracture healing simulations. The numerical model, was able to predict a realistic fracture healing process in a wide variety of radius geometries. Endochondral and mainly intramembranous ossification was predicted in the fractured area without callus formation. The model, therefore, appears appropriate to study metaphyseal bone healing under differing mechanical conditions and metaphyseal fractures in different bones and fracture types. Nevertheless, the outlined model was conducted in a simplified rotational symmetric case. Further studies may extend the model to a three-dimensional representation to investigate complex fracture shapes. This will help to optimize clinical treatments of radial fractures, medical implant design and foster biomechanical research in metaphyseal fracture healing.

1. Introduction

Accounting for 17% of all emergency department visits [1], distal radius fractures (DRFs) are the most common long-bone fracture. Its incidence is increasing because of an active lifestyle and an ageing society [1,2,3,4].
Typically, a DRF is a high-impact trauma that occurs in young patients during a fall, while in older patients it is characterized by a low-energy fracture occurring above the distal articular surface of the radius.
Particularly the tapering of the cortical bone, where the bone is reinforced by the trabecular structure, was found to be a critical point of failure [2]. DRFs can be categorized according to different classification schemes, depending on the type of defect and trauma [2,5,6]. DRFs can be treated either conservatively or by surgical intervention where the main objectives are repositioning, fracture reduction and immobilization to maintain the reduction [7,8].
The healing process of trabecular bone fractures differs from that of most diaphyseal fractures, in which healing occurs through endochondral ossification within a callus enclosing the fracture. In diaphyseal fractures, soft tissue is replaced by cartilage, which subsequently ossifies. In trabecular bone fractures, bone forms directly through intramembranous ossification with almost no intermediary cartilage [9,10]. However, in vivo studies by Claes et al. [11,12] and Sandberg et al. [13,14] demonstrated that cancellous bone healing follows similar biomechanical rules as diaphyseal healing even though the observed healing patterns differ.
The development of numerical models able to simulate such complex biological and mechanical processes has progressed significantly over recent decades [15,16,17,18,19,20,21,22,23,24]. But to date, most models have generally focused solely on the modelling of diaphyseal fractures in long bones. Shefelbine et al. [25] and Simon et al. [26] initiated one of the only attempts at modelling trabecular healing, later presented as the Ulm tissue-level bone healing model [21]. Thus, they predicted the sequence of events in metaphyseal fracture healing in an academic example with the use of a model able to predict diaphyseal and metaphyseal fracture healing.
As a part of a multinational collaboration over the course of the last few years, one of the most extensive clinical patient dataset tracking of in vivo fracture healing of the distal radius was acquired. Patients were inspected at six time points within the first year after fracture [27]. This allowed us for the first time to compare in vivo fracture data with the predictions of our fracture healing simulation.
Therefore, the aim of this study is to develop and adapt a pre-existing tissue-level bone-healing model to the metaphyseal fracture healing of the distal radius. To determine whether the Claes et al. [12] hypothesis that metaphyseal bone healing displays similar mechanobiological effects as diaphyseal healing appears reasonable because the developed simulation tool with its origin in simulating diaphyseal fractures provides an ideal possibility. Once we have established a reliable model, it may be useful as a tool for different research questions or even to provide support for optimizing clinical procedures.

2. Materials and Methods

2.1. Tissue Differentiation

The previously validated Ulm tissue-level bone healing model [21,25,28,29,30,31,32,33] underlies the biomechanical observations of Pauwels [34] and a further mechano-regulating hypothesis by Claes and Heigele [35]. From mechanical quantifications, Simon et al. [21] derived local distortional and dilatational strains as the primary mechanical stimuli for the formation and removal of the following tissue types: woven and lamellar bone, cartilage and connective tissue. In addition to the mechanics, biological stimuli involved in the bone healing process are represented. Properties like the locally enclosing tissue types and level of vascularization augment the differentiation rules, which are evaluated by a fuzzy logic controller. Through the iterative numerical procedure, as depicted in Figure 1, the tissue differentiation processes of fracture healing in an arbitrary setting can be predicted. More detailed insights can be found in the Supplementary Materials to this paper and publications by Simon et al. [21] and Niemeyer et al. [28].

2.2. Metaphyseal Healing

As Claes et al. [12] demonstrated, both metaphyseal and diaphyseal bone healing appear to be regulated by the same biomechanical stimuli (dilatational and distortional strain). They found that metaphyseal areas, exposed to low strains led to intramembranous bone formation, whereas higher strains additionally provoked endochondral ossification or fibrocartilage formation. As this coincides with assumptions of the diaphysis, all current theories and tissue differentiation processes as implemented in the most recent healing model [28] were used and not adapted or further calibrated.

2.3. Homogenized Material for Trabecular Structure

The present model assumes a continuum, that is, organ-level, representation of the involved materials. Therefore, by modelling a metaphyseal bone healing process, the complex trabecular geometrical structure is not required, but its mechanical influence can be approximated by the bone volume content of trabecular bone tissue over time. Studies by Arias-Moreno et al. [36], Hosseini et al. [37], MacNeil and Boyd [38] and Varga et al. [39] evaluated the differences between a detailed, high-resolution finite-element model of the trabecular structure (µFE) in the metaphyseal radius and a representative homogenized material finite-element model (hFE). In those studies, strain analysis revealed that a homogenization of the trabecular structure is a sufficient approximation. Considering that representing all of the fine trabeculae as in a µFE analysis would increase the computational effort for a fracture-healing simulation to an insoluble problem. By this approach, the trabecular structure is represented by a spatially homogenized mixture of lamellar bone and soft tissue.

2.4. Rule of Mixture

Because the elastic modulus scales according to the cube of the mineral density [40], we developed a rule of mixtures based on the experimental findings of Carter and Hayes [40], as a result of which the averaged material properties of the tissue compositions in each point could be calculated. As a linear calculation of the Young’s modulus is not appropriate, as demonstrated by Simon et al. [21], a cubic rule of mixture should be used. Therefore, the tissue composition of N-sorted tissues with elasticity parameters E1 < E2 < … < EN and the corresponding concentrations c1...cN are combined following Pietsch et al. [41]:
E e l = E 1 c 1 + i = 2 N ( E i E i 1 ) ( j = i N c j ) 3

2.5. Modelling of a Compression Fracture

Regarding the modelling of a compression fracture in a metaphyseal bone, the fractured area must be defined accurately. Therefore, we introduced here a new tissue type that we term “crushed trabecular bone” (CTB). As the name suggests, it should represent the fractured trabeculae bone structures, which are damaged and pressed into each other. In vivo and in vitro observations [11,14,42,43,44,45,46,47] showed that this damaged structure still contributes to the primary stability of the bone but is resorbed over time.
Modelling this densification is achieved by increasing the bone concentration in the fractured area according to an exponential function (2). This approach follows in vivo observations and studies by Thurner et al. [48] and Hambli [49], who investigated the compression of trabecular structures. The calculated concentration was assigned to the new designed material type CTB:
c CTB ( t 0 ) = 1 e 3 C bone
Thus, the fractured area consists of a certain amount of CTB and connective tissue.
Material properties with a Young’s modulus of 6000 MPa and a Poisson’s ratio of 0.3 were used for pure CTB [50]. The densification approach, as well as the material properties, followed in vivo observations. Beside resorption, this new tissue type does not participate in any other tissue differentiation processes as depicted in Figure 1.
Niemeyer et al. [28] implemented a method to consider the influence of adjacent tissue on the local differentiation process. This distance-weighted average tissue potential (µ) was used to formulate an influence of the surrounding substances on the resorption speed of CTB ( r CTB ) . From this, each tissue type was represented by a specific resorption potential p [ 0 , 1 ] . Tissues with high perfusion tended to contribute more to the resorption of the fractured bone structure. Consequently, higher resorption potentials ( r ) were assigned according to Table 1.
Δ c CTB = r CTB μ ,
with μ = ( ( p p ( c w ( , t ) + c l ( , t ) + c c ( , t ) ) + p s c s ( , t ) + p c t b c ctb ( , t ) ) G adj ) .
Here, the Gaussian kernel G adj is defined as introduced by Niemeyer et al. [28]. Resorption rates followed recent findings on the resorption speed by Atkins et al. [51] and Christen et al. [52]. Through parameter variations (see Section Parameter variation) using a different patient’s geometry than the ones reported in Figure 4 and using comparisons with in vivo observations, it was concluded that a maximal resorption speed of 11% per day ( r c t b = 0.11 ) was used in all simulations.

2.6. Application to Distal Radius Fractures

Because we have access to datasets of distal radius fractured patients, the newly developed method was applied to four patient specific geometries. Patients were scanned at six time points with a Scanco medical Xtreme CT-II-scan system in weeks 1, 3, 5, 12, 24 and 52 after trauma. All patients were treated conservatively and the cast was removed before the third scan. During the healing period, patients were urged to perform a daily hand gripping exercise using a hand-grip fitness tool. Approval was obtained from the ethical committee of the Medical University Innsbruck for this study (AN2014-0374 344/4.31).
To simulate this DRF treatment, a patient-specific simulation geometry, initial tissue concentrations and boundary conditions were defined according to Figure 3B.
Regarding computational efforts, two-dimensional (2D) rotational symmetric assumptions were used. Comparing the simplified 2D rotational symmetric simulation domain in the following patient-specific geometry application, it was important to note that the entire radius is represented by the 2D geometry. Only the radial styloid process cannot be modelled by this 2D approach (Figure 2). Inaccuracies in the articular surface are covered by the principle of Saint Venant. The geometrical patient specification was achieved through parameterized geometry generation. Parameters such as cortical thickness, radial diameters and fracture site were measured in the scans after training and assistance by clinicians. Details on the investigated patients are attached in the Supplementary Materials.
The Geometries are discretised by a finite element mesh with quadratic hexahedral elements (3724–4143 elements), and the mesh was checked by mesh convergence analysis. Boundary conditions were chosen to approximate the real situation within the constraints of a 2D rotational symmetric simulation: The bone was fixed to the ground and compressed onto the articular surface (Figure 3A). This artificial load case was intended to represent the daily mechanical stimulus to which the radius is exposed.
During the healing period, the patients used a hand-grip fitness tool (approx. 100 N compression force). The resulting axial compression force from the contracting muscles in the wrist served as the mechanical stimulus for the fracture healing process. The occurring force onto the radius was calculated with the aid of a recently developed musculoskeletal inverse dynamics model of the human hand [53,54,55]. Results showed that, on average, the radial joint surface was exposed to an axial load of 351 N without great influence of the patients’ anatomical variability.
The initial model had homogenous material properties: trabecular bone, cortical bone and marrow. To obtain a realistic baseline tissue distribution within the trabecular bone, the model underwent tissue differentiation until a homeostatic situation was obtained (Figure 3C). In this case, the apparent stiffness of the tissue concentrations was sufficient, and the implemented differentiation processes achieved no change in tissue concentrations. Through this process, the homeostatic situation of the specifically determined geometry and boundary condition was obtained as exemplary depicted in Figure 3C.
Through X-ray scans, an approximation of the patients fractured area was identified and assigned a CTB concentration as defined by formula (2) and depicted in Figure 3D in blue.
To simulate the conservative treatment until week 5, a rotational symmetric brace with specific material parameters [56] (Young’s modulus of 950 MPa and Poisson’s ratio of 0.3) was modelled around the radius geometry (Position (a) in Figure 4). Muscles, skin and other soft tissue in between are assumed to be pure connective tissue. To compare the simulation with six in vivo longitudinal X-ray scans, the simulation results were visualized by an emulated pseudo-X-ray view. Concentration values of lamellar and woven bone and CTB were multiplied by factors for their corresponding bone mineral density [57].
The entire simulation framework was implemented in Python (Version 3.7, Python Software Foundation, Beaverton, OR, USA) using the commercial finite-element solver ANSYS (Release 2020 R1; ANSYS, Inc., Canonsburg, PA, USA). Running the entire simulation procedure for one geometry on a local, well equipped workstation results in calculation times of approximately 180 min.

3. Results

The predicted fracture healing results of four patient-specific distal radius fractures are compared in Figure 4 by an emulated pseudo-X-ray view to longitudinal in vivo X-ray scans of the corresponding patient.
Clearly, artefacts of the conservatively treated patients’ braces were visible in the first two screenings before the cast was removed prior to the third scan. In the fractured areas, the trabecular structure was compressed. During weeks 5 to 12, this compressed structure was resorbed and new bone formed. In most cases, the fractured area was totally healed after one year. In Patient A, the higher bone density in the area of the compression fracture clearly decreased from day 21 to day 81. Simultaneously, new bone was formed until day 168. The simulated results displayed comparable results. The compression fracture, modelled through an increased concentration of CTB, was resorbed in the first 3 to 5 weeks and was simultaneously replaced by new bone.
In Position (a), the modelled brace around the radius was depicted and was present until week 5. Patient A at day 36 (Position b) displayed newly formed woven bone in the fractured area, but this bone was not evenly distributed. As the CTB structure was resorbed and laterally a higher concentration of CTB was still present, not as much woven bone was formed in this spot, but even here the woven bone formation concentration attained 53% compared to a maximal 86.4% in the medial area. This observation was also obtained in the corresponding X-ray, where specific spots of the fractured area were also more mineralized than the surrounding. When comparing the subjects to each other, it was noted that larger fractured areas tended to heal more slowly. In addition, the predicted homeostatic tissue distribution showed that smaller radii tended to exhibit higher bone density than larger radii.
The time course of tissue concentrations in the fractured area revealed the underlying healing processes in more detail: Figure 5 demonstrates the mean concentration values in the fractured compartment of Patient D. First, the CTB was rapidly resorbed over time. After 59 days, the resorption speed decreased because from days 29 to 65 woven bone was formed. This mostly originated from intramembranous ossification, but because only little cartilage formation was observed at the beginning some endochondral ossification also occurred. Clearly visible was the maturation process, whereby the newly formed woven bone became lamellar bone, which was, after 178 days, the only tissue presents apart from connective tissue.

Parameter variation

Formula (3) introduced a new parameter into the differentiation processes: the maximal resorption speed of CTB ( r c t b ). Variations of this parameter of between 5 and 23% (representing tissue loss per day) are depicted in Figure 6. A higher resorption rate resulted in drastically decreased bone volume. Bone formation initially increased in all cases and subsequently decayed to the same concentration level.

4. Discussion

In this study, a novel model to simulate metaphyseal bone healing based on the tissue differentiation hypothesis of Claes and Heigele [35] was introduced. The underlying Ulm fracture healing model for diaphyseal fracture healing [21,28] was used and extended for its application to the metaphysis. In addition to geometrical and boundary adaptations, the model had to be substantially extended to consider further materials and processes to simulate the observed in vivo effects.
The newly introduced material, CTB, represents the compressed and fractured trabecular tissue. Trabeculae which are broken and of which the blood supply is interrupted, cannot contribute to the accumulation of mesenchymal stem cells and the resulting bone formation processes. CTB comprises a passive component. Because its material properties are similar to those of trabecular bone, it stabilizes the fractured area and allows rapid bone formation by its mere existence. Furthermore, the modelled CTB had no direct biological influence, as proposed by Aspenberg and Sandberg [13], in that no obvious relationship between newly formed bone and the old trabeculae was observed. Despite the fact that biopsies [13] showed that osteoconduction on old trabeculae is not an important process, we maintained the consideration that its surface offers ideal osteoconductive properties [58].
The modelled CTB was resorbed through a single rate in all biological processes. This rate was dependent on the surrounding tissue type, which might be an appropriate assumption, but modelling more specific resorption processes could be promising.
Modelling the geometry of the radius as a rotationally symmetric representation on the one hand might be a good approach for limiting computational costs, but on the other hand it limits the ability to represent a patient-specific radius shape. The radial styloid process as well as the asymmetrical shape of the lunate and scaphoid facet cannot be considered in such an assumption. On top, the fractured area and its positioning are limited to a horizontal expansion and can only be personalized through its height and vertical positioning. However, using 3D computations would have increased the computational time by up to 2 weeks for each healing simulation, but that is unacceptable in the development process of a new method.
Nevertheless, the model also shows that in 2D the assumptions on metaphyseal fracture healing appear reasonable, and further studies with 3D models can be elaborated using the same framework and modelling techniques.
From a large dataset, the four presented patients (Figure 4) were selected because they covered a broad variety of different aspects. The geometrical shape as well as their sex and age varied widely. The selected patients’ X-ray scans were compared to the simulation outcomes through an emulated X-ray presentation. Therefore, it must be noted that the emulated image was closer to a computed-tomography scanned image slice than to an actual X-ray view because the 3D geometric radioscopy was not modelled in this view. We nevertheless believe that a comparison of the real and simulated results is clearer through such a presentation. Because of the long time-lag between X-ray scans, a detailed comparison of the healing process was difficult. A qualitative comparison showed the model to be capable of simulating the metaphyseal section, the compression fracture and its healing progress. The simulations were able to predict the healing patterns observed in the X-ray scans from a decrease in bone density in the fractured area to new bone formation over the same time frame.
Discrepancies in the bone concentration of Patients A and D around the rotational axis (Positions c and d) might have arisen from the artificial load case as well as the boundary conditions required to maintain axial rotation symmetry. The significant difference in the diameter between Patients B and D was reflected in the bone concentration (Positions e and f) in the proximal region and might be explained by the fact that the model first calculated a homeostatic situation based on hand gripping force, which triggered bone formation in the relatively small cortical shell of Patient C (see (e) in Figure 4).
De Jong et al. [43] observed an increase of total bone (11%) and trabecular bone (20%) density until weeks 6–8. In the simulated results, a similar behaviour was observed. Increases of total bone concentration of 11 to 17% from days 10 to 90 were predicted in our simulations.
Published quantified results of in vivo investigations on distal radius fractures by de Jong et al. [42,43] were compared to congeneric quantities obtained by the simulations (Figure 7). Same behaviours were obtained for the axial stiffness, with a decrease from day 21 to 42. Density comparisons showed that the increase of healthy woven and later lamellar bone slightly lagged behind the in vivo observations but showed the same behaviour.
Histological studies by Aspenberg and Sandberg [13], Chen et al. [47], Jarry and Uhthoff [9] and Claes et al. [12] showed that metaphyseal bone formation occurs mainly by two patterns. Either lamellar bone forms on existing trabeculae, or woven bone emerges directly in the marrow [13] between the old trabeculae without direct connection to pre-existing bone [47]. Cartilage was seldom observed [12,13], but islands of chondrogenic tissue indicated limited endochondral ossification [12]. All observations showed that direct intramembranous bone formation is the dominant process of bone formation.
Those findings strongly correspond to our simulated results. We detected little cartilage formation in the fractured area (Figure 5) and hence little endochondral ossification. Direct woven-bone formation in the fractured area indicated that the mechanical stimuli were within the ranges originally proposed for the diaphysis. The course of tissue formation in Figure 4 thus reflected the histological observations that new woven bone formed in between the old trabeculae (CTB) and matured within a plausible time to lamellar bone. Additionally, the amount of bone formed corresponded to the surrounding bone concentration of the uninvolved parts of the bone. X-ray scans one year after trauma showed a healed radius, and it was difficult to determine the previously fractured area. The same characteristics occurred in the simulated results.
The histological comparison showed that the simulation model represents the metaphyseal fracture healing in an accurate manner. Because the model assumes tissue differentiation processes and biomechanical findings originally formulated for diaphyseal fracture healing, it can be stated that this simulation model corroborated the hypothesis of Claes et al. [12]. Despite the fact that some parameters of the model might differ slightly between the diaphysis and metaphysis, the same mechanobiological rules controlling bone and cartilage formation and resorption are able to predict fracture healing patterns both for the diaphysis and the metaphysis.

5. Conclusions

This metaphyseal bone-healing model provided defined and adjustable biomechanical conditions for a wide range of radius geometry representations. The simulations confirm the hypothesis of Claes et al. [12] for similar mechanobiological rules in fracture healing of both, the metaphysis and diaphysis. Furthermore, the results demonstrated reasonable predictions of intramembranous, and little endochondral, ossification in the fractured area. Therefore, the model appears appropriate for studying fracture healing under differing mechanical conditions and varying metaphyseal bones and fracture types. Investigating comorbidities such as osteoporosis and the implicated reduced bone density, as well as the extension of the model to a 3D representation where patient-specific, non-symmetric fractured areas can be investigated, are part of future studies. Those fracture healing simulations could then help to understand the healing pattern better and derive clinical treatments as well as biomechanical knowledge of metaphyseal fracture healing.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/biomechanics1010003/s1, Information S1: The Ulm Fracture Healing Model, Information S2: Details on the investigated patients.

Author Contributions

Conceptualization: F.N., U.S. and A.I.; Methodology: L.E., F.N. and U.S.; Software: L.E. and F.N.; Validation: L.E.; Investigation: L.E.; Resources: F.N. and U.S.; Data curation: K.S. and M.B.; Writing—original draft preparation: L.E., F.N.; Writing—review and editing: P.C., R.M., M.B., K.U., A.I. and U.S.; Visualization: L.E.; Supervision: K.U., A.I. and U.S.; Project administration: P.C., R.M., M.B., A.I. and U.S.; Funding acquisition: P.C., R.M., M.B., K.U., A.I. and U.S. All authors have read and agreed to the published version of the manuscript.

Funding

Funding by the SNF (320030L_170205), DFG (SI 2196/2-1, IG 18/19-1) and FWF (I 3258-B27) for the DACHFX Project is gratefully acknowledged. Support from the AO Foundation under AO Clinical Priority Programs Fracture Fixation in Osteoporotic Bone (FFOB) is acknowledged.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of Medical University Innsbruck for this study (AN2014-0374 344/4.31).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bonafede, M.; Espindle, D.; Bower, A.G. The direct and indirect costs of long bone fractures in a working age US population. J. Med. Econ. 2013, 16, 169–178. [Google Scholar] [CrossRef]
  2. MacIntyre, N.J.; Dewan, N. Epidemiology of distal radius fractures and factors predicting risk and prognosis. J. Hand Ther. 2016, 29, 136–145. [Google Scholar] [CrossRef] [PubMed]
  3. Porrino, J.A.; Maloney, E.; Scherer, K.; Mulcahy, H.; Ha, A.S.; Allan, C. Fracture of the distal radius: Epidemiology and premanagement radiographic characterization. AJR Am. J. Roentgenol. 2014, 203, 551–559. [Google Scholar] [CrossRef] [PubMed]
  4. De Putter, C.E.; van Beeck, E.F.; Looman, C.W.N.; Toet, H.; Hovius, S.E.R.; Selles, R.W. Trends in Wrist Fractures in Children and Adolescents, 1997–2009. J. Hand Surg. Am. 2011, 36, 1810–1815.e2. [Google Scholar] [CrossRef]
  5. Jayakumar, P.; Teunis, T.; Giménez, B.B.; Verstreken, F.; Di Mascio, L.; Jupiter, J.B. AO Distal Radius Fracture Classification: Global Perspective on Observer Agreement. J. Wrist Surg. 2017, 6, 46–53. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Müller, M.E.; Koch, P.; Nazarian, S.; Schatzker, J. The Comprehensive Classification of Fractures of Long Bones; Springer: Berlin/Heidelberg, Germany, 1990; ISBN 978-3-540-18165-1. [Google Scholar]
  7. Lippisch, R.; Lucas, B.; Schüttrumpf, J.P.; Piatek, S.; Walcher, F. Distale Radiusfraktur. Trauma Berufskrankh 2016, 18, 413–420. [Google Scholar] [CrossRef] [Green Version]
  8. Meena, S.; Sharma, P.; Sambharia, A.K.; Dawar, A. Fractures of distal radius: An overview. J. Family Med. Prim. Care 2014, 3, 325–332. [Google Scholar] [CrossRef]
  9. Jarry, L.; Uhthoff, H.K. Differences in healing of metaphyseal and diaphyseal fractures. Can. J. Surg. 1971, 14, 127–135. [Google Scholar]
  10. Uhthoff, H.K.; Rahn, B.A. Healing patterns of metaphyseal fractures. Clin. Orthop. Relat. Res. 1981, 295–303. [Google Scholar] [CrossRef]
  11. Claes, L.; Veeser, A.; Göckelmann, M.; Simon, U.; Ignatius, A. A novel model to study metaphyseal bone healing under defined biomechanical conditions. Arch. Orthop. Trauma Surg. 2009, 129, 923–928. [Google Scholar] [CrossRef]
  12. Claes, L.; Reusch, M.; Göckelmann, M.; Ohnmacht, M.; Wehner, T.; Amling, M.; Beil, F.T.; Ignatius, A. Metaphyseal fracture healing follows similar biomechanical rules as diaphyseal healing. J. Orthop. Res. 2011, 29, 425–432. [Google Scholar] [CrossRef]
  13. Aspenberg, P.; Sandberg, O. Distal radial fractures heal by direct woven bone formation. Acta Orthop. 2013, 84, 297–300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Sandberg, O.; Aspenberg, P. Inter-trabecular bone formation: A specific mechanism for healing of cancellous bone. Acta Orthop. 2016, 87, 459–465. [Google Scholar] [CrossRef] [Green Version]
  15. Lacroix, D.; Prendergast, P.J. A mechano-regulation model for tissue differentiation during fracture healing: Analysis of gap size and loading. J. Biomech. 2002, 35, 1163–1171. [Google Scholar] [CrossRef]
  16. Doblaré, M.; García, J.M.; Gómez, M.J. Modelling bone tissue fracture and healing: A review. Eng. Fract. Mech. 2004, 71, 1809–1840. [Google Scholar] [CrossRef]
  17. Gómez-Benito, M.J.; García-Aznar, J.M.; Kuiper, J.H.; Doblaré, M. A 3D computational simulation of fracture callus formation: Influence of the stiffness of the external fixator. J. Biomech. Eng. 2006, 128, 290–299. [Google Scholar] [CrossRef]
  18. Isaksson, H.; van Donkelaar, C.C.; Huiskes, R.; Ito, K. Corroboration of mechanoregulatory algorithms for tissue differentiation during fracture healing: Comparison with in vivo results. J. Orthop. Res. 2006, 24, 898–907. [Google Scholar] [CrossRef] [PubMed]
  19. Geris, L.; Gerisch, A.; Sloten, J.V.; Weiner, R.; van Oosterwyck, H. Angiogenesis in bone fracture healing: A bioregulatory model. J. Theor. Biol. 2008, 251, 137–158. [Google Scholar] [CrossRef]
  20. Checa, S.; Prendergast, P.J. A mechanobiological model for tissue differentiation that includes angiogenesis: A lattice-based modeling approach. Ann. Biomed. Eng 2009, 37, 129–145. [Google Scholar] [CrossRef]
  21. Simon, U.; Augat, P.; Utz, M.; Claes, L. A numerical model of the fracture healing process that describes tissue development and revascularisation. Comput. Methods Biomech. Biomed. Eng. 2011, 14, 79–93. [Google Scholar] [CrossRef]
  22. Ghiasi, M.S.; Chen, J.; Vaziri, A.; Rodriguez, E.K.; Nazarian, A. Bone fracture healing in mechanobiological modeling: A review of principles and methods. Bone Rep. 2017, 6, 87–100. [Google Scholar] [CrossRef]
  23. Wilson, C.J.; Schütz, M.A.; Epari, D.R. Computational simulation of bone fracture healing under inverse dynamisation. Biomech. Model. Mechanobiol. 2017, 16, 5–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Ren, T.; Dailey, H.L. Mechanoregulation modeling of bone healing in realistic fracture geometries. Biomech. Modeling Mechanobiol. 2020. [Google Scholar] [CrossRef] [PubMed]
  25. Shefelbine, S.J.; Augat, P.; Claes, L.; Simon, U. Trabecular bone fracture healing simulation with finite element analysis and fuzzy logic. J. Biomech. 2005, 38, 2440–2450. [Google Scholar] [CrossRef] [PubMed]
  26. Simon, U.; Shefelbine, S.J.; Augat, P.; Claes, L. Simulation of fracture healing in metaphyseal bone. J. Biomech. 2006, 39, S456. [Google Scholar] [CrossRef]
  27. Stock, K.; Marano, G.; Blauth, M.; Arora, R.; Müller, R.; Christen, P. Local In Vivo Assessment of Human Distal Radius Fractures by Time-Lapse HR-pQCTradius Fractures. In Proceedings of the 28th Joint Annual Meeting of The Swiss Association against Osteoporosis (SVGO) and the Swiss Bone and Mineral Society (SBMS), Bern, Switzerland, 3 May 2018. [Google Scholar]
  28. Niemeyer, F.; Claes, L.; Ignatius, A.; Meyers, N.; Simon, U. Simulating lateral distraction osteogenesis. PLoS ONE 2018, 13, e0194500. [Google Scholar] [CrossRef] [PubMed]
  29. Steiner, M.; Claes, L.; Ignatius, A.; Niemeyer, F.; Simon, U.; Wehner, T. Prediction of fracture healing under axial loading, shear loading and bending is possible using distortional and dilatational strains as determining mechanical stimuli. J. R. Soc. Interface 2013, 10, 20130389. [Google Scholar] [CrossRef] [Green Version]
  30. Steiner, M.; Claes, L.; Ignatius, A.; Simon, U.; Wehner, T. Disadvantages of interfragmentary shear on fracture healing--mechanical insights through numerical simulation. J. Orthop. Res. 2014, 32, 865–872. [Google Scholar] [CrossRef]
  31. Steiner, M.; Claes, L.; Ignatius, A.; Simon, U.; Wehner, T. Numerical simulation of callus healing for optimization of fracture fixation stiffness. PLoS ONE 2014, 9, e101370. [Google Scholar] [CrossRef] [Green Version]
  32. Wehner, T.; Claes, L.; Niemeyer, F.; Nolte, D.; Simon, U. Influence of the fixation stability on the healing time--a numerical study of a patient-specific fracture healing process. Clin. Biomech. 2010, 25, 606–612. [Google Scholar] [CrossRef]
  33. Wehner, T.; Steiner, M.; Ignatius, A.; Claes, L. Prediction of the time course of callus stiffness as a function of mechanical parameters in experimental rat fracture healing studies—A numerical study. PLoS ONE 2014, 9, e115695. [Google Scholar] [CrossRef]
  34. Pauwels, F. Eine neue Theorie ber den Einflu mechanischer Reize auf die Differenzierung der Sttzgewebe. Z. Anat. Entwickl. Gesch. 1960, 121, 478–515. [Google Scholar] [CrossRef]
  35. Claes, L.; Heigele, C.A. Magnitudes of local stress and strain along bony surfaces predict the course and type of fracture healing. J. Biomech. 1999, 32, 255–266. [Google Scholar] [CrossRef]
  36. Arias-Moreno, A.J.; Hosseini, H.S.; Bevers, M.; Ito, K.; Zysset, P.; van Rietbergen, B. Validation of distal radius failure load predictions by homogenized- and micro-finite element analyses based on second-generation high-resolution peripheral quantitative CT images. Osteoporos. Int. 2019, 30, 1433–1443. [Google Scholar] [CrossRef] [Green Version]
  37. Hosseini, H.S.; Dünki, A.; Fabech, J.; Stauber, M.; Vilayphiou, N.; Pahr, D.; Pretterklieber, M.; Wandel, J.; van Rietbergen, B.; Zysset, P.K. Fast estimation of Colles’ fracture load of the distal section of the radius by homogenized finite element analysis based on HR-pQCT. Bone 2017, 97, 65–75. [Google Scholar] [CrossRef]
  38. MacNeil, J.A.; Boyd, S.K. Bone strength at the distal radius can be estimated from high-resolution peripheral quantitative computed tomography and the finite element method. Bone 2008, 42, 1203–1213. [Google Scholar] [CrossRef] [PubMed]
  39. Varga, P.; Dall’Ara, E.; Pahr, D.H.; Pretterklieber, M.; Zysset, P.K. Validation of an HR-pQCT-based homogenized finite element approach using mechanical testing of ultra-distal radius sections. Biomech. Modeling Mechanobiol. 2011, 10, 431–444. [Google Scholar] [CrossRef]
  40. Carter, D.R.; Hayes, W.C. The compressive behavior of bone as a two-phase porous structure. J. Bone Joint Surg. Am. 1977, 59, 954–962. [Google Scholar] [CrossRef] [PubMed]
  41. Pietsch, M.; Niemeyer, F.; Simon, U.; Ignatius, A.; Urban, K. Modelling the fracture-healing process as a moving-interface problem using an interface-capturing approach. Comput. Methods Biomech. Biomed. Eng. 2018, 21, 512–520. [Google Scholar] [CrossRef] [Green Version]
  42. de Jong, J.J.A.; Heyer, F.L.; Arts, J.J.C.; Poeze, M.; Keszei, A.P.; Willems, P.C.; van Rietbergen, B.; Geusens, P.P.; van den Bergh, J.P.W. Fracture Repair in the Distal Radius in Postmenopausal Women: A Follow-Up 2 Years Postfracture Using HRpQCT. J. Bone Miner. Res. 2016, 31, 1114–1122. [Google Scholar] [CrossRef] [PubMed]
  43. de Jong, J.J.A.; Willems, P.C.; Arts, J.J.; Bours, S.G.P.; Brink, P.R.G.; van Geel, T.A.C.M.; Poeze, M.; Geusens, P.P.; van Rietbergen, B.; van den Bergh, J.P.W. Assessment of the healing process in distal radius fractures by high resolution peripheral quantitative computed tomography. Bone 2014, 64, 65–74. [Google Scholar] [CrossRef] [Green Version]
  44. Meyer, U.; de Jong, J.J.; Bours, S.G.P.; Keszei, A.P.; Arts, J.J.; Brink, P.R.G.; Menheere, P.; van Geel, T.A.C.M.; van Rietbergen, B.; van den Bergh, J.P.W.; et al. Early changes in bone density, microarchitecture, bone resorption, and inflammation predict the clinical outcome 12 weeks after conservatively treated distal radius fractures: An exploratory study. J. Bone Miner. Res. 2014, 29, 2065–2073. [Google Scholar] [CrossRef]
  45. Loi, F.; Córdova, L.A.; Pajarinen, J.; Lin, T.-H.; Yao, Z.; Goodman, S.B. Inflammation, fracture and bone repair. Bone 2016, 86, 119–130. [Google Scholar] [CrossRef] [Green Version]
  46. Schlickewei, C.W.; Kleinertz, H.; Thiesen, D.M.; Mader, K.; Priemel, M.; Frosch, K.-H.; Keller, J. Current and Future Concepts for the Treatment of Impaired Fracture Healing. Int. J. Mol. Sci. 2019, 20, 5805. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Chen, W.T.; Han, D.C.; Zhang, P.X.; Han, N.; Kou, Y.H.; Yin, X.F.; Jiang, B.G. A special healing pattern in stable metaphyseal fractures. Acta Orthop. 2015, 86, 238–242. [Google Scholar] [CrossRef] [PubMed]
  48. Thurner, P.J.; Erickson, B.; Jungmann, R.; Schriock, Z.; Weaver, J.C.; Fantner, G.E.; Schitter, G.; Morse, D.E.; Hansma, P.K. High-speed photography of compressed human trabecular bone correlates whitening to microscopic damage. Eng. Fract. Mech. 2007, 74, 1928–1941. [Google Scholar] [CrossRef]
  49. Hambli, R. Evaluation of apparent stress and strain at fracture of human trabecular bone under quasi-static compression load using micro-finite element simulations, Chapter of Evaluation of Apparent Stress and Strain at the Fracture of the Human Trabecular Bone under Quasi-Static Compression Load Using Micro-CT Finite Element Simulation. In Trabecular BONE: Morphology, biomechanics and clinical implications; Nova Science Publishers: Hauppauge, NY, USA, 2013; pp. 77–114. [Google Scholar]
  50. Wu, D.; Isaksson, P.; Ferguson, S.J.; Persson, C. Young’s modulus of trabecular bone at the tissue level: A review. Acta Biomater. 2018, 78, 1–12. [Google Scholar] [CrossRef] [PubMed]
  51. Atkins, P.R.; Stock, K.; Horling, L.; Ohs, N.; Collins, C.J.; Christen, P.; Blauth, M.; Müller, R. Bone Resorption and Formation Rates are Reduced in Late Stage Healing of Dominant Arm Distal Radius Fractures; Orthopaedic Research Society (ORS): Phoenix, AZ, USA, 2020. [Google Scholar]
  52. Christen, P.; Boutroy, S.; Ellouz, R.; Chapurlat, R.; van Rietbergen, B. Least-detectable and age-related local in vivo bone remodelling assessed by time-lapse HR-pQCT. PLoS ONE 2018, 13, e0191369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Engelhardt, L.; Vogelaar, D.; Weresch, P.; Urban, K.; Simon, U. A Musculoskeletal Inverse-Dynamics Model of the Human Hand Including Intrinsic Muscles; WCB: Dublin, Ireland, 2018. [Google Scholar]
  54. Engelhardt, L.; Melzner, M.; Havelkova, L.; Christen, P.; Dendorfer, S.; Simon, U. A New Musculoskeletal Anybody™ Detailed Hand Model; Annual meeting of European Society of Biomechanics: Vienna, Austria, 2019. [Google Scholar]
  55. Engelhardt, L.; Melzner, M.; Havelkova, L.; Christen, P.; Dendorfer, S.; Simon, U. A new musculoskeletal AnyBody™ detailed hand model. J. Comput. Methods Biomech. Biomed. Eng. 2020. [Google Scholar] [CrossRef] [PubMed]
  56. Schmidt, V.E.; Somerset, J.H.; Porter, R.E. Mechanical properties of orthopedic plaster bandages. J. Biomech. 1973, 6, 173–185. [Google Scholar] [CrossRef]
  57. Banzer, D.; Fabian, C.; Andresen, R.; Banzer, J.P.; Felsenberg, D.; Reisinger, W.; Neher, K.M. Knochendichte von Spongiosa und Kortikalis der Lendenwirbelsäule. Beziehungen zu Geschlecht, Alter und Wirbeldeformitäten in einem regionalen Kollektiv der europäischen Studie zur vertebralen Osteoporose (EVOS). Med. Klin. 1998, 93 Suppl. 2, 56–62, 64–65. [Google Scholar] [CrossRef] [PubMed]
  58. Burchardt, H. The biology of bone graft repair. Clinical Orthopaedics and Related Research 1983, 174, 28–42. [Google Scholar] [CrossRef]
Figure 1. The Ulm tissue-level bone healing model. The numerical procedure of solving an initial value problem to describe the biological tissue concentrations in an arbitrary setting uses a fuzzy logic controller. This process uses the input of a structurally finite element analysis, the enclosing tissue types and further biological settings to determine the concentration changes of woven and lamellar bone, cartilage and connective tissue. Further, the extensions of this study, the fractured trabecular (crushed trabecular) bone is depicted. The numerical procedure is adapted from Niemeyer et al. [18] CC-BY. The tissue differentiation diagram is adapted with permission of Georg Thieme Verlag KG (Stuttgart).
Figure 1. The Ulm tissue-level bone healing model. The numerical procedure of solving an initial value problem to describe the biological tissue concentrations in an arbitrary setting uses a fuzzy logic controller. This process uses the input of a structurally finite element analysis, the enclosing tissue types and further biological settings to determine the concentration changes of woven and lamellar bone, cartilage and connective tissue. Further, the extensions of this study, the fractured trabecular (crushed trabecular) bone is depicted. The numerical procedure is adapted from Niemeyer et al. [18] CC-BY. The tissue differentiation diagram is adapted with permission of Georg Thieme Verlag KG (Stuttgart).
Biomechanics 01 00003 g001
Figure 2. Two-dimensional rotational symmetric representation of a right distal radius with fracture area from anterior to posterior. Cortical shell in solid color and fractured area are defined where visible in the X-ray scan.
Figure 2. Two-dimensional rotational symmetric representation of a right distal radius with fracture area from anterior to posterior. Cortical shell in solid color and fractured area are defined where visible in the X-ray scan.
Biomechanics 01 00003 g002
Figure 3. Entire simulation procedure: Geometry definition from patients X-ray with highlighted fractured area (A). Boundary conditions of the idealized rotational symmetric geometry and finite element mesh (B). Fixation on the bottom and compression with 351 N from the top. Colours represent the standardized initial tissue distribution. Homeostatic situation after remodelling processes of 200 days is shown by an X-ray emulator (C). After adding the fractured area and brace (D) the proper fracture healing simulation can be run to obtain the results (E).
Figure 3. Entire simulation procedure: Geometry definition from patients X-ray with highlighted fractured area (A). Boundary conditions of the idealized rotational symmetric geometry and finite element mesh (B). Fixation on the bottom and compression with 351 N from the top. Colours represent the standardized initial tissue distribution. Homeostatic situation after remodelling processes of 200 days is shown by an X-ray emulator (C). After adding the fractured area and brace (D) the proper fracture healing simulation can be run to obtain the results (E).
Biomechanics 01 00003 g003
Figure 4. Longitudinal anterior–posterior X-ray scans of four patients 1–4 (above) and corresponding patient-specific simulation results (below) depicted by a pseudo-X-ray emulation. Patients were scanned at varying days after trauma. Patient data: (A) female, 51 years; fracture type, A2.1; Radius; T-score, −1.00; Z-Score, −0.30; bone mineral density (BMD), 0.52. (B) female, 86 years; fracture type, C1.1; Radius; T-Score, −1.90; Z-Score, −0.30; BMD, 0.48. (C) male, 33 years; fracture type, A2.1; Radius; T-Score, −0.90; Z-Score, −0.70; BMD, 0.53. (D) male, 45 years; fracture type, A2.1; Radius; T-Score, −1.90; Z-Score, −1.60; BMD: 0.59.
Figure 4. Longitudinal anterior–posterior X-ray scans of four patients 1–4 (above) and corresponding patient-specific simulation results (below) depicted by a pseudo-X-ray emulation. Patients were scanned at varying days after trauma. Patient data: (A) female, 51 years; fracture type, A2.1; Radius; T-score, −1.00; Z-Score, −0.30; bone mineral density (BMD), 0.52. (B) female, 86 years; fracture type, C1.1; Radius; T-Score, −1.90; Z-Score, −0.30; BMD, 0.48. (C) male, 33 years; fracture type, A2.1; Radius; T-Score, −0.90; Z-Score, −0.70; BMD, 0.53. (D) male, 45 years; fracture type, A2.1; Radius; T-Score, −1.90; Z-Score, −1.60; BMD: 0.59.
Biomechanics 01 00003 g004
Figure 5. Mean tissue concentrations in the fractured area. “Crushed trabecular bone” (CTB) is resorbed over time and woven, and later lamellar bone is formed. Total bone describes the sum of woven and lamellar bone. A small amount of cartilage and subsequently endochondral ossification is observed. However, intramembranous ossification dominates the overall process.
Figure 5. Mean tissue concentrations in the fractured area. “Crushed trabecular bone” (CTB) is resorbed over time and woven, and later lamellar bone is formed. Total bone describes the sum of woven and lamellar bone. A small amount of cartilage and subsequently endochondral ossification is observed. However, intramembranous ossification dominates the overall process.
Biomechanics 01 00003 g005
Figure 6. Parameter variations of the “crushed trabecular bone” (CTB) resorption speed ( r c t b ) as a function of time. Bone concentration consists of woven, lamellar bone and CTB representing an overall bone concentration in the fractured area, depicted in blue.
Figure 6. Parameter variations of the “crushed trabecular bone” (CTB) resorption speed ( r c t b ) as a function of time. Bone concentration consists of woven, lamellar bone and CTB representing an overall bone concentration in the fractured area, depicted in blue.
Biomechanics 01 00003 g006
Figure 7. Left: Comparison of the bone density in distal radius fractures by de Jong et al. (2016) and the mean simulated values of pseudo-X-ray density of bone for the investigated patients in function of days after trauma. Right: Comparison of axial stiffness values, calculated by de Jong et al. (2016) and mean Young’s modulus values in the fractured area of the simulated patients in function of days after trauma.
Figure 7. Left: Comparison of the bone density in distal radius fractures by de Jong et al. (2016) and the mean simulated values of pseudo-X-ray density of bone for the investigated patients in function of days after trauma. Right: Comparison of axial stiffness values, calculated by de Jong et al. (2016) and mean Young’s modulus values in the fractured area of the simulated patients in function of days after trauma.
Biomechanics 01 00003 g007
Table 1. CTB resorption potentials and involved tissue type abbreviations.
Table 1. CTB resorption potentials and involved tissue type abbreviations.
Resorption Potential p
Perfused tissue (pp)0.8
Soft tissue (ps)1.0
CTB (pctb) 0.1
Tissue ConcentrationsVariable
Woven bone c w
Lamellar bone c l
Crushed trabecular bone c c t b
Cartilage c c
Connective/soft tissue c s
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Engelhardt, L.; Niemeyer, F.; Christen, P.; Müller, R.; Stock, K.; Blauth, M.; Urban, K.; Ignatius, A.; Simon, U. Simulating Metaphyseal Fracture Healing in the Distal Radius. Biomechanics 2021, 1, 29-42. https://0-doi-org.brum.beds.ac.uk/10.3390/biomechanics1010003

AMA Style

Engelhardt L, Niemeyer F, Christen P, Müller R, Stock K, Blauth M, Urban K, Ignatius A, Simon U. Simulating Metaphyseal Fracture Healing in the Distal Radius. Biomechanics. 2021; 1(1):29-42. https://0-doi-org.brum.beds.ac.uk/10.3390/biomechanics1010003

Chicago/Turabian Style

Engelhardt, Lucas, Frank Niemeyer, Patrik Christen, Ralph Müller, Kerstin Stock, Michael Blauth, Karsten Urban, Anita Ignatius, and Ulrich Simon. 2021. "Simulating Metaphyseal Fracture Healing in the Distal Radius" Biomechanics 1, no. 1: 29-42. https://0-doi-org.brum.beds.ac.uk/10.3390/biomechanics1010003

Article Metrics

Back to TopTop