Next Article in Journal
Photoluminescence Properties of Cyan-Emitting Lu3GaxAl5−xO12: Ce3+ Garnet Phosphors Synthesized in Nonreducing Atmosphere and at Different Temperature for High Quality w-LEDs
Previous Article in Journal
Numerical Investigation of the Temperature Field Effect on the Mechanical Responses of Conventional and Cool Pavements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Selected Problems of Random Free Vibrations of Rectangular Thin Plates with Viscoelastic Dampers

by
Marcin Kamiński
1,*,
Agnieszka Lenartowicz
2,
Michał Guminiak
3 and
Maciej Przychodzki
3
1
Department of Structural Mechanics, Faculty of Civil Engineering, Architecture & Environmental Engineering, Lodz University of Technology, Al. Politechniki 6, 90-924 Lodz, Poland
2
Doctoral School, Poznan University of Technology, Piotrowo 3 Street, 60-965 Poznan, Poland
3
Department of Structural Mechanics, Institute of Structural Analysis, Poznan University of Technology, 60-965 Poznan, Poland
*
Author to whom correspondence should be addressed.
Materials 2022, 15(19), 6811; https://doi.org/10.3390/ma15196811
Submission received: 3 August 2022 / Revised: 5 September 2022 / Accepted: 27 September 2022 / Published: 30 September 2022

Abstract

:
The main motivation of this work was to present a semi-analytical extension of the correspondence principle in stochastic dynamics. It is demonstrated for the stochastic structural free vibrations of Kirchhoff–Love elastic, isotropic and rectangular plates supported by viscoelastic generalized Maxwell dampers. The ambient temperature of the plate affects the dampers only and is included in a mathematical model using the frequency–temperature correspondence principle. The free vibration problem of the plate–viscoelastic damper system is solved using the continuation method and also the Finite Element Method (FEM). The stochastic approach begins with an initial deterministic sensitivity analysis to detect the most influential parameters and numerical FEM recovery of the polynomial representation for lower eigenfrequencies versus these parameters. A final symbolic integration leads to the first four basic probabilistic characteristics, all delivered as functions of the input uncertainties.

1. Introduction

Probabilistic mechanics is a topic that has been extensively studied, e.g., in [1,2], and one of its fundamental numerical methods—the Stochastic Finite Element Method (SFEM)—was invented and has been applied in the context of thin rectangular plate bending [3]. Considering multiple geometric scale uncertainties, the bending analysis of thin plates may also be performed using the Wavelet-based Stochastic Finite Element Method [4]. The SFEM was efficiently applied for the stochastic dynamic response analysis of graphite–epoxy composite plates [5]; other studies in this area can be found in [6,7]. An interesting scientific and engineering problem, the modeling of propeller blades, is presented in [8]. The authors did not use the typical random approach; however, they did perform an analysis of the deviation histogram of machining errors that may be random in nature.
The probabilistic structural response for the free damped vibrations of thin elastic and isotropic plates resting on viscoelastic supports is considered in this work. The nonlinear eigenproblem is solved here to determine the natural frequencies of the plate–viscoelastic damper system, and its solution is obtained thanks to the iterative continuation method presented by Lewandowski et al. [9,10,11]. The authors considered different types of viscoelastic dampers based on the generalized Maxwell model of a damper. An experimental study considering a generalized Maxwell model for nonlinear viscoelastic dampers was comprehensively performed by Lu et al. [12]. A comprehensive overview of some other deterministic and stochastic methods used in dynamics can be found in [13]. Several uncorrelated Gaussian random design variables are considered in this study, with an initial sensitivity study leading to the selection of the most influential parameters. The Least-Squares Method enables the determination of the random polynomials of eigenfrequencies, whose further integration with the Gaussian kernel finally returns the probabilistic characteristics. It should be underlined that the FEM experiments were entirely programmed in the Octave environment, whereas sensitivity and probabilistic analyses were all implemented in the computer algebra system MAPLE. The most important novelty of this work is the common application of a probabilistic numerical apparatus for a solution of an eigenproblem obtained using the continuation method.

2. Eigenvibration Analysis Methodology

