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

Introduction. The performance of total hip replacement (THR) devices can be affected by the quality of the tissues surrounding the joint or the mismatch of the component centres during hip replacement surgery. Experimental studies have shown that these factors can cause the separation of the two components during walking cycle (dynamic separation) and the contact of the femoral head with the rim of the acetabular liner (edge loading), which can lead to increased wear and shortened implant lifespan. 1. There is a need for flexible pre-clinical testing tools which allow THR devices to be assessed under these adverse conditions. In this work, a novel dynamic finite element model was developed that is able to generate dynamic separation as it occurs during the gait cycle. In addition, the ability to interrogate contact mechanics and material strain under separation conditions provides a unique means of assessing the severity of edge loading. This study demonstrates these model capabilities for a range of simulated surgical translational mismatch values, for ceramic-on-polyethylene implants. Methodology. The components of the THR were aligned and constrained as illustrated in Figure 1. CAD models of commercially available implant geometries were used (DePuy Synthes, Leeds, UK) modified for model simplicity by removing anti-rotation features. The polyethylene cup liner was given elastic-plastic behaviour. An axial load following the Paul cycle pattern (5 repetitive cycles) with maximum of 3KN and swing phase load of 0.3KN, was applied through the cup holder. The effect of translational mismatch was implemented by using a spring element connected to the cup unit on the lateral side. The spring was compressed by a fixed amount to replicate a degree of medial-lateral mismatch of the components. The instantaneous resultant force vector dictated the dynamic sliding behaviour of the cup against the head. In this study, translational medial-lateral mismatch values of 1, 2, 3 and 4mm were used and the medial-lateral dynamic separation, contact pressure maps and plastic strain were recorded. Results. The highest level of dynamic separation is achieved when the minimum axial load (during swing phase) is applied. The dynamic separation increases as the surgical translation mismatch increases (figure 2), with values over 0.5mm (radial clearance) representing cases where the head is in contact with the rim of the cup. Maximum separation occurred towards the end of the swing phase. Plots of the shape of the contact pressure at that point can be seen in Figure 3. Only the 4mm mismatch created substantial plastic deformation. Conclusion. The finite element model was able to predict medial-lateral separation as it occurred dynamically in the gait cycle, including cases where the femoral head was in contact with the rim of the cup. The increase in medial-lateral separation with increased translational mismatch was in broad agreement with existing experimental data. 2. Substantial plastic deformation was only seen in cases where the translational mismatch caused the femoral head to be in contact with the rim of the polyethylene cup

Background. It is known that severe cases of intervertebral disc (IVD) disease may lead to the loss of natural intervertebral height, which can cause radiating pain throughout the lower back and legs. To this point, surgeons perform lumbar fusion using interbody cages, posterior instrumentation and bone graft to fuse adjacent vertebrae together, thus restoring the intervertebral height and alleviating the pain. However, this surgical procedure greatly decreases the range of motion (ROM) of the treated segment, mainly caused by high cage stiffness. Additive manufacturing can be an interesting tool to reduce the cage's elastic modulus (E), by adding porosity (P) in its design. A porous cage may lead to an improved osteointegration since there is more volume in which bone can grow. This work aims to develop a finite element model (FEM) of the L4-L5 functional spinal unit (FSU) and investigate the loss of ROM induced by solid and porous cages. Materials and Methods. The Intact-FEM of L4-L5 was created, which considered the vertebrae, IVD and ligaments with their respective material properties. 1. The model was validated by comparing its ROM with that of other studies. Moments of 10 Nm were applied on top of L4 while the bottom of L5 was fixed to simulate flexion, extension, lateral bending and axial rotation. 2. The lumbar cages, posterior instrumentation and bone graft were then modelled to create the Cage-FEMs. Titanium was chosen for the instrumentation and cages. Cages with different stiffness were considered to represent porous structures. The solid cage had the highest modulus (E. 0. =110 GPa, P. 0. =0%) whereas the porous cages were simulated by lowering the modulus (E. 1. =32.8 GPa, P. 1. =55%; E. 2. =13.9 GPa, P. 2. =76%; E. 3. =5.52 GPa, P. 3. =89%; E. 4. =0.604 GPa, P. 4. =98%), following the literature. 3. The IVD was removed in Cage-FEMs to allow the implant's insertion [Fig. 1] and the previous loading scenarios were simulated to assess the effects of cage porosity on ROM. Results. The Intact-FEM presents acceptable ROM according to experimental and numerical studies, as shown by the red line in Figure 2. After insertion, lower ROM values in Cage-FEMs are measured for each physiological movement [Fig. 3]. In addition, highly porous cages have greater ROM, especially in axial rotation. Discussion. Significant reduction of ROM is expected after cage insertion because the main goal of interbody fusion is to allow bone growth. As such, the procedure's success is highly dependent on segmental stability, which is achieved by using cages in combination with bone graft and posterior instrumentation. Furthermore, higher cage porosities seem to affect the FSU. In fact, ROM increases more as the cage modulus approaches that of the cancellous bone (E. canc-bone. =0.2 GPa. 1. ). Next step will be to assess the effects of cage design on the L4-L5 FSU mechanical behavior and stress distribution. To conclude, additive manufacturing offers promising possibilities regarding implant optimization, being able to create porous cages, thus reducing their stiffness. For any figures or tables, please contact the authors directly

