Dynamic Stability Analysis in Hybrid Nanocomposite Polymer Beams Reinforced by Carbon Fibers and Carbon Nanotubes

The objective of this innovative research is assessment of dynamic stability for a hybrid nanocomposite polymer beam. The considered beam formed by multiphase nanocomposite, including polymer–carbon nanotubes (CNTs)–carbon fibers (CFs). Hence, as to compute the effective material characteristics related to multiphase nanocomposite layers, the Halpin–Tsai model, as well as micromechanics equations are employed. To model the structure realistically, exponential shear deformation beam theory (ESDBT) is applied and using energy methods, governing equations are achieved. Moreover, differential quadrature method (DQM) as well as Bolotin procedures are used for solving the obtained governing equations and the dynamic instability region (DIR) relative to the beam is determined. To extend this novel research, various parameters pinpointing the influences of CNT volume fraction, CFs volume percent, boundary edges as well as the structure’s geometric variables on the dynamic behavior of the beam are presented. The results were validated with the theoretical and experimental results of other published papers. The outcomes reveal that increment of volume fraction of CNT is able to shift DIR to more amounts of frequency. Further, rise of carbon fibers volume percent leads to increase the excitation frequency of this structure.


Introduction
When a structure is subjected to periodic and pulsatile loads, it is well known that the ordinary forced response will lead to dynamic instability. This subject is very important since the periodic loads may cause parametric vibrations, which may damage the structures. Nanocomposite plates can be used in different industries such as aircrafts and automobile which may be subjected to periodic and pulsatile loads. However, dynamic stability of nanocomposite plates is a novel topic which should be studied.
Different papers were brought into investigation regarding polymer nanocomposites (PNCs) special applications [1], device usages [2], its processing technologies for future usages [3][4][5][6] and safety analysis [7][8][9]. Pastoriza-Santos et al. [10] discussed the processes of manufacturing and different usages of plasmonic PNC ranging from light-harvesting improvement to tracing biologic as well as chemicals molecules. Leon et al. [11] brought up significant proficiency of PNC in terms of additive manufacturing usages. It is found that PNC are preferable due to their significant mechanical, chemical and thermal characteristics under harsh conditions. In the aspect of energy density, Liu et al. [12] analyzed dielectric PNC having multilayered structure. They claimed that the discharged energy density of this multilayered structure is the greatest discharged energy densities which were presented. Nevertheless, there are limitations in the conduction of polymers that Abbasi et al. [13] discussed and offered some advices. It is believed that there are constraints related to electrically conductive polymers and considering carbon-based nanoparticles can push back them. Hence, this research reviewed some improvements in this field.
Further, it should be noticed that hybrid composites refer to specific substances made up of organic and inorganic ingredients and nanocomposite term, points at hybrid inorganic-organic composites having components in nano scale. These materials are utilized in formations of structures because of their striking properties such as safety [14][15][16] and reliability [17][18][19]. Xia and Lo [20] offered a method by which ibuprofen (IBU)considered as organic pollutant in flowing waters, can be eliminated from the water using hybrid nanocomposites. It is shown that hybrid nanocomposite matrix can be increased the strength of organic pollutant. Filippov et al. [21] presented a novel method for characterizing hybrid nanocomposites. Ebrahimi and Dabbagh [22] carried out a research regarding vibration of hybrid nanocomposite plates. Using classical plate theory and Hamilton principle, the governing equations are obtained and solved analytically. It is observed that hybrid nanocomposites are able to resist high amounts of frequency compared to the ordinary composites. Dabbagh et al. [23] reported thermal buckling of hybrid nanocomposites beam. Using refined plate theory as well as Hamilton's principle, the governing equations are obtained and utilizing Galerkin's procedure, the motion equations are solved. It is shown that hybrid nanocomposite structures are tougher toward critical buckling loads than conventional composites. Moreover, creep as well as viscoelastic behaviors of hybrid glass/epoxy nanocomposites are studied experimentally by Salehi and colleague [24].
Without doubt, besides invented methods as to tightening structures in various domains, reinforcing structures via nanocomposites and CFs have been featured strongly since last decade. Nanocomposites such as CNTs have been preferable materials to be added into structures for this intention due to high fracture resistance and capability of vibration damping. Regardless of nanotube's types containing single walled carbon nanotubes (SWCNTs) and multi walled carbon nanotubes (MWCNTs) and their branches, they possess notable properties in different fields such as significant electrical and thermal conductivities, strength and elasticity and energy storage that have made them desirable. Moreover, investigations regarding CNTs and their composites are of interest between scientists to enhance damping properties of these materials, especially epoxy, in terms of loss modulus [25,26] or damping ratio [27,28]. There are researches exploring properties of nanocomposites Vinson [29] and Herman et al. [30]. Gul and Aydogdu [31] investigated wave dispersion of DWCNTs. They explained that, through Doublet Mechanics, the wave dispersion was studied which had precise outcomes for flexural as well as axial wave dispersion in nanotubes. Vibration analysis of SWCNTs was carried out by Avramov [32]. Using Sanders-Koiter shell methods-as well as nonlocal elasticity-the motion equations were derived and solved numerically via the Galerkin procedure. In another study, Bian and Wang [33] carried out a research as to buckling of DWCNTs in a thermal environment. The influence of finite temperature and various nonlocal effects have been studied. They revealed that buckling of DWCNT declines whereas the finite temperature rises. Natural frequency analysis of FG-CNTs curved shell panel which is shallow has been done by Mehar et al. [34]. In this study, higher-order theory is used for obtaining governing equations and finite element method (FEM) is utilized for solving equations. Further, the precision of this numerical method is compared with different numerical procedures. Zhang et al. [35] analyzed deflection of FG-CNTRC plate which is on the Pasternak medium. First-order shear deformation theory (FSDT) besides Ritz methods are utilized in order to obtain and solve governing equations, respectively. In this work, the influences relevant to CNT distributions and volume fraction and foundation up deflection are analyzed. Jiao et al. [36] discussed dynamic buckling behavior of FG-CNTRC shell. The structure is under time-varying load and FSDT is applied to derive motion equations. Galerkin method is used for solving the equations and the influence of dispersion of CNTs and Polymers 2021, 13, 106 3 of 15 their volume fraction and other parameters are taken into account. Mirzaei and Kiani [37] investigated buckling evaluation of FG-CNTRC plates in thermal environment. In this study, besides various parameters including different boundary conditions and aspect ratio, the influence of CNTs are discussed. They indicated that FG-X distribution type is the better choice. In another research, wave dispersion of FG nanocomposites which are reinforced through CNTs are investigated by Janghorban and Nami [38]. Dynamic analysis of structures induced by non-conservative loading and electromagnetic interactions has been studied in the literature. Based on the equivalent single-layer linear theory for laminated shells, Kolahchi et al. [39] presented dynamic buckling of sandwich nanoplate subjected to harmonic compressive load based on nonlocal elasticity theory. Mikhasev et al. [40] investigated free and forced vibrations of thin cylindrical sandwich panels with magnetorheological core. Hajmohammad et al. [41] analyzed dynamic response of sandwich plates under the blast load based on sinusoidal shear deformation theory. Malikan and Eremeyev [42] investigated the effect of flexoelectricity on dynamic response of piezoelectric nanobeam assuming internal viscoelasticity. Malikan et al. [43] studied torsional stability of a nanocomposite shell subjected to magnetic field utilizing nonlocal strain gradient theory. Keshtegar et al. [44] presented dynamic stability of nanocomposite-truncated sandwich conical shells based on differential cubature method.
To date, no research has been conducted studying and analyzing the dynamic behavior of a hybrid nanocomposite polymer beam reinforced via CFs as well as CNTs using ESDBT. In this study, we examine the beams composed by multiphase nanocomposites, including carbon/fiber/CNT/ polymer which the responses of beams are extracted using micromechanics-as well as Halpin-Tsai equations. With respect to DQM and Bolotin procedures, the obtained motion equations are solved, and DIR are attained, respectively. Distinct variables are investigated, such as the influences of CNT volume fraction, CF volume percent, boundary edges as well as the structure's geometric variables on the dynamic behavior of the beam. equations. Galerkin method is used for solving the equations and the influence of dispersion of CNTs and their volume fraction and other parameters are taken into account. Mirzaei and Kiani [37] investigated buckling evaluation of FG-CNTRC plates in thermal environment. In this study, besides various parameters including different boundary conditions and aspect ratio, the influence of CNTs are discussed. They indicated that FG-X distribution type is the better choice. In another research, wave dispersion of FG nanocomposites which are reinforced through CNTs are investigated by Janghorban and Nami [38]. Dynamic analysis of structures induced by non-conservative loading and electromagnetic interactions has been studied in the literature. Based on the equivalent single-layer linear theory for laminated shells, Kolahchi et al. [39] presented dynamic buckling of sandwich nanoplate subjected to harmonic compressive load based on nonlocal elasticity theory. Mikhasev et al. [40] investigated free and forced vibrations of thin cylindrical sandwich panels with magnetorheological core. Hajmohammad et al. [41] analyzed dynamic response of sandwich plates under the blast load based on sinusoidal shear deformation theory. Malikan and Eremeyev [42] investigated the effect of flexoelectricity on dynamic response of piezoelectric nanobeam assuming internal viscoelasticity. Malikan et al. [43] studied torsional stability of a nanocomposite shell subjected to magnetic field utilizing nonlocal strain gradient theory. Keshtegar et al. [44] presented dynamic stability of nanocomposite-truncated sandwich conical shells based on differential cubature method.

