The in-plane free vibration of sector, annular, and circular plates is investigated by isogeometric finite element approach on the basis of nonuniform rational B-spline (NURBS) basis functions. Under the current framework, both the displacement field and geometry of the sector, annular, and circular plates are modeled by NURBS basis functions to bridge the gap between the design of geometry and the analysis of variable field. The NURBS basic functions can not only preserve the exact geometry of sector, annular, and circular plates, but also provide higher continuity of basis function and its derivatives. The governing equations can be derived by employing the principle of virtual work and the desired solutions are obtained by using the Arnoldi Method. Several numerical examples of sector, annular, and circular plates are performed, and three refinement schemes (the h-, p-, and k-refinement strategies) are applied to demonstrate the convergence. Then the effectiveness and accuracy of the proposed approach are validated through comparisons with results obtained from the finite element analysis as well as open available literature. On this basis, some new numerical results of frequency parameters with mode shapes are shown and may be used as benchmark results for the vibration investigation in the future. In addition, the effects of sector angles and ratio of inside to outside radii on the in-plane vibration of sector, annular, and circular plates under different boundary conditions are fully demonstrated.
National Natural Science Foundation of China5177512551279035Fundamental Research Funds for the Central Universities of ChinaHEUCFG2017131. Introduction
The sector, annular, and circular plates are typical structural components used in engineering widely. The flexural vibration of sector, annular, and circular plates has been studied by researchers and design engineers with considerable interest. A number of investigations [1–5] have been devoted to flexural vibrations of structures, maybe because the flexural vibration has lower resonant frequencies and a decisive role in terms of fluid-structure coupling. However, the vibrations of structures also contain in-plane parts, which often appear in high frequency motions and large coupled structures. When the sound radiation and energy transmission of coupled structures are considered, the importance of in-plane vibration is nonnegligible. Consequently, it is of great significance to obtain deep in-plane vibration comprehensions of sector, annular, and circular plates.
Xing and Liu [6] applied a Rayleigh quotient variational principle to study the free in-plane vibration for rectangular plate. In the paper, all classical boundary restraints which included two various kinds of simple supports were taken into consideration. Bercin and Langley [7] analyzed the in-plane vibration problem of plate structures by adopting the dynamic stiffness technique and classical finite element assembled technique. Nefovska-Danilovic and Petronijevic [8] also used this method to investigate the in-plane free vibration of arbitrarily restrained isotropic rectangular plate. Gorman [9] used the method of superposition to solve the problem of free in-plane vibration of plates; then he used this method to solve in-plane vibration of Levy-type plates [10] which have a pair of simply supported opposite edges at least. Andrianov et al. [11] adopted homotopy perturbation method to study the in-plane vibration of rectangular plate. Chen et al. [12] applied Chebyshev-Lagrangian approach to study in- and out-of-plane vibration of plate with cutout. Mohazzab and Dozio [13] used the spectral collocation method to analyze in-plane vibration of isotropic skewed geometrical plates. Some other efforts [14–18] have also been devoted to the in-plane vibration problem of plate structures.
It is needed to point out that a great majority of existing investigations in aforementioned issues are restricted within rectangular plates. The related work about the in-plane vibration of sector, annular, and circular plates is very limited. Irie et al. [19] studied the natural in-plane vibration frequencies for annular and circular plates with clamped and free boundary conditions by using the transfer matrix method. The Rayleigh-Ritz method was employed by Bashmal et al. [20] to investigate the natural in-plane vibration characteristics for circular disk. Farag and Pan [21] assumed that the in-plane vibration mode shapes in circumferential and radial direction were, respectively, expressed in a series summation of trigonometric functions and Bessel functions, and then they solved the in-plane vibration equations for the clamped circular plates. It is well known that differential quadrature method (DQM) as one of point discretization methods can solve this vibration. With this DQM, Eftekhari [22] investigated in-plane vibration of circular arches. In this investigation variable thickness and a moving point load are taken into consideration. Singh and Muhammad [23] solved free in-plane vibration of isotropic rectangular, rhombic, and annular sector plates with the finite element analysis. In their paper, the natural coordinates with multiple points were used to map the geometry. Chan [24] calculated in-plane vibration frequency of circular plate structures by utilizing Hamilton’s principle. Shi et al. [25] calculated the natural in-plane vibration of annular sector plate with using the modified Fourier series method. This method was used, lately, by Wang et al. [26] to analyze the free in-plane vibration of orthotropic sector, annular, and circular plates.
The finite element method (FEM) is one of main numerical approaches to solve structural vibration. Under this method, the basis functions in traditional FEM are employed to describe the unknown structures’ fields. In engineering occasions, however, some drawbacks of FEM may be encountered when obtaining a discretized geometric model, such as time-consuming refinement and cumbersome mesh generation, as well as with geometrical net errors. To overcome the gap between FEM and Computer Aided Design (CAD), recently Hughes et al. [27] introduced an isogeometric analysis (IGA) based on the NURBS basis functions [28]. In the IGA method, the geometry and the variable fields of structures are constructed by using the NURBS basis functions. These functions can preserve exact geometry and provide higher order approximations than Lagrange functions used in standard FEM. What is more, the refinement and meshing process may be implemented without CAD geometry.
In recent decades, the IGA method has been also applied to solve several mechanics problems, such as solid, fluid, contact, and structural mechanics [29–33]. Cho and Ha [34] investigated the exact geometry and enhanced sensitivity of shape design optimization. Wall et al. [35] applied this approach to solve structural shape optimization about two-dimensional problems of linear elasticity. Based on patient-specific vascular NURBS modeling, Zhang et al. [36] used the isogeometric approach to analyze the flow of blood. Bazilevs et al. [37] analyzed the fluid-structure interaction problem about blood flow with the same theory. Kruse et al. [38] analyzed large deformation elasticity by isogeometric collocation methods; then they also solved frictional contact problems. Recently, the IGA was extended to study static and dynamic problems of various kinds of plate and shell structures. Isogeometric Kirchhoff-Love shell formulations were used by Kiendl et al. [39] to solve general hyperelastic materials. This method then was used by Tepole et al. [40] to analyze biological membranes. Benson et al. [41] used IGA method to study Reissner-Mindlin shell. What is more, the IGA method has been utilized in analyzing the vibration of rotation-free shells [42] and functionally graded and multilayered plates with various theories, such as higher order shear deformation [43–45], layerwise deformation [46], and new inverse trigonometric shear deformation theories [47]. Although a large number of studies on the basis of IGA method have been published, to the authors’ knowledge, the works regarding the use of this method for the in-plane vibration of sector, annular, and circular plates are still unavailable in the literature.
In this work, an approximate formulation based on NURBS functions is carried out for studying the in-plane vibration characteristics of the sector, annular, and circular plates, in which the CAD basis functions are utilized for approximating field variables. The principle of virtual work is used to derive the differential equations for free in-plane vibration of sector plate, annular plate, circular sector plate, and circular plate. The convergence study of the formulation is testified in several selected examples. Several numerical results of the natural in-plane frequency parameters for these plates are given. By comparing numerical results with data obtained by the FEM and existing literature, the convergence and reliability of the formulation are tested in several selected examples. In addition, the influence with regard to some effects on the frequency parameters of these plates and their associated mode shapes is also exerted, such as sector angles, the ratios between inner and outer radii, and different boundary conditions.
2. Theoretical Formulations2.1. A Brief of NURBS Functions
A group of nondecreasing values ξi∈R, i=1,2,3⋯m+p+1, are used to define the knot vector E=ξ1,ξ2⋯ξn+p+1. B-spline basis function is piecewise constant when the order of basis functions p is zero. (1)Ni,0ξ=1,ξi≤ξ≤ξi+10,ξi≥ξ,ξ≥ξi+1.For p>0, the basis functions can be calculated recursively as(2)Ni,pξ=ξ-ξiξi+p-ξiNi,p-1ξ+ξi+p+1-ξξi+p+1-ξi+1Ni+1,p-1ξ,where the assumption 0/0=0 is adopted. m and p, respectively, are the total number and polynomial order of basic functions. If there are p+1 same boundary knots in a knot vector, this vector belongs to the type of open knot vector. For the continuity of the basis function, it is easily found that the basis functions are Cp-k continuous at the knot having multiplicity k and Cp-1 continuous at nonrepeated knots. By employing the knot vector E=0,0,0,1,2,3,4,4,5,5,5, a series of quadratic B-spline basis functions can be easily obtained in Figure 1. It is obvious that these basis functions have C0 continuity at the location of ends of the knot vector with repeated knot ξ=0.8 and C1 continuity at the other nonrepeated knots.
Quadratic basis functions for the knot vector E=0,0,0,1,2,3,4,4,5,5,5.
However, it is insufficient to construct exactly arbitrary geometries only by B-spline basis functions, such as ellipse and circle. NURBS basis functions can overcome these limitations by introducing the weights to construct a great variety of geometrical entities. One-dimensional (1D) NURBS basis function of p order can be obtained: (3)Ri,pξ=ωiNi,pξ∑i=1nωiNi,pξ,where ωi is the weight referred to the ith B-spline basis function Ni,p(ξ).
Analogously, the 2D NURBS basis functions of orders p and q can be defined as(4)Ri,jp,qξ,η=ωi,jNi,pξMj,pη∑i=1m∑j=1nωi,jNi,pξMj,qη,where Mj,p(η) is the jth B-spline basis function in the η-direction. ωi,j is the corresponding weight. m and n are the numbers of total basic functions in the ξ- and η-direction, respectively.
Therefore, by giving a set of control points pi,j=xi,j,yi,j and 2D NURBS basis functions Ri,jp,q(ξ,η), a surface is written as(5)Sξ,η=∑i=1m∑j=1nRi,jp,qξ,ηpi,j.
For example, the exact geometries of annular sector plate, annular plate, circular sector plate, and circular plate with their corresponding control nets are shown in Figure 2.
The exact geometries and control nets for sector plates of various types.
Then, by employing the NURBS basis functions, the geometry coordinate (x and y) and the in-plane displacement fields (u and v) of structures can be described as follows: (6)xξ,η=∑a=1n×mRaξ,ηxa,yξ,η=∑a=1n×mRaξ,ηyauξ,η,t=∑a=1n×mRaξ,ηua,vξ,η,t=∑a=1n×mRaξ,ηva,where Ra(ξ,η)=Ri,jp,q(ξ,η) is the NURBS basis function, subscript a=(i-1)·n+j is used for convenience, xa and ya are the control point coordinates corresponding to relational basis function, and (ua, va) represents the displacement vector of the corresponding control points.
2.2. Description of the Models
The investigative geometry models for free in-plane vibration of the sector, annular, and circular plates in natural coordinates are presented in Figure 3. In following analysis, the materials of these sector plates are defined as isotropic with Poisson’s ratio ν, mass density ρ, and Young’s modulus E. The thickness, inner radius, and outer radius of the annular plate sector are h, R0, and R1, respectively. All annular plate, circular sector plate, and circular plate can be obtained from the annular sector plate. It is required that the sector angle ϕ can be set to 2π for the case of annular plate. Similarly, for the circular sector plate, it is required that the inner radius R0 can be set to zero, while for the circular plate the sector angle ϕ and inner radius R0 should be, respectively, set to 2π and zero. However, in the coupling edges ϕ=0, ϕ=2π and inner radius R0 = 0, the kinematic and physical compatibility conditions of the above described annular plate, circular sector plate, and circular plate should be taken into consideration. To satisfy this characteristic, a scheme for coupling coincident control points analogous to coupling nodes in the finite element analysis is used to couple control points and nodes in the model.
The stress-strain relations for the aforementioned isotropic sector, annular, and circular plates undergoing in-plane vibration can be obtained on the basis of generalized Hooke’s law, expressed in the matrix form as(7)σxσyτxy=Q11Q12Q16Q12Q22Q26Q16Q26Q66εxεyγxy,where σx,σy,τxy and εx,εy,γxy, respectively, are the stresses and strains. The elastic stiffness coefficients Qij (i,j=1,2,6) are given as follows: Q11=Q22=E/(1-ν2); Q12=Eν/(1-ν2); Q16=Q26=0; and Q66=E/(2(1+ν)).
Referring to the work of Singh and Muhammad [23], the in- and out-of-plane motions of these plates can be separated; therefore, the motions of these sector plates are assumed without out-of-plane part. Under this assumption the strain-displacement relations are written:(8)εx=∂u∂x,εy=∂v∂y,γxy=∂u∂y+∂v∂x.In this paper, only the in-plane displacement components (u and v), rather than the out-of-plane one (w), are involved in the formulation. Therefore, the current formulation is suitable for the plates composed of the symmetric materials with respect to the mid-plane, such as single-layer isotropic, transverse isotropic, and orthotropic plates, as well as symmetric multilayered composite plates.
2.4. Solution Procedure
The potential (U) energy of sector plate and the kinetic (T) energy of sector plate are written as(9)U=12∬SQ11εx2+2Q12εxεy+Q22εy2+Q66γxy2hdxdy(10)T=12∬Sρhu˙2+v˙2dxdy,where ρ, h, and S are, respectively, the mass density, thickness, and integral surface of sector plate, u˙=∂u/∂t and v˙=∂v/∂t, respectively, represent in-plane velocity field components along the x, y directions. The stress, strain components, and elastic stiffness coefficients are determined by (7) and (8).
By applying Hamilton’s principle [48], the principle of virtual work for sector plate is expressed as (11)δ∫t1t2U+Tdt=0.
Substituting (9)-(10) into (11), the weak statement of in-plane vibration can been obtained as follows: (12)∫t1t2∬SQ11∂u∂x+Q12∂v∂yδ∂u∂x+Q12∂u∂x+Q22∂v∂yδ∂v∂y+Q66∂u∂y+∂v∂xδ∂u∂y+Q66∂u∂y+∂v∂xδ∂v∂x-ρu˙δu˙+v˙δv˙dxdydt=0.
The in-plane vibration governing equations of sector plate can be written by variable separation approach:(13)Q11∂2u∂x2+Q12+Q66∂2v∂x∂y+Q66∂2u∂y2=ρu¨Q22∂2v∂y2+Q12+Q66∂2u∂x∂y+Q66∂2v∂x2=ρv¨.And the following definitions of stress resultants are used:(14)N-x=Q11∂u∂x+Q12∂v∂y,N-xy=Q66∂u∂y+∂v∂xat r=R0,R1N-y=Q12∂u∂x+Q22∂v∂y,N-yx=Q66∂u∂y+∂v∂xat ϕ=0,ϕ1,where N-x and N-xy are stress resultants in the x- and y-axis directions in the radial edges, respectively. Similarly, N-y and N-yx are stress resultants in the y- and x-axis directions in circumferential edges.
The boundary conditions are of the following form (see [49]): (15)N-x or uN-xy or vN-y or vN-yx or u.According to the study of Cottrell [50], the IGA is analogous to the traditional FEM for using the homogeneous Dirichlet boundary conditions on u, v. The following boundary conditions could be used in this paper:
As for the center of circular and circular sector plates and the circumferential edges of annular and circular plates, it is essential to consider their physical compatibility and kinematic conditions. Therefore, similar to standard FEM for coupling coincident nodes, IGA method can couple relevant control points. From Figure 2, there are six control points in coupling edges ϕ=0 of annular plate and four control points in the center of circular sector plate in original one. After coupling relevant control points, there are three control points and one point in coupling edges ϕ=0 of annular plate and the center of circular sector plate, respectively.
The NURBS functions substitute for interpolation functions which are utilized in FEM to construct the unknown fields. A typical element for two-dimensional structures is used in the isogeometric analysis. There are present space for calculating integral equation and physical space in conventional finite element implementations, only. The parametric space is unnecessary. Figure 4 shows the parent, parametric, and physical spaces and mapping variations. The current space is mapped to geometrical space by Xe=S∘ϕ~, where mappings ϕ~:Ω~→Ω^ and S:Ω~→Ω, respectively, are from current space Ω~ to parametric space Ω~ and from parametric space Ω~ to geometrical space Ω.
The mappings among parent space, parametric space, and physical space.
An example is presented to explain the mapping ϕ~:Ω~→Ω^ for the Ω^e=[ξi,ξi+1]⊗[ηj,ηj+1] parametric space of an element(18)ϕ~eξη=12ξi+1-ξiξ+12ξi+1+ξi12ηj+1-ηjη+12ηj+1+ηjand the corresponding Jacobian determinant(19)J~=14ξi+1-ξiηj+1-ηj.
According to (6), it is easy to calculate the mapping S:Ω^→Ω. And the Jacobian of the relevant transformation is given: (20)J^=∂x∂ξ∂x∂η∂y∂ξ∂y∂η.By using the first derivatives of NURBS basis functions, the Jacobian may be represented by the following expression:(21)J^=∑ane∂Ra∂ξxae∑ine∂Ra∂ηxae∑ane∂Ra∂ξyae∑ane∂Ra∂ηyae,where ne is the number of total control points of each element. Therefore, the Jacobian |J^| can be determined by (21).
Considering the above Jacobian determinant |J^| and |J~|, the Jacobian determinant can be written as(22)J=J^·J~.
In order to facilitate the analysis, the transformation formula which is used for interpretation of a quadrature rule from geometrical space to current space is written as(23)∬Sfx,ydxdy=∑eNel∬Ωefx,ydxdy=∑eNel∬Ω~efξ~,η~Jdξ~dη~,where Nel is the total number of elements. It is easy to see that transformation formula can be obtained with Jacobian determinant and the mappings in Figure 4.
By using the separation of variables u(x,t)=u~(x)ejωlt and v(x,t)=v~(x)ejωlt, the principle of virtual work for sector plate in (10) can be written as (24)U=12∬ShQ11∂u~∂x2+2Q12∂v~∂y∂u~∂x+Q22∂v~∂y2+Q66∂u~∂y+∂v~∂x2dxdyT=-12∬Sωl2ρhu~2+v~2dxdy,where ωl is the lth free vibration frequency of the sector plate.
The Lagrange function of sector plate can be written as(25)L=12∬ShQ11∂u~∂x2+2Q12∂v~∂y∂u~∂x+Q22∂v~∂y2+Q66∂u~∂y+∂v~∂x2-ωl2ρu~2+v~2dxdy.The variables and geometry coordinate of control points in (6) can be written as (26)xξ,η=∑a=1n×mRaξ,ηxa,yξ,η=∑a=1n×mRaξ,ηya,u~ξ,η=∑a=1n×mRaξ,ηu~a,vξ,η=∑a=1n×mRaξ,ηv~a.
Similarly, the geometry coordinate and the in-plane displacement fields of each element can be written as(27)xeξ,η=∑a=1neRaξ,ηxae,yξ,η=∑a=1neRaξ,ηyae,u~eξ,η=∑a=1neRaξ,ηuae,v~eξ,η=∑a=1neRaξ,ηv~ae,where ne is the number of total control points of each element; the superscript e denotes use for each element.
By substituting (26) into (25) and then applying Rayleigh-Ritz procedure to solve every u~a and v~a, (28)δL=0.
According to the element assembly operator in the standard FEM and quadrature rule in (23), (27) can be also written as (29)δ12∑e=1Nel∬ΩeQ11∂u~e∂x2+2Q12∂v~e∂y∂u~e∂x+Q22∂v~e∂y2+Q66∂u~e∂y+∂v~e∂x2-ωl2ρu~e2+v~e2dxdy=0.Therefore, the in-plane vibration governing expression can be written as(30)K-ωl2Mφ~l=0,where M and K are, respectively, the global mass and stiffness matrices for sector plate. ωl denotes the lth natural frequency of sector plate, and φ~l=u~,v~ denotes the corresponding eigenvector nodal displacements, in which displacement component vectors u~ and v~ are written as (31)u~=u~1,u~2,u~3⋯u~n·mTv~=v~1,v~2,v~3⋯v~n·mT.The eigenproblem in (30) can be solved by using Arnoldi Method [51, 52], and the mode shapes can be expressed by multiplying the NURBS basis functions with eigenvector.
To calculate the global stiffness matrix and mass matrix in (30), each local element mass matrix and stiffness matrix should be computed according to (29). Therefore, the local element stiffness matrix can be written as(32)Ke=K11K12K21K22.The explicit forms of the local element stiffness matrix Ke are given by (33)Kij11=∬ΩeQ11∂Ri∂x∂Rj∂x+Q66∂Ri∂y∂Rj∂ydxdyKij12=∬ΩeQ12∂Ri∂x∂Rj∂y+Q66∂Ri∂y∂Rj∂xdxdyKij21=∬ΩeQ12∂Ri∂y∂Rj∂x+Q66∂Ri∂x∂Rj∂ydxdyKij22=∬ΩeQ22∂Ri∂y∂Rj∂y+Q66∂Ri∂x∂Rj∂xdxdy.Each local element mass matrix for sector plates can be written as(34)Me=M11M12M21M22.The explicit forms of the local element mass matrix Me are given by (35)Mij11=∬ΩeρRiRjdxdyMij12=Mij21=0Mij22=∬ΩeρRiRjdxdy.By using integral transform which is determined in (23), the integral region of element can be transformed to numerical integral region. Then numerical integration can be obtained by applying the standard Gauss-Legendre quadrature.
The global stiffness and mass matrices of annular sector plate can be obtained by using (29)–(35). For the stiffness and mass matrices of the annular plate, circular sector plate, and circular plate, a coupling coincident control points scheme mentioned above is used. When the angle ϕ of the annular sector plate is equal to 2π and with coupling the edges ϕ=0, ϕ=2π, the annular sector plate can be regarded as an annular plate. When inner radius R0 of the annular sector plate is zero and with coupling the control points in the inner radius R0, the annular sector plate can be regarded as a circular sector plate. In a similar way, circular plate can be obtained when the angle ϕ of the circular sector plate is equal to 2π and with coupling the edges ϕ=0, ϕ=2π.
2.5. Refinement Algorithms
The accuracy of FEM computational solutions is greatly related to the number of elements and the order of functions used to describe unknown field variables. To capture fast convergence in the solution field, it is essential that the discretization of intricate geometries can be successively refined. The basic B-spline refinement mechanisms are different from the counterparts of the FEM. One of the main advantages of a discretization tool based on the B-spline basis functions is the ability to apply different refinements algorithms in a simple manner. Three types of refinement algorithms [53] are generally used in the IGA.
In the IGA, knot insertion is a strategy that inserts additional knots into the knot vectors. In this subdivision strategy, both the number and control points of element are changed but the order of basis functions is maintained. Degree elevation is a strategy that increases the order of B-spline basis functions. In this strategy, both the order and control points of basis B-spline functions are changed yet knot span is not increased. The knot insertion and degree elevation, respectively, are analogous to the h- and p-refinement of the standard FEA. In addition, the combination of knot insertion with degree elevation is the ability to refine discretization. A strategy that first elevates the order of NURBS basis functions and then inserted a unique knot into a nonzero knot span is said to be k-refinement. As an efficient and robust scheme, this strategy can increase the number of elements and provide higher order of basis functions.
3. Numerical Results and Discussions
In this part, the convergence of the current method is first tested. Then, several selected numerical examples of the in-plane vibration analysis for sector, annular, and circular plates are given to verify the validity and accuracy of this IGA method. Free (F) and clamped (C) boundary conditions are considered. For instance, CCCC and FFFF, respectively, represent all clamped and free sides, CCCF means the fourth side is free and other three sides are fixed. The sequence of boundary edges for the sector plates is defined in Figure 3. Then, some free in-plane vibration mode shapes of sector plates are plotted and relevant frequencies are presented. Finally, the effects of some parameters study examples are also given. In the following analysis, the isotropic annular, circular, and sector plates are considered as thin plates. The nondimensional frequency parameter Ω=2ωR1/πρ/E1-ν2 is used if no special statement is given; the material parameters of the different types of annular sector, annular, circular sector, and circular plates are E = 70 Gpa, ν = 0.3, and ρ = 2700 kg/m3.
3.1. Convergence and Validation
The annular sector plate as an example is applied to test the convergence of the present approach for free in-plane natural frequencies. For conciseness and convenience, three refinement schemes such as h-, p-, and k-refinements which elevate the polynomial order and insert knots with preserving the geometry and parameterization are employed. Then, the first eight nondimensional frequency parameters of the CCCC annular sector plate by using different orders and elements are shown in Table 1. The plates are characterized by the following geometry constants: R1 = 1 m, R0/R1 = 0.5, and ϕ=π/4. Obviously, the computed results in Table 1 show a great convergence trend.
Convergence of frequency parameters Ω=2ωR1/πρ/E1-ν2Ω of plates with CCCC boundaries and inner-to-outer radii ratios: R0/R1=0.5.
Order
Element
Mode number
1
2
3
4
5
6
7
8
p=q=2
8×8
4.4070
4.8359
5.5223
6.9141
7.0940
7.3009
8.0867
8.7015
16×16
4.4065
4.8353
5.5204
6.9101
7.0884
7.2908
8.0748
8.6908
32×32
4.4064
4.8352
5.5203
6.9098
7.0880
7.2902
8.0742
8.6901
p=q=3
8×8
4.4066
4.8355
5.5204
6.9101
7.0883
7.2908
8.0748
8.6908
16×16
4.4065
4.8353
5.5203
6.9099
7.0881
7.2903
8.0742
8.6902
32×32
4.4064
4.8352
5.5203
6.9098
7.0880
7.2902
8.0741
8.6901
p=q=4
8×8
4.4066
4.8355
5.5204
6.9099
7.0882
7.2905
8.0744
8.6903
16×16
4.4065
4.8354
5.5204
6.9100
7.0881
7.2903
8.0743
8.6903
32×32
4.4065
4.8353
5.5203
6.9099
7.0880
7.2902
8.0742
8.6902
p=q=5
8×8
4.4067
4.8356
5.5205
6.9099
7.0882
7.2905
8.0744
8.6903
16×16
4.4066
4.8354
5.5204
6.9100
7.0882
7.2904
8.0743
8.6903
32×32
4.4065
4.8353
5.5203
6.9099
7.0881
7.2903
8.0742
8.6902
The comparisons of frequency parameters of these plates including annular sector plate, annular plate, and circular sector plate are shown in Tables 2–4. The first six dimensionless natural frequencies Ω=ωR1ρ/E of annular sector plates are presented Table 2. In this table, the inner-to-outer radii ratios R0/R1=0.5, the sector angles ϕ=30°,90°, and different boundary conditions such as CCCC, CFCF, FFFF, and CCCF are taken into account. Some results in available literature reported by Singh and Muhammad [23] using traditional FEM and Shi et al. [25] using Ritz method are used in Table 2 to validate the accuracy of results of annular sector plate with the present method. From Table 2, a good agreement can be obtained by observing these results. The discrepancies of the results can be negligible.
Comparison of frequency parameters Ω=ωR1ρ/E for annular sector plates with various boundaries and inner-to-outer radii ratios: R0/R1=0.5.
ϕ (°)
Mode
CCCC
CFCF
FFFF
CCCF
Ref. [23]
Ref. [25]
Present
Ref. [23]
Ref. [25]
Present
Ref. [25]
Present
Ref. [25]
Present
30
1
8.1592
8.1500
8.1527
3.4489
3.4788
3.4774
4.8589
4.8777
5.3157
5.3112
2
8.8379
8.8085
8.8122
6.2770
6.2179
6.2184
5.6760
5.6938
6.8500
6.8399
3
10.1089
10.2116
10.2138
6.4037
6.4445
6.4441
5.7026
5.7199
7.4560
7.5122
4
12.7648
12.7775
12.7822
8.3676
8.3194
8.3189
6.1162
6.1370
10.2296
10.2655
5
13.2623
13.3090
13.3143
8.3358
8.4499
8.4478
7.6716
7.6990
10.7897
10.8673
6
13.2555
13.3195
13.3225
9.4104
9.4264
9.4264
8.5692
8.6995
11.5209
11.5250
90
1
4.7576
4.7508
4.7516
3.7834
3.8017
3.7986
1.5633
1.5641
4.0055
4.0121
2
6.3458
6.3216
6.3233
4.4797
4.4676
4.4636
2.6785
2.6702
5.4346
5.3987
3
6.9512
6.8564
6.8582
6.0711
5.9866
5.9864
2.9795
2.9792
6.5046
6.4447
4
6.9793
6.8921
6.8937
6.5393
6.5032
6.4992
4.4513
4.4293
6.6321
6.5915
5
7.5749
7.5573
7.5588
6.5952
6.5497
6.5504
4.4627
4.4754
7.0237
6.9826
6
8.0831
8.0973
8.0983
6.7786
6.7742
6.7739
4.6668
4.6575
7.3563
7.3541
Comparison of frequency parameters Ω=2ωR1/πρ/E1-ν2 for annular plates with various boundary conditions.
R0/R1
Mode
FF
FC
CC
Ref. [19]
Ref. [20]
Present
Ref. [19]
Ref. [20]
Present
Ref. [19]
Ref. [20]
Present
0.2
n=1
1.156
1.155
1.156
1.472
1.473
1.472
1.947
1.963
1.947
2.688
2.687
2.680
2.311
2.313
2.310
2.840
2.870
2.840
n=2
0.777
0.777
0.777
1.786
1.788
1.472
2.363
2.374
2.363
1.681
1.680
1.683
2.762
2.765
2.762
3.050
3.066
3.050
n=3
1.449
1.450
1.449
2.580
2.584
2.580
2.845
2.857
2.844
2.379
2.379
2.380
3.399
3.364
3.399
3.571
3.582
3.570
n=4
1.936
1.935
1.936
3.296
3.301
3.296
3.359
3.366
3.358
3.070
3.069
3.067
4.123
4.130
4.123
4.200
4.051
4.199
0.4
n=1
1.177
1.177
1.177
1.761
1.764
1.761
2.399
2.418
2.399
2.829
2.829
2.828
2.454
2.458
2.454
3.712
3.743
3.713
n=2
0.504
0.504
0.505
1.904
1.913
1.903
2.814
2.831
2.814
1.715
1.714
1.714
2.901
2.905
2.901
3.715
3.741
3.716
n=3
1.132
1.133
1.132
2.248
2.252
2.249
3.293
3.314
3.293
2.341
2.340
2.272
3.497
3.501
3.496
3.931
3.953
3.931
n=4
1.736
1.736
1.736
2.767
2.770
2.767
3.699
3.750
3.699
2.957
2.956
2.953
4.109
4.115
4.109
4.400
4.416
4.399
Comparison of the first nine frequency parameters Ω=2ωR1/πρ/E1-ν2 for circular sector plate with various boundary conditions.
BC
ϕ (°)
Method
Mode sequence number
1
2
3
4
5
6
7
8
9
CCC
60
ANSYS
3.4710
3.5239
4.2842
5.3930
5.5591
5.8322
5.8813
5.9433
7.0236
Present
3.4716
3.5246
4.2847
5.3941
5.5601
5.8331
5.8821
5.9443
7.0243
90
ANSYS
2.7035
2.9895
3.4421
4.1459
4.5198
4.5668
4.9401
5.1551
5.4300
Present
2.7039
2.9900
3.4423
4.1466
4.5202
4.5676
4.9407
5.1560
5.4306
180
ANSYS
1.8735
2.3630
2.5835
2.8205
3.1101
3.1957
3.6007
3.7763
3.8302
Present
1.8737
2.3632
2.5836
2.8207
3.1103
3.1961
3.6010
3.7768
3.8305
270
ANSYS
1.6512
1.9271
2.2385
2.5417
2.6073
2.6320
2.9348
3.0282
3.2372
Present
1.6512
1.9274
2.2387
2.5419
2.6074
2.6324
2.9351
3.0285
3.2375
CFC
60
ANSYS
2.0725
2.3418
2.9577
3.9833
4.4296
4.4502
4.7109
4.9565
5.2712
Present
2.0727
2.3420
2.9579
3.9835
4.4303
4.4507
4.7113
4.9568
5.2713
90
ANSYS
1.5530
1.6721
2.1934
2.7998
3.2069
3.4231
3.8270
4.0418
4.0815
Present
1.5532
1.6723
2.1934
2.8000
3.2073
3.4234
3.8274
4.0421
4.0817
180
ANSYS
0.7003
1.1850
1.4080
1.6112
1.8608
1.9677
2.4200
2.4398
2.6848
Present
0.7003
1.1851
1.4081
1.6112
1.8609
1.9677
2.4201
2.4400
2.6849
270
ANSYS
0.3983
0.7898
1.0803
1.3769
1.4061
1.4591
1.7509
1.7527
2.0124
Present
0.3983
0.7898
1.0803
1.3770
1.4061
1.4592
1.7510
1.7527
2.0125
FCF
60
ANSYS
0.8416
1.5705
1.6204
2.6790
3.1841
3.5447
3.7060
3.8662
4.0591
Present
0.8417
1.5707
1.6204
2.6791
3.1844
3.5449
3.7061
3.8663
4.0594
90
ANSYS
1.0527
1.5126
1.7130
2.5865
2.8986
3.1485
3.1677
3.6040
4.0099
Present
1.0528
1.5127
1.7131
2.5867
2.8987
3.1486
3.1679
3.6041
4.0100
180
ANSYS
1.3000
1.3378
1.6175
2.0817
2.3425
2.4921
2.6675
3.1048
3.1080
Present
1.3001
1.3379
1.6176
2.0818
2.3427
2.4922
2.6676
3.1050
3.1081
270
ANSYS
1.2704
1.4099
1.4850
1.7801
2.1763
2.2147
2.2991
2.5909
2.7198
Present
1.2705
1.4099
1.4850
1.7799
2.1765
2.2146
2.2992
2.5910
2.7200
Besides, annular and circular sector plates under this investigation are considered as special cases on the basis of the model of annular sector plate. As mentioned above, they are obtained by using a scheme for coupling coincident control points. Therefore, the results of frequency parameters of these plates should be investigated. The nondimensional frequency parameters results of the annular plate are presented in Table 3. The inner-to-outer radii ratios used for this table are 0.2 and 0.4. Different boundary conditions FF, FC, and CC are taken into consideration. Similarly, the natural frequencies of annular plate are compared with the results which are taken from Irie et al. [19] using transfer matrix method and Bashmal et al. [20] using Rayleigh-Ritz method. Table 4 presents the first nondimensional nine frequency parameters of the circular sector plate. The sector angles ϕ=30°,90° and boundary conditions CCC, CFC, and FCF are taken into consideration. To validate these results of the present method in Table 4, the natural frequency parameters computed by ANSYS based on the FEM are utilized. In this ANSYS computation, the type and the size of element are, respectively, shell82 and 0.02, and mapped meshing grid is adapted. Based on the above description of tables, it is obvious that this IGA method has a great accuracy and efficiency to evaluate the natural in-plane vibration frequencies for the annular sector, annular, and circular sector plates.
Figures 5–7 depict first four mode shapes and relevant frequency parameters of the free in-plane vibration for annular sector plate with different geometric parameters and boundary conditions. From these figures, it can be easily seen that in-plane vibration mode shapes are more complicated than their flexural vibration modes even though they belong to lower-order modes. For example, the deformation of extension and compression in one region can suddenly change into a shear state in another region as mode varies. One of the reasons is maybe that the in-plane vibration frequency is determined by many factors rather than individual variable of the plate. The deformation of the first mode extends and compresses in the circumferential direction under CCCC and CFCF boundary conditions but in the radial direction under FCFC boundary condition. The variations of mode shapes with sector angle also are obviously seen. Therefore, the boundary conditions and the geometry parameters play significant roles in their in-plane vibration mode shapes. Figures 8–10 present some mode shapes, corresponding to frequency parameters of annular, circular, and circular plates with different geometric parameters, respectively.
The first four mode shapes of CCCC annular sector plate (R0/R1=0.5) with different sector angles, that is, ϕ = 90°, 180°, and 300°.
The first four mode shapes of CFCF annular sector plate (R0/R1=0.5) with different sector angles, that is, ϕ = 90°, 180°, and 300°.
The first four mode shapes of FCFC annular sector plate (R0/R1=0.5) with different sector angles, that is, ϕ = 90°, 180°, and 300°.
The first seven mode shapes of CF annular plate with different inner-to-outer radii ratios, that is, R0/R1 = 0.2, 0.5, and 0.8.
The first seven mode shapes of CFC circular sector plate with different sector angles, that is, ϕ = 90°, 180°, and 300°.
The first four mode shapes of C circular plate with different boundary conditions.
Ω1 = 1.3692
Ω2 = Ω3 = 1.5859
Ω4 = Ω5 = 2.1321
Ω6 = Ω7 = 2.2237
3.2. The Effects of Parameters
In numerical studies, the effects of several key parameters are investigated in Figures 11–13. The order p=q=3 and the element 32 × 32 are used to calculate the frequency parameters. The change trends of the first three frequencies for annular sector plate and circular sector plate are shown in Figure 11. In this analysis, the inner-to-outer radii ratios R0/R1=0.5 for annular sector plate and R0/R1=0 for circular sector plate are used, and these sector annular angles gradually increase from 2° to 358°. For sector annular plates with CCCC, FCFC, and CCCF boundary conditions and circular sector plates with CCC, CFC, and CCF boundary conditions, the conclusion can be clearly obtained that the frequency parameters of the plate are in the monotonically decreasing state when the sector angle increases. The frequency parameters sharply decrease before the annular sector angle is small, that is, ϕ ≤ 10°, and the descending velocities of the frequency parameters become slower when the sector angle varies from 10° to 100°, and the frequency parameters remain stable after the annular sector angle ϕ ≥ 100°. For the cases of CFCF annular sector plate and FCF circular sector plate, the frequency parameters increase rapidly with beginning and then decease or increase slightly with sector angle increasing in general. The effects of inner-to-outer radii ratios on first three natural frequencies of annular sector plate and annular plate are illustrated in Figure 12. The sector angle 60° is considered for annular sector plate. For the case of sector annular plate under boundary conditions CCCC, CFCF, and CCCF and annular plates under boundary conditions CC, CF, and FC, the frequency parameters always increase and the ascending velocities of frequency parameters gradually become faster as the inner-to-outer radii ratios increase. For the FCFC sector annular plate, the frequency parameters always decrease. Figure 13 presents the variational rules of first three frequency parameters for the annular sector plate by changing the inner-to-outer radii ratios. The various sector angles and boundary conditions are considered in this analysis. It is obvious that the changes of natural frequencies for CCCC, CFCF, and CCCF boundary conditions also increase as the inner-to-outer radii ratios increase. But for FCFC boundary condition the variations of frequency parameter are more complicated; when the sector angle is larger than 10°, the first three frequency parameters decrease with the increasing inner-to-outer radii ratios. When the sector angle is smaller 10°, the variations of frequency parameters are different. For instance, the first frequency parameter monotonically decreases, the second frequency parameter first decreases and then increases, and there are almost few changes in the third frequency parameters.
The variations of some frequency parameters for sector plates with different boundary conditions and sector angles: (a) annular sector plate; (b) circular sector plate.
The variations of some frequency parameters for plates with different boundary conditions and inner-to-outer radii ratios: (a) annular sector plate; (b) annular plate.
The variation of first three frequency parameters for the annular sector plate with different boundary conditions, inner-to-outer radii ratios, and sector angles.
4. Conclusions
In this paper, the free in-plane vibration of annular sector, annular, circular sector, and circular plates is investigated by an effective approximate formulation based on isogeometric finite element analysis. By using the NURBS basic functions in this method, the polynomial order can be easily evaluated and the desired geometry can be obtained. The gap between CAD and FEA can be possibly narrowed by applying the identical functions to construct the unknown filed and the geometry. The in-plane vibration differential equations for sector plate, annular plate, and circular plate can be systematically derived by the principle of virtual work. Three refinement schemes (the h-, p-, and k-refinement strategies) are used to discretize the geometry and refine element, and the efficiency and rapid convergence of the k-refinement in IGA have been confirmed. Then several numerical results of the natural nondimensional frequencies for sector, annular, and circular plates are computed and compared with relevant data in open references or from finite element method. It can be easily obtained that the current formulation provides a reliability means for the in-plane vibration of these plates. Additionally, some new natural nondimensional frequency parameters and corresponding mode shapes for annular plate, circular plate, and sector plates are presented. The new calculated results can be adopted by other researchers as a benchmark for future researches. The effects of sector angle and ratio of inside to outside radii for these plates with different boundary conditions are investigated as well.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This paper is supported by the National Natural Science Foundation of China (nos. 51775125 and 51279035) and the Fundamental Research Funds for the Central Universities of China (no. HEUCFG201713).
YeT.JinG.SuZ.JiaX.A unified Chebyshev–Ritz formulation for vibration analysis of composite laminated deep open shells with arbitrary boundary conditions201484444147110.1007/s00419-013-0810-1JinG.SuZ.YeT.JiaX.Three-dimensional vibration analysis of isotropic and orthotropic conical shells with elastic boundary restraints2014892072212-s2.0-8490755400110.1016/j.ijmecsci.2014.09.005JinG. Y.YeT.MaX. L.ChenY.SuZ.XieX.A unified approach for the vibration analysis of moderately thick composite laminated cylindrical shells with arbitrary boundary conditions20137535737610.1016/j.ijmecsci.2013.08.0032-s2.0-84884655860YeT.JinG.SuZ.ChenY.A modified Fourier solution for vibration analysis of moderately thick laminated plates with general boundary restraints and internal line supports201480294610.1016/j.ijmecsci.2014.01.0012-s2.0-84893352855YeT.JinG.SuZ.Three-dimensional vibration analysis of laminated functionally graded spherical shells with general boundary conditions2014116157158810.1016/j.compstruct.2014.05.0462-s2.0-84904557880XingY. F.LiuB.Exact solutions for the free in-plane vibrations of rectangular plates20095132462552-s2.0-6154909528010.1016/j.ijmecsci.2008.12.009Zbl1264.74108BercinA. N.LangleyR. S.Application of the dynamic stiffness technique to the in-plane vibrations of plate structures19965958698752-s2.0-003056766410.1016/0045-7949(95)00308-8Nefovska-DanilovicM.PetronijevicM.In-plane free vibration and response analysis of isotropic rectangular plates using the dynamic stiffness method201515282952-s2.0-8492424116310.1016/j.compstruc.2015.02.001GormanD. J.Free in-plane vibration analysis of rectangular plates by the method of superposition20042723-58318512-s2.0-214276519310.1016/S0022-460X(03)00421-8GormanD. J.Exact solutions for the free in-plane vibration of rectangular plates with two opposite edges simply supported20062941-21311612-s2.0-3364649576210.1016/j.jsv.2005.10.023AndrianovI. V.AwrejcewiczJ.ChernetskyyV.Analysis of natural in-plane vibration of rectangular plates using homotopy perturbation approach2006200682059810.1155/MPE/2006/20598MR2278500ChenY.JinG.LiuZ.Flexural and in-plane vibration analysis of elastically restrained thin rectangular plate with cutout using Chebyshev-Lagrangian method2014892642782-s2.0-8490756318310.1016/j.ijmecsci.2014.09.006MohazzabA. H.DozioL.A spectral collocation solution for in-plane eigenvalue analysis of skew plates201594-951992102-s2.0-8492584274910.1016/j.ijmecsci.2015.03.008ZhouY.WangQ.ShiD.LiangQ.ZhangZ.Exact solutions for the free in-plane vibrations of rectangular plates with arbitrary boundary conditions201713011010.1016/j.ijmecsci.2017.06.004ShiZ.YaoX.PangF.WangQ.A semi-analytical solution for in-plane free vibration analysis of functionally graded carbon nanotube reinforced composite circular arches with elastic restraints20171824204342-s2.0-8503012983210.1016/j.compstruct.2017.09.045RostamiH.Rahbar RanjiA.Bakhtiari-NejadF.Free in-plane vibration analysis of rotating rectangular orthotropic cantilever plates2016115-1164384562-s2.0-8498994891610.1016/j.ijmecsci.2016.07.030LiH.PangF.WangX.LiS.Benchmark Solution for Free Vibration of Moderately Thick Functionally Graded Sandwich Sector Plates on Two-Parameter Elastic Foundation with General Boundary Conditions2017201735401862910.1155/2017/40186292-s2.0-85030530368HuangD. T.Effects of constraint, circular cutout and in-plane loading on vibration of rectangular plates2013681141242-s2.0-8487521561210.1016/j.ijmecsci.2013.01.005IrieT.YamadaG.MuramotoY.Natural frequencies of in-plane vibration of annular plates19849711711752-s2.0-024263810910.1016/0022-460X(84)90479-6BashmalS.BhatR.RakhejaS.In-plane free vibration of circular annular disks20093221-22162262-s2.0-6174908852610.1016/j.jsv.2008.11.024FaragN. H.PanJ.Modal characteristics of in-plane vibration of circular plates clamped at the outer edge20031134 I193519462-s2.0-003738119210.1121/1.1553456EftekhariS. A.Differential quadrature procedure for in-plane vibration analysis of variable thickness circular arches traversed by a moving point load2016407-84640466310.1016/j.apm.2015.11.046MR3478235SinghA. V.MuhammadT.Free in-plane vibration of isotropic non-rectangular plates20042731-22192312-s2.0-234255729710.1016/S0022-460X(03)00496-6ChanI. P.Frequency equation for the in-plane vibration of a clamped circular plate20083131-232533310.1016/j.jsv.2007.11.0342-s2.0-40949115127ShiX.LiW.ShiD.Free in-plane vibrations of annular sector plates with elastic boundary supportsProceedings of the 164th Meeting of the Acoustical Society of America 2012October 201210.1121/1.47738882-s2.0-84878957871WangQ.ShiD.LiangQ.e AhadF.A unified solution for free in-plane vibration of orthotropic circular, annular and sector plates with general boundary conditions20164021-229228925310.1016/j.apm.2016.06.005MR3549744HughesT. J.CottrellJ. A.BazilevsY.Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement200519439-414135419510.1016/j.cma.2004.10.008MR2152382PieglL. A.TillerW.1995Berlin, GermanySpringerZbl0828.68118CottrellJ. A.RealiA.BazilevsY.HughesT. .Isogeometric analysis of structural vibrations200619541-435257529610.1016/j.cma.2005.09.027MR2243311Zbl1119.740242-s2.0-33746838484YinS.HaleJ. S.YuT.BuiT. Q.BordasS. P. A.Isogeometric locking-free plate element: A simple first order shear deformation theory for functionally graded plates201411811211382-s2.0-8491773691010.1016/j.compstruct.2014.07.028HughesT. J.RealiA.SangalliG.Efficient quadrature for NURBS-based isogeometric analysis20101995-830131310.1016/j.cma.2008.12.004MR2576763KiendlJ.BletzingerK.-U.LinhardJ.WüchnerR.Isogeometric shell analysis with Kirchhoff-Love elements200919849-523902391410.1016/j.cma.2009.08.013MR2557478ShojaeeS.ValizadehN.IzadpanahE.BuiT.VuT.-V.Free vibration and buckling analysis of laminated composite plates using the NURBS-based isogeometric finite element method2012945167716932-s2.0-8485869844010.1016/j.compstruct.2012.01.012ChoS.HaS.-H.Isogeometric shape design optimization: exact geometry and enhanced sensitivity2009381537010.1007/s00158-008-0266-zMR2475898Zbl1274.742212-s2.0-59449101918WallW. A.FrenzelM. A.CyronC.Isogeometric structural shape optimization200819733-402976298810.1016/j.cma.2008.01.025MR24270962-s2.0-44749089360ZhangY.BazilevsY.GoswamiS.BajajC. L.HughesT. J.Patient-specific vascular NURBS modeling for isogeometric analysis of blood flow200719629-302943295910.1016/j.cma.2007.02.009MR2325400BazilevsY.CaloV. M.ZhangY.HughesT. J. R.Isogeometric fluid-structure interaction analysis with applications to arterial blood flow2006384-53103222-s2.0-3374567781310.1007/s00466-006-0084-3Zbl1161.74020KruseR.Nguyen-ThanhN.De LorenzisL.HughesT. J.Isogeometric collocation for large deformation elasticity and frictional contact problems20152967311210.1016/j.cma.2015.07.022MR33962922-s2.0-84939617131KiendlJ.HsuM.-C.WuM. C.RealiA.Isogeometric Kirchhoff-Love shell formulations for general hyperelastic materials201529128030310.1016/j.cma.2015.03.010MR33452512-s2.0-84928728099TepoleA. B.KabariaH.BletzingerK.-U.KuhlE.Isogeometric Kirchhoff-Love shell formulations for biological membranes201529332834710.1016/j.cma.2015.05.006MR33959112-s2.0-84930959329BensonD. J.BazilevsY.HsuM. C.HughesT. .Isogeometric shell analysis: the Reissner-Mindlin shell20101995-827628910.1016/j.cma.2009.05.011MR2576761Zbl1227.741072-s2.0-72549097830BensonD. J.BazilevsY.HsuM.-C.HughesT. .A large deformation, rotation-free, isogeometric shell201120013-161367137810.1016/j.cma.2010.12.003MR2774751Zbl1228.740772-s2.0-78751648559TranL. V.FerreiraA. J. M.Nguyen-XuanH.Isogeometric analysis of functionally graded plates using higher-order shear deformation theory201351436838310.1016/j.compositesb.2013.02.0452-s2.0-84877805644TranL. V.LyH. A.LeeJ.WahabM. A.Nguyen-XuanH.Vibration analysis of cracked FGM plates using higher-order shear deformation theory and extended isogeometric approach201596-97657810.1016/j.ijmecsci.2015.03.0032-s2.0-84927601982Nguyen-XuanH.TranL. V.ThaiC. H.KulasegaramS.BordasS. P. A.Isogeometric analysis of functionally graded plates using a refined plate theory2014642222342-s2.0-8490164886410.1016/j.compositesb.2014.04.001ThaiC. H.FerreiraA. J. M.CarreraE.Nguyen-XuanH.Isogeometric analysis of laminated composite and sandwich plates using a layerwise deformation theory201310419621410.1016/j.compstruct.2013.04.0022-s2.0-84879448094ThaiC. H.FerreiraA.BordasS.RabczukT.Nguyen-XuanH.Isogeometric analysis of laminated composite and sandwich plates using a new inverse trigonometric shear deformation theory20144318910810.1016/j.euromechsol.2013.09.001ReddyJ. N.1984New York, NY, USAJohn Wiley & Sons Inc.ReddyJ. N.A simple higher order theory for laminated composite plates198451474575210.1115/1.3167719Zbl0549.730622-s2.0-0021621010CottrellJ. A.HughesT. I. R.BazilevsY.2009New York, NY, USAJohn Wiley & Sons Inc.LehoucqR. B.SorensenD. C.Deflation techniques for an implicitly restarted Arnoldi iteration199617478982110.1137/S0895479895281484MR1410702Zbl0863.650162-s2.0-0030502142SorensenD. C.Implicit application of polynomial filters in a k-step Arnoldi method199213135738510.1137/0613025MR1146670CottrellJ. A.HughesT. J. R.RealiA.Studies of refinement and continuity in isogeometric structural analysis200719641-44416041832-s2.0-3454762976710.1016/j.cma.2007.04.007Zbl1173.74407