A fast and efficient numerical method for computing the stress concentration between closely located stiff inclusions of general shapes

Xiaofei Li,Shengqi Lin,Haojie Wang

2024-07-09

Abstract:When two stiff inclusions are closely located, the gradient of the solution to the Lamé system, in other words the stress, may become arbitrarily large as the distance between two inclusions tends to zero. To compute the gradient of the solution in the narrow region, extremely fine meshes are required. It is a challenging problem to numerically compute the stress near the narrow region between two inclusions of general shapes as their distance goes to zero. A recent study [15] has shown that the major singularity of the gradient can be extracted in an explicit way for two general shaped inclusions. Thus the complexity of the computation can be greatly reduced by removing the singular term and it suffices to compute the residual term only using regular meshes. The goal of this paper is to numerically compute the stress concentration in a fast and efficient way. In this paper, we compute the value of the stress concentration factor, which is the normalized magnitude of the stress concentration, for general shaped domain as the distance between two inclusions tends to zero. We also compute the solution for two closely located inclusions of general shapes and show the convergence of the solution. Only regular meshes are used in our numerical computation and the results clearly show that the characterization of the singular term method can be efficiently used for computation of the stress concentration between two closely located inclusions of general shapes.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to calculate the stress concentration in the narrow area between two irregular - shaped rigid inclusions quickly and efficiently when they are very close to each other.

Specifically, when the distance between two rigid inclusions approaches zero, the gradient of the solution of the Lamé system (i.e., stress) may become arbitrarily large. Therefore, calculating the stress gradient in the narrow area requires extremely fine meshing, which is a huge challenge in numerical calculations. This paper proposes a new numerical method. By removing the singular terms to simplify the computational complexity and using a regular grid to calculate the remaining terms, it can achieve fast and efficient calculation of stress concentration between inclusions of general shapes.

### Main contributions of the paper:
1. **Extracting singular terms**: Using existing research results, the main singularity of the stress gradient can be explicitly extracted, thereby greatly reducing the computational complexity.
2. **Efficient calculation**: By removing the singular terms, only the remaining terms need to be calculated with a regular grid, which significantly improves the computational efficiency.
3. **Convergence verification**: It proves the effectiveness of the proposed method in numerical calculations and shows the convergence of the calculation results as the number of grids increases.

### Mathematical description:
Let \(D_1\) and \(D_2\) be two strictly convex, simply - connected domains with smooth boundaries that are very close to each other. Define the distance \(\epsilon=\text{dist}(D_1, D_2)\) and assume that \(\epsilon\) is very small. Let \(z_1\in\partial D_1\) and \(z_2\in\partial D_2\) be the closest points, satisfying:

\[|z_1 - z_2|=\text{dist}(D_1, D_2)\]

The goal of the paper is to calculate the stress concentration factor, that is, the normalized amplitude of stress concentration, for a fast and efficient numerical calculation for domains of general shapes when the distance between the two inclusions approaches zero.

### Key formulas:
- The expression for the stress concentration factor \(\alpha_0\) is:
\[
\alpha_0=\frac{4\pi r_1 r_2}{r_1 + r_2}(\nu\cdot\nabla H)(z)
\]
where \(r_1\) and \(r_2\) are the radii of the two discs respectively, and \(H\) is the applied field.

- The definition of the singular function \(q(x)\):
\[
q(x)=\frac{1}{2\pi}(\ln|x - p_1|-\ln|x - p_2|)
\]
where \(p_1\) and \(p_2\) are the fixed points of the mixed reflection.

Through these methods, the paper achieves fast and efficient calculation of stress concentration between inclusions of general shapes and verifies the effectiveness and convergence of the method.