Buckling and free vibration analysis of thick rectangular plates resting on elastic foundation using mixed finite element and differential quadrature method

In this article, a combination of the finite element (FE) and differential quadrature (DQ) methods is used to solve the eigenvalue (buckling and free vibration) equations of rectangular thick plates resting on elastic foundations. The elastic foundation is described by the Pasternak (two-parameter) model. The three dimensional, linear and small strain theory of elasticity and energy principle are employed to derive the governing equations. The in-plane domain is discretized using two dimensional finite elements. The spatial derivatives of equations in the thickness direction are discretized in strong-form using DQM. Buckling and free vibration of rectangular thick plates of various thicknesses to width and aspect ratios with Pasternak elastic foundation are investigated using the proposed FE-DQ method. The results obtained by the mixed method have been verified by the few analytical solutions in the literature. It is concluded that the mixed FE-DQ method has good convergancy behavior; and acceptable accuracy can be obtained by the method with a reasonable degrees of freedom.     

Three-dimensional layerwise-finite element free vibration analysis of thick laminated annular plates on elastic foundation

Using a three-dimensional layerwise-finite element method, the free vibration of thick laminated circular and annular plates supported on the elastic foundation is studied. The Pasternak-type formulation is employed to model the interaction between the plate and the elastic foundation. The discretized governing equations are derived using the Hamilton’s principle in conjunction with the layerwise theory in the thickness direction, the finite element (FE) in the radial direction and trigonometric function in the circumferential direction, respectively. The fast rate of convergence of the method is demonstrated and to verify its accuracy, comparison studies with the available solutions in the literature are performed. The effects of the geometrical parameters, the material properties and the elastic foundation parameters on the natural frequency parameters of the laminated thick circular and annular plates subjected to various boundary conditions are presented.     

A mixed method for free and forced vibration of rectangular plates

This paper presents a very first combined application of Ritz method and differential quadrature (DQ) method to vibration problem of rectangular plates. In this study, the spatial partial derivatives with respect to a coordinate direction are first discretized using the Ritz method. The resulting system of partial differential equations and the related boundary conditions are then discretized in strong form using the DQ method. The mixed method combines the simplicity of the Ritz method and high accuracy and efficiency of the DQ method. The results are obtained for various types of boundary conditions. Comparisons are made with existing analytical and numerical solutions in the literature. Numerical results prove that the present method is very suitable for the problem considered due to its simplicity, efficiency, and high accuracy.     

DQM large amplitude vibration of composite beams on nonlinear elastic foundations with restrained edges

This paper presents an efficient and accurate differential quadrature (DQ) large amplitude free vibration analysis of laminated composite thin beams on nonlinear elastic foundation. Beams under consideration have elastically restrained against rotation and in-plane immovable edges. Elastic foundation has cubic nonlinearity with shearing layer. We impose the boundary conditions directly into the governing equations in spite of the conventional DQ method and without any extra efforts. A direct iterative method is used to solve the nonlinear eigenvalue system of equations after transforming the governing equations into the frequency domain. The fast rate of convergence of the method is shown and their accuracy is demonstrated by comparing the results with those for limit cases, i.e. beams with classical boundary conditions, available in the literature. Besides, we develop a finite element program to verify the results of the presented DQ approach and to show its high computational efficiency. The effects of different parameters on the ratio of nonlinear to linear natural frequency of beams are studied.     

Modified mixed Ritz-DQ formulation for free vibration of thick rectangular and skew plates with general boundary conditions

Recently, the present authors proposed a simple mixed Ritz-differential quadrature (DQ) methodology for free and forced vibration, and buckling analysis of rectangular plates. In this technique, the Ritz method is first used to discretize the spatial partial derivatives with respect to a coordinate direction of the plate. The DQ method is then employed to analogize the resulting system of ordinary or partial differential equations. The mixed method was shown to work well for vibration and buckling problems of rectangular plates with simple boundary conditions. But, due to the use of conventional Ritz method in one coordinate direction of the plate, the geometric boundary conditions of the problem can only be satisfied in that direction. Therefore, the conventional mixed Ritz-DQ methodology may encounter difficulties when dealing with rectangular plates involving adjacent free edges and skew plates. To overcome this difficulty, this paper presents a modified mixed Ritz-DQ formulation in which all the natural boundary conditions are exactly implemented. The versatility, accuracy and efficiency of the proposed method for free vibration analysis of thick rectangular and skew plates are tested against other solution procedures. It is revealed that the proposed method can produce highly accurate solutions for the natural frequencies of thick rectangular plates involving adjacent free edges and skew plates using a small number Ritz terms and DQ sampling points.     

Flexoelectric and surface effects on the electromechanical behavior of graphene-based nanobeams

