The polygonal scaled boundary thin plate element based on the discrete Kirchhoff theory

Chong-Jun Li,Ying Zhang,Yan-Mei Jia,Juan Chen

DOI: https://doi.org/10.1016/j.camwa.2021.05.036

2021-09-01

Abstract:<p>The scaled boundary finite element method (SBFEM) is a powerful method for solving elastostatics problems based on polygonal elements. In this paper, we firstly construct the quadratic polygonal scaled boundary element only depends on the boundary nodal displacements by transforming the additional degrees of freedom derived from the constant body loads to those by the boundary nodes. Further, combining with the discrete Kirchhoff theory, we construct the polygonal scaled boundary thin plate element, which can possess the second order completeness. The element stiffness matrix for the thin plate problem can be transformed by the stiffness matrix for the plane problem directly by avoiding to compute the shape functions of SBFEM. Numerical examples verify that the proposed polygonal scaled boundary thin plate element has good accuracy.</p>

mathematics, applied

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: in the polygonal thin - plate element, combine the discrete Kirchhoff theory and the scaled boundary finite element method (SBFEM) to construct a polygonal scaled - boundary thin - plate element with second - order completeness, and provide a method for directly calculating the element stiffness matrix without calculating the shape function.

Specifically, the paper mainly focuses on the following points:

1. **Constructing polygonal elements**: Based on the discrete Kirchhoff theory, construct polygonal scaled - boundary thin - plate elements. Different from traditional triangular or quadrilateral elements, polygonal elements provide greater flexibility in mesh generation and solution accuracy.

2. **Handling constant body forces**: By converting the extra degrees of freedom caused by constant body forces into degrees of freedom on boundary nodes, construct quadratic polygonal scaled - boundary elements. This enables the element to maintain second - order completeness.

3. **Avoiding shape function calculation**: Propose a method to directly derive the stiffness matrix of the thin - plate problem from the stiffness matrix of the plane problem without calculating the shape function. This is different from traditional methods based on given shape functions, such as DKT or DKQ elements.

4. **Numerical verification**: Verify through numerical examples that the proposed polygonal scaled - boundary thin - plate element has good accuracy.

The core contribution of the paper lies in combining SBFEM and the discrete Kirchhoff theory, solving the construction problem of polygonal thin - plate elements, and providing an efficient method for calculating the stiffness matrix, thereby enhancing the application potential of polygonal meshes in thin - plate bending problems.

### Formula summary

- **Scaled - boundary coordinate transformation**:
\[
\begin{cases}
x = x_0+\xi(x_b(\eta)-x_0)=x_0+\xi(N(\eta)x_b - x_0),\\
y = y_0+\xi(y_b(\eta)-y_0)=y_0+\xi(N(\eta)y_b - y_0),
\end{cases}
\]
where \(x_b = [x_i, x_{n + i}, x_{i+1}]^T\), \(y_b = [y_i, y_{n + i}, y_{i+1}]^T\), \(N(\eta)=[N_1(\eta), N_2(\eta), N_3(\eta)]^T\) is a one - dimensional quadratic Lagrange function vector.

- **Stiffness matrix calculation**:
\[
K = Y^T(W\circ G)Y,
\]
where \(W\) and \(G\) are respectively:
\[
W=\begin{bmatrix}
0&0&0&\cdots&0&0\\
0&0&0&\cdots&0&0\\
0&0&\frac{1}{2\lambda_3}&\cdots&\frac{1}{\lambda_3+\lambda_{4n}}&\frac{1}{\lambda_3 + 2}\\
\vdots&\vdots&\vdots&\ddots&\vdots&\vdots\\
0&0&\frac{1}{\lambda_3+\lambda_{4n}}&\cdots&\frac{1}{2\lambda_{4n}}&\frac{1}{\lambda_{4n}+2}\\
0&0&\frac{1}{\lambda_3+2}&\cdots&\frac{1}{\lambda_{4n}+2}&\frac{1}{4}\\
0&0&\frac{1}{\lambda_3+2}&\cdots&\frac{1}{\lambda_{4n}+2}&\frac{1}{4}
\end{bmatrix}