Bioengineering

Research

9 pages, 2729 KiB

Open AccessArticle

Determining the Relationship between Mechanical Properties and Quantitative Magnetic Resonance Imaging of Joint Soft Tissues Using Patient-Specific Templates

by Takehito Hananouchi, Shinji Satake, Kei Sakao, Hiroshi Katsuda, Nagakazu Shimada, Erik W. Dorthe and Darryl D. D’Lima

Bioengineering 2023, 10(9), 1050; https://doi.org/10.3390/bioengineering10091050 - 7 Sep 2023

Viewed by 1077

Abstract

To determine whether the mechanical properties of joint soft tissues such as cartilage can be calculated from quantitative magnetic resonance imaging (MRI) data, we investigated whether the mechanical properties of articular cartilage and meniscus scheduled to be resected during arthroplasty are correlated with [...] Read more.

To determine whether the mechanical properties of joint soft tissues such as cartilage can be calculated from quantitative magnetic resonance imaging (MRI) data, we investigated whether the mechanical properties of articular cartilage and meniscus scheduled to be resected during arthroplasty are correlated with the T2 relaxation time on quantitative MRI at the same location. Six patients who had undergone knee arthroplasty and seven who had undergone hip arthroplasty were examined. For the knee joint, the articular cartilage and lateral meniscus of the distal lateral condyle of the femur and proximal lateral tibia were examined, while for the hip joint, the articular cartilage above the femoral head was studied. We investigated the relationship between T2 relaxation time by quantitative MRI and stiffness using a hand-made compression tester at 235 locations. The patient-individualized template technique was used to align the two measurement sites. The results showed a negative correlation (from −0.30 to −0.35) in the less severely damaged articular cartilage and meniscus. This indicates that tissue mechanical properties can be calculated from T2 relaxation time, suggesting that quantitative MRI is useful in determining when to start loading after interventional surgery on cartilage tissue and in managing the health of elderly patients. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

12 pages, 3339 KiB

Open AccessArticle

3D Tortuosity and Diffusion Characterization in the Human Mineralized Collagen Fibril Using a Random Walk Model

by Fabiano Bini, Andrada Pica, Andrea Marinozzi and Franco Marinozzi

Bioengineering 2023, 10(5), 558; https://doi.org/10.3390/bioengineering10050558 - 7 May 2023

Viewed by 1695

Abstract

Bone tissue is mainly composed at the nanoscale of apatite minerals, collagen molecules and water that form the mineralized collagen fibril (MCF). In this work, we developed a 3D random walk model to investigate the influence of bone nanostructure on water diffusion. We [...] Read more.

Bone tissue is mainly composed at the nanoscale of apatite minerals, collagen molecules and water that form the mineralized collagen fibril (MCF). In this work, we developed a 3D random walk model to investigate the influence of bone nanostructure on water diffusion. We computed 1000 random walk trajectories of water molecules within the MCF geometric model. An important parameter to analyse transport behaviour in porous media is tortuosity, computed as the ratio between the effective path length and the straight-line distance between initial and final points. The diffusion coefficient is determined from the linear fit of the mean squared displacement of water molecules as a function of time. To achieve more insight into the diffusion phenomenon within MCF, we estimated the tortuosity and diffusivity at different quotes in the longitudinal direction of the model. Tortuosity is characterized by increasing values in the longitudinal direction. As expected, the diffusion coefficient decreases as tortuosity increases. Diffusivity outcomes confirm the findings achieved by experimental investigations. The computational model provides insights into the relation between the MCF structure and mass transport behaviour that may contribute to the improvement of bone-mimicking scaffolds. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

14 pages, 3642 KiB

Open AccessArticle

Finite Element Modeling of Residual Hearing after Cochlear Implant Surgery in Chinchillas

by Nicholas Castle, Junfeng Liang, Matthew Smith, Brett Petersen, Cayman Matson, Tara Eldridge, Ke Zhang, Chung-Hao Lee, Yingtao Liu and Chenkai Dai

Bioengineering 2023, 10(5), 539; https://doi.org/10.3390/bioengineering10050539 - 27 Apr 2023

Cited by 1 | Viewed by 1658

Abstract

Cochlear implant (CI) surgery is one of the most utilized treatments for severe hearing loss. However, the effects of a successful scala tympani insertion on the mechanics of hearing are not yet fully understood. This paper presents a finite element (FE) model of [...] Read more.

Cochlear implant (CI) surgery is one of the most utilized treatments for severe hearing loss. However, the effects of a successful scala tympani insertion on the mechanics of hearing are not yet fully understood. This paper presents a finite element (FE) model of the chinchilla inner ear for studying the interrelationship between the mechanical function and the insertion angle of a CI electrode. This FE model includes a three-chambered cochlea and full vestibular system, accomplished using µ-MRI and µ-CT scanning technologies. This model’s first application found minimal loss of residual hearing due to insertion angle after CI surgery, and this indicates that it is a reliable and helpful tool for future applications in CI design, surgical planning, and stimuli setup. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

18 pages, 3755 KiB

Open AccessArticle

Multi-Objective Optimisation of a Novel Bypass Graft with a Spiral Ridge