Problem Definition
To date, no research has been conducted studying and analyzing the dynamic behavior of a hybrid nanocomposite polymer beam reinforced via CFs as well as CNTs using ESDBT. In this study, we examine the beams composed by multiphase nanocomposites, including carbon/fiber/CNT/ polymer which the responses of beams are extracted using micromechanics-as well as Halpin-Tsai equations. With respect to DQM and Bolotin procedures, the obtained motion equations are solved, and DIR are attained, respectively. Distinct variables are investigated, such as the influences of CNT volume fraction, CF volume percent, boundary edges as well as the structure's geometric variables on the dynamic behavior of the beam.  In this section, in order to model the structure mathematically, ESDBT is employed and hence, the displacement vectors are given by [45]: In this section, in order to model the structure mathematically, ESDBT is employed and hence, the displacement vectors are given by [45]:

Problem Definition
Bx´ψ px, tq¯, u 2 px, z, tq " 0, u 3 px, z, tq " wpx, tq. In above relation, u 1 and u 3 respectively describe mid-plane displacements in longitudinal and thickness directions. Moreover, ψ defines rotation related to cross section and Φpzq represents shape function of the beam that can be written as in which h refers to thickness of structure. Furthermore, the strain components of this structure are described as It should be noted that stress-strain relations can be written as where C 11 and C 44 express the elastic constants achieved through Halpin-Tsai model.

