Abstract:We develop efficient and accurate numerical methods to solve a class of shallow shell problems of the von Karman type. The governing equations form a fourth-order coupled system of nonlinear biharnomic equations for the transverse deflection and Airy's stress function. A second-order finite difference discretization with three iterative methods (Picard, Newton and Trust-Region Dogleg) are proposed for the numerical solution of the nonlinear PDE system. Three simple boundary conditions and two application-motivated mixed boundary conditions are considered. Along with the nonlinearity of the system, boundary singularities that appear when mixed boundary conditions are specified are the main numerical challenges. Two approaches that use either a transition function or local corrections are developed to deal with these boundary singularities. All the proposed numerical methods are validated using carefully designed numerical tests, where expected orders of accuracy and rates of convergence are observed. A rough run-time performance comparison is also conducted to illustrate the efficiency of our methods. As an application of the methods, a snap-through thermal buckling problem is considered. The critical thermal loads of shell buckling with various boundary conditions are numerically calculated, and snap-through bifurcation curves are also obtained using our numerical methods together with a pseudo-arclength continuation method. Our results are consistent with previous studies.
What problem does this paper attempt to address?
This paper attempts to solve the problem of numerical solution of shallow - shell equations under thermal stress. Specifically, the article focuses on developing efficient and accurate numerical methods to solve a class of shallow - shell problems, which are described by the nonlinear bi - harmonic equations of von Karman type. The main challenges include:
1. **Non - linear system**: The shallow - shell equations are non - linear, which makes it difficult to obtain analytical solutions.
2. **Boundary singularity**: When it comes to mixed boundary conditions, singularities will occur at the boundaries, which challenges the accuracy of standard numerical methods.
To solve the above problems, the author proposes the following methods:
- **Finite - difference discretization**: A second - order finite - difference scheme is used to discretize the non - linear partial differential equations.
- **Iterative solution methods**: Three iterative methods (Picard, Newton and Trust - Region Dogleg) are introduced to solve the discretized non - linear equations.
- **Handling boundary singularities**: To deal with the boundary singularities under partially clamped mixed boundary conditions, two methods are proposed:
- **Transition function method**: By introducing a smooth transition function, the boundary conditions are made to transition smoothly from one state to another, thus eliminating discontinuities.
- **Local asymptotic solution method**: The local asymptotic analysis solution is used to correct the numerical scheme near the singular points.
In addition, the author also validates all the proposed numerical methods and observes the expected accuracy and convergence rate through carefully designed numerical tests. Finally, as an application example, the thermally - induced buckling problem is studied, the critical thermal loads for shell buckling under different boundary conditions are calculated, and the corresponding bifurcation curves are obtained.
### Mathematical formulas
The main mathematical models involved in the article are as follows:
- **Governing equations**:
\[
\nabla^4 \varphi' = -\frac{1}{2Eh} L[w', w'] - \frac{1}{Eh} L[w_0', w'] - \nabla^2 N(T')
\]
\[
D \nabla^4 w' = L[\varphi', w'] + L[\varphi', w_0'] - \frac{1}{1 - \nu} \nabla^2 M(T') + P'
\]
- **Dimensionless governing equations**:
\[
\nabla^4 \varphi = -\frac{1}{2} L[w, w] - L[w_0, w] - f_\varphi
\]
\[
\nabla^4 w = L[w, \varphi] + L[w_0, \varphi] + f_w
\]
- **Definition of bilinear operator \(L\)**:
\[
L[u, v] \equiv u_{xx} v_{yy} + u_{yy} v_{xx} - 2 u_{xy} v_{xy}
\]
- **Thermal force and thermal moment**:
\[
N(T') = E \alpha \int_{-h/2}^{h/2} T'(z) \, dz
\]
\[
M(T') = E \alpha \int_{-h/2}^{h/2} T'(z) z \, dz
\]
These formulas and methods together form the core content of this paper, aiming to provide an effective numerical tool to study the behavior of shallow shells under thermal stress.