Shock and Vibration

Shock and Vibration / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 1237674 | 20 pages | https://doi.org/10.1155/2019/1237674

A Semianalytical Three-Dimensional Elasticity Solution for Vibrations of Orthotropic Plates with Arbitrary Boundary Conditions

Academic Editor: Nerio Tullini
Received14 Apr 2019
Accepted09 Jun 2019
Published25 Jul 2019

Abstract

A semianalytical three-dimensional (3D) elasticity solution for the vibration of the orthotropic plate is presented under arbitrary boundary conditions. Three-dimensional (3D) elasticity theory provides the theoretical support for the energy function of orthotropic plates. The orthotropic plates which have the arbitrary boundary condition are realized by the way of arranging three sets of linear springs at the edges. With the aim of eliminating the nonsmooth phenomenon at the edges, the admissible displacement function of an orthotropic plate is expressed with a modified Fourier series solution. Under this framework, a change that occurs on the boundary conditions only needs to modify the boundary parameters of the orthotropic plate, without the need for new derivation, thus greatly saving the modeling time. The convergence and accuracy of the proposed method are better than those of the published literature. Lastly, the new vibration results and parametric research of thick orthotropic plates as well as the geometric parameter are also presented.

1. Introduction

As a basic building block, the orthotropic plates have been recognized and widely used by experts in various fields such as aerospace, architecture, shipbuilding, and other industries because of their good orthotropic mechanical properties such as strength to weight and high ratio of stiffness. In the above engineering applications, orthotropic rectangular plates as foundation components usually receive unknown loads. Thus, an intensive analysis about the vibration characteristics of the orthotropic rectangular plates has a beneficial effect on structural designers to carry out preliminary engineering design work. Compared with traditional isotropic materials, the vibration characteristics of orthotropic materials may be more complicated because of the material anisotropy. Traditional analytical and semianalytical numerical methods were used to solve the sandwich plate, laminated composite, and multilayered plates or functionally graded material plates with simply supported or clamped boundary conditions [13]. Therefore, it is very urgent and necessary to develop a method to handle the orthotropic plates under the conditions of arbitrary boundary.

Up to now, characteristic of vibration about the orthotropic rectangular plate has attracted much attention and has achieved fruitful results. According to two variable refined plate theories, Thai and Kim [46] proposed a Levy-type solution for free vibration and bucking analysis of orthotropic plates having two opposite edges simply supported and the other two edges having arbitrary classical boundary conditions simultaneously. Sepahvand et al. [7] concluded a new plate theory about a study of the orthotropic plate under the condition of stochastic-free vibration with a method of generalized polynomial chaotic expansion. And an analytical homogenization model for torsion of orthotropic sandwich panels was proposed by Abbès and Guo [8]. Ghugal and Sayyad [9] proposed a concept of triangular shear deformation of the thick orthotropic plate to study the free vibration of this plate under classical boundary conditions. The vibration mechanism of orthotropic plates on nonlinear elastic foundations is given by Hsu [10] through the differential quadrature method. Bahmyari and Rahbar-Ranji [11] researched the thickness variation of a rectangular orthotropic plate on the nonuniform elastic foundation in one or two directions in the way of applying the meshless Galerkin method. In view of the classical plate theory (CPT), Valizadeh et al. [12] developed the isogeometric analysis to research buckling and free and forced vibration analyses to research orthotropic plates under the conditions of classical boundary. On the basis separation of the variables method, Xing and Xu [13] calculated the exact solution when the rectangular thin plates with orthotropic properties are free to vibrate and apply them to various classical boundary conditions. It was precise because of the rational use of the mixed first-order shear deformation plate theory that Zenkour [14] studied on the bending response about the simply supported orthotropic plate. Bakhadda and his coworkers [1518] studied the free vibration, bending response, dynamic response of composite plates, and FGM plates by various valuable shear deformation theories. Papkov and Banerjee [19] proposed a new solution, which aims to solve the free vibration and buckling analysis of orthotropic rectangular plates under the different classical boundary conditions. And the characteristics of the orthotropic rectangular plate under free vibration are studied by Liu et al. [20]; in particular, the plates undergo free vibration in one or two directions. Through the Rayleigh–Ritz procedure, Shi et al. [21] studied a special plane vibration analysis and the special point is that the elastic boundary constraints of the orthotropic rectangular plates are nonuniform as well as internally supported. An effective hybrid method proposed by Jafari and Eftekhari [22] is suitable for the free vibration and buckling analysis about the orthotropic rectangular plate under various classical boundary conditions. Liu et al. [23, 24] proposed an accurate closed-form solution regarding orthotropic rectangular Mindlin plates to solve the free vibration problem under some classical boundary conditions. Chang and Tarn [25] researched the 3D analysis associated with orthotropic rectangular plates which supported external loads from the top and bottom surfaces by using the state space approach. In the aspect of 3D analysis, Teo and Liew [26] found that a method called differential quadrature (DQ) can effectively analyze rectangular orthotropic plates which have the boundary conditions of SSSS, SCCC, CSSS, and CCCC. Tian et al. [27, 28] proposed a new higher order analysis model for sandwich and FGM plates with flexible core and obtained the fundamental frequencies and static bending under sinusoidal and uniformly distributed loads. A semianalytical method is employed to analyze free vibration behaviors of composite laminated doubly curved, cylindrical, and spherical shells by Li et al. [2931]. Abdelaziz and his coworkers [3234] analyzed free vibration, bending, buckling, and postbuckling of FGM sandwich plates and nanobeams with various boundary conditions by various shear deformation theories. Cui et al. [35] investigated energy absorption of double layered and Sandwich plates under ballistic impact via experiment and numerical simulation. Liu et al. [36] assessed the material strain rate effects on steel plates under lateral impact loads.

