During computer navigated total knee arthroplasty, pin holes are drilled in the femur and tibia to allow the placement of navigation trackers, and fractures associated with these pin holes have recently been reported. We hypothesized that an increase in stress around the pin holes is one of the most relevant factors contributing to the fracture. In this study, we used finite element analysis to investigate the stresses around femoral pin holes with respect to the mode of pin penetration, the diameter of the pin holes, and the degree of osteoporosis. Our results indicate that increases in pin hole diameter and reduction in bone strength as a result of osteoporosis intensify the stresses around the pin holes, especially in cases of transcortical pin penetration.
1 January 2012 in Journal of Biomechanical Engineering
Determination of physiological loads in human lumbar spine is critical for understanding the mechanisms of lumbar diseases and for designing surgical treatments. Computational models have been used widely to estimate the physiological loads of the spine during simulated functional activities. However, various assumptions on physiological factors such as the intra-abdominal pressure (IAP), centers of mass (COMs) of the upper body and lumbar segments, and vertebral centers of rotation (CORs) have been made in modeling techniques. Systematic knowledge of how these assumptions will affect the predicted spinal biomechanics is important for improving the simulation accuracy. In this paper, we developed a 3D subject-specific numerical model of the lumbosacral spine including T12 and 90 muscles. The effects of the IAP magnitude and COMs locations on the COR of each motion segment and on the joint/muscle forces were investigated using a global convergence optimization procedure when the subject was in a weight bearing standing position. The data indicated that the line connecting the CORs showed a smaller curvature than the lordosis of the lumbar spine in standing posture when the IAP was 0 kPa and the COMs were 10 mm anterior to the geometric center of the T12 vertebra. Increasing the IAP from 0 kPa to 10 kPa shifted the location of CORs toward the posterior direction (from 1.4 ± 8.9 mm anterior to intervertebral disc (IVD) centers to 40.5 ± 3.1 mm posterior to the IVD centers) and reduced the average joint force (from 0.78 ± 0.11 Body weight (BW) to 0.31 ± 0.07 BW) and overall muscle force (from 349.3 ± 57.7 N to 221.5 ± 84.2 N). Anterior movement of the COMs from −30 mm to 70 mm relative to the geometric center of T12 vertebra caused an anterior shift of the CORs (from 25.1 ± 8.3 mm posterior to IVD centers to 7.8 ± 6.2 mm anterior to IVD centers) and increases of average joint forces (from 0.78 ± 0.1 BW to 0.93 ± 0.1 BW) and muscle force (from 348.9 ± 47.7 N to 452.9 ± 58.6 N). Therefore, it is important to consider the IAP and correct COMs in order to accurately simulate human spine biomechanics. The method and results of this study could be useful for designing prevention strategies of spinal injuries and recurrences, and for enhancing rehabilitation efficiency.
Although the clinical importance of extragraft bone formation (ExGBF) and bridging (ExGBB) has been reported, few studies have investigated the biomechanical influences of ExGBF on the motion segment. In this study, ExGBF was simulated at the C5-C6 motion segment after anterior cervical discectomy and fusion using a developed finite element model and a sequential bone-remodelling algorithm in flexion and extension. The computer simulation results showed that extragraft bone was primarily formed in the extension motion and grew to form ExGBB. A stepwise decrease in the intersegmental rotation angle, maximum von Mises stress and strain energy density on the trabecular bone with ExGBF were predicted in extension. When ExGBB was formed in the trabecular bone region, the intersegmental rotation angle slightly decreased with additional bone formation. However, the stress and strain energy density on the trabecular bone region decreased until ExGBB reached the peripheral cortical margin. The results offer a rationale supporting the hypothesis that mechanical stimuli influence ExGBF. ExGBF was helpful in increasing the stability of the motion segment and decreasing the fracture risk of trabecular bones, even in cases in which ExGBB was not formed. ExGBB can be classified as either soft or hard bridging based on a biomechanical point of view.
Mechanical seals, such as labyrinth seals, are typically installed at the turbine outlets to prevent oil leakage. However, these seals undergo deformation because of the vibrations of the rotor, even during normal turbine operating conditions, which may cause an increase in oil leakage. In this study, the oil leakage performance of three labyrinth seals with different types of seal teeth, narrow stainless teeth (Type 1), wide aluminum teeth fixed on the body (Type 2), and fixed wide aluminum movable teeth (Type 3), were evaluated using finite element (FE) and computational fluid dynamics (CFD) analyses. Three-dimensional FE models of the rotor and oil deflectors were developed, and the plastic deformation of the teeth of the labyrinth seals was predicted when the rotor impacted the sealing teeth during turbine operation. The oil leakage was predicted using CFD analysis. The results indicated that the Type 3 seal, including movable teeth, is beneficial in preventing leakage and tooth deformation compared with the other types. The Type 2 seal is advantageous because it results in a smaller increase in gap size and greater vena contracta effects than the Type 1 seal. The results of this study could be helpful when designing and selecting the teeth of a labyrinth seal.
The core inlet flow distribution in the APR1000 reactor is critical for ensuring the reactors safety and efficient operation by maintaining uniform coolant flow across fuel assemblies. Previous studies, though insightful, faced challenges in fully replicating reactor-scale flow conditions due to technical and economic constraints associated with scaled-down experimental models and the limited numerical validation methodologies. This study addresses these limitations by developing and validating a robust computational fluid dynamics (CFD) methodology to accurately analyze the core inlet flow distribution. A 1/5 scaled-down experimental model adhering to similarity laws was employed for validation. CFD analyses using ANSYS Fluent and CFX, combined with turbulence model evaluations and grid sensitivity studies, demonstrated that the SST and RNG k-ε turbulence models provided the most accurate predictions, with a high correlation to previous experimental data. Full-scale simulations revealed uniform coolant distribution at the core inlet, with peripheral assemblies exhibiting higher flow rates, consistent with previous experimental observations. Quantitative metrics such as the coefficient of variation (COV), relative error (RD), and root mean square error (RMSE) confirmed the superior performance of the SST model in CFX, achieving a COV of 7.993% (experimental COV: 5.694%) and an RD of 0.047. This methodology not only validates the CFD approach but also highlights its applicability to reactor design optimization and safety assessment. The findings of this study provide critical guidelines for analyzing complex thermal-fluid systems in nuclear reactor designs.
25 March 2021 in International Journal of Precision Engineering and Manufacturing