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.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

Comparative Studies between Picard’s and Taylor’s Methods of Numerical Solutions of First Ordinary Order Differential Equations Arising from Real-Life Problems

To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.

Adaptive multi-step piecewise interpolation reproducing kernel method for solving the nonlinear time-fractional partial differential equation arising from financial economics

This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.

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.

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.

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.

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.

Quasi-static and dynamical analyses of a thermoviscoelastic Timoshenko beam using the differential quadrature method

The quasi-static and dynamic responses of a thermoviscoelastic Timoshenko beam subject to thermal loads are analyzed. First, based on the small geometric deformation assumption and Boltzmann constitutive relation, the governing equations for the beam are presented. Second, an extended differential quadrature method(DQM)in the spatial domain and a differential method in the temporal domain are combined to transform the integro-partial-differential governing equations into the ordinary differential equations. Third, the accuracy of the present discrete method is verified by elastic/viscoelastic examples, and the effects of thermal load parameters, material and geometrical parameters on the quasi-static and dynamic responses of the beam are discussed. Numerical results show that the thermal function parameter has a great effect on quasi-static and dynamic responses of the beam. Compared with the thermal relaxation time, the initial vibrational responses of the beam are more sensitive to the mechanical relaxation time of the thermoviscoelastic material.

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.

ERROR ANALYSIS IN DIFFERENTIAL QUADRATURE METHOD

Error estimation for the differential quadrature (DQ)with various sets of grid spacing is presented. A general formula is given to compute the weighting coefficients directly. It is found that the maximum error of the Do method with roots of Chebyshev polynomials including two end POints (-1 and + 1 ) is the smallest among the several sets of grid spacing investigated herein.

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.

Buckling analysis of functionally graded plates partially resting on elastic foundation using the differential quadrature element method

We extend the differential quadrature element method (DQEM) to the buckling analysis of uniformly in-plane loaded functionally graded (FG) plates fully or partially resting on the Pasternak model of elastic support. Material properties of the FG plate are graded in the thickness direction and assumed to obey a power law distribution of the volume fraction of the constituents. To set up the global eigenvalue equation, the plate is divided into sub-domains or elements and the generalized differential quadrature procedure is applied to discretize the governing, boundary and compatibility equations. By assembling discrete equations at all nodal points, the weighting coefficient and force matrices are derived. To validate the accuracy of this method, the results are compared with those of the exact solution and the finite element method. At the end, the effects of different variables and local elastic support arrangements on the buckling load factor are investigated.

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.

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.

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.

Nonlinear Analysis of Organic Polymer Solar Cells Using Differential Quadrature Technique with Distinct and Unique Shape Function

Four numerical schemes are introduced for the analysis of photocurrent transients in organic photovoltaic devices.Themathematicalmodel for organic polymer solar cells contains a nonlinear diffusion-reaction partial differential equation system with electrostatic convection attached to a kinetic ordinary differential equation.To solve the problem,Polynomial-based differential quadrature,Sinc,and Discrete singular convolution are combined with block marching techniques.These schemes are employed to reduce the problem to a nonlinear algebraic system.The iterative quadrature technique is used to solve the reduced problem.The obtained results agreed with the previous exact one and the finite element method.Further,the effects of different times,different mobilities,different densities,different geminate pair distances,different geminate recombination rate constants,different generation efficiencies,and supporting conditions on photocurrent have been analyzed.The novelty of this paper is that these schemes for photocurrent transients in organic polymer solar cells have never been presented before,so the results may be useful for improving the performance of solar cells.

Structural dynamic responses analysis applying differential quadrature method

Unconditionally stable higher-order accurate time step integration algorithms based on the differential quadrature method (DQM) for second-order initial value problems were applied and the quadrature rules of DQM, computing of the weighting coefficients and choices of sampling grid points were discussed. Some numerical examples dealing with the heat transfer problem, the second-order differential equation of imposed vibration of linear single-degree-of-freedom systems and double-degree-of-freedom systems, the nonlinear move differential equation and a beam forced by a changing load were computed, respectively. The results indicated that the algorithm can produce highly accurate solutions with minimal time consumption, and that the system total energy can remain conservative in the numerical computation.

An inverse analysis of fluid flow through granular media using differentiable lattice Boltzmann method

This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.