Modeling CNT/Fiber/Polymer Multiphase Nanocomposite
Carbon nanotubes were grown directly on carbon fibers using chemical vapor deposition. When embedded in a polymer matrix, the change in length scale of carbon nanotubes relative to carbon fibers results in a multiscale composite [46]. Incorporation of micromechanics as well as Halpin-Tsai model [47] are utilized for reaching polymer beam's equivalent material properties in two phases. Orthotropic effective characteristics related to CNT-reinforced multi-phase laminates are given by in which G, E and ρ respectively express shear modulus, Young's modulus and mass density; V as well as ν represent volume fraction and Poisson's ratio, respectively. Further, the subscript and superscript MNC and F, respectively refer to matrix of nanocomposite and fibers. With respect to Halpin-Tsai equations, nanocomposite's elastic modulus is written as follow [48]: in which in which V M as well as E M respectively describe matrix's volume fraction and Young's modulus. Moreover V CN , E CN , t CN , d CN and CN respectively define volume fraction, Young's modulus, thickness, outer diameter and length correlative to CNTs. The CNT volume fraction is given by in which ρ m and ρ CN respectively hint at mass density of matrix and CNTs and w CN is mass fraction. Likewise, Poisson's ratio, mass density and shear modulus related to MNC are written as below: in which ν MNC as well as ν M respectively represent Poisson's ratio related to MNC and matrix. It is mentioned that because of the CNT's small extent, the Poisson's ratio related to matrix and MNC will be assumed analogous [49].

