Search results

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.

This paper aims to investigate the application of adaptive integration in element-free Galerkin methods for solving problems in structural and solid mechanics to obtain accurate reference solutions.

An adaptive quadrature algorithm which allows user control over integration accuracy, previously developed for integrating boundary value problems, is adapted to elasticity problems. The algorithm allows the development of a convergence study procedure that takes into account both integration and discretisation errors. The convergence procedure is demonstrated using an elasticity problem which has an analytical solution and is then applied to accurately solve a soft-tissue extension problem involving large deformations.

The developed convergence procedure, based on the presented adaptive integration scheme, allows the computation of accurate reference solutions for challenging problems which do not have an analytical or finite element solution.

This paper investigates the application of adaptive quadrature to solid mechanics problems in engineering analysis using the element-free Galerkin method to obtain accurate reference solutions. The proposed convergence procedure allows the user to independently examine and control the contribution of integration and discretisation errors to the overall solution error. This allows the computation of reference solutions for very challenging problems which do not have an analytical or even a finite element solution (such as very large deformation problems).

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 simulate numerical solutions of nonlinear Burgers' equation with two well‐known problems in order to verify the accuracy of the cubic B‐spline differential quadrature methods.

Cubic B‐spline differential quadrature methods have been used to discretize the Burgers' equation in space and the resultant ordinary equation system is integrated via Runge‐Kutta method of order four in time. Numerical results are compared with each other and some former results by calculating discrete root mean square and maximum error norms in each case. A matrix stability analysis is also performed by determining eigenvalues of the coefficient matrices numerically.

Numerical results show that differential quadrature methods based on cubic B‐splines generate acceptable solutions of nonlinear Burgers' equation. Constructing hybrid algorithms containing various basis to determine the weighting coefficients for higher order derivative approximations is also possible.

Nonlinear Burgers' equation is solved by cubic B‐spline differential quadrature methods.

The purpose of this paper is to develop two meshfree algorithms based on multiquadric radial basis functions (RBFs) and differential quadrature (DQ) technique for numerical simulation and to capture the shocks behavior of Burgers’ type problems.

The algorithms convert the problems into a system of ordinary differential equations which are solved by the Runge–Kutta method.

Two meshfree algorithms are developed and their stability is discussed. Numerical experiment is done to check the efficiency of the algorithms, and some shock behaviors of the problems are presented. The proposed algorithms are found to be accurate, simple and fast.

The present algorithms LRBF-DQM and GRBF-DQM are based on radial basis functions, which are new for Burgers’ type problems. It is concluded from the numerical experiments that LRBF-DQM is better than GRBF-DQM. The algorithms give better results than available literature.

The Convolution Quadrature Method (CQM) has been widely applied to solve transient elastodynamic problems because of its stability and generality. However, the CQM suffers from the problems of huge memory requirement in case of direct implementation in time domain or CPU time in case of its reformulation in Laplace domain. The purpose of this paper is to combine the CQM with the pseudo-initial condition method (PICM) to achieve a good balance between memory requirement and CPU time.

The combined methods first subdivide the whole analysis into a few sub-analyses, which is dealt with the PICM, namely, the results obtained by previous sub-analysis are used as the initial conditions for the next sub-analysis. In each sub-analysis, the time interval is further discretized into a number of sub-steps and dealt with the CQM. For non-zero initial conditions, the pseudo-force method is used to transform them into equivalent body forces. The boundary face method is employed in the numerical implementation. Three examples are analyzed. Results are compared with analytical solutions or FEM results and the results of reformulated CQM.

Results demonstrate that the computation time and the storage requirement can be reduced significantly as compared to the CQM, by using the combined approach.

The combined methods can be successfully applied to the problems of long-time dynamic response, which requires a large amount of computer memory when CQM is applied, while preserving the CQM stability. If the number of time steps is high, then the accuracy of the proposed approach can be deteriorated because of the pseudo-force method.

The purpose of the study is to provide a dataset about geometrical constructions of early Ottoman tombs for conservation studies. Thus, a proposal for the restitution phase of the damaged tombs aims to develop.

The study is composed of four phases. First, the representative plan and section drawings of early Ottoman Tombs were redrawn; second, a geometrical analysis was made, a proposal table was prepared for the restitution of the damaged tombs; and last, this table was applied to tomb examples and restitution drawings are verified with the original situation of the tombs.

Early Ottoman tombs may be interpreted through geometric shapes, including the square, circle, triangle, octagon, arsin grid and quadrature systems. The arsin grid system provides information about the position of the domes and the height of the drums and windows. Quadrature establishes the highest point of domes and entrances. The proposal table, developed from the obtained results, facilitated the identification of the original elements, including the dome, drum, window and portal. This information is crucial for conducting further studies on restitution.

The abundance and dispersed nature of tomb structures compared to other architectural designs pose challenges in their scholarly examination. The early Ottoman tombs, which experienced an increase in numbers following the Conquest of Istanbul, serve as the initial expressions and embodiments of novel architectural endeavors. Thus, the determination of design ideas of the early Ottoman tombs sheds light on Ottoman architectural practice, which has remained largely unknown and guided the conservation studies of the tombs that have lost their integrity and originality.

The purpose of this paper is to study the wave propagation in thermoelastic diffusive medium.

The present paper deals with the numerical study of wave propagation in coupled thermoelastic diffusive medium by using DQ method together with fourth‐order Runge‐Kutta method.

The paper finds solutions of displacements, temperature change and concentration.

The paper can be sued to solve non‐linear partial differential equations.

The solutions of displacements, temperature change and concentration are illustrated graphically. Numerical examples show that the method yields very good results.

Cantilever plates subject to axial flow can lose stability by flutter and properties such as viscoelasticity and laminar friction affect dynamic stability. The purpose of the present study is to investigate the dynamic stability of viscoelastic cantilever plates subject to axial flow by using the differential quadrature method.

Equation of motion of the viscoelastic plate is derived by implementing Kelvin-Voigt model of viscoelasticity and applying inverse Laplace transformation. The differential quadrature method is employed to discretize the equation of motion and the boundary conditions leading to a generalized eigenvalue problem. The solution is verified using the existing results in the literature and numerical results are given for critical flow velocities

It is observed that higher aspect ratios lead to imaginary part of third frequency becoming negative and causing single-mode flutter instability. It was found that flutter instability does not occur at low aspect ratios. Moreover the friction coefficient is found to affect the magnitude of critical flow velocity, however, its effect on the stability behaviour is minor.

The effects of various problem parameters on the dynamic stability of a viscoelastic plate subject to axial flow were established. It was shown that laminar friction coefficient of the flowing fluid increases the critical fluid velocity and higher aspect ratios lead to single-mode flutter instability. The effect of increasing damping of viscoelastic material on the flutter instability was quantified and it was found that increasing viscoelasticity can lead to divergence instability.

We study an efficient numerical method for pricing European options when the dynamics of the underlying asset are described by Levy processes. In this case. we can write a characteristic function solution for a specific Levy option model and then take its inversion numerically. Specifically we use Variance Gamma process as an example of Levy option model and consider various characteristic function representation forms of European option price such as Carr and Madan (1999), Bakshi and Madan (2000). and Lewis (2001). Fast Fourier Transform method is applied to solve the numerical inversion problem with parameters for the KOSPI 200 options data. After analysing the problems in the FFT method, we propose alternative numerical inversion method, Gaussian Quadrature. This paper reports that Gaussian Quadrature numerical inversion method with the representation form of Bakshi & Madan (2000) is more efficient and accurate than other alternatives considered in this paper.