The above results have been clearly described. It does not have trouble in finding that most of the current research results are confined to thin orthotropic plates and moderately thick orthotropic plates. In their research, they are based on the two-dimensional (2D) elasticity theory. From existing literature and materials, we can see that two-dimensional elasticity theories are based on assumptions. However, due to the existence of the assumption, the computation results are not credible at high modes. Based on this background, the three-dimensional theory is developed. And its results are generally accepted by researchers and used as contrastive data for reference in two-dimensional elasticity theory. However, for the greatest knowledge of the author, only two papers are currently committed to research the 3D vibration analysis of thick orthotropic plates that the conditions are constrained within the classical boundary.

A 3D vibration analysis of thick orthotropic plates under the general boundary condition is firstly performed in this paper. The theoretical formula is obtained by using this three-dimensional elastic theory. The arbitrary boundary conditions of the orthotropic plates are achieved through arranging three sets of linear springs at the edges. No matter how the boundary conditions are, the admissible displacement function of a thick orthotropic plate is uniformly expressed with improved Fourier series. Under this framework, when the boundary conditions of the orthotropic plates change, only the boundary parameters need to be modified, without the need for new derivation, thus greatly saving the modeling time. On this foundation, the astringency and exactitude of existing methods were verified by contrasting the results of existing literature and data. Not limited to this, some new numerical results are further displayed as reference data for future data validation.

2. Theoretical Formulations

2.1. Description of the Thick Orthotropic Plates

Figure 1 shows a thick orthotropic plate under the conditions of general boundary. The coordinate system (x, y, z) of the thick orthotropic plate is located at the bottom. a, b, and h represent length, width, and height successively. The symbols u, , and , respectively, represent displacement components in the length, width, and thickness directions. As shows in Figure 1(b), three kinds of linear springs, i.e., linear elastic ku in the direction of length, linear elastic in the direction of width, and linear elastic in the direction of thickness, are uniformly arranged at the structural boundaries. By simply setting the stiffness coefficient, the real boundary force can be simulated and the equivalent ideal boundary condition is available. Each boundary spring stiffness values are set to infinity (ku =  =  = 1016); for instance, it becomes very easy to obtain clamped boundary conditions.

2.2. Energy Expression of Thick Orthotropic Plates

On account of the 3D elastic theory, the linear strain-displacement regulations of thick orthotropic plates are given as follows:

The stress of a thick orthotropic plates constrained by linear elastic theory can be defined as follows:where stiffness coefficients Cij are obtained as below:

It must be emphasized here that when the above material parameters are simplified to isotropic materials, a thick isotropic plate can be obtained directly, and their stiffness coefficients can be simplified as follows:

Next, the U representing the strain energy of the thick orthotropic plates which is based on the elastic theory is as follows:

The T representing kinetic energy [37] about thick orthotropic plates is as follows:

Boundary conditions are efficiently simulated by artificial spring boundary technology. The actual operation is based on the introduction of three types of springs at the boundary to correspond to the actual boundary forces. Therefore, Usp representing the potential energy which is located in the boundary spring can be expressed as follows:

The Lagrangian function of the thick orthotropic plates obtained from the above expression is expressed as follows:

2.3. Admissible Displacement Functions

The main task of this chapter is to establish admissible displacement function of the thick plate displacement when using the Ritz method. According to the previous research, it is easy to find that the traditional displacement function has some shortcomings. Therefore, in this paper, we will make use of the modified Fourier series proposed by the author before. The principle and mathematical analysis process of the modified Fourier series solution are included in [3841]. Therefore, on the basis of the improved series solution principle and the motion equation, the displacement tolerance function of the thick orthotropic plate is uniformly described as a periodic function:where UΩ, VΩ, and WΩ are representative of the internal displacement distribution function about the thick orthotropic plate and , , and represent the complementary sequences of the displacement on the boundaries of the thick orthotropic plates. And the three 3D Fourier coefficients vector are Amnq, Bmnq, and Cmnq, respectively. The following formulas represent the specific elements of each vector symbol:

2.4. Solution Procedure

With the aim of extending Ritz method to the partial differential operation of the unknown coefficients about above displacement functions, the kinetic energy as described below can be obtained by the way of replacing equations (9a)–(9c) with the kinetic energy (T) function:

In addition, by substituting equations (9a)–(9c) into an expression about strain energy, the following kinetic energy can be obtained:

Substituting equation (7) for equations (9a)–(9c) under the same processing form can be obtained:

The Lagrangian equation [31] L of the vibration of the thick orthotropic plate structure can be obtained in the light of the above expression, and the equation is expressed in the following manner:

In addition, equations (11)–(13) are substituted into equation (14), and then the Lagrangian equations are partially deduced to zero according to the sealed Fourier coefficient vector, and its expression is as follows:

The matrix form can also be used to describe the equation:where K and M indicate the stiffness matrices and mass matrices of the FGP annular and sector plate. And X is the sealed Fourier coefficient vector that is described as . The two important parameters natural frequency and associated mode shape about thick orthotropic plates are acquired by the generalized eigenvalues of equation (16).

3. Numerical Results and Discussion

A series of numerical studies are carried out in this section, mainly including the following contents: (a) firstly, the convergence of this method is studied, the correctness is verified, and the definition of boundary conditions is studied; (b) secondly, accuracy of the method is validated by comparing with the relevant literature and datum; and (c) on this basis, the vibration performance of thick orthotropic plates which have arbitrary boundary conditions is further researched. In identifying the boundary conditions, these three symbols F, C, and S represent three boundary conditions, namely, free, clamping, and simple supported boundary conditions. And all the calculation codes are based on the MATLAB 8.2 platform. In addition, if no special declaration is made, the frequency dimensionless formula is defined as .

3.1. Convergence and Stiffness Value Study

As we expected, as an excellent numerical method, it depends on its convergence and numerical constancy. However, it is found that the convergence and numerical constancy of proposed solution depend on the truncation number of admissible displacement functions (equations (10a)–(10k)). Therefore, it is important to research the convergence property and numerical constancy about these admissible displacements functions (equations (10a)–(10k)). The convergence of frequency parameters Ω about the orthotropic plate, which own different censoring numbers, is plotted in Table 1. Three kinds of thickness, i.e., h/a = 0.1 (thin plate), h/a = 0.3 (moderately thick plate), and h/a = 0.5 (thick plate), and seven types of truncated number, i.e., MNQ = 4  4  4, 5  5  5, 6  6  6, 7  7  7, 8  8  8, 9  9  9, and 10  10  10, are examined in this table. The material property of this table is the aragonite crystal, and the specific material information is as follows: C22/C11 = 0.543103, C12/C11 = 0.23319, C33/C11 = 0.530172, C13/C11 = 0.010776, C23/C11 = 0.098276, C66/C11 = 0.262931, C55/C11 = 0.159914, and C44/C11 = 0.26681. The information in the table is available that no matter how thick the plate is, its natural frequency tends to stabilize rapidly when the number of truncation MNQ increases. Especially when the truncated number value is more than 7  7  7, the calculation result can keep the 4 digit effective number unchanged. Therefore, at this point, we can think that the calculation results have stabilized steadily. And the MNQ representing truncated values of the revised sequence under the continuation of numerical examples is defined as MNQ = 8  8  8.


MNQMode number
12345678

4  4  40.047640.103570.119000.161080.161080.169670.190380.21697
5  5  50.047450.103330.118830.161080.161080.169440.189500.21697
6  6  60.047450.103330.118830.161080.161080.169440.189040.21697
7  7  70.047420.103300.118800.161080.161080.169420.188940.21697
8  8  80.047420.103300.118800.161080.161080.169420.188870.21697
9  9  90.047420.103290.118800.161080.161080.169410.188850.21697
10  10  100.047420.103290.118800.161080.161080.169410.188830.21697
DQM0.047420.103290.118800.161090.161090.169410.188760.21697

4  4  40.332310.483230.483230.650210.650430.656650.898120.96645
5  5  50.332030.483230.483230.649800.650430.656180.897500.96645
6  6  60.332030.483230.483230.649790.650430.656180.897500.96645
7  7  70.332000.483230.483230.649740.650430.656110.897410.96645
8  8  80.332000.483230.483230.649740.650430.656110.897410.96645
9  9  90.331990.483230.483230.649730.650430.656100.897390.96645
10  10  100.331990.483230.483230.649730.650430.656100.897390.96645
DQM0.332000.483270.483270.649740.650430.656120.897420.96654

4  4  40.703910.805380.805381.082351.243941.305231.492271.61076
5  5  50.703440.805380.805381.082351.242591.304461.492271.61076
6  6  60.703430.805380.805381.082351.242571.304441.492271.61076
7  7  70.703370.805380.805381.082341.242371.304341.492271.61076
8  8  80.703370.805380.805381.082341.242371.304331.492271.61076
9  9  90.703360.805380.805381.082341.242321.304311.492271.61076
10  10  100.703360.805380.805381.082341.242321.304311.492271.61076
DQM0.703380.805450.805451.082401.242301.304301.492301.61090

In this paper, we aim to establish an analytical model for vibration of orthotropic plates under the conditions of the arbitrary boundary by the method of arranging three kinds of linear springs at the edges. In addition, the simulation of the clamped boundary conditions is to set the spring stiffness value to infinity, and in practical cases, numerical instability will occur when the stiffness value is set to infinity due to the limitation of computational resources. Therefore, we must use a finite value to replace the infinite stiffness value. Thus, the vibration performances on orthotropic plates are affected with the elastic boundary.

Frequency parameters Ω of the elastically constrained orthotropic plate are summarized in Tables 24. Three elastic boundary conditions are defined in the following ways: (1) elastically restrained in the x direction (Table 2), ; (2) elastically restrained in the y direction (Table 3), ; and (3) elastically restrained in the z direction (Table 4), . To simplify the research, the y = constant boundary is a free boundary, one end of orthotropic plates (x = a) is a clamped boundary condition, and the other edge x = 0 is elastically constrained. Through the comparison of three tabular data, the results show that the parameters have the greatest effect with the vibration behavior of orthotropic plates, followed by and finally the parameter. This phenomenon can also be found directly from the modal shapes in Figures 24. Based on the above research, we further define three conditions with the classical boundary and four conditions of the elastic boundary in the following example, and the detailed spring stiffness coefficient is shown in Table 5.


kuMode number
12345678

1020.187020.291620.294970.621040.744590.786040.809410.84366
1040.187020.291620.294970.621040.744590.786040.809410.84366
1060.187020.291620.294970.621050.744590.786040.809410.84366
1080.187050.291660.294970.621120.744860.786200.809420.84369
10100.190010.295130.295190.628160.771150.800770.809600.84620
10120.227820.298150.333920.753640.817250.887581.029901.09273
10140.233660.298930.338990.776240.823080.896421.061011.09704
10160.233720.298940.339040.776480.823170.896531.061351.09709
10180.233720.298940.339040.776480.823180.896531.061351.09709


kuMode number
12345678

1020.187020.291620.294970.621040.744590.786040.809410.84366
1040.187020.291620.294970.621040.744590.786040.809410.84366
1060.187020.291630.294970.621040.744590.786040.809420.84366
1080.187020.292590.295100.621050.744590.786160.809510.84367
10100.187020.306980.369850.621370.744700.798400.817900.84449
10120.187050.438160.623920.659140.745180.856570.908491.07874
10140.187060.459800.624220.667170.745210.861110.922341.10855
10160.187060.460100.624230.667250.745210.861210.922541.10886
10180.187060.460100.624230.667250.745210.861210.922541.10886


kuMode number
12345678

1020.187020.291620.294970.621040.744590.786040.809410.84366
1040.187020.291620.294970.621040.744590.786040.809410.84366
1060.187030.291620.294980.621050.744590.786040.809420.84366
1080.188750.291620.295620.621380.744590.786040.809590.84397
10100.291620.299120.347120.656550.744590.786040.825680.87358
10120.291630.481090.534660.744610.786050.942371.052151.06005
10140.291630.484540.543690.744610.786060.953551.052931.07903
10160.291630.484580.543790.744610.786060.953681.052941.07922
10180.291630.484580.543790.744610.786060.953681.052941.07922


Boundary conditionSpring stiffness
ku

Free edge (F)000
Simply supported edge (S)010161016
Clamped edge (C)101610161016
Elastic constraint 1 (E1)101010161016
Elastic constraint 2 (E2)101610101016
Elastic constraint 3 (E3)101610161010
Elastic constraint 4 (E4)101010101010

3.2. Vibration Analysis of Thick Orthotropic Plates with Various Boundary Conditions

As mentioned before, there are relatively few studies on this topic. Therefore, in order to dig deeper into the details of the research content, the author will further carry out the study of vibration analysis with these thick orthotropic plates of various boundary conditions. The numerical results of a thick plate with different boundary conditions compared with published literature [26] are shown in Table 6. It can be easy to find that the maximum error between the present method and literature data is no more than 0.15%. So, it can be concluded that the present method is accurate to solve the vibration behaviors of thick plates with various boundary conditions.


ModeCCCCSCCCCSSS
0.30.50.30.50.30.5
Reference [23]PresentReference [23]PresentReference [23]PresentReference [23]PresentReference [23]PresentReference [23]Present