by Antonios Xenakis, Andres Ruiz-Soler and Amir Keshmiri

Bioengineering 2023, 10(4), 489; https://doi.org/10.3390/bioengineering10040489 - 19 Apr 2023

Cited by 3 | Viewed by 1475

Abstract

The low long-term patency of bypass grafts is a major concern for cardiovascular treatments. Unfavourable haemodynamic conditions in the proximity of distal anastomosis are closely related to thrombus creation and lumen lesions. Modern graft designs address this unfavourable haemodynamic environment with the introduction [...] Read more.

The low long-term patency of bypass grafts is a major concern for cardiovascular treatments. Unfavourable haemodynamic conditions in the proximity of distal anastomosis are closely related to thrombus creation and lumen lesions. Modern graft designs address this unfavourable haemodynamic environment with the introduction of a helical component in the flow field, either by means of out-of-plane helicity graft geometry or a spiral ridge. While the latter has been found to lack in performance when compared to the out-of-plane helicity designs, recent findings support the idea that the existing spiral ridge grafts can be further improved in performance through optimising relevant design parameters. In the current study, robust multi-objective optimisation techniques are implemented, covering a wide range of possible designs coupled with proven and well validated computational fluid dynamics (CFD) algorithms. It is shown that the final set of suggested design parameters could significantly improve haemodynamic performance and therefore could be used to enhance the design of spiral ridge bypass grafts. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

13 pages, 4597 KiB

Open AccessArticle

Biomechanical Effects of a Novel Pedicle Screw W-Type Rod Fixation for Lumbar Spondylolysis: A Finite Element Analysis

by Jo-Hsi Pan, Chen-Sheng Chen, Chien-Lin Liu and Po-Hsin Chou

Bioengineering 2023, 10(4), 451; https://doi.org/10.3390/bioengineering10040451 - 7 Apr 2023

Cited by 2 | Viewed by 1745

Abstract

Lumbar spondylolysis involves anatomical defects of the pars interarticularis, which causes instability during motion. The instability can be addressed through instrumentation with posterolateral fusion (PLF). We developed a novel pedicle screw W-type rod fixation system and evaluated its biomechanical effects in comparison with [...] Read more.

Lumbar spondylolysis involves anatomical defects of the pars interarticularis, which causes instability during motion. The instability can be addressed through instrumentation with posterolateral fusion (PLF). We developed a novel pedicle screw W-type rod fixation system and evaluated its biomechanical effects in comparison with PLF and Dynesys stabilization for lumbar spondylolysis via finite element (FE) analysis. A validated lumbar spine model was built using ANSYS 14.5 software. Five FE models were established simulating the intact L1–L5 lumbar spine (INT), bilateral pars defect (Bipars), bilateral pars defect with PLF (Bipars_PLF), Dynesys stabilization (Bipars_Dyn), and W-type rod fixation (Bipars_Wtyp). The range of motion (ROM) of the affected segment, the disc stress (DS), and the facet contact force (FCF) of the cranial segment were compared. In the Bipars model, ROM increased in extension and rotation. Compared with the INT model, Bipars_PLF and Bipars_Dyn exhibited remarkably lower ROMs for the affected segment and imposed greater DS and FCF in the cranial segment. Bipars_Wtyp preserved more ROM and generated lower stress at the cranial segment than Bipars_PLF or Bipars_Dyn. The injury model indicates that this novel pedicle screw W-type rod for spondylolysis fixation could return ROM, DS, and FCF to levels similar to preinjury. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

12 pages, 6396 KiB

Open AccessArticle

A Symmetry-Based Superposition Method for Planning and Surgical Outcome Assessment

by Yu-Ching Hsiao and Jing-Jing Fang

Bioengineering 2023, 10(3), 335; https://doi.org/10.3390/bioengineering10030335 - 6 Mar 2023

Cited by 1 | Viewed by 1388

Abstract

Computer-aided surgical planning has been widely used to increase the safety and predictability of surgery. The validation of the target of surgical planning to surgical outcomes on a patient-specific model is an important issue. The aim of this research was to develop a [...] Read more.

Computer-aided surgical planning has been widely used to increase the safety and predictability of surgery. The validation of the target of surgical planning to surgical outcomes on a patient-specific model is an important issue. The aim of this research was to develop a robust superposition method to assess the deviation of planning and outcome by using the symmetrical characteristic of the affected target. The optimal symmetry plane (OSP) of an object is usually used to evaluate the degree of symmetry of an object. We proposed a refined OSP-based contouring method to transfer a complex three-dimensional superposition operation into two dimensions. We compared the typical iterative closest point (ICP) algorithm with the refined OSP-based contouring method and examined the differences between them. The results using the OSP-based method were much better than the traditional method. As for processing time, the OSP-based contouring method was 11 times faster than the ICP method overall. The proposed method was not affected by the metallic artifacts from medical imaging or geometric changes due to surgical intervention. This technique can be applied for post-operative assessment, such as quantifying the differences between surgical targets and outcomes as well as performing long-term medical follow-up. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

13 pages, 47857 KiB

