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

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

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

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

Objectives. In this prospective cohort study, we investigated whether patient-specific finite element (FE) models can identify patients at risk of a pathological femoral fracture resulting from metastatic bone disease, and compared these FE predictions with clinical assessments by experienced clinicians. Methods. A total of 39 patients with non-fractured femoral metastatic lesions who were irradiated for pain were included from three radiotherapy institutes. During follow-up, nine pathological fractures occurred in seven patients. Quantitative CT-based FE models were generated for all patients. Femoral failure load was calculated and compared between the fractured and non-fractured femurs. Due to inter-scanner differences, patients were analyzed separately for the three institutes. In addition, the FE-based predictions were compared with fracture risk assessments by experienced clinicians. Results. In institute 1, median failure load was significantly lower for patients who sustained a fracture than for patients with no fractures. In institutes 2 and 3, the number of patients with a fracture was too low to make a clear distinction. Fracture locations were well predicted by the FE model when compared with post-fracture radiographs. The FE model was more accurate in identifying patients with a high fracture risk compared with experienced clinicians, with a sensitivity of 89% versus 0% to 33% for clinical assessments. Specificity was 79% for the FE models versus 84% to 95% for clinical assessments. Conclusion. FE models can be a valuable tool to improve clinical fracture risk predictions in metastatic bone disease. Future work in a larger patient population should confirm the higher predictive power of FE models compared with current clinical guidelines. Cite this article: F. Eggermont, L. C. Derikx, N. Verdonschot, I. C. M. van der Geest, M. A. A. de Jong, A. Snyers, Y. M. van der Linden, E. Tanck. Can patient-specific finite element models better predict fractures in metastatic bone disease than experienced clinicians? Towards computational modelling in daily clinical practice. Bone Joint Res 2018;7:430–439. DOI: 10.1302/2046-3758.76.BJR-2017-0325.R2

Introduction. Density-modulus relationships are often used to map the mechanical properties of bone based on CT- intensity in finite element models (FEMs). Although these relationships are thought to be site-specific, relationships developed for alternative anatomic locations are often used regardless of bone being modeled. Six relationships are commonly used in finite element studies of the shoulder; however, the accuracy of these relationships have yet to be compared. This study compares each of these six relationships ability to predict apparent strain energy density (SED. app. ) in trabecular bone cores from the glenoid. Methods. Quantitative-CT (QCT) (0.625 mm isotropic voxels), and µ-CT scans (0.032 mm isotropic voxels) were obtained for fourteen cadaveric scapulae (7 male, 7 female). Micro finite element models (µ-FEMs) were created from 98 virtual ‘cores’ using direct conversion to hexahedral elements. Two µ-FEM cases were considered: homogeneous tissue modulus of 20 GPa, and heterogeneous tissue modulus scaled by CT intensity of the µ-CT images (196 models). Each µ- FEM model was compressively loaded to 0.5% apparent strain and apparent strain energy density (SED. app. ) was calculated. Additionally, each of the six density-modulus relationships were used to map heterogeneous material properties to co- registered QCT-derived models (588 models in total). The loading and boundary conditions were replicated in the QCT-FEMs and the SED. app. was calculated and compared to the µ-FEM SED. app. To account for more samples than donors, restricted maximum likelihood estimation (REML) linear regression compared µ-FEM SED. app. and QCT-FEM SED. app. for each relationship. Results. When considering comparisons between QCT-FEMs and µ-FEMs with a homogeneous tissue modulus, near absolute statistical agreement (Y=X) was observed between the µ-FEMs and the QCT-FEMs using the Morgan et al. (2003) pooled relationship. Not surprisingly, due to the similarity between the two relationships, the Gupta & Dan (2004) and Carter and Hayes (1977) models showed near identical REML linear regression fit parameters. All relationships other than the Morgan et al. (2003) pooled relationship, greatly underestimated the µ-FEM apparent strain energy density (SED. app. ) when considering a homogeneous tissue modulus in the µ-FEMs. The same result with the pooled relationship did not hold true when heterogeneous tissue modulus was considered in the µ-FEMs. The Büchler et al., (2002) relationship most accurately predicted the SED. app. for this comparison. Interestingly, the Gupta & Dan (2004) and Carter and Hayes (1977) relationships again showed near identical REML linear regression fit parameters. DISCUSSION. This study compared six common density-modulus relationships used to map mechanical properties of bone in shoulder FE studies. It was found that when considering a homogeneous tissue modulus for µ-FEMs, relationships pooled from alternative anatomic locations may accurately predict the mechanical properties of glenoid trabecular bone. However, when considering a heterogeneous tissue modulus, this did not hold true. Further studies to determine if these relationships can be translated to whole bones may provide insight into the predictive capabilities of using pooled density-modulus equations in the mapping of mechanical properties in future FEMs of the shoulder

