Objective To design customized titanium alloy lunate prosthesis, construct three-dimensional finite element model of wrist joint before and after replacement by finite element analysis, and observe the biomechanical changes of wrist joint after replacement, providing biomechanical basis for clinical application of prosthesis. Methods One fresh frozen human forearm was collected, and the maximum range of motions in flexion, extension, ulnar deviation, and radialis deviation tested by cortex motion capture system were 48.42°, 38.04°, 35.68°, and 26.41°, respectively. The wrist joint data was obtained by CT scan and imported into Mimics21.0 software and Magics21.0 software to construct a wrist joint three-dimensional model and design customized titanium alloy lunate prosthesis. Then Geomagic Studio 2017 software and Solidworks 2017 software were used to construct the three-dimensional finite element models of a normal wrist joint (normal model) and a wrist joint with lunate prosthesis after replacement (replacement model). The stress distribution and deformation of the wrist joint before and after replacement were analyzed for flexion at and 15°, 30°, 48.42°, extension at 15°, 30°, and 38.04°, ulnar deviation at 10°, 20°, and 35.68°, and radial deviation at 5°, 15°, and 26.41° by the ANSYS 17.0 finite element analysis software. And the stress distribution of lunate bone and lunate prosthesis were also observed. Results The three-dimensional finite element models of wrist joint before and after replacement were successfully constructed. At different range of motion of flexion, extension, ulnar deviation, and radial deviation, there were some differences in the number of nodes and units in the grid models. In the four directions of flexion, extension, ulnar deviation, and radial deviation, the maximum deformation of wrist joint in normal model and replacement model occurred in the radial side, and the values increased gradually with the increase of the range of motion. The maximum stress of the wrist joint increased gradually with the increase of the range of motion, and at maximum range of motion, the stress was concentrated on the proximal radius, showing an overall trend of moving from the radial wrist to the proximal radius. The maximum stress of normal lunate bone increased gradually with the increase of range of motion in different directions, and the stress position also changed. The maximum stress of lunate prosthesis was concentrated on the ulnar side of the prosthesis, which increased gradually with the increase of the range of motion in flexion, and decreased gradually with the increase of the range of motion in extension, ulnar deviation, and radialis deviation. The stress on prosthesis increased significantly when compared with that on normal lunate bone. Conclusion The customized titanium alloy lunate prosthesis does not change the wrist joint load transfer mode, which provided data support for the clinical application of the prosthesis.
ObjectiveTo explore the biomechanical stability of the medial column reconstructed with the exo-cortical placement of humeral calcar screw by three-dimensional finite element analysis. MethodsA 70-year-old female volunteer was selected for CT scan of the proximal humerus, and a wedge osteotomy was performed 5 mm medially inferior to the humeral head to form a three-dimensional finite element model of a 5 mm defect in the medial cortex. Then, the proximal humeral locking plate (PHILOS) was placed. According to distribution of 2 calcar screws, the study were divided into 3 groups: group A, in which 2 calcar screws were inserted into the lower quadrant of the humeral head in the normal direction for supporting the humeral head; group B, in which 1 calcar screw was inserted outside the cortex below the humeral head, and the other was inserted into the humeral head in the normal direction; group C, in which 2 calcar screws were inserted outside the cortex below the humeral head. The models were loaded with axial, shear, and rotational loadings, and the biomechanical stability of the 3 groups was compared by evaluating the peak von mises stress (PVMS) of the proximal humerus and the internal fixator, proximal humeral displacement, neck-shaft angle changes, and the rotational stability of the proximal humerus. Seven cases of proximal humeral fractures with comminuted medial cortex were retrospectively analyzed between January 2017 and December 2020. Locking proximal humeral plate surgery was performed, and one (5 cases) or two (2 cases) calcar screws were inserted into the inferior cortex of the humeral head during the operation, and the effectiveness was observed. Results Under axial and shear force, the PVMS of the proximal humerus in group B and group C was greater than that in group A, the PVMS of the internal fixator in group B and group C was less than that in group A, while the PVMS of the proximal humerus and internal fixator between group B and group C were similar. The displacement of the proximal humerus and the neck-shaft angle change among the 3 groups were similar under axial and shear force, respectively. Under the rotational torque, compared with group A, the rotation angle of humerus in group B and group C increased slightly, and the rotation stability decreased slightly. All the 7 patients were followed up 6-12 months. All the fractures healed, and the healing time was 8-14 weeks, with an average of 10.9 weeks; the neck-shaft angle changes (the difference between the last follow-up and the immediate postoperative neck-shaft angle) was (1.30±0.42)°, and the Constant score of shoulder joint function was 87.4±4.2; there was no complication such as humeral head varus collapse and screw penetrating the articular surface. ConclusionFor proximal humeral fractures with comminuted medial cortex, exo-cortical placement of 1 or 2 humeral calcar screw of the locking plate outside the inferior cortex of the humeral head can also effectively reconstruct medial column stability, providing an alternative approach for clinical practice.
Objective To establish the finite element model of Y-shaped patellar fracture fixed with titanium-alloy petal-shaped poly-axial locking plate and to implement the finite element mechanical analysis. Methods The three-dimensional model was created by software Mimics 19.0, Rhino 5.0, and 3-Matic 11.0. The finite element analysis was implemented by ANSYS Workbench 16.0 to calculate the Von-Mises stress and displacement. Before calculated, the upper and lower poles of the patella were constrained. The 2.0, 3.5, and 4.4 MPa compressive stresses were applied to the 1/3 patellofemoral joint surface of the lower, middle, and upper part of the patella respectively, and to simulated the force upon patella when knee flexion of 20, 45, and 90°. Results The number of nodes and elements of the finite element model obtained was 456 839 and 245 449, respectively. The max value of Von-Mises stress of all the three conditions simulated was 151.48 MPa under condition simulating the knee flexion of 90°, which was lower than the yield strength value of the titanium-alloy and patella. The max total displacement value was 0.092 8 mm under condition simulating knee flexion of 45°, which was acceptable according to clinical criterion. The stress concentrated around the non-vertical fracture line and near the area where the screws were sparse. Conclusion The titanium-alloy petal-shaped poly-axial locking plate have enough biomechanical stiffness to fix the Y-shaped patellar fracture, but the result need to be proved in future.
Objective The biomechanical characteristics of three internal fixation modes for femoral subtrochanteric spiral fracture in osteoporotic patients were compared and analyzed by finite element technology, so as to provide the basis for the optimization of fixation methods for femoral subtrochanteric spiral fracture. MethodsTen female patients with osteoporosis and femoral subtrochanteric spiral fractures caused by trauma, aged 65-75 years old, with a height of 160-170 cm and a body weight mass of 60-70 kg, were selected as the study subjects. The femur was scanned by spiral CT and a three-dimensional model of the femur was established by digital technology. The computer aided design models of proximal intramedullary nail (PFN), proximal femoral locking plate (PFLP), and the combination of the two (PFLP+PFN) were constructed under the condition of subtrochanteric fracture. Then the same load of 500 N was applied to the femoral head, and the stress distribution of the internal fixators, the stress distribution of the femur, and the displacement of femur after fracture fixation were compared and analyzed under the three finite element internal fixation modes, so as to evaluate the fixation effect. ResultsIn the PFLP fixation mode, the stress of the plate was mainly concentrated in the main screw channel, the stresses of the different part of the plate were not equal, and gradually decreased from the head to the tail. In the PFN fixation mode, the stress was concentrated in the upper part of the lateral middle segment. In the PFLP+PFN fixation mode, the maximum stress appeared between the first and the second screws in the lower segment, and the maximum stress appeared in the lateral part of the middle segment of the PFN. The maximum stress of PFLP+PFN fixation mode was significantly higher than that of PFLP fixation mode, but significantly lower than that of PFN fixation mode (P<0.05). In PFLP and PFN fixation modes, the maximum stress of femur appeared in the medial and lateral cortical bone of the middle femur and the lower side of the lowest screw. In PFLP+PFN fixation mode, the stress of femur concentrated in the medial and lateral of the middle femur. There was no significant difference in the maximum stress of femur among the three finite element fixation modes (P>0.05). The maximum displacement occurred at the femoral head after three finite element fixation modes were used to fix subtrochanteric femoral fractures. The maximum displacement of femur in PFLP fixation mode was the largest, followed by PFN, and PFLP+PFN was the minimum, with significant differences (P<0.05). ConclusionUnder static loading conditions, the PFLP+PFN fixation mode produces the smallest maximum displacement when compared with the single PFN and PFLP fixation modes, but its maximum plate stress is greater than the single PFN and PFLP fixation mode, suggesting that the combination mode has higher stability, but the plate load is greater, and the possibility of fixation failure is higher.
This study was aimed to compare the mechanical characteristics under different physiological load conditions with three-dimensional finite element model of rigid fixation and elastic fixation in the lumbar. We observed the stress distribution characteristics of a sample of healthy male volunteer modeling under vertical, flexion and extension torque situation. The outcomes showed that there existed 4-6 times pressure on the connecting rod of rigid fixation compared with the elastic fixations under different loads, and the stress peak and area of force on elastic fixation were much higher than that of the rigid fixations. The elastic fixation has more biomechanical advantages than rigid fixation in promoting interbody lumbar fusion after surgery.
ObjectiveTo investigate biomechanical effects of pseudo-patella baja on stress of patellofemoral joint after total knee arthroplasty (TKA) by using finite element analysis (FEA).MethodsA series of CT and MRI of the left knee joint of two healthy volunteers and three-dimensional (3D) scanned data of TKA prosthesis were taken, and the 3D models of knee before and after TKA were established. The finite element model of pseudo-patella baja, normal patella, and alta patella after TKA were constructed by Insall-Salvafi (IS) ratio and Blackburne-Peel (BP) ratio. The load was applied along the direction of quadriceps femoris. After testing the validity of the finite element model, the high contact stress of patellofemoral joint was measured on the von Mise stress nephogram of pseudo-patella baja, normal patella, and alta patella after TKA when the knee flexion was 30°, 60°, and 90°. The average contact area was calculated according to two volunteers’ data.ResultsOn the finite element model of the normal patella after TKA with knee flexion 30°, 475 N pressure was applied along the direction of quadriceps femoris. The contact stress of patellofemoral joint was (1.29±0.41) MPa, which was similar to the results reported previously. The finite element model was valid. The von Mise stress nephogram showed that the stress mainly focused on the medial patellofemoral articular surface during knee flexion, and the contact point gradually moved up with the knee flexion deepened. The stress on the medial and lateral patellofemoral articular surface increased with the knee flexion deepened but decreased with the increase of patellar height. The effects of patellar height and knee flexion on the high contact stress of patellofemoral joint were similar among the finite element models after TKA based on the data of two volunteers. The high contact stress of patellofemoral joint increased with the knee flexion deepened in the same patellar height models (P<0.05), but decreased with the increase of patellar height in the same knee flexion models (P<0.05). The high contact stress of patellofemoral joint of pseudo-patella baja model was significantly higher than normal and alta patella models (P<0.05). The average contact area of patellofemoral joint of pseudo-patella baja was bigger than normal and alta patella models with the knee flexion deepened.ConclusionThe pseudo-patella baja after TKA has an important effect on the biomechanics of patellofemoral joint. Reserving the joint line and avoiding the occurrence of pseudo-patella baja can decrease the risk of anterior knee pain, patellar arthritis, and other complications caused by the increasing of contact stress of patellofemoral joint.
The purpose of this study is to analyze the biomechanics of ankle cartilage and ligaments during a typical Tai Chi movement–Brush Knee and Twist Step (BKTS). The kinematic and kinetic data were acquired in one experienced male Tai Chi practitioner while performing BKTS and in normal walking. The measured parameters were used as loading and boundary conditions for further finite element analysis. This study showed that the contact stress of the ankle joint during BKTS was generally less than that during walking. However, the maximum tensile force of the anterior talofibular ligament, the calcaneofibular ligament and the posterior talofibular ligament during BKTS was 130 N, 169 N and 89 N, respectively, while it was only 57 N, 119 N and 48 N during walking. Therefore, patients with arthritis of the ankle can properly practice Tai Chi. Practitioners with sprained lateral ligaments of the ankle joint were suggested to properly reduce the ankle movement range during BKTS.
Cholangiocarcinoma is a highly malignant tumor. It is not sensitive to radiotherapy and chemotherapy and has a poor prognosis. At present, there is no effective treatment. As a new method for treating cancer, magnetic fluid hyperthermia has been clinically applied to a variety of cancers in recent years. This article introduces it to the cholangiocarcinoma model and systematically studies the effect of magnetic fluid hyperthermia on cholangiocarcinoma. Starting from the theory of magnetic fluid heating, the electromagnetic and heat transfer models were constructed in the finite element simulation software COMSOL using the Pennes biological heat transfer equation. The Helmholtz coil was used as an alternating magnetic field generating device. The relationship between the magnetic fluid-related properties and the heating power was analyzed according to Rosensweig’s theory. After the multiphysics coupling simulation was performed, the electromagnetic field and thermal field distribution in the hyperthermia region were obtained. The results showed that the magnetic field distribution in the treatment area was uniform, and the thermal field distribution met the requirements of hyperthermia. After the magnetic fluid injection, the cholangiocarcinoma tissue warmed up rapidly, and the temperature of tumor tissues could reach above 42 °C, but the surrounding healthy tissues did not heat up significantly. At the same time, it was verified that the large blood vessels around the bile duct, the overflow of the magnetic fluid, and the eddy current heat had little effect on thermotherapy. The results of this article can provide a reference for the clinical application of magnetic fluid hyperthermia for cholangiocarcinoma.
The current finite element analysis of vascular stent expansion does not take into account the effect of the stent release pose on the expansion results. In this study, stent and vessel model were established by Pro/E. Five kinds of finite element assembly models were constructed by ABAQUS, including 0 degree without eccentricity model, 3 degree without eccentricity model, 5 degree without eccentricity model, 0 degree axial eccentricity model and 0 degree radial eccentricity model. These models were divided into two groups of experiments for numerical simulation with respect to angle and eccentricity. The mechanical parameters such as foreshortening rate, radial recoil rate and dog boning rate were calculated. The influence of angle and eccentricity on the numerical simulation was obtained by comparative analysis. Calculation results showed that the residual stenosis rates were 38.3%, 38.4%, 38.4%, 35.7% and 38.2% respectively for the 5 models. The results indicate that the pose has less effect on the numerical simulation results so that it can be neglected when the accuracy of the result is not highly required, and the basic model as 0 degree without eccentricity model is feasible for numerical simulation.
Objective To investigate the impact of bone mass and volume of low-density zones beneath the tibial plateau on the maximum von Mises stresses experienced by the cartilage and meniscus in the knee joint. Methods The study included one healthy adult volunteer, from whom CT scans were obtained, and one patient diagnosed with knee osteoarthrisis (KOA), for whom X-ray films were acquired. A static model of the knee joint featuring a low-density zone was established based on a normal knee model. In the finite element analysis, axial loads of 1 000 N and 1 800 N were applied to the weight-bearing region of the upper surface of the femoral head for model validation and subsequent finite element studies, respectively. The maximum von Mises stresses in the femoral cartilage, as well as the medial and lateral tibial cartilage and menisci, were observed, and the stress percentage of the medial and lateral components were concurrently analyzed. Additionally, HE staining, as well as alkaline magenta staining, were performed on the pathological specimens of patients with KOA in various low-density regions. ResultsThe results of model validation indicated that the model was consistent with normal anatomical structures and correlated with previous calculations documented in the literature. Static analysis revealed that the maximum von Mises stress in the medial component of the normal knee was the lowest and increased with the advancement of the hypointensity zone. In contrast, the lateral component exhibited an opposing trend, with the maximum von Mises stress in the lateral component being the highest and decreasing as the hypointensity zone progressed. Additionally, the medial component experienced an increasing proportion of stress within the overall knee joint. HE staining demonstrated that the chondrocyte layer progressively deteriorated and may even disappear as the hypointensity zone expanded. Furthermore, alkaline magenta staining indicated that the severity of microfractures in the trabecular bone increased concurrently with the expansion of the hypointensity zone. Conclusion The presence of subtalar plateau low-density zone may aggravate joint degeneration. In clinical practice, it is necessary to pay attention to the changes in the subtalar plateau low-density zone and actively take effective measures to strengthen the bone status of the subtalar plateau low-density zone and restore the complete biomechanical function of the knee joint, in order to slow down or reverse the progression of osteoarthritis.