10.442530.442820.852760.853330.420550.420810.827400.827850.351290.351370.715420.71557
20.713150.713731.302501.303650.660430.660371.100801.100720.483270.483230.805450.80538
30.749960.750681.361501.361480.698530.698871.292501.293180.580910.580880.967750.96771
40.816760.816781.393301.394210.745410.746091.297101.297120.652290.652341.252901.25337
50.965680.966591.704001.703790.778160.778161.391001.391930.671450.671641.306701.30682
61.022301.022141.760401.761900.959120.959851.675301.675260.841990.841901.400001.39988
71.046601.047681.793501.793481.035401.036051.755101.756230.902300.902451.610901.61076
81.076601.076571.861901.864301.056201.056101.760001.759890.966540.966451.636201.63612

The free vibration experiments of the isotropic plate structure was carried out to verify the feasibility of the present method, and the material properties are chosen as E = 210 GPa, ρ = 7850 kg/m3, μ = 0.3, a = b = 1.2 m, and h = 0.008 m. The test was carried out by loading on multiple excitation points; meanwhile, the responding signals were collected from a single point. The comparison of natural frequencies which are obtained by the present method, experiment, and FEM for the isotropic plate is shown in Table 7. It is easy to find that the results achieved by three different approaches are almost infallible, and the maximum error is less than 2%. Figure 5 exhibits the test system and model. In order to better understand the vibration characteristics of the isotropic plate, the first three mode shapes obtained by experiment are shown in Figure 6. It can be concluded that the presented method can solve the vibration behaviors of isotropic plates.


Mode numberPresentExperimentalError (%)FEMError (%)

1329.1331.40.70329.80.21
2802.3811.51.15806.60.54
31538.91564.71.681552.80.90

Tables 810 provide frequency parameters Ω belonging to the first eight bits of the orthotropic plate under different boundary conditions, and these different boundary conditions mentioned above include four kinds of classical boundary conditions and four groups of elastic boundary conditions. In this above table, the material parameters and aspect ratios a/b = 1 of the orthotropic plate are consistent with those in Table 1 and the thickness ratios are h/a = 0.3 (Table 8), h/a = 0.6 (Table 9), and h/a = 1 (Table 10), respectively. As can be seen from the table above, frequency parameters Ω obtained under the boundary conditions of CCCC are the largest, while the frequency parameter obtained under the E4E4E4E4 boundary condition is the smallest. Therefore, it is reflected that the increase of boundary stiffness of the orthotropic plate will result in the increase of frequency parameters. For the sake of better show the advantages of this method, Figures 79 provide three-dimensional modal shapes with the orthotropic plate under different boundary conditions. Unlike other methods, only two-dimensional modal shapes can be drawn, so the vibration characteristics of the orthotropic plate can be more clearly displayed. Lastly, the influence of aspect ratio b/a and thickness ratio h/a is also reported. The change regulation about nondimensional frequency parameter Ω of the orthotropic rectangular plates under the different aspect ratio b/a is shown in Figure 10. There are four sets of boundary conditions CCCC, CFCF, E2E2E2E2, and E4E4E4E4 which are used in this example. From these pictures, it can be clearly detected that the frequency parameters of orthotropic rectangular plates gradually decrease as the increase of aspect ratio a/b, and it need to be noticed that the frequency parameters almost remain unchanged when the aspect ratio increases to a critical aspect ratio. Moreover, it is found that the critical aspect ratio is directly related to the boundary conditions. Figure 11 shows the relation between the nondimensional frequency parameter Ω and the thickness ratio h/a of orthotropic rectangular plates. From this diagram, it becomes apparent that the frequency parameters on orthotropic rectangular plates rise monotonously with the raise of thickness ratio h/a, regardless of the boundary conditions or modal order. Finally, the influence of material properties on structural vibration characteristics is carried out. Figure 12 shows the histogram of the variation of the first three modes frequency parameters of the structure with different orthogonal stiffness ratios E11/E22 under different boundary conditions. The material parameters are defined as follows: Ey = Ez = 8.963 GPa, Ex/Ey = 5, 10, 15, 20, 25, 30, 35, and 40; Gxy = 6.205 GPa, Gxz = 7.101 GPa, and Gyz = 6.205 GPa; μxy = 0.3, μyz = 0.49, and μxz = 0.3; and ρ = 1605 kg/m3. The geometric parameters are defined as follows: a = b = 1 m and h = 0.5 m. It can be clearly found from Figure 12 that the frequency parameters of the structure always increase with the increase of the orthotropic stiffness value no matter what the boundary conditions of the structure are, and the trend is more obvious at the stage of smaller stiffness. Through the above discussion and study, it is easy to discover that the vibration behavior of the orthotropic rectangular plates is related to both geometric parameters and boundary conditions.


Mode numberBoundary conditions
CCCCSSSSCFCFSFSFE1E1E1E1E2E2E2E2E3E3E3E3E4E4E4E4

10.442820.332000.298580.197590.170390.169870.163840.15148
20.713730.483230.329600.276210.393450.366940.305140.16813
30.750680.483230.446530.384660.526650.526420.311040.16823
40.816780.649740.503230.483260.664140.650420.477750.22784
50.966590.650430.607140.494370.676700.677370.619660.22816
61.022140.656110.659530.559850.734080.682600.624440.23910
71.047680.897410.767500.624220.943860.714010.745910.35537
81.076570.966450.813760.631070.964910.893510.772380.40865


