- 1Applied Mechanics, Saarland University, Saarbrücken, Germany
- 2Department of Diagnostic and Interventional Radiology, University Hospital Aachen, Aachen, Germany
- 3Faculty of Medicine, BG Hospital Tuebingen, University Hospital Tuebingen on Behalf of the Eberhard-Karls-University Tuebingen, Tuebingen, Germany
Purpose: The prescribed amount of weight-bearing after tibial plateau fractures is controversial because it affects osteosynthetic construct stability and fracture healing. We aim to introduce a simulation model that adequately predicts the effects of different weight-bearing amounts on stability and healing, based on the patient’s individual fracture pattern and treatment construct.
Methods: To safely test different amounts of weight-bearing limits, we first extracted knee joint forces for different weight-bearing limits from musculoskeletal simulation based on monitoring data of 22 uninjured participants. Correct loading was ensured with a force-measuring insole. We then tested three patients after tibial plateau fracture with their current weight-bearing level and constructed a simulation model determining implant stress, knee joint force, and fracture gap interfragmentary strain. The patient-specific weight-bearing level was then substituted for weight-normalized uninjured participant data to test different weight-bearing levels in the simulation model.
Results: The simulation model calculated individual construct stiffness and interfragmentary strain at different weight-bearing levels following the clinical course. When comparing the patient’s individual weight-bearing input with the weight-normalized input of the uninjured participants at the same level, comparable knee joint forces were extracted, showing the feasibility of this approach.
Conclusion: Using an adapted reference movement database, the model allows the determination of safe weight-bearing ranges concerning construct stability and fracture healing based on individual fracture morphology and treatment without exposing patients to excessive weight-bearing. Future studies can test this approach in more extensive patient-number studies and different treatment situations.
1 Introduction
Tibial plateau fractures represent about 1% of all fractures in adults and are considered among the most complex joint injuries to be treated (Court-Brown and Caesar, 2006; Delp et al., 2007). Resulting from high-energy mechanisms in younger patients, as well as low-energy trauma with associated poor bone quality in the elderly, this fracture entity is associated with a limited prognosis due to cartilage, as well as additional soft tissue injury (Arnold et al., 2017; Gaston et al., 2005; Neal et al., 2017). While the specific osteosynthetic construct, depending on the fracture configuration, is still debated, the mainstay of treatment is open reduction and internal fixation (Canton et al., 2022; Mcnamara et al., 2015). According to the AO recommendations for weight-bearing after tibial plateau fracture treatment, patients should be kept non-weight-bearing for 10–12 weeks. However, a recent survey analysis among Dutch trauma surgeons has shown that this recommendation is only inconsistently followed by surgeons and that the prescribed start of at least some weight-bearing mostly ranges between 0 and 6 weeks (van der Vusse et al., 2017). Furthermore, the general compliance of patients to non- or limited weight-bearing regimens is limited in several studies (Braun et al., 2017; Chiodo et al., 2016; Kammerlan et al., 2018). In addition, current studies have suggested that some degree of early or even permissive full weight bearing may be safe in select cases (Canton et al., 2022; The et al., 2017; Williamson et al., 2018).
Multiple studies have focused on determining a safe range for weight-bearing after these injuries, showing that at least a moderate amount of weight-bearing early on is not secure with standard open reduction internal fixation constructs, but might be beneficial for the patient regarding reported outcomes (Arnold et al., 2017; Canton et al., 2022; Iliopoulos and Galanis, 2020; Pishtiwan Hassan Shaker Kalmet et al., 2019; Williamson et al., 2018). Studies have gone so far as to recommend permissive weight-bearing, according to the patient’s comfort, but even then, patients reached full weight-bearing on average only after 14 weeks (Pishtiwan Hassan Shaker Kalmet et al., 2019). Recent randomized and comparative studies demonstrated that immediate or earlier weight bearing after internal fixation can lead to superior functional outcomes without increased radiological complications, and that earlier loading does not increase complication rates or compromise fracture union, further questioning the necessity of prolonged restriction (Erick et al., 2023; Ibrahim et al., 2025). Depending on the fracture type and treatment, there is limited evidence that higher average joint reaction forces at the early point can lead to an increase in intra-articular step-off formation, thus showing the clinical need for a reliable tool to determine safe ranges of weight-bearing on an individual level (Thewlis et al., 2015).
To accurately determine the safe range of weight-bearing after tibial plateau fractures, musculoskeletal simulation can help assess individual construct stability. This is a reliable method to assess implant stresses, joint reaction forces, and interfragmentary fracture strain in other lower extremity fracture entities (Braun et al., 2021; Marmor et al., 2020). While standard models of the proximal tibia have already been constructed to test different mechanical conditions in fracture models (Carrera et al., 2018; Carrera et al., 2016), to our knowledge, no model has been developed to calculate individual, case-specific construct stability and the resulting interfragmentary strain. Previous research often relies on generic or population-averaged models that do not incorporate patient-specific anatomical data or individualized loading conditions (Cunningham et al., 2017). These standard models typically apply simplified loading assumptions, uniform implant configurations, or generalized boundary conditions that may not accurately reflect the biomechanical environment of a healing fracture in a specific patient.
The aim of the current study was thus to construct a simulation model that adequately predicts the effect of different weight-bearing amounts on the stability and interfragmentary strain of fracture treatment, based on the patient’s individual fracture pattern and treatment construct.
2 Materials and methods
This study recruited 22 healthy subjects, including 11 women and 11 men with no history of lower extremity fracture, musculoskeletal disease, or neurological impairment. These individuals served exclusively as a normative reference group for musculoskeletal simulation input data. In contrast, the clinical part of the study involved three postoperative patients with tibial plateau fractures treated by open reduction and internal fixation, who were included to test the feasibility of the patient-specific simulation workflow. We used the Xsens™ system, known for its precision in recording human movement, to capture motion data. The participants actively participated in the Timed Up and Go (TUG) test, a standardized assessment that measures their mobility and functional capacity (Podsiadlo and Richardson, 1991). During the TUG test, individuals had to rise from a seated position, walk a short distance, turn around, and return to the seated position, with the primary outcome being their completion time, measured in seconds. Before conducting the test, we provided all participants with a comprehensive explanation of the TUG procedure and its significance. We alternated the partial load conditions between subjects, with one group (n = 11) of 20 kg and the other group (n = 11) an equivalent load of half their body weight as partial weight bearing. To ensure participants adhered to their assigned partial weight, we actively employed a Moticon™ insole system to monitor and record the load during the TUG test. We actively instructed participants in properly handling crutches using the 3-point gait method to ensure consistency and safety.
Furthermore, participants actively assessed their perception of the partial load using a household scale. Before the formal testing, all participants walked over the scale ten times to familiarize themselves with the sensation of the ordered partial load. We conducted the TUG test in two rounds: first, with a full load to establish reference values, and then with the arranged partial load, determined by the subject’s group assignment (20 kg or half body weight). This dual-round approach allowed for a comprehensive evaluation of each participant’s performance under full and partial load conditions, enabling an active comparison of their mobility and functional capacity in different load scenarios. The TUG test, a quick, reliable, and standardized assessment, actively evaluated participants’ mobility and functional capacity throughout this study (Izzah et al., 2019; Kristensen et al., 2007).
The workflow simultaneously measures the subject using the Xsens™ motion capturing system (Supplementary Figure S1), Moticon™ measurement insoles (Supplementary Figure S2), and Novel™ loadpads (Supplementary Figure S3). The Xsens™ motion data is exported as. bvh (BioVision Hierarchy) files and imported into AnyBody™ software (Supplementary Figure S4), a musculoskeletal simulation tool that works seamlessly with Xsens™. Data from the measurement insoles and loadpads are imported as. txt files into AnyBody™, where the musculoskeletal model is scaled to the subject’s anthropometric data. This data is then used to calculate the joint and muscle forces inversely. The interface between Xsens™ and AnyBody™ ensures seamless data integration. Xsens™ gait data is exported as. bvh files for each subject and processed in AnyBody™ using the internal AMMR (AnyBody™ Model Repository), an open collection of examples and models of the musculoskeletal system, to map the gait cycle accurately. After completing the motion data acquisition with the 22 healthy subjects using the Xsens™ system, the data were analyzed using AnyBody™. Based on the recorded data, AnyBody™ allows us to perform a detailed analysis of joints, muscles, and moments. This analysis is critical for evaluating the effects of partial loads, particularly on knee forces, as part of our study of tibial plateau fractures (Figure 1). The knee joint forces included for a controlled partial load (20 kg, or 50% of body weight) were derived from healthy individuals. Patient data were used to determine the force values for full weight-bearing. To allow for inter-subject comparison, all computed knee joint reaction forces from the AnyBody™ musculoskeletal model were normalized to each participant’s body weight (%BW). This was done by dividing the joint force values (in Newtons) by the individual’s body weight (in Newtons) and multiplying by 100. The resulting values are expressed as a percentage of body weight (%BW) and plotted across the gait cycle. The systems used for inertial motion capture (Xsens™), insole-based force measurement (Moticon™), and musculoskeletal modeling (AnyBody™) provide a validated and reproducible workflow, as demonstrated in previous studies (Burns et al., 2019; Andersen et al., 2009; Muro-de-la-Herran et al., 2014). Detailed instrumentation specifications are in the Supplementary Material.

