The $ AA $-iterative algorithm in hyperbolic spaces with applications to integral equations on time scales

Aynur Şahin; Zeynep Kalkan; Aynur Şahin; Zeynep Kalkan

doi:10.3934/math.20241192

AIMS Mathematics

2024, Volume 9, Issue 9: 24480-24506. doi: 10.3934/math.20241192

Previous Article Next Article

Research article Special Issues

The $AA$ -iterative algorithm in hyperbolic spaces with applications to integral equations on time scales

Aynur Şahin ^,,
Zeynep Kalkan

Department of Mathematics, Faculty of Sciences, Sakarya University, Sakarya 54050, Türkiye

Received: 18 March 2024 Revised: 08 August 2024 Accepted: 16 August 2024 Published: 21 August 2024
MSC : 45B05, 47H09, 47H10

We explored the $AA$ -iterative algorithm within the hyperbolic spaces (HSs), aiming to unveil a stability outcome for contraction maps and convergence outcomes for generalized $(\alpha, \beta)$ -nonexpansive ( $G\alpha \beta N$ ) maps in such spaces. Through this algorithm, we derived compelling outcomes for both strong and $\Delta$ -convergence and weak $w^2$ -stability. Furthermore, we provided an illustrative example of $G\alpha \beta N$ maps and conducted a comparative analysis of convergence rates against alternative iterative methods. Additionally, we demonstrated the practical relevance of our findings by applying them to solve the linear Fredholm integral equations (FIEs) and nonlinear Fredholm-Hammerstein integral equations (FHIEs) on time scales.

Keywords:

hyperbolic space,
fixed point,
generalized $(\alpha, \beta)$ -nonexpansive map,
the AA-iterative algorithm,
time scale,
integral equations,
weak $w^2$ -stability

Citation: Aynur Şahin, Zeynep Kalkan. The $AA$ -iterative algorithm in hyperbolic spaces with applications to integral equations on time scales[J]. AIMS Mathematics, 2024, 9(9): 24480-24506. doi: 10.3934/math.20241192

Related Papers:

[1]	Nabil Kerdid . Asymptotic analysis of stretching modes for a folded plate. AIMS Mathematics, 2023, 8(10): 23974-23988. doi: 10.3934/math.20231222
[2]	Javier de Lucas, Julia Lange, Xavier Rivas . A symplectic approach to Schrödinger equations in the infinite-dimensional unbounded setting. AIMS Mathematics, 2024, 9(10): 27998-28043. doi: 10.3934/math.20241359
[3]	Muhammad Afzal, Tariq Ismaeel, Azhar Iqbal Kashif Butt, Zahid Farooq, Riaz Ahmad, Ilyas Khan . On the reducibility of a class of almost-periodic linear Hamiltonian systems and its application in Schrödinger equation. AIMS Mathematics, 2023, 8(3): 7471-7489. doi: 10.3934/math.2023375
[4]	Maxime Krier, Julia Orlik . Solvability of a fluid-structure interaction problem with semigroup theory. AIMS Mathematics, 2023, 8(12): 29490-29516. doi: 10.3934/math.20231510
[5]	Maha M. Helmi, Ali M. Mubaraki, Rahmatullah Ibrahim Nuruddeen . Analysis of horizontally polarized shear waves on a highly inhomogeneous loaded bi-material plate. AIMS Mathematics, 2023, 8(1): 2117-2136. doi: 10.3934/math.2023109
[6]	Nabil Kerdid . Asymptotic analysis of high frequency modes for thin elastic plates. AIMS Mathematics, 2023, 8(8): 18618-18630. doi: 10.3934/math.2023948
[7]	Min Zhang, Yi Wang, Yan Li . Reducibility and quasi-periodic solutions for a two dimensional beam equation with quasi-periodic in time potential. AIMS Mathematics, 2021, 6(1): 643-674. doi: 10.3934/math.2021039
[8]	Jia Li, Xia Li, Chunpeng Zhu . Reducibility for a class of almost periodic Hamiltonian systems which are degenerate. AIMS Mathematics, 2023, 8(1): 2296-2307. doi: 10.3934/math.2023119
[9]	Syed T. R. Rizvi, Sana Ghafoor, Aly R. Seadawy, Ahmed H. Arnous, Hakim AL Garalleh, Nehad Ali Shah . Exploration of solitons and analytical solutions by sub-ODE and variational integrators to Klein-Gordon model. AIMS Mathematics, 2024, 9(8): 21144-21176. doi: 10.3934/math.20241027
[10]	Aziz Ur Rehman, Muhammad Bilal Riaz, Ilyas Khan, Abdullah Mohamed . Time fractional analysis of Casson fluid with application of novel hybrid fractional derivative operator. AIMS Mathematics, 2023, 8(4): 8185-8209. doi: 10.3934/math.2023414

Abstract

1. Introduction

The moderately thick plate is an important material that can adapt to high pressure, high temperature and strong radiation environments, often used in construction engineering, machinery manufacturing, container manufacturing and so on ^[1,2]. In recent years, many achievements have been made in numerical methods of moderately thick plate problems. There are some typical numerical methods, such as the meshless local radial point interpolation method ^[3], discrete singular convolution methods ^[4], spectral element method ^[5], meshless local Petrov-Galerkin method ^[6], differential quadrature element method ^[7], finite integral transform method ^[8] and chaotic radial basis function method ^[9]. Compared with numerical methods, the analytical method seems more accurate, but it is difficult to construct. Through the efforts of many scholars, the analytical method has made great progress. The common analytical methods are integral transformation ^[10,11], semi-inverse method ^[12], infinite series method ^[13] and Green's functions ^[14]. Traditionally, some analytical solutions always need to eliminate the unknown functions by various methods in the range of a class of variables so as to obtain a partial differential equation with higher order, which will make separating variables and other mathematical methods impossible. In order to avoid this phenomenon, academician Zhong put forward the SEA in the 1990s. The essence of SEA is the method of separation of variables based on Hamiltonian system, which does not need to assume the form of the solution in advance. The main feature of the approach is that the solution procedure is conducted in the symplectic space with dual variables rather than in the Euclidean space with one kind of variables, which is rational to solve the equations of elasticity ^[15,16]. Within the framework of the Hamiltonian system, the analytical solutions of the considered problem can be obtained by using the complete expansion of symplectic eigenvectors without any prior assumptions of the solution forms, which shows the distinctive advantage of the SEA.

SEA needs to transform the equations of elasticity into a proper Hamiltonian system and can uniformly solve the state vector, including as many unknown functions, such as the physical displacements and forces, as possible. Nowadays, the SEA has been extensively studied by some scholars. Xu et al. ^[17] presented the fracture analysis of fractional two-dimensional viscoelastic media by SEA. Zhao et al. ^[18] utilized SEA on the wave propagation problem for periodic structures with uncertainty. Qiao et al. ^[19] considered the plane elasticity problems of two-dimensional octagonal quasicrystals by SEA. Zhou et al. ^[20] studied new buckling solutions of moderately thick rectangular plates by the symplectic superposition method within the Hamiltonian system framework. Li et al. ^[21] explored new analytic free vibration solutions of doubly curved shallow shells. Su et al. ^[22] obtained an analytical free vibration solution of fully free orthotropic rectangular thin plates on two-parameter elastic foundations by the symplectic superposition method. In the last thirty years, the application of SEA has been extended continuously in the electromagnetic elastic solid problem ^[23], non-Lévy-type functionally graded rectangular plates ^[24], special mechanical problems in mine engineering ^[25], buckling of regular and auxetic honeycombs ^[26], free vibration of non-Lévy-type porous FGM rectangular plates ^[27], the free vibration of edge-cracked thick rectangular plates ^[28], the buckling problem of CNT reinforced composite rectangular plates ^[29] and so on.

In this paper, based on the idea of ^[30], the mechanical problems of rectangular moderately thick plates are analyzed and a unified framework for a class of rectangular moderately thick plates under the simply supported boundary conditions at opposite sides is given. The model considered in this paper contains free parameters, including the bending, buckling and free vibration problems of the plates. We start from the basic equations of the rectangular moderately thick plates, guide the equations to the Hamiltonian system by introducing appropriate state vector and solve the eigenvalue problem of the Hamiltonian operators. More importantly, we prove the completeness of the eigenfunction system, which lays a theoretical foundation for our general solution. For the applications of the model, based on the obtained general solutions, rectangular moderately thick plates under Winkler type foundation and Pasternak type two-parameter foundation are numerically simulated, which have not been solved by SEA so far. The free vibration problem of moderately thick plates with fully simply supported is also discussed. The satisfactory agreement further confirms the stability of the present method.

2. The model of the moderately thick plates and SEA

2.1. Basic equations and Hamiltonian system

For the problems of rectangular moderately thick plates, we propose the following general model:

$\begin{align} &\frac{\partial M_{x}}{\partial x}+\frac{\partial M_{xy}}{\partial y}-Q_{x} = q_{1}(x,y), \end{align}$

(2.1a)

$\begin{align} &\frac{\partial M_{y}}{\partial y}+\frac{\partial M_{xy}}{\partial x}-Q_{y} = q_{2}(x,y), \end{align}$

(2.1b)

$\begin{align} &\frac{\partial Q_{x}}{\partial x}+\frac{\partial Q_{y}}{\partial y}+\alpha_{1} \frac{\partial^{2} w}{\partial x^{2}}+\beta_{1} \frac{\partial^{2} w}{\partial y^{2}}+\gamma w = q_{3}(x,y), \end{align}$

(2.1c)

where $q_{1}$ , $q_{2}$ , $q_{3}$ are loads, $\alpha_{1}$ , $\beta_{1}$ , $\gamma$ are free parameters, $M_{x}$ , $M_{y}$ , $M_{xy}$ , $Q_{x}$ , $Q_{y}$ , $w$ are the bending moments, torsional moment, shear forces and the transverse modal displacement, respectively.

The internal forces of the plates can be presented as

$\begin{align} &M_{x} = -D(\frac{\partial \psi_{x}}{\partial x}+\nu\frac{\partial \psi_{y}}{\partial y}), \end{align}$

(2.2a)

$\begin{align} &M_{y} = -D(\frac{\partial \psi_{y}}{\partial y}+\nu\frac{\partial \psi_{x}}{\partial x}), \end{align}$

(2.2b)

$\begin{align} &M_{xy} = -\frac{D(1-\nu)}{2}(\frac{\partial \psi_{x}}{\partial y}+\frac{\partial \psi_{y}}{\partial x}), \end{align}$

(2.2c)

$\begin{align} &Q_{x} = C(\frac{\partial w}{\partial x}-\psi_{x}), \end{align}$

(2.2d)

$\begin{align} &Q_{y} = C(\frac{\partial w}{\partial y}-\psi_{y}), \end{align}$

(2.2e)

where $C = \frac{5Eh}{12(1+\nu)}, \ D = \frac{Eh^3}{12(1-\nu^2)}$ are the elastic constants, in which $E$ denotes the modulus of elasticity, $h$ represents the thickness of the plate and $\nu$ is Poisson's ratio. From (2.1a)–(2.2e), the model can be transformed into the following form:

$\begin{align} &( -D \frac{\partial^2}{\partial x^2}-\frac{D(1-\nu)}{2} \frac{\partial^2}{\partial y^2}+C)\psi_{x}-\frac{D(1+\nu)}{2}\frac{\partial^2}{\partial x \partial y}\psi_{y}-C\frac{\partial}{\partial x}w = q_{1}(x,y), \end{align}$

(2.3a)

$\begin{align} &( -D \frac{\partial^2}{\partial y^2}-\frac{D(1-\nu)}{2} \frac{\partial^2}{\partial x^2}+C)\psi_{y}-\frac{D(1+\nu)}{2}\frac{\partial^2}{\partial x \partial y}\psi_{x}-C\frac{\partial}{\partial y}w = q_{2}(x,y), \end{align}$

(2.3b)

$\begin{align} &-C \frac{\partial}{\partial x}\psi_{x}-C\frac{\partial}{\partial y}\psi_{y}+(\alpha\frac{\partial^2}{\partial x^2}+\beta\frac{\partial^2}{\partial y^2}+\gamma)w = q_{3}(x,y), \end{align}$

(2.3c)

where $\alpha = \alpha_{1}+C$ , $\beta = \beta_{1}+C$ .

Setting $V_x = \alpha\frac{\partial w}{\partial x} -C \psi _x$ and combing (2.2a)–(2.3c), we can obtain the following matrix equation:

$\begin{align} \frac{\partial {\textbf{U}}(x,y)}{\partial x} = {\textbf{H}} U(x,y)+{\textbf{f}}(x,y). \end{align}$

(2.4)

It is easy to verify that ${\textbf{H}} = \left(\begin{array}{cc} {\textbf{0}} & {\textbf{T}} \\ {\textbf{S}} & {\textbf{0}} \\ \end{array} \right)$ satisfies ${\textbf{H}}^* = {\textbf{JHJ}}$ in which ${\textbf{J}} = \left(\begin{array}{cc} {\textbf{0}} & {\textbf{I}} _3\\ -{\textbf{I}}_3 & {\textbf{0}} \\ \end{array} \right)$ , ${\textbf{I}}_3$ is a $3\times3$ unit matrix. So ${\textbf{H}}$ is a Hamiltonian operator matrix, in which

${\textbf{ T}} = \left( \begin{array}{ccc} -\frac{1}{D} & -\nu \frac{\partial}{\partial y} & 0 \\ \nu \frac{\partial}{\partial y} & D\left( \nu^2-1\right) \frac{\partial^2}{\partial y^2}+C & -C \frac{\partial}{\partial y} \\ 0 & C\frac{\partial}{\partial y} & -\beta\frac{\partial^2}{\partial y^2} -\gamma \\ \end{array} \right),\,\,\,\, {\textbf{S}} = \left(\begin{array}{ccc} \frac{C^2}{\alpha }-C & \frac{\partial}{\partial y} & \frac{C}{\alpha } \\ -\frac{\partial}{\partial y} & -\frac{2}{D (\nu-1)} & 0 \\ \frac{C}{\alpha } & 0 & \frac{1}{\alpha } \\ \end{array} \right),$

${\textbf{f}}(x, y) = (0, -q_{2}(x, y), q_{3}(x, y), q_{1}(x, y), 0, 0)^T$ , ${\textbf{U}}(x, y) = (\psi _x, -M_{xy}, V_x, M_x, \psi _y, w)^T$ .

Next, we continue to use the SEA to solve the Hamiltonian system (2.4).

2.2. Eigenvalues and eigenvectors of ${\textbf{H}}$

To solve the Hamiltonian system, the homogeneous system is first considered:

$\begin{align} \frac{\partial {\textbf{U}}(x,y)}{\partial x} = {\textbf{H}} {\textbf{U}}(x,y). \end{align}$

(2.5)

Using the separation of variables method, assuming $U(x, y) = X(x)Y(y)$ , we have

$\begin{align} \frac{{\mathrm{d}}X(x)}{{\mathrm{d}}x} = \mu X(x), \end{align}$

(2.6)

$\begin{align} {\textbf{H}}Y{\textbf{}}(y) = \mu {\textbf{Y}}(y), \end{align}$

(2.7)

