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 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.
As a representative part of the oral system and masticatory robot system, the modeling method of the dental model is an important factor influencing the accuracy of the multi-body dynamic model. Taking the right first molars of the masticatory robot as the research object, an equivalent model, point-contact higher kinematic pair composed of v-shaped surface and sphere surface, was proposed. Firstly, the finite element method was used to analyze the occlusal dynamics of the original model in three static contact cases (intrusive contact, centric occlusion, and extrusive contact) and one dynamic chewing case, and the expected bite force was obtained. Secondly, the Hertz contact model was adopted to establish the analytical expression of the bite force of the equivalent model in three static contact cases. The normal vectors and contact stiffness in the expression were designed according to the expected bite force. Finally, the bite force performance of the equivalent model in three static contact cases and one dynamic chewing case was evaluated. The results showed that the equivalent model could achieve the equivalent bite force of 8 expected items in the static contact cases. Meanwhile, the bite force in the early and late stages of the dynamic chewing case coincides well with the original model. In the middle stage, a certain degree of impact is introduced, but it can be weakened by subsequent trajectory planning. The equivalent modeling scheme of the dental model proposed in this paper further improves the accuracy of the dynamic model of the multi-body system. It provides a new idea for the dynamic modeling of other complex human contacts.
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.
We developed a three-dimensional finite element model of the pelvis. According to Letournel methods, we established a pelvis model of T-shaped fracture with its three different fixation systems, i.e. double column reconstruction plates, anterior column plate combined with posterior column screws and anterior column plate combined with quadrilateral area screws. It was found that the pelvic model was effective and could be used to simulate the mechanical behavior of the pelvis. Three fixation systems had great therapeutic effect on the T-shaped fracture. All fixation systems could increase the stiffness of the model, decrease the stress concentration level and decrease the displacement difference along the fracture line. The quadrilateral area screws, which were drilled into cortical bone, could generate beneficial effect on the T-type fracture. Therefore, the third fixation system mentioned above (i.e. the anterior column plate combined with quadrilateral area screws) has the best biomechanical stability to the T-type fracture.
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 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.
Craniofacial malformation caused by premature fusion of cranial suture of infants has a serious impact on their growth. The purpose of skull remodeling surgery for infants with craniosynostosis is to expand the skull and allow the brain to grow properly. There are no standardized treatments for skull remodeling surgery at the present, and the postoperative effect can be hardly assessed reasonably. Children with sagittal craniosynostosis were selected as the research objects. By analyzing the morphological characteristics of the patients, the point cloud registration of the skull distortion region with the ideal skull model was performed, and a plan of skull cutting and remodeling surgery was generated. A finite element model of the infant skull was used to predict the growth trend after remodeling surgery. Finally, an experimental study of surgery simulation was carried out with a child with a typical sagittal craniosynostosis. The evaluation results showed that the repositioning and stitching of bone plates effectively improved the morphology of the abnormal parts of the skull and had a normal growth trend. The child’s preoperative cephalic index was 65.31%, and became 71.50% after 9 months’ growth simulation. The simulation of the skull remodeling provides a reference for surgical plan design. The skull remodeling approach significantly improves postoperative effect, and it could be extended to the generation of cutting and remodeling plans and postoperative evaluations for treatment on other types of craniosynostosis.
In this paper, we established magnetic fluid hyperthermia (MFH) model for rat tumor using the finite element software COMSOL based on the linear response theory. By analyzing four kinds of magnetic medium within relaxation mechanism, such as Fe3O4、FeCo、fccFePt and L10FePt, we studied the influence of the change of magnetic medium radius on dissipation power and temperature field, respectively. At the same time, the optimization method for the parameters of several magnetic medium is proposed, and the applications of four kinds of magnetic medium are given as well. By increasing the dissipation power of the magnetic medium as much as possible, the dose of magnetic medium used in the treatment can be reduced, meanwhile, the adverse effects on health tissue surrounding the tumor will be minimized. The conclusions of this paper can provide reference for magnetic medium preparation applied to MFH.
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.