Dynamic Characteristics of Functionally Graded Timoshenko Beams by Improved Differential Quadrature Method

This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.

DIFFERENTIAL QUADRATURE METHOD TO STABILITY ANALYSIS OF PIPES CONVEYING FLUID WITH SPRING SUPPORT

It is a new attempt to extend the differential quadrature method(DQM) to stability analysis of the straight and curved centerlinepipes conveying fluid. Emphasis is placed on the study of theinfluences of several parameters on the critical flow velocity.Compared to other methods, this method can more easily deal with thepipe with spring support at its boundaries and asks for much lesscomputing effort while giving ac- ceptable precision in the numericalresults.

Differential Quadrature Method for Bending Problem of Plates with Transverse Shear Effects

A differential quadrature (DQ) method for orthotropic plates was proposed based on Reddy' s theory of plates with the effects of the higher-order transverse shear deformations. Wang-Bert's DQ approach was also further extended to handle the boundary conditions of plates. The computational convergence was studied, and the numerical results were obtained for different grid spacings and compared with the existing results. The results show that the DQ method is fairly reliable and effective.

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.

Size-dependent effect on biaxial and shear nonlinear buckling analysis of nonlocal isotropic and orthotropic micro-plate based on surface stress and modified couple stress theories using differential quadrature method

The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elasticity theories using the differential quadrature method （DQM） is presented. Main advantages of the MCST over the classical theory （CT） are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter. Based on the nonlinear von Karman assumption, the governing equations of equilibrium for the micro-classical plate consid- ering midplane displacements are derived based on the minimum principle of potential energy. Using the DQM, the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained. Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature. A parametric study is conducted to show the effects of the aspect ratio, the side-to-thickness ratio, Eringen＇s nonlocal parameter, the material length scale parameter, Young＇s modulus of the surface layer, the surface residual stress, the polymer matrix coefficients, and various boundary conditions on the dimensionless uniaxial, biaxial, and shear critical buckling loads. The results indicate that the critical buckling loads are strongly sensitive to Eringen＇s nonlocal parameter, the material length scale parameter, and the surface residual stress effects, while the effect of Young＇s modulus of the surface layer on the critical buckling load is negligible. Also, considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate. The results show that the critical biaxial buckling load increases with an increase in G12/E2 and vice versa for E1/E2. It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude. Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios, it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.

Application of Mixed Differential Quadrature Method for Solving the Coupled Two-Dimensional Incompressible Navier-Stokes Equation and Heat Equation

The traditional differential quadrature method was improved by using theupwind difference scheme for the convective terms to solve the coupled two-dimensionalincompressible Navier-stokes equations and heat equation. The new method was compared with theconventional differential quadrature method in the aspects of convergence and accuracy. The resultsshow that the new method is more accurate, and has better convergence than the conventionaldifferential quadrature method for numerically computing the steady-state solution.

DIFFERENTIAL QUADRATURE METHOD FOR BENDING OF ORTHOTROPIC PLATES WITH FINITE DEFORMATION AND TRANSVERSE SHEAR EFFECTS

Based on the Reddy's theory of plates with the effect of higher-order shear deformations, the governing equations for bending of orthotropic plates with finite deformations were established. The differential quadrature (DQ) method of nonlinear analysis to the problem was presented. New DQ approach, presented by Wang and Bert (DQWB), is extended to handle the multiple boundary conditions of plates. The techniques were also further extended to simplify nonlinear computations. The numerical convergence and comparison of solutions were studied. The results show that the DQ method presented is very reliable and valid. Moreover, the influences of geometric and material parameters as well as the transverse shear deformations on nonlinear bending were investigated. Numerical results show the influence of the shear deformation on the static bending of orthotropic moderately thick plate is significant.

Differential Quadrature Method for Steady Flow of an Incompressible Second-Order Viscoelastic Fluid and Heat Transfer Model

