1.Three-dimensional finite element feature analysis of the mandible and morphology and position of temporomandibular joint in patients with unilateral and bilateral molar scissor bite.
Tianhao CHU ; Xueying ZHANG ; Haocheng WANG ; Haojie MA ; Yuanyuan LIU
West China Journal of Stomatology 2025;43(1):114-125
OBJECTIVES:
The objective of this study is to measuring the morphology and position of bilateral temporomandibular joints in patients with unilateral and bilateral molar scissor bite and simulating the deformation of the mandible during occlusion, in order to provide thesis for the diagnosis of temporomandibular joint disease in patients with unilateral and bilateral molar scissor bite.
METHODS:
This study was a retrospective study. A total of 10 patients with unilateral molar scissor bite (the unilateral molar scissor bite group) and 10 patients with bilateral molar scissor bite (the bilateral molar scissor bite group) were selected as the experimental group, and 20 adult patients with classⅠ of angle classification of similar ages were selected as the control group. All patients underwent cone beam computed tomography scans, by measuring the width of the fossa, height of the fossa, articular eminence inclination, long axis of the condyle, minor axis of the condyle, horizontal angle of the condyle and the space of the temporomandibular joint, compare temporomandibular joint morphology and position. The three-dimensional finite element analysis of the mandible morphology was carried out to evaluate the force and deformation of the mandible by using software to simulate the occlusion of the patients. It was further explored the relationship between the force of the mandible morphology and the possible temporomandibular joint disorder symptoms of the patients.
RESULTS:
Intergroup comparisons for the unilateral molar scissor bite group and left sides of the other groups revealed that the superior articular space in the group with unilateral molar scissor bite was shorter than that in the control group (P<0.05); the long axis of the condyle in the unilateral and bilateral molar scissor bite group were both shorter than that of the control group (P<0.05); among which the unilateral group was larger than the bilateral group, and the minor axis of the condyle in bilateral molar scissor bite group was smaller than in the control group (P<0.05), and the unilateral and bilateral condylar groups were larger than the control group (P<0.05); and the condylar horizontal angle in the unilateral and bilateral groups were larger than that in the control group (P<0.05). The normal sides of the unilateral molar scissor bite group and right sides of the other groups had smaller superior articular space than the control group (P<0.05); and the condylar long-axis in bilateral group was smaller than the control group (P<0.05); and the normal side of the condylar short-axis unilateral group was larger than that of the bilateral condylar group. Three-dimensional finite element analysis: the condyle of patients with molar scissor bite was a concentrated area of deformation during the bite of the mandible, when the first molar occlusion of the scissors bite side was simulated, the maximum deformation was located in the condyle in the X-axis and Z-axis directions. The amount of deformation was greater than that of the scissor bite side in the X-axis direction, while in the Z-axis direction, the normal side was greater than the scissor bite side. The maximum sites of local deformation in the X-axis direction were located in anterior and posterior the transverse crest of scissor bite side, and the minimum sites of local deformation was at 1/3 of the anterior slope of the inner pole of the normal side, the maximum local deformation sites in the Z-axis direction were located in the outer pole and below the outer pole of the normal side. The X-axis deformation value was the largest in the molars occlusion on the normal side, the Y-axis deformation value was in the premolars occlusion on the normal side, and the Z-axis deformation value was the largest in the centric occlusion, the deformation value of the condyle was not most significant in molar scissor bite.
CONCLUSIONS
Unilateral and bilateral molar scissor bite resulting in a short condyle morphology, and the bilateral group had a shorter condylar morphology than the unilateral group. The condyle of the patient with molar scissor bite is a concentrated area of poor occlusal deformation, and the largest sites of deformation are distributed near the transverse ridge of the inner and outer poles of the condyle. Different occlusion conditions have an effect on condylar deformation values, but do not indicate whether there is a clear association between them.
Humans
;
Finite Element Analysis
;
Retrospective Studies
;
Temporomandibular Joint/pathology*
;
Cone-Beam Computed Tomography
;
Mandible/pathology*
;
Imaging, Three-Dimensional
;
Adult
;
Temporomandibular Joint Disorders/diagnostic imaging*
;
Mandibular Condyle/diagnostic imaging*
;
Female
;
Male
;
Molar
2.Biomechanical analysis of three kinds of rigid internal fixation methods for condylar head fractures.
Junhui SUN ; Duoduo LAN ; Dong WANG ; Yao XU ; Zeyu WANG ; Chenchen ZHANG ; Kai ZHANG ; Tao XU
West China Journal of Stomatology 2025;43(1):126-132
OBJECTIVES:
This study aims to analyze the biomechanics of three kinds of rigid internal fixation methods for condylar head fractures.
METHODS:
A three dimensional finite element model of the normal mandible was constructed. It was then used to prepare condylar head fracture finite element model and three kinds of rigid internal fixation finite element model (unilateral tension screw, bilateral tension screw, tension screw+titanium plate). The mechanical characteristics and changes of the mandible condyle under the same mechanical conditions were compared among the three different rigid internal fixation methods.
RESULTS:
The maximum equivalent stress and displacement of the non-free end of condyle under the rigid internal fixation method of unilateral tension screw were 71.03 MPa and 4.72 mm, respectively. The maximum equivalent stress and displacement of the free end of condyle were 78.45 MPa and 4.50 mm, respectively. The maximum stress of fracture suture was 3.27 MPa. The maximum equivalent stress and displacement of the non-free end of condyle under the rigid internal fixation method of bilateral tension screw were 70.52 MPa and 4.00 mm, respectively. The maximum equivalent stress and displacement of the free end of condyle were 72.49 MPa and 3.85 mm, respectively. The maximum stress of fracture suture was 2.33 MPa. The maximum equivalent stress and maximum displacement of the non-free end of condyle under the rigid internal fixation method of tension screw+titanium plate were 67.26 MPa and 2.66 mm, respectively. The maximum equivalent stress and maximum displacement of the free end of condyle were 69.66 MPa and 2.50 mm, respectively. The maximum stress of fracture suture was 2.18 MPa.
CONCLUSIONS
The tension screw+titanium plate rigid internal fixation method is the most conducive to biomechanical distribution for condylar head fractures.
Fracture Fixation, Internal/instrumentation*
;
Mandibular Condyle/surgery*
;
Biomechanical Phenomena
;
Bone Screws
;
Finite Element Analysis
;
Humans
;
Mandibular Fractures/surgery*
;
Bone Plates
;
Titanium
;
Stress, Mechanical
3.Biomechanical effects of medial and lateral translation deviations of femoral components in unicompartmental knee arthroplasty on tibial prosthesis fixation.
Jingting XU ; Jing ZHANG ; Bing ZHANG ; Wen CUI ; Weijie ZHANG ; Zhenxian CHEN
Journal of Biomedical Engineering 2025;42(1):105-112
Prosthesis loosening is the leading cause of postoperative revision in unicompartmental knee arthroplasty (UKA). The deviation of medial and lateral translational installation of the prosthesis during surgery is a common clinical phenomenon and an important factor in increasing the risk of prosthesis loosening. This study established a UKA finite element model and a bone-prosthesis fixation interface micromotion prediction model. The predicted medial contact force and joint motion of the knee joint from a patient-specific lower extremity musculoskeletal multibody dynamics model of UKA were used as boundary conditions. The effects of 9 femoral component medial and lateral translational installation deviations on the Von Mises stress of the proximal tibia, the contact stress, and the micro-motion of the bone prosthesis fixation interface were quantitatively studied. It was found that compared with the neutral position (a/A of 0.492), the lateral translational deviation of the femoral component significantly increased the tibial Von Mises stress and the bone-prosthesis fixation interface contact stress. The maximum Von Mises stress and the maximum contact stress of the fixation interface increased by 14.08% and 143.15%, respectively, when a/A was 0.361. The medial translational deviation of the femoral component significantly increased the bone-prosthesis fixation interface micro-motion. The maximum value of micromotion under the conditions of femoral neutral and medial translation deviation was in the range of 20-50 μm, which is suitable for osseointegration. Therefore, based on considerations such as the micromotion range suitable for osseointegration reported in the literature, the risk of reducing prosthesis loosening, and factors that may induce pain, it is recommended that clinicians control the mounting position of the femoral component during surgery within the safe range of 0-4 mm medial translation deviation.
Humans
;
Arthroplasty, Replacement, Knee/methods*
;
Finite Element Analysis
;
Biomechanical Phenomena
;
Knee Prosthesis
;
Tibia/surgery*
;
Femur/surgery*
;
Stress, Mechanical
;
Prosthesis Failure
;
Knee Joint/surgery*
;
Prosthesis Design
4.Quantitative analysis of transcranial temporal interference stimulation in rodents: A simulation study on electrode configurations.
Xiaoxi LIU ; Hongli YU ; Fushuai GOU ; Boai DU ; Pengyi LU ; Chunfang WANG
Journal of Biomedical Engineering 2025;42(2):280-287
Transcranial temporal interference stimulation (tTIS) is a novel non-invasive transcranial electrical stimulation technique that achieves deep brain stimulation through multiple electrodes applying electric fields of different frequencies. Current studies on the mechanism of tTIS effects are primarily based on rodents, but experimental outcomes are often significantly influenced by electrode configurations. To enhance the performance of tTIS within the limited cranial space of rodents, we proposed various electrode configurations for tTIS and conducted finite element simulations using a realistic mouse model. Results demonstrated that ventral-dorsal, four-channel bipolar, and two-channel configurations performed best in terms of focality, diffusion of activated brain regions, and scalp impact, respectively. Compared to traditional transcranial direct current stimulation (tDCS), these configurations improved by 94.83%, 50.59%, and 3 514.58% in the respective evaluation metrics. This study provides a reference for selecting electrode configurations in future tTIS research on rodents.
Animals
;
Transcranial Direct Current Stimulation/instrumentation*
;
Electrodes
;
Mice
;
Computer Simulation
;
Finite Element Analysis
;
Brain/physiology*
5.A simulation study of nerve fiber activation in the lumbar segment under kilohertz-frequency transcutaneously spinal cord stimulation.
Qi XU ; Xinru LI ; Zhixin LU ; Yongchao WU
Journal of Biomedical Engineering 2025;42(2):300-307
Clinical trials have demonstrated that kilohertz-frequency transcutaneous spinal cord stimulation (TSCS) can be used to facilitate the recovery of sensory-motor function for patients with spinal cord injury, whereas the neural mechanism of TSCS is still undetermined so that the choice of stimulation parameters is largely dependent on the clinical experience. In this paper, a finite element model of transcutaneous spinal cord stimulation was used to calculate the electric field distribution of human spinal cord segments T 12 to L 2, whereas the activation thresholds of spinal fibers were determined by using a double-cable neuron model. Then the variation of activation thresholds was obtained by varying the carrier waveform, the interphase delay, the modulating frequency, and the modulating pulse width. Compared with the sinusoidal carrier, the usage of square carrier could significantly reduce the activation threshold of dorsal root (DR) fibers. Moreover, the variation of activation thresholds was no more than 1 V due to the varied modulating frequency and decreases with the increased modulating pulse width. For a square carrier at 10 kHz modulated by rectangular pulse with the frequency of 50 Hz and the pulse width of 1 ms, the lowest activation thresholds of DR fibers and dorsal column fibers were 27.6 V and 55.8 V, respectively. An interphase delay of 5 μs was able to reduce the activation thresholds of the DR fibers to 20.1 V. The simulation results can lay a theoretical foundation on the selection of TSCS parameters in clinical trials.
Humans
;
Spinal Cord Stimulation/methods*
;
Nerve Fibers/physiology*
;
Finite Element Analysis
;
Spinal Cord/physiology*
;
Computer Simulation
;
Spinal Cord Injuries/physiopathology*
;
Lumbosacral Region
;
Lumbar Vertebrae
;
Transcutaneous Electric Nerve Stimulation/methods*
;
Models, Neurological
6.Biomechanical study on wing shaped titanium plate fixation of acetabular anterior column and posterior hemi-transverse fracture under multiple working conditions.
Jianwu ZHANG ; WURIKAIXI AIYITI ; Gang LYU ; MAIMAIAILI YUSHAN ; Zhiqiang MA ; Chao MA
Journal of Biomedical Engineering 2025;42(2):351-358
This article aims to compare and analyze the biomechanical differences between wing-shaped titanium plates and traditional titanium plates in fixing acetabular anterior column and posterior hemi-transverse (ACPHT) fracture under multiple working conditions using the finite element method. Firstly, four sets of internal fixation models for acetabular ACPHT fractures were established, and the hip joint stress under standing, sitting, forward extension, and abduction conditions was calculated through analysis software. Then, the stress of screws and titanium plates, as well as the stress and displacement of the fracture end face, were analyzed. Research has found that when using wing-shaped titanium plates to fix acetabular ACPHT fractures, the peak stress of screws decreases under all working conditions, while the peak stress of wing-shaped titanium plates decreases under standing and sitting conditions and increases under forward and outward extension conditions. The relative displacement and mean stress of the fracture end face decrease under all working conditions, but the values are higher under forward and outward extension conditions. Wing-shaped titanium plates can reduce the probability of screw fatigue failure when fixing acetabular ACPHT fractures and can bear greater loads under forward and outward extension conditions, improving the mechanical stability of the pelvis. Moreover, the stress on the fracture end surface is more conducive to stimulating fracture healing and promoting bone tissue growth. However, premature forward and outward extension rehabilitation exercises should not be performed.
Titanium
;
Bone Plates
;
Humans
;
Acetabulum/surgery*
;
Fracture Fixation, Internal/methods*
;
Biomechanical Phenomena
;
Finite Element Analysis
;
Bone Screws
;
Fractures, Bone/surgery*
;
Stress, Mechanical
;
Working Conditions
7.A study on the predictive model of porous hyperelastic properties of human alveolar bone based on computed tomography imaging.
Bin WU ; Mingna LI ; Fan YANG ; Le YUAN ; Yi LU ; Di JIANG ; Yang YI ; Bin YAN
Journal of Biomedical Engineering 2025;42(2):359-365
Alveolar bone reconstruction simulation is an effective means for quantifying orthodontics, but currently, it is not possible to directly obtain human alveolar bone material models for simulation. This study introduces a prediction method for the equivalent shear modulus of three-dimensional random porous materials, integrating the first-order Ogden hyperelastic model to construct a computed tomography (CT) based porous hyperelastic Ogden model (CT-PHO) for human alveolar bone. Model parameters are derived by combining results from micro-CT, nanoindentation experiments, and uniaxial compression tests. Compared to previous predictive models, the CT-PHO model shows a lower root mean square error (RMSE) under all bone density conditions. Simulation results using the CT-PHO model parameters in uniaxial compression experiments demonstrate more accurate prediction of the mechanical behavior of alveolar bone under compression. Further prediction and validation with different individual human alveolar bone samples yield accurate results, confirming the generality of the CT-PHO model. The study suggests that the CT-PHO model proposed in this paper can estimate the material properties of human alveolar bone and may eventually be used for bone reconstruction simulations to guide clinical treatment.
Humans
;
Tomography, X-Ray Computed/methods*
;
Porosity
;
Alveolar Process/physiology*
;
Bone Density
;
Computer Simulation
;
Elasticity
;
X-Ray Microtomography
;
Stress, Mechanical
;
Finite Element Analysis
;
Models, Biological
8.Modeling and finite element analysis of human trabecular meshwork outflow pathways.
Shiya BAO ; Qing SUN ; Si CHEN ; Xinyu CHEN ; Xiang PENG ; Jing ZHANG
Journal of Biomedical Engineering 2025;42(3):585-591
Glaucoma is the leading cause of irreversible blindness worldwide, with its primary risk factor arising from elevated intraocular pressure (IOP) due to an imbalance between aqueous humor production and outflow. This study aims to establish quantitative correlations among IOP, iris mechanical properties, channel microstructures, and aqueous humor dynamics through three-dimensional modeling and finite element analysis, overcoming the limitations of conventional experimental techniques in studying aqueous flow within the trabecular meshwork (TM) outflow pathway. A three-dimensional fluid-structure interaction (FSI) model incorporating the layered TM structure, Schlemm's canal (SC), iris, and other anterior segment tissues was developed based on human ocular anatomy. FSI simulations were performed to quantify the effects of IOP variations and iris Young's modulus on tissue morphology and aqueous humor dynamics parameters. The computational results demonstrated that axial iris deformation showed significant correlations with IOP and iris Young's modulus. Although elevated IOP exhibited minimal effects on hydrodynamic parameters in the anterior and posterior chambers, it markedly suppressed aqueous flow velocity in the TM region. Additionally, wall shear stress in SC and collector channels displayed high sensitivity to IOP variations. These findings reveal that the tissue mechanics-FSI mechanism modulates outflow resistance by regulating aqueous humor dynamics, offering valuable references for developing clinical therapies targeting IOP reduction in glaucoma management.
Humans
;
Trabecular Meshwork/anatomy & histology*
;
Finite Element Analysis
;
Aqueous Humor/metabolism*
;
Intraocular Pressure/physiology*
;
Glaucoma/physiopathology*
;
Iris/anatomy & histology*
;
Computer Simulation
;
Models, Biological
9.Development of a balloon sensor device for force-electrical coupling measurement of esophagus.
Peng RAN ; Ying ZHONG ; Yingbing LAI ; Lei LIU ; Yanhang ZHU ; Huantao ZHU
Journal of Biomedical Engineering 2025;42(3):610-619
To address the challenges of capturing micro-strains in detecting esophageal motility disorders and the limitations of existing high-resolution manometry and functional intraluminal imaging probes in directly measuring esophageal tissue electrical impedance, this study proposes a novel flexible balloon sensor structure that integrates a piezoelectric film assembly with a distributed impedance electrode array. Using the electrical analysis module in the finite element analysis (FEA) software, simulations of the forward problem for esophageal impedance detection were conducted to optimize the excitation source parameters, and a physical prototype was fabricated. Under a relative excitation mode with a voltage sensitivity of 2.059%, the voltage output characteristics of the impedance electrode array were analyzed during linear changes in the balloon filling volume. Based on the performance variation of the piezoelectric film assembly, 80% was selected as the optimal filling volume. Force-electric coupling tests were conducted on the balloon sensor using a pressure testing platform, revealing that both the piezoelectric film assembly inside the balloon and the impedance electrodes outside the balloon exhibited significant load differentiation characteristics as the force application point shifted. In summary, this balloon sensor facilitates the localization of force application while simultaneously analyzing esophageal tissue properties, offering a novel diagnostic approach and objective tool for esophageal disease detection.
Esophagus/physiology*
;
Electric Impedance
;
Humans
;
Finite Element Analysis
;
Manometry/methods*
;
Electrodes
;
Esophageal Motility Disorders/physiopathology*
;
Equipment Design
10.Effects of elastic modulus of the metal block on the condylar-constrained knee prosthesis tibial fixation stability.
Yuhan ZHANG ; Jing ZHANG ; Tianqi DONG ; Xuan ZHANG ; Weijie ZHANG ; Lei GUO ; Zhenxian CHEN
Journal of Biomedical Engineering 2025;42(4):782-789
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.
Knee Prosthesis
;
Humans
;
Finite Element Analysis
;
Tibia/surgery*
;
Elastic Modulus
;
Arthroplasty, Replacement, Knee/methods*
;
Stress, Mechanical
;
Metals
;
Prosthesis Design
;
Knee Joint/surgery*
;
Biomechanical Phenomena

Result Analysis
Print
Save
E-mail