Mode numberBoundary conditions
CCCCSSSSCFCFSFSFE1E1E1E1E2E2E2E2E3E3E3E3E4E4E4E4

11.059860.897410.680200.555750.340790.339750.332170.32414
21.600310.966450.764020.729930.995000.893880.730780.33626
31.633880.966450.893090.769291.053291.052850.764420.33647
41.707991.297241.255880.966511.301641.366001.121050.44272
52.044651.536831.338201.211431.329861.423221.403880.44290
62.151751.585011.491641.256281.570231.519771.576750.47822
72.157861.638411.526771.258131.639411.636271.619200.85949
82.272511.888701.715511.312461.700461.660791.662791.03227


Mode numberBoundary conditions
CCCCSSSSCFCFSFSFE1E1E1E1E2E2E2E2E3E3E3E3E4E4E4E4

11.881501.610761.184451.057600.567980.566270.555630.54909
22.725681.610761.361751.256281.378691.660371.303790.56044
32.776871.687791.488511.282181.755491.718531.448040.56080
42.787672.042732.310181.317111.811131.754761.997540.68732
52.893282.137442.312991.610862.158702.279322.352620.68739
63.376542.286412.330432.042812.221192.300822.370890.79706
73.380052.557732.385912.058952.575912.389932.443141.47435
83.387072.667662.403522.113172.629572.576782.527161.67874

4. Conclusions

This paper analyzed the free vibration of orthotropic plates under arbitrary boundary conditions by semianalytical three-dimensional elastic solution. And 3D elastic theory is adopted to achieve the energy functional of orthotropic plates. The real boundary conditions of orthotropic plates are simulated by uniformly arranging three groups of linear springs on the edges of plates. The admissible displacement functions of orthotropic plates are consistently expanded as a modified Fourier series that it is based on 3D Fourier cosine series and auxiliary function. Under this framework, when the boundary conditions of the orthotropic plates change, just only the boundary parameters need to be modified and do not need a new derivation, thus greatly saving the modeling time. The astringency, accuracy, and constancy of the method of this paper are validated by means of series of numerical examples. Moreover, some new results and parameters of orthotropic plates under different boundary conditions, thickness ratios, and aspect ratios are also presented which can be used as a benchmark to help subsequent scholars engage in such research. At the same time, the results of this paper can be used as reference data for future researches in this field.

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported financially by the National Natural Science Foundation of China (Grant nos. 51109101, 51779111, and 51622902), the Key University Science Research Project of Jiangsu Province (Grant no. 17KJA416003), and the Key Laboratory Foundation of Jiangsu Province (Grant no. CJ1601).

