Abstract. Objectives. The aim of this study was to develop an open-source finite element model of the ankle for identification of the best clinical treatment to restore stability to the ankle after injury. Methods. The ankle geometry was defined from the Visible Human Project Female CT dataset available from the National Library of Medicine, and segmented using Dragonfly software (Object Research Systems, 2020). The finite element model was created with FEBio (University of Utah, 2021) using the dynamic nonlinear implicit solver. Linear isotropic material properties were assigned to the bones (E=7300MPa, ν=0.3, ρ=1730kg/m. 3. ) and cartilage (E=10MPa, ν=0.4, ρ=1100kg/m. 3. ). Spring elements were used to represent the ligaments and material properties were taken from Mondal et al. [1]. Lagrangian contact was defined between the cartilaginous surfaces with μ=0.003. A standing load case was modelled, assuming even distribution of load between the feet. A reaction force of 344.3N was applied to the base of the foot, a muscle force of 252.2N, and the proximal ends of the tibia and fibula were fully constrained. Results. The von Mises stresses closely matched those reported by Mondal et al. for the fibula (Present study: 1.00MPa, Mondal: 1.30MPa) and the talus (Present study: 2.20MPa, Mondal: 2.39MPa). However stresses within the tibia were underpredicted (Present study: 1.08MPa, Mondal: 5.86MPa). This was because the present study modelled a shorter tibial length because of a limitation in the CT slices available, which reduced the bending force. Conclusions. This first step in producing an open source ankle model for the orthopaedics community has shown the potential of the model to generate results comparable with those found in the literature. Future work is underway to examine the robustness of the model under different loading and explore alternative open-source CT datasets. [1] Mondal, S., & Ghosh, R. (2017). J Orthopaedics, 14(3), 329–335. . https://doi.org/10.1016/j.jor.2017.05.003

Finite element models of the musculoskeletal system have the possibility of describing the in vivo situation to a greater extent than a single in vitro experimental study ever could. However these models and the assumptions made must be validated before they can be considered truly useful. The object of this study was to validate, using digital image correlation (DIC) and strain gauging, a novel free boundary condition finite element model of the femur. The femur was treated as a complete musculoskeletal construct without specific fixed restraint acting on the bone. Spring elements with defined force-displacement relationships were used to characterize all muscles and ligaments crossing the hip and knee joints. This model was subjected to a loading condition representing single leg stance. From the developed model muscle, ligament and joint reaction forces were extracted as well as displacement and strain plots. The muscles with the most influence were selected to be represented in the simplified experimental setup. To validate the finite element model a balanced in vitro experimental set up was designed. The femur was loaded proximally through a construct representative of the pelvis and balanced distally on a construct representing the tibio-femoral joint. Muscles were represented using a cabling system with glued attachments. Strains were recorded using DIC and strain gauging. DIC is an image analysis technique that enables non-contact measurement of strains across surfaces. The resulting strain distributions were compared to the finite element model. The finite element model produced hip and knee joint reaction forces comparable to in vivo data from instrumented implants. The experimental models produced strain data from both DIC and strain gauging; these were in good agreement with the finite element models. The DIC process was also shown to be a viable method for measuring strain on the surface of the specimen. In conclusion a novel approach to finite element modeling of the femur was validated, allowing greater confidence for the model to be further developed and used in clinical settings

