Abstract
The main aim of this paper is to propose a new boundary element algorithm for describing thermomechanical interactions in anisotropic soft tissues. The governing equations are studied based on the dual-phase lag bioheat transfer and Biot's theory. Due to the advantages of convolution quadrature boundary element method (CQBEM), such as low CPU usage, low memory usage and suitability for treatment of soft tissues that have complex shapes, it is a versatile and powerful method for modeling of bioheat distribution in anisotropic soft tissues and the related deformation. The resulting linear systems for bioheat and mechanical equations are solved by Transpose-free quasi-minimal residual (TFQMR) solver with a dual-threshold incomplete LU factorization technique (ILUT) preconditioner that reduces the iterations number and total CPU time. Numerical results demonstrate the validity, efficiency and accuracy of the proposed algorithm and technique.