with eigenvalue $\mu$ and eigenvector ${\textbf{Y}}(y)$ .

The considered model is solved in the rectangular region $\{(x, y)\mid 0 \leq x\leq a, 0 \leq y \leq b\}$ and the boundary conditions are simply supported at $y = 0$ and $y = b$ as

$\begin{align} w = \psi _x = M_y = 0. \end{align}$

(2.8)

The eigenvalue problems (2.6) and (2.7) are actually quite complicated to compute. For the sake of brevity below, let $\zeta_n = \frac{n \pi}{b}$ and

$\begin{equation*} P_{n} = \sqrt{\frac{C^4-2\alpha C^3+C^2 \left(\alpha ^2+2 D \zeta_n^2 (\alpha-\beta )+2\gamma D\right)+2\alpha C D \left(\zeta_n^2 (\alpha-\beta )+\gamma \right)+D^2 \left(\zeta_n^2 (\alpha-\beta )+\gamma \right)^2}{D^2}}, \end{equation*}$

$n\in \mathbb{N}$ , where $\mathbb{N}$ stands for the set of positive integers.

With the help of symbolic software, the corresponding results may be divided into the following three cases.

Case 1. $\mu_{1} = \frac{\sqrt{2} \sqrt{-C}}{\sqrt{D (\nu-1)}}, \mu_{2} = -\mu_{1}$ are eigenvalues of ${\textbf{H}}$ and the corresponding eigenvectors are

${\textbf{Y}}_{1} = (0, C, 0, 0, \mu_{1}, 0)^{T}, {\textbf{Y}}_{2} = (0, C, 0, 0, \mu_{2}, 0)^{T},$

where the superscript $T$ signifies the transpose of a vector.

Case 2. When $P_{n}\neq0$ , ${\textbf{H}}$ has six simple eigenvalues, which can be expressed as

$\mu_{n1} = \frac{\sqrt{D (\nu-1) \zeta _n^2-2 C}}{\sqrt{D (\nu-1)}},\, \mu_{n2} = -\mu _{n1};$

$\mu_{n3} = \sqrt{-\frac{C^2+D \left(\zeta _n^2 (-(\alpha +\beta ))+P_n+\gamma \right)-\alpha C}{2 \alpha D}},\, \mu_{n4} = -\mu_{n3}$

and

$\mu_{n5} = \sqrt{\frac{-C^2+D \left(\zeta _n^2 (\alpha +\beta )+P_n-\gamma \right)+\alpha C}{2 \alpha D}},\, \mu_{n6} = -\mu_{n5}\ (n\in \mathbb{N}).$

The corresponding eigenvectors with respect to $\mu_{ni}$ can be rewritten as

$\begin{align*} {\textbf{Y}}_{n, i} = \left( \begin{array}{c} \sin \left(y \zeta _n\right) \\ -\frac{D (\nu -1) \cos \left(y \zeta _n\right) \left(C^2 \left(\zeta _n^2+\mu_{ni}^2\right)-C \left(\gamma -\beta \zeta _n^2+\alpha \mu_{ni}^2\right)+D \left(\nu \zeta _n^2+\mu_{ni}^2\right) \left(\gamma -\beta \zeta _n^2+\alpha \mu_{ni}^2\right)\right)}{\zeta _n \left(2 C^2+D (\nu +1) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)} \\ \frac{C \sin \left(y \zeta _n\right) \left(-2 C^2+2 \alpha C-\gamma D (\nu +1)+D \zeta _n^2 (\alpha (-\nu )+\alpha +\beta \nu +\beta )-2 \alpha D \mu _{ni}^2\right)}{2 C^2+D (\nu +1) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)} \\ \frac{D \sin \left(y \zeta _n\right) \left(2 C^2 (\nu -1) \mu _{ni}^2-2 C \nu \left(\gamma -\beta \zeta _n^2+\alpha \mu_{ni}^2\right)+D (\nu -1) \left(\nu \zeta _n^2+\mu _{ni}^2\right) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)}{\mu_{ni} \left(2 C^2+D (\nu +1) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)} \\ \frac{\cos \left(y \zeta _n\right) \left(2 C^2 \mu_{ni}^2-2 C \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)+D \left((\nu -1) \zeta _n^2+2 \mu_{ni}^2\right) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)}{\zeta_n \mu_{ni} \left(2 C^2+D (\nu +1) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)} \\ \frac{C \sin \left(y \zeta _n\right) \left(2 C+D (\nu -1) \left(\mu _{ni}^2-\zeta _n^2\right)\right)}{\mu _{ni} \left(2 C^2+D (\nu +1) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)} \end{array} \right), \end{align*}$

( $i = 1, 2, \cdots, 6$ ; $n\in \mathbb{N}$ ).

Case 3. When $P_n = 0$ , ${\textbf{H}}$ contains two single eigenvalues:

$\mu_{n1} = \frac{\sqrt{D (\nu-1) \zeta _n^2-2 C}}{\sqrt{D (\nu-1)}}, \mu_{n2} = -\mu_{n1},$

and two double eigenvalues:

$\mu_{n3} = \sqrt{\frac{C (\alpha -C)+D \zeta _n ^2 (\alpha +\beta )-\gamma D}{2 \alpha D}}, \mu_{n4} = -\mu_{n3} \, (n\in \mathbb{N}).$

In this case, the basic eigenvectors corresponding to $\mu_{n1} \, (i = 1, 2, 3, 4; n\in \mathbb{N})$ are

$\begin{align*} {\textbf{Y}}_{n,i}^0 = \left( \begin{array}{c} \sin \left(y \zeta _n\right) \\ -\frac{D (\nu -1) \cos \left(y \zeta _n\right) \left(C^2 \left(\zeta _n^2+\mu_{ni}^2\right)-C \left(\gamma -\beta \zeta _n^2+\alpha \mu_{ni}^2\right)+D \left(\nu \zeta _n^2+\mu _{ni}^2\right) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)}{\zeta _n \left(2 C^2+D (\nu +1) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)} \\ \frac{C \sin \left(y \zeta _n\right) \left(-2 C^2+2 \alpha C-\gamma D (\nu +1)+D \zeta _n^2 (\alpha (-\nu )+\alpha +\beta \nu +\beta )-2 \alpha D \mu _{ni}^2\right)}{2 C^2+D (\nu +1) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)} \\ \frac{D \sin \left(y \zeta _n\right) \left(2 C^2 (\nu -1) \mu _{ni}^2-2 C \nu \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)+D (\nu -1) \left(\nu \zeta _n^2+\mu _{ni}^2\right) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)}{\mu _{ni} \left(2 C^2+D (\nu +1) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)} \\ \frac{\cos \left(y \zeta _n\right) \left(2 C^2 \mu_{ni}^2-2 C \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)+D\left((\nu -1) \zeta _n^2+2 \mu _{ni}^2\right) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)}{\zeta _n \mu _{ni} \left(2 C^2+D (\nu +1) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)} \\ \frac{C \sin \left(y \zeta _n\right) \left(2 C+D (\nu -1) \left(\mu _{ni}^2-\zeta _n^2\right)\right)}{\mu _{ni} \left(2 C^2+D (\nu +1) \left(\gamma -\beta \zeta _n^2+\alpha \mu _{ni}^2\right)\right)} \\ \end{array} \right). \end{align*}$

As for the situation of double eigenvalues, there exist first-Jordan form eigenvectors. Solving $\gamma$ from $P_n = 0$ yields that

$\gamma = \gamma_{n1} = \frac{-2 \sqrt{\alpha } C^{\frac{3}{2}}-C^2-\alpha C-\alpha D \zeta _n^2+\beta D \zeta _n^2}{D}$

and

$\gamma = \gamma_{n2} = \frac{2 \sqrt{\alpha } C^{\frac{3}{2}}-C^2-\alpha C-\alpha D \zeta _n^2+\beta D \zeta _n^2}{D}$

for further simplification. According to ${\textbf{H}}Y_{n, i}^1(y) = \mu_n {\textbf{Y}}_{n, i}^1(y)+{\textbf{Y}}_{n, i}^0(y), \ (i = 3, 4, n\in \mathbb{N}).$ The first-order Jordan eigenvectors are simplified as follows.

Case 3.1. If $\gamma = \gamma_{n1} = \frac{-2 \sqrt{\alpha } C^{\frac{3}{2}}-C^2-\alpha C-\alpha D \zeta _n^2+\beta D \zeta _n^2}{D},$ then