Motion Equations
In this part, the structure's strain energy is described as below: Likewise, by introducing Equations (3) and (4) into Equation (19), we have Therefore, stress resultants are given by In order to extend above equations, Equations (3)-(6) are substituted into Equations (21)- (24) and the outcomes are presented in Appendix A. For the next step, the kinetic energy for this structure are written as Substituting Equation (1) into Equation (25) yields in which ρ refers to the density of the beam. The inertia moments are expressed as follow Factoring in Equation (27), the Equation (26) is rewritten as Eventually, motion equations are achieved utilizing Hamilton's principle: By introducing Equations (20), (28) into Equation (29), the governing equations are expressed as follow δu : δw where P " αP cr`β P cr cospωtq in which α, β, P cr and ω are static load factor, dynamic load factor, critical load and frequency, respectively. For expansion of above equations, Equations (A1)-(A4) are substituted into Equations (30)-(32) and hence, they are expressed as below: At the end, for investigating different cases of this beam, various boundary edges are considered which can be computed as below

Solving Procedure
As mentioned, for solution of motion equations besides determination of DIR, DQM is employed. In the method, the various orders of the beam's differential equations are converted to a set of algebraic equations based on weighting coefficients. In other words, with respect to this precise procedure, one derivative of a function at one contemplated separate point can be computed, utilizing the sum of the extent of function at every separate point opted out in the scope of the solution. The approximation of the derivative function can be expressed in a general form as [50][51][52] in which f pxq refers to the function, N describes number of points, x i represents an instance point of the function scope, f i expresses the amount of the function at ith instance point and C ij denotes weighting coefficients. Moreover, selecting grid points, as well as weighting coefficients would be essential parameter in gaining precise consequences. Grid points can be contemplated through Chebyshev polynomials as With respect to the Chebyshev polynomials, grid points would be denser in the neighbor of boundaries. The weighting coefficients are achieved using following relation where With respect to DQM, motion equations in matrix are expressed form as below: , . - where rMs,rKs as well as rCs define mass matrix, stiffness matrix and damper matrix, respectively. Further, td b u and td d u define boundary points as well as domain points, respectively. With respect to Bolotin's method, elements related to tYu can be described in Fourier series having period 2T as tYu " In this part, by substituting Equation (44) into Equation (43) and setting the factors of cosine and sine-as well as the sum of constant terms to zero-we havěˇˇˇˆr

Ks´ˆα˘β
2˙P cr rKs G˙´r Ms Hence, as to gain variation of ω and DIR based on β, the mentioned relation can be solved according to eigenvalue problem.