Open AccessArticle

Modeling the Impact of Meniscal Tears on von Mises Stress of Knee Cartilage Tissue

by Oleg Ardatov, Viktorija Aleksiuk, Algirdas Maknickas, Rimantas Stonkus, Ilona Uzieliene, Raminta Vaiciuleviciute, Jolita Pachaleva, Giedrius Kvederas and Eiva Bernotiene

Bioengineering 2023, 10(3), 314; https://doi.org/10.3390/bioengineering10030314 - 1 Mar 2023

Cited by 5 | Viewed by 2312

Abstract

The present study aims to explore the stressed state of cartilage using various meniscal tear models. To perform this research, the anatomical model of the knee joint was developed and the nonlinear mechanical properties of the cartilage and meniscus were verified. The stress–strain [...] Read more.

The present study aims to explore the stressed state of cartilage using various meniscal tear models. To perform this research, the anatomical model of the knee joint was developed and the nonlinear mechanical properties of the cartilage and meniscus were verified. The stress–strain curve of the meniscus was obtained by testing fresh tissue specimens of the human meniscus using a compression machine. The results showed that the more deteriorated meniscus had greater stiffness, but its integrity had the greatest impact on the growth of cartilage stresses. To confirm this, cases of radial, longitudinal, and complex tears were examined. The methodology and results of the study can assist in medical diagnostics for meniscus treatment and replacement. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Graphical abstract

16 pages, 6055 KiB

Open AccessArticle

Hemodynamic Effects of Subaortic Stenosis on Blood Flow Characteristics of a Mechanical Heart Valve Based on OpenFOAM Simulation

by Aolin Chen, Adi Azriff Basri, Norzian Bin Ismail and Kamarul Arifin Ahmad

Bioengineering 2023, 10(3), 312; https://doi.org/10.3390/bioengineering10030312 - 1 Mar 2023

Cited by 2 | Viewed by 1944 | Correction

Abstract

Subaortic stenosis (SAS) is a common congenital heart disease that can cause significant morbidity and mortality if not treated promptly. Patients with heart valve disease are prone to complications after replacement surgery, and the existence of SAS can accelerates disease progression, so timely [...] Read more.

Subaortic stenosis (SAS) is a common congenital heart disease that can cause significant morbidity and mortality if not treated promptly. Patients with heart valve disease are prone to complications after replacement surgery, and the existence of SAS can accelerates disease progression, so timely diagnosis and treatment are required. However, the effects of subaortic stenosis on mechanical heart valves (MHV) are unknown. This study aimed to investigate flow characteristics in the presence of subaortic stenosis and computationally quantify the effects on the hemodynamics of MHV. Through the numerical simulation method, the flow characteristics and related parameters in the presence of SAS can be more intuitively observed. Based on its structure, there are three types of SAS: Tunnel-type SAS (TSS); Fibromuscular annulus SAS (FSS); Discrete SAS (DSS). The first numerical simulation study on different types of SAS found that there are obvious differences among them. Among them, the tunnel-type SAS formed a separated vortex structure on the tunnel-type narrow surface, which exhibits higher wall shear force at a low obstacle percentage. However, discrete SAS showed obvious differences when there was a high percentage of obstacles, forming high peak flow, high wall shear stress, and a high-intensity complex vortex. The presence of all three types of SAS results in the formation of high-velocity jets and complex vortices in front of the MHV, leading to increased shear stress and stagnation time. These hemodynamic changes significantly increase the risk of MHV dysfunction and the development of complications. Despite differences between the three types of SAS, the resultant effects on MHV hemodynamics are consistent. Therefore, early surgical intervention is warranted in SAS patients with implanted MHV. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

20 pages, 17249 KiB

Open AccessArticle

Development of a Hip Joint Socket by Finite-Element-Based Analysis for Mechanical Assessment

by Ana Karen González, Juvenal Rodríguez-Reséndiz, José Eli Eduardo Gonzalez-Durán, Juan Manuel Olivares Ramírez and Adyr A. Estévez-Bén

Bioengineering 2023, 10(2), 268; https://doi.org/10.3390/bioengineering10020268 - 18 Feb 2023

Cited by 1 | Viewed by 1914

Abstract

This article evaluates a hip joint socket design by finite element method (FEM). The study was based on the needs and characteristics of a patient with an oncological amputation; however, the solution and the presented method may be generalized for patients with similar [...] Read more.

This article evaluates a hip joint socket design by finite element method (FEM). The study was based on the needs and characteristics of a patient with an oncological amputation; however, the solution and the presented method may be generalized for patients with similar conditions. The research aimed to solve a generalized problem, taking a typical case from the study area as a reference. Data were collected on the use of the current improving prosthesis—specifically in interaction with its socket—to obtain information on the new approach design: this step constituted the work’s starting point, where the problems to be solved in conventional designs were revealed. Currently, the development of this type of support does not consider the functionality and comfort of the patient. Research has reported that 58% of patients with sockets have rejected their use, because they do not fit comfortably and functionally; therefore, patients’ low acceptance or rejection of the use of the prosthesis socket has been documented. In this study, different designs were evaluated, based on the FEM as scientific support for the results obtained, for the development of a new ergonomic fit with a 60% increase in patient compliance, that had correct gait performance when correcting postures, improved fit–user interaction, and that presented an esthetic fit that met the usability factor. The validation of the results was carried out through the physical construction of the prototype. The research showed how the finite element method improved the design, analyzing the structural behavioral, and that it could reduce cost and time instead of generating several prototypes. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

