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.
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.
Astronauts who are exposed to weightless environment in long-term spaceflight might encounter bone density and mass loss for the mechanical stimulus is smaller than normal value. This study built a three dimensional model of human femur to simulate the remodeling process of human femur during bed rest experiment based on finite element analysis (FEA). The remodeling parameters of this finite element model was validated after comparing experimental and numerical results. Then, the remodeling process of human femur in weightless environment was simulated, and the remodeling function of time was derived. The loading magnitude and loading cycle on human femur during weightless environment were increased to simulate the exercise against bone loss. Simulation results showed that increasing loading magnitude is more effective in diminishing bone loss than increasing loading cycles, which demonstrated that exercise of certain intensity could help resist bone loss during long-term spaceflight. At the end, this study simulated the bone recovery process after spaceflight. It was found that the bone absorption rate is larger than bone formation rate. We advise that astronauts should take exercise during spaceflight to resist bone loss.
The validated finite element head model (FEHM) of a 3-year-old child, a 6-year-old child and a 50th percentile adult were used to investigate the effects of head dimension and material parameters of brain tissues on the head rotational responses based on experimental design. Results showed that the effects of head dimension and directions of rotation on the head rotational responses were not significant under the same rotational loading condition, and the same results appeared in the viscoelastic material parameters of brain tissues. However, the head rotational responses were most sensitive to the shear modulus (G) of brain tissues relative to decay constant (β) and bulk modulus (K). Therefore, the selection of material parameters of brain tissues is most important to the accuracy of simulation results, especially in the study of brain injury criterion under the rotational loading conditions.
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.
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.
This study aims to overcome the shortcomings such as low efficiency, high cost and difficult to carry out multi-parameter research, which limited the optimization of infusion bag configuration and manufacture technique by experiment method. We put forward a fluid cavity based finite element method, and it could be used to simulate the stress distribution and deformation process of infusion bag under external load. In this paper, numerical models of infusion bag with different sizes was built, and the fluid-solid coupling deformation process was calculated using the fluid cavity method in software ABAQUS subject to the same boundary conditions with the burst test. The peeling strength which was obtained from the peeling adhesion test was used as failure criterion. The calculated resultant force which makes the computed peeling stress reach the peeling strength was compared with experiment data, and the stress distribution was analyzed compared with the rupture process of burst test. The results showed that considering the errors caused by the difference of weak welding and eccentric load, the flow cavity based finite element method can accurately model the stress distribution and deformation process of infusion bag. It could be useful for the optimization of multi chamber infusion bag configuration and manufacture technique, leading to cost reduction and study efficiency improvement.
Objective To establish a finite element model of the knee joint based on coronal plane alignment of the knee (CPAK) typing method, and analyze the biomechanical characteristics of different types of knee joints.Methods The finite element models of the knee joint were established based on CT scan data of 6 healthy volunteers. There were 5 males and 1 female with an average age of 24.2 years (range, 23-25 years). There were 3 left knees and 3 right knees. According to the CPAK typing method, the knees were rated as types Ⅰ to Ⅵ. Under the same material properties, boundary conditions, and axial loading, biomechanical simulations were performed on the finite element model of the knee joint. Based on the Von Mises stress nephogram and displacement nephogram, the peak stresses of the meniscus, femoral cartilage, and tibial cartilage, and the displacement of the meniscus were compared among different types of knee joints. Results The constructed finite element model of the knee joint was verified to be effective, and the stress and displacement results were consistent with previous literature. Under the axial load of 1 000 N, the stress nephogram showed that the stress distribution of the medial and lateral meniscus and tibial cartilage of CPAK type Ⅲ knee joint was the most uneven. The peak stresses of the lateral meniscus and tibial cartilage were 9.969 6 MPa and 2.602 7 MPa, which were 173% and 165% of the medial side, respectively. The difference of peak stress between the medial and lateral femoral cartilage was the largest in type Ⅳ knee joint, and the medial was 221% of the lateral. The displacement nephogram showed that the displacement of the medial meniscus was greater than that of the lateral meniscus except for types Ⅲ and Ⅵ knee joints. The difference between medial and lateral meniscus displacement of type Ⅲ knee joint was the largest, the lateral was 170% of the medial. Conclusion In the same type of joint line obliquity (JLO), the medial and lateral stress distribution of the knee was more uniform in varus and neutral positions than in valgus position. At the same time, the distal vertex of JLO subgroup can help to reduce the uneven medial and lateral stress distribution of varus knee, but increase the uneven distribution of valgus knee.
Males typically have high rates of morbidity of primary bladder neck obstruction, while the existing urodynamic examination is invasive and more likely to cause false diagnosis. To build a non-invasive biomechanical detecting system for the male lower urinary tract, a finite element model for male lower urinary tract based on the collodion slice images of normal male lower urinary tract was constructed, and the fluid-structure interaction of the lower urinary tract was simulated based on the real urination environment. The finite element model of the lower urinary tract was validated by comparing the clinical experiment data with the simulation result. The stress, flow rate and deformation of the lower urinary tract were analyzed, and the results showed that the Von Mises stress and the wall shear stress at the membrane sphincter in the normal male lower urinary tract model reached a peak, and there was nearly 1 s delay than in the bladder pressure, which helped to validate the model. This paper lays a foundation for further research on the urodynamic response mechanism of the bladder pressure and flow rate of the lower urinary tract obstruction model, which can provide a theoretical basis for the research of non-invasive biomechanical detecting system.
Based on the CT data and the structure characteristics of the femoral fractures during different healing stages, medical FE models of fractured femur treated with locking compression plate (LCP)were built.Under the physiological load of a standard body weight (70 kg) and the constraint condition,the stress distributions of LCP and fractured femur during healing were calculated by means of three-dimensional finite element analysis (3D-FEA).The results showed that the stress distribution in the LCP and the fractured femur was similar,during the initial stage which there was no newly formed bone or soft tissue in fracture site.The maximum von Mises stress (371.23,272.76 MPa) in the fractured femur was much higher than that in natural femur,and the intensive stress was concentrated mainly in the proximal area of the fractured femur.With the growth of bony callus bone in fracture site,the intensity of stress in proximal femur decreased.Contrasted to the two cases mentioned above,the value of the maximum von Mises stress (68.17 MPa) in bony callus bone stage decreased significantly,and was lower than the safe strength of natural bone.Therefore,appropriate training which is benefitial for the growth to new bone could be arranged for the better rehabilitation.