Introduction. Achieving an appropriate primary stability after implantation is a prerequisite for the long-term viability of a dental implant. Virtual testing of the bone-implant construct can be performed with finite element (FE) simulation to predict primary stability prior to implantation. In order to be translated to clinical practice, such FE modeling must be based on clinically available imaging methods. The aim of this study was to validate an FE model of dental implant primary stability using cone beam computed tomography (CBCT) with ex vivo mechanical testing. Method. Three cadaveric mandibles (male donors, 87-97 years old) were scanned by CBCT. Twenty-three bone samples were extracted from the bones and conventional dental implants (Ø4.0mm, 9.5mm length) were inserted in each. The implanted specimens were tested under quasi-static bending-compression load (cf. ISO 14801). Sample-specific homogenized FE (hFE) models were created from the CBCT images and meshed with hexahedral elements. A non-linear constitutive model with element-wise density-based material properties was used to simulate bone and the implant was considered rigid. The experimental loading conditions were replicated in the FE model and the ultimate force was evaluated. Result. The experimental ultimate force ranged between 67 N and 789 N. The simulated ultimate force correlated better with the experimental ultimate force (R. 2. =0.71) than the peri-implant bone density (R. 2. =0.30). Conclusion. The developed hFE model was demonstrated to provide stronger prediction of primary stability than peri-implant bone density. Therefore, hFE Simulations based on this clinically available low-radiation imaging modality, is a promising technology that could be used in future as a surgery planning tool to assist the clinician in evaluating the load-bearing capacity of an implantation site. Acknowledgements. Funding: EU's Horizon 2020 grant No: 953128 (I-SMarD). Dental implants: THOMMEN Medical AG

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

In severe cases of total knee arthroplasty which cannot be treated with off-the-shelf implants anymore custom-made knee implants may serve as one of the few remaining options to restore joint function or to prevent limb amputation. Custom-made implants are specifically designed and manufactured for one individual patient in a single-unit production, in which the surgeon is responsible for the implant design characteristics in consultation with the corresponding engineer. The mechanical performance of these custom-made implants is challenging to evaluate due to the unique design characteristics and the limited time until which the implant is needed. Nevertheless, the custom-made implant must comply with clinical and regulatory requirements. The design of custom-made implants is often based on a underlying reference implant with available biomechanical test results and well-known clinical performance. To support surgeons and engineers in their decision whether a specific implant design is suitable, a method is proposed to evaluate its mechanical performance. The method uses finite element analysis (FEA) and comprises six steps: (1) Identification of the main potential failure mechanism and its corresponding FEA quantity of interest. (2) Reproduction of the biomechanical test of the reference implant via FEA. (3) Identification of the maximum value of the corresponding FEA quantity of interest at the required load level. (4) Definition of this value as the acceptance criteria for the FEA of the custom-made implant. (5) Reproduction of the biomechanical test with the custom-made implant via FEA. (6) Conclusion whether the acceptance criteria is fulfilled or not. The method was applied to two exemplary cases of custom-made knee implants. The FEA acceptance criteria derived from the reference implants were fulfilled in both custom-made implants. Subsequent biomechanical tests verified the FEA results. This study suggests and applies a non-destructive and efficient method for pre-clinical testing of a single-unit custom-made knee implant to evaluate whether the design is mechanically suitable

