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.
Modified electroconvulsive therapy (MECT) and magnetic seizure therapy (MST) are effective treatments for severe major depression. MECT has better efficacy in the treatment than MST, but it has cognitive and memorial side effects while MST does not. To study the causes of these different outcomes, this study contrasted the electric filed strength and spatial distribution induced by MECT and MST in a realistic human head model. Electric field strength induced by MECT and MST are simulated by the finite element method, which was based on a realistic human head model obtained by magnetic resonance imaging. The electrode configuration of MECT is standard bifrontal stimulation configuration, and the coil configuration of MST is circular. Maps of the ratio of the electric field strength to neural activation threshold are obtained to evaluate the stimulation strength and stimulation focality in brain regions. The stimulation strength induced by MECT is stronger than MST, and the activated region is wider. MECT stimulation strength in gray matter is 17.817 times of that by MST, and MECT stimulation strength in white matter is 23.312 times of that by MST. As well, MECT stimulation strength in hippocampi is 35.162 times of that by MST. More than 99.999% of the brain volume is stimulated at suprathreshold by MECT. However, MST activated only 0.700% of the brain volume. The stimulation strength induced by MECT is stronger than MST, and the activated region is wider may be the reason that MECT has better effectiveness. Nevertheless, the stronger stimulation strength in hippocampi induced by MECT may be the reason that MECT is more likely to give rise to side effects. Based on the results of this study, it is expected that a more accurate clinical quantitative treatment scheme should be studied in the future.
Numerical simulation of stent deployment is very important to the surgical planning and risk assess of the interventional treatment for the cardio-cerebrovascular diseases. Our group developed a framework to deploy the braided stent and the stent graft virtually by finite element simulation. By using the framework, the whole process of the deployment of the flow diverter to treat a cerebral aneurysm was simulated, and the deformation of the parent artery and the distributions of the stress in the parent artery wall were investigated. The results provided some information to improve the intervention of cerebral aneurysm and optimize the design of the flow diverter. Furthermore, the whole process of the deployment of the stent graft to treat an aortic dissection was simulated, and the distributions of the stress in the aortic wall were investigated when the different oversize ratio of the stent graft was selected. The simulation results proved that the maximum stress located at the position where the bare metal ring touched the artery wall. The results also can be applied to improve the intervention of the aortic dissection and the design of the stent graft.
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.
This paper adopted UG8.0 to bulid the stent and blood vessel models. The models were then imported into the finite element analysis software ANSYS. The simulation results of ANSYS software showed that after endothelial stent implantation, the velocity of the blood was slow and the fluctuation of velocity was small, which meant the flow was relatively stable. When blood flowed through the endothelial stent, the pressure gradually became smaller, and the range of the pressure was not wide. The endothelial shear stress basically unchanged. In general, it can be concluded that the endothelial stents have little impact on the flow of blood and can fully realize its function.
The objective of this study was to determine the visco-hyperelastic constitutive law of brain tissue under dynamic impacts. A method combined by finite element simulations and optimization algorithm was employed for the determination of material variables. Firstly, finite element simulations of brain tissue dynamic uniaxial tension, with a maximum stretch rate of 1.3 and strain rates of 30 s–1 and 90 s–1, were developed referring to experimental data. Then, fitting errors between the engineering stress-strain curves predicted by simulations and experimental average curves were assigned as objective functions, and the multi-objective genetic algorithm was employed for the optimation solution. The results demonstrate that the brain tissue finite element models assigned with the novel obtained visco-hyperelastic material law could predict the brain tissue’s dynamic mechanical characteristic well at different loading rates. Meanwhile, the novel material law could also be applied in the human head finite element models for the improvement of the biofidelity under dynamic impact loadings.