A computing method for bending problem of thin plate on Pasternak foundation

Abstract

The governing equation of the bending problem of simply supported thin plate on Pasternak foundation is degraded into two coupled lower order differential equations using the intermediate variable, which are a Helmholtz equation and a Laplace equation. A new solution of two-dimensional Helmholtz operator is proposed as shown in Appendix 1. The R-function and basic solutions of two-dimensional Helmholtz operator and Laplace operator are used to construct the corresponding quasi-Green function. The quasi-Green’s functions satisfy the homogeneous boundary conditions of the problem. The Helmholtz equation and Laplace equation are transformed into integral equations applying corresponding Green’s formula, the fundamental solution of the operator, and the boundary condition. A new boundary normalization equation is constructed to ensure the continuity of the integral kernels. The integral equations are discretized into the nonhomogeneous linear algebraic equations to proceed with numerical computing. Some numerical examples are given to verify the validity of the proposed method in calculating the problem with simple boundary conditions and polygonal boundary conditions. The required results are obtained through MATLAB programming. The convergence of the method is discussed. The comparison with the analytic solution shows a good agreement, and it demonstrates the feasibility and efficiency of the method in this article.

Keywords

Plate bending problem Pasternak foundation Helmholtz equation Laplace equation R-function Green’s function

Introduction

In the bending problem of Pasternak foundation plate, its governing differential equation is a higher order equation and contains biharmonic operator.¹ Therefore, it is difficult to study the plate with complex boundary conditions in depth. R-function theory and Green’s function method can always get the corresponding boundary normalization equation for the thin plate with complex boundary conditions and transform the high-order control differential equation into the low-order equations which are a Helmholtz equation² and a Laplace equation³ and then get the effective corresponding solution. Scholars at home and abroad have conducted extensive research on the plate problem on Pasternak foundation, Helmholtz equation, and Laplace equation.

The solution of the bending problem of rectangular plate with free edge on Pasternak foundation was obtained by Cai⁴ using the Fourier series. The vibration characteristics of elastic rotation constrained laminated plates on Pasternak foundation was studied by Yang et al.⁵ using differential quadrature (DQ)–Galerkin semi-analytical method. The analytical solution of the bending problem of rectangular thin plates with four sides simply supported on Pasternak foundation was studied by Wang.⁶ The bending of rectangular plate with free edge on Pasternak foundation was studied by Liu and Peng⁷ using DQ method of control equation. The axisymmetric large deflection of circular plates on Pasternak foundation was studied by Liu⁸ using the perturbation method combined with the integral equation method. Vibration analysis of a thin functionally graded plate having an out-of-plane material inhomogeneity resting on Winkler–Pasternak foundation under different combinations of boundary conditions was researched by Piyush et al.⁹ The exact nonstationary responses of rectangular thin plate on Pasternak foundation excited by stochastic moving loads was researched by Chen et al.¹⁰ A Frobenius exact solution method for stability analysis of orthotropic rectangular thin plate under biaxial nonlinear in-plane loading resting on Pasternak foundation was researched by Hosein et al.¹¹ A semi-analytical closed-form solution on static behavior of thin skew plates on Winkler and Pasternak foundations was researched by Amin Joodaky and Iman Joodaky.¹² Dual reciprocity hybrid radial boundary node method for Winkler and Pasternak foundation thin plate was researched by Yan et al.¹³ Vibration analysis of thin plates resting on Pasternak foundations was researched by Ehsan et al.¹⁴ using element free Galerkin method. Three approximate analytical solutions for thermal buckling of clamped thin rectangular functionally graded material plates resting on Pasternak elastic foundation were researched by Kiani et al.¹⁵ The performance of a rotated 5-point Laplacian operator for computing the harmonic potentials was investigated by Dahalan et al.¹⁶ Koyunbakan et al.¹⁷ analyzed a nonlinear eigenvalue problem for the p-Laplacian operator with zero Dirichlet boundary conditions. A metaheuristic optimization method for discretization of fractional Laplacian without discretization operator was studied by Mahata et al.¹⁸ Tan and Li¹⁹ studied the solutions for nonlinear fractional differential equations with p-Laplacian operator nonlocal boundary value problem in a Banach space. The existence and uniqueness of solutions for fractional boundary value problems with p-Laplacian operator were studied by Bai.²⁰ The edge extraction of quantum image based on Laplace operator and zero-crossing method was studied by Fan et al.²¹ The several positive solutions of nonlinear mixed fractional differential equation with p-Laplacian operator were studied by Li.²² The stereo-matching algorithm of weighted guided image filter based on Gauss Laplace operator was studied by Bo et al.²³ The damage detection of building after disaster based on local Laplace operator was studied by Li et al.²⁴ Assari et al.² provided a numerical method for solving logarithmic Fredholm integral equations which occur as a reformulation of two-dimensional Helmholtz equations over the unit circle with the Robin boundary conditions. Kravcenko et al.³ proposed a distributed fast boundary element methods for Helmholtz problems. Tian et al.²⁵ devoted their work to study the error analysis and stability of the method of fundamental solutions for the case of the modified Helmholtz equation. A wavelet collocation method was proposed by Chen et al.²⁶ for solving the linear boundary integral equations reformulated from the modified Helmholtz equation with Robin boundary conditions. Wang et al.²⁷ made a first attempt to use a new localized method of fundamental solutions to accurately and stably solve the inverse Cauchy problems of two-dimensional Laplace and biharmonic equations in complex geometries.

