Aims. The presence of facet tropism has been correlated with an elevated susceptibility to lumbar disc pathology. Our objective was to evaluate the impact of facet tropism on chronic lumbosacral discogenic pain through the analysis of clinical data and finite element modelling (FEM). Methods. Retrospective analysis was conducted on clinical data, with a specific focus on the spinal units displaying facet tropism, utilizing FEM analysis for motion simulation. We studied 318 intervertebral levels in 156 patients who had undergone provocation discography. Significant predictors of clinical findings were identified by univariate and multivariate analyses. Loading conditions were applied in FEM simulations to mimic biomechanical effects on intervertebral discs, focusing on maximal displacement and intradiscal pressures, gauged through alterations in disc morphology and physical stress. Results. A total of 144 discs were categorized as ‘positive’ and 174 discs as ‘negative’ by the results of provocation discography. The presence of defined facet tropism (OR 3.451, 95% CI 1.944 to 6.126) and higher Adams classification (OR 2.172, 95% CI 1.523 to 3.097) were important predictive parameters for discography-‘positive’ discs. FEM simulations showcased uneven stress distribution and significant disc displacement in tropism-affected discs, where loading exacerbated stress on facets with greater angles. During varied positions, notably increased stress and displacement were observed in discs with tropism compared to those with normal facet structure. Conclusion. Our findings indicate that facet tropism can contribute to disc herniation and changes in intradiscal pressure, potentially exacerbating disc degeneration due to altered force distribution and increased mechanical stress. Cite this article: Bone Joint Res 2024;13(9):452–461

Aims

Precise implant positioning, tailored to individual spinopelvic biomechanics and phenotype, is paramount for stability in total hip arthroplasty (THA). Despite a few studies on instability prediction, there is a notable gap in research utilizing artificial intelligence (AI). The objective of our pilot study was to evaluate the feasibility of developing an AI algorithm tailored to individual spinopelvic mechanics and patient phenotype for predicting impingement.

Aims

The aim of this study was to determine the clinical outcomes and factors contributing to failure of transposition osteotomy of the acetabulum (TOA), a type of spherical periacetabular osteotomy, for advanced osteoarthritis secondary to hip dysplasia.

The February 2024 Wrist & Hand Roundup³⁶⁰ looks at: Occupational therapy for thumb carpometacarpal osteoarthritis?; Age and patient-reported benefits from operative management of intra-articular distal radius fractures: a meta-regression analysis; Long-term outcomes of nonsurgical treatment of thumb carpometacarpal osteoarthritis: a cohort study; Semi-occlusive dressing versus surgery in fingertip injuries: a randomized controlled trial; Re-fracture in partial union of the scaphoid waist?; The WALANT distal radius fracture: a systematic review; Endoscopic carpal tunnel release with or without hand therapy?; Ten-year trends in the level of evidence in hand surgery.

Adult Spine Deformity (ASD) is a degenerative condition of the adult spine leading to altered spine curvatures and mechanical balance. Computational approaches, like Finite Element (FE) Models have been proposed to explore the etiology or the treatment of ASD, through biomechanical simulations. However, while the personalization of the models is a cornerstone, personalized FE models are cumbersome to generate. To cover this need, we share a virtual cohort of 16807 thoracolumbar spine FE models with different spine morphologies, presented in an online user-interface platform (SpineView). To generate these models, EOS images are used, and 3D surface spine models are reconstructed. Then, a Statistical Shape Model (SSM), is built, to further adapt a FE structured mesh template for both the bone and the soft tissues of the spine, through mesh morphing. Eventually, the SSM deformation fields allow the personalization of the mean structured FE model, leading to generate FE meshes of thoracolumbar spines with different morphologies. Models can be selectively viewed and downloaded through SpineView, according to personalized user requests of specific morphologies characterized by the geometrical parameters: Pelvic Incidence; Pelvic Tilt; Sacral Slope; Lumbar Lordosis; Global Tilt; Cobb Angle; and GAP score. Data quality is assessed using visual aids, correlation analyses, heatmaps, network graphs, Anova and t-tests, and kernel density plots to compare spinopelvic parameter distributions and identify similarities and differences. Mesh quality and ranges of motion have been assessed to evaluate the quality of the FE models. This functional repository is unique to generate virtual patient cohorts in ASD.

