Search results

This paper aims to present a numerical solution of non‐linear Burger's equation using differential quadrature method based on sinc functions.

Sinc Differential Quadrature Method is used for space discretization and four stage Runge‐Kutta algorithm is used for time discretization. A rate of convergency analysis is also performed for shock‐like solution. Numerical stability analysis is performed.

Sinc Differential Quadrature Method generates more accurate solutions of Burgers' equation when compared with the other methods.

This combination, Sinc Differential Quadrature and Runge‐Kutta of order four, has not been used to obtain numerical solutions of Burgers' equation.

The purpose of this paper is to present a general formulation of the quadrature element method (QEM). The method is then used to investigate the free vibration of functionally graded (FG) beams with general boundary conditions and different variations of material properties.

The quadrature elements with arbitrary number of nodes and nodal distributions are established on the basis of two types of FG Timoshenko beam theories. One called TBT-1 takes the cross-sectional rotation as the unknown function and the other called TBT-2 uses the transverse shear strain as the unknown function. Explicit formulas are provided via the help of the differential quadrature (DQ) rule and thus the elements can be implemented adaptively with ease.

The suitability and computational efficiency of the proposed quadrature elements for the vibration analysis of FG beams are demonstrated. The convergence rate of the proposed method is high. The elements are shear-locking free and can yield accurate solutions with a small number of nodes for both thin and moderately thick beams. The performance of the element based on TBT-1 is better than the one based on TBT-2.

The present QEM is different from the existing one which exclusively uses Gauss–Lobatto–Legendre (GLL) nodes and GLL quadrature and thus is more general. The element nodes can be either the same or different from the integration points, making the selection of element nodes more flexible. Presented data are accurate and may be a reference for other researchers to develop new numerical methods. The QEM may be also useful in multi-scale modeling and in the analysis of civil infrastructures.

The purpose of the present paper is to investigate the nonconservative instability of a single-walled carbon nanotube (SWCNT) with an added mass through nonlocal theories. The governing equations are discretized by means of the differential quadrature (DQ) rules, as introduced by Bellman and Casti. DQ rules have been largely used in engineering and applied sciences. Recently, they were applied to enhance some numerical schemes, such as step-by-step integration schemes and Picard-like numerical schemes.

In the present paper, the DQ rules are used to investigate the nonconservative instability of a SWCNT through nonlocal theories.

To show the sensitivity of the SWCNT to the values of added mass and the influence of nonlocal parameter on the fundamental frequencies values, some numerical examples have been performed and discussed. Yet, the effect of the different boundary conditions on the instability behaviour has been investigated. The validity of the present model has been confirmed by comparing some results against the ones available in literature.

Applying the nonlocal elasticity theory, this paper presents a re-formulation of Hamilton’s principle for the free vibration analysis of a uniform Euler–Bernoulli nanobeam. The main purpose of this paper is to investigate the free vibration response of an SWCNT with attached mass and for various values of small scale effects.

This study aims to demonstrate the numerical application of differential quadrature (DQ) methods and show the experimental application of free vibration analysis of fiber-metal laminated composite (FML) plates with various boundary conditions.

The FMLs are hybrid structures consisting of fiber-reinforced polymer matrix composites such as carbon, glass, aramid and different metal sheets, and are currently widely used in the automobile, aircraft and aerospace industries. Thus, free vibration analysis of these hybrid materials is necessary for the design process. The governing equations of motion are derived based on the classical plate theory. The DQ, generalized DQ (GDQ) and harmonic DQ (HDQ) differential quadrature methods have been used to solve the governing equations of an FML composite plate numerically. The accuracy and convergence of the numerical model have been verified by comparing the results available in the published literature with the results obtained from these methods. Moreover, an experimental procedure has been performed in order to compare the results against those of the numerical methods.

It is noteworthy that a high degree of similarity and accuracy was observed between the numerical results obtained by the DQ methods and the experimental results. Thus, the present study validates the applicability of the DQ methods for designing the FML composite plates.

In this study, the advantages of the DQ methods have been demonstrated differently from previous studies on the vibration analysis of the FML plates.

The purpose of this paper is to present some modifications in the spline‐based differential quadrature method (DQM), in order to accelerate the convergence of the method. The improvements are explained and examined by the examples of the free vibration of conical shells. The composite laminated shell, as well as isotropic one, are taken under consideration.

To determine weighting coefficients for the DQM, the spline interpolation with non‐standard definitions of the end conditions is used. One of these definitions combines natural and not‐a‐knot end conditions, while the other one uses the boundary conditions for considered problem as the end conditions. The weighting coefficients can be determined by solving set of equations arising from spline interpolation.

It is shown that the proposed modifications significantly improve the convergence of the method, especially when the boundary conditions are introduced at the stage of the computation of the weighting coefficients. Unfortunately, the use of this approach is limited to some types of boundary conditions.

The paper describes development of the modified spline interpolation dedicated to DQM and examines the possibility of building boundary conditions into the weighting coefficients at the stage of the computation of these coefficients.