The two-dimensional steady flow of an incompressible second-order viscoelastic fluid between two parallel plates was studied in terms of vorticity, the stream function and temperature equations. The governing equations were expanded with respect to a snmll parameter to get the zeroth- and first-order approximate equations. By using the differenl2al quadrature method with only a few grid points, the high-accurate numerical results were obtained.

Buckling analysis of nanobeams with exponentially varying stiffness by differential quadrature method

We present the application of differential quadrature（DQ） method for the buckling analysis of nanobeams with exponentially varying stiffness based on four different beam theories of Euler-Bernoulli, Timoshenko, Reddy, and Levison.The formulation is based on the nonlocal elasticity theory of Eringen. New results are presented for the guided and simply supported guided boundary conditions. Numerical results are obtained to investigate the effects of the nonlocal parameter,length-to-height ratio, boundary condition, and nonuniform parameter on the critical buckling load parameter. It is observed that the critical buckling load decreases with increase in the nonlocal parameter while the critical buckling load parameter increases with increase in the length-to-height ratio.

Free vibration and critical speed of moderately thick rotating laminated composite conical shell using generalized differential quadrature method

The generalized differential quadrature method （GDQM） is employed to con- sider the free vibration and critical speed of moderately thick rotating laminated compos- ite conical shells with different boundary conditions developed from the first-order shear deformation theory （FSDT）. The equations of motion are obtained applying Hamilton＇s concept, which contain the influence of the centrifugal force, the Coriolis acceleration, and the preliminary hoop stress. In addition, the axial load is applied to the conical shell as a ratio of the global critical buckling load. The governing partial differential equations are given in the expressions of five components of displacement related to the points ly- ing on the reference surface of the shell. Afterward, the governing differential equations are converted into a group of algebraic equations by using the GDQM. The outcomes are achieved considering the effects of stacking sequences, thickness of the shell, rotating velocities, half-vertex cone angle, and boundary conditions. Furthermore, the outcomes indicate that the rate of the convergence of frequencies is swift, and the numerical tech- nique is superior stable. Three comparisons between the selected outcomes and those of other research are accomplished, and excellent agreement is achieved.

Static Analysis of Doubly-Curved Shell Structures of Smart Materials and Arbitrary Shape Subjected to General Loads Employing Higher Order Theories and Generalized Differential Quadrature Method

The article proposes an Equivalent Single Layer(ESL)formulation for the linear static analysis of arbitrarily-shaped shell structures subjected to general surface loads and boundary conditions.A parametrization of the physical domain is provided by employing a set of curvilinear principal coordinates.The generalized blendingmethodology accounts for a distortion of the structure so that disparate geometries can be considered.Each layer of the stacking sequence has an arbitrary orientation and is modelled as a generally anisotropic continuum.In addition,re-entrant auxetic three-dimensional honeycomb cells with soft-core behaviour are considered in the model.The unknown variables are described employing a generalized displacement field and pre-determined through-the-thickness functions assessed in a unified formulation.Then,a weak assessment of the structural problem accounts for shape functions defined with an isogeometric approach starting fromthe computational grid.Ageneralizedmethodology has been proposed to define two-dimensional distributions of static surface loads.In the same way,boundary conditions with three-dimensional features are implemented along the shell edges employing linear springs.The fundamental relations are obtained from the stationary configuration of the total potential energy,and they are numerically tackled by employing the Generalized Differential Quadrature(GDQ)method,accounting for nonuniform computational grids.In the post-processing stage,an equilibrium-based recovery procedure allows the determination of the three-dimensional dispersion of the kinematic and static quantities.Some case studies have been presented,and a successful benchmark of different structural responses has been performed with respect to various refined theories.

Investigation of the Free Vibrations of Radial Functionally Graded Circular Cylindrical Beams Based on Differential Quadrature Method

