Abstract
The principal objective of this work was to develop a semi-implicit hybrid boundary element method (HBEM) to describe the nonlinear fractional biomechanical interactions in functionally graded anisotropic (FGA) soft tissues. The local radial basis function collocation method (LRBFCM) and general boundary element method (GBEM) were used to solve the nonlinear fractional dual-phase-lag bioheat governing equation. The boundary element method (BEM) was then used to solve the poroelastic governing equation. To solve equations arising from boundary element discretization, an efficient partitioned semi-implicit coupling algorithm was implemented with the generalized modified shift-splitting (GMSS) preconditioners. The computational findings are presented graphically to display the influences of the graded parameter, fractional parameter, and anisotropic property on the bio-thermal stress. Different bioheat transfer models are presented to show the significant differences between the nonlinear bio-thermal stress distributions in functionally graded anisotropic biological tissues. Numerical findings verified the validity, accuracy, and efficiency of the proposed method.