Worldwide, osteoporosis, causes more than 8.9 million fractures annually, resulting in an osteoporotic fracture every 3 seconds, where 1 in every 3 women and 1 in every 5 men aged over 50 will experience osteoporotic fractures at least once in their lifetime. Vertebral fractures, estimated at 1.4 million/year are among the most common fractures, posing enormous health and socioeconomic challenges to the individual and society at large. Considering that the great majority of individuals at high risk (up to 80%), who have already had at least one osteoporotic fracture, are neither identified nor treated, prediction of the risk factors for vertebral fractures can be of great value for prevention/early diagnosis. Recent studies show that finite element analysis of computed tomography (CT) scans provides noninvasive means to assess fracture risk and has the potential to be clinically implemented upon proper validation. The objective of this study was to develop a voxel-based finite element model using quantitative computed tomography (QCT) images in conjunction with in-vitro experiments to evaluate the strength of the vertebral bodies and predict the fracture risk criteria. A total of 10 vertebrae were dissected from juvenile sheep lumbar spines. The attached soft tissues and posterior elements and facet joints were completely removed, and the upper and lower vertebral bodies were polished using glass paper to provide smooth surfaces. The specimens were wrapped in phosphate buffer saline (PBS) soaked gauze, sealed in plastic bags, and stored in a refrigerator at −22°C. QCT scans of the specimens were captured using a bone density calibration phantom (QRM Co., Moehrendorf, Germany) with three 18 mm cylindrical inserts, providing 0, 100 and 200 mg HA/ccm, respectively. All the specimens, preserved hydrated in PBS solution, were mechanically tested at room temperature using a mechanical testing apparatus (Zwick/Roell, Ulm-Germany). The QCT images were then used to reconstruct the voxel-based FE model employing a custom-developed heterogeneous material mapping code. Five different equations for the correlation of the density and the elastic modulus were used to validate the efficiency of the FE model as compared to the in-vitro experiments. The results of the voxel-based FE models matched well with the in-vitro experiments, with an average error of 11.38 (±4.09)% based on the power law equation. A failure criterion was embedded in the FE models and the initiation of fracture was successfully predicted for all specimens. Further, typical kyphoplasty treatment was simulated in the 5 models to evaluate the application of the validated algorithm in the estimation of the failure patterns. Our novel voxel-based FE model can be used in future studies to predict the outcome of different types of therapeutic modalities/surgeries and estimate fracture risk including postoperative fractures

Summary. A retrospective study on 98 patients shows that FE-based bone strength from CT data (using validated FE models) is a suitable candidate to discriminate fractured versus controls within a clinical cohort. Introduction. Subject-specific Finite element models (FEM) from CT data are a promising tool to non-invasively assess the bone strength and the risk of fracture of bones in vivo in individual patients. The current clinical indicators, based on the epidemiological models like the FRAX tool, give limitation estimation of the risk of femoral neck fracture and they do not account for the mechanical determinants of the fracture. Aim of the present study is to prove the better predictive accuracy of individualised computer models based a CT-FEM protocol, with the accuracy of a widely used standard of care, the FRAX risk indicator. Patients and Methods. This retrospective cohort is individually-matched case control study composed by 98 Caucasian women who were at least 5 years post menopause. The case group consisted of 49 patients who had sustained a hip fracture (36 intra-capsular and 13 extra-capsular fractures) within the previous 90 days due to low-energy trauma. The CT datasets were segmented (using the ITK-Snap software) in order to extract the periosteal bone surface. Unstructured meshes (10-node tetrahedral elements) were generated using ANSYS mesh morphing software. Each CT dataset was calibrated using the European Spine Phantom. The inhomogeneous material properties were mapped from CT datasets into the FEM with the BoneMat_V3 software. Bone strength was evaluated in quasi-axial loading conditions, for a set of 12 different configurations sampling the cone of recorded in vivo hip joint reactions, and was defined as the minimum load inducing on the femoral neck surface an elastic principal strain value greater than a limit value. Results. There were no statistically significant difference between the fracture and the control groups for age, height and weight (p<0.05). All indices of areal bone mineral density (aBMD) and the volumetric mineral density (vBMD) between fractured and controls showed on average a lower value for fractured respect of the controls, with similar mean difference (14% for aBMD and 13% for the vBMD). FEM-predicted strength differed between fractured and non-fractured on average for 20%. To evaluate its ability to identify patients at risk of hip fracture, FEM-based strength was compared to the FRAX predictor by computing for each predictor the Receiver Operating Characteristic (ROC) curve, and the Area Under the Curve (AUC). The individualised risk predictor based on FEM bone strength was found to perform significantly better (AUC=0.76) than FRAX (AUC=0.66). When the FEM-based strength indicator was combined with available clinical information in a logistic regression, the resulting predictor achieved in this retrospective study an excellent accuracy (AUC=0.82). Discussion. This study confirms that individualised, CT- FEM, when generated using to the state-of-the-art protocols, can provide a predictor of the risk of hip fracture more accurate than those based on clinical data alone. In the integrated workflow developed in the VPHOP Project (FP7-ICT-223865) CT-based risk prediction is requested only for those patients for whom the clinical decision is uncertain