Failure of reverse total shoulder arthroplasty (rTSA) due to loosening of the metaglene remains a concern. The metaglene is typically affixed to the glenoid via four peripheral bone screws, and the orientations of these screws can affect the stability of the metaglene. The purpose of this finite element analysis (FEA) study was to investigate whether screw orientations should be considered on a patient-specific basis to maximise early fixation. Three-dimensional geometries of four scapula specimens were obtained by segmenting from CT data in 3D Slicer. A metaglene and four rigidly attached 4.5 mm diameter, 18 mm long cylinders representing screws, were placed on each reamed glenoid. Each screw was placed at one of four orientations, 15° or 7.5° toward or away from the central axis of the metaglene face, while all others were held in the baseline (BL) configuration, where all screws were perpendicular to the metaglene face. Finite element models were created by meshing with linear tetrahedral elements. Material properties of titanium (E=113.8 GPa, v=0.34) were applied to the metaglene and screws. Cortical bone material properties were considered uniform (E=17.5 GPa, ν=0.3) while cancellous bone material properties were non-uniform and mapped on an element-by-element basis using CT attenuation data. The scapula was fully constrained, and a 252 N superiorly oriented shear force was applied to the inferior portion of the metaglene. Contact was modelled at bone-implant and bone-screw interfaces. Displacements of the metaglene with respect to the glenoid were measured. The orientations of each screw that minimised in-plane displacement were used for specimen-specific (SS) configurations. A global (GL) configuration was also defined based on the averages of SS orientations. FE model-predicted metaglene displacements of the SS, GL, and BL screw configurations were compared using paired t-tests. The average in-plane metaglene displacements for the SS, GL, and BL configurations were 4.8 ± 1.2, 6.5 ± 3.7, and 5.3 ± 1.5 um, respectively. SS configurations significantly decreased displacements by −0.4 ± 0.3 um (−8.5%, p = 0.024) when compared to BL, but the difference of −1.6 ± 3.1 um (25.3%, p = 0.187) was not significant when compared to the GL configuration. In general, the SS configurations resulted in smaller metaglene displacements than the GL configurations, however the difference was not statistically significant. In one specimen, the GL configuration resulted in abnormally large displacements. These results indicate that, while on average, patient-specific orientations won't yield significantly greater fixation than global configurations; non-patient-specific configurations can, in some cases, yield poor results. Therefore, to ensure optimal fixation for all patients, screw orientations should be considered on a patient-specific basis

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