Introduction

Fretting corrosion at the taper interface of modular connections can be studied using Finite Element (FE) analyses. However, the loading conditions in FE studies are often simplified, or based on generic activity patterns. Using musculoskeletal modeling, subject-specific muscle and joint forces can be calculated, which can then be applied to a FE model for wear predictions. The objective of the current study was to investigate the effect of incorporating more detailed activity patterns on fretting simulations of modular connections.

INTRODUCTION. In theory, Finite Element Analysis (FEA) is an attractive method for elucidating the mechanics of modular implant junctions, including variations in materials, designs, and modes of loading. However, the credence of any computational model can only be established through validation using experimental data. In this study we examine the validity of such a simulation validated by comparing values of interface motion predicted using FEA with values measured during experimental simulation of stair-climbing. MATERIALS and METHODS. Two finite element models (FEM) of a modular implant assembly were created for use in this study, consisting of a 36mm CoCr femoral head attached to a TiAlV rod with a 14/12 trunnion. Two head materials were modelled: CoCr alloy (118,706 10-noded tetrahedral elements), and alumina ceramic (124,710 10-noded tetrahedral elements). The quasi-static coefficients of friction (µ. s. ) of the CoCr-TiAlV and Ceramic-TiAlV interfaces were calculated from uniaxial assembly (2000N) and dis-assembly experiments performed in a mechanical testing machine (Bionix, MTS). Interface displacements during taper assembly and disassembly were measured using digital image correlation (DIC; Dantec Dynamics). The assembly process was also simulated using the computational model with the friction coefficient set to µ. s. and solved using the Siemens Nastran NX 11.0 Solver. The frictional conditions were then varied iteratively to find the value of µ providing the closest estimate to the experimental value of head displacement during assembly. To validate the FEA model, the relative motion between the head and the trunnion was measured during dynamic loading simulating stair-climbing. Each modular junction was assembled in a drop tower apparatus and then cyclically loaded from 230–4300N at 1 Hz for a total of 2,000 cycles. The applied load was oriented at 25° to the trunnion axis in the frontal plane and 10° in the sagittal plane. The displacement of the head relative to the trunnion during cyclic loading was measured by a three-camera digital image correlation (DIC) system. The same loading conditions were simulated using the FEA model using the optimal value of µ derived from the initial head assembly trials. RESULTS. For both head materials, the predicted values of axial displacement of the head on the trunnion closely approximated the measured values derived from DIC measurements, with differences of −0.17% to +6.5%, respectively. Larger differences were calculated for individual components of motion for the stair climbing activity. However, the predicted magnitude of interface motion was still within 10% of the observed values, ranging from −7% to −5%. CONCLUSIONS. Our simulations closely approximated physical testing using complex loading, coming within 7% of the target values. By generating a validated computational model of a modular junctions with varying head materials, we will be able to simulate additional activities of daily living to determine micromotion and areas of peak pressure and contact stresses generated. For any figures or tables, please contact authors directly

