An efficient mixed methodology for free vibration and buckling analysis of orthotropic rectangular plates

In this paper, a simple and efficient mixed Ritz-differential quadrature (DQ) method is presented for free vibration and buckling analysis of orthotropic rectangular plates. The mixed scheme combines the simplicity of the Ritz method and high accuracy and efficiency of the DQ method. The accuracy of the proposed method is demonstrated by comparing the calculated results with those available in the literature. It is shown that highly accurate results can be obtained using a small number of Ritz terms and DQ sample points. The proposed method is suitable for the problem considered due to its simplicity and potential for further development.     

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.     

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.     

矩形平板自由和强迫振动以及屈曲的有限元-微分求积混合法(FE-DQM)

首次提出有限元法(FEM)和微分求积法(DQM)的组合应用,分析矩形平板的振动和屈曲问题.混合法综合了FEM几何适应性强,以及DQM的高精度和高效率.与已有文献的计算结果比较,验证了该方法的正确性.研究表明,使用少量的有限单元和不多的DQM样本点,就可以得到高精度的结果.由于该方法简单且具备进一步发展的潜力,被认为适用于这类问题的求解.     

Bernstein Ritz‐Galerkin method for solving an initial‐boundary value problem that combines Neumann and integral condition for the wave equation

In this article, the Ritz‐Galerkin method in Bernstein polynomial basis is implemented to give an approximate solution of a hyperbolic partial differential equation with an integral condition. We will deal here with a type of nonlocal boundary value problem, that is, the solution of a hyperbolic partial differential equation with a nonlocal boundary specification. The nonlocal conditions arise mainly when the data on the boundary cannot be measured directly. The properties of Bernstein polynomial and Ritz‐Galerkin method are first presented, then Ritz‐Galerkin method is used to reduce the given hyperbolic partial differential equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique presented in this article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Coupling dual reciprocity BEM and differential quadrature method for time-dependent diffusion problems

This paper presents the very first combined application of dual reciprocity BEM (DRBEM) and differential quadrature (DQ) method to time-dependent diffusion problems. In this study, the DRBEM is employed to discretize the spatial partial derivatives. The DQ method is then applied to analogize temporal derivatives. The resulting algebraic formulation is the known Lyapunov matrix equation, which can be very efficiently solved by the Bartels–Stewart algorithms. The mixed scheme combines strong geometry flexibility and boundary-only feature of the BEM and high accuracy and efficiency of the DQ method. Its superiority is demonstrated through the solution of some benchmark diffusion problems. The DQ method is shown to be numerically accurate, stable and computationally efficient in computing dynamic problems. In particular, the present study reveals that the DRBEM is also very efficient for transient diffusion problems with Dirichlet boundary conditions by coupling the DQ method in time discretization.     

A mixed differential quadrature and finite element free vibration and buckling analysis of thick beams on two-parameter elastic foundations

As a first endeavor, a mixed differential quadrature (DQ) and finite element (FE) method for boundary value structural problems in the context of free vibration and buckling analysis of thick beams supported on two-parameter elastic foundations is presented. The formulations are based on the two-dimensional theory of elasticity. The problem domain along axial direction is discretized using finite elements. The resulting system of equations and the related boundary conditions are discretized in the thickness direction and in strong-form using DQM. The method benefits from low computational efforts of the DQ in conjunction with the effectiveness of the FE method in general geometry and systematic boundary treatment resulting in highly accurate and fast convergence behavior solution. The boundary conditions at the top and bottom surface of the beams are implemented accurately. The presented formulations provide an effective analysis tool for beams free of shear locking. Comparisons are made with results from elasticity solutions as well as higher-order beam theory.     

On the solution of the non-local parabolic partial differential equations via radial basis functions

In this paper, the problem of solving the one-dimensional parabolic partial differential equation subject to given initial and non-local boundary conditions is considered. The approximate solution is found using the radial basis functions collocation method. There are some difficulties in computing the solution of the time dependent partial differential equations using radial basis functions. If time and space are discretized using radial basis functions, the resulted coefficient matrix will be very ill-conditioned and so the corresponding linear system cannot be solved easily. As an alternative method for solution, we can use finite-difference methods for discretization of time and radial basis functions for discretization of space. Although this method is easy to use but an accurate solution cannot be provided. In this work an efficient collocation method is proposed for solving non-local parabolic partial differential equations using radial basis functions. Numerical results are presented and are compared with some existing methods.     

A Comparison Study of Deep Galerkin Method and Deep Ritz Method for Elliptic Problems with Different Boundary Conditions

Recent years have witnessed growing interests in solving partial differential equations by deep neural networks, especially in the high-dimensional case. Unlike classical numerical methods, such as finite difference method and finite element method, the enforcement of boundary conditions in deep neural networks is highly nontrivial. One general strategy is to use the penalty method. In the work, we conduct a comparison study for elliptic problems with four different boundary conditions, i.e., Dirichlet, Neumann, Robin, and periodic boundary conditions, using two representative methods: deep Galerkin method and deep Ritz method. In the former, the PDE residual is minimized in the least-squares sense while the corresponding variational problem is minimized in the latter. Therefore, it is reasonably expected that deep Galerkin method works better for smooth solutions while deep Ritz method works better for low-regularity solutions. However, by a number of examples, we observe that deep Ritz method can outperform deep Galerkin method with a clear dependence of dimensionality even for smooth solutions and deep Galerkin method can also outperform deep Ritz method for low-regularity solutions.Besides, in some cases, when the boundary condition can be implemented in an exact manner, we find that such a strategy not only provides a better approximate solution but also facilitates the training process.     

Galerkin approximations with embedded boundary conditions for retarded delay differential equations

Finite-dimensional approximations are developed for retarded delay differential equations (DDEs). The DDE system is equivalently posed as an initial-boundary value problem consisting of hyperbolic partial differential equations (PDEs). By exploiting the equivalence of partial derivatives in space and time, we develop a new PDE representation for the DDEs that is devoid of boundary conditions. The resulting boundary condition-free PDEs are discretized using the Galerkin method with Legendre polynomials as the basis functions, whereupon we obtain a system of ordinary differential equations (ODEs) that is a finite-dimensional approximation of the original DDE system. We present several numerical examples comparing the solution obtained using the approximate ODEs to the direct numerical simulation of the original non-linear DDEs. Stability charts developed using our method are compared to existing results for linear DDEs. The presented results clearly demonstrate that the equivalent boundary condition-free PDE formulation accurately captures the dynamic behaviour of the original DDE system and facilitates the application of control theory developed for systems governed by ODEs.     

Ritz–Galerkin method for solving a parabolic equation with non‐local and time‐dependent boundary conditions

The paper is devoted to the investigation of a parabolic partial differential equation with non‐local and time‐dependent boundary conditions arising from ductal carcinoma in situ model. Approximation solution of the present problem is implemented by the Ritz–Galerkin method, which is a first attempt at tackling parabolic equation with such non‐classical boundary conditions. In the process of dealing with the difficulty caused by integral term in non‐local boundary condition, we use a trick of introducing the transition function G(x,t) to convert non‐local boundary to another non‐classical boundary, which can be handled with the Ritz–Galerkin method. Illustrative examples are included to demonstrate the validity and applicability of the technique in this paper. Copyright © 2015 John Wiley & Sons, Ltd.     

A Cartesian‐grid collocation method based on radial‐basis‐function networks for solving PDEs in irregular domains

This paper reports a new Cartesian‐grid collocation method based on radial‐basis‐function networks (RBFNs) for numerically solving elliptic partial differential equations in irregular domains. The domain of interest is embedded in a Cartesian grid, and the governing equation is discretized by using a collocation approach. The new features here are (a) one‐dimensional integrated RBFNs are employed to represent the variable along each line of the grid, resulting in a significant improvement of computational efficiency, (b) the present method does not require complicated interpolation techniques for the treatment of Dirichlet boundary conditions in order to achieve a high level of accuracy, and (c) normal derivative boundary conditions are imposed by means of integration constants. The method is verified through the solution of second‐ and fourth‐order PDEs; accurate results and fast convergence rates are obtained. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Suboptimal control of a 2D molten carbonate fuel cell PDAE model

Numerical results for suboptimal boundary control of an integro partial differential algebraic equation system of dimension 28 are presented. The application is a molten carbonate fuel cell power plant. The technically and economically important fast tracking of the new stationary cell voltage during a load change is optimized. The nonstandard optimal control problem s.t. degenerated PDE is discretized by the method of lines yielding a very large DAE constrained standard optimal control problem. An index analysis is performed to identify consistent initial conditions.     

A hybrid DQ/LMQRBF-DQ approach for numerical solution of Poisson-type and Burger’s equations in irregular domain

Differential quadrature (DQ) is an efficient and accurate numerical method for solving partial differential equations (PDEs). However, it can only be used in regular domains in its conventional form. Local multiquadric radial basis function-based differential quadrature (LMQRBF-DQ) is a mesh free method being applicable to irregular geometry and allowing simple imposition of any complex boundary condition. Implementation of the latter numerical scheme imposes high computational cost due to the necessity of numerous matrix inversions. It also suffers from sensitivity to shape parameter(s). This paper presents a new method through coupling the conventional DQ and LMQRBF-DQ to solve PDEs. For this purpose, the computational domain is divided into a few rectangular shapes and some irregular shapes. In such a domain decomposition process, a high percentage of the computational domain will be covered by regular shapes thus taking advantage of conventional DQM eliminating the need to implement Local RBF-DQ over the entire domain but only on a portion of it. By this method, we have the advantages of DQ like simplicity, high accuracy, and low computational cost and the advantages of LMQRBF-DQ like mesh free and Dirac’s delta function properties. We demonstrate the effectiveness of the proposed methodology using Poisson and Burgers’ equations.     

Matrix displacement mappings in the numerical solution of functional and nonlinear differential equations with the tau method

The use of matrix displacement mappings reduces most matrix operations required in the construction of an approximate solution of a functional or differential equation by means of Ortiz' formulation of the Tau method to index shifts. The coefficient vector of the approximate solution is defined implicitly by a very sparse system of linear algebraic equations. The contributions of the differential or functional operator, and of the supplementary conditions of the problem (initial, boundary, or multipoint conditions) are treated with a single and versatile procedure of remarkable simplicity, which can be easily implemented in a computer. We give two nontrivial examples on the application of this approach: the first is a nonlinear boundary value problem with a continuous locus of singular points and multiple solutions, where stiffness is present, the second is a functional differential equation arising in analytic number theory. In both cases we obtain results of nigh 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.     

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.     

On the practical use of the spectral homotopy analysis method and local linearisation method for unsteady boundary-layer flows caused by an impulsively stretching plate

This paper expands the ideas of the spectral homotopy analysis method to apply them, for the first time, on non-linear partial differential equations. The spectral homotopy analysis method (SHAM) is a numerical version of the homotopy analysis method (HAM) which has only been previously used to solve non-linear ordinary differential equations. In this work, the modified version of the SHAM is used to solve a partial differential equation (PDE) that models the problem of unsteady boundary layer flow caused by an impulsively stretching plate. The robustness of the SHAM approach is demonstrated by its flexibility to allow linear operators that are partial derivatives with variable coefficients. This is seen to significantly improve the convergence and accuracy of the method. To validate accuracy of the the present SHAM results, the governing PDEs are also solved using a novel local linearisation technique coupled with an implicit finite difference approach. The two approaches are compared in terms of accuracy, speed of convergence and computational efficiency.     

Exact and approximation product solutions form of heat equation with nonlocal boundary conditions using Ritz–Galerkin method with Bernoulli polynomials basis

In this article, a new method is introduced for finding the exact solution of the product form of parabolic equation with nonlocal boundary conditions. Approximation solution of the present problem is implemented by the Ritz–Galerkin method in Bernoulli polynomials basis. The properties of Bernoulli polynomials are first presented, then Ritz–Galerkin method in Bernoulli polynomials is used to reduce the given differential equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the techniques presented in this article for finding the exact and approximation solutions. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1143–1158, 2017