Numerical Consequences
This part is presented to investigate and analyze the influence of diverse variables on the dynamic behavior of the beam with length of L = 2 m and thickness of h = 30 cm. As noted previously, the hypothesized beam is composed of epoxy and CFs. In addition, CNTs are considered the reinforcements and properties of these materials are assumed, according to [34]. Young's modulus, Poisson's ratio and density of epoxy are E M " 3.51 GPa, ν M " 0.3 and ρ M " 1200 Kg{m 3 , respectively. Further, material characteristics including Young's modulus, Shear modulus, Poisson's ratio and density related to the CFs, respectively are E F 11 " 233.05 GPa, E F 22 " 23.1 GPa, G F 12 " 8.96 GPa, ν F " 0.2 and ρ F " 1750 Kg{m 3 . Likewise, Young's modulus, Poisson's ratio, density, outer diameter, thickness and length of CNT, respectively are E CN " 640 GPa, ν CN " 0.27, ρ CN " 1350 Kg{m 3 , d CN " 1.4 nm, t CN " 0.34 nm and CN " 25ˆ10´6 m.

Validation
To date, in spite of proliferation of researches in mechanic and material fields, there is no work analyzing dynamic behavior of hybrid nanocomposite polymer beams reinforced via CFs and CNTs using ESDBT and determining DIR. Therefore, in order to demonstrate the precision-as well as validity of the considered theories and numerical methods-some characteristics containing CNTs, CFs and ESDBT are disregarded to fit in this beam with the research published by Joubaneh et al [53]. In this research, vibration of a three-layer beam is studied based on theoretical and experimental methods. For this purpose, a beam with shear modulus of 22.1 GPa, Poisson's ratio of 0.3, density of 60 kg/m 3 , thickness of 15 mm and length of 260 mm is assumed which is covered with two layers at the top and bottom of core with Young's modulus of 210 GPa, Poisson's ratio of 0.3, density of 7900 kg/m 3 and thickness of 1.9 mm. For experimental analyses, a clamped-free boundary condition of beam with three layers is installed on a VDL shaker (B & K V830-335-SPA16K) using a head plate and a fixture, both made of aluminum, as shown in Figure 2. Hence, as to gain variation of  and DIR based on  , the mentioned relation can be solved according to eigenvalue problem.

Numerical Consequences
This part is presented to investigate and analyze the influence of diverse variables on the dynamic behavior of the beam with length of L = 2 m and thickness of h = 30 cm. As noted previously, the hypothesized beam is composed of epoxy and CFs. In addition, CNTs are considered the reinforcements and properties of these materials are assumed, according to [34]

Validation
To date, in spite of proliferation of researches in mechanic and material fields, there is no work analyzing dynamic behavior of hybrid nanocomposite polymer beams reinforced via CFs and CNTs using ESDBT and determining DIR. Therefore, in order to demonstrate the precision-as well as validity of the considered theories and numerical methods-some characteristics containing CNTs, CFs and ESDBT are disregarded to fit in this beam with the research published by Joubaneh et al [53]. In this research, vibration of a three-layer beam is studied based on theoretical and experimental methods.  The theoretical and experimental results of [53] are shown in Table 1: the outcomes are validated with our numerical method. As can be seen, the results of this study are close to the theoretical and experimental results [53,54]. The theoretical and experimental results of [53] are shown in Table 1: the outcomes are validated with our numerical method. As can be seen, the results of this study are close to the theoretical and experimental results [53,54]. In another validation, the buckling of nanocomposite beam is studied. A Poly methyl methacrylate (PMMA) with Young's modulus, Poisson's ratio and density of epoxy are E M " 2.5 GPa, ν M " 0.3 and ρ M " 1190 Kg{m 3 , respectively is assumed which is reinforced with 0.12 carbon nanotube with Young's modulus, Poisson's ratio and the density of epoxy are E CN 11 " 600 GPa, E CN 22 " 10 GPa, ν CN " 0.19 and ρ CN " 1400 Kg{m 3 , respectively. The non-dimensional buckling load P " P{hE M for different boundary condition is presented in Table 2. It is shown that the result of this study are in good agreement with Asadi and Wang [55] and Yas and Samadi [56].  Table 3 presents the dimensionless frequency of the beam reinforced by CNTs, considering the material properties the same as those mentioned in before validation. The results are reported for different theories of first-order shear deformation theory (FSDT), third order shear deformation theory (TSDT), exponential order shear deformation theory (ESDT), higher order shear deformation theory (HSDT) and trigonometric shear deformation theory (TrSDT) which are presented by Wattanasakulpong and Ungbhakorn [57]. It is seen that the accuracy of the obtained results is good.

Parametric Study
In this section, the dynamic stability of the structure is studied for different parameters. All of figures show the excitation frequency against dynamic load factor (i.e., β) which represent the DIR of structure. In these figures, the regions inside and outside the boundary curves correspond to unstable (parametric resonance) and stable regions, respectively.
At first, so as to assess the precision as well as convergence of applied numerical method, Figure 3 is depicted. As observed, excitation frequencies were obtained for various grid points. It is clear that fifteen grid points can be factored in the precise outcomes in this numerical method.  Figure 4 presents the influence of CNT weight-percentage on DIR. Likewise, it is worth noting that the inside sections of triangular shapes represent unstable areas and consequently, the outside sections indicate stable areas of this beam. It is obvious that rise of CNT weight percentage triggers increase in excitation frequency of this structure. To put it differently, the DIR takes place at higher extents of excitation frequencies while increasing weight percentage of CNTs. It is justified with this fact that increase in CNT weight percentage leads to enhancement of stiffness as well as bending rigidity related to this structure. In this work, various boundary edges were factored in to extend this research for different cases. Therefore, Figure 5 shows the influence of different boundary edges containing clamped-clamped (CC), clamped-simply (CS) and simply-simply (SS) boundary edges on the dynamic behaviors of this system. As observed, boundary edges have an undisputed effect on the DIR of the beam, and it  Figure 4 presents the influence of CNT weight-percentage on DIR. Likewise, it is worth noting that the inside sections of triangular shapes represent unstable areas and consequently, the outside sections indicate stable areas of this beam. It is obvious that rise of CNT weight percentage triggers increase in excitation frequency of this structure. To put it differently, the DIR takes place at higher extents of excitation frequencies while increasing weight percentage of CNTs. It is justified with this fact that increase in CNT weight percentage leads to enhancement of stiffness as well as bending rigidity related to this structure.
Polymers 2020, X, x FOR PEER REVIEW 10 of 16 Figure 3. Precision and convergence of differential quadrature method (DQM). Figure 4 presents the influence of CNT weight-percentage on DIR. Likewise, it is worth noting that the inside sections of triangular shapes represent unstable areas and consequently, the outside sections indicate stable areas of this beam. It is obvious that rise of CNT weight percentage triggers increase in excitation frequency of this structure. To put it differently, the DIR takes place at higher extents of excitation frequencies while increasing weight percentage of CNTs. It is justified with this fact that increase in CNT weight percentage leads to enhancement of stiffness as well as bending rigidity related to this structure. In this work, various boundary edges were factored in to extend this research for different cases. Therefore, Figure 5 shows the influence of different boundary edges containing clamped-clamped (CC), clamped-simply (CS) and simply-simply (SS) boundary edges on the dynamic behaviors of this system. As observed, boundary edges have an undisputed effect on the DIR of the beam, and it In this work, various boundary edges were factored in to extend this research for different cases. Therefore, Figure 5 shows the influence of different boundary edges containing clamped-clamped (CC), clamped-simply (CS) and simply-simply (SS) boundary edges on the dynamic behaviors of this system. As observed, boundary edges have an undisputed effect on the DIR of the beam, and it is vivid that this structure with CC boundary edges contains higher DIR compared to CS as well as SS having lower quantities. In fact, it raises bending rigidity.
Polymers 2020, X, x FOR PEER REVIEW 11 of 16 is vivid that this structure with CC boundary edges contains higher DIR compared to CS as well as SS having lower quantities. In fact, it raises bending rigidity.

Figure 5.
Effect of distinct boundary edges on DIR. Figure 6 illustrates to indicate the significance related to presence of CFs and its volume percent and consequently its influence on the dynamic behaviors of the system. As expected, increase in volume percent of CFs results in shift of DIR to higher extents of excitation frequencies owing to rise of system's stiffness. Having been plotted below, Figures 7 and 8 are the exemplary states to ascertain the influence of geometric variables on the dynamic behavior of the beam. Figures 7 and 8, respectively evaluate the changes of length and thickness of structure. It is construed that rise of length parameter of polymeric beam results in alteration of DIR into the lower extents. In fact, this change can affect  Figure 6 illustrates to indicate the significance related to presence of CFs and its volume percent and consequently its influence on the dynamic behaviors of the system. As expected, increase in volume percent of CFs results in shift of DIR to higher extents of excitation frequencies owing to rise of system's stiffness.
Polymers 2020, X, x FOR PEER REVIEW 11 of 16 is vivid that this structure with CC boundary edges contains higher DIR compared to CS as well as SS having lower quantities. In fact, it raises bending rigidity.

Figure 5.
Effect of distinct boundary edges on DIR. Figure 6 illustrates to indicate the significance related to presence of CFs and its volume percent and consequently its influence on the dynamic behaviors of the system. As expected, increase in volume percent of CFs results in shift of DIR to higher extents of excitation frequencies owing to rise of system's stiffness. Having been plotted below, Figures 7 and 8 are the exemplary states to ascertain the influence of geometric variables on the dynamic behavior of the beam. Figures 7 and 8, respectively evaluate the changes of length and thickness of structure. It is construed that rise of length parameter of polymeric beam results in alteration of DIR into the lower extents. In fact, this change can affect Having been plotted below, Figures 7 and 8 are the exemplary states to ascertain the influence of geometric variables on the dynamic behavior of the beam. Figures 7 and 8, respectively evaluate the changes of length and thickness of structure. It is construed that rise of length parameter of polymeric beam results in alteration of DIR into the lower extents. In fact, this change can affect characteristics related to the structure and cause softness of whole structure. Inversely, increment of the thickness leads to improvement in the stiffness of system and move in DIR to the higher extents. characteristics related to the structure and cause softness of whole structure. Inversely, increment of the thickness leads to improvement in the stiffness of system and move in DIR to the higher extents.

Conclusions
This concerned the evaluation of dynamic instability related to hybrid nanocomposite polymer beam. The considered beam composed by multiphase nanocomposite, include polymer-CNTs-CFs. Halpin-Tsai as well as micromechanics equations were applied to compute effective material characteristics correlative to multiphase nanocomposite layers. Employing ESDBT besides Hamilton's principle, governing equations were obtained. Further, the DQM and Bolotin procedures were utilized to solve the governing equations and achieve dynamic behavior of polymer beam. Different parameters such as various boundary edges, geometric parameters, CNTs and CFs volume Polymers 2020, X, x FOR PEER REVIEW 12 of 16 characteristics related to the structure and cause softness of whole structure. Inversely, increment of the thickness leads to improvement in the stiffness of system and move in DIR to the higher extents.

Conclusions
This concerned the evaluation of dynamic instability related to hybrid nanocomposite polymer beam. The considered beam composed by multiphase nanocomposite, include polymer-CNTs-CFs. Halpin-Tsai as well as micromechanics equations were applied to compute effective material characteristics correlative to multiphase nanocomposite layers. Employing ESDBT besides Hamilton's principle, governing equations were obtained. Further, the DQM and Bolotin procedures were utilized to solve the governing equations and achieve dynamic behavior of polymer beam. Different parameters such as various boundary edges, geometric parameters, CNTs and CFs volume

Conclusions
This concerned the evaluation of dynamic instability related to hybrid nanocomposite polymer beam. The considered beam composed by multiphase nanocomposite, include polymer-CNTs-CFs. Halpin-Tsai as well as micromechanics equations were applied to compute effective material characteristics correlative to multiphase nanocomposite layers. Employing ESDBT besides Hamilton's principle, governing equations were obtained. Further, the DQM and Bolotin procedures were utilized to solve the governing equations and achieve dynamic behavior of polymer beam. Different parameters such as various boundary edges, geometric parameters, CNTs and CFs volume fractions were brought up and their influences were illustrated on DIR. The prominent outcomes are presented as following

‚
Increments of weight-percentage of CNTs can lead to shift of DIR to the right; ‚ Presence of CFs as well as CNTs play paramount role in dynamic behavior of hybrid polymer structure and raise the excitation frequency; Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.