Osteosynthesis aims to maintain fracture reduction until bone healing occurs, which is not achieved in case of mechanical fixation failure. One form of failure is plastic plate bending due to overloading, occurring in up to 17% of midshaft fracture cases and often necessitating reoperation. This study aimed to replicate in-vivo conditions in a cadaveric experiment and to validate a finite element (FE) simulation to predict plastic plate bending. Six cadaveric bones were used to replicate an established ovine tibial osteotomy model with locking plates in-vitro with two implant materials (titanium, steel) and three fracture gap sizes (30, 60, 80 mm). The constructs were tested monotonically until plastic plate deformation under axial compression. Specimen-specific FE models were created from CT images. Implant material properties were determined using uniaxial tensile testing of dog bone shaped samples. The experimental tests were replicated in the simulations. Stiffness, yield, and maximum loads were compared between the experiment and FE models. Implant material properties (Young's modulus and yield stress) for steel and titanium were 184 GPa and 875 MPa, and 105 GPa and 761 MPa, respectively. Yield and maximum loads of constructs ranged between 469–491 N and 652–683 N, and 759–995 N and 1252–1600 N for steel and titanium fixations, respectively. FE models accurately and quantitatively correctly predicted experimental results for stiffness (R2=0.96), yield (R2=0.97), and ultimate load (R2=0.97). FE simulations accurately predicted plastic plate bending in osteosynthesis constructs. Construct behavior was predominantly driven by the implant itself, highlighting the importance of modelling correct material properties of metal. The validated FE models could predict subject-specific load bearing capacity of osteosyntheses in vivo in preclinical or clinical studies. Acknowledgements: This study was supported by the AO Foundation via the AOTRAUMA Network (Grant No.: AR2021_03)

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: Contemporary surgical interventions for adolescent idiopathic scoliosis (AIS) include both anterior and posterior rod systems, in which a single or double rod construct provides curve correction and stability. This paper presents a methodology for development of patient-specific finite element methods to predict the biomechanical outcomes of scoliosis surgery pre-operatively, with the aim of optimising the performance of instrumentation constructs for anterior single rod AIS surgery. Methods: Geometry for each patient-specific finite element model is obtained from pre-operative thoracolumbar CT scans taken in the supine position using a low dose multi-slice imaging protocol. The finite element model incorporates vertebrae, intervertebral discs, and posterior processes with associated ligaments and zygapophysial joints. A custom pre-processor generates the entire model according to user-specified meshing parameters, providing rapid model generation once the geometric parameters have been extracted from each CT dataset. Material properties are currently based on published values. Simulated movements about axes corresponding to flexion/extension, left/right lateral bending, and trunk rotation are solved using the ABAQUS/Standard software, allowing assessment of predicted loads and stresses before and after addition of instrumentation. Results: The total time per patient required for model generation is currently about six hours, with manual measurement of spine geometry from the CT stack accounting for most of this time. Actual solution time for each finite element model is expected to be around four hours, making patient-specific pre-operative planning for endoscopic scoliosis surgery a feasible option at least in terms of processing time per patient. Discussion: A finite element methodology has been developed for patient-specific simulation of endoscopic scoliosis surgery. Issues to be addressed in future include prescription of patient-specific material properties, analysis of errors associated with geometry measurement from CT scans, and validation of the methodology by comparison of predicted and actual outcomes for scoliosis patients. Patient-specific simulation of scoliosis surgery has the potential to optimize surgical outcomes and reduce biomechanical complications associated with the use of endoscopic scoliosis instrumentation systems

Musculoskeletal models of the lower limb lend insight into muscle forces and joint mechanics during dynamic activities. However, traditional musculoskeletal modeling is based on rigid body assumptions, and frequently represents the knee as a hinge joint, neglecting the complex interactions between the patella, femur, and tibia. Implementation of the musculoskeletal modeling framework in an explicit finite element environment allows joint contact to be easily incorporated, as well as representation of any structure as rigid or fully deformable in order to evaluate, for example, implant stresses or bone strain. Prediction of these values is particularly valuable when evaluating implant mechanics after total knee replacement. A finite element, musculoskeletal model of an implanted right lower limb was constructed, including thirteen muscles crossing the knee joint. A Hill-type muscle model was developed to allow muscle activation within the explicit FE framework. Muscle forces were predicted by optimization of muscle activation patterns during flexion-extension and chair-rise activities. The effect of muscle path representation was investigated using two approaches: lines of action directly between the origin and insertion sites of the muscles, and lines of action along the centroid of the muscle bodies. Incorporating anatomic muscle paths into the model reduced the predicted peak quadriceps force during the chair-rise activity by 46%, and reduced the peak tibio-femoral contact pressure by 14%. In addition, bone strain was predicted during the activity for the implanted patella, and showed peak bone strain at the edge of the implant near the inferior pole. The muscle-activated models demonstrated the advantages of an explicit finite element framework, and allow rapid, rigid body simulation in addition to the full contact, deformable analyses when greater resolution is required

Introduction. Non-union is debilitating, costly and affects 2–8% of intramedullary fixed fractures. Clinical data suggest that percutaneous interfragmentary screws offer a less invasive alternative to exchange nailing. This study aimed to assess their efficiency with biomechanical analyses. Materials and Methods. A tibia was prepared for finite element analysis by creating a fracture of AO classification 42A2b, prior to reaming and insertion of an intramedullary nail. A callus was modelled as granulation tissue and gait loads were applied. The model was validated against published data and with sensitivity studies. The effects of weightbearing, fracture gap and angle, percutaneous screws and exchange nailing were compared through quantification of interfragmentary motion and strain, with the latter used to gauge healing performance via mechano-regulation theory. Results. Axial interfragmentary motion increased with increasing weightbearing, however, shear decreased at 25–50% weightbearing, leading to superior healing performance. Fracture gap had minimal effect on axial motion, but larger gaps gave greater shear, compromising healing. Elevated fracture obliquity culminated in more shear and inferior healing. Exchange nailing reduced axial motion by ∼30%, but had little effect on shear. Conversely, percutaneous screws had negligible effect on axial motion, but reduced shear by ∼15%, with three screws having a similar net effect on healing as exchange nailing from 10 to 11mm. Conclusions. This study provides new insight into fracture healing biomechanics and discovered that partial weightbearing, less oblique fractures and percutaneous screws all reduce shear, enhancing healing

Introduction. On the basis of a proposal by Noble, the marrow cavity form can be classified into three categories: stovepipe, normal, and champagne-fluted. In the present study, three typical finite element femoral models were created using CT data based on Noble's three categories. The purpose was to identify the relationship of stress distribution of the surrounding areas between femoral bone marrow cavity form and hip stem. The results shed light on whether the distribution of the high-stress area reflects the stem design concept. In order to improve the results of THA, researchers need to consider the instability of a stem design based on the pressure zone and give feedback on future stem selection. Methods. To develop finite element models, two parts (cortical bone and stem) were constructed using four-node tetrahedral elements. The model consisted of about 40,000 elements. The material characteristics were defined by the combination of mass density, elastic coefficient, and Poisson's ratio. Concerning the analysis system, HP Z800 Workstation(HP, Japan) was used as hardware and LS-DYNA Ver. 971 (Livermore Software Technology Corporation, USA) as software. The distal end of the femur was constrained in all directions. On the basis of ISO 7206 Part 4,8 that specifies a method of endurance testing for joint prostheses, the stem was tilted 10°, and a 500 N resultant force in the area around the hip joint was applied to the head at an angle of 25° with the long axis. Automatic contact with a consideration of slip was used. Von Mises stress during a 1.0 s period after loading was analyzed, and stress distribution in the stem and its maximum value were calculated. Result. The maximum stress at marrow cavity form of normal was shown to be 72 MPa. The stress of champagne-fluted was evenly distributed from proximal to distal, and the maximum stress was 67 MPa. For stovepipe, the maximum proximal stress was shown to be 120 MPa; moreover, stress concentration was observed. Discussion. The design concept for a Zweymüller-type stem can distribute load across a wide range of cortical bone from the middle position to the distal femur. It is determined using this concept that a wide range of stress was absorbed at the middle position and distal femur in the champagne-fluted and normal cases. On the other hand, the contact pressure zone of stovepipe could not meet the expected level at the distal femur. The method of this research involves controlling the stress conditions within the stem design. At this point, it is considered possible for the stability of various stem designs to be predicted and the stability to be assessed positively. On the basis of Noble's categories, three types of finite element model were made, and stress distribution measurement and finite element analyses were performed. The results indicate that Zweymüller stem has clinical validity for securing force in the champagne-fluted and stovepipe types from the stress distribution

Polyethylene wear of joint replacements can cause severe clinical complications, including; osteolysis, implant loosening, inflammation and pain. Wear simulator testing is often used to assess new designs, but it is expensive and time consuming. It is possible to predict the volume of polyethylene implant wear from finite element models using a modification of Archard's classic wear law [1–2]. Typically, linear elastic isotropic, or elasto-plastic material models are used to represent the polyethylene. The purpose of this study was to investigate whether use of a viscoelastic material model would significantly alter the predicted volumetric wear of a mobile-bearing unicompartmental knee replacement. Tensile creep-recovery experiments were performed to characterise the creep and relaxation behaviour of the polyethylene (moulded GUR 4150 samples machined to 180×20×1 mm). Samples were loaded to 3 MPa stress in 4 minutes, and then held for 6 hours, the tensile stress was removed and samples were left to relax for 6 hours. The mechanical test data was used fit to a validated three–dimensional fractional Maxwell viscoelastic constitutive material model [3]. An explicit finite element model of a mobile–bearing unicompartmental knee replacement was created, which has been described previously [4]. The medial knee replacement was loaded to 1200 N over a period of 0.2 s. The bearing was meshed using quadratic tetrahedral elements (1.5 mm seeding size based on results of a mesh convergence study), and the femoral component was represented as an analytical rigid body. Wear predictions were made from the contact stress and sliding distance using Archard's law, as has been described in the literature [1–2]. A wear factor of 5.24×10. −11. was used based upon the work by Netter et al. [2]. All models were created and solved using ABAQUS finite element software (version 6.14, Simulia, Dassault Systemes). The fractional viscoelastic material model predicted almost twice as much wear (0.119 mm. 3. /million cycles) compared to the elasto-plastic model (0.069 mm. 3. /million cycles). The higher wear prediction was due to both an increased sliding distance and higher contact pressures in the viscoelastic model. These preliminary findings indicate the simplified elasto-plastic polyethylene material representation can underestimate wear predictions from numerical simulations. Polyethylene is known to be a viscoelastic material which undergoes creep clinically, and it is not surprising that it is necessary to represent that viscoelastic behaviour to accurately predict implant wear. However, it does increase the complexity and run time of such computational studies, which may be prohibitive

Finite element analysis (FEA) has been applied for the biomechanical analysis of acetabular dysplasia, but not for biomechanical studies of periacetabular osteotomy (PAO) or those performing analysis taking into consideration the severity of acetabular dysplasia. This study aimed to perform biomechanical evaluation of changes in stress distribution following PAO and to determine the effect of the severity of developmental dysplasia of the hip (DDH) using three-dimensional FEA. A normal model was designed with a 25° center-edge (CE) angle and a 25° vertical-center-anterior margin (VCA) angle. DDH models were designed with CE and VCA angles each of 10, 0, or −10°. Post-PAO models were created by separating each DDH model and rotating the acetabular bone fragment in the anterolateral direction so that the femoral head was covered by the acetabular bone fragment, with CE and VCA angles each at 25°. Compared to the normal hip joint model, the DDH models showed stress concentration in the acetabular edge and contacting femoral head, and higher stress values; stress increased with decreasing CE and VCA angles. Compared to the DDH models, the post-PAO models showed near-normal patterns of stress distribution in the acetabulum and femoral head, with stress concentration areas shifted from the lateral to medial sides. Stress dispersion was especially apparent in the severe acetabular dysplasia models. PAO provided greater decreases in the maximum values of von Mises stress in the load-bearing area of the acetabulum and femoral head when applied to the DDH models of higher degrees of severity, although the values increased with increasing severity of DDH. PAO is expected to provide biomechanical improvement of the hip joint, although the results also suggest a limitation in the applicability of PAO for the patients with severe acetabular dysplasia

Introduction. On the basis of a proposal by Noble, the marrow cavity form can be classified into three categories: normal, champagne-fluted and stovepipe. In the present study, three typical finite element femoral models were created using CT data based on Noble's three categories. The purpose was to identify the relationship of stress distribution of the surrounding areas between femoral bone marrow cavity form and hip stems. The results shed light on whether the distribution of the high-stress area reflects the stem design concept. In order to improve the results of THA, researchers need to consider the instability of a stem design based on the stress distributioin and give feedback on future stem selection. Methods. As analyzing object, we selected SL-PLUS and BiCONTACT stems. To develop finite element models, two parts (cortical bone and stem) were constructed using four-node tetrahedral elements. The model consisted of about 60,000 elements. The material characteristics were defined by the combination of mass density, elastic coefficient, and Poisson's ratio. Concerning the analysis system, HP Z800 Workstation was used as hardware and LS-DYNA Ver. 971 as software. The distal end of the femur was constrained in all directions. On the basis of ISO 7206 Part 4,8 that specifies a method of endurance testing for joint prostheses, the stem was tilted 10°, and a 1500 N resultant force in the area around the hip joint was applied to the head at an angle of 25° with the long axis. Automatic contact with a consideration of slip was used. Result. The maximum stress on femur implanted a SL-PLUS with marrow cavity form of normal, champagne-fluted and stovepipe were shown to be 90MPa, 90MPa and 45MPa. The maximum stress on a BiCONTACT with marrow cavity form of normal, champagne-fluted and stovepipe were shown to be 45MPa, 90MPa and 15MPa. Discussion. The design concept for aZweymüller-type stem can distribute load across a wide range of cortical bone from the middle position to the distal femur. It is determined using this concept that a wide range of stress was absorbed at the middle position and distal femur in the champagne-fluted and normal cases. On the other hand, the contact pressure zone of stovepipe could not meet the expected level at the distal femur. The method of this research involves controlling the stress conditions within the stem design. At this point, it is considered possible for the stability of various stem designs to be predicted and the stability to be assessed positively. On the basis of Noble's categories, three types of finite element model were made, and stress distribution measurement and finite element analyses were performed. The results indicate that Zweymüller stem has clinical validity for securing force in the champagne-fluted and stovepipe types from the stress distribution