Li and Yuan^28–31 studied the thin plate and shell problems using the mathematical method of Laplace equation. In this article, R-function theory and quasi-Green function methods of Helmholtz equation as shown in Appendix 1 and Laplace equation are used to study the bending of thin plates on Pasternak foundation. First, the biharmonic operator of the plate bending problem on Pasternak foundation is reduced to Helmholtz equation and Laplace equation. Using Green’s function, the Helmholtz equation and Laplace equation are transformed into two integral equations, and a new boundary normalization equation is constructed to ensure the continuity of the integral kernels. Then the two integral equations can be discretized and programmed by MATLAB, and the bending result of thin plate can be obtained. The R-function theory guarantees that the corresponding boundary normalization equation can be found for any complex region. Thus, the original problem can be transformed into the second kind of Fredholm integral equation without singularity. Five numerical examples are given to verify the validity of the proposed method as shown in Appendix 1 which uses the quasi-Green function method of Helmholtz equation combing with the quasi-Green function method of Laplace equation.

Basic equation

According to the Kirchhoff hypothesis, under transverse load $Z (x)$ , the governing equation and boundary conditions for the bending problem of simply-supported thin plate on Pasternak foundation³² are

D \nabla^{4} w - G_{p} \nabla^{2} w + kw (x) = Z (x), x \in Ω

(1a)

w = 0, M_{n} = 0, x \in Γ

(1b)

in which $\nabla^{4} = (\partial^{2} / \partial x_{1}^{2} + \partial^{2} / \partial x_{2}^{2})^{2}$ is biharmonic operator, the bending stiffness is $D = E h^{3} / 12 (1 - ν^{2})$ , $\nabla^{2} = \partial^{2} / \partial x_{1} + \partial^{2} / \partial x_{2}$ is Laplace operator, $h$ is the thickness of the plate, $E$ and $ν$ are elastic modulus and Poisson’s ratio, respectively, $k$ is elastic foundation coefficient, $G_{p}$ is shear modulus of foundation, and $w$ is deflection. $x = (x_{1}, x_{2})$ , $Ω$ is the domain contained in the middle plane ( $x_{1} - x_{2}$ plane) of the plate in Cartesian coordinate system.

Introducing intermediate variable

M = \frac{1}{1 + ν} (M_{1} + M_{2})

(2)

where $M_{1}$ and $M_{2}$ are expressed as

\begin{matrix} M_{1} = - D (\frac{\partial^{2} w}{\partial {x_{1}}^{2}} + ν \frac{\partial^{2} w}{\partial {x_{2}}^{2}}) \\ M_{2} = - D (\frac{\partial^{2} w}{\partial {x_{2}}^{2}} + ν \frac{\partial^{2} w}{\partial {x_{1}}^{2}}) \end{matrix}

(3)