In this novel work, the electromechanical behavior of graphene-based nanocomposite (GNC) beams with flexoelectric and surface effects were investigated using size-dependent Euler-Bernoulli theory, linear piezoelectricity and Galerkin's weighted residual method along with modified strength of materials and finite element (FE) approaches. In addition, analytical and FE models were developed to study the static response of flexoelectric GNC nanobeams with various boundary conditions: cantilever, simply-supported and clamped-clamped. The developed models predict that the effective piezoelectric coefficients of GNC are responsible for the actuation capability of a graphene layer in the transverse direction due to the applied field in its axial direction and the predictions by both the models are found to be in good agreement. Results reveal that the flexoelectric and surface effects on the static response of GNC nanobeams are significant and should be taken into account. The electromechanical response of GNC nanobeams can be tailored to achieve the required coupled electromechanical characteristics of a vast range of NEMS using various boundary conditions and thickness of nanobeam as well as volume fraction of graphene. Our fundamental study sheds a light on the possibility of developing high-performance and lightweight graphene-based NEMS such as nanosensors, nanogenerators and nanoresonators using non-piezoelectric graphene.     

SBFEM elements for thin-walled composite beams with arbitrary layup

The scaled boundary finite element method (SBFEM) is a semi-analytical method in which only the boundary is discretized. The results on the boundary are scaled into the domain with respect to a scaling center which must be “visible” from the whole boundary. For beam-like problems the scaling center can be selected at infinity and only the cross-section is discretized. Two new elements for thin-walled beams have been developed on the basis of the first order shear deformation theory. The beam sections are considered to be multilayered laminate plates with arbitrary layup. The arbitrary cross-section is discretized with beam elements of Timoshenko type. Using the virtual work principle gives the SBFEM equation, which is a system of differential equations of a gyroscopic type. The solution is calculated using the matrix exponential function. The elements have been tested and compared with a finite element model and they give good results. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Accurate vibration analysis of skew plates by the new version of the differential quadrature method

An accurate free vibration analysis of skew plates is presented by using the new version of the differential quadrature method (DQM). Eight combinations of simply supported (S), clamped (C) and free (F) boundary conditions are considered. Detailed solution procedures are given and key points for success by using the DQM are emphasized. A way to simplifying the programming in using the DQM is proposed. Convergence study is made for the simply supported skew plate with a large skew angle. Good convergence of frequencies is observed. The DQ results agree very well with the existing first known accurate upper bound solutions, obtained by using Ritz method taking into considerations of the bending stress singularities occurred at corners having obtuse angles. Since slight discrepancy between the DQ data and the known accurate solutions is observed for plates with large skew angles, the DQ results are also compared with data obtained by using finite element method with very fine meshes to verify their accuracy.     

The differential quadrature element method irregular element torsion analysis model

The differential quadrature element method (DQEM) has been proposed. The element weighting coefficient matrices are generated by the differential quadrature (DQ) or generic differential quadrature (GDQ). By using the DQ or GDQ technique and the mapping procedure the governing differential or partial differential equations, the transition conditions of two adjacent elements and the boundary conditions can be discretized. A global algebraic equation system can be obtained by assembling all of the discretized equations. This method can convert a generic engineering or scientific problem having an arbitrary domain configuration into a computer algorithm. The DQEM irregular element torsion analysis model is developed.     

Primal mixed formulations for the coupling of FEM and BEM. Part I: linear problems

We apply the boundary integral equation method and a primal mixed finite element approach to study the weak solvability and Galerkin approximations of linear interior transmission problems arising in potential theory and elastostatics. The existence and uniqueness of solution of the resulting weak formulations and of the associated discrete schemes are derived by using the classical theory for variational problems with constraints. Suitable finite element subspaces of Lagrange type satisfying the compatibility conditions are utilized for defining the Galerkin scheme. The error analysis and corresponding rates of convergence are also provided.     

Analytical analysis of free vibration of non-uniform and non-homogenous beams: Asymptotic perturbation approach

This paper applies the asymptotic perturbation approach (APA) to obtain a simple analytical expression for the free vibration analysis of non-uniform and non-homogenous beams with different boundary conditions. A linear governing equation of non-uniform and non-homogeneous beams is obtained based on the Euler–Bernoulli beam theory. The perturbative theory is employed to derive an asymptotic solution of the natural frequency of the beam. Finally, numerical solutions based on the analytical method are illustrated, where the effect of a variable width ratio on the natural frequency is analyzed. To verify the accuracy of the present method, two examples, piezoelectric laminated trapezoidal beam and axially functionally graded tapered beam, are presented. The results are compared with those results obtained from the finite element method (FEM) simulation and the published literature, respectively, and a good agreement is observed for lower-order beam frequencies.     

Finite volume analysis of adaptive beams with piezoelectric sensors and actuators

This paper presents a finite volume (FV) formulation for the free vibration analysis and active vibration control of the smart beams with piezoelectric sensors and actuators. The governing equations based on Timoshenko beam theory are discretized using the finite volume method. For the purpose of forced vibration control of beam structures, the negative velocity feedback controller is designed for the single-input, single-output system. To achieve the best effect, the piezoelectric sensors and actuators are coupled with the host structure in different positions and then the performance of the designed control system is evaluated for each position. In the test examples, first the shear locking free feature of the present formulation is demonstrated. This has been performed by doing static and natural frequency analysis of some reference models. Then, the capability of the proposed method for the prediction of uncontrolled forced vibration response and active vibration control of a beam structure is studied.     