Introduction. Aligning the tibial tray is a critical step in total knee arthroplasty (TKA). Malalignment, (especially in varus) has been associated with failure and revision surgery. While the link between varus malalignment and failure has been attributed to increased medial compartmental loading and generation of shear stress, quantitative biomechanical evidence to directly support this mechanism is incomplete. We therefore constructed and validated a finite element model of knee arthroplasty to test the hypothesis that varus malalignment of the tibial tray would increase the risk of tray subsidence. Methods. Cadaver Testing. Fresh human knees (N = 4) were CT scanned and implanted with TKA cruciate-retaining tibial tray (Triathlon CR, Stryker Orthopaedics, New Jersey). The specimens were subjected to ISO-recommended knee wear simulation loading for up to 100,000 cycles. Micromotion sensors were mounted between the tray and underlying bone to measure micromotion. In two of the specimens, the application of vertical load was shifted medially to generate a load distribution ratio of 55:45 (medial:lateral) to represent neutral varus-valgus alignment. In the remaining two specimens, a load distribution ratio of 75:25 was generated to represent varus alignment. Finite element analysis. qCT scans of the tested knees were segmented using MIMICS (Materialise, Belgium). Material properties of bone were spatially assigned after converting bone density to elastic modulus. A finite element model of the tibia implanted with a tibial tray was constructed (Abaqus 6.8, Simulia, Dassault Syst`mes). Boundary conditions were applied to simulate experimental mounting conditions and the tray was subjected to a single load cycle representing that applied during cadaver loading. Results. The two cadaver specimens tested at 55:45 medial:lateral (M:L) force distribution survived the 100,000 cycle test, while both cadaver specimens tested at 75:25 M:L force distribution failed. The finite element model generated distinct differences in compressive strain distribution patterns in the proximal tibia. A threshold of 2000 microstrain was used for fatigue damage in bone under cyclic loading. Both specimens loaded under 75:25 M:L distribution demonstrated substantially larger cortical bone volumes in the proximal tibial cortex that were greater than this fatigue threshold. Discussion & Conclusion. We validated a finite element model of tibial loading after TKA. Local compressive strains directly correlated with subsidence and failure in cadaver testing. A significantly greater volume of proximal tibial cortical bone was compressed to a strain greater than the fatigue threshold in the varus alignment group, indicating an increased risk for fatigue damage. This model is extremely valuable in studying the effect of surgical alignment, loading, and activity on damage to proximal bone. Emerging techniques that customize tibial tray placement to the individual patient's pre-arthritic alignment run counter to the traditional recommendations for coronal alignment to the mechanical axis of the knee. A method that determines the risk of bone damage in a patient-specific manner can provide the surgeon with a safe range for component alignment and may even be applicable in preoperative planning

Introduction. The process of wear and corrosion at the head-neck junction of a total hip replacement is initiated when the femoral head and stem are joined together during surgery. To date, the effects of the surface topography of the femoral head and metal stem on the contact mechanics during assembly and thus on tribology and fretting corrosion during service life of the implant are not well understood. Therefore, the objective of this study was to investigate the influence of the surface topography of the metal stem taper on contact mechanics and wear during assembly of the head-neck junction using Finite Element models. Materials and Methods. 2D axisymmetric Finite Element models were developed consisting of a simplified head-neck junction incorporating the surface topography of a threaded stem taper to investigate axial assembly with 1 kN. Subsequently, a base model and three modifications of the base model in terms of profile peak height and plateau width of the stem taper topography and femoral head taper angle were calculated. To account for the wear process during assembly a law based on the Archard equation was implemented. Femoral head was modeled as ceramic (linear-elastic), taper material was either modeled as titanium, stainless steel or cobalt-chromium (all elastic-plastic). Wear volume, contact area, taper subsidence, equivalent plastic strain, von Mises stress, engagement length and crevice width was analyzed. Results. Titanium tapers showed largest wear volume throughout all simulations, followed by stainless steel and cobalt-chromium. A larger head taper angle resulted in an increase of the wear volume for all taper materials while the increase of the plateau width resulted in a decrease of the wear volume. Taper subsidence, von Mises stress and equivalent plastic strain followed the same trends. Contact area was largest for the models with a large plateau width for all taper materials. Other taper parameters had little effect on contact area. A pure increase of the angular mismatch (AM) resulted in the strongest decrease of the engagement length, while a combined increase of the AM and plateau width showed only a moderate decrease. The smallest effect concerning the engagement length was found when a combined increase of the profile peak height and AM was simulated. Crevice width was largest for a pure increase of the AM and for a combined increase of the AM and profile peak height for all taper materials. Discussion. This study showed that depending on the surface topography and material of the stem taper, wear and taper mechanics during assembly could be affected. For the examined surface topographies wear is distinctively elevated by increasing the AM and the profile peak height due to the resulting higher mechanical loading. More parameter studies under in vivo loading and the study of other taper surface parameters like the peak-to-peak distance have to be conducted to get a deeper insight into taper mechanics and wear effects. However, this study demonstrates the importance of good manufacturing practice of components for hip replacement systems to guarantee reproducible taper mechanics. For any figures or tables, please contact authors directly

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

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: The mechanobiology and response of bone formation to strain under physiological loading is well established, however investigation into exceedingly soft scaffolds relative to cancellous bone is limited. In this study we designed and 3D printed mechanically-optimised low-stiffness implants, targeting specific strain ranges inducing bone formation and assessed their biological performance in a pre-clinical in vivo load-bearing tibial tuberosity advancement (TTA) model. The TTA model provides an attractive pre-clinical framework to investigate implant osseointegration within an uneven loading environment due to the dominating patellar tendon force. A knee finite element model from ovine CT data was developed to determine physiological target strains from simulated TTA surgery. We 3D printed low-stiffness Ti wedge osteotomy implants with homogeneous stiffness of 0.8 GPa (Ti1), 0.6 GPa (Ti2) and a locally-optimised design with a 0.3 GPa cortex and soft 0.1 GPa core (Ti3), for implantation in a 12-week ovine tibial advancement osteotomy (9mm). We quantitatively assessed bone fusion, bone area, mineral apposition rate and bone formation rate. Optimised Ti3 implants exhibited evenly high strains throughout, despite uneven wedge osteotomy loading. We demonstrated that higher strains above 3.75%, led to greater bone formation. Histomorphometry showed uniform bone ingrowthin optimised Ti3 compared to homogeneous designs (Ti1 and Ti2), and greater bone-implant contact. The greatest bone formation scores were seen in Ti3, followed by Ti2 and Ti1. Results from our study indicate lower stiffness and higher strain ranges than normally achieved in Ti scaffolds stimulate early bone formation. By accounting for loading environments through rational design, implants can be optimised to improve uniform osseointegration. Design and 3D printing of exceedingly soft titanium orthopaedic implants enhance strain induced bone formation and have significant importance in future implant design for knee, hip arthroplasty and treatment of large load-bearing bone defects

Introduction. Geometric variations of the hip joint can give rise to repeated abnormal contact between the femur and acetabular rim, resulting in cartilage and labrum damage. Population-based geometric parameterisation can facilitate the flexible and automated in silico generation of a range of clinically relevant hip geometries, allowing the position and size of cams to be defined precisely in three dimensions. This is advantageous compared to alpha angles, which are unreliable for stratifying populations by cam type. Alpha angles provide an indication of cam size in a single two-dimensional view, and high alpha angles have been observed in asymptomatic individuals. Parametric geometries can be developed into finite element models to assess the potential effects of morphological variations in bone on soft tissue strains. The aim of this study was to demonstrate the capabilities of our parameterisation research tool by assessing impingement severity resulting from a range of parametrically varied femoral and acetabular geometries. Methods. Custom made MATLAB (MathWorks) and Python codes. [1]. were used to generate bone surfaces, which were developed into finite element models in Abaqus (SIMULIA). Parametric femoral surfaces were defined by a spherical proximal head and ellipse sections through the neck/cam region. This method produced surfaces that were well fitted to bone geometry segmented from CT scans of cam patients and capable of producing trends in results similar to those found using segmented models. A simplified spherical geometry, including the labrum and acetabular cartilage, represented the acetabulum. Femoral parameters were adjusted to define relevant variations in cam size and position. Two radii (small and large cams) and two positions (anterior and superior cams) were defined resulting in four models. Alpha angles of these parametric femurs were measured in an anterior-posterior view and a cross-table lateral view using ImageJ (NIH). A further model was developed using a femur with a medium cam size and position, and the level of acetabular coverage and labrum length were varied. Bones were modelled as rigid bodies and soft tissues were modelled as transversely isotropic linearly elastic materials. With the acetabulum fully constrained in all cases, the femurs were constrained in translation and rotated to simulate flexion followed by internal rotation to cause impingement against the labrum. Results and Discussion. Models generated using the parametric approach showed that potential for tissue damage, indicated through local strain, was not predicted by measured alpha angle, but resulted from cam extent and position as defined by the ellipses. When variations were made to the acetabular rim, an increase in bone coverage had the greatest effect on impingement severity, indicated by strain in the cartilage labral-junction. An increase in labral length increased labral displacement, but had less effect on cartilage-labral strain. Patient specific models currently require full image segmentation, but there is potential to further develop these parametric methods to assess likely impingement severity based on a series of measures of the neck and acetabulum when three-dimensional imaging of patients is available

The advantages of unicompartmental knee arthroplasty (UKA) include its bone preserving nature, lower relative cost and superior functional results. Some temporary pain has been reported clinically following this procedure. Could this be related to bone remodeling? A validated bone remodeling algorithm may have the answers…. A 3D geometry of an intact human cadaveric tibia was generated using CT images. An all poly unicompartmental implant geometry was positioned in an inlay and onlay configuration on the tibia and the post-operative models created. An adaptive bone remodeling algorithm was used with finite element modeling to predict the bone remodeling behavior surrounding the implant in both scenarios. Virtual DEXA images were generated from the model and bone mineral density (BMD) was measured in regions of interest in the AP and ML planes. BMD results were compared to clinical results. The bone remodelling algorithm predicted BMD growth in the proximal anterior regions of the tibia, with an inward tendency for both inlay and onlay models. Looking in the AP plane, a maximum of up to 7% BMD growth was predicted and in the ML plane this was as high as 16%. Minimal BMD loss was observed, which suggests minimal disturbance to the natural bone growth following UKA. Positron emission tomography (PET) scans showed active hot spots in the antero- medial regions of the tibia. These results were consistent with the finite element modeling results. Bone remodeling behavior was found to be sensitive to sizing and positioning of the implant. The adaptive bone remodeling algorithm predicted minimal BMD loss and some BMD growth in the anterior region of the tibia following UKA. This is consistent with patient complaint and PET scans

Distal radius fractures are the most common osteoporotic fractures among women. The treatment of these fractures has been shifting from a traditional non-operative approach to surgery, using volar locking plate (VLP) technology. Surgery, however, is not without risk, complications including failure to restore an anatomic reduction, fracture re-displacement, and tendon rupture. The VLP implant is also marked by bone loss due to stress-shielding related to its high stiffness relative to adjacent bone. Recently, a novel internal, composite-based implant, with a stiffness less than the VLP, was designed to eradicate the shortcomings associated with the VLP implant. It is unclear, however, what effect this less-stiff implant will have upon adjacent bone density distributions long-term. The objective of this study was to evaluate the long-term effects of the two implants (the novel surgical implant and the gold-standard VLP) by using subject-specific finite element (FE) models integrated with an adaptive bone formation/resorption algorithm. Specimen: One fresh-frozen human forearm specimen (female, age = 84 years old) was imaged using CT and was used to create a subject-specific FE model of the radius. Finite element modeling: In order to simulate a clinically relevant (unstable) fracture of the distal radius, a wedge of bone was removed from the model, which was approximately 10 mm wide and centered 20 mm proximal to the tip of the radial styloid. Bone remodeling algorithm: A strain-energy density (SED) based bone remodeling theory was used to account for bone remodeling. With this approach, bone density decreased linearly when SED per bone density was less than 67.5 µJ/g and increased when it was more than 232.5 µJ/g. When it was in the lazy zone (67.5 to 232.5 µJ/g), no changes in density occurred. Boundary conditions: A 180 N quasi-static force representing the scaphoid, and a 120 N quasi-static force representing the lunate was applied to the radius. The midshaft of the radius was constrained. FE outcomes: To examine the effects of stress shielding associated with each implant, the long-term changes of bone density within proximal transverse cross-sections of radius were inspected. The regional density analysis focused on three transverse cross-sections. The transverse cross-sections were positioned proximal to the subchondral plate, and were distanced 50 (cross-section A), 57 (cross-section B), and 64 mm (cross-section C) from the subchondral endplate. For both implants in all three cross-sections, cortical bone was reserved completely at the volar side. On the dorsal side, the cortical bone was completely resorbed in the VLP model. In all cross-sections, the averaged resultant density was higher for the “novel implant”. The difference ranged from 33% (cross-section A) to 36% (cross-section C) in favor of the “novel implant”. On average, the density values of the novel implant were 34% higher in transverse cross-sections (A, B, and C). This study showed that the novel implant offered higher density distributions compared to the VLP, which suggests that the novel implant may be superior to the VLP in terms of avoiding stress shielding

Revision surgeries for orthopaedic infections are done in two stages – one surgery to implant an antibiotic spacer to clear the infection and another to install a permanent implant. A permanent porous implant, that can be loaded with antibiotics and allow for single-stage revision surgery, will benefit patients and save healthcare resources. Gyroid structures can be constructed with high porosity, without stress concentrations that can develop in other period porous structures [1] [2]. The purpose of this research is to compare the resulting bone and prosthesis stress distributions when porous versus solid stems are implanted into three proximal humeri with varying bone densities, using finite element models (FEM). Porous humeral stems were constructed in a gyroid structure at porosities of 60%, 70%, and 80% using computer-aided design (CAD) software. These CAD models were analyzed using FEM (Abaqus) to look at the stress distributions within the proximal humerus and the stem components with loads and boundary conditions representing the arm actively maintained at 120˚ of flexion. The stem was assumed to be made of titanium (Ti6Al4V). Three different bone densities were investigated, representing a healthy, an osteopenic, and an osteoporotic humerus, with an average bone shape created using a statistical shape and density model (SSDM) based on 75 cadaveric shoulders (57 males and 18 females, 73 12 years) [3]. The Young's moduli (E) of the cortical and trabecular bones were defined on an element-by-element basis, with a minimum allowable E of 15 MPa. The Von Mises stress distributions in the bone and the stems were compared between different stem scenarios for each bone density model. A preliminary analysis shows an increase in stress values at the proximal-lateral region of the humerus when using the porous stems compared to the solid stem, which becomes more prominent as bone density decreases. With the exception of a few mesh dependent singularities, all three porous stems show stress distributions below the fatigue strength of Ti-6Al-4V (410 MPa) for this loading scenario when employed in the osteopenic and osteoporotic humeri [4]. The 80% porosity stem had a single strut exceeding the fatigue strength when employed in the healthy bone. The results of this study indicate that the more compliant nature of the porous stem geometries may allow for better load transmission through the proximal humeral bone, better matching the stress distributions of the intact bone and possibly mitigating stress-shielding effects. Importantly, this study also indicates that these porous stems have adequate strength for long-term use, as none were predicted to have catastrophic failure under the physiologically-relevant loads. Although these results are limited to a single boney geometry, it is based on the average shape of 75 shoulders and different bone densities are considered. Future work could leverage the shape model for probabilistic models that could explore the effect of stem porosity across a broader population. The development of these models are instrumental in determining if these structures are a viable solution to combatting orthopaedic implant infections

Introduction. Relative motion at the modular head-neck junction of hip prostheses can lead to severe surface damage through mechanically-assisted corrosion. One factor affecting the mechanical performance of modular junctions is the frictional resistance of the mating surfaces to relative motion. Low friction increasing forces normal to the head-neck interface, leading to a lower threshold for slipping during weight-bearing. Conversely, a high friction coefficient is expected to limit interface stresses but may also allow uncoupling of the interface in service. This study was performed to examine this trade-off using finite element models of the modular head-neck junction. Methods. A finite element model (FEM) of the trunnion/ head assembly of a total hip prosthesis was initially created and experimentally validated. CAD models of a stem trunnion (taper size: 12/14mm) and a prosthetic femoral head (diameter: 28mm) were discretized into elements for finite element analysis (FEA). The trunnion (Ti6Al4V) was modelled with a hexahedral mesh (33,648 elements) and the femoral head (CoCrMo) with a tetrahedral mesh (51,182 elements). A friction-based sliding contact interface was defined between the mating surfaces. The model was loaded in 2 stages: (i) an assembly load of 4000N applied along the trunnion axis, and (ii) 500N applied along the trunnion axis in combination with a torque of 10Nm. A linear static solution was set up using Siemens NX-Nastran solver. Multiple simulations were executed by modulating the frictional coefficient at the taper-bore interface from 0.05 to 0.15 in increments of 0.01, the coefficient of 0.1 serving as the control case (Swaminathan and Gilbert, 2012). Results. The vertical and tangential displacements of the nodes on the taper of the trunnion relative to the femoral head demonstrated a strong inverse dependence upon the coefficient of friction at the interface (Fig. 1). A similar trend was observed with respect to the peak interface pressure (Fig. 2). The peak von Mises stress, however, increases with increasing coefficient of friction (Fig. 2). A Fisher's R to Z correlation test was performed on each output variable to determine its correlation with coefficient of friction. The coefficient of friction correlated significantly (p<0.0001) with both tangential displacement (r = −0.990) and vertical displacement (r = −0.974). Peak von Mises stress (r = 0.995) and peak contact pressure (r = −0.984) were also found to be significantly (p<0.0001) correlated to the coefficient of friction. Discussion. A higher coefficient of friction at the taper-bore interface led to lower contact pressure and sliding at the modular junction. However, higher coefficients of friction also led to increased von Mises stresses within the bore and the trunnion increasing the risk of yielding and fatigue failure. The current results strongly indicate that factors affecting the frictional coefficient at the interface likely influence the occurrence of and severity of mechanically-assisted corrosion in THA. Significance. The results from this study will help us set tolerances for the interlocking mechanism, identifying the minimum frictional coefficient required to obtain stable implant mechanics

Mechanical wear and corrosion lead to the release of metal particulate debris and subsequent release of metal ions at the trunnion-taper surface. In order to quantify the amount of volume loss to ultimate locations in the surrounding joint space, finite element analysis of the modular head-stem junction is being carried out. The key purpose being to determine a set of optimum design changes that offer the least material loss at the taper-trunnion junction using optimization algorithms such as the gradient based local search (Sequential Quadratic Programming–SQP) and global search (Non-Dominated Sorting Genetic Algorithm-II–NSGA-II). In a broader sense, the principal goal is to work toward the minimization of wear debris produced in the hip joint, thereby resulting in a longer prosthetic lifetime. A numerical approach that simulates wear in modular hip prostheses with due consideration to the taper-trunnion junction on metal-on-metal contacts is proposed. A quasi-static analysis is performed considering realistic loading stages in the gait cycle, and nonlinear contact analysis is to be employed. The technique incorporates a measured wear rate as an input to the finite element model. The simulation of wear is performed by progressively changing nodal coordinates to simulate the wear loss that occurs during surface interaction. The geometry of the worn surface is updated under gait loading. With a given geometry and gait loading, the linear and volumetric wear increases with the number of gait cycles. The continuous wear propagation is discretized and an approximation scheme known as surrogates is to be developed using Artificial Neural Networks (ANN) to reduce the expensive computational simulations during optimization. The model is employed in the optimization schemes coded in MATLAB and linked to the finite element model developed in ANSYS batch mode. The objective function of the optimization problem is to minimize the volumetric wear at taper-trunnion interface under some constraints. By minimizing the volumetric wear, the chance of failure of modular hip implants is also minimized. The FE model developed to reproduce fretting wear is validated through in vitro wear simulations. The important taper design variables considered to have impact on the fretting corrosion performance include; medial-lateral offset, neck length, taper head diameter, trunnion length and diameter, included angle for the head/neck tapers, angle of mismatch or variation in taper trunnion angle, etc. It is expected from clinical outcomes that increased offset and large taper diameter has serious implications in the fretting corrosion behavior primarily because these variables control the bending stresses and strains along the length of the taper. During cyclic loading of the taper, the higher the strain range, the higher will be the relative micromotion at the point of engagement between the stem and head tapers. This research is carried out with the objective to optimize the effects of these geometrical factors at the mating taper interfaces. The developed models have great potentiality for accurate assessment of wear in a range of metal-on-metal (MoM) hip prostheses at the femoral head taper-trunnion junction while substantially reducing the wear and failure rate of prostheses

Introduction:. The mechanical stresses and strains surrounding orthopaedic implants can influence bone resorption and formation, micro-fracture, and consequently implant fixation or loosening. Experimental measurement of these internal parameters is generally not feasible. Computational predictions by finite element modeling are promising, but until recently have been limited to assuming the surrounding cancellous bone as a continuous volume, without modeling individual trabeculae. A recent study demonstrated errors in bone-implant stiffness exceeding 100% when using this continuum assumption [1]. Conversely, recently micro-finite element computer models have been built from high resolution imaging of trabecular bone. In the present study we developed such models of central pegs cemented into cadaveric glenoids. We hypothesized that additional applied cement would lead to stronger implant fixation, but less physiologic strains in the trabeculae. Methods:. Two cadaveric specimens were implanted, with the applied cement volume in the Specimen 2 approximately double that of Specimen 1. The specimens were imaged by micro-computed tomography (vivaCT 40, Scanco, Switzerland) with a resolution of 12 microns. Images were filtered and resampled, then imported in Mimics (Materialise, Belgium) for semi-automated segmentation and 3D reconstruction based on our laboratory's published methods. Finite element models containing 1.7 to 1.8 million elements having sides of 0.1 mm were generated by a direct image voxel-to-element approach [2] (Fig. 1). The material properties of cement and bone were assumed linear elastic (bone: E = 3.5 GPa, cement: E = 3.0 GPa, and implant (UHMWPE): E = 1.3 GPa), and interfaces were assumed fully bonded. All outer walls of the bone were fixed, and a downward force of 250 N was applied to the implant peg. Simulations were run using Abaqus (Simulia, Pawtucket RI) on a 32-core, 1 TB-memory server at PSU's High Performance Computing Systems. Results:. Specimen 1 had 254 mm. 3. cement measured in the model, whereas Specimen 2 had 535 mm. 3. Strain energy density was less for Specimen 2 for bone underneath the implant, but similar between specimens for bone around the implant sides (Figs 2 and 3), providing initial indication of complex effects of cement volume on peri-implant strains. In Specimen 2 a slightly larger volume of cement (8.6 vs. 6.8 mm. 3. ) was exposed to von Mises stresses exceeding 10 MPa. Discussion:. This study is novel in its prediction of stresses and strains down to the level of individual glenoid trabeculae surrounding a cemented implant. In this pilot investigation we found that bone embedded in the cement mantle is subject to low strains, whereas the bone immediately surrounding the cement mantle is subject to abnormally high strains, with both cement technique and trabecular architecture likely influencing results. The study is limited by the lack of application of more complex loads and boundary conditions. Future work includes modeling of additional specimens and statistical analyses, and investigation of the roles of cement stiffness and peg design in dictating peri-implant bone strains

Accurate in vivo knee joint contact forces are required for joint simulator protocols and finite element models during the development and testing of total knee replacements (Varadarajan et al., 2008.) More accurate knowledge of knee joint contact forces during high flexion activities may lead to safer high flexion implant designs, better understanding of wear mechanisms, and prevention of complications such as aseptic loosening (Komistek et al., 2005.) High flexion is essential for lifestyle and cultural activities in the developing world, as well as in Western cultures, including ground-level tasks and chores, prayer, leisure, and toileting (Hemmerich et al., 2006.) In vivo tibial loads have been reported while kneeling; but only while the subject was at rest in the kneeling position (Zhao et al., 2007), meaning that the loads were submaximal due to muscle relaxation and thigh-calf contact support. The objective of this study was to report the in vivo loads experienced during high flexion activities and to determine how closely the measured axial joint contact forces can be estimated using a simple, non-invasive model. It provides unique data to better interpret non-invasively determined joint-contact forces, as well as directly measured tiobiofemoral joint contact force data for two subjects. Two subjects with instrumented tibial implants performed kneeling and deep knee bend activities. Two sets of trials were carried out for each activity. During the first set, an electromagnetic tracking system and two force plates were used to record lower limb kinematics and ground reaction forces under the foot and under the knee when it was on the ground. In the second set, three-dimensional joint contact forces were directly measured in vivo via instrumented tibial implants (Heinlein et al., 2007.) The measured axial joint contact forces were compared to estimates from a non-invasive joint contact force model (Smith et al., 2008.). The maximum mean axial forces measured during the deep knee bend were 24.2 N/kg at 78.2° flexion (subject A) and 31.1 N/kg at 63.5° flexion (subject B) during the deep knee bend (Figure 1.) During the kneeling activity, the maximum mean axial force measured was 29.8 N/kg at 86.8° flexion (subject B.) While the general shapes of the model-estimated curves were similar to the directly measured curves, the axial joint contact force model underestimated the measured contact forces by 7.0 N/kg on average (Figure 2.) The most likely contributor to this underestimation is the lack of co-contraction in the model. The study protocol was limited in that data could not be simultaneously collected due to electromagnetic interference between the motion tracking system and the inductively powered instrumented tibial component. Because skin-mounted markers were used, kinematics may be affected by skin motion artefacts. Despite these limitations, this study presents valuable information that will advance the development of high flexion total knee replacements. The study provides in vivo measurements and non-invasive estimates of joint contact forces during high flexion activities that can be used for joint simulator protocols and finite element modeling

Introduction. Successful designs of total hip replacement need to be robust to surgery-related variability. Until recently, only simple parametric studies have explored the influence of surgical variability [1]. This study presents a systematic method for quantifying the effect of variability in positioning on the primary stability of femoral stems using finite element (FE) models. Methods. Patient specific finite element models were generated of two femurs, one male and one female. An automated algorithm positioned and sized a Corail stem (DePuy Synthes, Warsaw) into each of the femurs to achieve maximum fill of the medullary canal without breaching into the cortical bone boundaries.. Peak joint contact and muscle forces associated with level gait were applied[2] and scaled to the body mass of each subject, whilst the distal femur was rigidly constrained. The space prone to surgical variation was defined by the “gap” between the stem and the inner boundary of the cortical bone. The anterior/posterior and the varus/valgus alignment of the stem within this “gap” was controlled by varying the location of the points defining the shaft axis. The points were taken at 20% and 80% of the stem length (Figure 1). The anteversion angle as well as the vertical and the medial position of the stem were controlled by changing the location of the head centre within the femoral head radius. The location of these points was varied using Latin Hypercube sampling to generate 200 models per femur, each with a unique stem position. The risk of failure was evaluated based on stem micromotion, equivalent strains, and percentage of the bone-prosthesis contact area experiencing more than 7000 µstrains [3]. Results. The range of positions covered in this study adhered to the anatomy of the subjects (Table 1) and none of the stem positions breached into the cortical bone of the femur. The 90th percentile peri-prosthetic strains were between 1770 – 4792 µstrains for the male subject, and 2710 – 11260µstrains for the female subject. The 90th percentile micromotion was between (15.6 – 47) µm for the male subject, and (42.4 – 102.4) µm for the female subject. The percentage of the contact area experiencing more than 7000 µstrains was between (0% – 0.33%) for the male subject, and (0% – 12%) for the female subject. Discussion. A systematic method for studying the effect of surgical-related variation on primary stability was presented its applicability demonstrated on two femurs. The study found that variation in stem position may result in large variation (up to 1.5 times the baseline position) in strains and micromotions. The magnitude Up to three times the magnitudes for the ideal stem position. This method can be applied to larger samples to understand the influence of different alignment parameters on the primary stability of femoral stems

Introduction. Varus alignment in total knee replacement (TKR) results in a larger portion of the joint load carried by the medial compartment. [1]. Increased burden on the medial compartment could negatively impact the implant fixation, especially for cementless TKR that requires bone ingrowth. Our aim was to quantify the effect varus alignment on the bone-implant interaction of cementless tibial baseplates. To this end, we evaluated the bone-implant micromotion and the amount of bone at risk of failure. [2,3]. Methods. Finite element models (Fig.1) were developed from pre-operative CT scans of the tibiae of 11 female patients with osteoarthritis (age: 58–77 years). We sought to compare two loading conditions from Smith et al.;. [1]. these corresponded to a mechanically aligned knee and a knee with 4° of varus. Consequently, we virtually implanted each model with a two-peg cementless baseplate following two tibial alignment strategies: mechanical alignment (i.e., perpendicular to the tibial mechanical axis) and 2° tibial varus alignment (the femoral resection accounts for additional 2° varus). The baseplate was modeled as solid titanium (E=114.3 GPa; v=0.33). The pegs and a 1.2 mm layer on the bone-contact surface were modeled as 3D-printed porous titanium (E=1.1 GPa; v=0.3). Bone material properties were non-homogeneous, determined from the CT scans using relationships specific to the proximal tibia. [2,4]. The bone-implant interface was modelled as frictional with friction coefficients for solid and porous titanium of 0.6 and 1.1, respectively. The tibia was fixed 77 mm distal to the resection. For mechanical alignment, instrumented TKR loads previously measured in vivo. [5]. were applied to the top of the baseplate throughout level gait in 2% intervals (Fig.1a). For varus alignment, the varus/valgus moment was modified to match the ratio of medial-lateral force distribution from Smith et al. [1]. (Fig.1b). Results. For both alignments and all bones, the largest micromotion and amount of bone at risk of failure occurred during mid stance, at 16% of gait (Figs.2,3). Peak micromotion, located at the antero-lateral edge of the baseplate, was 153±32 µm and 273±48 µm for mechanical and varus alignment, respectively. The area of the baseplate with micromotion above 40 µm (the threshold for bone ingrowth. [3]. ) was 28±5% and 41±4% for mechanical and varus alignment, respectively. The amount of bone at risk of failure at the bone-implant interface was 0.5±0.3% and 0.8±0.3% for the mechanical and varus alignment, respectively. Discussion. The peak micromotion and the baseplate area with micromotion above 40 µm increased with varus alignment compared to mechanical alignment. Furthermore, the amount of bone at risk of failure, although small for both alignments, was greater for varus alignment. These results suggest that varus alignment, consisting of a combination of femoral and tibial alignment, may negatively impact bone ingrowth and increase the risk of bone failure for cementless tibial baseplates of this TKR design