By substituting equation (2) into differential equation (1), the differential equation for bending problem of thin plates on Pasternak foundation can be degraded into Helmholtz equation and Laplace equation

\nabla^{2} M - \frac{G_{p}}{D} M = kw - Z

(4)

\nabla^{2} w = - \frac{M}{D}

(5)

where $\nabla^{2} = \partial^{2} / \partial x_{1}^{2} + \partial^{2} / \partial x_{2}^{2}$ is a Laplace operator and $x \in Ω$ .

Boundary constraints

The deflection on the plate boundary should be zero and the bending moment should be equal to zero under simply supported constraint boundary conditions

w = 0, x \in Γ

(6)

M = 0, x \in Γ

(7)

The integral equations

Using the solution of two-dimensional Helmholtz operator, the boundary value problems (equations (4) and (7)) are integrated as (according to Appendix 1)

M (x) = M_{0} (x) + \int_{Ω} M (ξ) K_{M} (x, ξ) d_{ξ} Ω

(8)

M_{0} (x) = \frac{- k}{2 π} \int_{Ω} G_{M} (x, ξ) w (ξ) d_{ξ} Ω + \frac{Z}{2 π} \int_{Ω} G_{M} (x, ξ) d_{ξ} Ω

(9)

K_{M} (x, ξ) = - \frac{1}{2 π} (\nabla_{ξ}^{2} - λ) e_{M} (x, ξ)

(10)

Using the quasi-Green’s function method of Laplace operator in references,^28–31 the boundary value problems (equations (5) and (6)) are integrated as

w (x) = w_{0} (x) + \int_{Ω} w (ξ) K_{w} (x, ξ) d_{ξ} Ω

(11)

w_{0} (x) = - \frac{1}{D} \int_{Ω} G_{w} (x, ξ) M (ξ) d_{ξ} Ω

(12)

K_{w} (x, ξ) = (\frac{\partial^{2}}{\partial ξ_{1}^{2}} + \frac{\partial^{2}}{\partial ξ_{2}^{2}}) e_{w} (x, ξ)

(13)

According to the corresponding formulas in references^28–31 and Appendix 1, the following is obtained

\begin{matrix} K_{M} (x, ξ) = - \frac{1}{2 π} \sqrt{λ} k_{1} (\sqrt{λ} R) \\ {- \frac{2}{R} [1 + ω (x) \nabla^{2} ω] + \frac{2}{R^{3}} {[\vec{r} + 2 ω (x) \nabla ω]}^{2}} \\ - \frac{1}{2 π} λ k_{0} (\sqrt{λ} R) {\frac{1}{R^{2}} {[\vec{r} + 2 ω (x) \nabla ω]}^{2} - 1} \end{matrix}

(14)

\begin{matrix} K_{w} (x, ξ) = \\ \frac{R^{2} ω (x) \nabla^{2} ω + 4 ω (x) ω - 4 (\vec{r} \cdot \nabla ω) ω (x) - 4 ω^{2} (x) {(\nabla ω)}^{2}}{π R^{4}} \end{matrix}

(15)

In the formulas, $k_{1}$ is the second kind of first-order deformed Bessel function, $ω = ω (ξ)$ , $\nabla = \nabla_{ξ}$ . The singularity of integral kernels is treated by references,^28–31 which guarantees the continuity of integral kernels $K_{M} (x, ξ)$ and $K_{w} (x, ξ)$ .

The integral equations (8) and (11) are discretized by the method of references.^28–31 The natural frequencies and corresponding modes can be obtained by programming of MATLAB.

Numerical examples

Example 1

Simply supported thin plate on Pasternak foundation model is shown in Figure 1. We set a = b = 1, Poisson’s ratio v = 0.3, thickness t = 0.1, Elastic modulus E = 4 × 10⁶, the coefficient of elastic foundation k = 4 × 10⁴, k = 5 × 10⁴, and k = 6 × 10⁴, respectively, Shear modulus of foundation G_P = 1 × 10⁴, and uniform load subjected Z = 1000. According to R-function theory,^33–35 the following is obtained

ω_{0} = ω_{1} + ω_{2} - \sqrt{ω_{1}^{2} + ω_{2}^{2}}