3D free vibration analysis of laminated cylindrical shell integrated piezoelectric layers using the differential quadrature method

This paper addresses the free vibration problem of multilayered shells with embedded piezoelectric layers. Based on the three-dimensional theory of elasticity, an approach combining the state space method and the differential quadrature method (DQM) is used. The shell has arbitrary end boundary conditions. For the simply supported boundary conditions closed-form solution is given by making the use of Fourier series expansion. Applying the differential quadrature method to the state space formulations along the axial direction, new state equations about state variables at discrete points are obtained for the other cases such as clamped or free end conditions. Natural frequencies of the hybrid laminated shell are presented by solving the eigenfrequency equation which can be obtained by using edges boundary condition in this state equation. Accuracy and convergence of the present approach is verified by comparing the natural frequencies with the results obtained in the literatures. Finally, the effect of edges conditions, mid-radius to thickness ratio, length to mid-radius ratio and the piezoelectric thickness on vibration behaviour of shell are investigated.     

Large-amplitude free vibrations of functionally graded beams by means of a finite element formulation

The large-amplitude free vibration analysis of functionally graded beams is investigated by means of a finite element formulation. The Von-Karman type nonlinear strain–displacement relationships are employed where the ends of the beam are constrained to move axially. The effects of the transverse shear deformation and rotary inertia are included based upon the Timoshenko beam theory. The material properties are assumed to be graded in the thickness direction according to the power-law distribution. A statically exact beam element which devoid the shear locking effect with displacement fields based on the first order shear deformation theory is used to study the geometric nonlinear effects on the vibrational characteristics of functionally graded beams. The finite element method is employed to discretize the nonlinear governing equations, which are then solved by the direct numerical integration technique in order to obtain the nonlinear vibration frequencies of functionally graded beams with different boundary conditions. The influences of power-law exponent, vibration amplitude, beam geometrical parameters and end supports on the free vibration frequencies are studied. The present numerical results compare very well with the results available from the literature where possible. Some new results for the nonlinear natural frequencies are presented in both tabular and graphical forms which can be used for future references.     

分数阶Burgers方程的有限元计算(英文)

建立了一维和二维分数阶Burgers方程的有限元格式.时间分数阶导数使用L1方法离散,空间方向使用有限元方法离散.通过选择合适的基函数,将离散后的方程转化成一个非线性代数方程组,并应用牛顿迭代方法求解.数值实验显示出了方法的有效性.     

任意边界条件下圆柱厚壳自由振动分析的空间变量变换法

本文提出一种一般解析方法——空间变量变换法,用以求解任意边界条件下圆柱厚壳自由振动问题.运用本文方法对悬臂圆柱厚壳的自振特性作了计算,计算结果与薄壳理论相应结果及试验值作了比较.理论分析和计算结果表明,本文方法具有很好的收敛性和精确性,可以推广用于分析梁、板、壳的自由振动.     

On optimal shape design of systems governed by mixed Dirichlet-Signorini boundary value problems

A problem for finding optimal shape for systems governed by the mixed unilateral boundary value problem of Dirichlet-Signorini-type is considered. Conditions for the solvability of the problem are stated when a variational inequality formulation and when a penalty method is used for solving the state problem in question. The asymptotic relation of design problems based on these two formulations is presented. The optimal shape design problem is discretized by means of finite element method. The convergence results for the approximation are proved. The discretized versions are then formulated as a non-linear programming problem. Results of practical computations of the problem in question are reported.     

Formulation and Application of the Initially Rigid Cohesive Finite Element Method

The presented procedure for cohesive crack propagation is based on an adaptive finite element (FE) implementation, which enables the introduction of cohesive surfaces in dependence on the current crack state. In contrast to already existing formulations, the focus of the present model lies on failure processes that can be described at quasi-static conditions within an implicit framework. Furthermore, an extension for mesh independent crack propagation in terms of an additional mesh adaptive formulation is presented. By the evaluation of the failure criterion considering the preferred crack direction, a new crack tip coordinate is computed and the discretization is accordingly adjusted. The remaining mesh is modified for the new boundary representation. The application of the proposed method is shown by the numerical investigation of a concrete fracture specimen from an experimental research project. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A two-dimensional layerwise-differential quadrature static analysis of thick laminated composite circular arches

In this paper, a high accuracy and rapid convergence hybrid approach is developed for two-dimensional static analyses of circular arches with different boundary conditions. The method essentially consists of a layerwise technique in the thickness direction in conjunction with differential quadrature method (DQM) in the axial direction. Hence, the high accuracy and fast convergence of DQM with generality of layerwise formulations for modeling the transverse deformations of arbitrary laminated composite thick arches are combined. This results in superior accuracy with fewer degrees of freedom than conventional finite element method (FEM) or finite difference method (FDM). The convergence behavior of the method is shown and to verify its accuracy, the results are compared with those obtained based on the first order shear deformation Reissner–Naghdi type shell theory and also higher order shear deformation theory. The effects of opening angles, ply angle, boundary conditions, and thickness-to-length ratio on the stress and displacement components are studied.