Acknowledgements: European Commission (MSCA-TN-ETN-2020-Disc4All-955735, ERC-2021-CoG-O-Health-101044828)

Intervertebral discs (IVD) provide flexibility to the back and ensure functional distributions of the spinal loads. They are avascular, and internal diffusion-dependent metabolic transport is vital to supply nutrients to disc cells1, but interactions with personalized IVD shapes and mechanics remain poorly explored. Poromechanical finite element models of seven personalized lumbar IVD geometries, with mean heights ranging from 8 to 16 mm were coupled with a reactive oxygen, glucose and lactate transport model linked with tissue deformations and osmosis . In previous studies, reduced formulations of the divergence of the solute flux (∇ .J = ∇ . (D∇ C) = ∇ D. ∇ C +D∇ 2C) ignored the dependence of the diffusion on the deformation gradients, ∇ D. ∇C. We simulated this phenomenon to explore its significance in mechano-metabolic -transport couplings, in the different geometries, over 24h of simulated rest (8h) and physical activity (16h). ∇ D. ∇ C affected the daily variations of glucose concentrations in IVD thinner than 12 mm but with neglectable variation ranges, while not considering ∇ D. ∇ C in taller discs only slightly overestimated the glucose concentration. Most importantly, tall IVD had nearly 60% less glucose than thin IVD, with local drops below the concentration of 0.5 mM, considered to be critical for disc cells3, in the anterior nucleus pulposus. On the one hand, previous reduced formulations for mechanometabolic-transport models of the IVD seem acceptable, even for patient-specific modelling. On the other hand, tall IVD might suffer from unfortunate combinations of deformation-dependent solute diffusion and large diffusion distances, which may favor early. Acknowledgements: Catalan Government and European Commission (2020 BP 00282; ERC-2021-CoG-O-Health-101044828)

Although 3D-printed porous dental implants may possess improved osseointegration potential, they must exhibit appropriate fatigue strength. Finite element analysis (FEA) has the potential to predict the fatigue life of implants and accelerate their development. This work aimed at developing and validating an FEA-based tool to predict the fatigue behavior of porous dental implants. Test samples mimicking dental implants were designed as 4.5 mm-diameter cylinders with a fully porous section around bone level. Three porosity levels (50%, 60% and 70%) and two unit cell types (Schwarz Primitive (SP) and Schwarz W (SW)) were combined to generate six designs that were split between calibration (60SP, 70SP, 60SW, 70SW) and validation (50SP, 50SW) sets. Twenty-eight samples per design were additively manufactured from titanium powder (Ti6Al4V). The samples were tested under bending compression loading (ISO 14801) monotonically (N=4/design) to determine ultimate load (F. ult. ) (Instron 5866) and cyclically at six load levels between 50% and 10% of F. ult. (N=4/design/load level) (DYNA5dent). Failure force results were fitted to F/F. ult. = a(N. f. ). b. (Eq1) with N. f. being the number of cycles to failure, to identify parameters a and b. The endurance limit (F. e. ) was evaluated at N. f. = 5M cycles. Finite element models were built to predict the yield load (F. yield. ) of each design. Combining a linear correlation between FEA-based F. yield. and experimental F. ult. with equation Eq1 enabled FEA-based prediction of F. e. . For all designs, F. e. was comprised between 10% (all four samples surviving) and 15% (at least one failure) of F. ult. The FEA-based tool predicted F. e. values of 11.7% and 12.0% of F. ult. for the validation sets of 50SP and 50SW, respectively. Thus, the developed FEA-based workflow could accurately predict endurance limit for different implant designs and therefore could be used in future to aid the development of novel porous implants. Acknowledgements: This study was funded by EU's Horizon 2020 grant No. 953128 (I-SMarD). We gratefully acknowledge the expert advice of Prof. Philippe Zysset