Iterative finite element (FE) models are used to simulate bone remodelling that takes place due to the surgical insertion of an implant or to simulate fracture healing. In such simulations element material properties are calculated after each iteration of solving the model. New material properties are calculated based on the results derived by the model during the last iteration. Once the FE model has gone through a number of such iterations it is often necessary to assess the remodelling that has taken place. The method widely used to do this is to analyse element Young's modulus plots taken at particular sections through the model. Although this method gives relevant information which is often helpful when comparing different implants, the information is rather abstract and is difficult to compare with patient data which is commonly in the form of radiographs. The authors suggest a simple technique that can be used to generate synthetic radiograph images from FE models. These images allow relatively easy comparisons of FE derived information with patient radiographs. Another clear advantage of this technique is that clinicians (who are familiar with reading radiographs) are able to understand and interpret them readily. To demonstrate the technique a three dimensional (3D) model of the proximal tibia implanted with an Oxford Unicompartmental Knee replacement was created based on CT data obtained from a cadaveric tibia. The model's initial element material properties were calculated from the same CT data set using a relationship between radiographic density and Young's modulus. The model was subject to simplified loading conditions and solved over 365 iterations representing one year of in vivo remodelling. After each iteration the element material properties were recalculated based on previously published remodelling rules. Next, synthetic anteroposterior radiographs were generated by back calculating radiographic densities from material properties of the model after 365 iterations. A 3D rectangular grid of sampling points which encapsulated the model was defined. For each of the elements in the FE model radiographic densities were back calculated based on the same relationships used to calculate material properties from radiographic densities. The radiographic density of each element was assigned to all the sampling grid points within the element. The 3D array of radiographic densities was summed in the anteroposterior direction thereby creating a 2D array of radiographic densities. This 2D array was plotted giving an image analogous to anteroposterior patient radiographs. Similar to a patient radiograph denser material appeared lighter while less dense material appeared darker. The resulting synthetic radiographs were compared to patient radiographs and found to have similar patterns of dark and light regions. The synthetic radiographs were relatively easy to produce based on the FE model results, represented FE results in a manner easily comparable to patient radiographs, and represented FE results in a clinician friendly manner

