An efficient active-stress electromechanical isogeometric shell model for muscular thin film simulations

Michele Torre,Simone Morganti,Alessandro Nitti,Marco Donato de Tullio,Josef Kiendl,Francesco Silvio Pasqualini,Alessandro Reali

DOI: https://doi.org/10.1016/j.mechmat.2024.105046

2024-06-11

Abstract:We propose an isogeometric approach to model the deformation of active thin films using layered, nonlinear, Kirchhoff Love shells. Isogeometric Collocation and Galerkin formulations are employed to discretize the electrophysiological and mechanical sub-problems, respectively, with the possibility to adopt different element and time-step sizes. Numerical tests illustrate the capabilities of the active stress based approach to effectively simulate the contraction of thin films in both quasi-static and dynamic conditions.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: how to efficiently simulate the mechanical behavior of Muscular Thin Films (MTFs) under electrical stimulation. Specifically, the author proposes a method based on Isogeometric Analysis to model and simulate the deformation behavior of a thin - film structure composed of a cardiomyocyte layer and an elastic substrate during activation.

### Problem Background
1. **Challenges in Biomedical Applications**:
- It is very complicated to directly measure the force generated by tissues during contraction.
- Researchers usually indirectly evaluate force generation by measuring the deformation of artificially constructed objects such as muscle thin films.

2. **Limitations of Existing Methods**:
- Although the traditional finite element method can be used to simulate these constructs, it has a high computational cost when dealing with thin structures with a high aspect ratio.
- Existing electromechanical coupling models are less efficient in dealing with complex three - dimensional problems.

### Main Contributions of the Paper
1. **Proposing a New Isogeometric Shell Model**:
- The Kirchhoff - Love shell theory is used to describe the mechanical behavior of the thin film.
- The isogeometric Collocation and Galerkin methods are combined to discretize the electrophysiological and mechanical sub - problems respectively, which improves the computational efficiency.

2. **Introducing the Active Stress Formula**:
- A formula based on active stress is proposed to simulate the contraction caused by cell activation, instead of the commonly used active strain method.
- This enables the model to more accurately describe the mechanical response of the thin film during activation.

3. **Solving the Multi - Physics Coupling Problem**:
- The electrophysiological activity is coupled with mechanical deformation to simulate the dynamic behavior of muscle thin films under different conditions.

### Formula Summary
- **Monodomain Equation**:
\[
\begin{cases}
C_m \frac{\partial v}{\partial t}=\nabla\cdot(D\nabla v)-I_{ion}+I_{app} & \text{in } A_e\times[0, T] \\
n_e\cdot D\nabla v = I_n & \text{on }\partial A_e\times[0, T] \\
v(r_e, 0)=v_0 & \text{in }\bar{A}_e \text{ for }t = 0
\end{cases}
\]

- **Green - Lagrange Strain Tensor**:
\[
E_{\alpha\beta}=\frac{1}{2}(C_{\alpha\beta}-\bar{g}_{\alpha\beta})
\]

- **Active Stress Tensor**:
\[
S_a=\sigma_a\frac{\bar{f}_0\otimes\bar{f}_0}{\lambda_f^2}
\]

Through these improvements, the paper provides an efficient and accurate simulation method, which is suitable for studying the mechanical behavior of muscle thin films under electrical stimulation.