In the current research,an effective differential quadrature method(DQM)has been developed to solve natural frequency and vibration modal functions of circular section beams along radial functional gradient.Based on the high-order theory of transverse vibration of circular cross-section beams,lateral displacement equation was reconstructed neglecting circumferential shear stress.Two equations coupled with deflection and rotation angles were derived based on elastic mechanics theory and further simplified into a constant coefficient differential equation with natural frequency as eigenvalue.Then,differential quadrature method was applied to transform the eigenvalue problem of the derived differential equation into a set of algebraic equation eigenvalue problems.Natural frequencies of the free vibrations of cylindrical beams with circular cross-sections were calculated at one time,and corresponding modal functions were solved together.The obtained numerical results indicated that the natural frequencies of functionally graded(FG)circular cylindrical beams obtained using differential quadrature method agreed with the results reported in related literatures.In addition,influences of varying gradient parameters on the modal shapes of circular cylindrical beams were found to be strongly consistent with previous studies.Numerical results further validated the feasibility and accuracy of the developed differential quadrature method in solving the transverse vibration of FG circular cross-section beams.

A MIXTURE DIFFERENTIAL QUADRATURE METHOD FOR SOLVING TWO-DIMENSIONAL INCOMPRESSIBLE NAVIER_STOKES EQUATIONS

Differential quadrature method (DQM) is able to obtain highly accurate numerical solutions of differential equations just using a few grid points. But using purely differential quadrature method, good numerical solutions of two_dimensional incompressible Navier_Stokes equations can be obtained only for low Reynolds number flow and numerical solutions will not be convergent for high Reynolds number flow. For this reason, in this paper a combinative predicting_correcting numerical scheme for solving two_dimensional incompressible Navier_Stokes equations is presented by mixing upwind difference method into differential quadrature one. Using this scheme and pseudo_time_dependent algorithm, numerical solutions of high Reynolds number flow are obtained with only a few grid points. For example, 1∶1 and 1∶2 driven cavity flows are calculated and good numerical solutions are obtained.

Mechanical quadrature methods and extrapolation for solving nonlinear boundary Helmholtz integral equations

This paper presents mechanical quadrature methods （MQMs） for solving nonlinear boundary Helmholtz integral equations. The methods have high accuracy of order O（h3） and low computation complexity. Moreover, the mechanical quadrature methods are simple without computing any singular integration. A nonlinear system is constructed by discretizing the nonlinear boundary integral equations. The stability and convergence of the system are proved based on an asymptotical compact theory and the Stepleman theorem. Using the h3-Richardson extrapolation algorithms （EAs）, the accuracy to the order of O（h5） is improved. To slove the nonlinear system, the Newton iteration is discussed extensively by using the Ostrowski fixed point theorem. The efficiency of the algorithms is illustrated by numerical examples.

NUMERICAL SIMULATION OF FLOWS FOR SECOND-ORDER VISCOELASTIC FLUID COUPLED WITH HEAT TRANSFER BY DIFFERENTIAL QUADRATURE METHOD

The problem of two-dimensional steady flow of an incompressible second-order viscoelastic fluid coupled with heat transfer between parallel plates was considered. A viscous dissipation function was included in the energy equation. When the elastic property of the fluid is weaker, the zeroth-order and first-order approximate governing equations were obtained by means of the perturbation method. To understand the behavior of flow near the tube wall, the half-domain was divided into two sub-domains, in which one is a thin layer near the wall called the inner domain and the remainder is called the outer domain. The governing equations in the inner domain and in the outer domain were discretized respectively by using the Differential Quadrature Method (DQM). The matching conditions at the interface between the inner and outer domains were presented. An iterative method for solving these discretized equations was given in this paper. The numerical results obtained agree with existing results.

Numerical Simulation of MHD Peristaltic Flow with Variable Electrical Conductivity and Joule Dissipation Using Generalized Differential Quadrature Method

