Objective To compare the biomechanical differences among the five internal fixation modes in treatment of Day type Ⅱ crescent fracture dislocation of pelvis (CFDP), and find an internal fixation mode which was the most consistent with mechanical principles. Methods Based on the pelvic CT data of a healthy adult male volunteer, a Day type Ⅱ CFDP finite element model was established by using Mimics 17.0, ANSYS 12.0-ICEM, Abaqus 2020, and SolidWorks 2012 softwares. After verifying the validity of the finite element model by comparing the anatomical parameters with the three-dimensional reconstruction model and the mechanical validity verification, the fracture and dislocated joint of models were fixed with S1 sacroiliac screw combined with 1 LC-Ⅱ screw (S1+LC-Ⅱ group), S1 sacroiliac screw combined with 2 LC-Ⅱ screws (S1+2LC-Ⅱ group), S1 sacroiliac screw combined with 2 posterior iliac screws (S1+2PIS group), S1 and S2 sacroiliac screws combined with 1 LC-Ⅱ screw (S1+S2+LC-Ⅱ group), S2-alar-iliac (S2AI) screw combined with 1 LC-Ⅱ screw (S2AI+LC-Ⅱ group), respectively. After each internal fixation model was loaded with a force of 600 N in the standing position, the maximum displacement of the crescent fracture fragments, the maximum stress of the internal fixation (the maximum stress of the screw at the ilium fracture and the maximum stress of the screw at the sacroiliac joint), sacroiliac joint displacement, and bone stress distribution around internal fixation were observed in 5 groups. Results The finite element model in this study has been verified to be effective. After loading 600 N stress, there was a certain displacement of the crescent fracture of pelvis in each internal fixation model, among which the S1+LC-Ⅱ group was the largest, the S1+2LC-Ⅱ group and the S1+2PIS group were the smallest. The maximum stress of the internal fixation mainly concentrated at the sacroiliac joint and the fracture line of crescent fracture. The maximum stress of the screw at the sacroiliac joint was the largest in the S1+LC-Ⅱ group and the smallest in the S2AI+LC-Ⅱ group. The maximum stress of the screw at the ilium fracture was the largest in the S1+2PIS group and the smallest in the S1+2LC-Ⅱ group. The displacement of the sacroiliac joint was the largest in the S1+LC-Ⅱ group and the smallest in the S1+S2+LC-Ⅱ group. In each internal fixation model, the maximum stress around the sacroiliac screws concentrated on the contact surface between the screw and the cortical bone, the maximum stress around the screws at the iliac bone concentrated on the cancellous bone of the fracture line, and the maximum stress around the S2AI screw concentrated on the cancellous bone on the iliac side. The maximum bone stress around the screws at the sacroiliac joint was the largest in the S1+LC-Ⅱ group and the smallest in the S2AI+LC-Ⅱ group. The maximum bone stress around the screws at the ilium was the largest in the S1+2PIS group and the smallest in the S1+LC-Ⅱ group. Conclusion For the treatment of Day type Ⅱ CFDP, it is recommended to choose S1 sacroiliac screw combined with 1 LC-Ⅱ screw for internal fixation, which can achieve a firm fixation effect without increasing the number of screws.
Aiming at the problem of scaffold degradation in bone tissue engineering, we studied the feasibility that controlls bone defect repair effect with the inhomogeneous structure of scaffold. The prediction model of bone defect repair which contains governing equations for bone formation and scaffold degradation was constructed on the basis of analyzing the process and main influence factors of bone repair in bone tissue engineering. The process of bone defect repair and bone structure after repairing can be predicted by combining the model with finite element method (FEM). Bone defect repair effects with homogenous and inhomogeneous scaffold were simulated respectively by using the above method. The simulation results illustrated that repair effect could be impacted by scaffold structure obviously and it can also be controlled via the inhomogeneous structure of scaffold with some feasibility.
In the study of oral orthodontics, the dental tissue models play an important role in finite element analysis results. Currently, the commonly used alveolar bone models mainly have two kinds: the uniform and the non-uniform models. The material of the uniform model was defined with the whole alveolar bone, and each mesh element has a uniform mechanical property. While the material of the elements in non-uniform model was differently determined by the Hounsfield unit (HU) value of computed tomography (CT) images where the element was located. To investigate the effects of different alveolar bone models on the biomechanical responses of periodontal ligament (PDL), a clinical patient was chosen as the research object, his mandibular canine, PDL and two kinds of alveolar bone models were constructed, and intrusive force of 1 N and moment of 2 Nmm were exerted on the canine along its root direction, respectively, which were used to analyze the hydrostatic stress and the maximal logarithmic principal strain of PDL under different loads. Research results indicated that the mechanical responses of PDL had been affected by alveolar bone models, no matter the canine translation or rotation. Compared to the uniform model, if the alveolar bone was defined as the non-uniform model, the maximal stress and strain of PDL were decreased by 13.13% and 35.57%, respectively, when the canine translation along its root direction; while the maximal stress and strain of PDL were decreased by 19.55% and 35.64%, respectively, when the canine rotation along its root direction. The uniform alveolar bone model will induce orthodontists to choose a smaller orthodontic force. The non-uniform alveolar bone model can better reflect the differences of bone characteristics in the real alveolar bone, and more conducive to obtain accurate analysis results.
In order to check the neck response and injury during motor vehicle accidents, we developed a detailed finite element model for human cervical spine C4-C6. This model consisted of cortical bone, cancellous bone, annulus, nucleus, ligaments and articular facet, and it also set up contact in the contacting parts for simulating the movement perfectly under frontal impact. This model could be used for stress and strain distribution after the frontal impact load was applied on this model. During the process of frontal impact, the most displacement simulated data were in the interval range of experimental data. The experimental results showed that this model for the human cervical spine C4-C6 simulated the movement under the frontal impact with fidelity, and reflected the impact dynamics response on the whole.
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.
A three-dimensional finite element model of premaxillary bone and anterior teeth was established with ANSYS 13.0. The anterior teeth were fixed with strong stainless labial archwire and lingual frame. In the horizontal loading experiments, a horizontal retraction force of 1.5 N was applied bilaterally to the segment through hooks at the same height between 7 and 21 mm from the incisal edge of central incisor; in vertical loading experiments, a vertical intrusion force of 1.5 N was applied at the midline of lingual frame with distance between 4 and 16 mm from the incisal edge of central incisor. After loading, solution was done and displacement and maximum principle stress were calculated. After horizontal loading, lingual displacement and stress in periodontal membrane (PDM) was most homogeneous when the traction force was 14 mm from the edge of central incisor; after vertical loading, intrusive displacement and stress in PDM were most homogeneous when the traction force was 12 mm from the incisal edge of central incisor. The results of this study suggested that the location of center of resistance (CRe) of six maxillary anterior teeth is about 14 mm gingivally and 12 mm lingually to incisal edge of central incisor. The location can provide evidence for theoretical and clinical study in orthodontics.
We developed a three-dimensional finite element model of development dysplasia of hip (DDH) of a patient. And then we performed virtual Bernese periacetabular osteotomy (PAO) by rotating the acetabular bone with different angle so as to increase femoral head coverage and distribute the contact pressure over the cartilage surface. Using finite element analysis method, we analyzed contact area, contact pressure, and von Mises stress in the acetabular cartilage to determine the effect of various rotation angle. We also built a normal hip joint model. Compared to the normal hip joint model, the DDH models showed stress concentration in the acetabular edge, and higher stress values. Compared to the DDH models, the post-PAO models showed decreases in the maximum values of von Mises stress and contact pressure while we increased the contact area. An optimal position could be achieved for the acetabulum that maximizes the contact area while minimizing the contact pressure and von Mises stress in the acetabular cartilage. These would provide theoretical bases to pre-operative planning.
ObjectiveTo explore the biomechanical characteristics and clinical application effects of three-dimensional (3D) printed osteotomy guide plate combined with Ilizarov technique in the treatment of rigid clubfoot. Methods A retrospective analysis was performed on the clinical data of 11 patients with rigid clubfoot who met the inclusion criteria and were admitted between January 2019 and December 2024. There were 6 males and 5 females, aged 21-60 years with an average of 43.2 years. Among them, 5 cases were untreated congenital rigid clubfoot, 4 cases were recurrent rigid clubfoot after previous treatment, and 2 cases were rigid clubfoot due to disease sequelae. All 11 patients first received slow distraction using Ilizarov technique combined with circular external fixator until the force lines of the foot and ankle joint were basically normal. Then, 1 male patient aged 24 years was selected, and CT scanning was used to obtain imaging data of the ankle joint and foot. A 3D finite element model was established and validated using the plantar stress distribution nephogram of the patient. After validation, the biomechanical changes of the tibiotalar joint under the same load were simulated after triple arthrodesis and fixation. The optimal correction angle of the hindfoot was determined to fabricate 3D-printed osteotomy guide plates, and all 11 patients underwent triple arthrodesis using these guide plates. The functional recovery was evaluated by comparing the American Orthopaedic Foot and Ankle Society (AOFAS) score, International Clubfoot Study Group (ICFSG) score, and 36-Item Short Form Survey (SF-36) score before and after operation. Results Finite element analysis showed that the maximum peak von Mises stress of the tibiotalar joint was at hindfoot varus 3° and the minimum at valgus 6°; the maximum peak von Mises stress of the 3 naviculocuneiform joints under various conditions appeared at lateral naviculocuneiform joint before operation, and the minimum appeared at lateral naviculocuneiform joint at neutral position 0°; the maximum peak von Mises stress of the 5 tarsometatarsal joints under various conditions appeared at the 2nd tarsometatarsal joint at hindfoot neutral position 0°, and the minimum appeared at the 1st tarsometatarsal joint at valgus 6°. Clinical application results showed that the characteristics of clubfoot deformity observed during operation were consistent with the preoperative 3D reconstruction model. All 11 patients were followed up 8-24 months with an average of 13.1 months. One patient had postoperative incision exudation, which healed after dressing change; the remaining patients had good incision healing. All patients achieved good healing of the osteotomy segments, with a healing time of 3-6 months and an average of 4.1 months. At last follow-up, the AOFAS score, SF-36 score, and ICFSG score significantly improved when compared with those before operation (P<0.05). ConclusionThe 3D-printed osteotomy guide plate combined with Ilizarov technique has favorable biomechanical advantages in the treatment of rigid clubfoot, with significant clinical application effects. It can effectively improve the foot function of patients and achieve precise and personalized treatment.
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 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.