15 pages, 4377 KiB

Open AccessArticle

Three-Dimensional Computational Model Simulating the Initial Callus Growth during Fracture Healing in Long Bones: Application to Different Fracture Types

by José M. Naveiro, Luis Gracia, Jorge Roces, Jorge Albareda and Sergio Puértolas

Bioengineering 2023, 10(2), 190; https://doi.org/10.3390/bioengineering10020190 - 2 Feb 2023

Cited by 1 | Viewed by 2231

Abstract

Bone fractures are among the most common and potentially serious injuries to the skeleton, femoral shaft fractures being especially severe. Thanks to recent advances in the area of in silico analysis, several approximations of the bone healing process have been achieved. In this [...] Read more.

Bone fractures are among the most common and potentially serious injuries to the skeleton, femoral shaft fractures being especially severe. Thanks to recent advances in the area of in silico analysis, several approximations of the bone healing process have been achieved. In this context, the objective of this work was to simulate the initial phase of callus formation in long bones, without a pre-meshed domain in the 3D space. A finite element approach was computationally implemented to obtain the values of the cell concentrations along the whole domain and evaluate the areas where the biological quantities reached the thresholds necessary to trigger callus growth. A voxel model was used to obtain the 3D domain of the bone fragments and callus. A mesh growth algorithm controlled the addition of new elements to the domain at each step of the iterative procedure until complete callus formation. The implemented approach is able to reproduce the generation of the primary callus, which corresponds to the initial phase of fracture healing, independently of the fracture type and complexity, even in the case of several bone fragments. The proposed approach can be applied to the most complex bone fractures such as oblique, severely comminuted or spiral-type fractures, whose simulation remains hardly possible by means of the different existing approaches available to date. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Graphical abstract

12 pages, 8773 KiB

Open AccessArticle

The Role of Multifidus in the Biomechanics of Lumbar Spine: A Musculoskeletal Modeling Study

by Kuan Wang, Zhen Deng, Xinpeng Chen, Jiang Shao, Lulu Qiu, Chenghua Jiang and Wenxin Niu

Bioengineering 2023, 10(1), 67; https://doi.org/10.3390/bioengineering10010067 - 4 Jan 2023

Cited by 22 | Viewed by 3458

Abstract

Background: The role of multifidus in the biomechanics of lumbar spine remained unclear. Purpose: This study aimed to investigate the role of multifidus in the modeling of lumbar spine and the influence of asymmetric multifidus atrophy on the biomechanics of lumbar spine. Methods: [...] Read more.

Background: The role of multifidus in the biomechanics of lumbar spine remained unclear. Purpose: This study aimed to investigate the role of multifidus in the modeling of lumbar spine and the influence of asymmetric multifidus atrophy on the biomechanics of lumbar spine. Methods: This study considered five different multifidus conditions in the trunk musculoskeletal models: group 1 (with entire multifidus), group 2 (without multifidus), group 3 (multifidus with half of maximum isometric force), group 4 (asymmetric multifidus atrophy on L5/S1 level), and group 5 (asymmetric multifidus atrophy on L4/L5 level). In order to test how different multifidus situations would affect the lumbar spine, four trunk flexional angles (0°, 30°, 60°, and 90°) were simulated. The calculation of muscle activation and muscle force was done using static optimization function in OpenSim. Then, joint reaction forces of L5/S1 and L4/L5 levels were calculated and compared among the groups. Results: The models without multifidus had the highest normalized compressive forces on the L4/L5 level in trunk flexion tasks. In extreme cases produced by group 2 models, the normalized compressive forces on L4/L5 level were 444% (30° flexion), 568% (60° flexion), and 576% (90° flexion) of upper body weight, which were 1.82 times, 1.63 times, and 1.13 times as large as the values computed by the corresponding models in group 1. In 90° flexion, the success rate of simulation in group 2 was 49.6%, followed by group 3 (84.4%), group 4 (89.6%), group 5 (92.8%), and group 1 (92.8%). Conclusions: The results demonstrate that incorporating multifidus in the musculoskeletal model is important for increasing the success rate of simulation and decreasing the incidence of overestimation of compressive load on the lumbar spine. Asymmetric multifidus atrophy has negligible effect on the lower lumbar spine in the trunk flexion posture. The results highlighted the fine-tuning ability of multifidus in equilibrating the loads on the lower back and the necessity of incorporating multifidus in trunk musculoskeletal modeling. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Graphical abstract

19 pages, 6673 KiB

Open AccessArticle

Hierarchical Structure and Properties of the Bone at Nano Level

by Farah Hamandi and Tarun Goswami

Bioengineering 2022, 9(11), 677; https://doi.org/10.3390/bioengineering9110677 - 10 Nov 2022