This study investigates the relationships between Intervertebral Disc (IVD) morphology and biomechanics using patient-specific (PS) finite element (FE) models and poromechanical simulations. 169 3D lumbar IVD shapes from the European project MySpine (FP7-269909), spanning healthy to Pfirrmann grade 4 degeneration, were obtained from MRIs. A Bayesian Coherent Point Drift algorithm aligned meshes to a previously validated structural FE mesh of the IVD. After mesh quality analyses and Hausdorff distance measurements, mechanical simulations were performed: 8 and 16 hours of sleep and daytime, respectively, applying 0.11 and 0.54 MPa of pressure on the upper cartilage endplate (CEP). Simulation results were extracted from the anterior (ATZ) and posterior regions (PTZ) and the center of the nucleus pulposus (CNP). Data mining was performed using Linear Regression, Support Vector Machine, and eXtreme Gradient Boosting techniques. Mechanical variables of interest in DD, such as pore fluid velocity (FLVEL), water content, and swelling pressure, were examined. The morphological variables of the simulated discs were used as input features. Local morphological variables significantly impacted the local mechanical response. The local disc heights, respectively in the mid (mh), anterior (ah), and posterior (ph) regions, were key factors in general. Additionally, fluid transport, reflected by FLVEL, was greatly influenced (r2 0.69) by the shape of the upper and lower cartilage endplates (CEPs). This study suggests that disc morphology affects Mechanical variables of interest in DD. Attention should be paid to the antero-posterior distribution and local effects of disc heights. Surprisingly, the CEP morphology remotely affected the fluid transport in NP volumes around mid-height, and mechanobiological implications shall be explored. In conclusion, patient-specific IVD modeling has strong potential to unravel important correlations between IVD phenotypes and local tissue regulation. Acknowledgments: European Commission: Disc4All-MSCA-2020-ITN-ETN GA: 955735; O-Health-ERC-CoG-2021-101044828

It is known that the gait dynamics of elderly substantially differs from that of young people. However, it has not been well studied how this age-related gait dynamics affects the knee biomechanics, e.g., cartilage mechanical response. In this study, we investigated how aging affects knee biomechanics in a female population using subject-specific computational models. Two female subjects (ages of 23 and 69) with no musculoskeletal disorders were recruited. Korea National Institute for Bioethics Policy Review Board approved the study. Participants walked at a self-selected speed (SWS), 110% of SWS, and 120% of SWS on 10 m flat ground. Three-dimensional marker trajectories and ground reaction forces (Motion Analysis, USA), and lower limbs’ muscle activities were measured (EMG, Noraxon USA). Knee cartilage and menisci geometries were obtained from subjects’ magnetic resonance images (3T, GE Health Care). An EMG-assisted musculoskeletal finite element modeling workflow was used to estimate knee cartilage tissue mechanics in walking trials. Knee cartilage and menisci were modeled using a transversely isotropic poroviscoelastic material model. Walking speed in SWS, 110%, and 120% of SWS were 1.38 m/s, 1.51 m/s, and 1.65 m/s for the young, and 1.21 m/s, 1.34 m/s and 1.46 m/s for the elderly, respectively. The maximum tensile stress in the elderly tibial cartilage was ~25%, ~33%, and ~32% lower than the young at SWS, 110%, and 120% of SWS, respectively. These preliminary results suggest that the cartilage in the elderly may not have enough stimulation even at 20% increases in walking speed, which may be one reason for tissue degeneration. To enhance these findings, further study with more subjects and different genders will investigate how age-related gait dynamics affects knee biomechanics. Acknowledgments: Australian NHMRC Ideas Grant (APP2001734), KITECH (JE220006)

