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.
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.
In the present study, a finite element model of L4-5 lumbar motion segment was established based on the CT images and a combination with image processing software, and the analysis of lumbar biomechanical characteristics was conducted on the proposed model according to different cases of flexion, extension, lateral bending and axial rotation. Firstly, the CT images of lumbar segment L4 to L5 from a healthy volunteer were selected for a three dimensional model establishment which was consisted of cortical bone, cancellous bone, posterior structure, annulus, nucleus pulposus, cartilage endplate, ligament and facet joint. The biomechanical analysis was then conducted according to different cases of flexion, extension, lateral bending and axial rotation. The results showed that the established finite element model of L4-5 lumbar segment was realistic and effective. The axial displacement of the proposed model was 0.23, 0.47, 0.76 and 1.02 mm, respectively under the pressure of 500, 1 000, 1 500 and 2 000 N, which was similar to the previous studies in vitro experiments and finite element analysis of other people under the same condition. The stress distribution of the lumbar spine and intervertebral disc accorded with the biomechanical properties of the lumbar spine under various conditions. The established finite element model has been proved to be effective in simulating the biomechanical properties of lumbar spine, and therefore laid a good foundation for the research of the implants of biomechanical properties of lumbar spine.
Objective To investigate the validity of improving the femur’s mechanical characteristics by implanting calcium phosphate ceramic screws after removing dynamic hip screw (DHS). Methods The three dimensional finite element model of the femur was built based on the CT scanning of a normal male volunteer. Then the models of the femur with and without DHS were established. According to calcium phosphate ceramic screws with porosity and apparent elastic modulus, 80% and 0.1 GPa were set as group A, 50% and 1.0 GPa as group B, and 30% and 1.5 GPa as group C. Von Mises stress distribution and maximum stress were recorded when the joint was maximally loaded in a gait cycle. Results The Von Mises in normal femoral shaft was uniform; no phenomena of stress concentration was observed and the maximum stress located at the joint load-bearing site of the proximal femur. The stress concentration was observed in the femur without DHS, and the maximum stress located at the distal femur around the screw hole. By comparing several different calcium phosphate ceramic screws, the stress distribution of group B was similar to normal femur model, and the maximum stress located at the joint load-bearing site. The other screws of groups A and C showed varying degrees of stress concentration. Conclusion Implanting calcium phosphate ceramic screw can improve the mechanical characteristics of the femur after removing dynamic hip screw, and the calcium phosphate ceramic screw with 50% porosity and 1.0 GPa apparent elastic modulus is suitable for implanting.
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.
Stress distribution of denture is an important criterion to evaluate the reasonableness of technological parameters, and the bite force derived from the antagonist is the critical load condition for the calculation of stress distribution. In order to improve the accuracy of stress distribution as much as possible, all-ceramic crown of the mandibular first molar with centric occlusion was taken as the research object, and a bite force loading method reflecting the actual occlusal situation was adopted. Firstly, raster scanning and three dimensional reconstruction of the occlusal surface of molars in the standard dental model were carried out. Meanwhile, the surface modeling of the bonding surface was carried out according to the preparation process. Secondly, the parametric occlusal analysis program was developed with the help of OFA function library, and the genetic algorithm was used to optimize the mandibular centric position. Finally, both the optimized case of the mesh model based on the results of occlusal optimization and the referenced case according to the cusp-fossa contact characteristics were designed. The stress distribution was analyzed and compared by using Abaqus software. The results showed that the genetic algorithm was suitable for solving the occlusal optimization problem. Compared with the reference case, the optimized case had smaller maximum stress and more uniform stress distribution characteristics. The proposed method further improves the stress accuracy of the prosthesis in the finite element model. Also, it provides a new idea for stress analysis of other joints in human body.
The risk of vertebral cortical shell fracture increases with aging. However, it remains unclear how aging contributes to cortex fracture at present. The aim of this study is to make understanding of the mechanism of how the spinal aging influences the cortical shell strain. Two finite element (FE) models of spinal segments (mildly and fully aged) were created, and then were compared to the FE models of the healthy spinal segment. The FE models of the aged spinal segments were generated by modifying both the geometry of the intervertebral disc (IVD) and the material properties of the spinal components. To find out under which case the cortical shell strain was influenced more, we created two types of FE model comparison methods: one with changes only in the spinal material properties and the other with changes only in the IVD geometry. The results showed that the cortical shell strains increased with aging and that compared to the changes of IVD geometry, the changes of spinal material property have a higher influence on the cortical shell strains. This study may suggest that for the prevention and treatment of vertebral cortex fracture, the augmentation of the vertebral body is a more effective treatment.
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.
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.
The human spine injury and various lumbar spine diseases caused by vibration have attracted extensive attention at home and abroad. To explore the biomechanical characteristics of different approaches for lumbar interbody fusion surgery combined with an interspinous internal fixator, device for intervertebral assisted motion (DIAM), finite element models of anterior lumbar interbody fusion (ALIF), transforaminal lumbar interbody fusion (TLIF) and lateral lumbar interbody fusion (LLIF) are created by simulating clinical operation based on a three-dimensional finite element model of normal human whole lumbar spine. The fusion level is at L4–L5, and the DIAM is implanted between spinous process of L4 and L5. Transient dynamic analysis is conducted on the ALIF, TLIF and LLIF models, respectively, to compute and compare their stress responses to an axial cyclic load. The results show that compared with those in ALIF and TILF models, contact forces between endplate and cage are higher in LLIF model, where the von-Mises stress in endplate and DIAM is lower. This implies that the LLIF have a better biomechanical performance under vibration. After bony fusion between vertebrae, the endplate and DIAM stresses for all the three surgical models are decreased. It is expected that this study can provide references for selection of surgical approaches in the fusion surgery and vibration protection for the postsurgical lumbar spine.