Published: 26 November 2020

Bending of a circular flexible plate taking into account the initial curvature

G. U. Mamatova1
А. K. Sugirbekova2
N. M. Bekaulov3
1Department of Chemical Physics and Materials Science, Al-Farabi Kazakh National University, Almaty 050040, Kazakhstan
2Department of Chair of Mathematics, Physics and Informatics Teaching Methods, Abai Kazakh National Pedagogical University, Almaty 050010, Kazakhstan
3Department of Computer Engineering and Information Security, International IT University, Almaty 050000, Kazakhstan
Corresponding Author:
G. U. Mamatova
Editor's Pick
Views 302
Reads 136
Downloads 938

Abstract

In this article, an analytical solution of the problem of bending a round plate with an initial curvature is obtained. A round plate with an initial curvature as a structural element is widely used, the calculation of which leads to many questions related to the design of round Foundation plates, turbine disks, flexible shaft connections and other. The nonlinear theory of round flexible plates provides the key to explaining the process of loss of stability, which often leads to complete destruction of the structure. Therefore, there is a high need for analytical methods for solving problems of calculating the stress-strain state of round flexible plates, taking into account the initial curvature. These problems are mathematically reduced to differential equations with variable coefficients, the exact solution of which, as a rule, does not exist. Therefore, the construction of analytical solutions to these problems is very relevant.

Bending of a circular flexible plate taking into account the initial curvature

1. Introduction

The problem of bending elastic round plates of initial curvature is one of the actual problems of the technical theory of elasticity. A round plate with an initial curvature as a structural element is widely used, which is calculated by many questions related to the design of round Foundation plates, turbine disks, flexible shaft connections, hydraulic machine blades, disk springs, etc. [1], The development of modern practice requires researchers and designers to create methods for solving a large number of strength problems related to the variability of thickness, elastic modulus, Poisson's ratio, and the presence of initial curvature in aggregates. There is a high need for analytical and approximate analytical methods for solving problems of calculating the stress-strain state of inhomogeneous plates, taking into account the initial curvature. Existing studies of symmetrical bending of round plates with initial curvature are of a particular nature. Asymmetric bending of round plates of variable thickness has been studied even less; here, in addition to solutions for hyperbolic profile plates, there are practically no solutions for any other profiles of practical significance. These problems are mathematically reduced to differential equations with variable coefficients, the exact solution of which, as a rule, does not exist. Therefore, the construction of analytical solutions to these problems is very relevant [2]. This work substantially fulfills the specified gap. Using the method of partial discretization of differential equations, based on the theory of generalized functions, an analytical solution to the problem of bending elastic round plates with account for the initial curvature is obtained.

2. Materials and methods

Let’s consider a round plate that is fixed along the contour in a certain way. We will assume that the plate is ring-shaped with external and internal radii rb, ra and thickness h.

Fig. 1Basic designations when considering a round plate

Basic designations when considering a round plate

In this case, the deflections are symmetrical with respect to the axis z and the deflection, as well as all other values that characterize the stressed and deformed state, will be functions only r. When constructing solutions to the solving equations, we use the partial discretization method, which allows us to construct an analytical solution to the problem of bending a round flexible plate with an initial deflection, fixed along the contour and subjected to a uniformly distributed transverse load of intensity q. In this case, the deflection, as well as all other values that characterize the stressed and deformed state, will be functions of only r.

The basic system of differential equations of the theory of a flexible round plate with initial deflection has the following form [3]:

1
Dddr2ω=ψ+hrdΦdrdωdr+dωindr,ddr2Φ=-Еr12dωdr2+dωindrdωdr,

where ψ is the load function equal to ψ=1r0rqrdr, Φ is the voltage function entered by the formulas:

σr=1rdΦdr, σφ=d2Φdr2,

and E is the elastic modulus, ω is the deflection, σr, σφ are radial and tangential stresses.

Using the method of partial discretization of differential equation [3], we write the second equation of the system Eq. (1) as:

2
d3Φdr3+1rd2Φdr2-1r2dΦdr=-12k=1nrk+rk+1Erk12dωkdr2+dωinkdrdωkdrδr-rk
-Erk+112dωk+1dr2+dωink+1drdωk+1drδr-rk+1.

The general solution of Eq. (2) will be:

dΦdr=C1r+C21r-14k=1nrk+rk+1Erk12dωkdr2+dωinkdrdωkdrr-rk2rHr-rk
-Erk+112dωk+1dr2+dωink+1drdωk+1drr-rk+12rHr-rk+1.

Consider a plate with a hole that is pivotally supported along the contour, then the boundary conditions are Eq. [4]:

3
1rdΦdrr=ra=0,
4
Mrr=rb=Dd2ωdr2+μrdωdrr=rb=0,
5
ωr=rb=0.

Due to the fact that the ra is not subject to a smallness constraint, we will assume an arbitrary constant C2 equal to zero. By virtue of condition Eq. (3), the solution of problem Eq. (2, 3) will be:

6
dΦdr=-E4k=1nrk+rk+11rk12dωkdr2+dωinkdrdωkdrr-rk2rHr-rk
-1rk+112dωk+1dr2+dωink+1drdωk+1drr-rk+12rHr-rk+1.

Substitute the resulting expression in the first equation of the system Eq. (1) and we will perform further discretization of the multiplier dω/dr in the right part of this equation. Then the first equation of the system Eq. (1) will take the form:

7
d3ωdr3+1rd2ωdr2-1r2dωdr=ψD-E4hDrk=1nrk+rk+11rk12dωkdr2+dωinkdrdωkdr
×r-rk2rHr-rk
-1rk+112dωk+1dr2+dωink+1drdωk+1drr-rk+12rHr-rk+1dωindr
-E8hDrdωdrk=1nrk+rk+11rk12dωkdr2+dωinkdrdωkdrr-rk2rHr-rk
-1rk+112dωk+1dr2+dωink+1drdωk+1drr-rk+12rHr-rk+1k=1nrk+rk+1δr-rk-δr-rk+1.

After the corresponding transformations, taking into account the properties of generalized functions and conditions Eq. (4), we get:

8
dωdr=r2Dψdr-12Drr2ψdr-r2Dr=rbψdr+1-μ1+μ1rb2r=rbr2ψdr
+Ehr8Dk=1nrk+rk+112dωkdr2+dωinkdrdωkdr1r2r2rk-rkdωindrHr-rkdr
-1rk-rkr2dωindrHr-rkdr-μ-1rb21+μr=rbr2rk-rkdωindrHr-rkdr
+r=rb1rk-rkr2dωindrHr-rkdr-12dωk+1dr2+dωink+1drdωk+1dr
×1r2r2rk+1-rk+1dωindrHr-rk+1dr-1rk+1-rk+1r2dωindrHr-rk+1dr
-μ-1rb21+μr=rbr2rk+1-rk+1dωindrHr-rk+1dr+r=rb1rk+1-rk+1r2dωindr
×H(r-rk+1+Ehr16D1+μk=2nrk+1-rk-11+rk2rb2+μ1-rk2rb2
-1+μ1-rk2r2Hr-rki=1k-1ri+ri+112dωidr2+dωinidrdωidr1ri-ri rk2
-12dωi+1dr2+dωini+1drdωi+1dr1ri+1-ri+1 rk2dωkdr-Ehr16D1+μrn-rn+1
×1+rn+12rb2+μ1-rn+12rb2-1+μ1-rn+12r 2Hr-rn+1
i=1nri+ri+112dωidr2++dωinidrdωidr1ri-ri rn+12
-12dωi+1dr2+dωini+1drdωi+1dr1ri+1-ri+1 rn+12dωn+1dr.

In this case q=const, the analytical expression of the rotation angle at points rk for has the expression:

9
dω1dr=qr1316D-qc23+μr116D1+μ,
10
dω2dr=qr216Dr22-c23+μ1+μ1-Ehr2r3-r132D1+μ1+r22c2+μ1-r22c2r1+r2dω1dr21r11-r12r22,

for j=3,n¯.

11
dωjdr=qrj16Drj2-c23+μ1+μ+Ehrj32D1+μk=2j-1rk+1-rk-11+rk2c2+μ1-rk2c2
-1+μ1-rk2rj2i=1k-1ri+ri+1dωidr21ri-rirk2-dωi+1dr21ri+1-ri+1 rk2dωkdr
/1-Ehrj32D1+μrj+1-rj-11+rj2c2+μ1-rj2c2
×i=1j-1ri+ri+1dωidr21ri-rirj2-dωi+1dr21ri+1-ri+1rj2.

Integrating Eq. (8) and taking into account condition Eq. (5), we write the solution of the first differential equation of the system Eq. (1) as:

12
ω=q64Dr4-rb4+qrb23+μ32D1+μrb2-r2+Eh8Dk=1nrk+rk+112dωkdr2
+dωinkdrdωkdr1rr2rk-rkdωindrHr-rkdr
-r1rk-rkr2dωindrHr-rkdrdr
-r=rb1rr2rk-rkdωindrHr-rkdr-r1rk-rkr2dωindrHr-rkdrdr
+1-μ21+μr2rb2-1r=rbr2rk-rkdωindrHr-rkdr+r22-rb22r=rb1rk-rkr2
×dωindrHr-rkdr-12dωk+1dr2+dωink+1drdωk+1dr1rr2rk+1-rk+1
×dωindrHr-rk+1dr-r1rk+1-rk+1r2dωindrHr-rk+1drdr
-r=rb1rr2rk+1-rk+1×dωindrHr-rk+1dr
-r1rk+1-rk+1r2dωindrHr-rk+1drdr+1-μ21+μr2rb2-1
×r=rbr2rk+1-rk+1dωindrHr-rk+1dr
+r22-rb22r=rb1rk+1-rk+1r2dωindrHr-rk+1dr
+Eh16Dk=2nrk+1-rk-11+1-μ1+μrk2c2r22-rk21+μ-1+rk2c2-r22-11-rk2r2Hr-rk+i=1k-1ri+ri+112Vi2+ViniVi1ri-rirk2
-12Vi+12+Vini+1Vi+11ri+1-ri+1rk2Vk.

Fig. 2 shows graphs of deflection changes for a flexible round plate with different initial deflections.

3. Results and discussion

To obtain a solution of the system of resolving differential equations of the problem under consideration, A. S. Volmir sets the deflection form and then uses the method of Bubnov-Galerkin. A system of equations is often solved by a small parameter method, when the multiplier before the nonlinear terms is small. If there is no such restriction on this multiplier, then getting a solution is problematic. Therefore, we can use the method of partial discretization of differential equations, the accuracy of the resulting solution for which has been repeatedly checked when solving problems related to obtaining plate bends.

The Fig. 2 shows the bending curve of the flexible plate for ωin=fin1-r2/rb22 individual initial values of the deflection boom fin.

Fig. 2Deflection change curves for: fin=0,2h; fin=0,5h; fin=0,8h; h= 0,02 m; q= 104 N/m2; E= 2⋅1011 N/m2; μ= 0,3.

Deflection change curves for: fin=0,2h; fin=0,5h; fin=0,8h;  h= 0,02 m; q= 104 N/m2; E= 2⋅1011 N/m2; μ= 0,3.

4. Conclusions

The existing method of finding solutions to the nonlinear differential equations of Bubnov-Galerkin [5], although it is very useful, but in the General formulation of the problem and the tendency of the number of linear combinations of a given linear independent system to infinity, cannot guarantee even a weak convergence of the approximate solution to the exact one. In this regard, the application of the partial discretization method to the system of nonlinear equations under consideration is very appropriate.

References

  • Reddy J. N. Theory and Analysis of Elastic Plates and Shells. CRC Press, Taylor and Francis, 2007.
  • Szilard R. Theories and Applications of Plate Analysis: Classical, Numerical and Engineering Methods. John Wiley & Sons, 2004.
  • Volmir A. S. Flexible Plates and Shells. Moscow, State Technical Publishing House, 1956.
  • Tyurekhodgaev, A. N., Mamatova G. U., Kalzhanova G. K. Bending of Inhomogeneous and Nonlinear Plates and Shells. LAP LAMBERT Academic Publishing is a trademark of: OmniScriptum GmbH & Co, Germany, 2015.
  • Smolina I. Y., Puteeva L. E., Shilnikov S. M. Calculation of Plate Bending by the Method of Bubnov-Galerkin. Tomsk: Publishing House of Tomsk State University of Architecture and Building Engineering, 2009.

About this article

Received
02 October 2020
Accepted
26 October 2020
Published
26 November 2020
SUBJECTS
Mathematical models in engineering
Keywords
the bending of the plate
stress-strain state
nonlinear differential equations
partial discretization method
class of generalized functions