News Detail
Stress Analysis of Tibial Bone Using Three Different Materials for Bone Fixation PlatesAbstractStress shielding is a problem for traditional metal bone fixation plates made of magnesium and titanium alloys. This problem can be solved by using composite materials with a low elastic modulus. This study analyzed the effect of carbon fiber reinforced PEEK (CFRP) composites on stress shielding under static loading using finite element simulations. Callus formation times relative to the healing period were gradually imposed according to the elapsed time, considering 1% and 75% as healing stages. The Inventor? 3D CAD 2024 software was used for modeling, and the ANSYS? FEA R2023 software was used for analysis. The results showed that metal fixation plates made of titanium and magnesium alloys transferred less stress to the bone than the CFRP fixation plate. In particular, the use of the CFRP fixation plate resulted in a higher peak stress and a more uniform stress field in the bone, especially in the bone-plate contact area, where the risk of stress shielding is higher in the 1% and 75% healing phases. Keywords: human tibia fractures; fixation plates; metallic alloys; CFRP; stress shielding; finite element analysis1. IntroductionHuman tibia fractures are common bone injuries due to falls or traffic accidents [1,2,3]. Depending on the fracture type and severity, different healing methods are used [4]. For high energy fractures involving significant bone debris, indirect healing (i.e., plastering) is often adopted. For simpler (low energy) fractures of long bones, direct healing, involving the direct application of internal fixation devices to the broken bones, is usually preferred. The healing of bone fractures may be primary or secondary [5]. Primary healing requires the precise anatomical realignment of fractured extremities without gaps, supporting direct and stable contact. For unstable and indirect contact, fractured bones heal by secondary healing, combining intramembranous and endochondral ossification [5]. Mechanical stimuli, such as stress and deformation, play a crucial role in secondary healing. Indeed, mechanical stress applied to fractured bone stimulates bone cells to produce new bone tissue, promoting bone callus formation and accelerating the healing process [6,7,8]. In addition, applying controlled stress levels to bone during remodeling contributes to increased bone density and the formation of more resistant bone tissue. Without sufficient surrounding stress, bone becomes fragile and unstable and may prematurely re-fracture after removal of the bone plate [8]. Bone plates are an important type of internal fixation device for stabilizing human bone fractures and are typically made of stainless steel and titanium alloys. Although metals have high strength, ductility, and wear resistance, they may undergo corrosion in physiological conditions, increased thermal conductivity, as well as be sensitive to electro-magnetic fields [9,10,11,12]. Ti-6Al-4V has a Young’s modulus five times higher than human cortical bone [13]. The significantly higher stiffness of metallic plates for bone fixation with respect to cortical bones leads to the increased transfer of higher loads through plates rather than through bones. This phenomenon, referred to as the “stress shielding effect” or “stress protection”, retards the natural bone healing process [14]. However, the presence of these metal plates in the human body for a long period of time may cause pain, infection, tissue irritation, etc. Therefore, they must be removed from the body after a certain period of time through secondary surgery. Other metal alloys, such as magnesium alloys, can naturally degrade in the human body. For example, the ZM21 alloy easily decomposes over time, thus eliminating the need for plate removal via secondary surgery. In the presence of chloride and at neutral pH, as in bodily fluids, Mg hydroxide can be transformed into the completely soluble Mg chloride (Equation (1)): Magnesium is an essential element for the human body, acting as co-factor for many enzymes. However, a high concentration of magnesium (>1.05 mmol/L) in bodily fluid may not be tolerated. Witte et al. [15] and Yamasaki et al. [16] have shown that magnesium plays a key role in bone remodeling. In addition, the mechanical properties of magnesium are closer to that of natural bone than other metal alloys and biodegradable polymers, which may reduce the stress shielding effect commonly associated with metallic materials. The main problem with using magnesium biodegradable implants is that degradation rate of pure Mg in physiological conditions is very high. Such a high degradation rate may cause the structural integrity of Mg implants to vanish before completing the bone healing process. Recently, several magnesium alloys have been investigated for use in biodegradable implants [17,18,19]. The AZ series alloys (aluminum and zinc) are the most investigated, mainly because of their commercial availability, showing less degradation than pure magnesium. The addition of calcium to magnesium alloys has improved degradation behavior. In addition, various coating technologies, such as polymeric, ceramic, and even metallic coating, have been developed for magnesium alloys in order to reduce initial degradation [20]. Degradation tests have shown, for example, that hydroxyapatite particles stabilize the degradation rate and present more uniform degradation in artificial seawater and cell solutions [21]. Another way to minimize stress shielding is to change the geometry of the bone fixation plate. The equivalent stiffness of metal parts and the material consumption in manufacturing processes may be adjusted by integrating topology optimization and additive manufacturing [22]. An alternative solution to stress shielding is to find innovative materials for fixation plates: the basic requirement is that the elastic modulus of fixation materials should be as close as possible to the bone modulus [22]. In the last decade, composite materials have had a significant impact in the field of prosthetic implants, with a particular focus on bone fixation plates [23,24,25]. Recent studies [26,27] have been actively exploring the use of composites in bone fixation plates due to their ability to improve properties such as strength, stiffness, and corrosion resistance. In clinical applications, carbon fiber reinforced PEEK (CF/PEEK) plates are radiolucent and do not obstruct radiographic analysis; hence, their CT images have much less artifacts than those recorded for metal alloy fixation plates [28]. In addition, as shown above, metal ions released from metal implants can cause allergic reactions in tissues, as the body may react with hypersensitivity to such materials. Replacing metal with CF/PEEK can prevent this potential hypersensitivity, because this composite does not release metal ions, thus reducing the risk of allergic reactions in surrounding tissues [29]. For the purpose of developing high-performance composite bone plates with mechanical properties that promote bone fracture healing, several research studies have been conducted. Fujihara et al. [30] investigated the manufacturing process of carbon/PEEK thermoplastic bone plates; mechanical performance was partially verified by comparing experimental data gathered from 4-point bending tests with analytical results. Kim et al. [31] compared the mechanical performance of a composite bone plate with a commercially available stainless steel bone plate then determined the best conditions for the composite plate. In evaluating the best performance of composite bone plates, the callus stress distributions generated between the fracture surfaces were considered an important criterion. Following this rationale, the present study utilized the finite element method (FEM) to analyze and evaluate the effect of stress shielding, in terms of stress transmission to bone, on a tibial fracture fixed by a titanium alloy plate, a magnesium alloy plate, and a composite material (carbon fiber reinforced PEEK, CFRP) plate. The manuscript is organized as follows: Section 2 describes the geometric modeling and the finite element modeling/analysis of the fixation plate. Section 3 presents the results of FE simulations relative to the different materials of the fixation plates. Section 4 discusses the results while in view of the limitations of the study. Section 5 summarizes the main findings of the study and outlines future research directions. 2. Materials and Methods2.1. Model GeometryThis study focuses on the fixation of tibial fractures with composite material plates. The tibia, one of the main load-bearing bones of the human leg, is shown in Figure 1. It transfers loads between the foot and femur during standing or walking. Similar to other bones, the outer part of bone is made up of cortical bone, while the core is made up of trabecular bone. Cortical bone is composed of hard, dense tissue and is involved in load transfer, such as flexion and compression. In contrast, the trabecular part is composed of spongy tissue to reduce structural mass [31,32]. Figure 1. 3D view of the tibia bone and details of the bone cross-section highlighting the cortical and trabecular regions. To simplify the geometry, the tibia was modeled as a hollow cylinder. The geometric specifications of the tibia (averaged from anatomical databases) were 10 mm inner diameter, 25 mm outer diameter, 103 mm total length (see Figure 2) [32]. This model was used to perform finite element simulations for evaluating the structural behavior of the tibia and the fixation devices used during the healing process. The fracture site was assumed to have 1 mm thickness [32], with the outer and inner diameters corresponding to those of the tibia. The plate was modeled with six holes that crossed both the plate itself and the thickness of the bone. The plate dimensions were defined as length 68 mm, width 10 mm and thickness 3 mm. All three-dimensional models were created using Autodesk Inventor 2024 software. 2.2. Material PropertiesBones are typically composed of cortical bone in the outer region and trabecular bone in the inner region [32], which have an anisotropic behavior. An anisotropic material is described by 21 elastic constants, but bone has been shown to exhibit transverse isotropic behavior in the YZ and XZ planes. Therefore, it can be described as an orthotropic model with 9 elastic constants (Ex, Ey, Ez, Gxy, Gyz, Gxz). Two stages of healing in the fractured debris (1% and 75%) were considered for analysis. It was assumed that the callus acts as a bridge at the fracture site, and that the callus material is homogeneous and isotropic [32]. Regarding bone plate materials, the following fixation plate materials were considered in this study.
On the other hand, CFRP mechanical behavior depends on the fibers’ mechanical properties as well as on the fibers’ orientation in the different layers. In this study, configuration I was selected, with fibers oriented at 0°, +45°, ?45°, and 90°, where the 0° orientation corresponds to the longitudinal direction of the tibia. This configuration allows the optimal distribution of the applied load and results in a homogeneous mechanical behavior of the system [33,34,35,36,37]. The selected orientations of fibers help optimally distribute and transmit the applied load to the bone. In addition, mechanical behavior is transversely isotropic in this configuration. An approximately 60% fiber volume fraction was selected, as this is the critical volume below which mechanical properties of composite material are inferior to those of the matrix [38]. The mechanical properties of all materials used in the FEA of fixed bone are listed in Table 1 [38]. 2.3. Loading and Kinematic Boundary ConditionsThe osteogenic ability of bone is influenced by the level of principal stress generated by applied loading modes. Stress varies if the bone is subjected to compression, flexion, and torsion. It has been observed that the maximum principal stress under compression and flexion can favor osteogenesis more than torsion [32,38,39,40,41]. Therefore, the bone and fixation plate assembly were modeled here subject to compression and flexion, excluding the torsion case, as shown in Figure 3. Figure 3. Loading and kinematic constraint conditions: The bone’s left extremity was fixed, while a pressure of 1.1 MPa and a bending moment of 1 Nmm was applied on the bone’s right extremity. For a 90 kg patient standing on both legs, the magnitude of the load acting on the tibia cross-sectional area is approximately (Equation (2)) 2.4. Finite Element ModelThe numerical simulation of the mechanical behavior of the plate-fractured bone assembly was performed with the commercial finite element code ANSYS Workbench? Version R2023. The geometric models were imported from Inventor 2024 and a mesh was then generated for all parts considered in the study. The SOLID 226 hexahedral element was selected for the static analysis [22,32,41,42]. To ensure the accuracy of results, convergence analysis was carried out to find the optimal mesh size. A coarser element size was initially selected and then gradually reduced until the average variation of nodal displacements between two successive iterations was less than 0.01%. The final mesh size was set to 1.5 mm for bone, fixation plate, and screws and 0.5 mm for the fractured bone area [32]. In total, the FE model included 336,257 nodes and 196,421 elements (see Figure 5). To calculate behavior at the interfaces, two contact surfaces were defined between (i) the bone fixation plate and the tibial bone, (ii) screws and bone, and (iii) screws and bone fixation plate. A coefficient of friction of 0.4 was chosen between the plate and the bone [31]. Following Perren et al. [43], the appropriate traction force to stabilize fractures with the fixation plate in place ranges between 2000 and 3000 N. Higher screw fixation forces lead to excessive interface stress, causing blood flow problems in contact areas, while lower forces lead to excessive movement between bone plate and bone, resulting in ineffective bone stabilization. Therefore, a preload of 2000 N was applied to the screw. All screws were considered to be fixed in the bone. 3. ResultsThe von Mises equivalent stress distributions developed in the healing phases at 1% and 75% of the fractured area are shown in Figure 6, highlighting the influence of various bone fixation plate materials. Figure 6. Distribution of Von Mises stress in fractured bone cross-sectional areas for different plate materials: 1% and 75% healing level. The von Mises stress gradually increases from the fixation plate side to the opposite side of bone. As expected, the minimum equivalent stress is localized on the fixation plate side that carries most of the load. In the beginning to 1% stage of fracture healing (approximately one week after surgery), there is no substantial hard tissue in the fracture gap. Therefore, the fracture site is characterized primarily by the presence of hematoma and fluid. These soft tissues are unable to withstand mechanical loads due to their relatively low Young’s modulus. Consequently, stresses in the fractured area are relatively low for all bone plate materials, and the stress shielding caused by the fixation plate is at its highest level at the 1% healing stage compared to the 75% healing stage. Figure 6 shows how, at the 75% healing stage, more callus is formed: the callus becomes hard enough to support additional external loads. In fact, the stress distribution at this stage tends to be higher and more uniform in the bone. To analyze the stress distribution in the fractured area in detail, three control points were defined, as shown in Figure 7. Specifically, point A is localized right below the fixation plate; point B is localized in the middle diameter of the bone; and point C is localized in the distal part of point A. Figure 7. Positions of the three control points chosen for analyzing stress distribution in fractured bone. Figure 8 shows the values of the von Mises stress at the control points A, B, and C, considering a healing phase of 1%. Figure 8. Von Mises stress values in fractured bone, considering different materials for the bone fixation plate and a healing stage of 1%. For all fixation plates, the minimum stress occurs at point A, i.e., at the contact between the plate and the bone, with 1.12 MPa for the titanium alloy Ti-6Al-4V, 4.21 MPa for magnesium alloy ZM 21, and 7.35 MPa for the CFRP composite. It can also be observed that the stress increases as the control points move away from the plate. This can be explained because, looking at Figure 4, the stress changes linearly across the fractured bone section, with a point of minimum at the plate–bone contact, due to the phenomenon of stress shielding. In the 1% healing phase, the Ti-6Al-4V plate has a lower average stress in the bone than the other two plates. Since bone healing begins at the 1% healing stage, the generated tissues are not strong enough, and stress shielding is more evident. Figure 9 shows the values of von Mises stress developed in the fractured bone at the 75% healing stage. Von Mises stress values at all control points of the fractured bone are higher than in the healing phase 1% for the titanium alloy and composite material. This can be explained with the argument that bone gains strength during healing due to external mechanical stimulation and can therefore absorb more stress and become more resistant. Using the CFRP fixation plate results in higher stresses in the bone than in the case of titanium and magnesium alloy plates (see Table 2). Figure 9. Von Mises stress values in fractured bone, considering different materials for the bone fixation plate and a healing stage of 75%. Table 2. Von Mises stress values for different materials at the position (A, B, C) of the tibia, for a healing phase of 75%. Figure 6 also shows that the stress distribution in the bone is more uniform with the CFRP tread at 75% healing. Thus, the effect of stress shielding is less significant for CFRP than for titanium and magnesium alloys. This can be explained by the fact that the stiffness of CFRP is closer to that of bone than that of the two metallic materials (see Table 1. From the point of view of mechanical behavior, stress analysis within the bone plate, which serves as a structural component to support the fractured area, is of considerable importance. In the 1% healing phase, the fractured area is weak and has the least capacity to support external loads. Most of the load is carried by the bone fixation plate, resulting in the highest stress within the plate at this stage. Figure 10 shows stress distributions in fixation plates in the 1% healing phase. Figure 10. Von Mises stress distributions in the three bone fixation plates at the 1% healing stage. It can be seen that stress concentration is highest in the two central screw holes near fracture site. This is because the fractured area is too weak to support external stress, and most stress is retained from the bone plate and screw, concentrating at the edge of the screw hole. From the cross-sectional view, it is possible to analyze the stress through the cross-section. Figure 10 shows that the most stressed plate is Ti-6Al-4V with 32 MPa of stress, followed by magnesium alloy ZM 21 with 25 MPa, and CFRP with 18.76 MPa. This shows that a stiffer material such as titanium absorbs most of the stress by shielding the bone from stress, therefore the stress shielding phenomenon is less evident when materials with mechanical properties closer to that of bone are used. 4. DiscussionBone plates are the most used implants for the internal fixation of bone fractures. They present numerous advantages, including adequate stability and resistance to tension, compression, creep, torsion, and bending [44,45]. These are important characteristics for plates. They provide the support and strength necessary for proper fracture healing. In terms of materials, metal alloys such as 316 L stainless steel have been used in the manufacture of bone plates since 1895 [46]. However, due to high corrosion rates in human physiological environments and poor resistance to stress and wear, they were eliminated. Subsequently, titanium alloys such as Ti-6Al-4V were selected as candidates to replace steels. These alloys exhibit lower modulus, higher specific strength, superior biocompatibility, and improved corrosion resistance compared to stainless steels [47]. Resistance to corrosion from body fluids is provided by the presence of a passive TiO2 film on the surface, which has the property of immediately reforming at body temperature and in physiological fluids, ensuring the excellent corrosion resistance of titanium alloy implants in vivo [48]. However, these materials present complex problems related to wear resistance and fatigue strength. These problems can lead to potential accidents such, as fractures and plate loosening [49]. The large difference in stiffness between metallic materials and bone causes stress shielding. This phenomenon occurs when a rigid orthopedic implant absorbs most of the load, preventing the surrounding bone tissue from performing its natural role of support and loading [50]. This can lead to problems such as reduced bone density and the risk of secondary fractures [51]. To overcome this problem, a material with similar strength and stiffness to bone must be chosen. Polymeric materials, such as PEEK, have been considered as an alternative to metallic materials for orthopedic implants [52]. Polymers are lightweight, easy to form, and can have a Young’s modulus closer to that of bone tissue than metallic materials [52]. The inferior mechanical properties of polymeric fixation plates, compared to conventional titanium alloy plates, have been indicated as the main limitation of polymeric plates in several studies [52,53,54,55]. Particularly in skeletal sites subject to loading, lower mechanical strength may have a negative impact on fracture stability. In fact, in a study by Schliemann et al. [56], the lower stiffness of the PEEK plate, together with the increased flexion and movement between bone and implant, contributed to increased movement in the fracture zone, resulting in a reduced load bearing capacity and, consequently, an earlier failure of fracture fixation. Recently, magnesium (Mg) and its alloys proved to be promising biomaterials for bone plates due to their excellent mechanical properties, biocompatibility, and degradability in physiological fluids [57]. The Young’s modulus of Mg is similar to that of cortical bone, which reduces stress shielding and promotes fracture healing [58]. Several studies have investigated Mg-based metallic materials for internal fracture fixation, with research focusing on reducing degradation rates and increasing corrosion resistance through methods such as CaP coating [59,60,61,62]. Looking at bone tissue, it is mainly composed of organic materials, such as collagen, and inorganic materials, such as hydroxyapatite. The organic materials are responsible for flexibility and strength, while the inorganic materials are responsible for hardness and compressive strength [63]. Composite materials inspired by this natural composition may achieve optimal mechanical properties [64,65,66]. In a study conducted by Sentihil Maharaj et al. [67], femoral bone fixation plates were developed using carbon fabric interwoven with flax fabric and bio-epoxy resin. Such a composite material was shown to have good strength and stiffness, like that of femur bone. This minimizes stress shielding. For example, the tensile strengths of stainless steel (515 MPa) and titanium (950 MPa) are much higher than that of femoral bone (139 MPa). A major advantage of polymer matrix composites is the possibility of tailoring mechanical performance based on several factors, including matrix type, reinforcement type, reinforcing phase particle size distribution, and reinforcing phase quantity. This flexibility allows composite properties to be tailored to the specific application [68]. Furthermore, Petersen et al. [69] found that carbon fiber-reinforced composites stimulate osseointegration more effectively than titanium, suggesting a potential improvement in biocompatibility and bone integration of implants. Furthermore, carbon fibers have conductive properties that may positively influence cellular activity. It is hypothesized that they may act as conductive micro-biocircuits, helping to remove damaging free radicals and improve the environment for cell growth. The finite element method (FEM) can be used to study the stress distribution in bone tissues. The aim is to determine the possible causes of the stress shielding problem and how it can be reduced by modifying the material and geometry of the prosthetic element. For example, Ceddia et al. [22] carried out topology optimization of a femoral prosthesis made of CF/PEEK composite. Remarkably, the total mass of the prosthesis was reduced by 30% while maximum stress decreased from 987 MPa to about 810 MPa. However, the results shown in this study confirmed the theory that a material of similar stiffness to bone can equalize the stress distribution in bone by reducing the phenomenon of stress shielding. This explains why the magnesium alloy fixation plate reached its stress peak at the 1% healing phase. Conversely, the composite CFRP fixation showed more uniform and consistent stress distributions in the bone at the 1% and 75% healing stages. Such a behavior makes CFRP an excellent first choice material to use for bone fixating plates. Since this study included limitations such as the simplified numerical and geometric modeling of the tibia and the use of a limited set of loading conditions, FE results should be compared with experimental data gathered from in vitro/vivo tests to validate the numerical model proposed here. In addition, the implementation of dynamic loads could help to better understand the fatigue life of the materials studied in this work, providing clinicians and designers of orthopedic prostheses with optimized selection in reference to the patient’s characteristics. 5. ConclusionsThe study used finite element simulations to analyze the influence of different types of materials: titanium alloy Ti-6Al-4V, magnesium alloy ZM21, and the composite material CFRP. Simulation results carried out on a simplified but realistic model showed that the CFRP bone fixation plate had a higher bone loading capacity than both metallic plates. This was evident not only from stress values at the fracture site near the plate but also from the average bone stress values throughout the fracture site, which were highest for the CFRP fixation. Therefore, in view of the required physical and chemical properties for fracture healing, this study indicates the CFRP fixation plate as the most suitable material for the fracture healing process. In fact, CFRP allows a better and more uniform stress distribution to be achieved in the fracture area. The stress shielding phenomenon is more pronounced in the post-operative phase, when the bone is in the 1% healing phase, compared to the metal bone fixation plates currently used in biomedical and healthcare applications. Author ContributionsConceptualization, M.C. and B.T.; methodology, M.C.; software, M.C.; validation, B.T. and G.S.; formal analysis, M.C., L.L. and B.T.; investigation, M.C.; resources, M.C. and B.T.; data curation, M.C.; writing—original draft preparation, M.C. and B.T.; writing—review and editing, M.C., B.T., L.L. and M.T.; visualization, B.T. and L.L.; supervision, G.S. and B.T.; project administration, G.S. and B.T. All authors have read and agreed to the published version of the manuscript. FundingThis research received no external funding. Institutional Review Board StatementNot applicable. Informed Consent StatementNot applicable. Data Availability StatementAll experimental data to support the findings of this study are available by contacting the corresponding author upon request. Conflicts of InterestThe authors declare no conflicts of interest. References
|