Clinical investigations show that the cervical spine presents wide inter-individual variability, where its motion patterns and load sharing strongly depend on the anatomy. The magnitude and scope of cervical diseases, including disc degeneration, stenosis, and spondylolisthesis, constitute serious health and socioeconomic challenges that continue to increase along with the world”s growing aging population. Although complex exact finite element (FE) modeling is feasible and reliable for biomechanical studies, its clinical application has been limited as it is time-consuming and constrained to the input geometry, typically based on one or few subjects. The objective of this study was twofold: first to develop a validated parametric subject-specific FE model that automatically updates the geometry of the lower cervical spine based on different individuals; and second to investigate the motion patterns and biomechanics associated with typical cervical spine diseases. Six healthy volunteers participated in this study upon informed consent. 26 parameters were identified and measured for each vertebra in the lower cervical spine from Lateral and AP radiographs in neutral, flexion and extension viewpoints in the standing position. The lower cervical FE model was developed including the typical vertebrae (C3-C7), intervertebral discs, facet joints, and ligaments using ANSYS (PA, USA). In order to validate the FE model, the bottom surface of C7 was fixed, and a 73.6N preload together with a 1.8 N.m pure moment were input into the model in both flexion and extension. The results were compared to experimental studies from literature. Disc degeneration disease (DDD) was used as an example, where the geometry of C5-C6 disc was changed in the model to simulate 3 different grades of disc degeneration (mimicking grades 1 to 3), and the resulting biomechanical responses were evaluated. The average ranges of motion (ROM) were found to be 4.84 (±0.73) degrees and 5.36 (±0.68) degrees for flexion and extension for C5-C6 functional unit, respectively, in alignment with literature. The total ROM of the model with disc generation grades 2 and 3 was found to have decreased significantly as compared to the intact model. In contrast, the axial stresses on the degenerated discs were significantly higher than the intact discs for all 3 degeneration grades. Our preliminary results show that this novel validated subject-specific FE model provides a potential valuable tool for noninvasive time and cost effective analyses of cervical spine biomechanical (kinematic and kinetic) changes associated with various diseases. The model also provides an opportunity for clinicians to use quantitative data towards subject-specific informed therapy and surgical planning. Ongoing and future work includes expanding the studied population to investigate individuals with different cervical spine afflictions

Summary. At the clinical CT image resolution level, there is no influence of the image voxel size on the derived finite element human cancellous bone models. Introduction. Computed tomography (CT)-based finite element (FE) models have been proved to provide a better prediction of vertebral strength than dual-energy x-ray absorptiometry [1]. FE models based on µCTs are able to provide the golden standard results [2], but due to the sample size restriction of the µCT and the XtremeCT machines, the clinical CT-based FE models is still the most promising tool for the in vivo prediction of vertebrae's strength. It has been found [3] that FE predicted Young's modulus of human cancellous bone increases as the image voxel size increases at the µCT resolution level [3]. However, it is still not clear whether the image voxel size in the clinical range has an impact on the predicted mechanical behavior of cancellous bone. This study is designed to answer this question. Methods. For this study, 6 thoracolumbar vertebrae (Th12) obtained from the female donors were scanned in the non-dissected cadavers under 2 different resolutions – group A: 120 kVp, 100 mAs, with a resolution of 0.29×0.29×1.3 mm. 3. ; group B: 120 kVp, 360 mAs, with a resolution of 0.18×0.18×0.6 mm. 3. A solid calibration phantom (QRM-BDC) was placed beneath the cadavers during the scans. Cuboids with the size of 12.3×12.3×14.3 mm. 3. were cropped from the center of each vertebral body. The FE model was created by converting each image voxel into hexahedron (C3D8). Inhomogeneous material property was defined for the cuboid [4], i.e. the image greyscale value were firstly calibrated into the bone mineral density (BMD), then the Young's modulus and yield stress were calculated from the BMD [5] for each element. Statistical analysis was performed to compare the FE predicted mechanical properties between the groups and the significance level was set to 95% (α=0.05). Results. The trabecular structure is more clearly mimicked in the models from group B than those from group A. The modulus (mean ± SE) in group A is 5.9% higher than that in group B (193.33 ± 31.67 MPa vs. 182.50 ± 27.07 MPa). The yield strength (mean ± SE) in group A is 6.4% higher than that in group B (0.99 ± 0.21MPa vs. 0.93 ± 0.17MPa). However, the paired t-test shows there is no significant difference of the mechanical properties in the two groups (p=0.109 for the modulus and p=0.234 for the yield strength). Discussion. This study shows that there is no influence of the voxel size on the clinical CT derived FE cancellous bone models. This finding can help choose a better, less invasive CT protocol for the patient when creating a clinical CT image based FE model. Acknowledgements. This study is financially supported by the Federal Ministry of Education and Research and the state of Hamburg, Germany

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)