$\begin{align*} {\textbf{Y}}_{n,3}^1 = \left( \begin{array}{c} \frac{\sin \left(y \zeta _n\right) \left(-C \left(2 \sqrt{\alpha }+\sqrt{C}\right)-2 \sqrt{\alpha } D \zeta _n^2\right)}{C^{\frac{3}{2}} \sqrt{\frac{\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\\ \frac{D (\nu -1) \zeta _n \cos \left(y \zeta _n\right) \left(C^{\frac{3}{2}}+2 \sqrt{\alpha } C+2 \sqrt{\alpha } D \zeta _n^2\right)}{C^{\frac{3}{2}} \sqrt{\frac{\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\\ \frac{\sin \left(y \zeta _n\right) \left(C^{\frac{3}{2}}+\sqrt{\alpha } C+2 \sqrt{\alpha } D \zeta _n^2\right)}{\sqrt{C} \sqrt{\frac{\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\\ \frac{2 \sqrt{\alpha } D \sin \left(y \zeta _n\right) \left(C+(D-D \nu ) \zeta _n^2\right)}{C^{\frac{3}{2}}}\\ -\frac{2 \sqrt{\alpha } D \zeta _n \cos \left(y \zeta _n\right)}{C^{\frac{3}{2}}}\\ 0\\ \end{array} \right), {\textbf{Y}}_{n,4}^1 = \left( \begin{array}{c} \frac{\sin \left(y \zeta _n\right) \left(C^{\frac{3}{2}}+2 \sqrt{\alpha } C+2 \sqrt{\alpha } D \zeta _n^2\right)}{C^{\frac{3}{2}} \sqrt{\frac{\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\\ -\frac{D (\nu -1) \zeta _n \cos \left(y \zeta _n\right) \left(C^{\frac{3}{2}}+2 \sqrt{\alpha } C+2 \sqrt{\alpha } D \zeta _n^2\right)}{C^{\frac{3}{2}} \sqrt{\frac{\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\\ \frac{\sin \left(y \zeta _n\right) \left(-C \left(\sqrt{\alpha }+\sqrt{C}\right)-2 \sqrt{\alpha } D \zeta _n^2\right)}{\sqrt{C} \sqrt{\frac{\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\\ \frac{2 \sqrt{\alpha } D \sin \left(y \zeta _n\right) \left(C+(D-D \nu ) \zeta _n^2\right)}{C^{\frac{3}{2}}}\\ -\frac{2 \sqrt{\alpha }D \zeta _n \cos \left(y \zeta _n\right)}{C^{\frac{3}{2}}}\\ 0\\ \end{array} \right). \end{align*}$

Case 3.2. If $\gamma = \gamma_{n2} = \frac{2 \sqrt{\alpha } C^{\frac{3}{2}}-C^2-\alpha C-\alpha D \zeta _n^2+\beta D \zeta _n^2}{D},$ we have

$\begin{align*} {\textbf{Y}}_{n,3}^1 = \left( \begin{array}{c} \frac{\sin \left(y \zeta _n\right) \left(-C \left(-2 \sqrt{\alpha }+\sqrt{C}\right)+2 \sqrt{\alpha } D \zeta _n^2\right)}{C^{\frac{3}{2}} \sqrt{\frac{-\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\\ \frac{D (1-\nu) \zeta _n \cos \left(y \zeta _n\right) \left(-C^{\frac{3}{2}}+2 \sqrt{\alpha } C+2 \sqrt{\alpha } D \zeta _n^2\right)}{C^{\frac{3}{2}} \sqrt{\frac{-\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\\ \frac{\sin \left(y \zeta _n\right) \left(C^{\frac{3}{2}}-\sqrt{\alpha } C-2 \sqrt{\alpha } D \zeta _n^2\right)}{\sqrt{C} \sqrt{\frac{-\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\\ -\frac{2 \sqrt{\alpha } D \sin \left(y \zeta _n\right) \left(C+(D-D \nu ) \zeta _n^2\right)}{C^{\frac{3}{2}}}\\ \frac{2 \sqrt{\alpha } D \zeta _n \cos \left(y \zeta _n\right)}{C^{\frac{3}{2}}}\\ 0\\ \end{array} \right), {\textbf{Y}}_{n,4}^1 = \left( \begin{array}{c} \frac{\sin \left(y \zeta _n\right) \left(C^{\frac{3}{2}}-2 \sqrt{\alpha } C-2 \sqrt{\alpha } D \zeta _n^2\right)}{C^{\frac{3}{2}} \sqrt{\frac{-\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{d}}}\\ \frac{D (\nu -1) \zeta _n \cos \left(y \zeta _n\right) \left(-C^{\frac{3}{2}}+2 \sqrt{\alpha } C+2 \sqrt{\alpha } D \zeta _n^2\right)}{C^{\frac{3}{2}} \sqrt{\frac{-\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\\ \frac{\sin \left(y \zeta _n\right) \left(-C\left(-\sqrt{\alpha }+\sqrt{C}\right)+2 \sqrt{\alpha } D \zeta _n^2\right)}{\sqrt{C} \sqrt{\frac{-\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\\ -\frac{2 \sqrt{\alpha } D \sin \left(y \zeta _n\right) \left(C+(D-D \nu ) \zeta _n^2\right)}{C^{\frac{3}{2}}}\\ \frac{2 \sqrt{\alpha } D \zeta _n \cos \left(y \zeta _n\right)}{C^{\frac{3}{2}}}\\ 0\\ \end{array} \right). \end{align*}$

2.3. Properties of eigenvectors

After calculating the eigenvalues and eigenvectors, we continue to study the properties of the eigenvectors. And for the analytical solution of this kind of moderately thick rectangular plates, the choice of space is very important, and here we choose the Hilbert space $\mathcal{Y}: = L^{2}(0, b)\oplus L^{2}(0, b)\oplus L^{2}(0, b) \oplus L^{2}(0, b)\oplus L^{2}(0, b)\oplus L^{2}(0, b)$ , where $L^{2}(0, b)$ is the space of all square integrable functions. For $L(y), R(y)\in \mathcal{Y}$ , the symplectic inner product is defined as

$\left\langle L(y), R(y) \right \rangle = \int_{0}^{b}L(y)^{T}{\textbf{J}} R(y)\, {\mathrm{d}}y,$

where ${\textbf{J}}$ is the $6\times 6$ symplectic unit matrix.

2.3.1. When the eigenvalues are simple

Symplectic orthogonality of eigenvectors lays a foundation for the expansion of eigenvectors. The following lemma can be obtained via straightforward calculation.

Lemma 1. The system of eigenvectors $\{{\boldsymbol{Y}}_{1}(y), {\boldsymbol{Y}}_{2}(y)\}\cup\{{\boldsymbol{Y}}_{n, i}(y)\mid n\in \mathbb{N}, i = 1, 2, \cdots, 6\}$ of ${\boldsymbol{H}}$ satisfies the following symplectic orthogonality relations:

$\left\langle {\boldsymbol{Y}}_{1}(y), {\boldsymbol{Y}}_{2}(y)\right \rangle = \left\{\begin{array}{cc} \frac{2 \sqrt{2} b (-C)^{\frac{3}{2}}}{\sqrt{D (\nu-1)}},&n = m,\\ 0,& n\neq m; \end{array}\right.$

$\left\langle {\boldsymbol{Y}}_{n,3}(y), {\boldsymbol{Y}}_{m,4}(y)\right \rangle = \left\{\begin{array}{cc} A_{34,n}, &n = m, \\ 0, &n\neq m; \end{array}\right.$

$\left\langle {\boldsymbol{Y}}_{n,5}(y), {\boldsymbol{Y}}_{m,6}(y)\right \rangle = \left\{\begin{array}{cc} A_{56,n}, &n = m, \\ 0, &n\neq m; \end{array}\right.$

$\left\langle {\boldsymbol{Y}}_{n,1}(y), {\boldsymbol{Y}}_{m,2}(y)\right \rangle = \left\{\begin{array}{cc} \frac{-b C \sqrt{D (\nu -1) \zeta _n^2-2 C}}{\sqrt{D (\nu -1)} \zeta _n^2},&n = m,\\ 0,& n\neq m, \end{array}\right.$

where the constants $A_{34, n}$ and $A_{56, n}$ are presented in Appendix A. Other eigenvectors satisfy the following symplectic orthogonal relation:

$\begin{align*} &\left\langle {\boldsymbol{Y}}_{n,1}(y), {\boldsymbol{Y}}_{m,3}(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,1}(y), {\boldsymbol{Y}}_{m,4}(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,1}(y), {\boldsymbol{Y}}_{m,5}(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,1}(y), {\boldsymbol{Y}}_{m,6}(y) \right \rangle\\ & = \left\langle {\boldsymbol{Y}}_{n,2}(y), {\boldsymbol{Y}}_{m,3}(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,2}(y), {\boldsymbol{Y}}_{m,4}(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,2}(y), {\boldsymbol{Y}}_{m,5}(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,2}(y), {\boldsymbol{Y}}_{m,6}(y) \right \rangle\\ & = \left\langle {\boldsymbol{Y}}_{n,3}(y), {\boldsymbol{Y}}_{m,5}(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,3}(y), {\boldsymbol{Y}}_{m,6}(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,4}(y), {\boldsymbol{Y}}_{m,5}(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,4}(y), {\boldsymbol{Y}}_{m,6}(y) \right \rangle\\ & = \left\langle {\boldsymbol{Y}}_{n,i}(y), {\boldsymbol{Y}}_{m,i}(y) \right \rangle = 0, n, m\in \mathbb{N}. \end{align*}$

Based on the conjugate symplectic orthogonality of eigenvectors, we can obtain the following completeness theorem of eigenvectors, which provides the theoretical basis for the expression of general solutions.

Theorem 1. The system of eigenvectors $\{{\boldsymbol{Y}}_{1}(y), {\boldsymbol{Y}}_{2}(y)\}\cup\{{\boldsymbol{Y}}_{n, i}(y)\}_{n = 1}^{+\infty}$ $(i = 1, 2, \cdots, 6)$ of ${\boldsymbol{H}}$ is complete in $\mathcal{Y}$ .

Proof. For any ${\textbf{F}}(y) = (f_{1}(y), f_{2}(y), f_{3}(y), f_{4}(y), f_{5}(y), f_{6}(y))^{T}\in \mathcal{Y}$ , it is sufficient to prove that there exist constant sequences $\{C_{n, i}\}_{n = 1}^{+\infty}$ , $(i = 1, 2, \cdots, 6)$ , $C_{1}$ and $C_{2}$ such that

$\begin{align} {\textbf{F}}(y) = C_1{\textbf{Y}}_1+C_2{\textbf{Y}}_2+\sum\limits_{n = 1}^{+\infty}(\sum\limits_{i = 1}^{6} C_{n,i}{\textbf{Y}}_{n,i}(y)). \end{align}$

(2.9)

Applying Lemma 1, we can take

$\begin{align} \begin{split} &C_{n,1} = \frac{\left\langle {\textbf{Y}}_{n,2}(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{n,2}(y), {\textbf{Y}}_{n,1}(y) \right \rangle}, C_{n,2} = \frac{\left\langle {\textbf{Y}}_{n,1}(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{n,1}(y), {\textbf{Y}}_{n,2}(y) \right \rangle}, C_{n,3} = \frac{\left\langle {\textbf{Y}}_{n,4}(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{n,4}(y), {\textbf{Y}}_{n,3}(y) \right \rangle},\\ &C_{n,4} = \frac{\left\langle {\textbf{Y}}_{n,3}(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{n,3}(y), {\textbf{Y}}_{n,4}(y) \right \rangle}, C_{n,5} = \frac{\left\langle {\textbf{Y}}_{n,6}(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{n,6}(y), {\textbf{Y}}_{n,5}(y) \right \rangle}, C_{n,6} = \frac{\left\langle {\textbf{Y}}_{n,5}(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{n,5}(y), {\textbf{Y}}_{n,6}(y) \right \rangle},\\ &C_{1} = \frac{\left\langle {\textbf{Y}}_{2}(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{2}(y), {\textbf{Y}}_{1}(y) \right \rangle}, C_{2} = \frac{\left\langle {\textbf{Y}}_{1}(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{1}(y), {\textbf{Y}}_{2}(y) \right \rangle}, \end{split} \end{align}$

(2.10)

and then

$\begin{align*} &C_1{\textbf{Y}}_1+C_2{\textbf{Y}}_2+\sum\limits_{n = 1}^{+\infty}(\sum\limits_{i = 1}^{6} C_{n,i}{\textbf{Y}}_{n,i}(y))\\ & = \left( \begin{array}{c} 0 \\ \frac{1}{b}\int_{0}^{b}f_2(y){\mathrm{d}}y \\ 0 \\ 0 \\ \frac{1}{b}\int_{0}^{b}f_5(y){\mathrm{d}}y \\ 0 \\ \end{array} \right)+\sum\limits_{n = 1}^{+\infty}\left( \begin{array}{c} \frac{2}{b}(\int_{0}^{b}f_{1}(y)\sin\frac{n\pi y}{b}{\mathrm{d}}y)\sin\frac{n\pi y}{b}\\ \frac{2}{b}(\int_{0}^{b}f_{2}(y)\cos\frac{n\pi y}{b}{\mathrm{d}}y)\cos\frac{n\pi y}{b}\\ \frac{2}{b}(\int_{0}^{b}f_{3}(y)\sin\frac{n\pi y}{b}{\mathrm{d}}y)\sin\frac{n\pi y}{b}\\ \frac{2}{b}(\int_{0}^{b}f_{4}(y)\sin\frac{n\pi y}{b}{\mathrm{d}}y)\sin\frac{n\pi y}{b}\\ \frac{2}{b}(\int_{0}^{b}f_{5}(y)\cos\frac{n\pi y}{b}{\mathrm{d}}y)\cos\frac{n\pi y}{b}\\ \frac{2}{b}(\int_{0}^{b}f_{6}(y)\sin\frac{n\pi y}{b}{\mathrm{d}}y)\sin\frac{n\pi y}{b} \end{array} \right). \end{align*}$

Clearly, \, $\sum\limits_{n = 1}^{+\infty} \frac{2}{b}(\int_{0}^{b}f_{i}(y)\sin\frac{n\pi y}{b}{\mathrm{d}}y)\sin\frac{n\pi y}{b} \ (i = 1, 3, 4, 6)$ and

$\frac{1}{b}\int_{0}^{b}f_i(y){\mathrm{d}}y+\sum\limits_{n = 1}^{+\infty} \frac{2}{b}(\int_{0}^{b}f_{i}(y)\cos\frac{n\pi y}{b}{\mathrm{d}}y)\cos\frac{n\pi y}{b} \, (i = 2, 5)$

are the corresponding Fourier series expansions of $f_{i}(y)$ based on the orthogonal systems $\{\sin\frac{n\pi y}{b}\}_{n = 1}^{+\infty}$ and $\{\cos\frac{n\pi y}{b}\}_{n = 0}^{+\infty}$ in $L^{2}(0, b)$ , respectively. Thus, ${\textbf{F}}(y)$ can be expressed by (2.9) and the proof is completed. □

2.3.2. When the eigenvalues are double

Similar to the case where the eigenvalues are single, we first discuss the orthogonal conjugations of eigenvectors. These relations play a crucial role in proving the symplectic eigenvectors expansion theorem.

Lemma 2. The system of eigenvectors $\{{\boldsymbol{Y}}_1(y), {\boldsymbol{Y}}_2(y)\}\cup\{{\boldsymbol{Y}}_{n, i}^0 (y), {\boldsymbol{Y}}_{n, j}^1 (y)\mid n\in \mathbb{N}; i = 1, 2, 3, 4; j = 3, 4\}$ of ${\boldsymbol{H}}$ have the following symplectic orthogonality relations:

$\begin{align*} &\left\langle {\boldsymbol{Y}}_{n,1}^0(y), {\boldsymbol{Y}}_{m,3}^0(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,1}^0(y), {\boldsymbol{Y}}_{m,3}^1(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,1}^0(y), {\boldsymbol{Y}}_{m,4}^0(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,1}^0(y), {\boldsymbol{Y}}_{m,4}^1(y) \right \rangle\\ & = \left\langle {\boldsymbol{Y}}_{n,2}^0(y), {\boldsymbol{Y}}_{m,3}^0(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,2}^0(y), {\boldsymbol{Y}}_{m,3}^1(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,2}^0(y), {\boldsymbol{Y}}_{m,4}^0(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,2}^0(y), {\boldsymbol{Y}}_{m,4}^1(y) \right \rangle \\ & = \left\langle {\boldsymbol{Y}}_{n,3}^0(y), {\boldsymbol{Y}}_{m,3}^1(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,3}^0(y), {\boldsymbol{Y}}_{m,4}^0(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,4}^0(y), {\boldsymbol{Y}}_{m,4}^1(y) \right \rangle = \left\langle {\boldsymbol{Y}}_{n,i}^0(y), {\boldsymbol{Y}}_{m,i}^0(y) \right \rangle \\ & = 0,\, \, \,\, \, n, m\in \mathbb{N}. \end{align*}$

When $\gamma = \gamma_{n1} = \frac{-2 \sqrt{\alpha } C^{\frac{3}{2}}-C^2-\alpha C-\alpha D \zeta _n^2+\beta D \zeta _n^2}{D}$ , we have

$\left\langle {\boldsymbol{Y}}_{1}(y), {\boldsymbol{Y}}_{2}(y) \right \rangle = \left\{\begin{array}{cc} \frac{2 \sqrt{2} b {(-C)}^{\frac{3}{2}} }{\sqrt{D (\nu-1)}}, &n = m,\\ 0, &n\neq m; \end{array}\right.$

$\left\langle {\boldsymbol{Y}}_{n,3}^0(y), {\boldsymbol{Y}}_{m,4}^1(y) \right \rangle = \left\{\begin{array}{cc} \frac{2 \sqrt{\alpha } b D}{\sqrt{C}}, & n = m, \\ 0, &n\neq m; \end{array}\right.$

$\left\langle {\boldsymbol{Y}}_{n,3}^1(y), {\boldsymbol{Y}}_{m,4}^0(y) \right \rangle = \left\{\begin{array}{cc} -\frac{2 \sqrt{\alpha } b D}{\sqrt{C}}, & n = m,\\ 0,& n\neq m; \end{array}\right.$

$\left\langle {\boldsymbol{Y}}_{n,1}^0(y), {\boldsymbol{Y}}_{m,2}^0(y) \right \rangle = \left\{\begin{array}{cc} \frac{-b C \sqrt{D (\nu -1) \zeta _n^2-2 C}}{\sqrt{D (\nu -1)} \zeta _n^2},& n = m,\\ 0,& n\neq m; \end{array}\right.$

$\left\langle {\boldsymbol{Y}}_{n,3}^1(y), {\boldsymbol{Y}}_{m,4}^1(y) \right \rangle = \left\{\begin{array}{cc} -\frac{2 b \left(\sqrt{\alpha } C^{\frac{3}{2}} D+2 \alpha C D+2 \alpha D^2 \zeta _n^2\right)}{C^2 \sqrt{\frac{\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}, &n = m, \\ 0, &n\neq m. \end{array}\right.$

If $\gamma = \gamma_{n2} = \frac{2 \sqrt{\alpha } C^{\frac{3}{2}}-C^2-\alpha C-\alpha D \zeta _n^2+\beta D \zeta _n^2}{D}$ , then

$\left\langle {\boldsymbol{Y}}_{1}(y), {\boldsymbol{Y}}_{2}(y) \right \rangle = \left\{\begin{array}{cc} \frac{2 \sqrt{2} b {(-C)}^{\frac{3}{2}}}{\sqrt{D (\nu-1)}},&n = m,\\ 0,&n\neq m; \end{array}\right.$

$\left\langle {\boldsymbol{Y}}_{n,3}^0(y), {\boldsymbol{Y}}_{m,4}^1(y) \right \rangle = \left\{\begin{array}{cc} -\frac{2 \sqrt{\alpha } b D}{\sqrt{C}}, &n = m, \\ 0, &n\neq m; \end{array}\right.$

$\left\langle {\boldsymbol{Y}}_{n,3}^1(y), {\boldsymbol{Y}}_{m,4}^0(y) \right \rangle = \left\{\begin{array}{cc} \frac{2 \sqrt{\alpha } b D}{\sqrt{C}},&n = m,\\ 0,&n\neq m; \end{array}\right.$

$\left\langle {\boldsymbol{Y}}_{n,3}^1(y), {\boldsymbol{Y}}_{m,4}^1(y) \right \rangle = \left\{\begin{array}{cc} \frac{2 b \left(\sqrt{\alpha } C^{\frac{3}{2}} D-2 \alpha C D-2 \alpha D^2 \zeta _n^2\right)}{C^2 \sqrt{\frac{-\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}, &n = m, \\ 0, &n\neq m. \end{array}\right.$

Based on Lemma 2, the following completeness theorem is obtained.

Theorem 2. The system of eigenvectors $\{{\boldsymbol{Y}}_{1}(y), {\boldsymbol{Y}}_{2}(y)\}\cup\{{\boldsymbol{Y}}_{n, i}^0(y)\}\cup\{{\boldsymbol{Y}}_{n, j}^1(y)\}$ , $(i = 1, 2, 3, 4; j = 3, 4; n\in \mathbb{N})$ of ${\boldsymbol{H}}$ is complete in $\mathcal{Y}$ .

Proof. It is similar to the proof of Theorem 1 that we just have to take the sequences:

$\begin{align} \begin{split} &C_{n,1}^0 = \frac{\left\langle {\textbf{Y}}_{n,2}^0(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{n,2}^0(y), {\textbf{Y}}_{n,1}^0(y) \right \rangle}, C_{n,2}^0 = \frac{\left\langle {\textbf{Y}}_{n,1}^0(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{n,1}^0(y), {\textbf{Y}}_{n,2}^0(y) \right \rangle}, C_{n,3}^1 = \frac{\left\langle {\textbf{Y}}_{n,4}^0(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{n,4}^0(y), {\textbf{Y}}_{n,3}^1(y) \right \rangle},\\ &C_{n,4}^1 = \frac{\left\langle {\textbf{Y}}_{n,3}^0(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{n,3}^0(y), {\textbf{Y}}_{n,4}^1(y) \right \rangle}, C_{1} = \frac{\left\langle {\textbf{Y}}_{2}(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{2}(y), {\textbf{Y}}_{1}(y) \right \rangle}, C_{2} = \frac{\left\langle {\textbf{Y}}_{1}(y),{\textbf{F}}(y) \right \rangle}{\left\langle {\textbf{Y}}_{1}(y), {\textbf{Y}}_{2}(y) \right \rangle},\\ &C_{n,3}^0 = \frac{\left\langle {\textbf{Y}}_{n,4}^1(y),{\textbf{F}}(y) \right \rangle-C_{n,3}^1\left\langle {\textbf{Y}}_{n,4}^1(y),{\textbf{Y}}_{n,3}^1(y) \right \rangle}{\left\langle {\textbf{Y}}_{n,4}^1(y), {\textbf{Y}}_{n,3}^0(y) \right \rangle}, C_{n,4}^0 = \frac{\left\langle {\textbf{Y}}_{n,3}^1(y),{\textbf{F}}(y) \right \rangle-C_{n,4}^1\left\langle {\textbf{Y}}_{n,3}^1(y),{\textbf{Y}}_{n,4}^1(y) \right \rangle}{\left\langle {\textbf{Y}}_{n,3}^1(y), {\textbf{Y}}_{n,4}^0(y) \right \rangle}. \end{split} \end{align}$

(2.11)

It is easy to calculate that

$\begin{align*} &C_1{\textbf{Y}}_1+C_2{\textbf{Y}}_2+\sum\limits_{n = 1}^{+\infty}(\sum\limits_{i = 1}^{4} C_{n,i}^0{\textbf{Y}}_{n,i}^0(y)+\sum\limits_{j = 3}^{4} C_{n,j}^1{\textbf{Y}}_{n,j}^1(y))\\ & = \left( \begin{array}{c} 0 \\ \frac{1}{b}\int_{0}^{b}f_2(y){\mathrm{d}}y \\ 0 \\ 0 \\ \frac{1}{b}\int_{0}^{b}f_5(y){\mathrm{d}}y \\ 0 \\ \end{array} \right)+\sum\limits_{n = 1}^{+\infty}\left( \begin{array}{c} \frac{2}{b}(\int_{0}^{b}f_{1}(y)\sin\frac{n\pi y}{b}{\mathrm{d}}y)\sin\frac{n\pi y}{b}\\ \frac{2}{b}(\int_{0}^{b}f_{2}(y)\cos\frac{n\pi y}{b}{\mathrm{d}}y)\cos\frac{n\pi y}{b}\\ \frac{2}{b}(\int_{0}^{b}f_{3}(y)\sin\frac{n\pi y}{b}{\mathrm{d}}y)\sin\frac{n\pi y}{b}\\ \frac{2}{b}(\int_{0}^{b}f_{4}(y)\sin\frac{n\pi y}{b}{\mathrm{d}}y)\sin\frac{n\pi y}{b}\\ \frac{2}{b}(\int_{0}^{b}f_{5}(y)\cos\frac{n\pi y}{b}{\mathrm{d}}y)\cos\frac{n\pi y}{b}\\ \frac{2}{b}(\int_{0}^{b}f_{6}(y)\sin\frac{n\pi y}{b}{\mathrm{d}}y)\sin\frac{n\pi y}{b} \end{array} \right). \end{align*}$

As a result,

$\begin{align} {\textbf{F}}(y) = C_1{\textbf{Y}}_1+C_2{\textbf{Y}}_2+\sum\limits_{n = 1}^{+\infty}(\sum\limits_{i = 1}^{4} C_{n,i}^0{\textbf{Y}}_{n,i}^0(y)+\sum\limits_{j = 3}^{4} C_{n,j}^1{\textbf{Y}}_{n,j}^1(y)). \end{align}$

(2.12)

□

2.4. General solution

This section will give the general solutions of the Hamiltonian system under simply supported boundary conditions at $y = 0$ and $y = b$ based on the proven completeness theorems.

2.4.1. When eigenvalues are simple

The general solution of (2.5) can be represented in the form:

$\begin{align} {\textbf{U}}(x,y) = X_1(x){\textbf{Y}}_1(y)+X_2(x){\textbf{Y}}_2(y)+\sum\limits_{n = 1}^{+\infty}(\sum\limits_{i = 1}^{6}( X_{n,i}(x){\textbf{Y}}_{n,i}(y))). \end{align}$

(2.13)

The nonhomogeneous term ${\bf{f}}(x, y)$ expands to

$\begin{align} {\bf{f}}(x, y) = \theta_1(x){\textbf{Y}}_1(y)+\theta_2(x){\textbf{Y}}_2(y)+\sum\limits_{n = 1}^{+\infty}(\sum\limits_{i = 1}^{6}( \theta_{n,i}(x){\textbf{Y}}_{n,i}(y))). \end{align}$

(2.14)

From the above formulas and (2.6), we obtain

$\begin{align} &X_{1}(x) = \Psi_{1} e^{\mu_{1} x}+\int_0^x\theta_{1}(t)e^{\mu_{1}(x-t)}{\mathrm{d}}t, \\ &X_{2}(x) = \Psi_{2} e^{\mu_{2} x}+\int_0^x\theta_{2}(t)e^{\mu_{2}(x-t)}{\mathrm{d}}t, \\ &X_{n,i}(x) = \Psi_{ni}e^{\mu_{ni} x}+\int_{0}^{x}\theta_{n,i}(t)e^{\mu_{ni}(x-t)}{\mathrm{d}}t, \end{align}$

(2.15)

where $\Psi_{1}, \Psi_{2}, \Psi_{ni}$ are undetermined constants, $\theta_{1}, \theta_{2}, \theta_{n, i}$ can be expressed in a form similar to (2.10) via the symplectic inner product of ${\bf{f}}(x, y)$ and $Y$ . In particular, if $q_{k}(x, y)$ $(k = 1, 2, 3)$ is the uniformly distributed load acting on the plate, which means that it is a constant independent of $x$ and of $y$ . If $q_{k}(x, y)$ $(k = 1, 2, 3)$ represents a concentrated load of strength $B$ acting on the plate, then it can be expressed as $B\delta(x-x_0, y-y_0)$ in which $\delta(x-x_0, y-y_0)$ denotes the $\delta$ function.

The general solutions of $\psi_{y}$ and $w$ can be obtained by taking the fifth and sixth components of (2.13), respectively.

$\begin{align} \psi_{y} = &\left[\left(\Psi_1 e^{\mu_1 x}+\int _0^x\int _0^y\frac{q_2(x,y)}{2 b C}e^{\mu_1 (x-t)}{\mathrm{d}}y{\mathrm{d}}t\right)-\left(\Psi_2 e^{\mu_2 x}+\int _0^x\int _0^y\frac{q_2(x,y)}{2 b C}e^{\mu_2 (x-t)}{\mathrm{d}}y{\mathrm{d}}t\right)\right]\mu_1\\ &+\sum\limits_{n = 1}^{+\infty}\left\{\left[\left(\int_0^x e^{\mu_{n1} (x-t)} \theta _{n,1}(t) \, {\mathrm{d}}t+\Psi_{n1} e^{\mu_{n1} x}\right)-\left(\int_0^x e^{-\mu_{n1} (x-t)} \theta _{n,2}(t) \, {\mathrm{d}}t+\Psi_{n2} e^{-\mu_{n1} x}\right)\right]\right.\\ &\frac{\cos \left(y \zeta _n\right) \left[2 C^2 \mu_{n1}^2-2 C \left(\alpha \mu_{n1}^2+\gamma -\beta \zeta _n^2\right)+D \left(2 \mu_{n1}^2+(\nu -1) \zeta _n^2\right) \left(\alpha \mu_{n1}^2+\gamma -\beta \zeta _n^2\right)\right]}{ \mu_{n1} \zeta _n \left(2 C^2+D (\nu +1) \left(\alpha \mu_{n1}^2+\gamma -\beta \zeta _n^2\right)\right)}\\ &+\left[\left(\int_0^x e^{\mu_{n3} (x-t)} \theta _{n,3}(t) \, {\mathrm{d}}t+e^{\mu_{n3} x} \Psi_{n3}\right)-\left(\int_0^x e^{-\mu_{n3} (x-t)} \theta _{n,4}(t) \, {\mathrm{d}}t+\Psi_{n4} e^{-\mu_{n3} x}\right)\right]\\ &\frac{\cos \left(y \zeta _n\right) \left[2 C^2 \mu_{n3}^2-2 C \left(\alpha \mu_{n3}^2+\gamma -\beta \zeta _n^2\right)+D \left(2 \mu_{n3}^2+(\nu -1) \zeta _n^2\right) \left(\alpha \mu_{n3}^2+\gamma -\beta \zeta _n^2\right)\right]}{\mu_{n3} \zeta _n \left(2 C^2+D (\nu +1) \left(\alpha \mu_{n3}^2+\gamma -\beta \zeta _n^2\right)\right)}\\ &+\left[\left(\int_0^x e^{\mu_{n5} (x-t)} \theta _{n,5}(t) \, {\mathrm{d}}t+\Psi_{n5} e^{\mu_{n5} x}\right)-\left(\int_0^x e^{-\mu_{n5} (x-t)} \theta _{n,6}(t) \, {\mathrm{d}}t+\Psi_{n6} e^{-\mu_{n5} x}\right)\right]\\ &\left.\frac{\cos \left(y \zeta _n\right) \left[2 C^2 \mu_{n5}^2-2 C \left(\alpha \mu_{n5}^2+\gamma -\beta \zeta _n^2\right)+D \left(2 \mu_{n5}^2+(\nu -1) \zeta _n^2\right) \left(\alpha \mu_{n5}^2+\gamma -\beta \zeta _n^2\right)\right]}{\mu_{n5} \zeta _n \left(2 C^2+D (\nu +1) \left(\alpha \mu_{n5}^2+\gamma -\beta \zeta _n^2\right)\right)}\right\}, \end{align}$

(2.16)

$\begin{align} w = &\sum\limits_{n = 1}^{+\infty}\left\{\left[\left(\int_0^x e^{\mu_{n3}(x-t)} \Psi_{n,3}(t) \, dt+\Psi_{n3} e^{\mu_{n3} x}\right)-\left(\int_0^x e^{-\mu_{n3} (x-t)} \theta _{n,4}(t) \, {\mathrm{d}}t+\Psi_{n4} e^{-\mu_{n3} x}\right)\right]\right.\\ &\frac{C \sin \left(y \zeta _n\right) \left(2 C+D (\nu-1) \left(\mu_{n3}^2-\zeta _n^2\right)\right)}{\mu_{n3} \left(2 C^2+D (\nu+1) \left(\alpha \mu_{n3}^2+\gamma -\beta \zeta _n^2\right)\right)}+\frac{C \sin \left(y \zeta _n\right) \left(2 C+D (\nu-1) \left(\mu_{n5}^2-\zeta _n^2\right)\right)}{\mu_{n5} \left(2 C^2+D (\nu+1) \left(\alpha \mu_{n5}^2+\gamma -\beta \zeta _n^2\right)\right)}\times\\ &\left.\left[\left(\int_0^x e^{\mu_{n5} (x-t)} \theta _{n,5}(t) \, dt+\Psi_{n5} e^{\mu_{n5} x}\right)-\left(\int_0^x e^{-\mu_{n5} (x-t)} \theta _{n,6}(t) \, {\mathrm{d}}t+\Psi_{n6} e^{-\mu_{n5} x}\right)\right]\right\}. \end{align}$

(2.17)

2.4.2. When eigenvalues are double

Compared to the case where the eigenvalues are simple, we combine (2.6) with (2.12) to obtain:

$\begin{align*} &X_{1}(x) = Q_{1} e^{\mu_{1} x}+\int_0^x\Theta_{1}(t)e^{\mu_{1}(x-t)}{\mathrm{d}}t, \\ &X_{2}(x) = Q_{2} e^{\mu_{2} x}+\int_0^x\Theta_{2}(t)e^{\mu_{2}(x-t)}{\mathrm{d}}t,\\ &X_{n,1}^0(x) = Q_{n,1}^0 e^{\mu_{n1}(x)}+\int_{0}^{x}\Theta_{n,1}^0(t)e^{\mu_{n1}(x-t)}{\mathrm{d}}t,\\ &X_{n,2}^0(x) = Q_{n,2}^0 e^{\mu_{n2}(x)}+\int_{0}^{x}\Theta_{n,2}^0(t)e^{\mu_{n2}(x-t)}{\mathrm{d}}t,\\ &X_{n,3}^1(x) = Q_{n,3}^1 e^{\mu_{n3}(x)}+\int_{0}^{x}\Theta_{n,3}^0(t)e^{\mu_{n3}(x-t)}{\mathrm{d}}t,\\ &X_{n,4}^1(x) = Q_{n,4}^1 e^{\mu_{n4}(x)}+\int_{0}^{x}\Theta_{n,4}^0(t)e^{\mu_{n4}(x-t)}{\mathrm{d}}t,\\ &X_{n,3}^0(x) = (Q_{n,3}^0+Q_{n,3}^1x)e^{\mu_{n3}(x)}+\int_{0}^{x}\Theta_{n,3}^0(t)e^{\mu_{n3}(x-t)}{\mathrm{d}}t+\int_{0}^{x}\int_{0}^{t}\Theta_{n,3}^1(t)e^{\mu_{n3}(x-\tau)}{\mathrm{d}}\tau{\mathrm{d}}t,\nonumber\\ &X_{n,4}^0(x) = (Q_{n,4}^0+Q_{n,4}^1x)e^{\mu_{n4}(x)}+\int_{0}^{x}\Theta_{n,4}^0(t)e^{\mu_{n4}(x-t)}{\mathrm{d}}t+\int_{0}^{x}\int_{0}^{t}\Theta_{n,4}^1(t)e^{\mu_{n4}(x-\tau)}{\mathrm{d}}\tau{\mathrm{d}}t, \end{align*}$

where $Q_{1}, Q_{2}, Q_{n, i}^0, Q_{n, j}^1$ $(i = 1, 2, 3, 4; \ j = 3, 4)$ are undetermined constants, $\Theta_{1, 2}, \Theta_{n, i}^0, \Theta_{n, j}^1$ can be computed similarly as (2.11) by symplectic inner product between ${\bf{f}}(x, y)$ and the eigenvectors.

Then, ${\textbf{U}}(x, y)$ can be given by

$\begin{align} {\textbf{U}}(x,y) = X_1(x){\textbf{Y}}_1(y)+X_2(x){\textbf{Y}}_2(y)+ \sum\limits_{n = 1}^{+\infty}(\sum\limits_{i = 1}^{4}( X_{n,i}^0(x){\textbf{Y}}_{n,i}^0(y))+\sum\limits_{j = 3}^{4}( X_{n,j}^1(x){\textbf{Y}}_{n,j}^1(y))). \end{align}$

(2.18)

As for the case of double eigenvalues, in order to obtain a simple form of $\psi_{y}$ and $w$ , a classification discussion is carried out.

(i) When $\gamma = \gamma_{n1} = \frac{-2 \sqrt{\alpha } C^{\frac{3}{2}}-C^2-\alpha C-\alpha D \zeta _n^2+\beta D \zeta _n^2}{D},$ we have

$\begin{align*} \psi_{y} = &\left[\left(Q_1 e^{\mu_1 x}+\int _0^x\Theta _{1}(t)e^{\mu_1 (x-t)}{\mathrm{d}}t\right)-\left(Q_2 e^{\mu_2 x}+\int _0^x\Theta _{2}(t)e^{-\mu_1 (x-t)}{\mathrm{d}}t\right)\right]\frac{\sqrt{-2 C}}{\sqrt{D (\nu -1)}}\nonumber\\ &+\sum\limits_{n = 1}^{+\infty}\left\{\left[e^{\mu_{n1} x} Q_{n,1}^0+\int_0^x e^{\mu_{n1} (x-t)} \Theta _{n,1}^0(t) \, {\mathrm{d}}t-\left(e^{\mu_{n2} x} Q_{n,2}^0+\int_0^x e^{\mu_{n2} (x-t)} \Theta _{n,2}^0(t) \, {\mathrm{d}}t\right)\right]\right.\nonumber\\ &\left(\frac{\cos \left(y \zeta _n\right) \sqrt{D (\nu -1) \zeta _n^2-2 C}}{\sqrt{D (\nu -1)} \zeta _n}\right)+\left[\left(e^{\mu_{n3} x} \left(x Q_{n,3}^1+Q_{n,3}^0\right)+\int_0^x e^{\mu_{n3} (x-t)} \Theta _{n,3}^0(t) \, {\mathrm{d}}t+\right.\right.\nonumber\\ &\left.\int _0^x\int _0^te^{\mu_{n3} (x-\tau )} \Theta _{n,3}^1(t){\mathrm{d}}\tau {\mathrm{d}}t\right)-\left(e^{\mu_{n4} x} \left(x Q_{n,4}^1+Q_{n,4}^0\right)+\int_0^x \int _0^te^{\mu_{n4}(x-\tau )} \Theta _{n,4}^1(t) \,{\mathrm{d}}\tau {\mathrm{d}}t\right.\nonumber\\ &+\left.\left.\int _0^x e^{\mu_{n4}(x-t) } \Theta _{n,4}^9(t){\mathrm{d}}t\right)\right]\left(\frac{\zeta _n \cos \left(y \zeta _n\right)}{\sqrt{\frac{\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\right)+\left[(e^{\mu_{n3} x} Q_{n,3}^1+\int_0^x e^{\mu_{n3} (x-t)} \Theta _{n,3}^0(t) \, {\mathrm{d}}t)\right.\nonumber\\ &\left.\left.+(e^{\mu_{n4} x} Q_{n,4}^1+\int_0^x e^{\mu_{n4} (x-t)} \Theta _{n,4}^0(t) \, {\mathrm{d}}t)\right]\left(-\frac{2 \sqrt{\alpha } D \zeta _n \cos \left(y \zeta _n\right)}{C^{\frac{3}{2}}}\right)\right\}, \end{align*}$

$\begin{align*} w = &\sum\limits_{n = 1}^{+\infty}\left\{\left[e^{\mu_{n3} x} \left(x Q_{n,3}^1+Q_{n,3}^0\right)+\int _0^x\int _0^te^{\mu_{n3} (x-\tau )} \Theta _{n,3}^1(t){\mathrm{d}}\tau {\mathrm{d}}t+\int_0^x e^{\mu_{n3} (x-t)} \Theta _{n,3}^0(t) \, {\mathrm{d}}t\right]\right.\nonumber\\ &-\left[e^{\mu_{n4} x} \left(x Q_{n,4}^1+Q_{n,4}^0\right)+\int _0^x\int _0^te^{\mu_{n4} (x-\tau )} \Theta _{n,4}^1(t){\mathrm{d}}\tau {\mathrm{d}}t+\int_0^x e^{\mu_{n4} (x-t)} \Theta _{n,4}^0(t) \, {\mathrm{d}}t\right]\\ &\times\nonumber\left.\left(-\frac{\sqrt{C} \sin \left(y \zeta _n\right)}{\sqrt{\alpha } \sqrt{\frac{\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\right)\right\}. \end{align*}$

(ii) When $\gamma = \gamma_{n2} = \frac{2 \sqrt{\alpha } C^{\frac{3}{2}}-C^2-\alpha C-\alpha D \zeta _n^2+\beta D \zeta _n^2}{D},$ we can see that

$\begin{align*} \psi_{y} = &\left[\left(Q_1 e^{\mu_1 x}+\int_0^x \Theta _1(t) e^{\mu_1 (x-t)} \, dt\right)-\left(Q_2 e^{\mu_2 x}+\int_0^x \Theta _2(t) e^{\mu_2 (x-t)} \, dt\right)\right]\frac{\sqrt{-2 C}}{\sqrt{D (\nu -1)}}\nonumber\\ &+\sum\limits_{n = 1}^{+\infty}\left\{\left[e^{\mu_{n1} x} Q_{n,1}^0+\int_0^x e^{\mu_{n1} (x-t)} \Theta _{n,1}^0(t) \, {\mathrm{d}}t-\left(e^{\mu_{n2} x} Q_{n,2}^0+\int_0^x e^{\mu_{n2} (x-t)} \Theta _{n,2}^0(t) \, {\mathrm{d}}t\right)\right]\right.\nonumber\\ &\left(\frac{\cos \left(y \zeta _n\right) \sqrt{D (\nu -1) \zeta _n^2-2 C}}{\sqrt{D (\nu -1)} \zeta _n}\right)+\left[\left(e^{\mu_{n3} x} \left(x Q_{n,3}^1+Q_{n,3}^0\right)+\int_0^x e^{\mu_{n3} (x-t)} \Theta _{n,3}^0(t) \, {\mathrm{d}}t\right.\right.\nonumber\\ &\left.+\int _0^x\int _0^te^{\mu_{n3} (x-\tau )} \Theta _{n,3}^1(t){\mathrm{d}}\tau {\mathrm{d}}t\right)-\left(e^{\mu_{n4} x} \left(x Q_{n,4}^1+Q_{n,4}^0\right)+\int_0^x e^{\mu_{n4} (x-t)} \Theta _{n,4}^0(t) \, {\mathrm{d}}t\right.\nonumber\\ &\left.\left.+\int _0^x\int _0^te^{\mu_{n4} (x-\tau )} \Theta _{n,4}^1(t){\mathrm{d}}\tau {\mathrm{d}}t\right)\right]\left(\frac{\zeta _n \cos \left(y \zeta _n\right)}{\sqrt{\frac{-\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\right)+\left(\frac{2 \sqrt{\alpha } D \zeta _n \cos \left(y \zeta _n\right)}{C^{\frac{3}{2}}}\right)\\ &\times\nonumber\left.\left[\left(e^{\mu_{n3} x} Q_{n,3}^1+\int_0^x e^{\mu_{n3} (x-t)} \Theta _{n,3}^0(t) \, {\mathrm{d}}t\right)+\left(e^{\mu_{n4} x} Q_{n,4}^1+\int_0^x e^{\mu_{n4} (x-t)} \Theta _{n,4}^0(t) \, {\mathrm{d}}t\right)\right]\right\}, \end{align*}$

$\begin{align*} w = &\sum\limits_{n = 1}^{+\infty}\left\{\left[\left(e^{\mu_{n3} x} \left(x Q_{n,3}^1+Q_{n,3}^0\right)+\int _0^x\int _0^te^{\mu_{n3} (x-\tau )} \Theta _{n,3}^1(t)d\tau dt+\int_0^x e^{\mu_{n3} (x-t)} \Theta _{n,3}^0(t) \, {\mathrm{d}}t\right)\right.\right.\nonumber\\ &\left.-\left(e^{\mu_{n4} x} \left(x Q_{n,4}^1+Q_{n,4}^0\right)+\int _0^x\int _0^te^{\mu_{n4} (x-\tau )} \Theta _{n,4}^1(t){\mathrm{d}}\tau {\mathrm{d}}t+\int_0^x e^{\mu_{n4} (x-t)} \Theta _{n,4}^0(t) \, {\mathrm{d}}t\right)\right]\\ &\times\nonumber\left.\left(\frac{\sqrt{C} \sin \left(y \zeta _n\right)}{\sqrt{\alpha } \sqrt{\frac{-\frac{C^{\frac{3}{2}}}{\sqrt{\alpha }}+C+D \zeta _n^2}{D}}}\right)\right\}. \end{align*}$

3. Applications

The general solution of a kind of rectangular moderately thick plates has been obtained, which can be applied to some specific mechanical problems, such as bending, buckling, free vibration, etc. For example, when $\alpha = \beta = C, \gamma = q_1(x, y) = q_2(x, y) = 0, q_3(x, y) = -q$ , the model reduces to the bending problem of moderately thick plates with two opposite edges simply supported ^[31]; when $\alpha = C+N_{x}, \beta = C+N_{y}, \gamma = q_1(x, y) = q_2(x, y) = q_3(x, y) = 0$ , it represents the document ^[20] in discussing the buckling problem of the moderately thick rectangular plates; the reference [] discusses a free vibration problem of rectangular moderately thick plates, corresponding to $\alpha = \beta = C, \, \gamma = \rho h \varpi^2$ ( $\varpi$ represents circular frequency) and $q_1(x, y) = q_2(x, y) = q_3(x, y) = 0$ . For these specific problems, although they may correspond to different boundary conditions, most of them can be obtained by the general solution in this paper.

So far, SEA has not been applied to the problems of moderately thick rectangular plates with elastic foundations. In this section, the solutions to two elastic foundation problems are discussed first based on the general solutions of the considered model and the numerical results are compared with those in the literature.

Example 1. Consider a uniformly loaded and simply supported square plate resting on a Winkler-type elastic foundation with side length $a$ . In this case, $q_{1}(x, y) = q_{2}(x, y) = 0$ , $q_{3}(x, y) = q$ , $\alpha_1 = \beta_1 = 0$ and $\gamma = k$ , where $q$ denotes the uniform load and $k$ is the subgrade reaction operator for Winkler-type foundation ^[33,34,35], the corresponding equilibrium has the following form,

$\begin{align} \left\{ \begin{array}{cc} \frac{\partial M_{x}}{\partial x}+\frac{\partial M_{xy}}{\partial y}-Q_{x} = 0,\\ \frac{\partial M_{y}}{\partial y}+\frac{\partial M_{xy}}{\partial x}-Q_{y} = 0,\\ \frac{\partial Q_{x}}{\partial x}+\frac{\partial Q_{y}}{\partial y}+kw = q. \end{array}\right. \end{align}$

(3.1)

According to (2.17), the form of $w$ is given in detail:

$\begin{array}{l} w = \lim\limits_{N\rightarrow +\infty}\sum\limits_{n = 1}^{N}\left[\frac{e^{\mu_{n3} (-x)} \left(\left(e^{\mu_{n3} x}-1\right)^2 \theta _{n,3}+\mu_{n3}\Psi_{n3} e^{2 \mu_{n3} x}-\mu_{n3}\Psi_{n4}\right)}{\mu_{n3}}\right.\\\times\frac{\sqrt{2} \sin \left(y \zeta _n\right) \left(4 C^2-D (\nu -1) (k+P_1)\right)}{\left(4 C^2+D (\nu +1) (k-P_1)\right) \sqrt{\frac{2 C \zeta _n^2-k-P_1}{C}}}\nonumber\\ \left.+\frac{e^{\mu_{n5} (-x)} \left(\left(e^{\mu_{n5} x}-1\right)^2 \theta _{n,5}+\mu_{n5} \Psi_{n5} e^{2 \mu_{n5} x}-\mu_{n5}\Psi_{n6}\right)}{\mu_{n5}} \\\times \frac{\sqrt{2} \sin \left(y \zeta _n\right) \left(4 C^2-D (\nu -1) (k-P_1)\right)}{\left(4 C^2+D (\nu +1) (k+P_1)\right) \sqrt{\frac{2 C \zeta _n^2-k+P_1}{C}}}\right], \end{array}$

where $\mu_{n3} = \sqrt{-\frac{-2 C \zeta _n^2+k+P_1}{2C}}, \mu_{n5} = \sqrt{\frac{2 C \zeta _n^2-k+P_1}{2C}}, P_1 = \sqrt{k \left(\frac{4 C^2}{D}+k\right)}$ ,

$\begin{align*} \theta_{n,3} = &-\left[2 C^2 q \zeta _n \left(\cos \left(b \zeta _n\right)-1\right) \left(4 C^2-D (\nu -1) (k+P_1)\right) \left(4 C^2+D (\nu +1) (k-P_1)\right)\right]/\nonumber\\ &\left\{bD\left[32 C^5 (k-P_1) \zeta _n^2+16 C^4 k^2 (\nu -1)-8 C^3 D \zeta _n^2 \left(k^2 \left(\nu ^2-2\right)+2 k P_1+(1-2 \nu ) P_1^2\right)\right.\right.\nonumber\\ &\left.\left.+8 C^2 D k(\nu -1) \left(k^2-P_1^2\right)+D^2 (\nu -1) \left(k^2-P_1^2\right)^2-2 C D^2 \left(\nu ^2-1\right) (k-P_1)^2 (k+P_1) \zeta _n^2\right]\right\},\nonumber\\ \theta_{n,5} = &-\left[2 C^2 q \zeta _n \left(\cos \left(b \zeta _n\right)-1\right) \left(4 C^2-D (\nu -1) (k-P_1)\right) \left(4 C^2+D (\nu +1) (k+P_1)\right)\right]/\nonumber\\ &\left\{b D\left[32 C^5 (k+P_1) \zeta _n^2+16 C^4k^2 (\nu -1)-8 C^3 D \zeta _n^2 \left(k^2 \left(\nu ^2-2\right)-2 k P_1+(1-2 \nu ) P_1^2\right)\right.\right.\nonumber\\ &\left.\left.+8 C^2 D k (\nu -1) \left(k^2-P_1^2\right)+D^2 (\nu -1) \left(k^2-P_1^2\right)^2-2 C D^2 \left(\nu ^2-1\right) (k-P_1) (k+P_1)^2 \zeta _n^2\right]\right\}. \end{align*}$

Herein, $\Psi_{n, i}$ are determined by the simply supported conditions at $x = 0$ and $x = a$ . The constants $E = 2.06\times10^{11}pa, \ \nu = 0.3, \ q = 6.89\times10^{7} pa, \ a = 1.016m$ are taken from ^[33,34,35]. The corresponding results are shown in Table 1.

Table 1. Central deflection

$w$ of SEA and comparison with others.

$k$	$h/a$	$w$ (mm)
$(10^{5}N/m^{2})$		^[33]	^[34]	^[35]	Present
$54.25$	0.5	–	–	–	0.2758
	0.3	0.8153	0.8154	0.8588	0.8157
	0.2	2.2896	2.2896	2.3546	2.2858
	0.1	15.9992	15.987	16.1275	15.8836
	0.05	117.463	118.32	117.94	118.2312
$542.5$	0.5	–	–	–	0.2759
	0.3	0.8125	0.8129	0.8555	0.8162
	0.2	2.2670	2.2688	2.3307	2.2696
	0.1	14.9535	15.023	15.0655	15.0431
	0.05	77.6112	78.535	77.999	78.0131
$5425$	0.5	–	–	–	0.2763
	0.3	0.7846	0.7858	0.8247	0.8195
	0.2	2.0625	2.0615	2.1150	2.0843
	0.1	9.0424	9.0472	9.0833	9.0525
	0.05	17.668	17.921	17.703	17.8997

| Show Table

DownLoad: CSV

In Table 1, we only take the first 20 items of the series to achieve a better fitting effect with the data in the reference [33,,], which indicates that it is satisfactory to use the SEA to solve this moderately thick rectangular plate with single parameter foundation. presents the data convergence analysis for the central deflection $w$ when $\frac{h}{a} = 0.3$ with the different parameters $k$ . It can be seen that convergence of the series is satisfactory, and it confirms the completeness theorem of the eigenvector system discussed above in this paper.

Table 2. Convergence analysis of central deflection

$w$ when

$h/a = 0.3$ .

$k$				$w$ (mm)
$(10^{5}N/m^{2})$	$N$	1	5	8	10	15	20	25	30
$54.25$		0.8362	0.8169	0.8155	0.8162	0.8159	0.8157	0.8157	0.8157
$542.5$		0.8362	0.8173	0.8158	0.8165	0.8162	0.8162	0.8162	0.8162
$5425$		0.8397	0.8205	0.8191	0.8197	0.8195	0.8195	0.8195	0.8195

| Show Table

DownLoad: CSV

Example 2. The equilibrium equations of the bending problem of rectangular moderately thick plates with Pasternak type two-parameter elastic foundation are^[36,37]:

$\begin{align*} \left\{ \begin{array}{cc} \frac{\partial M_{x}}{\partial x}+\frac{\partial M_{xy}}{\partial y}-Q_{x} = 0,\\ \frac{\partial M_{y}}{\partial y}+\frac{\partial M_{xy}}{\partial x}-Q_{y} = 0,\\ \frac{\partial Q_{x}}{\partial x}+\frac{\partial Q_{y}}{\partial y}+G_f(\frac{\partial^2 w}{\partial x^2}+\frac{\partial^2 w}{\partial y^2})-k_fw = -q. \end{array}\right. \end{align*}$

For comparison with the data in ^[36,37], a square plate with the side length of $a$ is considered under fully simply supported conditions. This is the special case of our model when $q_{1}(x, y) = q_{2}(x, y) = 0$ , $q_{3}(x, y) = -q$ , $\alpha = \beta = C+G_f$ , $\gamma = -k_f$ . Furthermore, we take the general solution of the fifth and sixth components of ${\textbf{U}}(x, y)$ , and the specific forms of the deflection $w$ and the bending moment $M_x$ are as follows:

$\begin{align*} w = -\sum\limits_{n = 1}^{+\infty}&\left[\frac{e^{\mu_{n3} (-x)} \left(\left(e^{\mu_{n3} x}-1\right)^2 \theta _{n,3}+\mu_{n3} \Psi_{n3} e^{2 \mu_{n3} x}-\mu_{n3} \Psi_{n4}\right)}{\mu_{n3}}\right.\times\\ &\frac{C \sin \left(y \zeta _n\right) \left(4 C^2+C (\nu +3) G_f+D (\nu -1) \left(k_f-P_{2,n}\right)\right)}{\left(C+G_f\right) \left(4 C^2+C (\nu +1) G_f-D (\nu +1) \left(k_f+P_{2,n}\right)\right) \sqrt{\frac{C G_f+D k_f-D P_{2,n}}{2 C D+2 D G_f}+\zeta _n^2}}\\ &+\frac{e^{\mu_{n5} (-x)} \left(\left(e^{\mu_{n5} x}-1\right)^2 \theta _{n,5}+\mu_{n5} \Psi_{n5} e^{2 \mu_{n5} x}-\mu_{n5} \Psi_{n6}\right)}{\mu_{n5}}\times\\ &\left.\frac{C \sin \left(y \zeta _n\right) \left(4 C^2+C (\nu +3) G_f+D (\nu -1) \left(k_f+P_{2,n}\right)\right)}{\left(C+G_f\right) \left(4 C^2+C (\nu +1) G_f-D (\nu +1) \left(k_f-P_{2,n}\right)\right) \sqrt{\frac{C G_f+D k_f+D P_{2,n}}{2 C D+2 D G_f}+\zeta _n^2}}\right], \end{align*}$

$\begin{align*} M_x = -\sum\limits_{n = 1}^{+\infty}&\left\{\frac{\sin \left(y \zeta _n\right)(\Psi_{n1} D e^{\mu_{n1} x}-\Psi_{n2} D e^{-\mu_{n1}x})}{\mu_{n1} \left(2 C^2+D (\nu +1) \left(\left(C+G_f\right) \left(-\left(\zeta _n^2-\mu_{n1}^2\right)\right)-k_f\right)\right)}\times\left[2 C^2 \mu_{n1}^2 (\nu -1)+2C\nu\times\right.\right.\\ &\left.\left(\left(C+G_f\right)(\zeta _n^2-\mu_{n1}^2)-k_f\right)+D (\nu -1) \left(\mu_{n1}^2+\nu \zeta _n^2\right) \left(\left(C+G_f\right)\left(-\left(\zeta _n^2-\mu_{n1}^2\right)\right)-k_f\right)\right]\\ &+\frac{\sin \left(y \zeta _n\right)D e^{-\mu_{n3} x}\left(\left(e^{\mu_{n3} x}-1\right)^2 \theta_{n,3}+\mu_{n3}\left(\Psi_{n3} e^{2 \mu_{n3} x}-\Psi_{n4}\right)\right)}{\mu_{n3}^2 \left(2 C^2+D (\nu +1) \left(\left(C+G_f\right) \left(-\left(\zeta _n^2-\mu_{n3}^2\right)\right)-k_f\right)\right)}\times\left[2 C^2\mu_{n3}^2 (\nu -1)\right.\\ &\left.-2 C \nu \left(\left(C+G_f\right) \left(-\left(\zeta _n^2-\mu_{n3}^2\right)\right)-k_f\right)-D (\nu -1)\left(\mu_{n3}^2+\nu \zeta _n^2\right)\left(\left(C+G_f\right)\left(\zeta _n^2-\mu_{n3}^2\right)\right.\right.\\ &\left.\left.+k_f\right)\right]+\frac{\sin \left(y \zeta _n\right) D e^{-\mu_{n5}x} \left(\left(e^{\mu_{n5} x}-1\right)^2 \theta _{n,5}+\mu_{n5} \left(\Psi_{n5} e^{2 \mu_{n5} x}-\Psi_{n6}\right)\right)}{\mu_{n5}^2 \left(2 C^2+D (\nu +1) \left(\left(C+G_f\right) \left(-\left(\zeta _n^2-\mu_{n5}^2\right)\right)-k_f\right)\right)}\times\left[2 C^2\mu_{n5}^2\times \right.\\ &(\nu -1)-2 C \nu \left(\left(C+G_f\right) \left(-\left(\zeta _n^2-\mu_{n5}^2\right)\right)-k_f\right)-D (\nu -1) \left(\mu_{n5}^2+\nu \zeta _n^2\right)\left(\left(C+G_f\right)\times\right.\\ &\left.\left.\left(\left(\zeta _n^2-\mu_{n5}^2\right)+k_f\right)\right]\right\}, \end{align*}$

in which

$\begin{align*} &P_{2,n} = \sqrt{\frac{C^2 \left(G_f^2-4 D k_f\right)-2 C D G_f k_f+D^2 k_f^2}{D^2}},\\ &\mu_{n3} = \sqrt{\frac{C G_f+D k_f-D P_{2,n}}{2 C D+2 D G_f}+\zeta _n^2},\\ &\mu_{n5} = \sqrt{\frac{C G_f+D k_f+D P_{2,n}}{2 C D+2 D G_f}+\zeta _n^2}, \end{align*}$

$\theta_{n, 3}$ and $\theta_{n, 5}$ are displayed in Appendix B.

Taking the first 20 terms of the above formula, we calculate the center deflection $w$ and bending moment $M_x$ corresponding to different values of $h/a$ and $G_f$ . The relevant results are listed in and compared with the data in the references. The reference [] uses the grid method of $19 \times 19$ , and the reference [37] uses the analytic trigonometric series method. Compared with them, it is shown that the SEA in the presented paper is effective.

Table 3. The center deflection

$w$ and bending moment

$M_x$ of fully simply supported square plate under two-parameter foundation, where

$k_f = \frac{200D}{a^{4}}$ ,

$\nu = 0.25$ .

		$G_f=\frac{5D}{a^{2}}$		$G_f=\frac{20D}{a^{2}}$
$h/a$		$w \times(\frac{10^3 D }{q a^{4}})$	$M_x \times(\frac{10^2 }{q a^{2}})$	$w \times(\frac{10^3 D }{q a^{4}})$	$M_x \times(\frac{10^2 }{q a^{2}})$
$0.005$	Present	2.264014	2.418051	1.567611	1.612293
	^[36]	2.264013	2.417782	1.567607	1.612873
	^[37]	2.263888	2.417870	1.567556	1.612893
$0.100$	Present	2.313681	2.361804	1.587336	1.569424
	^[36]	2.313650	2.361704	1.587351	1.569177
	^[37]	2.303082	2.352117	1.581634	1.563309
$0.200$	Present	2.449645	2.207744	1.641966	1.449790
	^[36]	2.449532	2.207587	1.641869	1.450101
	^[37]	2.440886	2.198786	1.637296	1.444638

| Show Table

DownLoad: CSV

Example 3. For the free vibration of fully simply supported moderately thick rectangular plates, its governing equations are as follows^[38,39]:

$\begin{align*} \left\{ \begin{array}{cc} \frac{\partial M_{x}}{\partial x}+\frac{\partial M_{xy}}{\partial y}-Q_{x} = 0,\\ \frac{\partial M_{y}}{\partial y}+\frac{\partial M_{xy}}{\partial x}-Q_{y} = 0,\\ \frac{\partial Q_{x}}{\partial x}+\frac{\partial Q_{y}}{\partial y}+\rho \omega ^2 w = 0, \end{array}\right. \end{align*}$

where $\rho$ is the mass per unit area of the plate and $\omega$ is the natural frequency of the plates.

By analyzing the boundary conditions and letting the determinant of the coefficients be equal to zero, the frequency equation of the free vibration problem is given by

$\begin{align} &\frac{1}{4 \mu_{n5}^2}\left(C^2 D^2 \left(e^{2 a \mu_{n1}}-1\right) \left(e^{2 a \mu_{n3}}-1\right) \left(e^{2 a \mu_{n5}}-1\right) \left(\mu_{n1}^2-\mu_{n3}^2\right)^2 e^{-a (\mu_{n1}+\mu_{n3}+\mu_{n5})}\right.\\ &\left.\left(D (\nu -1) \zeta _n^2-2 C\right){}^2 \left(2 C \mu_{n5}^2+D (\nu -1) \zeta _n^2 \left(\zeta _n^2-\mu_{n5}^2\right)\right)\right) = 0. \end{align}$

(3.2)

To justify the frequency equation, we give the numerical results for the lowest six natural frequency parameters $\bar{k} = \frac{\rho \omega^2 b^4}{D \pi^4}$ in Table 4. Compared with ^[38] and ^[39], the satisfactory agreement further confirms the stability of the present model and method.

Table 4. The lowest six natural frequency parameters

$\bar{k}$ of fully simply supported moderately thick rectangular plates, where

$\delta_b$ =

$\frac{ \pi^2 D}{b^2 C}$ ,

$\vartheta$ =

$\frac{b}{a}$ and

$\nu$ = 0.3.

	$\vartheta$		0.2			0.6			1			2
$\delta_b$	$Modes$	^[38]	^[39]	$Present$	^[38]	^[39]	$Present$	^[38]	^[39]	$Present$	^[38]	^[39]	$Present$
$0.05$	1	1.028	1.028	1.0281	1.732	1.732	1.7318	3.636	3.636	3.6364	20.00	20.00	20.0000
	2	1.272	1.272	1.2718	5.306	5.306	5.3062	20.00	20.00	20.0000	45.71	45.71	45.7143
	3	1.732	1.732	1.7318	14.83	14.83	14.8330	45.71	45.71	45.7143	102.4	102.4	102.4240
	4	-	-	2.4858	-	-	15.6072	-	-	66.6667	-	-	156.2160
	5	-	-	3.6364	-	-	23.2654	-	-	102.4240	-	-	200.0000
	6	-	-	5.3062	-	-	34.1537	-	-	156.2160	-	-	277.7780
$0.10$	1	0.980	0.980	0.9797	1.628	1.628	1.6282	3.333	3.333	3.3333	16.67	16.67	16.6667
	2	1.206	1.206	1.2057	4.786	4.786	4.7858	16.67	16.67	16.6667	35.56	35.56	35.5556
	3	1.628	1.628	1.6282	12.63	12.63	12.6247	35.56	35.56	35.5556	73.48	73.48	73.4783
	4	-	-	2.3107	-	-	13.2379	-	-	50.0000	-	-	107.0370
	5	-	-	3.3333	-	-	19.1668	-	-	73.4783	-	-	133.3330
	6	-	-	4.7859	-	-	27.2657	-	-	107.0370	-	-	178.5710
$0.20$	1	0.895	0.895	0.89536	1.454	1.454	1.4541	2.857	2.857	2.8571	12.50	12.50	12.5000
	2	1.092	1.092	1.0922	4.001	4.001	4.0011	12.50	12.50	12.5000	24.62	24.62	24.6154
	3	1.454	1.454	1.4541	9.728	9.728	9.7281	24.62	24.62	24.6154	46.94	46.94	46.9444
	4	-	-	2.0253	-	-	10.1547	-	-	33.3333	-	-	65.6818
	5	-	-	2.8571	-	-	14.1732	-	-	46.9444	-	-	80.0000
	6	-	-	4.0011	-	-	19.4293	-	-	65.6818	-	-	104.1670
$0.40$	1	0.764	0.764	0.7638	1.198	1.198	1.1979	2.222	2.222	2.2222	8.333	8.333	8.3333
	2	0.919	0.919	0.9191	3.013	3.013	3.0130	8.333	8.333	8.3333	15.24	15.24	15.2381
	3	1.198	1.198	1.1979	6.668	6.668	6.6683	15.24	15.24	15.2381	27.26	27.26	27.2581
	4	-	-	1.6242	-	-	6.9277	-	-	20.0000	-	-	37.0513
	5	-	-	2.2222	-	-	9.3179	-	-	27.2581	-	-	44.4444
	6	-	-	3.0130	-	-	12.3374	-	-	37.0513	-	-	56.8182
$0.60$	1	0.666	0.666	0.6660	1.019	1.019	1.0185	1.818	1.818	1.8182	6.250	6.250	6.2500
	2	0.793	0.793	0.7934	2.416	2,416	2.4162	6.250	6.250	6.2500	11.03	11.03	11.0345
	3	1.019	1.019	1.0185	5.073	5.073	5.0727	11.03	11.03	11.0345	19.21	19.21	19.2045
	4	-	-	1.3557	-	-	5.2571	-	-	14.2857	-	-	25.8036
	5	-	-	1.8182	-	-	6.9403	-	-	19.2045	-	-	30.7692
	6	-	-	2.4162	-	-	9.0383	-	-	25.8036	-	-	39.0625
$0.80$	1	0.590	0.590	0.5904	0.886	0.886	0.8858	1.538	1.538	1.5385	5.000	5.000	5.0000
	2	0.698	0.698	0.6979	2.017	2.017	2.0168	5.000	5.000	5.0000	8.649	8.649	8.6487
	3	0.886	0.886	0.8858	4.093	4.093	4.0933	8.649	8.649	8.6487	14.83	14.83	14.8246
	4	-	-	1.1633	-	-	4.2357	-	-	11.1111	-	-	19.7945
	5	-	-	1.5385	-	-	5.5295	-	-	14.8246	-	-	23.5294
	6	-	-	2.0168	-	-	7.1313	-	-	19.7945	-	-	29.7619
$1.00$	1	0.530	0.530	0.5302	0.784	0.784	0.7837	1.333	1.333	1.3333	4.167	4.167	4.1667
	2	0.623	0.623	0.6230	1.731	1.731	1.7307	4.167	4.167	4.1667	7.111	7.111	7.1111
	3	0.784	0.784	0.7837	3.431	3.431	3.4308	7.111	7.111	7.1111	12.07	12.07	12.0714
	4	-	-	1.0188	-	-	3.5466	-	-	9.0909	-	-	16.0556
	5	-	-	1.3333	-	-	4.5953	-	-	12.0710	-	-	19.0476
	6	-	-	1.7307	-	-	5.8889	-	-	16.0560	-	-	24.0385

| Show Table

DownLoad: CSV

4. Conclusions

In this paper, the model for a class of moderately thick rectangular plates is considered, which includes the bending, free vibration and buckling problems of plates and is expected to be widely used in functionally graded material problems. Based on the equilibrium equations of the model, the relationship between force and displacement, the Hamiltonian system and the corresponding Hamiltonian operator are derived, which provides the possibility for the applications of SEA and the symplectic superposition method. Based on the calculation of eigenvalues and the analysis of the properties of symplectic eigenvectors, the general solution under the simply supported boundary conditions of opposite edges is given by using the proved completeness theorem. In addition, the general solution is applied to two kinds of elastic foundation problems of moderately thick rectangular plates, which have not been solved by the SEA before. At the same time, the higher order and more accurate frequency parameters are solved for the free vibration problem of moderately thick rectangular plates with fully simply supported. The satisfactory numerical results confirm the effectiveness and universality of our model. It should be noted that the buckling problem of the plates is solved similarly to the free vibration problem except that the buckling factor is included in the values of $\alpha_{1}$ and $\beta_{1}$ in the model.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

This work was supported by the Natural Science Foundation of Inner Mongolia (No.2021MS01004)

Conflict of interest

All authors declare no conflicts of interest in this paper.

Appendix

Appendix A

$\begin{align*} A_{34,n} = &-\left(\left(b\left((1-\nu) \left(C^4-2 \alpha C^3+C^2 \left(\alpha ^2+2 \gamma D+D \zeta _n^2 (\alpha (\nu -1)-2 \beta )\right)+C D\alpha (2 \gamma+ \right.\right.\right.\right.\\ &\left.\left.\zeta _n^2 (\alpha (-\nu )+\alpha -2 \beta )\right)-D^2 \left(-\gamma +(\beta -\alpha ) \zeta _n^2+P_n\right) \left(\gamma +\zeta _n^2 (-(\alpha \nu +\beta ))+P_{n}\right)\right)\\ &\left(C^4-2 \alpha C^3+C^2 \left(\alpha ^2+2 \gamma D+2 D \zeta _n^2 (\alpha (\nu -2)-\beta )\right)+2 \alpha C D \left(\gamma -\zeta _n^2 (\alpha \nu +\beta )\right)\right.\\ &\left.-D^2 \left(-\gamma +(\beta -\alpha ) \zeta _n^2+P_{n}\right) \left(\gamma +\zeta _n^2 (-(2 \alpha \nu +\alpha +\beta ))+P_{n}\right)\right)+\alpha D \zeta _n^2\left(C^2 (\nu- \right.\\ &3)\left.+D (\nu +1) \left(-\gamma +(\beta -\alpha ) \zeta _n^2+P_{n}\right)-\alpha C (\nu +1)\right)\left(3 C^4 (\nu -1)+\alpha C^3 (2-6 \nu )\right.\\ &+D^2 (\nu -1) \left(-\gamma +(\beta -\alpha ) \zeta _n^2+P_{n}\right) \left(-\gamma +\zeta _n^2 (2 \alpha \nu +\alpha +\beta )-P_{n}\right)-2 \alpha C D(P_{n}\\ &\left.-2 \gamma \nu +\zeta _n^2 (\alpha ((\nu -2) \nu -1)+2 \beta \nu )+\nu P_{n}\right)+C^2 \left(\alpha ^2 (3 \nu +1)+2 D (\nu -1)(2 \gamma +\right.\\ &\left.\left.\zeta _n^2 (\alpha (\nu -1)-2 \beta )+P_{n}\right)\right)+4 C^2 D \alpha \zeta _n^2 \left(C^2 (\nu -1)-\alpha C (\nu +3)+D (\nu -1)(\gamma + \right.\\ &\left.\left.\left.(\alpha -\beta ) \zeta _n^2+P_{n}\right))(C^2-\alpha C+D \left(\nu \left(\gamma +(\alpha -\beta ) \zeta _n^2\right)-P_{n}\right)\right)\right)\left/\left(\sqrt{2} \alpha ^2 D \zeta _n^2\left(C^2\right.\right.\right.\\ &\left.(\nu -3)-\alpha C (\nu +1)+D (\nu +1) \left(-\gamma +(\beta -\alpha ) \zeta _n^2+P_{n}\right)\right)^{2}(C (\alpha -C)-D(P_{n}+\\ &\left.\alpha +C (\nu +1)+D (\nu +1) \left(\gamma +(\alpha -\beta ) \zeta _n^2+P_{n}\right)\right)^{2}, \end{align*}$

$\begin{align*} A_{56,n} = &\left(b(\left(\nu -1)\left(-C^4+2 \alpha C^3+C^2 \left(-\alpha ^2-2 \gamma D+2 D \zeta _n^2 (\beta -\alpha (\nu -2))\right)\right.+2 \alpha C D\right.\right.\\ &\left.\left(\zeta _n^2 (\alpha \nu +\beta )-\gamma \right)+D^2 \left(\gamma +(\alpha -\beta ) \zeta _n^2+P_{n}\right) \left(-\gamma +\zeta _n^2 (2 \alpha \nu +\alpha +\beta )+P_{n}\right)\right)\\ &\left(-C^4+2 \alpha C^3-C^2 \left(\alpha ^2+2 \gamma D+D \zeta _n^2 (\alpha (\nu -1)-2 \beta )\right)+D^2 \left(\gamma +(\alpha -\beta ) \zeta _n^2\right.\right.\\ &\left.+P_{n})(-\gamma +\zeta _n^2 (\alpha \nu +\beta )+P_{n})+\alpha C D \left(\zeta _n^2 (\alpha (\nu -1)+2 \beta )-2 \gamma \right)\right)-4 C^2 \alpha \left(C^2\right.\\ &\left.(\nu -1)-\alpha C (\nu +3)-D (\nu -1) \left(-\gamma +(\beta -\alpha ) \zeta _n^2+P_{n}\right)\right)\left(C^2-\alpha C+D\right.\\ &\left.(\nu \left(\gamma +(\alpha -\beta ) \zeta _n^2\right)+P_{n})\right)+\alpha \left(-C^2 (\nu -3)+\alpha C (\nu +1)+D(\nu +1)(P_{n}+\gamma +\right.\\ &\left.\left.(\alpha -\beta ) \zeta _n^2\right)\right)\left(3 C^4 (\nu -1)+\alpha C^3 (2-6 \nu )-D^2 (\nu -1) \left(\gamma +(\alpha -\beta ) \zeta _n^2+P_{n}\right)\right.\\ &(-\gamma +\zeta _n^2 (2 \alpha \nu +\alpha +\beta )+P_{n})+2 C D \alpha \left(2 \gamma \nu +\zeta _n^2 (\alpha (1-(\nu -2) \nu )-2 \beta \nu )+P_{n}\right.\\ &\left.\left.\left.\nu +P_{n})+C^2 \left(\alpha ^2 (3 \nu +1)-2 D (\nu -1) \left(-2 \gamma +\zeta _n^2 (\alpha (-\nu )+\alpha +2 \beta )+P_{n}\right)\right)\right)\right)\right))\\ &\left/\left(\sqrt{\frac{-C^2+\alpha C+D \left(-\gamma +(\alpha +\beta ) \zeta _n^2+P_{n}\right)}{\alpha D}}\sqrt{2}D\zeta _n^2\alpha ^2\left(C (\nu +1)-C^2 (\nu -3)\right.\right.\right.\\ &\left.\alpha +C (\nu +1)+D (\nu +1) \left(\gamma +(\alpha -\beta ) \zeta _n^2+P_{n}\right)\right)^{2}. \end{align*}$

Appendix B

$\begin{align*} \theta_{n,3} = &-\left(C D(C+G_f)q(4 C^2+C G_f(\nu +1)-D (\nu +1) (k_f+P_{2,n}))\left(4 C^2+D \left(k_f-P_{2,n}\right)\right.\right.\\ &\left.\left.C (\nu +3) G_f+(\nu -1)\right)\zeta _n\sqrt{\frac{2 \left(C G_f+D \left(k_f-P_{2,n}\right)\right)}{D \left(C+G_f\right)}+4 \zeta _n^2} (\cos(b \zeta _n)-1)\right)/ \\ &\left(b\sqrt{\frac{C G_f+D k_f-D P_{2,n}}{2 C D+2 D G_f}+\zeta _n^2}\left(4 C^2 D\left(4 C^2+C (\nu +3) G_f+D (\nu -1) \left(k_f-P_{2,n}\right)\right)\right.\right.\\ &\left(C+G_f\right)\zeta _n^2 \left(C G_f+D \left(\nu k_f+P_{2,n}\right)\right)+(1-\nu)\left(-4 C^2 D k_f+C^2 G_f^2-2 C D G_f k_f+\right.\\ &D^2 (k_f^2-P_{2,n}^2)\left.+D \zeta _n^2 \left(C+G_f\right) \left(-4 C^2-C (\nu +1) G_f+D (\nu +1) \left(k_f+P_{2,n}\right)\right)\right)\left(C^2 G_f^2\right.\\ &+D^2 (k_f^2- P_{2,n}^2)-4 C^2 D k_f-2 C D G_f k_f+2 D (C+G_f)\left(-4 C^2-C (\nu +1) G_f+ D\times\right.\\ &\left.\left.(\nu +1)\left(k_f+P_{2,n}\right)\right)\zeta _n^2\right)+D \zeta _n^2 \left(C+G_f\right) \left(-4 C^2-C (\nu +1) G_f+D (\nu +1) \left(k_f+P_{2,n}\right)\right)\\ &\left(C^3 (2-6 \nu ) \left(C+G_f\right)+3 C^4 (\nu -1)+D^2(k_f+P_{2,n})(\nu -1)\left(2 \left(C+G_f\right)+k_f-P_{2,n}\times\right.\right.\\ &\left.(\nu +1) \zeta _n^2\right)-2 C D \left(C+G_f\right) \left(\left(\nu ^2-1\right) \zeta _n^2 \left(C+G_f\right)+2 \nu k_f+\nu P_{2,n}+P_{2,n}\right)+\left(\left(C+G_f\right){}^2\right.\\ &\left.\left.\left.\left.2 D (\nu -1) \left((\nu -3) \zeta _n^2 \left(C+G_f\right)-2 k_f+P_{2,n}\right)+(3 \nu +1)\right)C^2\right)\right)\right), \end{align*}$

$\begin{align*} \theta_{n,5} = &\left(C (C+G_f)q(4 C^2+C G_f(\nu +1)-D (\nu +1) (k_f+P_{2,n}))\left(4 C^2+D \left(k_f+P_{2,n}\right)\right.\right.\\ &\left.\left.C (\nu +3) G_f+(\nu -1)\right)\zeta _n\sqrt{\frac{2 \left(C G_f+D \left(k_f+P_{2,n}\right)\right)}{D \left(C+G_f\right)}+4 \zeta _n^2} (\cos(b \zeta _n)-1)\right)/ \\ &\left(b\zeta _n\sqrt{\frac{C G_f+D k_f+D P_{2,n}}{2 C D+2 D G_f}+\zeta _n^2}\left(-4 C^2 \left(C+G_f\right)(C G_f+D \nu k_f-D P_{2,n})\right.\right.\left(4 C^2+\right.\\ +1 &\left.C (\nu +3) G_f+D (\nu -1) \left(k_f+P_{2,n}\right)\right)+\frac{1}{D \zeta _n^2}(\nu -1)\left(4 C^2 D k_f-C^2 G_f^2+2 C D G_f k_f\right.\\ &+D^2P_{2,n}^2-D^2k_f^2\left.+D \zeta _n^2 \left(C+G_f\right) \left(4 C^2+C (\nu +1) G_f-D (\nu +1) \left(k_f-P_{2,n}\right)\right)\right)\times\\ &\left(-C^2 G_f^2+4 C^2 D k_f+2 C D G_f k_f-D^2 k_f^2+D^2 P_{2,n}^2+2 D (C+G_f)\left(C (\nu +1) G_f\right.\right.+\\ &4 C^2-\left.\left.D (\nu +1) \left(k_f-P_{2,n}\right)\right)\zeta _n^2\right)+\left(C+G_f\right) \left(4 C^2+C (\nu) G_f-D (\nu +1)\left(k_f-P_{2,n}\right)\right) \\ &\left(4 C^3 \left(G_f-2 D (\nu -1) \zeta _n^2\right)+D^2 (\nu -1) \left(k_f-P_{2,n}\right) \left(2 (\nu +1) G_f \zeta _n^2+k_f+P_{2,n}\right)+C^2 \right.\\ &\left(-2 D (\nu -1) (\nu +5) G_f \zeta _n^2+4 D \left(-2 \nu k_f+k_f+P_{2,n}\right)+(3 \nu +1) G_f^2\right)+2 C DG_f \times\\ &\left.\left.\left.\left(-2 \nu k_f+\nu P_{2,n}+P_{2,n}\right)-\left(\nu ^2-1\right) \zeta _n^2 \left(D \left(P_{2,n}-k_f\right)+G_f^2\right)\right)\right)\right). \end{align*}$

References

[1]	K. Arrow, Social coice and individual values, New Haven: Yale University Press, 2012. http://dx.doi.org/10.12987/9780300186987
[2]	A. Turing, On computable numbers, with an application to the entscheidungsproblem, Proc. Lond. Math. Soc., 42 (1937), 230–265. http://dx.doi.org/10.1112/plms/s2-42.1.230 doi: 10.1112/plms/s2-42.1.230
[3]	E. Witten, Dynamical breaking of supersymmetry, Nucl. Phys. B, 188 (1981), 513–554. http://dx.doi.org/10.1016/0550-3213(81)90006-7 doi: 10.1016/0550-3213(81)90006-7
[4]	B. Zitova, J. Flusser, Image registration methods: a survey, Image Vision Comput., 21 (2003), 977–1000. http://dx.doi.org/10.1016/S0262-8856(03)00137-9 doi: 10.1016/S0262-8856(03)00137-9
[5]	J. Nash, Equilibrium points in $n$ -person games, PNAS, 36 (1950), 48–49. http://dx.doi.org/10.1073/pnas.36.1.48 doi: 10.1073/pnas.36.1.48
[6]	F. Black, M. Scholes, The pricing of options and corporate liabilities, J. Polit. Econ., 81 (1973), 637–654. http://dx.doi.org/10.1086/260062 doi: 10.1086/260062
[7]	S. Wasserman, K. Faust, Social network analysis: methods and applications, Cambridge: Cambridge University Press, 1994. http://dx.doi.org/10.1017/CBO9780511815478
[8]	J. Hutchinson, Fractals and self-similarity, Indiana U. Math. J., 30 (1981), 713–747.
[9]	M. Barnsley, A. Vince, The chaos game on a general iterated function system, Ergod. Theor. Dyn. Syst., 31 (2011), 1073–1079. http://dx.doi.org/10.1017/S0143385710000428 doi: 10.1017/S0143385710000428
[10]	L. Brouwer, Über abbildung von mannigfaltigkeiten, Math. Ann., 71 (1911), 97–115. http://dx.doi.org/10.1007/BF01456931 doi: 10.1007/BF01456931
[11]	J. Schauder, Der fixpunktsatz in funktionalraümen, Stud. Math., 2 (1930), 171–180.
[12]	E. Picard, Mémoire sur la théorie des équations aux dérivés partielles et la méthode des approximations successives, J. Math. Pure. Appl., 6 (1890), 145–210.
[13]	S. Banach, Sur les opérations dans les ensembles abstraites et leurs applications, Fund. Math., 3 (1922), 133–181.
[14]	W. Mann, Mean value methods in iteration, Proc. Amer. Math. Soc., 4 (1953), 506–510. http://dx.doi.org/10.2307/2032162 doi: 10.2307/2032162
[15]	S. Ishikawa, Fixed points by a new iteration method, Proc. Amer. Math. Soc., 44 (1974), 147–150. http://dx.doi.org/10.1090/S0002-9939-1974-0336469-5 doi: 10.1090/S0002-9939-1974-0336469-5
[16]	M. Noor, New approximation schemes for general variational inequalities, J. Math. Anal. Appl., 251 (2000), 217–229. http://dx.doi.org/10.1006/jmaa.2000.7042 doi: 10.1006/jmaa.2000.7042
[17]	R. Agarwal, D. O'Regan, D. Sahu, Iterative construction of fixed points of nearly asymptotically nonexpansive mappings, J. Nonlinear Convex Anal., 8 (2007), 61–79.
[18]	M. Abbas, T. Nazir, A new faster iteration process applied to constrained minimization and feasibility problems, Mat. Vesnik, 66 (2014), 223–234.
[19]	B. Thakur, D. Thakur, M. Postolache, A new iteration scheme for approximating fixed points of nonexpansive mappings, Filomat, 30 (2016), 2711–2720. http://dx.doi.org/10.2298/FIL1610711T doi: 10.2298/FIL1610711T
[20]	K. Ullah, M. Arshad, Numerical reckoning fixed points for Suzuki's generalized nonexpansive mappings via new iteration process, Filomat, 32 (2018), 187–196. http://dx.doi.org/10.2298/FIL1801187U doi: 10.2298/FIL1801187U
[21]	J. Ali, F. Ali, A new iterative scheme to approximating fixed points and the solution of a delay differential equation, J. Nonlinear Convex Anal., 21 (2020), 2151–2163.
[22]	M. Abbas, M. Asghar, M. De la Sen, Approximation of the solution of delay fractional differential equation using AA-iterative scheme, Mathematics, 10 (2022), 273. http://dx.doi.org/10.3390/math10020273 doi: 10.3390/math10020273
[23]	I. Beg, M. Abbas, M. Asghar, Convergence of $AA$ -iterative algorithm for generalized $\alpha$ -nonexpansive mappings with an application, Mathematics, 10 (2022), 4375. http://dx.doi.org/10.3390/math10224375 doi: 10.3390/math10224375
[24]	M. Asghar, M. Abbas, C. Eyni, M. Omaba, Iterative approximation of fixed points of generalized $\alpha_m$ -nonexpansive mappings in modular spaces, AIMS Mathematics, 8 (2023), 26922–26944. http://dx.doi.org/10.3934/math.20231378 doi: 10.3934/math.20231378
[25]	C. Suanoom, A. Gebrie, T. Grace, The convergence of $AA$ -iterative algorithm for generalized $AK$ - $\alpha$ -nonexpansive mappings in Banach spaces, Science and Technology Asia, 28 (2023), 82–90. http://dx.doi.org/10.14456/scitechasia.2023.47 doi: 10.14456/scitechasia.2023.47
[26]	M. Asghar, M. Abbas, B. Rouhani, The $AA$ -viscosity algorithm for fixed point, generalized equilibrium and variational inclusion problems, Axioms, 13 (2024), 38. http://dx.doi.org/10.3390/axioms13010038 doi: 10.3390/axioms13010038
[27]	M. Abbas, C. Ciobanescu, M. Asghar, A. Omame, Solution approximation of fractional boundary value problems and convergence analysis using $AA$ -iterative scheme, AIMS Mathematics, 9 (2024), 13129–13158. http://dx.doi.org/10.3934/math.2024641 doi: 10.3934/math.2024641
[28]	U. Kohlenbach, Some logical metatherems with applications in functional analysis, Trans. Amer. Math. Soc., 357 (2005), 89–128. http://dx.doi.org/10.2307/3845213 doi: 10.2307/3845213
[29]	T. Suzuki, Fixed point theorems and convergence theorems for some generalized nonexpansive mappings, J. Math. Anal. Appl., 340 (2008), 1088–1095. http://dx.doi.org/10.1016/j.jmaa.2007.09.023 doi: 10.1016/j.jmaa.2007.09.023
[30]	R. Pant, R. Shukla, Approximating fixed points of generalized $\alpha$ -nonexpansive mappings in Banach spaces, Numer. Func. Anal. Opt., 38 (2017), 248–266. http://dx.doi.org/10.1080/01630563.2016.1276075 doi: 10.1080/01630563.2016.1276075
[31]	R. Pant, R. Pandey, Existence and convergence results for a class of nonexpansive type mappings in hyperbolic spaces, Appl. Gen. Topol., 20 (2019), 281–295. http://dx.doi.org/10.4995/agt.2019.11057 doi: 10.4995/agt.2019.11057
[32]	A. Şahin, Some new results of $M$ -iteration process in hyperbolic spaces, Carpathian J. Math., 35 (2019), 221–232. http://dx.doi.org/10.37193/CJM.2019.02.10 doi: 10.37193/CJM.2019.02.10
[33]	A. Şahin, Some results of the Picard-Krasnoselskii hybrid iterative process, Filomat, 33 (2019), 359–365. http://dx.doi.org/10.2298/FIL1902359S doi: 10.2298/FIL1902359S
[34]	S. Khatoon, I. Uddin, M. Başarır, A modified proximal point algorithm for a nearly asymptotically quasi-nonexpansive mapping with an application, Comp. Appl. Math., 40 (2021), 250. http://dx.doi.org/10.1007/s40314-021-01646-9 doi: 10.1007/s40314-021-01646-9
[35]	H. Piri, B. Daraby, S. Rahrovi, M. Ghasemi, Approximating fixed points of generalized $\alpha$ -nonexpansive mappings in Banach spaces by new faster iteration process, Numer. Algor., 81 (2019), 1129–1148. http://dx.doi.org/10.1007/s11075-018-0588-x doi: 10.1007/s11075-018-0588-x
[36]	K. Ullah, J. Ahmad, M. De la Sen, On generalized nonexpansive maps in Banach spaces, Computation, 8 (2020), 61. http://dx.doi.org/10.3390/computation8030061 doi: 10.3390/computation8030061
[37]	S. Reich, I. Shafrir, Nonexpansive iterations in hyperbolic spaces, Nonlinear Anal.-Theor., 15 (1990), 537–558. http://dx.doi.org/10.1016/0362-546X(90)90058-O doi: 10.1016/0362-546X(90)90058-O
[38]	W. Takahashi, A convexity in metric space and nonexpansive mappings, Kodai Math. Sem. Rep., 22 (1970), 142–149. http://dx.doi.org/10.2996/kmj/1138846111 doi: 10.2996/kmj/1138846111
[39]	M. Bridson, A. Haefliger, Metric spaces of non-positive curvature, Berlin: Springer, 2013. http://dx.doi.org/10.1007/978-3-662-12494-9
[40]	A. Khan, H. Fukhar-ud-din, M. Khan, An implicit algorithm for two finite families of nonexpansive maps in hyperbolic spaces, Fixed Point Theory Appl., 2012 (2012), 54. http://dx.doi.org/10.1186/1687-1812-2012-54 doi: 10.1186/1687-1812-2012-54
[41]	L. Leuştean, Nonexpansive iterations in uniformly convex W-hyperbolic spaces, In: Nonlinear analysis and optimization I: nonlinear analysis, Ramat-Gan: American Mathematical Society, 2010,193–209. http://dx.doi.org/10.1090/conm/513/10084
[42]	T. Lim, Remarks on some fixed point theorems, Proc. Amer. Math. Soc., 60 (1976), 179–182. http://dx.doi.org/10.1090/S0002-9939-1976-0423139-X doi: 10.1090/S0002-9939-1976-0423139-X
[43]	T. Cardinali, P. Rubbioni, A generalization of the Caristi fixed point theorem in metric spaces, Fixed Point Theory, 11 (2010), 3–10.
[44]	I. Timiş, On the weak stability of Picard iteration for some contractive type mappings, Ann. Univ. Craiova-Mat., 37 (2010), 106–114.
[45]	S. Hilger, Analysis on measure chains a unified approach to continuous and discrete calculus, Results Math., 18 (1990), 18–56. http://dx.doi.org/10.1007/BF03323153 doi: 10.1007/BF03323153
[46]	M. Bohner, A. Peterson, Dynamic equations on time scales: an introduction with applications, Boston: Birkhauser, 2001. http://dx.doi.org/10.1007/978-1-4612-0201-1
[47]	G. Guseinov, Integration on time scales, J. Math. Anal. Appl., 285 (2003), 107–127. http://dx.doi.org/10.1016/S0022-247X(03)00361-5 doi: 10.1016/S0022-247X(03)00361-5
[48]	G. Adomian, A review of the decomposition method in applied mathematics, J. Math. Anal. Appl., 135 (1988), 501–544. http://dx.doi.org/10.1016/0022-247X(88)90170-9 doi: 10.1016/0022-247X(88)90170-9
[49]	A. Şahin, E. Öztürk, G. Aggarwal, Some fixed-point results for the $KF$ -iteration process in hyperbolic metric spaces, Symmetry, 15 (2023), 1360. http://dx.doi.org/10.3390/sym15071360 doi: 10.3390/sym15071360
[50]	H. Senter, W. Dotson, Approximating fixed points of nonexpansive mappings, Proc. Amer. Math. Soc., 44 (1974), 375–380.
[51]	R. Farengo, Y. Lee, P. Guzdar, An electromagnetic integral equation: application to microtearing modes, Phys. Fluids, 26 (1983), 3515–3523. http://dx.doi.org/10.1063/1.864112 doi: 10.1063/1.864112
[52]	A. Manzhirov, On a method of solving two-dimensional integral equations of axisymmetric contact problems for bodies with complex rheology, J. Appl. Math. Mech., 49 (1985), 777–782. http://dx.doi.org/10.1016/0021-8928(85)90016-4 doi: 10.1016/0021-8928(85)90016-4
[53]	M. Mirkin, A. Bard, Multidimensional integral equations: part 1. a new approach to solving microelectrode diffusion problems, J. Electroanal. Chem., 323 (1992), 1–27. http://dx.doi.org/10.1016/0022-0728(92)80001-K doi: 10.1016/0022-0728(92)80001-K
[54]	Z. Kalkan, A. Şahin, A. Aloqaily, N. Mlaiki, Some fixed point and stability results in $b$ -metric-like spaces with an application to integral equations on time scales, AIMS Mathematics, 9 (2024), 11335–11351. http://dx.doi.org/10.3934/math.2024556 doi: 10.3934/math.2024556

This article has been cited by:

1.	Ke Yuan, Chang Gao, Yujin Lin, Haidong Yu, A new analytical model for predicting the welding deformation of large welded plates with equivalent moments, 2024, 1362-1718, 10.1177/13621718241288693
2.	Mengmeng Zhang, Eburilitu Bai, Jinglong Wang, Symplectic superposition solution for the buckling problem of orthotropic rectangular plates with four clamped edges, 2025, 95, 0939-1533, 10.1007/s00419-024-02724-0

Reader Comments

Your name:*

Email:*
© 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Mathematics

1.8 3.4

Metrics

Article views(864) PDF downloads(72) Cited by(0)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

AIMS Mathematics

The AA AA -iterative algorithm in hyperbolic spaces with applications to integral equations on time scales

Related Papers:

Abstract

1. Introduction

2. The model of the moderately thick plates and SEA

2.1. Basic equations and Hamiltonian system

2.2. Eigenvalues and eigenvectors of H {\textbf{H}}

2.3. Properties of eigenvectors

2.3.1. When the eigenvalues are simple

2.3.2. When the eigenvalues are double

2.4. General solution

2.4.1. When eigenvalues are simple

2.4.2. When eigenvalues are double

3. Applications

4. Conclusions

Use of AI tools declaration

Acknowledgments

Conflict of interest

Appendix

Appendix A

Appendix B

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog

Abstract

1. Introduction

2. The model of the moderately thick plates and SEA

2.1. Basic equations and Hamiltonian system

2.2. Eigenvalues and eigenvectors of H {\textbf{H}}

2.3. Properties of eigenvectors

2.3.1. When the eigenvalues are simple

2.3.2. When the eigenvalues are double

2.4. General solution

2.4.1. When eigenvalues are simple

2.4.2. When eigenvalues are double

3. Applications

4. Conclusions

Use of AI tools declaration

Acknowledgments

Conflict of interest

Appendix

Appendix A

Appendix B

References

The $AA$ -iterative algorithm in hyperbolic spaces with applications to integral equations on time scales

2.2. Eigenvalues and eigenvectors of ${\textbf{H}}$

2.2. Eigenvalues and eigenvectors of ${\textbf{H}}$