The purpose of this paper is to develop a general mathematical model for the evaluation of the bending and vibration responses of the skew sandwich composite plate using higher-order shear deformation theory. The sandwich structural components are highly preferable in modern engineering application because of their desirable structural advantages despite the manufacturing and analysis complexities. The present model is developed to solve the bending and vibration problem of the skew sandwich composite plate with adequate accuracy numerically in the absence of the experimental analysis.

The skew sandwich composite plate structure is modelled in the present analysis by considering laminated face sheet in conjunction with isotropic and/or orthotropic core numerically with the help of the higher-order mathematical model. Further, the responses are computed numerically with the help of in-house computer code developed in matrix laboratory (MATLAB) environment in conjunction with finite element (FE) steps. The system governing equations are derived via variational technique for the computation of the static and the frequency responses.

The skew sandwich composite plate is investigated using the higher-order kinematic model where the transverse displacement through the thickness is considered to be linear. The convergence and the validation study of the bending and the frequency values of the sandwich structure indicate the necessary accuracy. Further, the current model has been used to highlight the applicability of the higher-order kinematics for the evaluation of the sandwich structural responses (frequency and static deflections) for different design parameters.

In the present paper, the bending and the vibration responses of the skew sandwich composite plate are analysed numerically using the equivalent single-layer higher-order kinematic theory for the isotropic and the orthotropic core numerically with the help of isoparametric FE steps. Finally, it is understood that the present model is capable of solving the sandwich structural responses with less computation cost and adequate accuracy.

The purpose of this paper is to present a computational technique based on Newton–Cotes quadrature rule for solving fractional order differential equation.

The numerical method reduces initial value problem into a system of algebraic equations. The method presented here is also applicable to non-linear differential equations. To deal with non-linear equations, a recursive sequence of approximations is developed using quasi-linearization technique.

The method is tested on several benchmark problems from the literature. Comparison shows the supremacy of proposed method in terms of robust accuracy and swift convergence. Method can work on several similar types of problems.

It has been demonstrated that many physical systems are modelled more accurately by fractional differential equations rather than classical differential equations. Therefore, it is vital to propose some efficient numerical method. The computational technique presented in this paper is based on Newton–Cotes quadrature rule and quasi-linearization. The key feature of the method is that it works efficiently for non-linear problems.

Various time integration methods and time finite element methods have been developed to obtain the responses of structural dynamic problems, but the accuracy and computational efficiency of them are sometimes not satisfactory. The purpose of this paper is to present a more accurate and efficient formulation on the basis of the weak form quadrature element method to solve linear structural dynamic problems.

A variational principle for linear structural dynamics, which is inspired by Noble's work, is proposed to develop the weak form temporal quadrature element formulation. With Lobatto quadrature rule and the differential quadrature analog, a system of linear equations is obtained to solve the responses at sampling time points simultaneously. Computation for multi-elements can be carried out by a time-marching technique, using the end point results of the last element as the initial conditions for the next.

The weak form temporal quadrature element formulation is conditionally stable. The relation between the normalized length of element and the suggested number of integration points in one element is given by a simple formula. Results show that the present formulation is much more accurate than other time integration methods and its dissipative property is also illustrated.

The weak form temporal quadrature element formulation provides a choice with high accuracy and efficiency for solution of linear structural dynamic problems.

Solutions for the earth return mutual impedance play an important role in analyzing couplings of multi-conductor systems. Generally, the mutual impedance is approximated by Pollaczek integrals. The purpose of this paper is devising fast algorithms for calculation of this kind of improper integrals and its applications.

According to singular points, the Pollaczek integral is divided into two parts: the finite integral and the infinite integral. The finite part is computed by combining an efficient Levin method, which is implemented with a Chebyshev differential matrix. By transforming the integration path, the tail integral is calculated with help of a transformed Clenshaw–Curtis quadrature rule.

Numerical tests show that this new method is robust to high oscillation and nearly singularities. Thus, it is suitable for evaluating Pollaczek integrals. Furthermore, compared with existing method, the presented algorithm gives high-order approaches for the earth return mutual impedance between conductors over a multilayered soil with wide ranges of parameters.

An efficient truncation strategy is proposed to accelerate numerical calculation of Pollaczek integral. Compared with existing algorithms, this method is easier to be applied to computation of similar improper integrals, such as Sommerfeld integral.

The global methods of generalized differential quadrature (GDQ) and generalized integral quadrature (GIQ) are applied to solve three‐dimensional, incompressible, laminar boundary layer equations. The streamwise and crosswise velocity components are taken as the dependent variables. The normal velocity is obtained by integrating the continuity equation along the normal direction where the integral is approximated by GIQ approach with high order of accuracy. All the spatial derivatives are discretized by a GDQ scheme. After spatial discretization, the resultant ordinary differential equations are solved by the 4‐stage Runge—Katta scheme. Application of GDQ—GIQ approach to a test problem demonstrated that accurate numerical results can be obtained using just a few grid points.