Using a poroelastodynamic model to investigate the dynamic behaviour of articular cartilage

Dean Chou; Yun-Di Li; Zartasha Mustansar; Chen-Yuan Chung

doi:10.1016/j.cmpb.2023.107481

Using a poroelastodynamic model to investigate the dynamic behaviour of articular cartilage

Comput Methods Programs Biomed. 2023 May:233:107481. doi: 10.1016/j.cmpb.2023.107481. Epub 2023 Mar 12.

Authors

Dean Chou¹, Yun-Di Li², Zartasha Mustansar³, Chen-Yuan Chung⁴

Affiliations

¹ Department of Biomedical Engineering, National Cheng Kung University, Tainan City, Taiwan; Medical Device Innovation Center, National Cheng Kung University, Tainan City, Taiwan. Electronic address: dean@gs.ncku.edu.tw.
² Department of Mechanical Engineering, National Central University, Taoyuan City, Taiwan.
³ Department of National University of Science and Technology, Islamabad, Pakistan. Electronic address: zmustansar@sines.nust.edu.pk.
⁴ Department of Mechanical Engineering, National Central University, Taoyuan City, Taiwan. Electronic address: cychung@ncu.edu.tw.

PMID: 36963330
DOI: 10.1016/j.cmpb.2023.107481

Abstract

Background and objective: There is still a few studies about the poroelastic model that performed dynamic behaviour, especially for the case of the poroelastic cartilage model. Therefore, this study is aimed to use the poroelastodynamic model to simulate the dynamic behaviour of cartilage.

Methods: The governing equations of the poroelastodynamic model is firstly established. The validation of the model is initialised by modifying the equations into the static poroelastic model. The modified equations are then discretised using the finite element method. Mandel's problem is used to validate the discretised equations. The numerical solution calculated using FreeFEM++ is validated with the analytical solution for the quasi-static state and compared with the results generated using COMSOL Multiphysics software. Finally, the quasi-static solution is compared with the dynamic solution to discuss the difference in pore pressure and displacement variations of the poroelastic cartilage model.

Results: The dynamic solution showed transient behaviour at the beginning of the excitation. When the compressive force acts on the cartilage, there are obvious fluctuations during the initial stage and then the dynamic numerical solution gradually approaches the quasi-static value over a period of time. The deduced results of the analytical solution were approximately the same as the numerical simulation results.

Conclusion: This study was able to use the poroelastodynamics equation to simulate the dynamic behaviour of the poroelastic cartilage model. The comparison between the result coming from poroelastodynamics equation with that of the validated numerical solution was satisfactorily compared. The approximate similarity between the results of quasi-static and dynamic solutions underscored the importance of performing the dynamic solution for a more realistic simulation. This dynamic solution can be further used for the analysis of vibration or stress waves in future research.

Keywords: Cartilage; Mandel's problem; Poroelasticity; Poroelastodynamic.

MeSH terms

Cartilage, Articular*
Computer Simulation
Elasticity
Finite Element Analysis
Models, Biological
Pressure
Stress, Mechanical