INTRODUCTION. Growth-guidance constructs are an alternative to growing rods for the surgical treatment of early onset scoliosis (EOS). In growth-guidance systems, free-sliding anchors preserve longitudinal spinal growth, thereby eliminating the need for surgical lengthening procedures. Non-segmental constructs containing ultra-high molecular weight polyethylene (UHMWPE) sublaminar wires have been proposed as an improvement to the traditional Luque trolley. In such a construct, UHMWPE sublaminar wires, secured by means of a knot, serve as sliding anchors at the proximal and distal ends of a construct, while pedicle screws at the apex prevent rod migration and enable curve derotation. Ideally, a construct with the optimal UHMWPE sublaminar wire density, offering the best balance between providing adequate spinal fixation and minimizing surgical exposure, is designed preoperatively for each individual patient. In a previous study, we developed a parametric finite element (FE) model that potentially enables preoperative patient-specific planning of this type of spinal surgery. The objective of this study is to investigate if this model can capture the decrease in range of motion (ROM) after spinal fixation as measured in an experimental study. MATERIALS AND METHODS. In a previous in vitro study, the ROM of an 8-segment porcine spine was measured before and after instrumentation, using different instrumentation constructs with a sequentally decreasing number of wire fixation points. In the current study, the parametric FE model of the thoracolumbar spine was first validated relative to ROM values reported in the literature. The rods, screws, and sublaminar wires were implemented, and the model was subsequently used to replicate the in vitro tests. The experimental and simulated ROM”s for the different instrumentation conditions were compared. RESULTS. Good agreement between in vitro biomechanical tests and FE simulations was observed in terms of the decrease in ROM for the complete construct with wires at each level. The stepwise increase in total ROM with decreasing number of wires at the construct ends was less prominent in silico in comparison to in vitro. CONCLUSION. Important first steps in the implementation and validation of a growth-guidance construct for EOS patients in a patient-specific FE model of the spine have been made in this study. The parametric nature of the FE model allows for rapid personalization. Although further improvements to the model will be necessary to better distinguish between different spinal instrumentation constructs, we conclude that the model can well capture essential aspects of spinal motion and the overall effect of instrumentation

Summary Statement. A coupled finite element - analytical model is presented to predict and to elucidate a clinical healing scenario where bone regenerates in a critical-sized femoral defect, bounded by periosteum or a periosteum substitute implant and stabilised via an intramedullary nail. Introduction. Bone regeneration and maintenance processes are intrinsically linked to mechanical environment. However, the cellular and subcellular mechanisms of mechanically-modulated bone (re-) generation are not fully understood. Recent studies with periosteum osteoprogenitor cells exhibit their mechanosensitivity in vitro and in situ. In addtion, while a variety of growth factors are implicated in bone healing processes, bone morphogenetic protein-2 (BMP-2) is recognised to be involved in all stages of bone regeneration. Furthermore, periosteal injuries heal predominantly via endochondral ossification mechanisms. With this background in mind, the current study aims to understand the role of mechanical environment on BMP-2 production and periosteally-mediated bone regeneration. The one-stage bone transport model [1] provides a clinically relevant experimental platform on which to model the mechanobiological process of periosteum-mediated bone regeneration in a critical-sized defect. Here we develop a model framework to study the cellular-, extracellular- and mechanically-modulated process of defect infilling, governed by the mechanically-modulated production of BMP-2 by osteoprogenitor cells located in the periosteum. Methods. Material properties of the healing callus and periosteum contribute to the strain stimulus sensed by osteoprogenitor cells therein. Using a mechanical finite element model, periosteal surface strains are first predicted as a function of callus properties. Strains are then input to a mechanistic mathematical model, where mechanical regulation of BMP-2 production mediates rates of cellular proliferation, differentiation and extracellular matrix (ECM) production, to predict healing outcomes. A parametric approach enables the spatial and temporal prediction of tissue regeneration via endochondral ossification. Predictions are compared with experimental, micro-computed tomographic and histologic, measures of cartilage and mineralised bone tissue regenerates. Model Predictions in Light of Experimental Case Studies: A validated baseline model predicts defect healing via cellular egression, extracellular matrix production and endochondral ossification, using parameters optimised to mimic experimental outcome measures at initial and final stages of healing. To elucidate which predictive model paramenters result in the intrinsic differences in experimental outcomes between defects bounded by either periosteum in situ or a periosteum substitute implant, model parameters are then varied by orders of magnitude to determine which factors exert dominant influence on achievement of experimentally relevant ECM area outcomes. Considering the complete set of parameters relevant to healing, the rate of osteoprogenitor to osteoblast differentiation, as well as rates of chondrocyte and osteoblast proliferation must be reduced and ECM production by chondrocytes must be increased from baseline, to achieve healing outcomes analogous to those observed in experiments. Discussion/Conclusion. The novel model framework presented here integrates a mechanistic feedback system, based on the mechanosensitivity of periosteal osteoprogenitor cells, which allows for modeling and prediction of tissue regeneration on multiple length and time scales