Figure 1. Integrating the Xsens™ motion capturing system with AnyBody™ uses input data from captured movements and the subsequent realistic boundary conditions to simulate knee joint reactions.
Furthermore, our research goes beyond healthy subjects as we actively measure actual patients with their prescribed partial weight-bearing. We follow a similar measurement workflow to healthy subjects and perform these measurements in the clinic using the Xsens™ motion measurement system. As with healthy subjects, the data obtained from the patients is also analyzed using AnyBody™ software. Figure 2 illustrates an example of knee joint force curves for a patient, using his motion-capturing data. Additionally, we utilize the patients’ CT images to generate patient-specific 3D models through segmentation and downstream image processing, which are then used in finite element (FE) simulations. This approach enables us to create a database of partial loads, providing insight into how various loads impact the patient’s musculoskeletal dynamics, with a particular focus on tibial plateau fractures and their treatment.
The first step in the FE analysis methodology is to create geometric models using the individual DICOM image stack of the patient’s postoperative CT scan and a six-rod bone density calibration phantom (QRM-BDC/6, QRM GmbH, Moehrendorf, Germany) (Braun et al., 2021). Therefore, the images are segmented into masks using an adaptive thresholding procedure based on the calibration phantom. The calibration phantom with known density values establishes a relationship between Hounsfield Units and bone mineral density. After completing the segmentation, the treating trauma surgeon visually controlled the segmentation results of the fracture gap and made corrections if necessary. Then, we create high-resolution adaptive finite element meshes using the FE Mesher from Simpleware™ (Simpleware™, Synopsys, Mountain View, United States), employing C3D10 elements and defining node sets at the distal and proximal ends of the bone. Specify the homogeneous materials with standard properties within Simpleware™, sourcing values for the implant and fracture gap from the literature (Claes and Heigele, 1999; Imam and Fraker, 1996). To establish the relationship between elasticity and bone density, map the grayscale values of the CT data to the Hounsfield scale and mechanical local bone properties (Cattaneo et al., 2001; Hvid et al., 1989; Rho et al., 1995; Zannoni et al., 1999). The exact workflow is outlined in the Supplementary Material, including additional information about the mesh (Supplementary Figure S5). Following previous studies, assume an isotropic heterogeneous material with a varying Young’s modulus and a fixed Poisson ratio (Nir et al., 2009; Yosibash et al., 2007). Define cortical and trabecular bone mapping based on local ash density and equivalent mineral density, as described in prior research (Braun et al., 2021; Brent Edwards et al., 2013; Helgason et al., 2008; Knowles et al., 2016). Pass all material properties to the FE meshes, storing them in the nodes and elements of the corresponding masks. After setting up node sets, elements, and material properties in Simpleware™, create the mesh and export it as a. inp file.
Only descriptive statistics (means, standard deviations) were reported; no inferential testing was performed due to the small sample size.
2.1 Simulation workflow
Finite element analysis was conducted using Abaqus/CAE 2020, Abaqus™ (Dassault Systems, Velizy-Villacoublay, France). The input file generated in Simpleware™, which contained the finite element mesh, material properties, and predefined node sets (distal and proximal), was imported into the software. To simulate fixation, a displacement/rotation boundary condition was applied to the distal node set, fully constraining it in all three translational directions (U1, U2, and U3). Time-dependent knee joint forces derived from a musculoskeletal simulation in AnyBody™ were implemented in Abaqus using tabular amplitude definitions. Separate amplitude curves were created for each force component (Fx, Fy, and Fz). These forces were applied as concentrated loads to the proximal node set, each aligned with the corresponding spatial direction. Next, a static analysis step was defined with a time period matching that of the applied amplitudes. Nonlinear effects were disabled, and appropriate incrementation settings were selected to ensure stable convergence. For postprocessing, field output requests included the stress tensor (S), von Mises stress (MISES), strain tensor (E), displacement (U), reaction force (RF), and integration point volume (IVOL). Finally, a simulation job was created and submitted for execution. The simulation follows linear elastic assumptions.
For the simulation of Interfragmentary movement, we focus on a mechanobiological approach based on mechanical stimuli, as discussed in the works of Claes et al. (1998), Shefelbine et al. (2005). Specifically, we examine the interaction between octahedral shear strain, which is derived from the deviatoric component of the strain tensor related to shape distortion, and volumetric strain, which is associated with volume changes and hydrostatic pressure (Ghiasi et al., 2017). This approach emphasizes the significance of mechanical deviatoric strains in influencing cell differentiation and tissue formation. This notion aligns well with experimental findings (Ito et al., 2003; Doblaré et al., 2004; Richardson et al., 2003). To precisely analyze the mechanics within the fracture gap, we read the strain tensor for each tetrahedral element in that region. Both octahedral shear strain and volumetric strain were calculated and assessed according to the parameters established by Claes et al. (1998), Shefelbine et al. (2005). A MATLAB script is used to convert ODB (Output Database) files, produced by Abaqus, to store simulation results, into VTK (Visualization Toolkit) files, enabling improved visualization and analysis through advanced post-processing tools. Further information about the simulation process can be found in the Supplementary Material.
2.2 Case data
2.2.1 Use case 1
The first case is a 57-year-old patient suffering from a closed proximal tibia fracture (AO/OTA 41 C1.2; Schatzker IV) after falling from a ladder. Initially treated at an external institution, the patient was introduced to our department with a posteromedial malreduction, resulting in an unreduced knee joint with posteromedial dislocation and a bone non-union. The patient underwent revision with an open reduction and fixation using a medial, posteromedial, and anterolateral locking plate, along with a medial tricortical allograft. Postoperatively, a 20 kg partial weight-bearing order was given. The patient experienced an uneventful healing course, regaining full weight-bearing capacity and demonstrating timely fracture consolidation, with high satisfaction.
2.2.2 Use case 2
The second case is of a 72-year-old patient with atrial fibrillation, hypertension, and hypothyroidism, who sustained a closed proximal tibia fracture in a fall from standing height (AO 41 B3.3h; Schatzker IV). She was treated with an open reduction and medial locking plate osteosynthetic fixation. Postoperatively, a 20 kg partial weight-bearing order was given. Initially, the patient exhibited delayed healing, with early screw-back failure at the 6-week follow-up. Upon further review of the history, while the patient attempted to comply with the weight-bearing recommendation, she reported difficulties maintaining balance, which may have contributed to an overload of the fracture situation. At this point, she was instructed to limit weight-bearing for another 6 weeks with only gradual increases under strict physical therapy guidance. The fracture then showed slow but steady consolidation over the next 6 months. The screw progressively backed out and was removed 3 months after the initial treatment.
2.2.3 Use case 3
Case three is a 52-year-old patient who sustained a closed proximal tibia fracture in a skiing accident (AO 41 C3.1e; Schatzker IV). She was treated with a medial, posteromedial, and anterolateral locking plate with an open reduction utilizing a lateral epicondylar osteotomy. Postoperatively, a 20 kg partial weight-bearing order was given. The patient exhibited an uneventful healing course, with timely fracture consolidation and a favorable clinical outcome. Accordingly, the simulation predicted an uneventful healing course during all weight-bearing levels.
For all three patients, the clinical prescription was an immediate postoperative partial weight-bearing of 20 kg. Higher load levels (50% body weight and full weight-bearing) were not prescribed clinically but were tested only virtually in the simulation model.
3 Results
This section first presents the results of the diagrams of the ground reaction forces during the TUG test, focusing on compliance with partial weight-bearing recommendations in healthy test subjects. The derived mean force curve of the knee joint, based on the analyzed AnyBody™ data, is then examined. This assessment facilitates a deeper understanding of the impact of partial weight-bearing on joint forces and, consequently, its influence on interfragmentary motion in patients. We then examine patient outcomes and the effects of different loading scenarios on design stiffness and interfragmentary motion. Here, the analysis focuses on the mechanical effects of different partial weight-bearing approaches and evaluates their effectiveness in promoting healing. This analysis contributes to a more nuanced understanding of the effectiveness of orthopedic treatments and bridges the gap between theoretical modeling and clinical application.
3.1 Healthy subjects
A total of 22 participants underwent testing, comprising 11 women and 11 men. The average age of the participants was 26 ± 4 years. The mean weight is 71 ± 15 kg, and the average height is 176.5 ± 11 cm. The testing protocol involved each participant undergoing 11 walking trials with 50% of their body weight, followed by ten repetitions with a partial load of 20 kg. Figure 3 illustrates a representative instance of the ground reaction forces of the left and right feet obtained from insole measurements during sessions with healthy test subjects. The results of the mean knee joint forces, which were standardized to the body weight of the test subjects for comparison purposes, are illustrated in detail in Figure 4. It shows the mean knee joint force as a percentage of body weight during a gait cycle, initiated and completed by heel strike, with a partial load of 20 kg or 50% of the subjects’ body weight. The testing commenced with recording the TUG task during normal walking (Figure 4A). Subsequently, the subject walked with crutches and underwent testing for partial weight-bearing using a household scale. The subject underwent repetitions of partial weight-bearing trials of 20 kg on the right foot, consistent with standard clinical protocols (Figure 4B). Following these trials, the TUG task was repeated under the prescribed partial weight-bearing conditions, evaluating the ground reaction forces associated with this training regimen (Figure 4C). A more detailed representation of the knee joint forces can be seen in Supplementary Figure S6. Therefore, partial weight-bearing was implemented in the AnyBody™ system, and inverse dynamics were performed to infer the respective joint forces. The implementation utilized the measured ground reaction force from Moticon™ in the gait model, along with the corresponding gait data from the Xsens™ measurement. This detailed analysis reveals the nuanced impact of partial load on knee joint forces. However, the differences between the test subjects, particularly at 50% of body weight, demonstrate the individual nature of the effects of partial loading on knee joint forces. An example is the stark contrast between the impact of 50 percent partial weight-bearing on subjects with a body weight of 50 kg and those with a body weight of 90 kg. Furthermore, it is crucial to recognize that a partial load of 20 kg can increase knee joint forces to a peak of 100% of a person’s body weight, highlighting the significant impact that even relatively small external loads can have on knee joint biomechanical loads.

