jinaojakezhang / FEDFEMBioheatDeform

Fast computation of soft tissue thermal response under deformation based on fast explicit dynamics finite element algorithm for surgical simulation

Fast computation of soft tissue thermal response under deformation based on fast explicit dynamics finite element algorithm for surgical simulation (MIT License)

This is the source repository for the paper:

Zhang, J., & Chauhan, S. (2020). Fast computation of soft tissue thermal response under deformation based on fast explicit dynamics finite element algorithm for surgical simulation. Computer Methods and Programs in Biomedicine, 187, 105244. doi:10.1016/j.cmpb.2019.105244.

which is based on the works of

[1] Zhang, J., & Chauhan, S. (2019). Fast explicit dynamics finite element algorithm for transient heat transfer. International Journal of Thermal Sciences, 139, 160-175. doi:10.1016/j.ijthermalsci.2019.01.030 [Code].
[2] Zhang, J., & Chauhan, S. (2019). Real-time computation of bio-heat transfer in the fast explicit dynamics finite element algorithm (FED-FEM) framework. Numerical Heat Transfer, Part B: Fundamentals, 75(4), 217-238. doi:10.1080/10407790.2019.1627812 [Code].

Please cite the above paper if you use this code for your research.

If this code is helpful in your projects, please help to ⭐ this repo or recommend it to your friends. Thanks😊

Environment:

• Windows 10
• Visual Studio 2017
• OpenMP

How to build:

1. Download the source repository.
2. Visual Studio 2017->Create New Project (Empty Project)->Project->Add Existing Item->BioheatDeform.cpp.
3. Project->Properties->C/C++->Language->OpenMP Support->Yes (/openmp).
4. Build Solution (Release/x64).

How to use:

1. (cmd)Command Prompt->build path>project_name.exe input.txt. Example:

   

2. Output: T.vtk, U.vtk, and Undeformed.vtk

How to visualize:

1. Open T.vtk and U.vtk. (such as using ParaView)

How to make input.txt:

1. Liver_Iso.inp (Abaqus input) is provided in the “models”, which was used to create Liver_Iso_n1.txt.

Material types:

1. Isotropic, orthotropic, and anisotropic thermal conductivities.
2. Neo-Hookean hyperelastic material

Boundary conditions (BCs):

1. Node index: Disp, FixP, HFlux, FixT.
2. Element index: Perfu, BodyHFlux.
3. All Elements: Gravity, Metabo.

Notes:

1. Node and Element index can start at 0, 1, or any but must be consistent in a file.
2. Index starts at 0: *.txt.
3. Index starts at 1: *_n1.txt.

Feedback:

Please send an email to [e-mail]. Thanks for your valuable feedback and suggestions.