in which $ω_{1} = (a^{2} - x^{2}) / (2 a) \geq 0$ and $ω_{2} = (b^{2} - y^{2}) / (2 b) \geq 0$ . Then $ω_{0} = 0$ is the first-order boundary normalized equation of the square plate; $ω_{1}$ and $ω_{2}$ are the sides of the square plate.

Figure 1.

Simply supported plates on Pasternak foundation.

25 (5 × 5), 49 (7 × 7), 81 (9 × 9), 121 (11 × 11), and 169 (13 × 13) grid layout schemes are adopted for square plates, respectively. The results are shown in Table 1 and Figure 2.

Table 1.

Maximum deflection w of square plate.

Method		The coefficient of elastic foundation
		k = 4 × 10⁴		k = 5 × 10⁴		k = 6 × 10⁴
		Maximum deflection	Error (%)	Maximum deflection	Error (%)	Maximum deflection	Error (%)
The method in this paper	25	0.01283	12.90	0.01174	11.329	0.01081	9.992
	49	0.01386	5.906	0.01257	5.060	0.01149	4.330
	81	0.01424	3.327	0.01287	2.795	0.01172	2.415
	121	0.01441	2.172	0.01300	1.813	0.01183	1.499
	169	0.01451	1.494	0.01308	1.208	0.01188	1.082
Analytic solutions		0.01473	—	0.01324	—	0.01201	—

Figure 2.

Deflection of square plate x = 0.

Deflection diagrams of square plate deflection of x = 0 with different elastic foundation coefficients k using 169 (13 × 13) grids are shown in Figure 2.

Example 2

Simply supported thin plate on Pasternak foundation model is shown in Figure 1. We set a = 1.25, b = 1; other parameters applied are shown in example 1. According to R-function theory,^33–35 the following is obtained

ω_{0} = ω_{1} + ω_{2} - \sqrt{ω_{1}^{2} + ω_{2}^{2}}

in which $ω_{1} = (a^{2} - x^{2}) / (2 a) \geq 0$ and $ω_{2} = (b^{2} - y^{2}) / (2 b) \geq 0$ . Then $ω_{0} = 0$ is the first-order boundary normalized equation of the rectangular plate; $ω_{1}$ and $ω_{2}$ are the sides of the rectangular plate.

The results are shown in Table 2 and Figure 3.

Table 2.

Maximum deflection w of rectangular plate.

Method		The coefficient of elastic foundation
		k = 4 × 10⁴		k = 5 × 10⁴		k = 6 × 10⁴
		Maximum deflection	Error (%)	Maximum deflection	Error (%)	Maximum deflection	Error (%)
The method in this article	25	0.01416	13.606	0.01281	11.777	0.01168	10.223
	49	0.01538	6.162	0.01377	5.165	0.01244	4.381
	81	0.01582	3.478	0.01410	2.893	0.01270	2.383
	121	0.01603	2.196	0.01425	1.860	0.01281	1.537
	169	0.01614	1.525	0.01434	1.240	0.01287	1.076
Analytic solutions		0.01639	—	0.01452	—	0.01301	—

Figure 3.

Rectangular deflection of x = 0.

Deflection diagrams of rectangular plate deflection of x = 0 with different elastic foundation coefficients k using 169 (13 × 13) grids are shown in Figure 3.

Example 3

Simply supported trapezoidal thin plates Pasternak foundation model is shown in Figure 4. We set a = 1.5, b = 1, c = 2. Other parameters applied are shown in Example 1. According to R-function theory,^33–35 the following is obtained

\begin{matrix} ω_{0} = ω_{1} + ω_{2} + ω_{3} - \sqrt{ω_{1}^{2} + ω_{2}^{2}} - \sqrt{ω_{1}^{2} + ω_{3}^{2}} \\ - \sqrt{ω_{2}^{2} + ω_{3}^{2}} + \sqrt{ω_{1}^{2} + ω_{2}^{2} + ω_{3}^{2}} \end{matrix}

in which

$ω_{2} = \frac{1}{\sqrt{1 + {(\frac{c}{a - b})}^{2}}} (\frac{ac}{a - b} + \frac{c}{a - b} x - y)$