Figure 3. Ground reaction force curves for (A) full weight-bearing, (B) Training with partial weight-bearing and crutches, and (C) partial weight-bearing reprehensive for one subject.

Figure 4. Mean knee joint force curve over all participants (n = 22) for (A) 20 kg weight-bearing, (B) 50% body weight, and (C) full weight-bearing and natural walking.
3.2 Patient-specific simulations
To investigate the effects of partial weight-bearing on different fracture scenarios and treatment modalities, digital twins representing the three patients were subjected to finite element simulation. These simulations replicate the biomechanical effects of various loading conditions on the patient’s lower extremities, simulating a gait cycle with partial loading of 20 kg, 50% of the patient’s body weight, and a full weight-bearing scenario. Clinically, all patients were limited to 20 kg partial weight-bearing postoperatively; the 50% body weight and full weight-bearing conditions were included solely as virtual test scenarios in the simulation. The methodology for applying these loads graphically, as shown in Figure 1, provides a clear visual reference for the biomechanical conditions applied. With the boundary conditions of the actual patient measurement, comparable results are obtained; however, these are limited by the variable testing of the partial load. Figure 5 presents the outcome of these simulations, with the data for the three patients arranged horizontally from use cases 1 to 3. In particular, the analysis considers the stresses and strains in the fracture gap. It interprets these variables through the bone healing window (Braun et al., 2021; Lutz, 2021; Lutz, 2017; Shefelbine et al., 2005) by linking local stresses and strains in the fracture gap, as described by Claes and Heigele (1999) (Figure 5D) - a concept that relates the mechanical environment of a fracture site to the potential for bone healing. Reducing interfragmentary motion at full weight-bearing improves fracture stability. In addition, the simulation output von Mises stress analyzed stresses within the plate construction. In all three patient cases, the maximum von Mises stresses were below the material-specific yield strength of medical titanium aluminum alloys (Geetha et al., 2009), indicating sufficient construct stability.
• Case #1: The simulation results predicted an uneventful healing course after the revision of all weight-bearing levels. The maximum value is below 400 MPa with uniform stress distribution, indicating increased stability when supplied with three plates.
• Case #2: High local shear stresses outside the defined healing window under full weight-bearing. The full load caused locally increased stresses (>600 MPa) in the area of the distal screws of the medial plate.
• Case #3: The simulation predicted an uneventful healing course at all weight-bearing levels. The stress values in all plates were consistently below 500 MPa.