The main purpose of this analysis is to determine the first few natural frequencies of thin rectangular elastic and isotropic plates supported by viscoelastic dampers and also their first four probabilistic characteristics. The numerical analysis of this problem is based upon the Finite Element Method (FEM) with a regular discretization including 4-node-plate finite elements with linear approximation functions, as shown in Figure 1 and Figure 2 [14]. The deformation vector w i e of the i -th node within the finite element e can be written as
w i e = [ w i φ i x φ i y ] T = [ w i w i y w i x ] T ; i = 1 ,   2 ,   3 ,   4 .
The displacement field within the element e is expressed as a linear combination of shape functions N k e ( x , y ) :
w e ( x , y ) = N e w e ,
where w e = [ w 1 e ,   w 2 e ,   w 3 e ,   w 4 e ] T and N e = [ N 1 e ,   N 2 e ,   N 3 e ,   N 12 e ] . The element stiffness matrix Ke and the consistent mass matrix M e are defined in the traditional way. The stiffness matrix is derived here analytically, and the mass matrix is derived numerically using 16-point Gaussian quadrature. Further, it is known that the equation of motion of a structure with viscoelastic dampers can be written in the following form [9,10]:
M q ¨ ( t ) + C q ˙ ( t ) + K q ( t ) = f ( t ) .
where C denotes the global plate-damping matrix.
The application of the Laplace transform with zero initial conditions leads to the following transform of Equation (3):
( s 2 M + s C + K ) q ¯ ( s ) = f ¯ ( s ) ,
where q ¯ ( s ) is the -transform of q ( t ) , and f ¯ ( s ) can be expressed as
f ¯ ( s ) = r = 1 n d ( K r + G r ( s ) ) L r q ¯ ( s ) .
In this formula, n d is the total number of dampers attached to the plate at selected nodes of a finite element mesh, and L r is a global matrix indicating the location of the r -th damper in the plate. A viscoelastic damper is represented graphically in Figure 3, and it consists of m spring-dashpot elements and an additional spring element. Each of the Maxwell elements contains a viscous part with the constant c j and an elastic part with the constant k j , where j = 1 ,   2 , , m .
All quantities appearing in Equation (5) can be expressed as follows:
K r = k 0 ;   G r ( s ) = j = 1 m k j s ν j + s
where ν j = k j / c j is the quotient of the stiffness and damping coefficients of the j -th Maxwell element. Obviously, the stiffness and damping parameters, k j and c j , of the individual elements constituting the viscoelastic damper attached to the structure additionally depend upon the temperature. Let the damper parameters be known for a certain reference temperature T 0 and have the values k j and c j . The damper model constants at temperature T can be expressed using the frequency–temperature correspondence principle as [9,10]
k ¯ j = k j ; j = 0 ,   1 ,   2 , , m ,
c ¯ j = c j α T ; j = 1 ,   2 , , m .
The shift factor α T is a function of temperature T and can be expressed by the William–Landel–Ferry formula as [15,16]
log 10 α T = C 1 ( T T 0 ) C 2 + T T 0 ,
where C 1 and C 2 are experimental constants. After substituting Equation (5) into Equation (4), the -transform of Equation (3) of the motion of a plate with viscoelastic dampers takes the following form:
( s 2 M + s C + K + K d + G d ( s ) ) q ¯ ( s ) = 0 ,
where
K d = r = 1 n d K r L r ,
G d ( s ) = r = 1 n d G r ( s ) L r .
Equation (10) represents a nonlinear eigenproblem that is solved for the eigenvalue s and the eigenvector q ¯ ( s ) using the continuation method comprehensively described by Lewandowski in, e.g., [9,10]. Below, the foundations of this method are quoted.
In the case of Equation (10), the components containing the variable s in the first power are multiplied by the parameter κ [ 0 ; 1 ] . The equation can then be re-written as
h 1 ( q ¯ , s ) = D ( s ) q ¯ ( s ) = 0 ,
where
D ( s ) = s 2 M + κ s C + K + K d + κ G d ( s )
In order for the elements of the eigenvector q ¯ corresponding to the eigenvalue s to be determined unambiguously, an additional normalizing equation of the following form is introduced into the matrix in Equation (13):
h 2 ( q ¯ , s ) = 1 2 q ¯ ( s ) T D ( s ) s q ¯ ( s ) a = 0 ,
where a has a given value.
In the first step of the continuation method, in Equation (13), the parameter κ 1 = 0 is assumed, and the generalized eigenproblem is solved.
( s 2 M + K + K d ) q ¯ ( s ) = 0
This problem was solved in the Octave program using the built-in command ‘eig’, which allows both standard and generalized eigenproblems to be solved. As a result of solving this problem, the first approximations of eigenvalues s 1 ( 1 ) , s 2 ( 1 ) , , s 3 n ( 1 ) and eigenvectors q ¯ 1 ( 1 ) , q ¯ 2 ( 1 ) , , q ¯ 3 n ( 1 ) are obtained. On their basis, the parameter a j ( 1 ) = s j ( 1 ) ( q ¯ j ( 1 ) ) T M q ¯ j ( 1 ) , where j = 1 ,   2 , ,   3 n , is determined.
In the l-th step ( l = 2 ,   3 ,   4 , ), the increment κ l is assumed, and the Newton method is used to solve the system of Equation (13) with the additional Equation (15). For this purpose, the system of incremental equations of the Newton method is solved using κ l = κ l 1 + κ l , s j ( k 1 ) , q ¯ j ( k 1 ) and a j ( k 1 ) . This system of equations takes the following form:
{ h 1 q ¯ δ q ¯ + h 1 s δ s = h 1 , h 2 q ¯ δ q ¯ + h 2 s δ s = h 2 ,
where
h 1 q ¯ = D ( s ) = s 2 M + κ s C + K + K d + κ G d ( s ) ,
h 1 s = ( 2 s M + κ C + κ G d ( s ) s ) q ¯ ,
h 2 q ¯ = q ¯ T ( 2 s M + κ C + κ G d ( s ) s ) ,
h 2 s = 1 2 q ¯ T ( 2 M + κ 2 G d ( s ) s 2 )
The derivatives in Equations (18b)–(18d) are calculated as follows:
G d ( s ) s = r = 1 n d j = 1 m k j ν j ( ν j + s ) 2 L r ,
2 G d ( s ) s 2 = r = 1 n d j = 1 m 2 k j ν j ( ν j + s ) 3 L r .
The increments δ q ¯ and δ s are obtained from the system of Equation (17), and the following are calculated:
s j ( k ) = s j ( k 1 ) + δ s ,
q ¯ j ( k ) = q ¯ j ( k 1 ) + δ q ¯ ,
a j ( k ) = 1 2 ( q ¯ j ( k ) ) T D ( s ) s q ¯ j ( k ) .
Successive approximations of the j -th eigenvalue and the j -th eigenvector in the l -th step of the algorithm are calculated until the desired accuracies ε 1 and ε 2 , of the final results are achieved, that is, until the following inequalities are satisfied:
δ s < ε 1 | s j ( k ) | ,
δ q ¯ < ε 2 q ¯ j ( k ) .
The final values of s j ( k ) , q ¯ j ( k ) and a j ( k ) obtained in the l -th step are taken as starting values for step l + 1 and the new parameter κ l + 1 = κ l + κ l + 1 .
The procedure described above is carried out up to the value of the parameter κ = 1 , when the final eigenvalues and eigenvectors for the nonlinear eigenproblem in Equation (10) are obtained.
The obtained eigenvalues of the problem in Equation (10) are complex numbers of the form s j = μ j + i η j . On this basis, the j -th natural frequency ω j of the structure and the non-dimensional damping ratio γ j of the j -th mode of vibration are determined from the formulas:
ω j 2 = μ j 2 + η j 2 ;   γ j = μ j ω j .
The continuation method makes it possible to calculate the first few natural frequencies and the corresponding non-dimensional damping ratios without having to solve the entire nonlinear eigenproblem. Using proprietary algorithms written in the Octave programming language, the natural frequencies of a plate equipped with viscoelastic vibration dampers are determined. In the program, the user can independently select the positions of the selected number of dampers in the FEM mesh nodes. With the help of proprietary software, the matrices occurring in the nonlinear eigenproblem in Equation (10) are determined. Using the built-in command, the program solves the generalized eigenproblem in Equation (16). Then, the program uses the above-described iterative algorithm of the continuation method, which allows the unknown natural frequencies of a plate equipped with viscoelastic vibration dampers to be determined.

3. Sensitivity and Uncertainty Analyses

The so-called normalized sensitivity gradients are determined using the following standard definition [1]:
Δ ω i Δ v j ( ω i v j ) v ¯ j · v ¯ j ω ¯ i
where v j ¯ denotes the mean value of the given parameter v j . The following design parameters affecting the natural frequencies of the rectangular plate are checked at the initial stage: (i) geometric dimensions of the plate l x × l y × H , (ii) material constants of the plate E , ν p and ρ p , (iii) damper parameters k 0 , k 1 and c 1 , (iv) the ambient temperature of the plate, namely, T , and (v) the reference temperature T 0 . Further probabilistic analysis is carried out for two parameters exhibiting the highest positive and negative sensitivity coefficients. The eigenfrequencies ω i of the plate under consideration are all found via the polynomial basis
ω i = j = 1 n C i j v j
via the Least-Squares Method fittings made on the basis of several FEM solutions for varying values of the parameter v [1,2]. Statistical optimization of this basis order is employed through the common minimization of the fitting variance and the maximization of the correlation factor. Finally, the basic probabilistic characteristics, i.e., expected values, standard deviations, coefficients of variation, skewness and kurtosis, are computed. The following integral definitions are applied:
E [ ω i ] = + j = 1 n C i j v j p v ( x ) d x ,   σ ( ω i ) = { + ( j = 1 n C i j v j E [ ω i ] ) 2 p v ( x ) d x } 1 2
α ( ω i ) = | σ ( ω i ) E [ ω i ] | ,   β ( ω i ) = μ 3 ( ω i ) σ 3 ( ω i ) ,   κ ( ω i ) = μ 4 ( ω i ) σ 4 ( ω i )

4. Numerical Experiment

A square isotropic plate fixed on one edge was discretized using the 14 × 14 plate rectangular finite element mesh, whose material properties are E = 205 GPa, νp = 0.3 and ρp = 7850 kg/m3. The plate dimensions are equal to l x × l y × H = ( 2.0 × 2.0 × 0.01 )   m . Three viscoelastic dampers are attached in the middle and at both ends of the free edge of the plate (see Figure 4).
These dampers contain a single spring element and also a Maxwell element (Figure 5) with the following parameters at T 0 = 0.2   : k 0 = 108.56   N / m ;   k 1 = 19968.09   N / m ;   c 1 = 229.63   Ns / m .
The influence of temperature on the values of the above-mentioned parameters is taken into account by applying the frequency–temperature correspondence principle. In order to calculate the value of the shift function from Equation (9), the values of the constants C 1 = 19.5 and C 2 = 80.2 were adopted. The initial sensitivity analysis results were computed analytically using polynomial responses and are compared in Table 1 below. Quite expectedly, the two most influential parameters for the given plate are its thickness (the minimum value of the gradient) and the edge length (the maximum gradient).
So, these two parameters were further selected for stochastic analysis. They were treated as Gaussian variables having expected values equal to the mean values given above and a coefficient of variation belonging to the interval α ( v ) [ 0.00 ,   0.20 ] . Such a wide interval was assumed to check theoretical variations in all characteristics being computed, and it includes the statistical scattering of all possible measurement techniques. It should be noted that the time effort and computer power required for the semi-analytical probabilistic solutions were only a little bit greater than the deterministic origin. Polynomial approximation was used to describe the response function. The degree of the polynomial was assumed to be 4 or 5 depending on the necessary accuracy of matching the response function. Examples of polynomials obtained for a random plate side length and a random plate thickness are given in Chapter 5 by the relations in (28) and (29), respectively.
Figure 6 and Figure 7 show the graphs of the dependence of the expected value, variance, skewness and kurtosis on the coefficient of variance when the random variable is the plate side length and its thickness, respectively; each of the graphs in Figure 6a–d and Figure 7a–d shows these characteristics for the first five natural frequencies, ω 1 ω 5 . These probabilistic coefficients were all found from their integral definitions, and they can be treated as exact in the probabilistic context (no convergence studies are necessary).
It is seen in Figure 6 that Gaussian uncertainty in the plate length causes some nonlinear increases even for the expected values of the fundamental free vibrations. This is in contradiction to the case illustrated in Figure 7, where they are simply constant. Moreover, the plate length randomness greatly amplifies the uncertainty in this problem, because its output-to-input ratio equals almost 3; the plate thickness shows a direct interrelation between the input and output CoVs (Figure 7b). Interestingly, the largest statistical dispersion is always associated with the first eigenvalue. Finally, it is seen that both the skewness and kurtosis monotonously increase together with an additional increase in the input CoV in Figure 6. Hence, the positive non-symmetry and concentration about the expected values remarkably increase together with the input uncertainty level. One notices that the differences between probabilistic characteristics for various eigenfrequencies are rather small. Figure 7c,d show that higher-order probabilistic coefficients relevant to the plate thickness oscillate about (for the 1st and the 2nd) or are almost equal to 0 (for higher eigenfrequencies).

5. Comparative Analysis—Results Validation

In order to better illustrate the previously obtained results, calculations were performed using three probabilistic approaches: the Semi-Analytical Method (SAM), the Stochastic Perturbation Technique (SPT) and Monte Carlo simulations (MCSs).
These analyses were applied for a plate exhibiting Gaussian uncertainty in the first natural frequency, uniquely defined by its mean value and the specific range of its coefficient of variation, i.e., α ( b ) [ 0.00 ,   0.025 ] . The global response function for random plate side length lx = ly = l was obtained in the form of the following fourth-order polynomial:
ω 1 ( l ) = 367.570938927759 555.768532731197 · l + 336.646138548979 · l 2   92.8342487373826 · l 3 + 9.66595643939497 · l 4
The choice of the degree of the approximating polynomial to the random quantity was dictated by the sufficient accuracy of matching the response function. The number of trials for the Monte Carlo simulations was equal to 150,000. The expected values E ( l ) , coefficients of variation α ( l ) , skewness β ( l ) and kurtosis κ ( l ) of the first natural frequency are presented in turn in Figure 8. Considering the large variations in the resulting statistics, the expected values and coefficients of variation for the first natural frequency are presented in Table 2 and Table 3.
A very good convergence of the results of the SAM and SPT approaches can be observed. On the other hand, MCS results are in good agreement only for the coefficient of variations α ( l ) .
Another random parameter for which the validation of the calculation results was performed is the plate thickness. In this case, the global response function for a random plate thickness h was obtained in the form of the following fifth-order polynomial:
ω 1 ( h ) = 8.35555733333377 7824.96232261055 · h   + 2.41798719696965 · 10 6 · h 2 2.55138948717944 · 10 8 · h 3   + 1.18596617132865 · 10 10 · h 4 2.02512179487178 · 10 11 · h 5
Similar to the above, the expected values E ( h ) , coefficients of variation α ( h ) , skewness β ( h ) and kurtosis κ ( h ) of the first natural frequency are presented in turn in Figure 9. Additionally, the expected values and coefficients of variation for the first natural frequency are presented in Table 4 and Table 5.
Similar to the previous case, a very good convergence of the results of the SAM and SPT approaches can be observed. The MCS results are in quite good agreement only for kurtosis and α ( h ) values from 0.0 to about 0.1.
The Monte Carlo simulation technique is the most time-consuming in terms of numerical calculations when performed on a typical PC computer using the MAPLE v.21 computational package. The computation time ratio can be expressed simply by the quotient t SAM t MCS or t SPT t MCS , which ranges from 1/50 to 1/100, depending mainly on the number of MCS trials.

6. Conclusions

(1) The theoretical and computational studies presented in this work clearly show that the common application of the semi-analytical stochastic approach with continuation methods allows for the fast and accurate determination of the probabilistic coefficients of free vibrations. It is demonstrated that the output randomness in rectangular elastic thin isotropic plate eigenfrequencies is usually not larger than the input statistical scattering of their design parameters. The only exception is in the plate dimension statistics; however, successful measuring techniques are so accurate now that despite the huge sensitivity to this parameter, the realistic coefficient of variation is less than a few percent. The fact that the largest resulting statistical dispersion is obtained for the first eigenfrequency may be important in the reliability assessment of various dynamical systems. This is due to the fact that a limit function, whose probability serves as the basis of the reliability index, is introduced as the difference between the induced vibrations and the lowest eigenfrequency.
(2) A comparative analysis for the first eigenfrequency was performed using three probabilistic approaches: the Semi-Analytical Method (SAM), the Stochastic Perturbation Technique (SPT) and Monte Carlo simulations (MCSs). For the first two methods, results with very high accuracy were obtained. The Monte Carlo simulation showed convergence only for selected random moments—the coefficient of variation for the random plate side length and kurtosis for the random plate thickness.
(3) A continuation of this research can include stochastic extensions of the continuation method in the numerical analysis of forced vibrations, possibly with the use of non-Gaussian random design parameters, too. In the case of any mathematical difficulties with computer algebra integration, the iterative generalized stochastic perturbation technique implemented as the SFEM is recommended. Further uncertainty analysis may be alternatively completed by the application of probabilistic entropy or its relative version presented recently in the literature [17,18].

Author Contributions

Conceptualization, M.K. and M.G.; Data curation, A.L.; Formal analysis, M.K. and M.G.; Funding acquisition, M.K.; Investigation, A.L.; Methodology, M.K., M.G. and M.P.; Project administration, M.K.; Software, M.K., A.L. and M.P.; Validation, M.K. and M.G.; Visualization, A.L.; Writing—original draft, A.L.; Writing—review & editing, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This paper was written in the framework of the research grant OPUS no. 2021/41/B/ST8/02432, entitled “Probabilistic entropy in engineering computations”, and sponsored by The National Science Center in Poland.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kamiński, M.M. The Stochastic Perturbation Method for Computational Mechanics, 1st ed.; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2013; ISBN 978-0-470-77082-5. [Google Scholar]
  2. Kamiński, M. On iterative scheme in determination of the probabilistic moments of the structural response in the stochastic perturbation-based Finite Element Method. Int. J. Numer. Methods Eng. 2015, 104, 1038–1060. [Google Scholar] [CrossRef]
  3. Arregui-Mena, J.D.; Margetts, L.; Mummery, P.M. Practical Application of the Stochastic Finite Element Method. Arch. Computat. Method. Eng. 2014, 23, 171–190. [Google Scholar] [CrossRef]
  4. Han, J.-G.; Ren, W.-X.; Huang, Y. A wavelet-based stochastic finite element method of thin plate bending. Appl. Math. Model. 2007, 31, 181–193. [Google Scholar] [CrossRef]
  5. Chandra, S.; Sepahvand, K.; Matsagar, V.A.; Marburg, S. Stochastic dynamic analysis of composite plate with random temperature increment. Compos. Struct. 2019, 226, 111159. [Google Scholar] [CrossRef]
  6. Hoi, C.-K.; Noh, H.-C. Stochastic finite element analysis of plate structures by weighted integral method. Struct. Eng. Mech. 1996, 4, 703–715. [Google Scholar]
  7. Mestrovic, M. Stochastic Finite Element Analysis of Plates, Theories of Plates and Shells Critical Review and New Application. In Proceedings of the Euromech Colloqium 444, Bremen, Germany, 22–25 September 2002. [Google Scholar]
  8. Zhu, L.; Yan, B.; Wang, Y.; Dun, Y.; Ma, J.; Li, C. Inspection of blade profile and machining deviation analysis based on sample points optimization and NURBS knot insertion. Thin-Walled Struct. 2021, 162, 107540. [Google Scholar] [CrossRef]
  9. Lewandowski, R. Non-Linear Free Vibrations of Beams by the Finite Element and Continuation Methods. J. Sound Vib. 1994, 170, 577–593. [Google Scholar] [CrossRef]
  10. Lewandowski, R.; Pawlak, Z. Dynamic analysis of frames with viscoelastic dampers modeled by rheological models with fractional derivatives. J. Sound Vib. 2011, 330, 923–936. [Google Scholar] [CrossRef]
  11. Lewandowski, R.; Bartkowiak, A.; Maciejewski, H. Dynamic analysis of frames with viscoelastic dampers: A comparison of damper models. Struct. Eng. Mech. 2012, 41, 113–137. [Google Scholar] [CrossRef]
  12. Lu, L.-Y.; Lin, G.-L.; Shih, M.-H. An experimental study on a generalized Maxwell model for nonlinear viscoelastic dampers used in seismic isolation. Eng. Struct. 2012, 34, 111–123. [Google Scholar] [CrossRef]
  13. Clough, R.W.; Penzien, J. Dynamics of Structures; Computers & Structures, Inc.: Berkeley, CA, USA, 1995. [Google Scholar]
  14. Hughes, T.J.R. The Finite Element Method. Linear Static and Dynamic Finite Element Analysis; Prentice Hall: Hoboken, NJ, USA, 1987. [Google Scholar]
  15. García-Barruetabeña, J.; Cortés, F.; Abete, J.M.; Fernández, P.; Lamela, M.J.; Fernández-Canteli, A. Experimental characterization and modelization of the relaxation and complex moduli of a flexible adhesive. Mater. Des. 2011, 32, 2783–2796. [Google Scholar] [CrossRef]
  16. Lewandowski, R.; Przychodzki, M. Influence of temperature on dynamic properties of frames with viscoelastic dampers. J. Civ. Eng. Environ. Archit. 2016, 33, 431–438. (In Polish) [Google Scholar]
  17. Kamiński, M. On Shannon entropy computations in selected plasticity problems. Int. J. Numer. Methods Eng. 2021, 122, 5128–5143. [Google Scholar] [CrossRef]
  18. Bredow, R.; Kamiński, M. Structural safety of the steel hall under dynamic excitation using the relative probabilistic entropy concept. Materials 2022, 15, 3587. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Plate discretization with rectangular plate finite elements.
Figure 1. Plate discretization with rectangular plate finite elements.
Materials 15 06811 g001
Figure 2. Finite element type plQ4 used for discretization of the tested plate: (a) node numbering; (b) active degrees of freedom. The numbers 1–4 are the element node numbers.
Figure 2. Finite element type plQ4 used for discretization of the tested plate: (a) node numbering; (b) active degrees of freedom. The numbers 1–4 are the element node numbers.
Materials 15 06811 g002
Figure 3. Generalized Maxwell model of viscoelastic damper.
Figure 3. Generalized Maxwell model of viscoelastic damper.
Materials 15 06811 g003
Figure 4. Plate fixed on one edge with three viscoelastic dampers.
Figure 4. Plate fixed on one edge with three viscoelastic dampers.
Materials 15 06811 g004
Figure 5. Model of a viscoelastic damper attached to the tested plate.
Figure 5. Model of a viscoelastic damper attached to the tested plate.
Materials 15 06811 g005
Figure 6. Expected values (a), coefficients of variation (b), skewness (c) and kurtosis (d) for a random length of the plate.
Figure 6. Expected values (a), coefficients of variation (b), skewness (c) and kurtosis (d) for a random length of the plate.
Materials 15 06811 g006
Figure 7. Expected values (a), coefficients of variation (b), skewness (c) and kurtosis (d) for a random thickness of the plate.
Figure 7. Expected values (a), coefficients of variation (b), skewness (c) and kurtosis (d) for a random thickness of the plate.
Materials 15 06811 g007
Figure 8. Expected values (a), coefficients of variation (b), skewness (c) and kurtosis (d) for a random plate side length.
Figure 8. Expected values (a), coefficients of variation (b), skewness (c) and kurtosis (d) for a random plate side length.
Materials 15 06811 g008aMaterials 15 06811 g008b
Figure 9. Expected values (a), coefficients of variation (b), skewness (c) and kurtosis (d) for a random plate thickness.
Figure 9. Expected values (a), coefficients of variation (b), skewness (c) and kurtosis (d) for a random plate thickness.
Materials 15 06811 g009
Table 1. Values of the normalized sensitivity gradient for successive design variables on which the value of the first natural frequency of the tested plate depends.
Table 1. Values of the normalized sensitivity gradient for successive design variables on which the value of the first natural frequency of the tested plate depends.
Design Variable
Name
Design VariableNormalized Sensitivity
Gradient
Plate length l x = l y [m]–1.908109
Thickness of the plate H [m]1.136478
Young’s modulus E [N/m2]0.473791
Poisson ratio ν p [–]0.0669623
Density ρ p [kg/m3]–0.584767
Stiffness parameter of the damper’s Kelvin element k 0 [N/m]0.0109382
Stiffness parameter of the damper’s Maxwell element k 1 [N/m]–0.0708819
Viscosity parameter of the damper’s Maxwell element c 1 [Ns/m]0.164465
Reference temperature T 0 [°C]0.0189089
Ambient temperature T [°C]–0.0295390
Table 2. Expected values of the first circular frequency for a random plate side length.
Table 2. Expected values of the first circular frequency for a random plate side length.
α ( b ) E ( l )
SAMSPTMCS
0.025362.046846366682362.046846366682362.047466435600
0.050362.046908966730362.046908966730362.046431550698
0.075362.047013300146362.047013300146362.046932076133
0.100362.047159366940362.047159366940362.046112284009
0.125362.047347167119362.047347167119362.046692036114
0.150362.047576700696362.047576700696362.049805913200
0.175362.047847967686362.047847967686362.048943295710
0.200362.048160968106362.048160968106362.050809955170
0.225362.048515701978362.048515701978362.048477124862
0.250362.048912169324362.048912169324362.051187451775
Table 3. Coefficients of variation of the first circular frequency for a random plate side length.
Table 3. Coefficients of variation of the first circular frequency for a random plate side length.
α ( b ) α ( l )
SAMSPTMCS
0.0250.000379138402320590.000379138402257680.00037820668516146
0.0500.000758276797916070.000758276797790300.00075890901680434
0.0750.001137415180061430.001137415179872770.00113579512680720
0.1000.001516553542031540.001516553541780010.00151584644876523
0.1250.001895691877101270.001895691876786860.00189866159120374
0.1500.002274830178545410.002274830178168110.00227385578919360
0.1750.002653968439638680.002653968439198510.00265356938565002
0.2000.003033106653655750.003033106653152690.00303502740231922
0.2250.003412244813871140.003412244813305230.00340377236370069
0.2500.003791382913559340.003791382912930500.00378941244520895
Table 4. Expected values of the first circular frequency for random plate thickness.
Table 4. Expected values of the first circular frequency for random plate thickness.
α ( b ) E ( h )
SAMSPTMCS
0.02515.102030393868015.102030393867955.8307084247803
0.05015.075052826638115.075052826635956.6189485257520
0.07515.030903051808615.0309030517976–182.896206992036
0.10014.970800325210014.970800325210059.6193578514476
0.12514.896451605004914.896451604919461.9339027954595
0.15014.810051551688114.810051551510964.6258515558037
0.17514.714282528086814.714282527758568.1357033688504
0.20014.612314599360014.612314598800072.1942549445554
0.22514.507805532999314.507805532102276.9222494795175
0.25014.404900798828114.404900797460982.4482857143314
Table 5. Coefficients of variation of the first circular frequency for random plate thickness.
Table 5. Coefficients of variation of the first circular frequency for random plate thickness.
α ( b ) α ( h )
SAMSPTMCS
0.0250.0216902067352790.0216902066460100.09677329520225
0.0500.0438605697533050.0438605695711280.19522123484153
0.0750.0669517571808770.066951756898434–0.10707504128592
0.1000.0913262534885260.0913262530945810.39417790477171
0.1250.1172313931798190.1172313926590760.49988129568701
0.1500.1447656158060370.1447656151389530.60551075944427
0.1750.1738508453200130.1738508444827070.71407533547204
0.2000.2042167136859340.2042167126503110.82487834596053
0.2250.2354072860784690.2354072848128280.93330578277003
0.2500.2668290370727050.2668290355432921.05206621395503
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kamiński, M.; Lenartowicz, A.; Guminiak, M.; Przychodzki, M. Selected Problems of Random Free Vibrations of Rectangular Thin Plates with Viscoelastic Dampers. Materials 2022, 15, 6811. https://doi.org/10.3390/ma15196811

AMA Style

Kamiński M, Lenartowicz A, Guminiak M, Przychodzki M. Selected Problems of Random Free Vibrations of Rectangular Thin Plates with Viscoelastic Dampers. Materials. 2022; 15(19):6811. https://doi.org/10.3390/ma15196811

Chicago/Turabian Style

Kamiński, Marcin, Agnieszka Lenartowicz, Michał Guminiak, and Maciej Przychodzki. 2022. "Selected Problems of Random Free Vibrations of Rectangular Thin Plates with Viscoelastic Dampers" Materials 15, no. 19: 6811. https://doi.org/10.3390/ma15196811

APA Style

Kamiński, M., Lenartowicz, A., Guminiak, M., & Przychodzki, M. (2022). Selected Problems of Random Free Vibrations of Rectangular Thin Plates with Viscoelastic Dampers. Materials, 15(19), 6811. https://doi.org/10.3390/ma15196811

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop