A mixed modal-differential quadrature method for free and forced vibration of beams in contact with fluid

A simple and accurate mixed modal-differential quadrature formulation is proposed to study the dynamic behavior of beams in contact with fluid. Both free and forced vibration problems are considered. The proposed mixed methodology uses the modal technique for the structural domain while it applies the differential quadrature method (DQM) to the fluid domain. Thus, the governing partial differential equations of the beam and fluid are reduced to a set of ordinary differential equations in time. In the case of forced vibration, the Newmark time integration scheme is employed to solve the resulting system of ordinary differential equations. The proposed formulation, in general, combines the simplicity of the modal method and high accuracy and efficiency of the DQM. Its application is shown by solving some beam-fluid interaction problems. Comparisons with analytical solutions show that the present method is very accurate and reliable. To demonstrate its efficiency, the test problems are also solved using the finite element method (FEM). It is found that the proposed method can produce better accuracy than the FEM using less computational time. The technique presented in this investigation is general and can be used to solve various fluid-structure interaction problems.     

Mixed finite element and differential quadrature method for free and forced vibration and buckling analysis of rectangular plates

This paper presents a combined application of the finite element method (FEM) and the differential quadrature method (DQM) to vibration and buckling problems of rectangular plates. The proposed scheme combines the geometry flexibility of the FEM and the high accuracy and efficiency of the DQM. The accuracy of the present method is demonstrated by comparing the obtained results with those available in the literature. It is shown that highly accurate results can be obtained by using a small number of finite elements and DQM sample points. The proposed method is suitable for the problems considered due to its simplicity and potential for further development.     

基于高阶剪切变形理论的四边形求积元板单元及其应用

近年来由各类新型复合材料或功能梯度材料构成的板结构在工程领域得到了广泛应用,其显著特点是材料性能沿板厚变化.为合理考虑横向剪切应变,许多学者基于Reddy高阶剪切变形理论,构建了不同的有限元单元对该类板结构进行分析,但其中满足$C^{1}$连续条件的单元相对较少.本文基于Reddy高阶剪切变形理论,采用求积元方法,建立了$C^{1}$连续的四边形板单元.利用该单元对均质材料、复合材料、功能梯度材料构成的等厚度矩形板、变厚度矩形板及等厚度斜板的线弹性弯曲和自由振动问题进行了计算分析,并与现有文献中的相应计算结果进行了对比.研究表明:基于高阶剪切变形理论的四边形求积元板单元具有较高的计算效率和良好的适应性,文中各类材料构成的等变厚度矩形板及等厚度斜板均只需1个单元即可得到理想的计算结果.对于等/变厚度矩形板,可仅使用9$\times$9个积分点,而对于等厚度斜板,随着斜角的增大,所需积分点的数目逐渐增多至15$\times $15.该四边形求积元板单元可进一步用于新型复合材料板的非线性分析.      

DSC-Ritz element method for vibration analysis of rectangular Mindlin plates with mixed edge supports

Based on Mindlin's first-order shear deformable plate theory, a DSC-Ritz element method is developed for the free vibration analysis of moderately thick rectangular plates with mixed supporting edges. The rationale of the present approach is not only to apply the discrete singular convolution (DSC) delta type wavelet kernel as a trial function with the Ritz method, but also to incorporate the method in finite elements in order to handle the mixed boundary constraints. The approach is novel and flexible as it passes through a bottleneck of the global DSC-Ritz method in treating the kinematic supporting edges with assorted discontinuities. A series of numerical simulations for rectangular Mindlin plates with various edge support discontinuities, plate thicknesses and aspect ratios are presented. For verification, the vibration frequencies thus established are directly compared with those reported in the open literature. New sets of numerical results for several other cases of moderately thick plates with mixed simply supported, clamped and free edges are presented and discussed in detail.     

Free vibration of composite plates with mixed boundary conditions based on higher-order shear deformation theory

A global higher-order shear deformation theory is devised to obtain the governing equations of composite plates under dynamic excitation. The time-harmonic solution leads to an eigenvalue problem for the natural frequencies of plates. The eigenvalue problem for rectangular plates is converted to a set of homogenous algebraic equations using differential quadrature method. The formulation of the problem allows direct application of various boundary conditions. Therefore, rectangular plates with mixed boundary conditions are also considered. To show the validity of results, the fundamental natural frequencies of composite plates with different boundary conditions and those of isotropic plates with mixed boundary conditions are compared against the results available in the literature.     

Buckling analysis of functionally graded arbitrary straight-sided quadrilateral plates on elastic foundations

As a first endeavor, the buckling analysis of functionally graded (FG) arbitrary straight-sided quadrilateral plates rested on two-parameter elastic foundation under in-plane loads is presented. The formulation is based on the first order shear deformation theory (FSDT). The material properties are assumed to be graded in the thickness direction. The solution procedure is composed of transforming the governing equations from physical domain to computational domain and then discretization of the spatial derivatives by employing the differential quadrature method (DQM) as an efficient and accurate numerical tool. After studying the convergence of the method, its accuracy is demonstrated by comparing the obtained solutions with the existing results in literature for isotropic skew and FG rectangular plates. Then, the effects of thickness-to-length ratio, elastic foundation parameters, volume fraction index, geometrical shape and the boundary conditions on the critical buckling load parameter of the FG plates are studied.     

Differential quadrature element method for static analysis of Reissner–Mindlin polar plates

In this paper, a new numerical technique, the differential quadrature element method (DQEM) , has been developed for static analysis of the two-dimensional polar Reissner–Mindlin plate in the polar coordinate system by integrating the domain decomposition method (DDM) with the differential quadrature method (DQM) . The detailed formulations for the sectorial DQEM plate bending element and the compatibility conditions between each element are presented. The convergence properties and the accuracy of the DQEM for bending of thick polar plates are investigated through a number of numerical computations. Consequently, the DQEM has been successfully applied to analyze several annular sector plates with discontinuous loading and boundary conditions and cutouts to illustrate the simplicity and flexibility of this method for solving Reissner–Mindlin plates in polar coordinate system which are not solvable directly using the differential quadrature method. The numerical results are verified by the existing exact solutions or the FEM solutions obtained using the software package ANSYS (Version 5.3) .     

Nonlinear vibration of symmetrically laminated composite skew plates by finite element method

Here, the large amplitude free flexural vibration behavior of symmetrically laminated composite skew plates is investigated using the finite element method. The formulation includes the effects of shear deformation, in-plane and rotary inertia. The geometric non-linearity based on von Kármán's assumptions is introduced. The nonlinear matrix amplitude equation obtained by employing Galerkin's method is solved by direct iteration technique. Time history for the nonlinear free vibration of composite skew plate is also obtained using Newmark's time integration technique to examine the accuracy of matrix amplitude equation. The variation of nonlinear frequency ratios with amplitudes is brought out considering different parameters such as skew angle, fiber orientation and boundary condition.     

Three-dimensional free vibration analysis of multi-directional functionally graded piezoelectric annular plates on elastic foundations via state space based differential quadrature method

The three-dimensional free vibration analysis of a multi-directional functionally graded piezoelectric(FGP) annular plate resting on two parameter(Pasternak)elastic foundations is investigated under different boundary conditions. The material properties are assumed to vary continuously along the radial and thickness directions and have exponent-law distribution. A semi-analytical approach named the state space based differential quadrature method(SSDQM) is used to provide an analytical solution along the thickness using the state space method(SSM) and an approximate solution along the radial direction using the one-dimensional differential quadrature method(DQM).The influence of the Winkler and shear stiffness of the foundation, the material property graded variations, and the circumferential wave number on the non-dimensional natural frequency of multi-directional FGP annular plates is studied.     

Large amplitude free flexural vibrations of laminated composite skew plates

Here, the large amplitude free flexural vibration behaviors of thin laminated composite skew plates are investigated using finite element approach. The formulation includes the effects of shear deformation, in-plane and rotary inertia. The geometric non-linearity based on von Karman's assumptions is introduced. The non-linear governing equations obtained employing Lagrange's equations of motion are solved using the direct iteration technique. The variation of non-linear frequency ratios with amplitudes is brought out considering different parameters such as skew angle, number of layers, fiber orientation, boundary condition and aspect ratio. The influence of higher vibration modes on the non-linear dynamic behavior of laminated skew plates is also highlighted. The present study reveals the redistribution of vibrating mode shape at certain amplitude of vibration depending on geometric and lamination parameters of the plate. Also, the degree of hardening behavior increases with the skew angle and its rate of change depends on the level of amplitude of vibration.     

Static and free vibration analysis of graphene platelets reinforced composite truncated conical shell,cylindrical shell,and annular plate using theory of elasticity and DQM

Abstract

In this paper, three-dimensional static and free vibration analysis of functionally graded graphene platelets-reinforced composite (FG-GPLRC) truncated conical shells, cylindrical shells and annular plates with various boundary conditions is carried out within the framework of elasticity theory. The main contribution of the present work is that formulation for free vibration and bending behavior of the FG-GPLRC truncated conical shell based on theory of elasticity has not yet been reported. Additionally, formulation and solution for cylindrical shell and annular plate are derived by changing the semi vertex angle in formulation and solution of FG-GPLRC truncated conical shell. A semi-analytical solution is proposed base on employing differential quadrature method (DQM) together with state-space technique. Validity of current approach is assessed by comparing its numerical results with those available in the literature. An especial attention is drawn to the role of GPLs weight fraction, patterns of GPLs distribution through the thickness direction, geometrical parameters such as semi-vertex angle, length to mid-radius ratio on natural frequencies and bending characteristics. Numerical results reveal that desirable static and free vibration response (such as lower radial deflection and higher natural frequencies) can be achieved by locating more square shaped GPLs near inner and outer surfaces.     

弹性约束边界圆环板面内自由振动的二维弹性解

基于二维线弹性理论,应用哈密顿原理导出弹性约束边界圆环板面内自由振动的控制微分方程。采用微分求积法(DQM)数值研究了弹性约束边界圆环板面内自由振动的频率特性。通过设置弹性刚度系数为0或∞,问题退化为四种典型边界圆环板的面内自由振动,与已有文献的计算数值结果进行比较,证实本文的分析求解方法行之有效。最后全面考虑了圆环板边界条件、几何系数及刚度系数对自振频率的影响。     

Vibration analysis of cracked plates using the extended finite element method

In the present paper, the extended finite element method (X-FEM) is adopted to analyze vibrations of cracked plates. Mindlin’s plate theory taking into account the effects of shear deformation and rotatory inertia is included in the development of the model. First, conventional FEM without any discontinuity is carried out, then the enrichment proposed by Moës et al. (Int J Numer Methods Eng 46, 131–150, 1999) of nodal elements containing cracks is added to the FEM formulation. Numerical implementation of enriched elements by discontinuous functions is performed, and thus dynamic equations (stiffness and mass matrices) are established. A FORTRAN computer code based on the X-FEM formulation is hence developed. Rectangular and square plates containing through-edge and central cracks with different boundary conditions are considered. The subspace iteration method is used to solve the eigenvalue problem. Natural frequencies as well as the corresponding eigenfunctions are consequently calculated as a function of the crack length. The obtained results show that the X-FEM is an efficient method in the dynamic analysis of plates containing discontinuities.     

NON-AXISYMMETRIC MIXED BOUNDARY VALUE PROBLEM FOR DYNAMIC INTERACTION BETWEEN SATURATED POROUS FOUNDATION AND RECTANGULAR PLATE$^{\bf 1)}$

在同一界面的不同区域具有多种边界条件, 称之为混合边界, 这是一个熟知的力学问题. 对这类问题进行精确分析时, 必须要进行混合边值问题的求解. 而对于一般的三维非轴对称情形, 混合边值问题的求解往往存在数学困难. 本文利用Hilbert定理和双重Fourier变换, 给出了一种求解三维非轴对称混合边值问题的解析方法, 利用该方法对具有混合透水边界的饱和多孔地基上矩形板的振动弯曲进行了解析研究(板与地基接触面为不透水边界, 其余为透水边界). 首先, 基于Kirchhoff理论和Biot多孔介质理论建立矩形板与饱和多孔地基的动力控制方程, 进行耦合求解. 针对板土接触面和非接触面的混合边值问题, 采用双重Fourier变换构造出两对二维对偶积分方程, 以接触应力和接触面孔隙压力为基本未知量, 用Jacobi正交多项式将未知量展开, 再利用Schmidt法对二维对偶积分方程完成求解, 最终推导出板土系统在动力作用下的位移和应力解析式. 通过将本文计算模型退化为单一弹性地基, 与已有研究结果进行对比, 验证了本文方法的正确性和有效性. 最后, 通过数值算例, 对饱和多孔地基上矩形板的动力响应及参数影响做出分析和讨论. 此外, 本文提出的解析法具有一般性, 可广泛应用于复杂接触问题和多场耦合问题的求解.     

饱和多孔地基与矩形板动力相互作用的非轴对称混合边值问题

在同一界面的不同区域具有多种边界条件, 称之为混合边界, 这是一个熟知的力学问题. 对这类问题进行精确分析时, 必须要进行混合边值问题的求解. 而对于一般的三维非轴对称情形, 混合边值问题的求解往往存在数学困难. 本文利用Hilbert定理和双重Fourier变换, 给出了一种求解三维非轴对称混合边值问题的解析方法, 利用该方法对具有混合透水边界的饱和多孔地基上矩形板的振动弯曲进行了解析研究(板与地基接触面为不透水边界, 其余为透水边界). 首先, 基于Kirchhoff理论和Biot多孔介质理论建立矩形板与饱和多孔地基的动力控制方程, 进行耦合求解. 针对板土接触面和非接触面的混合边值问题, 采用双重Fourier变换构造出两对二维对偶积分方程, 以接触应力和接触面孔隙压力为基本未知量, 用Jacobi正交多项式将未知量展开, 再利用Schmidt法对二维对偶积分方程完成求解, 最终推导出板土系统在动力作用下的位移和应力解析式. 通过将本文计算模型退化为单一弹性地基, 与已有研究结果进行对比, 验证了本文方法的正确性和有效性. 最后, 通过数值算例, 对饱和多孔地基上矩形板的动力响应及参数影响做出分析和讨论. 此外, 本文提出的解析法具有一般性, 可广泛应用于复杂接触问题和多场耦合问题的求解.      

Mechanical buckling and free vibration of thick functionally graded plates resting on elastic foundation using the higher order B-spline finite strip method

In this study, the mechanical buckling and free vibration of thick rectangular plates made of functionally graded materials (FGMs) resting on elastic foundation subjected to in-plane loading is considered. The third order shear deformation theory (TSDT) is employed to derive the governing equations. It is assumed that the material properties of FGM plates vary smoothly by distribution of power law across the plate thickness. The elastic foundation is modeled by the Winkler and two-parameter Pasternak type of elastic foundation. Based on the spline finite strip method, the fundamental equations for functionally graded plates are obtained by discretizing the plate into some finite strips. The results are achieved by the minimization of the total potential energy and solving the corresponding eigenvalue problem. The governing equations are solved for FGM plates buckling analysis and free vibration, separately. In addition, numerical results for FGM plates with different boundary conditions have been verified by comparing to the analytical solutions in the literature. Furthermore, the effects of different values of the foundation stiffness parameters on the response of the FGM plates are determined and discussed.     

