Published: 26 November 2020

Mathematical modeling of the solution to the problem of plate bending given various fastenings along the contour and various loads

Serikbay Toibayev1
Assemgul Uderbayeva2
Zhanar Bidakhmet3
Ayaulym Rakhmatulina4
Algazy Zhauyt5
1, 2, 3, 4Department of Engineering Graphics and Applied Mechanics, Almaty Technological University, Almaty, 050012, Kazakhstan
5Department of Electronics and Robotic, Almaty University of Power Engineering and Telecommunications, Almaty, 050013, Kazakhstan
Corresponding Author:
Serikbay Toibayev
Views 86
Reads 26
Downloads 954

Abstract

The classical equations of plate bending in various fastenings along the contour and various loads are considered. Using the decomposition method, three auxiliary problems are formulated to determine the deflection of the plate. The task is to determine the individual terms of the decomposition coefficients anmi from the individual terms of the load series. Then the system of algebraic equations is obtained to determine anmi. If the plate contour is pivotally fixed along the contour, then the system of equations takes a known form and the expression for the deflection function is explicitly found. Approximate analytical solutions are obtained for the functions of deflections, internal force factors and stresses under various boundary conditions of the plate.

1. Introduction

In our century, with the increasing complexity of the forms of building structures, the emergence of aircraft manufacturing, and the diverse demands of mechanical engineering, the role of elasticity theory methods has changed dramatically [1]. Now they form the basis for constructing practical methods for calculating deformable bodies and systems of bodies of various shapes [2]. At the same time, modern calculations take into account not only the complexity of the body shape and the variety of effects (power, temperature, etc.), but also the specificity of the physical properties of the materials from which the bodies are made [3]. The fact is that in modern designs along with traditional materials (steel, wood, concrete, etc.), new materials are widely used, in particular composites with a number of specific properties [4]. Thus, reinforcing polymers with fibers from high-strength materials makes it possible to obtain new lightweight structural material having high strength properties that exceed even the strength of modern steels [5]. But the presence of a polymer base gives such a composite material, viscous properties in addition to elasticity, which must be taken into account in the calculations [6]. Exact solutions in the analytical form of the equations of the elasticity theory subject to boundary conditions, which constitutes the so-called boundary value problem, are possible only in some special cases of loading bodies and the conditions for their fastening [7]. Therefore, for engineering practice, approximate, but fairly general methods for solving problems of the applied theory of elasticity are especially important [8].

2. Materials and methods

Consider a rectangular plate of constant thickness, having geometric dimensions in plan 0xl1, 0yl2. We assume that the plate along the contour in the general case is fixed elastically. To determine the bends, we have the differential equation describing the deformations:

1
D4Wx4+24Wx2y2+4Wy2=qx,y,

where D=E1δ312=Eδ3121-μ2 bending stiffness, qx,y external load. Boundary conditions at elastic fastening of the plate along the contour:

2
W=0, l1k2Wx2+1-kWx=0, x=0, x=l1,W=0, l2k2Wy2+1-kWy=0, y=0, y=l2.

Then boundary value problem Eq. (1) and Eq. (2) transform to the form (λ=l1/l2):

3
4να4+2λ24να2β2+λ44νβ4=qD.

The purpose of the work is to determine the bending of the plate. In accordance with the decomposition method Eq. (3), we formulate three auxiliary problems:

Problem 1:

4
4ν1α4=f1α,β, ν1=0, πk2ν1α2+1-kν1α=0, α=0, α=π.

Problem 2:

5
λ44ν2β4=f2α,β, ν2=0, πk2ν2β2+1-kν2β=0, β=0, β=π.

Problem 3:

6
2λ24ν3α2β2=qD-f1α,β-f2α,β.

Conditions:

7
ν1=ν2=ν3=ν,

needs to perform approximately. Let’s decompose RHS of Eqs. (4-6) to the Fourier series:

8
fiα,β=n=1m=1anmisinnαsinmβ, i=1,2,

where anmi unknown constants, approximating unknown functions of the decomposition method.

Let’s decompose the load function q to the Fourier series:

9
qα,β=n=1m=1qnmsinnαsinmβ,

where qnm are coefficients of the load decomposition, and are determined using formula:

10
qnm=4π20π0πqα,βsinnαsinmβdαdβ.

Then the auxiliary problems Eqs. (6-8) are written as follows.

Problem 1:

11
4ν1α4=n=1m=1anm1sinnαsinmβ, ν1=0, πk2ν1α2+1-kν1α=0,
α=0, α=π.

Problem 2:

12
λ44ν2β4=n=1m=1anm2sinnαsinmβ, ν2=0, πk2ν2β2+1-kν2β=0,
β=0, β=π.

Problem 3:

13
2λ24ν3α2β2=1Dn=1m=1qnmsinnαsinmβ
-n=1m=1anm1sinnαsinmβ-n=1m=1anm2sinnαsinmβ.

And meeting the boundary conditions Eq. (11), we find:

14
ν1=n=1m=1anm1n31nsinnα+4α3kk-1π21-2k-11k2+α24k-5k2+1π1-2k-11k2
+αk2-11-2k-11k2sinmβ,

and solution of the boundary value problem Eq. (12) will be similar:

15
ν2=n=1m=1anm2λ4m31msinmα+4β3kk-1π21-2k-11k2+β24k-5k2+1π1-2k-11k2
+βk2-11-2k-11k2sinnα,

where anmi bend surface decomposition coefficients to be determined. For this, the term of the series of n, m number for the load Eqs. (9) and the same term of the series of bend Eqs. (14-15). According to Eqs. (7) assume:

16
ν3π2,π2=12ν1π2,π2+ν2π2,π2, ν1π2,π2=ν2π2,π2.

Now the problem is to determine individual terms of decomposition coefficients anmi from individual terms of load series Eq. (9). Then to determine anmi from Eq. (13) and Eq. (16) we will derive a system of algebraic equations. So, functions Eqs. (14-15) are the solution to the problem posed, since they meet the conditions on the plate contour and, when choosing the coefficients, derived by solving the system of algebraic equations. We introduce the notations:

3. Results and discussion

Fig. 1Computed plot of the coordinates α and β

Computed plot of the coordinates α and β
Computed plot of the coordinates α and β
Computed plot of the coordinates α and β
Computed plot of the coordinates α and β
Computed plot of the coordinates α and β
Computed plot of the coordinates α and β

Under the action of a uniform load of the entire area we have:

Gnm=16qDnmπ2sinnπ2sinmπ2.

Fig. 2Computed plot of the evenly distributed load

Computed plot of the evenly distributed load
Computed plot of the evenly distributed load
Computed plot of the evenly distributed load
Computed plot of the evenly distributed load
Computed plot of the evenly distributed load
Computed plot of the evenly distributed load

Under the action of a load changing along a parabola in the directions of the axis α of the entire area we have:

Gnm=32qDnm3π41--1nsinnπ2sinmπ2.

Fig. 3Computed plot of the changing parabola

Computed plot of the changing parabola
Computed plot of the changing parabola
Computed plot of the changing parabola
Computed plot of the changing parabola
Computed plot of the changing parabola
Computed plot of the changing parabola

When the action is changing along an isosceles triangle in the directions of the axis α of the entire area we have:

Gnm=16qDnm2π3sin2nπ2sinmπ2.

Fig. 4Computed plot of the isosceles triangle

Computed plot of the isosceles triangle
Computed plot of the isosceles triangle
Computed plot of the isosceles triangle
Computed plot of the isosceles triangle
Computed plot of the isosceles triangle
Computed plot of the isosceles triangle

4. Conclusions

The paper obtained equations of plate bending, internal force factors and stresses, and the following new scientific results:

1) approximate analytical solutions were obtained for the functions of bends, the formulas of internal force factors and stresses given various boundary conditions of the plate and load;

2) the results of a numerical experiment for the plate were obtained given various fastenings of the edges and load;

3) mechanical effects, conditions of fastening along the contour and load were revealed.

References

  • Khairnasov K. Z. Mathematical modeling of dynamic stability of shell structures under large elastoplastic deformations. IOP Conference Series: Journal of Physics, Vol. 991, 2018, p. 012042.
  • Omarov T., Tulegenova K., Bekenov Y., Abdraimova G., Zhauyt A., Ibadullayev M. Determination of reduced mass and stiffness of flexural vibrating cantilever beam. Journal of Measurements in Engineering, Vol. 6, Issue 1, 2018, p. 1-9.
  • Mogilevich L. I., Popov V. S., Popova A. A., Christoforova A. V. Mathematical modeling of hydroelastic walls oscillations of the channel on Winkler foundation under vibrations. Vibroengineering Procedia, Vol. 8, 2016, p. 294-299.
  • Kondratov D. V., Kalinina A. V., Mogilevich L. I., Popova A. A., Kondratova Y. N. Mathematical model of elastic ribbed shill dynamics interaction with viscous liquid under vibrations. Vibroengineering Procedia, Vol. 8, 2016, p. 300-305.
  • Andreev V., Chepurnenko A. On the bending of a thin plate at nonlinear creep. Advanced Materials Research, Vol. 900, 2014, p. 707-710.
  • Dickey R. W. Nonlinear bending of circular plates. Journal of Applied Mathematics, Vol. 30, Issue 1, 1976, p. 1-9.
  • Reissner E. On the theory of bending of elastic plates. Journal of Mathematics and Physics, Vol. 23, 1944, p. 184-191.
  • Charles Chinwuba Ike Mathematical solutions for the flexural analysis of Mindlin’s first order shear deformable circular plates. Mathematical Models in Engineering, Vol. 4, Issue 2, 2018, p. 50-72.

Cited by

2021 IEEE International Symposium on the Physical and Failure Analysis of Integrated Circuits (IPFA)
Bao Yazhou | Yang Gang | Jia Xingge | Song Yong | Yu Hongjun | Che Chuncheng | Xue Hailin
(2021)

About this article

Received
20 February 2020
Accepted
26 May 2020
Published
26 November 2020
SUBJECTS
Mathematical models in engineering
Erratum
The institutional affiliation was incorrect for Bidakhmet Zhanar in the paper finally approved (after the acceptance) by the authors. For more information read Editor's Note.
Keywords
mathematical modeling
high-strength
composite material
contour
solving problems