$ω_{3} = \frac{1}{\sqrt{1 + {(\frac{c}{a - b})}^{2}}} (\frac{ac}{a - b} - \frac{c}{a - b} x - y)$

and

$ω_{1} = \frac{(c - y) y}{c}$

Then $ω_{0} = 0$ is the first-order boundary normalized equation of trapezoidal plates. $ω_{1} = 0$ , $ω_{2} = 0$ , and $ω_{3} = 0$ are each side of the trapezoidal plate.

Figure 4.

Simply supported trapezoidal plates on Pasternak foundation.

Deflection diagrams of trapezoidal plate deflection of x = 0 with different elastic foundation coefficients k using 169 (13 × 13) grids are shown in Figure 5.

Figure 5.

Trapezoidal plate deflection of x = 0.

Example 4

Simply supported L-shaped thin plates on Pasternak foundation model are shown in Figure 6. We set a = 1.5, b = 2.0, c = 0.75, d = 1.0, uniform load subjected Z = 800. Other parameters applied are shown in Example 1. According to R-function theory,^33–35 the following is obtained

ω_{0} = ω_{5} + ω_{6} - \sqrt{ω_{5}^{2} + ω_{6}^{2}}

in which $ω_{5} = ω_{1} + ω_{2} - \sqrt{ω_{1}^{2} + ω_{2}^{2}}$ and $ω_{6} = ω_{3} + ω_{4} + \sqrt{ω_{3}^{2} + ω_{4}^{2}}$ , and $ω_{1} = (a^{2} - x^{2}) / (2 a) \geq 0$ , $ω_{2} = (b^{2} - y^{2}) / (2 b) \geq 0$ , $ω_{3} = (c - x) \geq 0$ , $ω_{4} = (d - y) \geq 0$ . Then $ω_{0} = 0$ is the first-order boundary normalized equation for L-shaped plates; $ω_{1} = 0$ , $ω_{2} = 0$ , $ω_{3} = 0$ , and $ω_{4} = 0$ are the sides of the L-shaped plate.

Figure 6.

Simply supported L-shaped plates on Pasternak foundation.

Deflection diagrams of L-shaped plate deflection of x = −0.0536 with different elastic foundation coefficients k using 147 (49 × 3) grids are shown in Figure 7.

Figure 7.

L-shaped plate deflection of x = −0.0536.

Example 5

A simply supported I-shaped thin plate on Pasternak foundation model is shown in Figure 8. We set a = 1.2, b = 1.5, c = 1, d = 1, uniform load subjected Z = 800. Other parameters applied are shown in example 1. According to R-function theory,^33–35 the following is obtained

ω_{0} = ω_{5} + ω_{6} - \sqrt{ω_{5}^{2} + ω_{6}^{2}}

in which $ω_{5} = ω_{1} + ω_{2} - \sqrt{ω_{1}^{2} + ω_{2}^{2}}$ and $ω_{6} = ω_{3} + ω_{4} + \sqrt{ω_{3}^{2} + ω_{4}^{2}}$ , and $ω_{1} = (a^{2} - x^{2}) / (2 a) \geq 0$ , $ω_{2} = (b^{2} - y^{2}) / (2 b) \geq 0$ , $ω_{3} = (c^{2} - x^{2}) / (2 c) \geq 0$ , $ω_{4} = (y^{2} - d^{2}) / (2 d) \geq 0$ . Then $ω_{0} = 0$ is the first-order boundary normalized equation of I-shaped plates. $ω_{1} = 0$ , $ω_{2} = 0$ , $ω_{3} = 0$ , and $ω_{4} = 0$ are each side of the I-shaped plate.

Figure 8.

Simply supported I-shaped plates on Pasternak foundation.

Deflection diagrams of I-shaped plate deflection of x = 0 with different elastic foundation coefficients k using 147 (49 × 3) grids are shown in Figure 9.

Figure 9.

I-shaped plate deflection of x = 0.

Conclusion