Postbuckling of long unsymmetrically laminated composite plates under axial compression

An approximate solution is given for the postbuckling of infinitely long and unsymmetrically laminated composite plates. This solution is obtained by superposing a polynomial transverse displacement given by bending due to unsymmetric laminate configurations and a simple functional representation for the buckling mode in conjunction with the Galerkin method. Nondimensional parameters are used to express the approximate solution in a very simple and clear formulation. The results given by this solution for axial compression in the longitudinal direction are compared with the results given by the nonlinear finite element method (FEM) for finite length rectangular long plates. The influence of the boundary conditions on postbuckling response is also studied. For the FEM analysis, two different simply supported boundary conditions on the long edges of the plate are considered. It is found that these two sets of boundary conditions give different results for the buckling and postbuckling finite element analysis. In most cases the FEM analysis overestimate and, respectively, underestimate the approximate closed form solution, depending on the type of simply supported boundary condition considered. Thus, the approximate solution appears useful for design purposes as an averaged quantity between the two FEM analyses. Also, it is found that the reduced bending stiffness method can be successfully used for determining the approximate solution.     

基于边光滑有限元法的加筋板静力和自由振动分析

运用边光滑有限元法,研究分析了加筋板结构的静力和自由振动问题。在边光滑有限元法中,将基于边的应变光滑技术用于对原来的应变场进行光滑操作;由于应变光滑技术能够适当地软化原来过刚的有限元模型,从而能够得到更加接近于系统准确刚度的光滑有限元模型;鉴于三角形单元良好的适用性,选用三角形单元对模型进行网格划分;同时,为了解决低阶Reissner-Mindlin板单元弯曲过程中的横向剪切自锁问题,采用了一种新型的离散剪切间隙技术。算例的数值计算结果表明,与传统的有限元法相比,边光滑有限元法能够得到精度更高的计算结果,且收敛更快,计算效率更佳。     

Higher-order shear deformable theories for flexure of sandwich plates—Finite element evaluations

A simple isoparametric finite element formulation based on a higher-order displacement model for flexure analysis of multilayer symmetric sandwich plates is presented. The assumed displacement model accounts for non-linear variation of inplane displacements and constant variation of transverse displacement through the plate thickness. Further, the present formulation does not require the fictitious shear correction coefficient(s) generally associated with the first-order shear deformable theories. Two sandwich plate theories are developed: one in which the free shear stress conditions on the top and bottom bounding planes are imposed and another, in which such conditions are not imposed. The validity of the present development(s) is established through, numerical evaluations for deflections/stresses/stress-resultants and their comparisons with the available three-dimensional analyses/closed-form/other finite element solutions. Comparison of results from thin plate. Mindlin and present analyses with the exact three-dimensional analyses yields some important conclusions regarding the effects of the assumptions made in the CPT and Mindlin type theories. The comparative study further establishes the necessity of a higher-order shear deformable theory incorporating warping of the cross-section particularly for sandwich plates.     

Nonlinear vibration analysis of laminated composite Mindlin micro/nano-plates resting on orthotropic Pasternak medium using DQM

The nonlocal nonlinear vibration analysis of embedded laminated microplates resting on an elastic matrix as an orthotropic Pasternak medium is investigated. The small size effects of micro/nano-plate are considered based on the Eringen nonlocal theory. Based on the orthotropic Mindlin plate theory along with the von Kármán geometric nonlinearity and Hamilton's principle, the governing equations are derived. The differential quadrature method (DQM) is applied for obtaining the nonlinear frequency of system. The effects of different parameters such as nonlocal parameters, elastic media, aspect ratios, and boundary conditions are considered on the nonlinear vibration of the micro-plate. Results show that considering elastic medium increases the nonlinear frequency of system. Furthermore, the effect of boundary conditions becomes lower at higher nonlocal parameters.