Even minor lesions in articular cartilage (AC) can cause underlying bone damage creating an osteochondral (OC) defect. OC defects can cause pain, impaired mobility and can develop to osteoarthritis (OA). OA is a disease that affects nearly 10% of the population worldwide. [1]. , and represents a significant economic burden to patients and society. [2]. While significant progress has been made in this field, realising an efficacious therapeutic option for unresolved OA remains elusive and is considered one of the greatest challenges in the field of orthopaedic regenerative medicine. [3]. Therefore, there is a societal need to develop new strategies for AC regeneration. In recent years there has been increased interest in the use of tissue-specific aligned porous freeze-dried extracellular matrix (ECM) scaffolds as an off-the-shelf approach for AC repair, as they allow for cell infiltration, provide biological cues to direct target-tissue repair and permit aligned tissue deposition, desired in AC repair. [4]. However, most ECM-scaffolds lack the appropriate mechanical properties to withstand the loads passing through the joint. [5]. One solution to this problem is to reinforce the ECM with a stiffer framework made of synthetic materials, such as polylactic acid (PLA). [6]. Such framework can be 3D printed to produce anatomically accurate implants. [7]. , attractive in personalized medicine. However, typical 3D prints are static, their design is not optimized for soft-hard interfaces (OC interface), and they may not adapt to the cyclic loading passing through our joints, thus risking implant failure. To tackle this limitation, more compliant or dynamic designs can be printed, such as coil-shaped structures. [8]. Thus, in this study we use finite element modelling to create different designs that mimic the mechanical properties of AC and prototype them in PLA, using polyvinyl alcohol as support. The optimal design will be combined with an ECM scaffold containing a tailored microarchitecture mimicking aspects of native AC. Acknowledgments: This project has received funding from the European Union's Horizon Europe research and innovation MSCA PF programme under grant agreement No. 101110000

The December 2023 Spine Roundup³⁶⁰ looks at: Does size matter in adolescent pedicle screws?; Effect of lumbar fusion and pelvic fixation rigidity on hip joint stress: a finite element analysis; Utility of ultrasonography in the diagnosis of lumbar spondylolysis in adolescent patients; Rett syndrome-associated scoliosis a national picture.

Abstract

Aims. Focal knee arthroplasty is an attractive alternative to knee arthroplasty for young patients because it allows preservation of a large amount of bone for potential revisions. However, the mechanical behaviour of cartilage has not yet been investigated because it is challenging to evaluate in vivo contact areas, pressure, and deformations from metal implants. Therefore, this study aimed to determine the contact pressure in the tibiofemoral joint with a focal knee arthroplasty using a finite element model. Methods. The mechanical behaviour of the cartilage surrounding a metal implant was evaluated using finite element analysis. We modelled focal knee arthroplasty with placement flush, 0.5 mm deep, or protruding 0.5 mm with regard to the level of the surrounding cartilage. We compared contact stress and pressure for bone, implant, and cartilage under static loading conditions. Results. Contact stress on medial and lateral femoral and tibial cartilages increased and decreased, respectively, the most and the least in the protruding model compared to the intact model. The deep model exhibited the closest tibiofemoral contact stress to the intact model. In addition, the deep model demonstrated load sharing between the bone and the implant, while the protruding and flush model showed stress shielding. The data revealed that resurfacing with a focal knee arthroplasty does not cause increased contact pressure with deep implantation. However, protruding implantation leads to increased contact pressure, decreased bone stress, and biomechanical disadvantage in an in vivo application. Conclusion. These results show that it is preferable to leave an edge slightly deep rather than flush and protruding. Cite this article: Bone Joint Res 2023;12(8):497–503

The objective of this study was to use patient-specific finite element modeling to measure the 3D interfragmentary strain environment in clinically realistic fractures. The hypothesis was that in the early post-operative period, the tissues in and around the fracture gap can tolerate a state of strain in excess of 10%, the classical limit proposed in the Perren strain theory. Eight patients (6 males, 2 females; ages 22–95 years) with distal femur fractures (OTA/AO 33-A/B/C) treated in a Level I trauma center were retrospectively identified. All were treated with lateral bridge plating. Preoperative computed tomography scans and post-operative X-rays were used to create the reduced fracture models. Patient-specific materials properties and loading conditions (20%, 60%, and 100% body weight (BW)) were applied following our published method.[1]. Elements with von Mises strains >10% are shown in the 100% BW loading condition. For all three loading scenarios, as the bridge span increased, so did the maximum von Mises strain within the strain visualization region. The average gap closing (Perren) strain (mean ± SD) for all patient-specific models at each body weight (20%, 60%, and 100%) was 8.6% ± 3.9%, 25.8% ± 33.9%, and 39.3% ± 33.9%, while the corresponding max von Mises strains were 42.0% ± 29%, 110.7% ± 32.7%, and 168.4% ± 31.9%. Strains in and around the fracture gap stayed in the 2–10% range only for the lowest load application level (20% BW). Moderate loading of 60% BW and above caused gap strains that far exceeded the upper limit of the classical strain rule (<10% strain for bone healing). Since all of the included patients achieved successful unions, these findings suggest that healing of distal femur fractures may be robust to localized strains greater than 10%