References

  1. C.-P. Wu and Y.-C. Liu, “A review of semi-analytical numerical methods for laminated composite and multilayered functionally graded elastic/piezoelectric plates and shells,” Composite Structures, vol. 147, pp. 1–15, 2016. View at: Publisher Site | Google Scholar
  2. A. S. Sayyad and Y. M. Ghugal, “On the free vibration analysis of laminated composite and sandwich plates: a review of recent literature with some numerical results,” Composite Structures, vol. 129, pp. 177–201, 2015. View at: Publisher Site | Google Scholar
  3. A. S. Sayyad and Y. M. Ghugal, “Bending, buckling and free vibration of laminated composite and sandwich beams: a critical review of literature,” Composite Structures, vol. 171, pp. 486–504, 2017. View at: Publisher Site | Google Scholar
  4. H.-T. Thai and S.-E. Kim, “Levy-type solution for free vibration analysis of orthotropic plates based on two variable refined plate theory,” Applied Mathematical Modelling, vol. 36, no. 8, pp. 3870–3882, 2012. View at: Publisher Site | Google Scholar
  5. H.-T. Thai and S.-E. Kim, “Levy-type solution for buckling analysis of orthotropic plates based on two variable refined plate theory,” Composite Structures, vol. 93, no. 7, pp. 1738–1746, 2011. View at: Publisher Site | Google Scholar
  6. H.-T. Thai and S.-E. Kim, “Analytical solution of a two variable refined plate theory for bending analysis of orthotropic levy-type plates,” International Journal of Mechanical Sciences, vol. 54, no. 1, pp. 269–276, 2012. View at: Publisher Site | Google Scholar
  7. K. Sepahvand, S. Marburg, and H.-J. Hardtke, “Stochastic free vibration of orthotropic plates using generalized polynomial chaos expansion,” Journal of Sound and Vibration, vol. 331, no. 1, pp. 167–179, 2012. View at: Publisher Site | Google Scholar
  8. B. Abbès and Y. Q. Guo, “Analytic homogenization for torsion of orthotropic sandwich plates: application to corrugated cardboard,” Composite Structures, vol. 92, no. 3, pp. 699–706, 2010. View at: Publisher Site | Google Scholar
  9. Y. M. Ghugal and A. S. Sayyad, “Free vibration of thick orthotropic plates using trigonometric shear deformation theory,” Latin American Journal of Solids and Structures, vol. 8, no. 3, pp. 229–243, 2011. View at: Publisher Site | Google Scholar
  10. M.-H. Hsu, “Vibration analysis of orthotropic rectangular plates on elastic foundations,” Composite Structures, vol. 92, no. 4, pp. 844–852, 2010. View at: Publisher Site | Google Scholar
  11. E. Bahmyari and A. Rahbar-Ranji, “Free vibration analysis of orthotropic plates with variable thickness resting on non-uniform elastic foundation by element free Galerkin method,” Journal of Mechanical Science and Technology, vol. 26, no. 9, pp. 2685–2694, 2012. View at: Publisher Site | Google Scholar
  12. N. Valizadeh, T. Q. Bui, V. T. Vu, H. T. Thai, and M. N. Nguyen, “Isogeometric simulation for buckling, free and forced vibration of orthotropic plates,” International Journal of Applied Mechanics, vol. 5, no. 2, Article ID 1350017, 2013. View at: Publisher Site | Google Scholar
  13. Y. F. Xing and T. F. Xu, “Solution methods of exact solutions for free vibration of rectangular orthotropic thin plates with classical boundary conditions,” Composite Structures, vol. 104, pp. 187–195, 2013. View at: Publisher Site | Google Scholar
  14. A. M. Zenkour, “Bending of orthotropic plates resting on Pasternak’s foundations using mixed shear deformation theory,” Acta Mechanica Sinica, vol. 27, no. 6, pp. 956–962, 2011. View at: Publisher Site | Google Scholar
  15. B. Bakhadda, M. B. Bouiadjra, F. Bourada et al., “Dynamic and bending analysis of carbon nanotube-reinforced composite plates with elastic foundation,” Wind and Structures, vol. 27, no. 5, pp. 311–324, 2018. View at: Google Scholar
  16. F. Bourada, A. A. Bousahla, M. Bourada et al., “Dynamic investigation of porous functionally graded beam using a sinusoidal shear deformation theory,” Wind and Structures, vol. 28, no. 1, pp. 19–30, 2019. View at: Google Scholar
  17. A. Boukhari, H. A. Atmane, A. Tounsi, E. A. Adda Bedia, and S. R. Mahmoud, “An efficient shear deformation theory for wave propagation of functionally graded material plates,” Structural Engineering and Mechanics, vol. 57, no. 5, pp. 837–859, 2016. View at: Publisher Site | Google Scholar
  18. H. Fourn, H. A. Atmane, M. Bourada et al., “A novel four variable refined plate theory for wave propagation in functionally graded material plates,” Steel and Composite Structures, vol. 27, no. 1, pp. 109–122, 2018. View at: Google Scholar
  19. S. O. Papkov and J. R. Banerjee, “A new method for free vibration and buckling analysis of rectangular orthotropic plates,” Journal of Sound and Vibration, vol. 339, pp. 342–358, 2015. View at: Publisher Site | Google Scholar
  20. M.-F. Liu, T.-P. Chang, and Y.-H. Wang, “Free vibration analysis of orthotropic rectangular plates with tapered varying thickness and winkler spring foundation,” Mechanics Based Design of Structures and Machines, vol. 39, no. 3, pp. 320–333, 2011. View at: Publisher Site | Google Scholar
  21. D. Shi, Q. Wang, X. Shi, and F. Pang, “A series solution for the in-plane vibration analysis of orthotropic rectangular plates with non-uniform elastic boundary constraints and internal line supports,” Archive of Applied Mechanics, vol. 85, no. 1, pp. 51–73, 2015. View at: Publisher Site | Google Scholar
  22. A. A. Jafari and S. A. Eftekhari, “An efficient mixed methodology for free vibration and buckling analysis of orthotropic rectangular plates,” Applied Mathematics and Computation, vol. 218, no. 6, pp. 2670–2692, 2011. View at: Publisher Site | Google Scholar
  23. B. Liu and Y. Xing, “Exact solutions for free vibrations of orthotropic rectangular Mindlin plates,” Composite Structures, vol. 93, no. 7, pp. 1664–1672, 2011. View at: Publisher Site | Google Scholar
  24. B. Liu and Y. Xing, “Comprehensive exact solutions for free in-plane vibrations of orthotropic rectangular plates,” European Journal of Mechanics—A/Solids, vol. 30, no. 3, pp. 383–395, 2011. View at: Publisher Site | Google Scholar
  25. H.-H. Chang and J.-Q. Tarn, “Three-dimensional elasticity solutions for rectangular orthotropic plates,” Journal of Elasticity, vol. 108, no. 1, pp. 49–66, 2012. View at: Publisher Site | Google Scholar
  26. T. M. Teo and K. M. Liew, “Three-dimensional elasticity solutions to some orthotropic plate problems,” International Journal of Solids and Structures, vol. 36, no. 34, pp. 5301–5326, 1999. View at: Publisher Site | Google Scholar
  27. A. Tian, R. Ye, and Y. Chen, “A new higher order analysis model for sandwich plates with flexible core,” Journal of Composite Materials, vol. 50, no. 7, pp. 949–961, 2016. View at: Publisher Site | Google Scholar
  28. A. Tian, R. Ye, P. Ren et al., “New higher-order models for sandwich plates with a flexible core and their accuracy assessment,” International Journal of Structural Stability and Dynamics, vol. 19, no. 3, Article ID 1950024, 2019. View at: Publisher Site | Google Scholar
  29. H. Li, F. Pang, X. Miao, S. Gao, and F. Liu, “A semi analytical method for free vibration analysis of composite laminated cylindrical and spherical shells with complex boundary conditions,” Thin-Walled Structures, vol. 136, pp. 200–220, 2019. View at: Publisher Site | Google Scholar
  30. H. Li, F. Pang, and H. Chen, “A semi-analytical approach to analyze vibration characteristics of uniform and stepped annular-spherical shells with general boundary conditions,” European Journal of Mechanics—A/Solids, vol. 74, pp. 48–65, 2019. View at: Publisher Site | Google Scholar
  31. H. Li, F. Pang, H. Chen, and Y. Du, “Vibration analysis of functionally graded porous cylindrical shell with arbitrary boundary restraints by using a semi analytical method,” Composites Part B: Engineering, vol. 164, pp. 249–264, 2019. View at: Publisher Site | Google Scholar
  32. H. H. Abdelaziz, M. A. A. Meziane, A. A. Bousahla et al., “An efficient hyperbolic shear deformation theory for bending, buckling and free vibration of FGM sandwich plates with various boundary conditions,” Steel and Composite Structures, vol. 25, no. 6, pp. 693–704, 2017. View at: Google Scholar
  33. B. Karami, M. Janghorban, and A. Tounsi, “Galerkin’s approach for buckling analysis of functionally graded anisotropic nanoplates/different boundary conditions,” Engineering with Computers, 2018. View at: Publisher Site | Google Scholar
  34. M. A. A. Meziane, H. H. Abdelaziz, and A. Tounsi, “An efficient and simple refined theory for buckling and free vibration of exponentially graded sandwich plates under various boundary conditions,” Journal of Sandwich Structures & Materials, vol. 16, no. 3, pp. 293–318, 2014. View at: Publisher Site | Google Scholar
  35. J. Cui, R. Ye, N. Zhao, J. Wu, and M. Wang, “Assessment on energy absorption of double layered and sandwich plates under ballistic impact,” Thin-Walled Structures, vol. 130, pp. 520–534, 2018. View at: Publisher Site | Google Scholar
  36. K. Liu, B. Liu, R. Villavicencio, Z. Wang, and C. Guedes Soares, “Assessment of material strain rate effects on square steel plates under lateral dynamic impact loads,” Ships and Offshore Structures, vol. 13, no. 2, pp. 217–225, 2018. View at: Publisher Site | Google Scholar
  37. H. Li, F. Pang, X. Miao, and Y. Li, “Jacobi-Ritz method for free vibration analysis of uniform and stepped circular cylindrical shells with arbitrary boundary conditions: a unified formulation,” Computers & Mathematics with Applications, vol. 77, no. 2, pp. 427–440, 2019. View at: Publisher Site | Google Scholar
  38. H. Zhang, D. Shi, S. Zha, and Q. Wang, “Vibro-acoustic analysis of the thin laminated rectangular plate-cavity coupling system,” Composite Structures, vol. 189, pp. 570–585, 2018. View at: Publisher Site | Google Scholar
  39. H. Zhang, D. Shi, S. Zha, and Q. Wang, “Parameterization study on the moderately thick laminated rectangular plate-cavity coupling system with uniform or non-uniform boundary conditions,” Composite Structures, vol. 194, pp. 537–554, 2018. View at: Publisher Site | Google Scholar
  40. H. Zhang, D. Shi, S. Zha, and Q. Wang, “A simple first-order shear deformation theory for vibro-acoustic analysis of the laminated rectangular fluid-structure coupling system,” Composite Structures, vol. 201, pp. 647–663, 2018. View at: Publisher Site | Google Scholar
  41. J. Zhao, F. Xie, A. Wang, C. Shuai, J. Tang, and Q. Wang, “Vibration behavior of the functionally graded porous (FGP) doubly-curved panels and shells of revolution by using a semi-analytical method,” Composites Part B: Engineering, vol. 157, pp. 219–238, 2019. View at: Publisher Site | Google Scholar

Copyright © 2019 Jie Cui et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

338 Views | 251 Downloads | 1 Citation
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.