Cited by 4 | Viewed by 3360

Abstract

Bone is a highly hierarchical complex structure that consists of organic and mineral components represented by collagen molecules (CM) and hydroxyapatite crystals (HAC), respectively. The nanostructure of bone can significantly affect its mechanical properties. There is a lack of understanding how collagen fibrils [...] Read more.

Bone is a highly hierarchical complex structure that consists of organic and mineral components represented by collagen molecules (CM) and hydroxyapatite crystals (HAC), respectively. The nanostructure of bone can significantly affect its mechanical properties. There is a lack of understanding how collagen fibrils (CF) in different orientations may affect the mechanical properties of the bone. The objective of this study is to investigate the effect of interaction, orientation, and hydration on atomic models of the bone composed of collagen helix (CH) and HAC, using molecular dynamics simulations and therefrom bone-related disease origins. The results demonstrate that the mechanical properties of the bone are affected significantly by the orientation of the CF attributed to contact areas at 0° and 90° models. The molecular dynamics simulation illustrated that there is significant difference (p < 0.005) in the ultimate tensile strength and toughness with respect to the orientation of the hydrated and un-hydrated CF. Additionally, the results indicated that having the force in a longitudinal direction (0°) provides more strength compared with the CF in the perpendicular direction (90°). Furthermore, the results show that substituting glycine (GLY) with any other amino acid affects the mechanical properties and strength of the CH, collagen–hydroxyapatite interface, and eventually affects the HAC. Generally, hydration dramatically influences bone tissue elastic properties, and any change in the orientation or any abnormality in the atomic structure of either the CM or the HAC would be the main reason of the fragility in the bone, affecting bone pathology. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

17 pages, 7716 KiB

Open AccessArticle

The Effect of Intraocular Pressure Load Boundary on the Biomechanics of the Human Conventional Aqueous Outflow Pathway

by Alireza Karimi, Reza Razaghi, Seyed Mohammadali Rahmati, J. Crawford Downs, Ted S. Acott, Mary J. Kelley, Ruikang K. Wang and Murray Johnstone

Bioengineering 2022, 9(11), 672; https://doi.org/10.3390/bioengineering9110672 - 10 Nov 2022

Cited by 10 | Viewed by 2205

Abstract

Background: Aqueous humor outflow resistance in the trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and Schlemm’s canal (SC) endothelium of the conventional outflow pathway actively contribute to intraocular pressure (IOP) regulation. Outflow resistance is actively affected by the dynamic outflow pressure gradient across [...] Read more.

Background: Aqueous humor outflow resistance in the trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and Schlemm’s canal (SC) endothelium of the conventional outflow pathway actively contribute to intraocular pressure (IOP) regulation. Outflow resistance is actively affected by the dynamic outflow pressure gradient across the TM, JCT, and SC inner wall tissues. The resistance effect implies the presence of a fluid–structure interaction (FSI) coupling between the outflow tissues and the aqueous humor. However, the biomechanical interactions between viscoelastic outflow tissues and aqueous humor dynamics are largely unknown. Methods: A 3D microstructural finite element (FE) model of a healthy human eye TM/JCT/SC complex was constructed with elastic and viscoelastic material properties for the bulk extracellular matrix and embedded elastic cable elements. The FE models were subjected to both idealized and a physiologic IOP load boundary using the FSI method. Results: The elastic material model for both the idealized and physiologic IOP load boundary at equal IOPs showed similar stresses and strains in the outflow tissues as well as pressure in the aqueous humor. However, outflow tissues with viscoelastic material properties were sensitive to the IOP load rate, resulting in different mechanical and hydrodynamic responses in the tissues and aqueous humor. Conclusions: Transient IOP fluctuations may cause a relatively large IOP difference of ~20 mmHg in a very short time frame of ~0.1 s, resulting in a rate stiffening in the outflow tissues. Rate stiffening reduces strains and causes a rate-dependent pressure gradient across the outflow tissues. Thus, the results suggest it is necessary to use a viscoelastic material model in outflow tissues that includes the important role of IOP load rate. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

13 pages, 29970 KiB

Open AccessArticle

Topology Optimization-Based Localized Bone Microstructure Reconstruction for Image Resolution Enhancement: Accuracy and Efficiency

by Jisun Kim and Jung Jin Kim

Bioengineering 2022, 9(11), 644; https://doi.org/10.3390/bioengineering9110644 - 3 Nov 2022

Cited by 2 | Viewed by 1892

Abstract

Topology optimization is currently the only way to provide bone microstructure information by enhancing a 600 μm low-resolution image into a 50 μm high-resolution image. Particularly, the recently proposed localized reconstruction method for the region of interest has received much attention because it [...] Read more.

Topology optimization is currently the only way to provide bone microstructure information by enhancing a 600 μm low-resolution image into a 50 μm high-resolution image. Particularly, the recently proposed localized reconstruction method for the region of interest has received much attention because it has a high possibility to overcome inefficiency such as iterative large-scale problems of the conventional reconstruction. Despite the great potential, the localized method should be thoroughly validated for clinical application. This study aims to quantitatively validate the topology optimization-based localized bone microstructure reconstruction method in terms of accuracy and efficiency by comparing the conventional method. For this purpose, this study re-constructed bone microstructure for three regions of interest in the proximal femur by localized and conventional methods, respectively. In the comparison, the dramatically reduced total progress time by at least 88.2% (20.1 h) as well as computational resources by more than 95.9% (54.0 gigabytes) were found. Moreover, very high reconstruction accuracy in the trabecular alignment (up to 99.6%) and morphometric indices (up to 2.71%) was also found. These results indicated that the localized method could reconstruct bone microstructure, much more effectively preserving the originality of the conventional method. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Graphical abstract

11 pages, 3071 KiB

Open AccessArticle

Sensitivity Analysis of Cardiac Alternans and Tachyarrhythmia to Ion Channel Conductance Using Population Modeling

by Da Un Jeong, Aroli Marcellinus and Ki Moo Lim

Bioengineering 2022, 9(11), 628; https://doi.org/10.3390/bioengineering9110628 - 1 Nov 2022

Cited by 1 | Viewed by 1686

Abstract

Action potential duration (APD) alternans, an alternating phenomenon between action potentials in cardiomyocytes, causes heart arrhythmia when the heart rate is high. However, some of the APD alternans observed in clinical trials occurs under slow heart rate conditions of 100 to 120 bpm, [...] Read more.

Action potential duration (APD) alternans, an alternating phenomenon between action potentials in cardiomyocytes, causes heart arrhythmia when the heart rate is high. However, some of the APD alternans observed in clinical trials occurs under slow heart rate conditions of 100 to 120 bpm, increasing the likelihood of heart arrhythmias such as atrial fibrillation. Advanced studies have identified the occurrence of this type of APD alternans in terms of electrophysiological ion channel currents in cells. However, they only identified physiological phenomena, such as action potential due to random changes in a particular ion channel’s conductivity through ion models specializing in specific ion channel currents. In this study, we performed parameter sensitivity analysis via population modeling using a validated human ventricular physiology model to check the sensitivity of APD alternans to ion channel conductances. Through population modeling, we expressed the changes in alternans onset cycle length (AOCL) and mean APD in AOCL (AO meanAPD) according to the variations in ion channel conductance. Finally, we identified the ion channel that maximally affected the occurrence of APD alternans. AOCL and AO meanAPD were sensitive to changes in the plateau Ca²⁺ current. Accordingly, it was expected that APD alternans would be vulnerable to changes in intracellular calcium concentration. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Graphical abstract

49 pages, 19903 KiB

Open AccessArticle

Traumatic Brain Injury Biomarkers, Simulations and Kinetics

by Celeste Hicks, Akshima Dhiman, Chauntel Barrymore and Tarun Goswami

Bioengineering 2022, 9(11), 612; https://doi.org/10.3390/bioengineering9110612 - 25 Oct 2022

Cited by 10 | Viewed by 3452

Abstract

This paper reviews the predictive capabilities of blood-based biomarkers to quantify traumatic brain injury (TBI). Biomarkers for concussive conditions also known as mild, to moderate and severe TBI identified along with post-traumatic stress disorder (PTSD) and chronic traumatic encephalopathy (CTE) that occur due [...] Read more.

This paper reviews the predictive capabilities of blood-based biomarkers to quantify traumatic brain injury (TBI). Biomarkers for concussive conditions also known as mild, to moderate and severe TBI identified along with post-traumatic stress disorder (PTSD) and chronic traumatic encephalopathy (CTE) that occur due to repeated blows to the head during one’s lifetime. Since the pathways of these biomarkers into the blood are not fully understood whether there is disruption in the blood–brain barrier (BBB) and the time it takes after injury for the expression of the biomarkers to be able to predict the injury effectively, there is a need to understand the protein biomarker structure and other physical properties. The injury events in terms of brain and mechanics are a result of external force with or without the shrapnel, in the wake of a wave result in local tissue damage. Thus, these mechanisms express specific biomarkers kinetics of which reaches half-life within a few hours after injury to few days. Therefore, there is a need to determine the concentration levels that follow injury. Even though current diagnostics linking biomarkers with TBI severity are not fully developed, there is a need to quantify protein structures and their viability after injury. This research was conducted to fully understand the structures of 12 biomarkers by performing molecular dynamics simulations involving atomic movement and energies of forming hydrogen bonds. Molecular dynamics software, NAMD and VMD were used to determine and compare the approximate thermodynamic stabilities of the biomarkers and their bonding energies. Five biomarkers used clinically were S100B, GFAP, UCHL1, NF-L and tau, the kinetics obtained from literature show that the concentration values abruptly change with time after injury. For a given protein length, associated number of hydrogen bonds and bond energy describe a lower bound region where proteins self-dissolve and do not have long enough half-life to be detected in the fluids. However, above this lower bound, involving higher number of bonds and energy, we hypothesize that biomarkers will be viable to disrupt the BBB and stay longer to be modeled for kinetics for diagnosis and therefore may help in the discoveries of new biomarkers. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

11 pages, 1797 KiB

Open AccessArticle

