The lumbar intervertebral disc exhibits a complex physiological structure with interactions between various segments, and its components are extremely complex. The material properties of different components in the lumbar intervertebral disc, especially the water content (undergoing dynamic change as influenced by age, degeneration, mechanical loading, and proteoglycan content) - critically determine its mechanical properties. When the lumbar intervertebral disc is under continuous pressure, water seeps out, and after the pressure is removed, water re-infiltrates. This dynamic fluid exchange process directly affects the mechanical properties of the lumbar intervertebral disc, while previous isotropic modeling methods have been unable to accurately reflect such solid-liquid phase behaviors. To explore the load-bearing mechanism of the lumbar intervertebral disc and establish a more realistic mechanical model of the lumbar intervertebral disc, this study developed a solid-liquid biphasic, fiber-reinforced finite element model. This model was used to simulate the four movements of the human lumbar spine in daily life, namely flexion, extension, axial rotation, and lateral bending. The fluid pressure, effective solid stress, and liquid pressure-bearing ratio of the annulus fibrosus and nucleus pulposus of different lumbar intervertebral discs were compared and analyzed under the movements. Under all the movements, the fluid pressure distribution was closer to the nucleus pulposus, while the effective solid stress distribution was more concentrated in the outer annulus fibrosus. In terms of fluid pressure, the maximum fluid pressure of the lumbar intervertebral disc during lateral bending was 1.95 MPa, significantly higher than the maximum fluid pressure under other movements. Meanwhile, the maximum effective solid stress of the lumbar intervertebral disc during flexion was 2.43 MPa, markedly higher than the maximum effective solid stress under other movements. Overall, the liquid pressure-bearing ratio under axial rotation was smaller than that under other movements. Based on the solid-liquid biphasic modeling method, this study more accurately revealed the dominant role of the liquid phase in the daily load-bearing process of the lumbar intervertebral disc and the solid-phase mechanical mechanism of the annulus fibrosus load-bearing, and more effectively predicted the solid-liquid phase co-load-bearing mechanism of the lumbar intervertebral disc in daily life.
In unicompartmental replacement surgery, there are a wide variety of commercially available unicompartmental prostheses, and the consistency of the contact surface between the common liner and the femoral prosthesis could impact the stress distribution in the knee after replacement in different ways. Medial tibial plateau fracture and liner dislocation are two common forms of failure after unicompartmental replacement. One of the reasons is the mismatch in the mounting position of the unicompartmental prosthesis in the knee joint, which may lead to failure. Therefore, this paper focuses on the influence of the shape of the contact surface between the liner and the femoral prosthesis and the mounting position of the unicompartmental prosthesis on the stress distribution in the knee joint after replacement. Firstly, a finite element model of the normal human knee joint was established, and the validity of the model was verified by both stress and displacement. Secondly, two different shapes of padded knee prosthesis models (type A and type B) were developed to simulate and analyze the stress distribution in the knee joint under single-leg stance with five internal or external rotation mounting positions of the two pads. The results showed that under a 1 kN axial load, the peak contact pressure of the liner, the peak ACL equivalent force, and the peak contact pressure of the lateral meniscus were smaller for type A than for type B. The liner displacement, peak contact pressure of the liner, peak tibial equivalent force, and peak ACL equivalent force were the smallest for type A at 3° of internal rotation in all five internal or external rotation mounting positions. For unicompartmental replacement, it is recommended that the choice of type A or type B liner for prosthetic internal rotation up to 6° should be combined with other factors of the patient for comprehensive analysis. In conclusion, the results of this paper may reduce the risk of liner dislocation and medial tibial plateau fracture after unicompartmental replacement, providing a biomechanical reference for unicompartmental prosthesis design.
Objective To discuss the method of constructingbiomechanical model of rabbit femur.Methods The sample of rabbit femur was prepared as follows:firstly，femur section images were obtained，then the image wasput into the computer and processed to get the boundary contour line; secondly, through programming the contour line coordinate for modeling was obtained, then the data were put into the model software to find the threedimensional entity model. Results Whole three-dimensional model of rabbit femur was constructed. It simulated actually dissection form of femur. Conclusion The establishment of the model lays a foundation for ascertaining optimal parameter of vibration improving bone minerydensity by finite element analysis.
For the transportation process of rescuing wounded personnel on naval vessels, a new type of shoulder type exoskeleton stretcher for individual soldier was designed in this paper. The three-dimensional model of the shoulder type exoskeleton stretcher for individual soldier was constructed using three dimensional modeling software. Finite element analysis technique was employed to conduct statics simulation, modal analysis, and transient dynamics analysis on the designed exoskeleton stretcher. The results show that the maximum stress of the exoskeleton stretcher for walking on flat ground is 265.55 MPa, which is lower than the allowable strength of the fabrication material. Furthermore, the overall deformation of the structure is small. Modal analysis reveals that the natural frequency range of the exoskeleton stretcher under different gait conditions is 1.96 Hz to 28.70 Hz, which differs significantly from the swing frequency of 1 Hz during walking. This indicates that the designed structure can effectively avoid resonance. The transient dynamics analysis results show that the maximum deformation and stress of exoskeleton stretcher remain within the safety range, which meets the expected performance requirements. In summary, the shoulder type exoskeleton stretcher for individual soldier designed in this study can solve the problem of requiring more than 2 people to carry for the existing stretcher, especially suitable for narrow spaces of naval vessels. The research results of this paper can provide a new solution for the rescue of wounded personnel on naval vessels.
Although metal blocks have been widely used for reconstructing uncontained tibial bone defects, the influence of their elastic modulus on the stability of tibial prosthesis fixation remains unclear. Based on this, a finite element model incorporating constrained condylar knee (CCK) prosthesis, tibia, and metal block was established. Considering the influence of the post-restraint structure of the prosthesis, the effects of variations in the elastic modulus of the block on the von Mises stress distribution in the tibia and the block, as well as on the micromotion at the bone-prosthesis fixation interface, were investigated. Results demonstrated that collision between the insert post and femoral prosthesis during tibial internal rotation increased tibial von Mises stress, significantly influencing the prediction of block elastic modulus variation. A decrease in the elastic modulus of the metal block resulted in increased von Mises stress in the proximal tibia, significantly reduced von Mises stress in the distal tibia, decreased von Mises stress of the block, and increased micromotion at the bone-prosthesis fixation interface. When the elastic modulus of the metal block fell below that of bone cement, inadequate block support substantially increased the risk of stress shielding in the distal tibia and fixation interface loosening. Therefore, this study recommends that biomechanical investigations of CCK prostheses must consider the post-constraint effect, and the elastic modulus of metal blocks for bone reconstruction should not be lower than 3 600 MPa.
Objective To review recent advance in the application and research of three-dimensional digital knee model. Methods The recent original articles about three-dimensional digital knee model were extensively reviewed and analyzed. Results The digital three-dimensional knee model can simulate the knee complex anatomical structure very well. Based on this, there are some developments of new software and techniques, and good clinical results are achieved. Conclusion With the development of computer techniques and software, the knee repair and reconstruction procedure has been improved, the operation will be more simple and its accuracy will be further improved.
Objective To investigate the effect of first to third metatarsus defect and the effect of reconstruction with ilium on foot function. Methods The first to third metatarsus defect was simulated in a 3D foot model and rebuilt by ilium. The maximal displacement and stress calculated by the method of finite elements were used as the index of estimation. Five cases treated from Mar. 1996 to Jan. 2003 with metatarsus defect rebuilding by free vascular iliac bone incorporating free flapwere evaluated. Results Foot function was affected largely by the defect of the first to third metatarsus. Compared with the normal foot, the maximal displacement was increased by 2.15 times and the maximal stress was increased by 2.12 times in 100% defected foot, and in 50%-defected foot maximal displacement and stress were increased by 1.65 times and 2.05 times respectively. Follow-up had been conducted for 1 to 2 years. All bones and flaps of the 5 cases survived (2 excellent, 2 good, and 1 passable) by function evaluation. Conclusion The first to third metatarsus defect should be repaired, and the method of transplanting iliac bone added by flap is effective.
Objective To investigate the stability and the stress distributions of L3-5 fused with three different approaches (interbody, posterolateral and circumferential fusions) and to investigate degeneration of thesegment adjacent to the fused functional spinal unit. Methods A detailed L3-5 three-dimensional nonlinear finite element model of a normal man aged 32 was established and validated. Based on the model, the destabilized model, the interbody, posterolateral and circumferential fusions models of L4-5 were established. After the loadings were placed on all the models, we recorded the angular motions of the fused segment and the Von Mises stress of the adjacent intervertebral disc. Results The circumferential fusion was most stable than the others, and the interbody fusion was more stable than the posterolateral fusion. The maximal Von Mises stress of the adjacent L3,4 intervertebral disc in all the models was ranked descendingly as flexion,lateral bending,torsion and extension. For the three kinds of fusions, the stress increment of the L3,4 intervertebral disc was ranked ascendingly as interbody fusion,posterolateral fusion and circumferential fusion. Conclusion After destabilization of the L4,5 segment, the stability of the circumferential fusionis better than that of the others, particularly under the flexional or extensional loading. The stability of the interbody fusion is better than that of the posterolateral fusion, except for under the flexional loading. The feasibility of adjacent segment degeneration can be ranked descendingly as: circumferential fusion,posterolateral fusion and interbody fusion.
The finite element method is a new method to study the mechanism of brain injury caused by blunt instruments. But it is not easy to be applied because of its technology barrier of time-consuming and strong professionalism. In this study, a rapid and quantitative evaluation method was investigated to analyze the craniocerebral injury induced by blunt sticks based on convolutional neural network and finite element method. The velocity curve of stick struck and the maximum principal strain of brain tissue (cerebrum, corpus callosum, cerebellum and brainstem) from the finite element simulation were used as the input and output parameters of the convolutional neural network The convolutional neural network was trained and optimized by using the 10-fold cross-validation method. The Mean Absolute Error (MAE), Mean Square Error (MSE), and Goodness of Fit (R2) of the finally selected convolutional neural network model for the prediction of the maximum principal strain of the cerebrum were 0.084, 0.014, and 0.92, respectively. The predicted results of the maximum principal strain of the corpus callosum were 0.062, 0.007, 0.90, respectively. The predicted results of the maximum principal strain of the cerebellum and brainstem were 0.075, 0.011, and 0.94, respectively. These results show that the research and development of the deep convolutional neural network can quickly and accurately assess the local brain injury caused by the sticks blow, and have important application value for understanding the quantitative evaluation and the brain injury caused by the sticks struck. At the same time, this technology improves the computational efficiency and can provide a basis reference for transforming the current acceleration-based brain injury research into a focus on local brain injury research.
Objective To compare the biomechanical properties of the anterior transpedicular screw-artificial vertebral body (AVB) and conventional anterior screw plate system (AP) in lower cervical spine by finite element study. Methods CT images (C1-T1) were obtained from a 38-year-old female volunteer. The models of intact C3-7 (intact group), AP fixation (AP group), and AVB fixation (AVB group) were established and analyzed by Mimics 14.0, Geomagic Studio 2013, and ANSYS 14.0 softwares. The axial force of 74 N and moment couple of 1 N·m were loaded on the upper surface and upper facet joint surfaces of C3. Under conditions of flexion, extension, lateral bending, and rotation, the Von Mises stress distribution regularity and maximum equivalent stree of AP and AVB groups were recorded, and the range of motion (ROM) was also analyzed of 3 groups. Results The intact model of lower cervical spine (C3-7) was established, consisting of 286 382 elements and 414 522 nodes, and it was successfully validated with the previously reported cadaveric experimental data of Panjabi and Kallemeyn. The stress concentrated on the connection between plate and screw in AP group, while it distributed evenly in AVB group. Between AP and AVB groups, there was significant difference in maximum equivalent stress values under conditions of 74 N axial force, flexion, extension, and rotation. AVB group had smaller ROM of fixed segments and larger ROM of adjacent segments than AP group. Compared with intact group, whole ROM of the lower cervical spine decreased about 3°, but ROM of C3, 4 and C6, 7 segments increased nearly 5° in both AP and AVB groups. Conclusion As a new reconstruction method of lower cervical spine, AVB fixation provides better stability and lower risk of failure than AP fixation.