Autologous osteochondral grafting has demonstrated positive outcomes for treating articular cartilage defects by replacing the damaged region with a cylindrical graft consisting of bone with a layer of cartilage, taken from a non-loadbearing region of the knee. Despite positive clinical use, factors that cause graft subsidence or poor integration are relatively unknown. The aim of this study was to develop finite element (FE) models of osteochondral grafts within a tibiofemoral joint and to investigate parameters affecting osteochondral graft stability. Initial experimental tests on cadaveric femurs were performed to calibrate the bone properties and graft-bone frictional forces for use in corresponding FE models, generated from µCT scan data. The effects of cartilage defects and osteochondral graft repair were measured by examining contact pressure changes using in vitro tests on a single cadaveric human tibiofemoral joint. Six defects were created in the femoral condyles which were subsequently treated with osteochondral autografts or metal pins. Matching µCT scan-based FE models were created, and the contact patches were compared. Sensitivity to graft bone properties was investigated. The bone material properties and graft-bone frictional forces were successfully calibrated from the initial tests with good resulting levels of agreement (CCC=0.87). The tibiofemoral joint experiment provided a range of cases to model. These cases were well captured experimentally and represented accurately in the FE models. Graft properties relative to host bone had large effects on immediate graft stability despite limited changes to resultant cartilage contact pressure. Model confidence was built through extensive validation and sensitivity testing, and demonstrated that specimen-specific properties were required to accurately represent graft behaviour. The results indicate that graft bone properties affect the immediate stability, which is important for the selection of allografts and design of future synthetic grafts. Acknowledgements. Supported by the EPSRC-EP/P001076

Virtual mechanical testing is a method for measuring bone healing using finite element models built from computed tomography (CT) scans. Previously, we validated a dual-zone material model for ovine fracture callus that differentiates between mineralized woven bone and soft tissue based on radiodensity. 1. The objective of this study was to translate the dual-zone material model from sheep to two important clinical scenarios: human tibial fractures in early-stage healing and late-stage nonunions. CT scans for N = 19 tibial shaft fractures were obtained prospectively at 12 weeks post-op. A second group of N = 33 tibial nonunions with CT scans were retrospectively identified. The modeling techniques were based on our published method. 2. The dual-zone material model was implemented for humans by performing a cutoff sweep for both the 12-week and nonunion groups. Virtual torsional rigidity (VTR) was calculated as VTR = ML/φ [N-m. 2. /°], where M is the moment reaction, L is the diaphyseal segment length, and φ is the angle of twist. As the soft tissue cutoff was increased, the rigidity of the clinical fractures decreased and soft tissue located within the fracture gaps produced higher strains that are not predicted without the dual zone approach. The structural integrity of the nonunions varied, ranging from very low rigidities in atrophic cases to very high rigidities in highly calcified hypertrophic cases, even with dual-zone material modeling. Human fracture calluses are heterogeneous, comprising of woven bone and interstitial soft tissue. Use of a dual-zone callus material model may be instrumental in identifying delayed unions during early healing when callus formation is minimal and/or predominantly fibrous with little mineralization. ACKNOWLEDGEMENTS:. This work was supported by the National Science Foundation (NSF) grant CMMI-1943287