Finite Element Analysis and Experimental Validation of the Anterior Cruciate Ligament and Implications for the Injury Mechanism

by Shuang Ren, Huijuan Shi, Zhenlong Liu, Jiahao Zhang, Hanjun Li, Hongshi Huang and Yingfang Ao

Bioengineering 2022, 9(10), 590; https://doi.org/10.3390/bioengineering9100590 - 21 Oct 2022

Cited by 12 | Viewed by 2926

Abstract

This study aimed to establish a finite element model that vividly reflected the anterior cruciate ligament (ACL) geometry and investigated the ACL stress distribution under different loading conditions. The ACL’s three-dimensional finite element model was based on a human cadaveric knee. Simulations of [...] Read more.

This study aimed to establish a finite element model that vividly reflected the anterior cruciate ligament (ACL) geometry and investigated the ACL stress distribution under different loading conditions. The ACL’s three-dimensional finite element model was based on a human cadaveric knee. Simulations of three loading conditions (134 N anterior tibial load, 5 Nm external tibial torque, 5 Nm internal tibial torque) on the knee model were performed. Experiments were performed on a knee specimen using a robotic universal force/moment sensor testing system to validate the model. The simulation results of the established model were in good agreement with the experimental results. Under the anterior tibial load, the highest maximal principal stresses (14.884 MPa) were localized at the femoral insertion of the ACL. Under the external and internal tibial torque, the highest maximal principal stresses (0.815 MPa and 0.933 MPa, respectively) were mainly concentrated in the mid-substance of the ACL and near the tibial insertion site, respectively. Combining the location of maximum stress and the location of common clinical ACL rupture, the most dangerous load during ACL injury may be the anterior tibial load. ACL injuries were more frequently loaded by external tibial than internal tibial torque. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

17 pages, 3716 KiB

Open AccessArticle

Layer-Specific Damage Modeling of Porcine Large Intestine under Biaxial Tension

by Aroj Bhattarai, Charlotte Anabell May, Manfred Staat, Wojciech Kowalczyk and Thanh Ngoc Tran

Bioengineering 2022, 9(10), 528; https://doi.org/10.3390/bioengineering9100528 - 6 Oct 2022

Cited by 4 | Viewed by 2183

Abstract

The mechanical behavior of the large intestine beyond the ultimate stress has never been investigated. Stretching beyond the ultimate stress may drastically impair the tissue microstructure, which consequently weakens its healthy state functions of absorption, temporary storage, and transportation for defecation. Due to [...] Read more.

The mechanical behavior of the large intestine beyond the ultimate stress has never been investigated. Stretching beyond the ultimate stress may drastically impair the tissue microstructure, which consequently weakens its healthy state functions of absorption, temporary storage, and transportation for defecation. Due to closely similar microstructure and function with humans, biaxial tensile experiments on the porcine large intestine have been performed in this study. In this paper, we report hyperelastic characterization of the large intestine based on experiments in 102 specimens. We also report the theoretical analysis of the experimental results, including an exponential damage evolution function. The fracture energies and the threshold stresses are set as damage material parameters for the longitudinal muscular, the circumferential muscular and the submucosal collagenous layers. A biaxial tensile simulation of a linear brick element has been performed to validate the applicability of the estimated material parameters. The model successfully simulates the biomechanical response of the large intestine under physiological and non-physiological loads. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

15 pages, 5409 KiB

Open AccessArticle

Finite Element Analysis of a New Non-Engaging Abutment System for Three-Unit Implant-Supported Fixed Dental Prostheses

by Soo-Hwan Byun, Joung-Hwa Seo, Ran-Yeong Cho, Sang-Min Yi, Lee-Kyong Kim, Hyun-Sook Han, Sung-Woon On, Won-Hyeon Kim, Hyun-Wook An and Byoung-Eun Yang

Bioengineering 2022, 9(10), 483; https://doi.org/10.3390/bioengineering9100483 - 20 Sep 2022

Cited by 4 | Viewed by 5505

Abstract

(1) Background: The stability of implants plays a significant role in the success of osseointegration. The stability of the connection between the fixture and the abutment is one of the critical factors affecting osseointegration. When restoring multiple, non-parallel, and splinted implants, achieving a [...] Read more.