Abstract

Abstract. OBJECTIVES. Bone health deterioration is a major public health issue. General guidelines for the limitation of bone loss prescribe a healthy lifestyle and a minimum level of physical activity. However, there is no specific recommendation regarding targeted activities that can effectively maintain lumbar spine bone health. To provide a better understanding of such influencing activities, a new predictive modelling framework was developed to study bone remodelling under various loading conditions. METHODS. The approach is based on a full-body subject-specific musculoskeletal model [1] combined with structural finite element models of the lumbar vertebrae. Using activities recorded with the subject, musculoskeletal simulations provide physiological loading conditions to the finite element models which simulate bone remodelling using a strain-driven optimisation algorithm [2]. With a combination of daily living activities representative of a healthy lifestyle including locomotion activities (walking, stair ascent and descent, sitting down and standing up) and spine-focused activities involving twisting and reaching, this modelling framework generates a healthy bone architecture in the lumbar vertebrae. The influence of spine-focused tasks was studied by adapting healthy vertebrae to an altered loading scenario where only locomotion activities were performed. RESULTS. The spine-focused activities were responsible for 57% of the overall bone mechanical stimulus of the five lumbar vertebrae. Cortical bone maintenance was more influenced by these activities in the superior vertebrae than in the inferior ones, with a stimulus degradation of 74% in L1 against 24% in L5 when adapted to the altered loading scenario. Trabecular bone stimulus degradation varied between 53% and 68%. CONCLUSION. The study suggests that locomotion activities are insufficient to maintain lumbar spine bone health. When appropriate, larger spine movements should be recommended as part of the minimum daily physical activities. Declaration of Interest. (b) declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported:I declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research project

Introduction. The ankle cartilage has an important function in walking movements, mainly in sports; for active young people, between 20 and 30 years old, the incidence of osteochondral lesions is more frequent. They are also more frequent in men, affecting around 21,000 patients per year in USA with 6.5% of ankle injuries generating osteochondral lesions. The lesion is a result of ankle sprain and is most frequently found in the medial location, in 53% of cases. The main objective of this work was to develop an experimental and finite element models to study the effect of the ankle osteochondral lesion on the cartilage behavior. Materials and Methods. The right ankle joint was reconstructed from an axial CT scan presenting an osteochondral lesion in the medial position with 8mm diameter in size. An experimental model was developed, to analyze the strains and influence of lesion size and location similar to the patient. The experimental model includes two cartilages constructed by Polyjet™ 3D printing from rubber material (young modulus similar to cartilage) and bone structures from a rigid polymer. The cartilage was instrumented with two rosettes in the medial and lateral regions, near the osteochondral region. The fluid considered was water at room temperature and the experimental test was run at 1mm/s. The Finite element model (FE) includes all the components considered in the experimental apparatus and was assigned the material properties of bone as isotropic and linear elastic materials; and the cartilage the same properties of rubber material. The fluid was simulated as hyper-elastic one with a Mooney-Rivlin behavior, with constants c1=0.07506 and c2=0.00834MPa. The load applied was 680N in three positions, 15º extension, neutral and 10º flexion. Results. The experimental strain measured in the cartilage in the rosettes presents similar behavior in all experiments and repetitions. The maximum value observed near the osteochondral lesion was 3014(±5.6)µε in comparison with the intact condition it was 468 (±1.95)µε. The osteochondral lesion increases the strains around 6.5 times and the synovial liquid reduces the intensity of strain distribution. The numerical model presents a good correlation with the experiments (R2 0.944), but the FE model underestimates the values. Discussion and conclusion. As a first conclusion, the size of the osteochondral lesion is important for the strains developed in cartilage. The size of lesion greater than 10mm is critical for the strains concentration. The synovial fluid present an important aspect in the strains measured, it reduces the strains in the external surface of cartilage and induces an increase in the lower part. This phenomenon should be addressed in more studies to evaluate this effect