Maximum deflections w of square plate and rectangular plate shown in Tables 1 and 2 are convergent with the increasing number of grids, and the results are compared with those of analytic solution, which shows the convergence and high accuracy of the proposed method. Maximum deflections w of square plate and rectangular plate shown in Tables 1 and 2 are decreasing with the increasing of the elastic foundation coefficient k. Deflection diagrams of square plate deflection of x = 0, rectangular plate deflection of x = 0, trapezoidal plate deflection of x = 0, L-shaped plate deflection of x = −0.0536, and I-shaped plate deflection of x = 0 shown in Figures 2, 3, 5, 7, and 9 demonstrate that deflections are bigger in middle and become smaller in the sides. Deflection diagrams also show that they are decreasing with the increasing of the elastic foundation coefficient k. By describing the shape of square plate, rectangular plate, trapezoidal plate, L-shaped plate, and I-shaped plate, it demonstrates R-function theory is effective to express complicated shape. The advantages of the present solutions are to solve irregular shape plate problem such as L-shaped thin plate in Example 4 which is a six-sided shape and I-shaped thin plate in Example 5 which is a 12-sided shape, and the irregular shape plates have some extension of the boundary extends into it.

In this article, the convergence of the method is discussed. The comparison with the analytic solution shows a good agreement, and it demonstrates the feasibility and efficiency of the method in this article. The proposed method as shown in Appendix 1 applied to solve the elastic bending problems of simply supported thin plates on Pasternak foundation with various shapes is effective. R-function theory method is a very effective method to solve boundary value problems, which needs to be further developed in the field of cross-science.

Footnotes

Appendix 1

Handling Editor: James Baldwin

Data availability statement

The data used to support the findings of this study are available from the corresponding author upon request.

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The authors gratefully acknowledge the financial support provided by the Science and Technology Scheme of Guangzhou City (no. 201904010141).

ORCID iD

Shanqing Li

References

Dou

Z-C

Chen

, et al. Analysis on the method of fundamental solutions for biharmonic equations. Appl Math Comput 2018: 339: 346–366.

Assari

Asadi-Mehregan

Cuomob

A numerical scheme for solving a class of logarithmic integral equations arisen from two-dimensional Helmholtz equations using local thin plate splines. Appl Math Comput 2019; 356: 157–172.

Kravcenko

Merta

Zapletal

Distributed fast boundary element methods for Helmholtz problems. Appl Math Comput 2019; 362: 124503.

Cai

Fourier series solution for bending problem of the rectangular plates with free edges on Pasternak foundation. J Guizhou Univ Technol Nat Sci 1998; 27: 21–31 (in Chinese).

Yang

Shen

H-S

Zhao

Vibration of laminated plates with elastic rotational constraints and resting on Pasternak foundations. J Shang Jiao Tong Univ 2001; 35: 496–498 (in Chinese).

Wang

Z-W.

Solution to bending of rectangular thin plates simple supported on four edges on Pasternak model of foundation. J Luoyang Univ Soc Sci 2004; 19: 74–77 (in Chinese).

Liu

Peng

J-S.

Differential quadrature method for bending problem of rectangular plates with free edges on Pasternak foundation. J Chongq Inst Technol Nat Sci 2009; 23: 44–49 (in Chinese).

Liu

Axisymmetric large deflection of circular plate on two parameter elastic foundation. Chin J Appl Mech 1995; 12: 96–99 (in Chinese).

Piyush

Mohammad Sikandar

Vinayak

Vibration analysis of a thin functionally graded plate having an out of plane material inhomogeneity resting on Winkler–Pasternak foundation under different combinations of boundary conditions. Proc IMechE 2019; 233: 2636–2662.

10.

Chen

Meng

Yang

Exact nonstationary responses of rectangular thin plate on Pasternak foundation excited by stochastic moving loads. J Sound Vib 2018; 412: 166–183.

11.

Hosein

Fatemeh

Mohammad

An exact solution for stability analysis of orthotropic rectangular thin plate under biaxial nonlinear in-plane loading resting on Pasternak foundation. J Brazil Soc Mech Sci Eng 2015; 39: 953–967.

12.

Joodaky

A semi-analytical study on static behavior of thin skew plates on Winkler and Pasternak foundations. Int J Mech Sci 2015; 100: 322–327.

13.

Yan

Feng

X-T

Zhou

Dual reciprocity hybrid radial boundary node method for Winkler and Pasternak foundation thin plate. Arch Appl Mech 2013; 83: 225–239.