In this paper, the MHD peristaltic flow inside wavy walls of an asymmetric channel is investigated, where the walls of the channel are moving with peristaltic wave velocity along the channel length. During this investigation,the electrical conductivity both in Lorentz force and Joule heating is taken to be temperature dependent. Also, the long wavelength and low Reynolds number assumptions are utilized to reduce the governing partial differential equations into a set of coupled nonlinear ordinary differential equations. The new set of obtained equations is then numerically solved using the generalized differential quadrature method(GDQM). This is the first attempt to solve the nonlinear equations arising in the peristaltic flows using this method in combination with the Newton-Raphson technique. Moreover, in order to check the accuracy of the proposed numerical method, our results are compared with the results of built-in Mathematica command NDSolve. Taking Joule heating and viscous dissipation into account, the effects of various parameters appearing in the problem are used to discuss the fluid flow characteristics and heat transfer in the electrically conducting fluids graphically. In presence of variable electrical conductivity, velocity and temperature profiles are highly decreasing in nature when the intensity of the electrical conductivity parameter is strengthened.

SPLITTING EXTRAPOLATIONS FOR SOLVING BOUNDARY INTEGRAL EQUATIONS OF LINEAR ELASTICITY DIRICHLET PROBLEMS ON POLYGONS BY MECHANICAL QUADRATURE METHODS

Taking hm as the mesh width of a curved edge Гm （m = 1, ..., d ） of polygons and using quadrature rules for weakly singular integrals, this paper presents mechanical quadrature methods for solving BIES of the first kind of plane elasticity Dirichlet problems on curved polygons, which possess high accuracy O（h0^3） and low computing complexities. Since multivariate asymptotic expansions of approximate errors with power hi^3 （i = 1, 2, ..., d） are shown, by means of the splitting extrapolations high precision approximations and a posteriori estimate are obtained.

Free Vibration Analysis of Laminated Composite Beams Using Differential Quadrature Method

A higher-order theory for laminated composite beams is used to study the free vibration of laminated composite beams, and the differential quadrature method is employed to obtain the numerical solution of the governing differential equations. Free vibration analysis of beams with rectangular cross-section for various combinations of end conditions is studied. The results show that the differential quadrature method is reliable and accurate compared with other available results.

Stability Analysis and Order Improvement for Time Domain Differential Quadrature Method

The differential quadrature method has been widely used in scientific and engineering computation.However,for the basic characteristics of time domain differential quadrature method,such as numerical stability and calculation accuracy or order,it is still lack of systematic analysis conclusions.In this paper,according to the principle of differential quadrature method,it has been derived and proved that the weighting coefficients matrix of differential quadrature method meets the important V-transformation feature.Through the equivalence of the differential quadrature method and the implicit Runge-Kutta method,it has been proved that the differential quadrature method is A-stable and s-stage s-order method.On this basis,in order to further improve the accuracy of the time domain differential quadrature method,a class of improved differential quadrature method of s-stage 2s-order has been proposed by using undetermined coefficients method and Pad´e approximations.The numerical results show that the proposed differential quadrature method is more precise than the traditional differential quadrature method.

Surface Effect on Vibration of Timoshenko Nanobeam Based on Generalized Differential Quadrature Method and Molecular Dynamics Simulation

Nanobeams have promising applications in areas such as sensors,actuators,and resonators in nanoelectromechanical systems(NEMS).Considering the effects of gyration inertia,surface layer mass,surface residual stress,and surface Young's modulus,this study develops the vibration equations of the Timoshenko nanobeam.The generalized differential quadrature(GDQ)method and molecular dynamics(MD)simulation are used to study the surface effect on vibration.For a rectangular cross section,surface residual stress and surface Young's modulus are all affected by the height of the cross section rather than by the length-height ratio.If surface layer mass is considered,then the first three natural frequencies all decrease relative to their counterparts in the case in which surface layer mass is ignored.Results show that the effect of gyration inertia on resonance frequency is negligible.Longitudinal vibration does not easily occur relative to the bending and rotation vibrations of nanobeams.In addition,the results obtained by the GDQ method fit those obtained by MD simulation for beams with length-height ratios of 4-8.This study provides insights into the mechanism of the vibration of short and deep nanobeams and sheds light on the quantitative design of the elements in NEMSs.