Smith's fractures generally occur when falling on a flexed wrist; however, orthopedic trauma surgeons often encounter distal radius fractures with volar displacement in patients who have allegedly fallen on the palm of their hands. This study aimed to reveal both the basic and clinical pathogenesis of Smith's fracture through a step-by-step investigation. We enrolled 17 patients with Smith's fractures, of which 71% fell on the palm and only 6% on the dorsum of the hand. First, we interviewed the outpatients to determine the mechanics of the injury and the position of their arm during injury. Second, we created a three-dimensional (3D) finite element model to predict the arm's position when the Smith's fracture occurred, which finite element analysis revealed as a 30° angle between the long axis of the forearm and the ground in the sagittal plane. Third, using this predicted position, we conducted experiments on 10 fresh frozen cadavers to prove the possibility of causing a Smith's fracture by falling on the palm of the hand. The results showed Smith-type fractures in seven of 10 wrists, whereas Colles-type fractures did not occur. Finally, we analyzed stress distribution in the distal radius when a Smith's fracture occurs using the 3D finite element model. In conclusion, this study demonstrates that Smith's fractures can also occur by falling on the palm of the hand

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)

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

Variations in component positioning of total hip replacements can lead to edge loading of the liner, and potentially affect device longevity. These effects are evaluated using ISO 14242:4 edge loading test results in a dynamic system. Mediolateral translation of one of the components during testing is caused by a compressed spring, and therefore the kinematics will depend on the spring stiffness and damping coefficient, and the mass of the translating component and fixture. This study aims to describe the sensitivity of the liner plastic strain to these variables, to better understand how tests using different simulator designs might produce different amounts of liner rim deformation. A dynamic explicit deformable finite element model with 36mm Pinnacle metal-on-polyethylene bearing geometry (DePuy Synthes, Leeds, UK) was used with material properties for conventional UHMWPE. Setup was 65° clinical inclination, 4mm mismatch, 70N swing phase load, and 100N/mm spring. Fixture mass was varied from 0.5-5kg, spring damping coefficient was varied from 0-2Ns/mm. They were changed independently, and in combination. Maximum separation values were relatively insensitive to changes in the mass, damping coefficient, or both. The sensitivity of peak plastic strain, to this range of inputs, was similar to changing the swing phase load from 70N to approximately 150N – 200N. Increasing the fixture mass and/or damping coefficient increased the peak plastic strain, with values from 0.15-0.19. Liner plastic deformation was sensitive to the spring damping and fixture mass, which may explain some of the differences in fatigue and deformation results in UHMWPE liners tested on different machines or with modified fixtures. These values should be described when reporting the results of ISO14242:4 testing. Acknowledgements. Funded by EPSRC grant EP/N02480X/1; CAD supplied by DePuy Synthes

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

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%

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)