14.

Ehsan

Mohammad Mahdi

Mohammad

Vibration analysis of thin plates resting on Pasternak foundations by element free Galerkin method. Shock Vib 2013; 20: 309–326.

15.

Kiani

Bagherizadeh

Eslami

MR.

Thermal buckling of clamped thin rectangular FGM plates resting on Pasternak elastic foundation (Three approximate analytical solutions). Zeits Angewandte Math Mech 2011; 91: 581–593.

16.

Dahalan

Saudi

Sulaiman

, et al. Autonomous navigation in static indoor environment via rotated Laplacian operator. AIP Conf Proc 2018; 1974:020035.

17.

Koyunbakan

Pablo Pinasco

Scarola

Energy dependent potential problems for the one dimensional p-Laplacian operator. Nonlin Analy Real World Appl 2019; 45: 285–298.

18.

Mahata

Saha

Kar

, et al. A metaheuristic optimization approach to discretize the fractional order Laplacian operator without employing a discretization operator. Swarm Evolut Comput 2019; 44: 534–545.

19.

Tan

Solutions of fractional differential equations with p- Laplacian operator in Banach spaces. Bound Value Prob 2018; 2018: 15.

20.

Bai

Existence and uniqueness of solutions for fractional boundary value problems with p-Laplacian operator. Adv Diff Equat 2018; 2018: 4.

21.

Fan

Zhou

R-G

, et al. Quantum image edge extraction based on Laplacian operator and zero-cross method. Quant Inform Process 2019; 18: 27.

22.

Multiple positive solutions for nonlinear mixed fractional differential equation with p-Laplacian operator. Adv Diff Equat 2019; 2019: 112.

23.

Ling

Wei

Stereo-Matching algorithm using weighted guided image filtering based on Laplacian of Gaussian operator. Laser Optoelectron Progress 2019; 56: 226–232 (in Chinese).

24.

Tao

Liang

, et al. Damage buildings detection of post-disaster based on local Laplacian operators. Geomat Spat Inform Techol 2018; 41: 21–24 (in Chinese).

25.

Tian

Z-C

Huang

H-T

, et al. Analysis of the method of fundamental solutions for the modified Helmholtz equation. Appl Math Comput 2017; 305: 262–281.

26.

Chen

Xie

Luo

A wavelet collocation method for boundary integral equations of the modified Helmholtz equation. Appl Math Comput 2018; 321: 300–312.

27.

Wang

Fan

C-M

Hua

, et al. Localized MFS for the inverse Cauchy problems of two-dimensional Laplace and biharmonic equations. Appl Math Comput 2020; 364: 124658.

28.

S-Q

Yuan

Quasi-Green’s function method for free vibration of clamped thin plates on Winkler foundation. Appl Math Mech 2011; 32: 265–276.

29.

S-Q

Yuan

Quasi-Green’s function method for free vibration of simply- supported trapezoidal shallow spherical shell. Appl Math Mech 2010; 31: 635–642.

30.

Yuan

Green quasifunction method for free vibration of simply-supported trapezoidal shallow spherical shell on Winkler foundation. Acta Mechanica Solida Sinica 2010; 23: 370–376.

31.

Yuan

Green quasifunction method for free vibration of clamped thin plates. Acta Mechanica Solida Sinica 2012; 25: 37–45.

32.

Ahmad

Jang

Carrasco

, et al. An efficient iterative method for computing deflections of Bernoulli-Euler-von Karman beams on a nonlinear elastic foundation. Appl Math Comput 2018; 334: 269–287.

33.

Kurpa

Shmatko

Timchenko

Free vibration analysis of laminated shallow shells with complex shape using the R-functions method. Comp Struct 2010; 93: 225–233.

34.

Awrejcewicz

Kurpa

Osetrov

Investigation of the stress-strain state of the laminated shallow shells by R-functions method combined with spline-approximation. Zeits Angewandte Math Mech 2011; 91: 458–467.

35.

Kurpa

Shmatko

Awrejcewicz

Vibration analysis of laminated functionally graded shallow shells with clamped cutout of the complex form by the Ritz method and the R-functions theory. Latin Am J Solids Struct 2019; 16: 1–16.