CC BY
25Magazine of Civil Engineering. 2023. 118(1). Article No. 11811
Magazine of Civil Engineering issn
2712-8172
journal homepage: http://engstroy.spbstu.ru/
Research article UDC 69.04
DOI: 10.34910/MCE.118.11
Nonlinear vibrations and dynamic stability of viscoelastic anisotropic fiber reinforced plates
B.Kh. Eshmatov1 , R.A. Abdikarimov2 , M. Amabili , N.I. Vatin4
1 Almalik Branch, National University of Science and Technology MISiS, Almalik, Uzbekistan
2 Tashkent Institute of Architecture and Civil Engineering, Tashkent, Uzbekistan
3 McGill University, Montreal, Canada
4 Peter the Great St. Petersburg Polytechnic University, St.Petersburg, Russia [email protected]
Keywords: nonlinear vibrations, dynamic stability, viscoelastic anisotropic fiber-reinforced plate, integro-differential equations, weakly-singular Koltunov-Rzhanitsyn kernel
Abstract. Fiber-reinforced plastic composites are one of the most widely used composite materials because they balance well between properties and cost. Despite their widespread use in the aviation and automotive industries, there is currently a lack of effective mathematical models for their calculation under various dynamic loads. The research object of this work is an anisotropic viscoelastic fiber-reinforced simply supported rectangular plate. Two dynamic problems are considered: vibrations of the plate under the influence of a uniformly distributed static load; stability of the plate compressed in one direction. Within the Kirchhoff-Love hypothesis framework, a mathematical model was built in a geometrically nonlinear formulation, taking into account the tangential forces of inertia. By the Bubnov-Galerkin method, based on a polynomial approximation of the deflection and displacement, the problem was reduced to solving systems of nonlinear ordinary integro-differential equations. With a weakly singular Koltunov-Rzhanitsyn kernel with variable coefficients, the resulting system was solved by a numerical method based on quadrature formulas. By using experimental studies, considering the directions of the fibers, the values of the rheological parameters of some plastic materials (KAST-V and EDF) were obtained. The plate's dynamical behavior was investigated depending on the plate's geometric parameters, viscoelastic and inhomogeneous material properties. Results show the importance of taking into account the viscoelastic properties of the material when solving dynamic problems of anisotropic reinforced plates made of composite materials. In particular, when studying the problem of dynamic stability of an anisotropic reinforced plate made of KAST-V, the results obtained in elastic and viscoelastic formulations in some cases differ from each other by more than 20 %.
Funding: This research was supported by a grant from the Russian Science Foundation No. 21-19-00324, https://rscf.ru/project/21-19-00324/
Citation: Eshmatov, B.Kh., Abdikarimov, R.A., Amabili, M., Vatin, N. Nonlinear vibrations and dynamic stability of viscoelastic anisotropic fiber reinforced plates. Magazine of Civil Engineering. 2023. 118(2). Article no. 11811. DOI: 10.34910/MCE.118.11
1. Introduction
Mechanics of composite materials has made great progress in recent years. The application of composite materials in engineering structures considerably improved structures' operational characteristics and, in several cases, create structures that could not be realized with traditional materials. Polymer
© Eshmatov, B.Kh., Abdikarimov, R.A., Amabili, M., Vatin, N., 2023. Published by Peter the Great St. Petersburg Polytechnic University.
composite materials are widely used in diverse engineering fields. A large class of engineering structures relies on the use of composite materials and ceramics (for example, operating under conditions of high temperatures, contact pressures, etc.). Therefore, developing efficient algorithms for solving nonlinear problems of dynamis stability of shells, panels, and plates made of composite materials is important in today's industry.
The development of high-modulus and high-strength boron and graphite fibers joined together in plastic or metal binders arouses great interest in structural elements made of layers, each of which has different mechanical properties. In most cases, these composite materials include many layers of fibers with specific stacking sequences. The dimensions of the structural elements are such that each layer can be adequately modeled as a homogeneous anisotropic layer. Due to the unique characteristics of modern composite materials, interest in the problems of studying the dynamic behavior of anisotropic reinforced layered structures has not weakened in recent years.
An overview analysis of the using various natural fibers to form composites with various matrix materials and the effect of fiber length, chemical processing, and compositions on the dynamic mechanical characteristics of composites is available in [1].
An overview of works related to the study of typical functional capabilities of composites reinforced with natural fibers, including thermal insulation properties, biodegradable properties, and vibration damping properties, are given in [2]. This article concludes that composite materials that fully utilize their functionality rather than mechanical characteristics are important for the future development of composites reinforced with natural fibers. The analysis of some published literature on vibration damping applied to composites reinforced with natural fibers by experimental and numerical methods was carried out here. The effect of reinforcement with natural fibers on the thermosetting and thermoplastic polymer matrix vibration response was considered.
The article [3] investigates the effect of transverse normal and shear deformations on fiber-reinforced viscoelastic beams resting on two-parameter (Pasternak's) elastic foundations. The results were obtained using the refined theory of the sinusoidal shear deformation beam theory and compared with the results obtained using the simple theory of the sinusoidal shear deformation beam theory, the Timoshenko firstorder shear deformation beam theory, and the classical Euler-Bernoulli beam theory.
The characteristics of free vibration of rectangular plates reinforced with continuous grading fiber (CGFR) based on elastic bases were studied based on the three-dimensional linear theory of elasticity and the theory of small deformations [4]. The foundation was described by the Pasternak model or a two-parameter model. The CGFR plate was simply supported at the edges and was assumed to have an arbitrary change in the volume fraction of the fiber in the direction of thickness. Displacement functions that satisfy simply supported boundary conditions were used to reduce the equilibrium equations to a set of coupled ordinary differential equations with variable coefficients, which in turn have been solved by differential quadrature (DQM) to obtain natural frequencies. Studies have shown that this method has a high convergence rate, stable numerical performance, and accuracy. The results for a CGFR plate with an arbitrary change in the volume fraction of the fiber in the thickness of the plate were compared with a discrete layered composite plate. They also showed that the use of CGFR plates with graduated volume fractions of fibers has advantages over traditional discretely laminated plates.
The static response of an inhomogeneous fiber-reinforced viscoelastic multilayer plate was investigated in [5] by using the first-order shear deformation theory. Several types of sandwich plates were considered, taking into account the symmetry of the plate and the thickness of each layer. In addition, two cases were considered depending on the viscoelastic material that enters the core or the surface of multilayer plates. To solve equations determining the bending of simply supported inhomogeneous fiber-reinforced viscoelastic multilayer plates, the effective modules method and Ilyushin's approximation method were used. Numerical calculations were carried out to study the effect of the time parameter on deflections and stresses at different values of the aspect ratio, the ratio of sides to the thickness, and the constitutive parameter.
Amabili [6] investigated the dynamic stability of orthotropic thick plates subjected to periodic uniaxial stress and bending stress, considering the inertia of rotation and transverse stress. Based on Bolotin's method, the regions of dynamic instability of plates reinforced with graphite and fiberglass were estimated by solving eigenvalue problems. The influence of material properties and load parameters on the instability region and the index of dynamic instability of orthotropic plates was discussed.
The distorted modes in geometrically nonlinear forced oscillations of angular layered circular cylindrical shells were investigated using the highest order Amabili-Reddy's shear deformation theory [7]. Harmonic forces were applied in the radial direction and simply supported boundary conditions were assumed. The equations of motion were obtained using an energy approach based on Lagrange equations,
which preserves dissipation. Numerical results were obtained using the method of pseudo-arc continuation and bifurcation analysis.
Anisotropic composite plates were evaluated with nanofiber-reinforced matrices (NFRM) [8]. Seven different loading conditions were considered: three for uniaxial loading, three for paired combined loading, and one with three combined loads. The anisotropy had a complete stiffness matrix of 6 by 6, which was satisfied and solved using the Galerkin algorithm.
The influence of some geometric and material parameters on the free vibrations of FML plates (fibrous metal laminates) was studied in [9]. For the analytical solution of the governing equations of the composite plate, the first-order shear deformation theory (FSDT) and the Fourier series method were used.
In combination with the elastic constitutive model and fiber direction, a plastic square plate with two materials (orthotropic and anisotropic) was simulated [10]. Anisotropic characteristics of the thermal field, pressure field, and vibration characteristics of fiber-reinforced PA66 were analyzed using simulation comparisons.
A new nonlinear model of vibrations of reinforced composite thin plates in a thermal environment was proposed [11] by introducing nonlinear thermal and amplitude fitting coefficients. Based on the classical laminated plate theory, the complex module approach, the power function, and the Ritz methods, dynamic governing equations in high-temperature environments for solving nonlinear natural frequencies, vibration characteristics, and damping parameters were obtained. Experimental tests of thin plates made of carbon-epoxy composite CF130 were carried out to test this model.
The article [12] presents a new nonlinear model of vibrations of composite plate structures with amplitude-dependent properties, taking into account the nonlinear characteristics of stiffness and damping by introducing the nonlinearity of the material into the dynamic modeling of fiber-reinforced composite thin plates. In this model, the elastic modulus and loss coefficients were expressed as a function of the strain energy density based on the nonlinear Jones-Nelson material model. These elastic moduli and loss coefficients were characterized as a function of the maximum dimensionless strain energy density using specific parameters at different excitation amplitudes. The power function fitting method was then used to determine the model's nonlinear stiffness and damping parameters. The nonlinear natural frequencies, vibration characteristics, and damping coefficients of the TC300 carbon-epoxy composite thin plate were calculated and measured in a case study.
A Fourier series method based on Mindlin theory and Hamilton's variational principle was proposed [13] for modeling and analyzing the vibration of composite multilayer plates with arbitrary boundary conditions, in which vibrational displacements were sought as a linear combination of a double series of Fourier cosines and auxiliary series functions. Three types of springs with constraints were introduced, and a vibration model was created by combining the Hamilton energy principle to create a general structural model of composite layered materials. The influence of boundary conditions, the laying angle, and the laying layer on the vibration characteristics were analyzed.
The nonlinear dynamic model in a thermal environment was established for modeling a reinforced composite structure made of thin plates [14]. This model uses Hamilton's principle in combination with the classical theory of layered plates, the complex modulus method, and the strain energy method. Nonlinear dependences between elasticity modules, Poisson coefficients, loss coefficients, and temperature change were expressed by the polynomial method. Dynamic equations in a high-temperature environment were derived to solve the inherent characteristics, dynamic characteristics, and damping parameters, considering temperature-dependent properties. The principle of identification of the corresponding fitting coefficients in the theoretical model was illustrated. An experimental test of a thin plate made of carbon-epoxy composite TC500 was carried out to substantiate the correctness of the obtained results. The results of the developed model and experimental tests showed good consistency.
The buckling of a multilayer composite reinforced plate, with porosity effects, was analyzed according to the theory of first shear deformation plates [15]. The properties of the plate material were considered orthotropic properties. Three different models were used considering the effect of porosity in a multilayer composite plate. The Navier procedure was used to solve the problem for a simply supported plate. The influence of porosity coefficients, porosity models, fiber orientation angles, and layer sequences on critical buckling loads was discussed.
Othman, Abouelregal, and Said [16] analytically studied the influence of variable thermal conductivity and initial stress on a reinforced transversely isotropic thick plate. A model of the problem of generalized thermoelasticity with phase delays in an isotropic elastic medium with temperature-dependent mechanical properties was developed. It was assumed that the upper surface of the plate was thermally insulated with a given surface load, and the lower surface of the plate rested on a rigid base. Normal mode analysis was used to obtain analytical expressions of displacement components, power voltage, and temperature
distribution. The results obtained in the presence and absence of initial stress and reinforcement were compared.
The nonlinear Jones-Nelson theory combined with the classical theory of layered plates by Yadav, Amabili, et al. [17]. The polynomial fitting method, the strain energy method, and a damping model of a fiber-reinforced composite thin plate were used. The methods and model were created considering the amplitude-dependent properties. The elastic modulus was expressed as the strain energy density function in this model. The loss coefficients in the longitudinal, transverse, and shear directions were expressed as functions of the excitation amplitude. Three plates made of carbon-epoxy composite TC300 were taken as research objects for the case study. One of them was used to determine the amplitude-dependent coefficients of loss coefficients in fiber-reinforced composites by combining the least-squares method with the polynomial fitting method. The other two plates were used to verify the correctness of the theoretical model. The developed model results, taking into account the amplitude dependence and experimental verification, showed good consistency.
The strict polynomial of a higher order in the thickness coordinate was introduced to develop a theory for doubly curved multilayer composite shells, considering the thickness and shear deformations [18]. By applying the condition of zero transverse normal and tangential stresses on the upper and lower surfaces of the shells, the theory of thickness deformation and shear of the third order with six kinematic parameters was derived. At the same time, the nonlinear terms in all kinematic parameters were preserved. The accuracy of the proposed theory has been tested for static and dynamic reference cases. The results obtained based on the proposed theory were compared with the results obtained using a more complex and computationally time-consuming nine-parameter theory. Isotropic and cross-layered circular cylindrical shells under the action of radial forces and pressure and nonlinear forced oscillations of the cross-layered shell under harmonic radial excitation were considered.
The vibrations of circular cylindrical shells made of a carbon nanotube (CNT) fiber-reinforced composite (CNT-FRC) were studied in a geometrically nonlinear formulation [19]. Vibrations were created due to radial harmonic force, and the process itself was considered viscous structural damping. Elastic properties of randomly distributed CNTs in a polymer matrix (i.e., hybrid matrix) were calculated according to the Eshelby-Mori-Tanaka/Voigt scheme to account for the effect of CNT agglomeration in a hybrid matrix. The resulting hybrid matrix was reinforced with aligned fibers to prepare the CNT-FRC shell plate; its effective properties were evaluated using the Halpin Tsai homogenization approach. The CNT-FRC shell was modeled using von Karman geometric nonlinearity and first-order shear deformation theory (FSDT). According to the Hamilton principle, nonlinear control partial differential equations (PDEs) of CNT-FRC shells were derived. These PDEs were discretized into ordinary differential equations (ODES) using the Galerkin method. The ODES were solved by the incremental harmonic balance (IHB) method combined with the arc length extension method to obtain the frequency-amplitude characteristic of the shell. The influence of various types of CNT agglomeration models, CNT mass fraction, agglomeration parameters, and the sequence of laying layered materials on frequency-amplitude curves corresponding to forced and free nonlinear oscillations of the CNT-FRC shell was studied in detail.
The main purpose of the study [20] was to determine the effect of fiber orientation on the reduction of applied voltage for modeling the structure of a dielectric elastomer and improving its mechanical characteristics. Based on the nonlinear continuum mechanics and large inelastic deformations, the defining relations and equations determining the behavior of viscoelastic dielectric elastomers under harmonic electrical load were extracted and analyzed in various states. The obtained numerical results, such as phase diagrams, frequency amplitudes, and oscillations, illustrate the dynamic behavior of an anisotropic dielectric elastomer with different fiber orientations. Hyperelastic and rheological models were used to consider the viscoelastic properties of the material in combination with electrical coupling.
The analysis of many experimental and fundamental studies shows that most composite materials have clearly pronounced viscoelastic properties [21-23].
Most researchers who have attacked the class of problems mentioned above have considered the solution of problems with such a mathematical statement in an elastic case. In these works, only some properties of construction materials were taken into consideration: the problems were solved either in a linear context or in elastic wave propagation was not taken into account. Even if the problems were solved using a viscoelastic formulation, in many cases, the viscoelastic characteristics of the material were only taken into account in a restricted context. In these cases, viscoelastic properties were addressed by employing the Voight model or by using exponential relaxation kernels. However, mathematical models of problems of viscoelastic systems based on these assumptions cannot describe real processes occurring in shell constructions in early times [21, 23]. The choice of an exponential kernel in calculations was not incidental. The systems of integro-differential equations obtained from the calculations were, by way of differentiation, reduced to the solution of ordinary differential equations, which in most cases used to be solved by the well-known Runge-Kutta's numerical method. To the present day, none of the existing
methods has allowed one to solve such problems with weakly singular kernels of the type of Koltunov's, Rzhanitsyn's, Abel's, Rabotnov's, and others [21, 23].
Badalov, Eshmatov, and Yusupov [24] developed the numerical method based on the use of quadrature rules, which makes it possible to solve a system of nonlinear integro-differential equations with weakly singular kernels of the type of Koltunov-Rzhanitsyn's, Abel's, and Rabotnov's. Verlan, Abdikarimov, and Eshmatov [25] modified this method. The method provides reasonably high accuracy of results, is universal, enables one to solve a wide class of dynamic problems in viscoelasticity, and is economical from computer time consideration. Based on this method, many numerical results have been obtained by Eshmatov, Mukherjee [26], Mirsaidov, Abdikarimov, Khodzhaev [27], and other researchers [28-31] that are in a general agreement with evaluable experiments' data.
This work aims to study the influences of the viscoelastic and anisotropic properties of a material on the dynamical behavior of the plate. For the first time, an integral model is used to consider the viscoelastic properties of materials of an anisotropic structure. The integral mode correctly describes the rheological processes occurring in the studied structure during the entire time. The presented mathematical model makes it possible to obtain sufficiently accurate solutions that are well combined with experimental results.
2. Materials and Methods
The classical Kirchhoff-Love theory is used to construct a mathematical model of the dynamic behavior of a plate made of a material having anisotropic viscoelastic properties under the influence of various external loads (Fig. 1).
X
Figure 1. Anisotropic fiber-reinforced rectangular plate.
In this case, following [32-34], the stress and moment resultants have the form:
Nx = x + 4*2^ y + A1*6 Y xy + B1*1X x + B1*2X y + B1*6X xy Ny = Al* x + y + 4*6 Y xy + B1*2X x + B22X y + B*6X x, Nxy = A1*6S x + A26^ y + A66 Y xy + B1*6X x + B26X y + B6*6X x Mx = Bn£ x + B\2^ y + B1*6 Y xy + A*1X x + D1*2X y + D1*6X x My = B1*2S x + B22Z y + B26Y xy + D1*2X x + D22X y + D26X x Mxy = B\6^ x + B26t y + B66 Y xy + D1*6X x + D26X y + D6*6X; where A*■, B*, D*, i, j = 1, 2, 6 are the operators having the following form:
ф = I (Qijф)к (zk- zk-1), Bijф = 2 I (Q*ф)к (z2 - zL )
к=1
■к=1
D ф=|J; (Q ф)к (zk-zL),
QllФ=
Q12Ф =
Q11 cos4 0 + 2 (Q12 + 2Q66 ) sin2 20 + Q22 sin4 0
Q11 = VJE—, 1 (Q11 + Q22 - 4Q66 ) sin2 20 + Q12 f 1 -1 sin2 20
(1 -r*)
Ф,
(1 -r*)
Ф,
Q22 = E2
Q22Ф = Qi>= 026Ф =
1 - M-12^21
Q11 sin4 0 +1 (Q12 + 2Q66 ) sin2 20 + Q22 cos4 0
(1 -r* )ф, (1 -г* )ф,
Q11 sin3 0 cos 0 +1 (Q12 + 2Q66) sin 40 - Q22 sin 0 cos3 0 (1 - Г*) ф,
Q11 sin 0 cos3 0 -1 (Q12 + 2Q66) sin 40 - Q22 sin3 0 cos 0
£>66Ф_
4 (Q11 - 2Q12 + Q22 ) sin2 20 + Q66 cos2 20
(1 -г* )ф,
(2)
Q12 = = Q66 = G12, Г*Ф = jГ(t-т)ф(т)dт.
1 ^12^21 1 ^12^21 0
Here K is the number of plate layers, Eb E2 are the elastic moduli, G12 is the shear modulus, 212, 221 are the Poisson ratios, 0 is the angle characterizing the direction of the fibers relative to the
axis OX, Г* is the integral operator with the relaxation kernel Г(t).
The relationship between the strains in the median surface sx, Sy, yxy, Xx, Xy, Xxy and displacements u, v, w in directions x, y, z, written by considering the von Karman type of geometric nonlinearity, in the form [32-35]:
_du 1 f dw S x _ dx + 2 ^сХ У dy ' 2 Idy
dv 1 f dw ^ s y _ — + —
X x _
d 2 w
X y
V4T д 2 w
x dx2 ' y dy2
_du dv dw dw
, Y xy _ + + ~ ^ , Xxy _ -2"
xy dy dx dx dy d 2 w
(3)
dxdy
By substituting (1) and (3) into the equations of motion
dNx dN
xy
dx dy
+ Px = P
d 2u
dNy dNy
dt1
dx dy
+ Py =P
d 2v
d 2Mx dx2
d 2M,
dy2
d_ dy
- + 2-
d 2 M
xy
d
(
dxdy dx
Ar dw dw Nx — + Nxy—
Ji A. y \
dx dy
dt2
A
(4)
f
dw
dw
A
Nxy— + Ny —
dx dy
+q = p
d 2 w
y
dt2
the system of nonlinear integro-differential equations in partial derivatives that satisfies the boundary conditions of the problem (the edges are hinged) is obtained. The solution of this system is sought in the form:
M N
u ( x y, t)= Z I um
m=1 n=1
M N
v ( x y, t)= I I vm,
m=1 n=1
M N
w ( x ^ t )= II wm
m=1 n=1
, v mnx . nny
,( t ) cos-Sin——
w a b
nny
/ . mnx ( t ) sin-cos
a
/ \ . mnx . (t ) sin-sin
a
b
nny
(5)
where umn (t), vmn (t) and wmn (t), m, n = 1, 2, 3, ... are the unknown functions of time. Note that
eq.(5) gives just an approximation of the in-plane displacements u and v for the given boundary conditions since additional functions should be retained. Substituting the approximating functions (5) into the resulting system of equations and performing the procedure of the Bubnov-Galerkin method, a system of nonlinear ordinary integro-differential equations is obtained. Further, this system is integrated using the numerical method based on the use of quadrature formulas [24, 25].
3. Results and Discussion
The calculations use the simplest and, at the same time, quite common weakly-singular Koltunov-
Rzhanitsyn kernel having the form r(t) = Ae~^ta-1 (0 <a<l) as the relaxation kernel, where A, a,
P are the rheological viscosity parameters determined from the experiments [21, 23]. Experimental studies have shown that weak singular functions most accurately describe the rates of relaxation processes [21, 23]. One of the difficulties of using singular kernels in numerical calculations is their singularity at t = 0. A detailed description of the numerical method based on the use of quadrature formulas is given in [24, 25]. The method allows solving systems of nonlinear integro-differential equations, preliminarily transforming singular kernels into regular ones.
In the calculations, the following parameters for plastics: KAST-V ( Ei = 25.5 GPa, E2 = 14.91 GPa, G12 = 4.41 GPa, ^12 = 0.2, p = 1900 kg/m3) and EDF ( E1 = 24.5 GPa, E2 = 18.04 GPa, G12 = 4.1 GPa, p,12 = 0.18, p = 1830 kg/m3) have been used. To determine the
rheological parameters of the viscosity of these materials for different directions of reinforced fibers, the results of experiments given in [36] have been used. In this paper, based on the experiments, the rheological parameters of the viscosity for some glass-reinforced plastics, such as Textolite, KAST-V, SVAM, and EDF have been obtained. The parameters given in the paper correspond to the samples with directions of reinforced fibers of 0, 45, and 90 degrees relative to the OX axis. The values of the rheological viscosity parameters for samples with an arbitrary fiber orientation were obtained in [37] from experimental results by using interpolation formulas.
The problem of buckling of the anisotropic fiber-reinforced rectangular plate of thickness h with the sides a and b, subjected to dynamic compression along with one of the sides by force P(t) = Pot
(Po is the loading rate) was considered as follows. A problem with a similar mathematical formuation (elastic problem) was considered in [29].
Hereinafter, the following parameters of the plate are used in the calculations (unless otherwise specified): a = b = 0.5m, h = 0.5sm, 6 = 45°, P0 = 2MPa / s.
As a criterion determining the critical time, it is assumed that the sag of the deflection should not exceed a value equal to the thickness of the plate [35]. In shell structures, the greater the critical time, the more stable it is to dynamic loads.
The following graphs correspond to the results obtained for the midpoint of the simply supported plate. On the graphs, m (meter) is the dimension for the deflection, and s (second) is for time.
Fig. 2 shows the influence of the material's viscoelastic properties on the behavior of anisotropic reinforced plates made from KAST-V and EDF. In general, KAST-V and EDF are materials with different viscoelastic properties. It is evident from the figure that considering viscoelastic properties results in a decrease of the critical time. Especially, it is apparent for plates from KAST-V (A = 0.0208), which have more viscoelastic properties than EDF. The difference in critical times for plates from elastic and viscoelastic KAST-V reaches up to 21 %, whereas this difference is about 2.4 % for plates from EDF. This result shows that it is crucial to consider viscoelastic properties in solving this type of problem.
4 ^
; 2
1 fillip 7 /
0 0.72 1.44 2.16 2.88 3.6
Figure 2. Effect of viscoelastic properties on the behavior of anisotropic reinforced plate from KAST-V (1 - viscoelastic, 2 - elastic) and EDF (3 - viscoelastic, 4 - elastic).
Fig. 3 shows a graph of the changes in the deflections of the midpoints of the plates of various thicknesses. The results show that an increase in plate rigidity due to increased plate thickness leads to a proportional increase in the critical time value.
L-- 1 _J
0 0.6 1.2 1.8 2.4 3
Figure 3. Dependence of the deflection on time for various values of the thicknesses of the plate: 1 - h = 0.3 sm; 2 - h = 0.4 sm; 3 - h = 0.5 sm.
Fig. 4 shows the analogous results for X = 1; 1.2; 1.4. Here X is the ratio of plate edges. Therefore, if X = 1, then the plate is square-shaped. As it is clear from the graph, the increase in one of the edges of the plate increases the critical time.
u-.105r
Figure 4. Dependence of the deflection on time for various values of the parameter X(X = a/b) 1 - X = 1; 2 - X = 1.2; 3 - X = 1.4.
The various curves in Fig. 5 correspond to cases of changes in the deflections of the midpoint of a reinforced rectangular plate at different loading speeds. It should be noted here that in all cases, at the initial moments of time, the changes in the deflections are oscillations that are harmonic in shape, which begin to increase rapidly at certain points in time.
w.10®
3——_L7 i , i
0.6
1.2
1.8
2.4
Figure 5. Dependence of the deflection on time for various values of the velocities of loading 1 - Po = 2 MPals; 2 - Po = 2.5 MPals; 3 - Po = 3 MPals.
Fig. 6 shows the influence of changes in the direction of the fibers of the reinforced plate on the dynamic process. As the angle of direction of the fibers increases from 0 to 45 degrees, an increase in the critical time is observed. The difference between the critical time values for single-layer plates with fiber directions of 0 and 45 degrees is 12 %.
0 0.56 1.12 1.68 2.24 2.8
Figure 6. Dependence of the deflection on time for plates with different fiber orientations
1 - 0°; 2 - 15°; 3 - 30°; 4 - 45°.
As mentioned earlier, reinforced composites are a set (composition) of several reinforced layers, each with its own mechanical properties. Thus, by changing the composite structure, it is possible to create contractures, the behavior of which can be predicted in advance.
In this regard, the study of the behavior of laminated reinforced plates with different fiber orientations under the influence of axial compression is of particular interest. Fig. 7 shows the changes in the deflections of the midpoints of laminated reinforced plates made of KAST-V. Moreover, although all these plates have different fiber directions, their thickness is the same. The results show that for two-layer plates with fibers located at an angle of -45 degrees relative to the OX axis in one layer and 45 degrees in another, the critical time values are higher than the others. The layered fiber plate, parallel and perpendicular to the OX axis, has a lower critical time (i.e., it is less stable) than other plates with similar mechanical properties. The difference between the critical time values for the above two-layer plates is 13.3 %.
w ■ 10\--—,
0 0.56 1.12 1.68 2,24 2.8
Figure 7. Dependence of the deflection on time for sandwich plates with different fiber-orientated layers 1 - 0°/90°; 2 - 15°/-15°; 3 - 30°/-30°; 4 - 45°/-45°.
The results of studies of the behavior of reinforced plates for a wide range of changes in their mechanical, physical, and geometric parameters under dynamic compression of one of their sides are shown in Table 1.
Table 1. Critical time fiber-reinforced plates from
values for various geometric and physical parameters of anisotropic KAST-V.
Geometrical parameters of the
Physical parameters
Critical time values
Number
Fiber
a, m b, m h, sm q, Pa v, MPa/s of layers orientations Elastic problem Viscoelastic problem differ. (in %)
1 0.5 0.5 0.5 100 2 1 45° 3.2798 2.7117 21
2 0.6 0.5 0.5 100 2 1 45° 3.3358 2.7338 22
3 0.7 0.5 0.5 100 2 1 45° 3.5262 2.8604 23.3
4 0.5 0.5 0.4 100 2 1 45° 2.0238 1.6522 22.5
5 0.5 0.5 0.3 100 2 1 45° 0.9192 0.7680 19.7
6 0.5 0.5 0.5 200 2 1 45° 3.2110 2.5986 23.6
7 0.5 0.5 0.5 300 2 1 45° 3.1422 2.4982 25.8
8 0.5 0.5 0.5 100 2.5 1 45° 2.6268 2.1806 20.5
9 0.5 0.5 0.5 100 3 1 45° 2.1858 1.8248 19.8
10 0.5 0.5 0.5 100 2 1 0° 2.5984 2.3884 8.8
11 0.5 0.5 0.5 100 2 1 15° 2.7640 2.4038 15
12 0.5 0.5 0.5 100 2 1 30° 3.1046 2.6044 19.2
13 0.5 0.5 0.5 100 2 2 0°/90° 2.5984 2.3830 9
14 0.5 0.5 0.5 100 2 2 15°/-15° 2.7860 2.4152 15.4
15 0.5 0.5 0.5 100 2 2 30°/-30° 3.1396 2.6344 19.2
16 0.5 0.5 0.5 100 2 2 45°/-45° 3.3242 2.7496 20.9
17 0.5 0.5 0.5 100 2 3 45°/-45°/45° 3.2900 2.7168 21.1
The numerical results in the table show that taking into account the viscoelastic properties of structural materials leads to a significant decrease in the critical time value. In some cases, the differences in the results obtained by solving elastic and viscoelastic problems differ from each other by more than 20 %. It is also shown that the direction of the fibers in the structures also significantly affects the value of the critical time. For example, in two-layer plates, the fibers of which have directions at an angle of 45 and -45 degrees, the critical time values are greater than in a single-layer plate of the same thickness, the fibers of which are located at an angle of 45 degrees. This means that by changing the structure of the composite material, it is possible to create structures more resistant to dynamic loads.
4. Conclusions
The study of the dynamic stability of viscoelastic anisotropic reinforced plates, subjected to dynamic compression along with one of their sides, shows:
- It is important to consider the viscoelastic properties of the construction material. The results indicate that the difference in the critical time of elastic and viscoelastic problems is 21 % for plates made from KAST-V.
- The critical time values mainly depend on the direction of the reinforced fibers in each layer. In single-layer and double-layer plates, the difference in the critical time values depending on the direction of the reinforced fibers in places is 12 % and 13.3 %, respectively. An analysis of the results shows that the most resistant to these loads are double-layer plates with fibers located at an angle of -45 degrees relative to the OX axis in one layer and 45 degrees in another.
The results obtained in the work and the conclusions drawn on their basis allow us to accurately predict the dynamic behavior of plates made of reinforced plastics.
References
1. Ashok, R.B., Srinivasa, C.V, Basavaraju, B. Dynamic mechanical properties of natural fiber composites - a review. Advanced Composites and Hybrid Materials. 2019. 2 (4). Pp. 586-607. DOI: 10.1007/s42114-019-00121-8
2. Takagi, H. Review of Functional Properties of Natural Fiber-Reinforced Polymer Composites: Thermal Insulation, Biodegradation and Vibration Damping Properties. Advanced Composite Materials. 2019. 28 (5). Pp. 525-543. DOI: 10.1080/09243046.2019.1617093
3. Zenkour, A.M., Allam, M.N.M., Sobhy, M. Effect of Tranverse Normal and Shear Deformation on a Fiver-Reinforced Viscoelastic Beam Resting on Two-Parameter Elastic Foundations. International Journal of Applied Mechanics. 2010. 02 (01). Pp. 87-115. DOI: 10.1142/S1758825110000482
4. Yas, M.H., Aragh, B.S. Free vibration analysis of continuous grading fiber reinforced plates on elastic foundation. International Journal of Engineering Science. 2010. 48 (12). Pp. 1881-1895. DOI: https://doi.org/10.1016/j.ijengsci.2010.06.015
5. Allam, M.N.M., Zenkour, A.M., El-Mekawy, H.F. Bending response of inhomogeneous fiber-reinforced viscoelastic sandwich plates. Acta Mechanica. 2010. 209 (3). Pp. 231-248. DOI: 10.1007/s00707-009-0157-4
6. Chen, T.-J., Chen, C.-S., Chen, C.-W. Dynamic response of fiber-reinforced composite plates. Mechanics of Composite Materials. 2011. 47 (5). Pp. 549-560. DOI: 10.1007/s11029-011-9233-7
7. Amabili, M. Nonlinear vibrations of angle-ply laminated circular cylindrical shells: Skewed modes. Composite Structures. 2012. 94(12). Pp. 3697-3709. DOI: https://doi.org/10.1016/j.compstruct.2012.05.019
8. Chamis, C.C., Murthy, P.L.N. Vibration of Fiber Reinforced Anisotropic Composite Plates With Nanofiber Reinforced Matrices11-2010.
9. Ghasemi, F.A., Paknejad, R., Fard, K.M. Effects of geometrical and material parameters on free vibration analysis of fiber metal laminated plates. Mechanics & Industry. 2013. 14 (4). Pp. 229-238.
10. Wang, J., Zhang, J., Li, W., Xu, Z. Performance analysis of the plastic square plate based on the fiber reinforced PA66. MATEC Web Conf. 2016. 77.
11. Zu, X., Wu, H., Lv, H., Zheng, Y., Li, H. An Amplitude- and Temperature-Dependent Vibration Model of Fiber-Reinforced Composite Thin Plates in a Thermal Environment. 13 (7) 2020.
12. Li, H., Xue, P., Guan, Z., Han, Q., Wen, B. A new nonlinear vibration model of fiber-reinforced composite thin plate with amplitude-dependent property. Nonlinear Dynamics. 2018. 94 (3). Pp. 2219-2241. DOI: 10.1007/s11071-018-4486-5
13. Xue, Z.C., Li, Q.H., Wang, J.F., Lan, Z.X. Vibration Analysis of Fiber Reinforced Composite Laminated Plates with Arbitrary Boundary Conditions. Key Engineering Materials. 2019. 818. Pp. 104-112. DOI: 10.4028/www.scientific.net/KEM.818.104
14. Li, H., Wu, H., Zhang, T., Wen, B., Guan, Z. A nonlinear dynamic model of fiber-reinforced composite thin plate with temperature dependence in thermal environment. Composites Part B: Engineering. 2019. 162. Pp. 206-218. DOI: https://doi.org/10.1016/j.compositesb.2018.10.070
15. Yüksel, Y.Z., Akbaç, §.D. Buckling Analysis of a Fiber Reinforced Laminated Composite Plate with Porosity. Journal of Computational Applied Mechanics. 2019. 50 (2). Pp. 375-380. DOI: 10.22059/jcamech.2019.291967.448
16. Othman, M.I.A., Abouelregal, A.E., Said, S.M. The effect of variable thermal conductivity on an infinite fiber-reinforced thick plate under initial stress. Journal of Mechanics of Materials and Structures. 2019. 14 (2). Pp. 277-293. DOI: 10.2140/jomms.2019.14.277
17. Li, H., Niu, Y., Li, Z., Xu, Z., Han, Q. Modeling of amplitude-dependent damping characteristics of fiber reinforced composite thin plate. Applied Mathematical Modelling. 2020. 80. Pp. 394-407. DOI: https://doi.org/10.1016/j.apm.2019.11.048
18. Amabili, M., Reddy, J.N. The nonlinear, third-order thickness and shear deformation theory for statics and dynamics of laminated composite shells. Composite Structures. 2020. 244. Pp. 112265. DOI: https://doi.org/10.1016/j.compstruct.2020.112265
19. Yadav, A., Amabili, M., Panda, S.K., Dey, T., Kumar, R. Nonlinear damped vibrations of three-phase CNT-FRC circular cylindrical shell. Composite Structures. 2021. 255. Pp. 112939. DOI: https://doi.org/10.1016/j.compstruct.2020.112939
20. Allahyari, E., Asgari, M. Nonlinear dynamic analysis of anisotropic fiber-reinforced dielectric elastomers: A mathematical approach. Journal of Intelligent Material Systems and Structures. 2021. 32 (18-19). Pp. 2300-2324. DOI: 10.1177/1045389X21995879
21. Rabotnov, Y.N. Elements of Hereditary Solid Mechanics. First Thus edition . Moscow, MIR Publishers, 1980. 387 p. ISBN: 9780714714578.
22. Bogdanovich, A.E. Non-linear Dynamic Problems for Composite Cylindrical Shells. New York, Elsevier, 1993. 312 p.
23. Koltunov, M.A. Choice of kernels in solving problems involving creep and relaxation. Polymer Mechanics. 1969. 2 (4).
Pp. 303-311. DOI: 10.1007/BF00860064. URL: http://link.springer.com/10.1007/BF00860064 (date of application: 15.04.2021).
24. Badalov, F.B., Éshmatov, Kh., Yusupov, M. On certain methods of solving systems of integrodifferential equations encountered in viscoelasticity problems. Journal of Applied Mathematics and Mechanics. 1987. 51. Pp. 683-686.
25. Verlan, A.F., Abdikarimov, R.A., Eshmatov, Kh. Numerical modeling of nonlinear problems of the dynamics of viscoelastic systems with variable rigidity. Electronic modeling. 2010. 32 (2). Pp. 3-14.
26. Eshmatov, B., Mukherjee, S. Nonlinear Vibrations of Viscoelastic Composite Cylindrical Panels. Journal of Vibration and Acoustics. 2007. 129 (3). Pp. 285-296. DOI: 10.1115/1.2730532
27. Mirsaidov, M.M., Abdikarimov, R.A., Khodzhaev, D.A. Dynamics of a viscoelastic plate carrying concentrated mass with account of physical nonlinearity of material. Part 1. mathematical model, solution method and computational algorithm. PNRPU Mechanics Bulletin. 2019. 2019 (2). Pp. 143-153. DOI: 10.15593/perm.mech/2019.2.11
28. Eshmatov, B.Kh. Nonlinear vibrations and dynamic stability of viscoelastic orthotropic rectangular plates. Journal of Sound and Vibration. 2007. 300 (3-5). Pp. 709-726. DOI: 10.1016/j.jsv.2006.08.024.
29. Eshmatov, B., Abdikarimov, R., Komilova, Kh., Safarbayeva, N. Dynamics stability of anisotropic fiber-reinforced plate. E3S Web of Conferences 264. 05016 (2021). DOI: 10.1051/e3sconf/202126405016
30. Khodzhaev, D.A., Abdikarimov, R.A., Mirsaidov, M.M. Dynamics of a physically nonlinear viscoelastic cylindrical shell with a concentrated mass. Magazine of Civil Engineering. 2019. 91 (7). Pp. 39-48. DOI: 10.18720/MCE.91.4
31. Abdikarimov, R., Amabili, M., Vatin, N.I., Khodzhaev, D. Dynamic Stability of Orthotropic Viscoelastic Rectangular Plate of an Arbitrarily Varying Thickness. 11 (13) 2021.
32. Ashton, J.E., Whitney, J.M. Theory of laminated plates. Lancaster, PA, Technomic Publishing Co., Inc., 1970. 158 p.
33. Jones, R.M. Mechanics of Composite Materials. New York, McGraw-Hill Book Co., 1975. 450 p.
34. Whitney, J.M. Structural Analysis Of Laminated Anisotropic Plates. United Kingdom, Routledge, 2018. 356 p.
35. Volmir, A.S. The nonlinear dynamics of plates and shells. USA, Ohio, Foreign Technology Division Wright-Patterson Air Force, 1974. 543 p.
36. Tyuneeva, I.M. Relaxation characteristics of glass-reinforced plastics. Polymer Mechanics. 1970. 6. Pp. 492-494.
37. Eshmatov, B.Kh. Nonlinear Oscillations of a Viscoelastic Anisotropic Reinforced Plate. Mechanics of Solids. 2018. 53 (5). Pp. 568-572. DOI: 10.3103/S0025654418080101
Information about authors:
Bakhtiyor Eshmatov, PhD
ORCID: https://orcid. ora/0000-0003-0198-6679 E-mail: [email protected]
Rustamkhan Abdikarimov, Doctor of Physical and Mathematical Sciences ORCID: https://orcid.ora/0000-0001-8114-1187 E-mail: [email protected]
Marco Amabili, PhD
ORCID: https://orcid.ora/0000-0001-9340-4474 E-mail: marco. amabili@mcaill. ca
Nikolai Vatin, Doctor of Technical Sciences ORCID: https://orcid. ora/0000-0002-1196-8004 E-mail: [email protected]
Received 02.10.2022. Approved after reviewina 28.11.2022. Accepted 28.11.2022.