Non-optimal clinical alignment of components in total hip replacements (THRs) may lead to edge loading of the acetabular cup liner. This has the potential to cause changes to the liner rim not accounted for in standard wear models. A greater understanding of the material behaviours could be beneficial to design and surgical guidance for THR devices. The aim of this research was to combine finite element (FE) modelling and experimental simulation with microstructural assessment to examine material behaviour changes during edge loading. A dynamic deformable FE model, matching the experimental conditions, was created to simulate the stress strain environment within liners. Five liners were tested for 4Mc (million cycles) of standard loading (ISO14242:1) followed by 3Mc of edge loading with dynamic separation (ISO14242:4) in a hip simulator. Microstructural measurements by Raman spectroscopy were taken at unloaded and highly loaded rim locations informed by FE results. Gravimetric and geometric measurements were taken every 1Mc cycles. Under edge loading, peak Mises stress and plastic deformation occur below the surface of the rim during heel strike. After 7Mc, microstructural analysis determined edge loaded regions had an increased crystalline mass fraction compared to unloaded regions (p<0.05). Gravimetric wear rates of 12.5mm. 3. /Mc and 22.3mm. 3. /Mc were measured for standard and edge loading respectively. A liner penetration of 0.37mm was measured after 7Mc. Edge loading led to an increase in gravimetric wear rate indicating a different wear mechanism is occurring. FE and Raman results suggest that changes to material behaviour at the rim could be possible. These methods will now be used to assess more liners and over a larger number of cycles. They have potential to explore the impact of edge loading on different surgical and patient variables

Occlusal loading and muscle forces during mastication aids in assessment of dental restorations and implants and jaw implant design; however, three-dimensional bite forces cannot be measured with conventional transducers, which obstruct the native occlusion. The aim of this study was to combine accurate jaw kinematics measurements, together with subject-specific computational modelling, to estimate subject-specific occlusal loading and muscle forces during mastication. Motion experiments were performed on one male participant (age: 39yrs, weight: 82kg) with healthy dentition. Two low-profile magnetic sensors were fixed to the participant's teeth and the two dental arches digitised using an intra-oral scanner. The participant performed ten continuous of chewing on a polyurethane rubber sample of known material properties, followed by maximal compression (clenching). This was repeated at the molars, premolars of both the left and right sides, and central incisors. Jaw motion was simultaneously recorded from the sensors, and finite element modelling used to estimate bite force. Specifically, simulations of chewing and biting were performed by driving the model using the measured kinematics, and bite force magnitude and direction quantified. Muscle forces were then evaluated using a rigid-body musculoskeletal model of the patient's jaw. The first molars generated the largest bite forces during chewing (left: 309 N, right: 311 N) and maximum-force biting (left: 496 N, right: 495 N). The incisors generated the smallest bite forces during chewing (75 N) and maximum-force biting (114 N). The anterior temporalis and superficial masseter muscles had the largest contribution to maximum bite force, followed by the posterior temporalis and medial pterygoid muscles. This study presents a new method for estimating dynamic occlusal loading and muscle forces during mastication. These techniques provide new knowledge of jaw biomechanics, including muscle and occlusal loading, which will be useful in surgical planning and jaw implant design

Orthopaedic impaction-instruments are used to drive implants into the bone of the patient. Pre-clinical experimental testing protocols and computer models of those are used to assess robustness and functional efficiency of such instruments. This generally involves impaction of the instrument mounted on a substrate that should represent the mechanics of the patient. In this study, the effects of the substrate on stressing of the impaction-instruments were investigated using dynamic finite element analysis. Model results were compared with experimental data from lab protocols, which have been derived to recreate the mechanics of cadaveric implantations, which represent clinical conditions. FEA models of selected experimental protocols were created in which a simplified instrument was impacted on substrates with varying material properties and boundary conditions. After impaction, the instrument settled into a modal vibration which then decayed over time. The resulting axial strain data from the computational model was compared to strain-gauge data collected from experimental measurements. Strain signal amplitude, frequency and decay were compared. The damping-ratio was derived from the decay of the strain signal. The computational model slightly over-predicted the initial experimental strain amplitudes in all cases, but the frequency of the cyclic strain signals matched. However, the model underestimated the experimentally measured rate of signal decay. Inclusion of implant seating and soft-tissue conditions had little effect on decay. Clinical failures of impaction-instruments may be related to multiple fatigue cycles for each impaction and should be modelled accurately to allow failure prediction. Any soft substrate results in an impedance mismatch at the instrument interface, which reflects the pressure wave and causes vibration with a frequency related to the speed-of-sound in the instrument, and its geometry. While this could be accurately modelled computationally, signal decay was underestimated. Further experimental quantification of energy losses will be important to understand vibration decay