(1) Background: The stability of implants plays a significant role in the success of osseointegration. The stability of the connection between the fixture and the abutment is one of the critical factors affecting osseointegration. When restoring multiple, non-parallel, and splinted implants, achieving a passive fit can be complicated and challenging. A new EZ post non-engaging abutment system of the BlueDiamond^® (BD) implant allows a wide connection angle while achieving a passive prosthesis fit. This study aimed to confirm the new abutment system’s clinical applicability by evaluating its biomechanical characteristics using finite element analysis (FEA). (2) Methods: The implant-supported fixed three-unit dental prostheses model was reproduced for two groups of AnyOne^® (AO) and BD implants using FEA. The loading conditions were a preload of 200 N in the first step and loads of 100 N (axial), 100 N (15°), or 30 N (45°) in the second step. (3) Results: The peak Von Mises stress (PVMS) value of the fixture in the BD group was more than twice that in the AO group. In contrast, the PVMS values of the abutment and abutment screws were lower in the BD group than in the AO group. The AO group revealed higher maximal principal stress (MPS) values than that of the BD group in the cortical bone, cancellous bone, and crown. The average stress of the outer surface of the abutment was lower in the AO group than in the BD group. The stress distribution for the inner surface of the fixture confirmed that the BD group displayed a lower stress distribution than the AO group under axial and 15° loads; however, the average stress was 1.5 times higher at the 45° load. The stress values of the entire surface where the cortical and cancellous bone were in contact with the fixture were measured. The AO group showed a higher stress value than the BD group in both cortical and cancellous bone. (4) Conclusions: In the AO group, the PVMS value of the fixture and the stress distribution at the contact surface between the fixture and the abutment were lower than those of the BD group, suggesting that the stability of the fixture would be high. However, due to the high stress in the fastening area of the abutment and abutment screw, the risk of abutment fracture in the AO group is higher than that of the BD group. Therefore, the new EZ post non-engaging abutment of the BD implant can be used without any problems in clinics, similar to the non-engaging abutment of the AO implant, which has been widely used in clinical practice. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Graphical abstract

10 pages, 2335 KiB

Open AccessArticle

Understanding Sex-Based Kinematic and Kinetic Differences of Chasse-Step in Elite Table Tennis Athletes

by Xiaoyi Yang, Qichang Mei, Shirui Shao, Wenjing Gu, Yuqi He, Ruizhe Zhu and Yaodong Gu

Bioengineering 2022, 9(6), 246; https://doi.org/10.3390/bioengineering9060246 - 4 Jun 2022

Cited by 4 | Viewed by 3126

Abstract

With the progress and innovation of table tennis technology, individualized training programs may deserve special attention. This study aimed to analyze elite table tennis athletes in chasse-step, with a particular focus on sex-based biomechanical differences. A total of 36 (18 males and 18 [...] Read more.

With the progress and innovation of table tennis technology, individualized training programs may deserve special attention. This study aimed to analyze elite table tennis athletes in chasse-step, with a particular focus on sex-based biomechanical differences. A total of 36 (18 males and 18 females) elite table tennis athletes performed topspin forehand of chasse-step. Angles and moments of hip, knee, and ankle joints were calculated using OpenSim (v4.2) with marker trajectories and ground reaction forces were measured via Vicon motion capture system and AMTI in-ground force platform. Males had greater hip and knee flexion angles during the entire motion phase and greater internal rotation angles of the hip during the forward swing phase. The joint stiffness of knee in males was greater than females in the frontal plane. Females in the forward swing phase showed greater hip flexion, adduction, and internal rotation moments than males. It was suggested that the difference may be due to the limitation of anatomical structures in sexes. Male table tennis athletes should strengthen lower extremity muscle groups to improve performance, while female table tennis athletes should focus on hip joint groups to avoid injury. The sex differences presented in this study could help coaches and athletes to develop individualized training programs for table tennis. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures

Figure 1

Review

Jump to: Research, Other

23 pages, 8026 KiB

Open AccessReview

Computational Modeling of Microwave Tumor Ablation

by Marija Radmilović-Radjenović, Nikola Bošković and Branislav Radjenović

Bioengineering 2022, 9(11), 656; https://doi.org/10.3390/bioengineering9110656 - 5 Nov 2022

Cited by 7 | Viewed by 2857

Abstract

Microwave ablation is recognized as a minimally invasive, fast-recovery treatment for destroying cancer cells using the heat generated by microwave energy. Despite the unquestionable benefits of microwave ablation, the interaction of the microwave applicator with the tissue may result in localized heating and [...] Read more.

Microwave ablation is recognized as a minimally invasive, fast-recovery treatment for destroying cancer cells using the heat generated by microwave energy. Despite the unquestionable benefits of microwave ablation, the interaction of the microwave applicator with the tissue may result in localized heating and damage to the surrounding tissue. The majority of the tissue damage can be removed by clarifying the conditions for their development. In addition to experimental methods, computer modeling has proven to be an effective tool for optimizing the performance of microwave ablation. Furthermore, because the thermal spread in biological tissue is difficult to measure, developing a predictive model from procedural planning to execution may have a substantial influence on patient care. The comprehension of heat transport in biological tissues plays a significant role in gaining insight into the mechanisms underlying microwave ablation. Numerical methods that enable ablation size control are required to guarantee tumor destruction and minimize damage to healthy tissues. Various values of input power and ablation time correspond to different tumor shapes ensuring the preservation of healthy tissues. The optimal conditions can be estimated by performing full three-dimensional simulations. This topical review recapitulates numerous computational studies on microwave tumor ablation. Novel areas emerging in treatment planning that exploit the advantages of numerical methods are also discussed. As an illustration, the results of the three-dimensional simulations of real liver tumors in the 3D-IRCADb-01 database are presented and analyzed. The simulation results confirm that numerical methods are very useful tools for modeling microwave tumor ablation with minimal invasiveness and collateral damage. Full article

(This article belongs to the Special Issue Computational Biomechanics)

► Show Figures