Figure 5. Results of the biomechanical simulation of digital twins for use cases 1, 2, and 3 under three different loading conditions: (A) 20 kg partial weight-bearing, (B) 50% body weight, and (C) full weight-bearing. Additionally, the implant stress in MPa is included for the different loading conditions. (D) Color pattern of the healing window, illustrating the quantity and precise position of the elements in the various healing states.
3.3 Clinical outcomes
Clinical outcomes aligned with the simulation predictions:
• Case #1: Healed uneventfully, consistent with patient-specific simulation outcomes.
• Case #2: Clinically observed retraction of a screw indicates a mechanically critical situation, which clinicians addressed by adjusting the post-treatment. The course went on to fracture consolidation only after adaptation to her weight-bearing and physical therapy training.
• Case #3: Healing occurred without major complications. The results of the limit simulation under higher loads corresponded to the cautious clinical progression of the load required for this case.
4 Discussion
The presented results demonstrate a simulation approach for determining a safe range of postoperative weight-bearing in individual patients after tibial plateau fractures, considering construct stability and interfragmentary strain that influence fracture healing. Notably, this model integrates the patient’s individual fracture and treatment situation with weight-adjusted movement norm data from the uninjured reference group, enabling the safe testing of different weight-bearing scenarios without exposing the patient to potentially harmful situations. Furthermore, construct stability is ensured for every patient by analyzing the von Mises stress, which remains below the critical material limits (i.e., yield strength of ∼880 MPa) for titanium aluminum alloys used in medical-grade implants (Geetha et al., 2009). This guarantees that no threshold values are exceeded and the components are not subjected to excessive stress. In addition, the analysis shows that by increasing the number of implants, patients 1 and 3 had a mechanical stimulus in the lower part of the ranges expected for uneventful bone healing, indicating a favorable biomechanical environment for fracture healing. This improvement in stability is particularly beneficial for bone healing, as it reduces the mechanical stimulus at the fracture site and brings it into the optimal range of the healing window.
The feasibility and validity of the presented simulation approach are grounded in a series of previous methodological and clinical studies that have established reliable workflows for biomechanical fracture simulation. Braun et al. (2021) demonstrated the integration of motion capture in tibial pseudarthrosis treatment, introducing a patient-specific approach to assess fracture mechanics. Similarly, the automated CT-based finite element pipeline by Dahmen et al. (2015) enabled individualized simulation of bone-implant constructs and laid the groundwork for reproducible, image-based simulation methods. These technical principles were further refined by Roland et al. (2024), who proposed a direct algorithmic strategy for simulating bone healing based on CT imaging, and by Orth et al. (2023), who used simulations to assess the effects of motion and weight-bearing on fracture healing in the lower leg. The applicability of the measurement systems used in cases of tibial fractures was 21 demonstrated by Andres et al. (2025a) through weekly measurements during the healing phase. The current study extends these approaches by incorporating individualized musculoskeletal loading, derived from inertial motion capture and force-sensing insoles, and fracture-specific finite element modeling to create digital twins of postoperative patients. This enables the realistic simulation of weight-bearing scenarios without exposing patients to biomechanical risk, aligning with earlier feasibility studies that utilized dynamic in vivo load measurements and simulation-based analysis of fracture healing (Andres et al., 2025b; Braun et al., 2022). Concerning validation, the predictive power of the simulation model is supported by agreement with clinical outcomes in all three use cases, and notably, in case #2, simulation results under full weight-bearing conditions predicted localized overloading and interfragmentary shear strain outside the safe healing window, which corresponded to screw loosening and delayed consolidation observed in the clinical course. We acknowledge that delayed consolidation in this patient was multifactorial, influenced by age, comorbidities, and surgical construct, in addition to biomechanical factors. Nonetheless, the simulation highlighted overload at the fracture site consistent with the clinical finding of screw loosening, indicating that while the model cannot exclude non-mechanical contributors, it reliably identified a key mechanical risk factor. These findings are consistent with the mechanobiological principles of Claes and Heigele (1999) and Shefelbine et al. (2005), which establish threshold values for strain-driven tissue differentiation during bone healing. Moreover, recent experimental validation of interfragmentary motion in bone-implant systems, such as by Wickert et al. (2024), reinforces the technical accuracy of the FEM-based predictions. In addition, the use of validated software environments (e.g., AnyBody™ and Abaqus™) and standardized imaging-to-simulation pipelines (Augat et al., 2021) further enhances the reproducibility and transferability of the approach. The current model builds upon and extends a robust methodological foundation, demonstrating technical feasibility and clinical validity in predicting mechanical outcomes relevant to postoperative fracture care.
The recording of joint forces in healthy test subjects enables a differentiated understanding of the biomechanical effects of partial loading during walking. Detailed measurements of knee joint forces can be found in the OrthoLoad™ database (Heinlein et al., 2009). Nine patients were measured using an instrumented knee joint implant while performing various activities. However, the database only provides two data sets, k5r and k2l, during a 3-point walk. This innovative approach circumvents the limitations associated with the need for measurable implants, giving valuable insights into the effects of loading during various walking scenarios.
No commercially available intramedullary implants exist to determine the simulated parameters, making it impossible to validate the determined mechanical conditions experimentally. However, the general capability of fracture condition simulation has been previously demonstrated (Andres et al., 2025c; Ghiasi et al., 2017), and clinical validation through the healed fracture supports its efficacy. Furthermore, the application of a postoperative, simulated biomechanical healing control has been described in a case series and subsequently validated through animal experiments (Dailey et al., 2019; Schwarzenberg et al., 2021). Additionally, research has shown that finite element simulations of bone-implant systems can be effectively validated using experimental data, reinforcing the accuracy of these simulations in replicating physical behaviors. Our recently published paper validated our simulation through biomechanical experiments on human specimens, and the results are consistent with those from the experiments (Wickert et al., 2024).
In the aftercare of tibial plateau fractures, a primary focus with potential implications for the clinical and radiographic treatment outcomes is postoperative weight-bearing (Arnold et al., 2017; Canton et al., 2022; Thewlis et al., 2015; Williamson et al., 2018). Traditionally, long periods of non-weight bearing have been prescribed for fear of secondary fracture dislocation, intra-articular step of formation, and, ultimately, a higher risk of posttraumatic arthritis. The current surveys, as well as clinical analysis, have shown, however, that these strict non-weight-bearing regimes are increasingly being abandoned in clinical practice, and that early, permissive weight-bearing in select cases is not associated with a higher degree of complications (Arnold et al., 2017; Thewlis et al., 2015). However, a low-powered study directly measuring postoperative joint subsidence through radio stereometrics analysis has shown that higher average knee joint loading during the stance phase of gait is not associated with increased postoperative fracture fragment migration beyond 3 mm (Thewlis et al., 2015). The clinical significance of these increases in step-off formation, particularly in relation to long-term outcomes, remains unclear. The relevance of an individual simulation approach in determining this risk for secondary fracture migration and loss of reduction is evident in the second patient with medial-sided plating. This is the only patient simulated in this analysis with postoperative fracture subsidence and subsequent loosening of a screw in the plate apex, experiencing an outward pushing force running almost parallel to the fracture plane. This patient reported being permissively weight-bearing within the first weeks after the fracture and used no aids for walking, thus exposing her to at least her full body weight. Accordingly, the simulation of the case under full weight-bearing conditions showed the highest amount of fracture area, above safe ranges of interfragmentary motion according to the healing window (Lutz, 2017), confirming the simulation’s principal capability to assess situations at risk, already seen in studies (Braun et al., 2021; Braun et al., 2019). Interestingly, the fracture in this patient consolidated in this position and has led to a satisfactory result for the patient after screw removal and prescription of partial weight-bearing for 6 weeks.
The other simulated fracture situations with triple plate osteosynthesis have shown sufficient stability and safe ranges of interfragmentary motion throughout all weight-bearing levels. Accordingly, these fractures exhibited uneventful postoperative treatment courses, with no secondary loss of reduction or hardware failures. Both patients, however, required more than 8 weeks to gradually wean themselves off their walking aids. This aligns with current clinical studies and review analyses focusing on permissive weight-bearing. While it is safe to allow permissive weight-bearing, patients often require considerable time to regain full weight-bearing (Pishtiwan Hassan Shaker Kalmet et al., 2019). Simulation workflows, such as the one introduced here, can help identify patients who are safe for higher degrees of weight-bearing early during the treatment course and potentially increase their physical therapy training within their simulated safe range, guiding them toward earlier functional recovery (Braun et al., 2021; Braun et al., 2017).
The approach demonstrates the value of a simulation workflow in testing the weight-bearing associated effects on construct stability, joint reaction forces, and fracture healing. Identifying patients at risk for treatment failure early in the fracture course is a potentially significant factor in reducing the patient and socioeconomic burden of disease. On the patient’s side, early interventions, ranging from adapting the aftercare regimen to surgical treatment, are possible. Aftercare adaptations can reduce weight-bearing in situations with increased interfragmentary motion or the risk of fixation failure. In contrast, the frequency and intensity of full weight-bearing can be addressed in understimulated situations. Furthermore, individual, simulation-adjusted adaptation of ROM braces could prevent certain critical motions. If the simulation determines that healing cannot be achieved by modulation of the weight-bearing input, the stiffness of the osteosynthetic construct could be adjusted, i.e., through augmentative plating, as well as further biologic measures (autologous bone graft) in mechanically sound situations without clinical healing. By establishing this simulation as a rehabilitation and screening tool, the added socioeconomic effects can be expected through better healing rates through adapted therapy and earlier intervention in critical cases. Research on joint reaction forces in arthroplasty has shown that it is not only gait and associated weight-bearing that are putting stress on the proximal tibia, but also different movement situations that can exert high levels of force on a treated joint (Berg et al., 2014). The presented workflow and the database of reference movements can be extended to include more scenarios, such as getting out of bed, walking at different speeds, or rising from a chair, to provide a more comprehensive simulation of individual patients’ movement patterns. Furthermore, this would allow us to test the particular patient in more detail and at a higher level of injury specificity, such as high-degree flexion testing in flexion-type injuries, to adapt postoperative range of motion protocols.
Simulation approaches, such as the one introduced in this article, offer a unique opportunity to enhance our understanding of the effects of postoperative weight-bearing on osteosynthetic construct stability, fracture healing, and, ultimately, patient outcomes. Over the past 30 years, the discussion on potential weight-bearing recommendations following operative fracture treatment of the proximal tibia has remained active and controversial (Canton et al., 2022; Kubiak et al., 2013; Whitelaw et al., 1993). Determining the biomechanical effects of weight-bearing on joint mechanics without computer simulation requires both controlled weight-bearing and precise radiographic measures of joint subsidence and fracture movement early during the fracture course. Limited patient number studies combining both these techniques have been performed showing that increased amounts of weight-bearing after proximal tibia fractures can be allowed without negative radiographic and clinical consequences (Thewlis et al., 2015; The et al., 2017), but to get to a comprehensive recommendation on all types of proximal tibia fractures, while considering the many treatment and hardware placement options is not feasible with these highly complex and research-oriented techniques. Clinical studies recommending permissive weight-bearing in these patients have shown that this can be associated with a decreased short-term recovery time, but also with limitations to the amount of weight-bearing, especially at the study endpoint, when the patient converts to a total knee arthroplasty (Pishtiwan Hassan Shaker Kalmet et al., 2019). To safely test these limitations early during the fracture treatment on the individual treatment situation, but safely with simulated, different weight-bearing levels, thus not only offers the opportunity to increase our understanding of safe boundaries during fracture rehabilitation but actually can allow for individualized aftercare concerning weight-bearing with wearable techniques, in situations that would otherwise be hard to assess and control (Braun et al., 2017; Marmor et al., 2022). The clinical feasibility of our approach has been demonstrated, and it can now be further refined in larger-scale clinical studies. The model’s outcome predictions enable the development of tailored post-operative care plans, ensuring high-risk patients receive the necessary attention to prevent complications. Hospitals can allocate resources more efficiently, reducing healthcare costs by preventing complications and readmissions. Predictive insights enable clinicians to implement early intervention strategies, thereby improving patient outcomes. Additionally, the findings can inform enhanced recovery after surgery protocols, optimizing recovery times and reducing hospital stays.
5 Limitations
While the study provides insightful results, it also offers opportunities for further research and refinement. Including a relatively small number of patients provides a focused initial investigation, paving the way for future studies with larger cohorts to enhance the robustness and applicability of the results. While the simulations do not capture the full biological context of the injury, including concomitant soft tissue damage and the effects of the ligamentous injury on joint stability, they open the door to more comprehensive modeling that incorporates these factors. In addition, the study highlights an important area for further investigation: determining the precise level of mechanical stimuli within the fracture zone that is optimal for healing.
The inclusion of a relatively small number of patients enables an initial analysis and serves as preparation for future studies with larger cohorts in order to improve the validity and applicability of the results. In addition, no formal power analysis was performed, as the study was designed as a feasibility and proof-of-concept investigation. The cohort of the uninjured reference group consisted of young individuals without comorbidities, which does not correspond to the typical patient group with tibial plateau fractures.
Although the simulation predicted biomechanical overload in Case #2 consistent with screw loosening, delayed healing in this patient was likely multifactorial. Age, comorbidities, balance deficits, and construct choice all contributed, and the model cannot fully isolate mechanical loading as the sole determinant. Moreover, clinical validation was limited to retrospective case comparison, which may invite confirmation bias; prospective validation in larger patient populations is required.
No sensitivity analysis was performed beyond this. While the simulations do not capture the full biological context of the injury—including concomitant soft tissue damage, vascularity, or the effects of ligamentous injury on joint stability—they open the door to more comprehensive modeling that incorporates these factors. This is a feasibility study to demonstrate the technical plausibility and clinical relevance of combining motion capture, musculoskeletal simulation, and finite element modeling.
6 Conclusion
This study introduces a patient-specific simulation framework for evaluating postoperative weight-bearing after tibial plateau fractures. By integrating patient imaging with reference motion data, the approach enables safe, non-invasive testing of individual loading scenarios and has the potential to inform personalized rehabilitation protocols. Beyond walking, the model can be extended to incorporate other everyday movements such as stair climbing or rising from a chair, broadening its clinical utility. Future research with larger cohorts will be essential to refine predictive accuracy and establish its value as a clinical decision-support tool. Ultimately, the goal is to define safe postoperative boundary conditions for loading and range of motion, supporting fracture healing while preserving implant integrity.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by Institutional Ethics Committee of the University of Tuebingen (Project Number 318/2022BO2). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
AA: Conceptualization, Data curation, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review and editing. MR: Investigation, Project administration, Supervision, Writing – review and editing. KW: Data curation, Writing – review and editing. SD: Funding acquisition, Project administration, Resources, Writing – review and editing. DT: Writing – review and editing. TH: Funding acquisition, Project administration, Resources, Writing – review and editing. BB: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Writing – original draft, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This study was funded by the German Federal Ministry for Science and Education (BMBF) under the grant “VirtuS” (FKZ: 13GW0572).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2025.1659029/full#supplementary-material
References
Andersen, M. S., Damsgaard, M., and Rasmussen, J. (2009). Kinematic analysis of over-determinate biomechanical systems. Comput. Methods Biomech. Biomed. Engin. 12 (4), 371–384. doi:10.1080/10255840802459412
Andres, A., Roland, M., Orth, M., and Diebels, S. (2025a). From injury to full recovery: monitoring patient progress through advanced sensor and motion capture technology. Sensors (Basel). 1, 3853–11. doi:10.3390/s25133853
Andres, A., Roland, M., Wickert, K., Diebels, S., Stöckl, J., Herrmann, S., et al. (2025b). Advantages of digital twin technology in orthopedic trauma surgery – exploring different clinical use cases. Sci. Rep. 15 (1), 19987–13. doi:10.1038/s41598-025-04792-w
Andres, A., Roland, M., Wickert, K., Pohlemann, T., Orth, M., Diebels, S., et al. (2025c). Individual postoperative and preoperative workflow for patients with fractures of the lower extremities. Clin. Biomech., 106503. doi:10.1016/j.clinbiomech.2025.106503
Arnold, J. B., Tu, C. G., Phan, T. M., Rickman, M., Varghese, V. D., Thewlis, D., et al. (2017). Characteristics of postoperative weight bearing and management protocols for tibial Plateau fractures: findings from a scoping review. Injury 48 (12), 2634–2642. doi:10.1016/j.injury.2017.10.040
Augat, P., Hast, M. W., Schemitsch, G., Heyland, M., Trepczynski, A., Borgiani, E., et al. (2021). Biomechanical models: key considerations in study design. OTA Int. 4 (2), E099. doi:10.1097/OI9.0000000000000099
Bergmann, G., Bender, A., Graichen, F., Dymke, J., Rohlmann, A., Trepczynski, A., et al. (2014). Standardized loads acting in knee implants. PLoS One 9, e86035. doi:10.1371/journal.pone.0086035
Braun, B. J., Veith, N. T., Rollmann, M., Orth, M., Fritz, T., Herath, S. C., et al. (2017). Weight-bearing recommendations after operative fracture treatment—fact or fiction? Gait results with and feasibility of a dynamic, continuous pedobarography insole. Int. Orthop. 41 (8), 1507–1512. doi:10.1007/s00264-017-3481-7
Braun, B. J., Pohlemann, T., Herath, S. C., Klein, M., Rollmann, M. F., Derr, R., et al. (2019). An individualized simulation model based on continuous, independent, ground force measurements after intramedullary stabilization of a tibia fracture. Arch. Appl. Mech. 89 (11), 2351–2360. doi:10.1007/s00419-019-01582-5
Braun, B. J., Orth, M., Diebels, S., Wickert, K., Andres, A., Gawlitza, J., et al. (2021). Individualized determination of the mechanical fracture environment after tibial exchange nailing—A simulation-based feasibility study. Front. Surg. 8 (September), 749209–9. doi:10.3389/fsurg.2021.749209
Braun, B. J., Histing, T., Herath, S. C., Rollmann, M. F. R., Reumann, M., Menger, M. M., et al. (2022). Bewegungsanalyse und muskuloskeletale Simulation in der Pseudarthrosentherapie – Erfahrungen und erste klinische Ergebnisse. Die Unfallchirurgie 125 (8), 619–627. doi:10.1007/s00113-022-01208-6
Brent Edwards, W., Schnitzer, T. J., and Troy, K. L. (2013). Torsional stiffness and strength of the proximal tibia are better predicted by finite element models than DXA or QCT. J. Biomech. 46 (10), 1655–1662. doi:10.1016/j.jbiomech.2013.04.016
Burns, G. T., Zendler, J. D., and Zernicke, R. F. (2019). Validation of a wireless shoe insole for ground reaction force measurement. J. Sports Sci. 37 (10), 1129–1138. doi:10.1080/02640414.2018.1545515
Canton, G., Sborgia, A., Dussi, M., Rasio, N., and Murena, L. (2022). Early weight bearing in tibial Plateau fractures treated with ORIF: a systematic review of literature. J. Orthop. Surg. Res. 17 (1), 261–268. doi:10.1186/s13018-022-03156-8
Carrera, I., Gelber, P. E., Chary, G., González-Ballester, M. A., Monllau, J. C., and Noailly, J. (2016). Fixation of a split fracture of the lateral tibial plateau with a locking screw plate instead of cannulated screws would allow early weight bearing: a computational exploration. Int. Orthop. 40 (10), 2163–2169. doi:10.1007/s00264-015-3106-y
Carrera, I., Gelber, P. E., Chary, G., Gomez Masdeu, M., González Ballester, M. A., Monllau, J. C., et al. (2018). An intact fibula may contribute to allow early weight bearing in surgically treated tibial plateau fractures. Knee Surg. Sport. Traumatol. Arthrosc. 26 (3), 756–761. doi:10.1007/s00167-017-4428-7
Cattaneo, P. M., Dalstra, M., and Frich, L. H. (2001). A three-dimensional finite element model from computed tomography data: a semi-automated method. Proc. Inst. Mech. Eng. Part H. J. Eng. Med. 215 (2), 203–212. doi:10.1243/0954411011533760
Chiodo, C. P., Macaulay, A. A., Palms, D. A., Smith, J. T., and Bluman, E. M. (2016). Patient compliance with postoperative lower-extremity non-weight-bearing restrictions. J. Bone Jt. Surg. - Am. 98 (18), 1563–1567. doi:10.2106/JBJS.15.01054
Claes, L. E., and Heigele, C. A. (1999). Magnitudes of local stress and strain along bony surfaces predict the course and type of fracture healing. J. Biomech. 32 (3), 255–266. doi:10.1016/S0021-9290(98)00153-5
Claes, L. E., Heigele, C. A., Neidlinger-Wilke, C., Kaspar, D., Seidl, W., Margevicius, K. J., et al. (1998). Effects of mechanical factors on the fracture healing process. Clin. Orthop. Relat. Res. 355, S132–S147. doi:10.1097/00003086-199810001-00015
Court-Brown, C. M., and Caesar, B. (2006). Epidemiology of adult fractures: a review. Injury 37 (8), 691–697. doi:10.1016/j.injury.2006.04.130
Cunningham, B. P., Brazina, S., Morshed, S., and Miclau, T. (2017). Fracture healing: a review of clinical, imaging and laboratory diagnostic options. Injury 48, S69–S75. doi:10.1016/j.injury.2017.04.020
Dahmen, T., Roland, M., Tjardes, T., Bouillon, B., Slusallek, P., and Diebels, S. (2015). An automated workflow for the biomechanical simulation of a tibia with implant using computed tomography and the finite element method. Comput. Math. Appl. 70 (5), 903–916. doi:10.1016/j.camwa.2015.06.009
Dailey, H. L., Schwarzenberg, P., Daly, C. J., Boran, S. A. M., Maher, M. M., and Harty, J. A. (2019). Virtual mechanical testing based on low-dose computed tomography scans for tibial fracture: a pilot study of prediction of time to union and comparison with subjective outcomes scoring. J. Bone Jt. Surg. 101 (13), 1193–1202. doi:10.2106/JBJS.18.01139
Delp, S. L., Anderson, F. C., Arnold, A. S., Loan, P., Habib, A., John, C. T., et al. (2007). OpenSim: open-source software to create and analyze dynamic simulations of movement. IEEE Trans. Biomed. Eng. 54 (11), 1940–1950. doi:10.1109/TBME.2007.901024
Doblaré, M., García, J. M., and Gómez, M. J. (2004). Modelling bone tissue fracture and healing: a review. Eng. Fract. Mech. 71, 1809–1840. doi:10.1016/j.engfracmech.2003.08.003
Erick, M. H. L. G., Tang, A., Vialonga, M., Jankowski, J. M., Liporace, F. A., Yoon, R. S., et al. (2023). Early versus standard weight bearing following operative treatment of tibial Plateau fractures: do we really have to wait So long? J. Knee Surg. 37 (05), 402–408. doi:10.1055/s-0043-1772232
Gaston, P., Will, E. M., and Keating, J. F. (2005). Recovery of knee function following fracture of the tibial plateau. J. Bone Jt. Surg. - Ser. B 87 (9), 1233–1236. doi:10.1302/0301-620X.87B9.16276
Geetha, M., Singh, A. K., Asokamani, R., and Gogia, A. K. (2009). Ti based biomaterials, the ultimate choice for orthopaedic implants - a review. Prog. Mater. Sci. 54 (3), 397–425. doi:10.1016/j.pmatsci.2008.06.004
Ghiasi, M. S., Chen, J., Vaziri, A., Rodriguez, E. K., and Nazarian, A. (2017). Bone fracture healing in mechanobiological modeling: a review of principles and methods. Bone Rep. 6, 87–100. doi:10.1016/j.bonr.2017.03.002
Heinlein, B., Kutzner, I., Graichen, F., Bender, A., Rohlmann, A., Halder, A. M., et al. (2009). ESB clinical biomechanics award 2008: complete data of total knee replacement loading for level walking and stair climbing measured in vivo with a follow-up of 6-10 months. Clin. Biomech. 24 (4), 315–326. doi:10.1016/j.clinbiomech.2009.01.011
Helgason, B., Perilli, E., Schileo, E., Taddei, F., Brynjólfsson, S., and Viceconti, M. (2008). Mathematical relationships between bone density and mechanical properties: a literature review. Clin. Biomech. 23 (2), 135–146. doi:10.1016/j.clinbiomech.2007.08.024
Hvid, I., Bentzen, S. M., Linde, F., Mosekilde, L., and Pongsoipetch, B. (1989). X-ray quantitative computed tomography: the relations to physical properties of proximal tibial trabecular bone specimens. J. Biomech. 22 (8–9), 837–844. doi:10.1016/0021-9290(89)90067-5
Ibrahim, M. A.-a., Moustafa, M. M. A., Michel Brismée, J., Farouk, O., Kamal Abdelnasser, M., Said, H. G., et al. (2025). Immediate weight-bearing after tibial plateau fractures internal fixation results in better clinical outcomes with similar radiological outcomes: a randomized clinical trial. Int. Orthop. 49 (5), 1245–1253. doi:10.1007/s00264-025-06443-1
Iliopoulos, E., and Galanis, N. (2020). Physiotherapy after tibial plateau fracture fixation: a systematic review of the literature. SAGE Open Med. 8, 2050312120965316. doi:10.1177/2050312120965316
Imam, A. C., and Fraker, M. A. (1996). Titanium alloys as implant materials. Mater. Sci., 3–16. doi:10.1520/STP16066S
Ito, K., Bishop, N. E., Tami, I., and Schneider, E. (2003). In vivo comparison of early fracture healing under deviatoric and volumetric deformationsitle. Acta Bioeng. Biomech. 4, 754–755.
Ibrahim, N. I, Ahmad, M. S., Zulfarina, M. S., Zaris, S. N. A. S. M., Zainal, A. N. A., Isa, N. M., et al. (2019). Physical function assessment of older adults with lower body fractures at 3 months post-discharge from hospital. Ther. Clin. Risk Manag. 15, 201–210. doi:10.2147/TCRM.S189748
Kammerlander, C., Pfeufer, D., Lisitano, L. A., Mehaffey, S., Böcker, W., and Neuerburg, C. (2018). Inability of older adult patients with hip fracture to maintain postoperative weight-bearing restrictions. J. Bone Jt. Surg. - Am. 100 (11), 936–941. doi:10.2106/JBJS.17.01222
Knowles, N. K., Reeves, J. M., and Ferreira, L. M. (2016). Quantitative computed tomography (QCT) derived bone mineral density (BMD) in finite element studies: a review of the literature. J. Exp. Orthop. 3, 36. doi:10.1186/s40634-016-0072-2
Kristensen, M. T., Foss, N. B., and Kehlet, H. (2007). Timed “Up and Go” test as a predictor of falls within 6 months after hip. Fract. Surg. 87, 1. doi:10.2522/ptj.20050271
Kubiak, E. N., Beebe, M. J., North, K., Hitchcock, R., and Potter, M. Q. (2013). Early weight bearing after lower extremity fractures in adults. J. Am. Acad. Orthop. Surg. 21 (12), 727–738. doi:10.5435/JAAOS-21-12-727
Lutz, E. C. (2017). “Mechanobiology of fracture healing,” in From basic science to clinical application. Available online at: https://link.springer.com/bookseries/10280.
Lutz, C. (2021). Improvement of clinical fracture healing – what can be learned from mechano-biological research? J. Biomech. 115, 110148. doi:10.1016/j.jbiomech.2020.110148
Marmor, M. T., Dailey, H., Marcucio, R., and Hunt, A. C. (2020). Biomedical research models in the science of fracture healing - pitfalls and promises. Injury 51 (10), 2118–2128. doi:10.1016/j.injury.2020.06.025
Marmor, M. T., Grimm, B., Hanflik, A. M., Richter, P. H., Sivananthan, S., Yarboro, S. R., et al. (2022). Use of wearable technology to measure activity in orthopaedic trauma patients: a systematic review. Indian J. Orthop. 56 (7), 1112–1122. doi:10.1007/s43465-022-00629-0
Mcnamara, I. R., Smith, T. O., Shepherd, K. L., Clark, A. B., Nielsen, D. M., Simon, D., et al. (2015). Surgical fixation methods for tibial Plateau fractures. Cochrane Database Syst. Rev. 2015, CD009679. doi:10.1002/14651858.CD009679.pub2
Muro-de-la-Herran, A., García-Zapirain, B., and Méndez-Zorrilla, A. (2014). Gait analysis methods: an overview of wearable and non-wearable systems, highlighting clinical applications. Sensors Switz. 14 (2), 3362–3394. doi:10.3390/s140203362
Neal, S., Sahakian, V., and Dawson, M. (2017). Outcome after tibial Plateau fracture: how important is restoration of articular congruity? J. Orthop. Trauma 31 (3), 158–163. doi:10.1097/BOT.0000000000000762
Nir, T., Yosibash, Z., and Milgrom, C. (2009). Validation of subject-specific automated p-FE analysis of the proximal femur. J. Biomech. 42 (3), 234–241. doi:10.1016/j.jbiomech.2008.10.039
Orth, M., Ganse, B., Andres, A., Wickert, K., Warmerdam, E., Müller, M., et al. (2023). Simulation-based prediction of bone healing and treatment recommendations for lower leg fractures: effects of motion, weight-bearing and fibular mechanics. Front. Bioeng. Biotechnol. 11 (February), 1067845–13. doi:10.3389/fbioe.2023.1067845
Pishtiwan Hassan Shaker Kalmet, Van Horn, Y. Y., Sanduleanu, S., Seelen, H. A. M., Brink, P. R. G., and Poeze, M. (2019). Patient-reported quality of life and pain after permissive weight bearing in surgically treated trauma patients with tibial plateau fractures: a retrospective cohort study. Arch. Orthop. Trauma Surg. 139 (4), 483–488. doi:10.1007/s00402-018-3088-5
Podsiadlo, S., and Richardson, D. (1991). The timed “up and go”: a test of basic functional mobility for frail elderly persons. J. Am. Geriatr. Soc. 39 (2), 142–148. doi:10.1111/j.1532-5415.1991.tb01616.x
Rho, J. Y., Hobatho, M. C., and Ashman, R. B. (1995). Relations of mechanical properties to density and CT numbers in human bone. Med. Eng. Phys. 17 (5), 347–355. doi:10.1016/1350-4533(95)97314-F
Richardson, J. B., Garcia, J. M., Kuiper, J. H., and Doblare, M. (2003). A numerical model to study the mechanical influences on bone fracture healing. Acta Bioeng. Biomech. 4, 394–395.
Roland, M., Diebels, S., Wickert, K., Pohlemann, T., and Ganse, B. (2024). Finite element simulations of smart fracture plates capable of cyclic shortening and lengthening: which stroke for which fracture? Front. Bioeng. Biotechnol. 12 (July), 1420047–14. doi:10.3389/fbioe.2024.1420047
Schwarzenberg, P., Klein, K., Ferguson, S. J., von Rechenberg, B., Darwiche, S., and Dailey, H. L. (2021). Virtual mechanical tests out-perform morphometric measures for assessment of mechanical stability of fracture healing in vivo. J. Orthop. Res. 39 (4), 727–738. doi:10.1002/jor.24866
Shefelbine, S. J., Augat, P., Claes, L., and Ulrich, S. (2005). Trabecular bone fracture healing simulation with finite element analysis and fuzzy logic. J. Biomech. 38 (12), 2440–2450. doi:10.1016/j.jbiomech.2004.10.019
Thewlis, D., Fraysse, F., Callary, S. A., Verghese, V. D., Jones, C. F., Findlay, D. M., et al. (2017). Postoperative weight bearing and patient reported outcomes at one year following tibial plateau fractures. Injury 48 (7), 1650–1656. doi:10.1016/j.injury.2017.05.024
Thewlis, D., Callary, S. A., Fraysse, F., and Solomon, L. B. (2015). Peak loading during walking is not associated with fracture migration following tibial plateau fracture: a preliminary case series. J. Orthop. Res. 33 (9), 1398–1406. doi:10.1002/jor.22905
van der Vusse, M., Kalmet, P. H. S., Bastiaenen, C. H. G., van Horn, Y. Y., Brink, P. R. G., and Seelen, H. A. M. (2017). Is the AO guideline for postoperative treatment of tibial plateau fractures still decisive? A survey among orthopaedic surgeons and trauma surgeons in the Netherlands. Arch. Orthop. Trauma Surg. 137 (8), 1071–1075. doi:10.1007/s00402-017-2718-7
Whitelaw, G. P., Segal, D., Mallik, A. R., Wetzler, M. J., and Franchi, A. V. (1993)). Early weight bearing of lateral tibial plateau fractures. Clin. Orthop. Relat. Res. 294, 232–237. doi:10.1097/00003086-199309000-00032
Wickert, K., Roland, M., Andres, A., Diebels, S., Ganse, B., Kerner, D., et al. (2024). Experimental and virtual testing of bone-implant systems equipped with the AO Fracture monitor with regard to interfragmentary movement. Front. Bioeng. Biotechnol. 12 (March), 1370837–20. doi:10.3389/fbioe.2024.1370837
Williamson, M., Iliopoulos, E., Jain, A., Ebied, W., and Trompeter, A. (2018). Immediate weight bearing after plate fixation of fractures of the tibial Plateau. Injury 49 (10), 1886–1890. doi:10.1016/j.injury.2018.06.039
Yosibash, Z., Padan, R., Joskowicz, L., and Milgrom, C. (2007). A CT-based high-order finite element analysis of the human proximal femur compared to in-vitro experiments. J. Biomech. Eng. 129 (3), 297–309. doi:10.1115/1.2720906
Keywords: partial weight bearing, interfragmentary movement, construct stability, musculoskeletal simulation, proximal tibia fracture, motion capturing
Citation: Andres A, Roland M, Wickert K, Diebels S, Truhn D, Histing T and Braun B (2025) Predicting the effect of individual weight-bearing on tibial load and fracture healing after tibial plateau fractures–introduction of a biomechanical simulation model. Front. Bioeng. Biotechnol. 13:1659029. doi: 10.3389/fbioe.2025.1659029
Received: 04 July 2025; Accepted: 03 September 2025;
Published: 25 September 2025.
Edited by:
Alessandro Ruggiero, University of Salerno, ItalyReviewed by:
Ruta Jakusonoka, Riga Stradins Universitate Medicinas Fakultate, LatviaMariam Ibrahim, Assiut University Hospital, Egypt
Copyright © 2025 Andres, Roland, Wickert, Diebels, Truhn, Histing and Braun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Annchristin Andres, YW5uY2hyaXN0aW4uYW5kcmVzQHVuaS1zYWFybGFuZC5kZQ==