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.

The purpose of this paper is to highlight the implementation of a recently developed weak form quadrature element method for nonlinear free vibrations of Timoshenko beams subjected to three different boundary conditions.

The design of the paper is based on considering the geometrically nonlinear effects of axial strain, bending curvature, and shear strain. Then the quadrature element formulation of the beam is introduced.

The efficiency of the method is demonstrated by a convergence study. Ratios of the nonlinear fundamental frequencies to the corresponding linear frequencies are extracted. Their variations with the ratio of amplitude to radius of gyration and the slenderness ratio are examined. The effects of the nonlinearity on higher order frequencies and mode shapes are also investigated.

The computed results show fast convergence and compare well with available literature results.

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 this paper is to parametrically investigate the vibration and damping characteristics of a functionally graded (FG) inhomogeneous and porous curved sandwich beam with a frequency-dependent viscoelastic core.

The FG material properties in this study are assumed to vary through the beam thickness by power law distribution. Additionally, FG layers have porosities, which are analyzed individually in terms of even and uneven distributions. First, the equations of motion for the free vibration of the FG curved sandwich beam were derived by Hamilton's principle. Then, the generalized differential quadrature method (GDQM) was used to solve the resulting equations in the frequency domain. Validation of the proposed FG curved beam model and the reliability of the GDQ solution was provided via comparison with the results that already exist in the literature.

A series of studies are carried out to understand the effects on the natural frequencies and modal loss factors of system parameters, i.e. beam thickness, porosity distribution, power law exponent and curvature on the vibration characteristics of an FG curved sandwich beam with a ten-parameter fractional derivative viscoelastic core material model.

This paper focuses on the vibration and damping characteristics of FG inhomogeneous and porous curved sandwich beam with frequency dependent viscoelastic core by GDQM – for the first time, to the best of the authors' knowledge. Moreover, it serves as a reference for future studies, especially as it shows that the effect of porosity distribution on the modal loss factor needs further investigation. GDQM can be useful in dynamic analysis of sandwich structures used in aerospace, automobile, marine and civil engineering applications.

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.

Accurate numerical differentiation of approximate data by methods based on Green's second identity often involves singular or nearly singular integrals over domains or their boundaries. This paper applies the finite part integration concept to evaluate such integrals and to generate suitable quadrature formulae. The weak singularity involved in first derivatives is removable; the strong singularities encountered in computing higher derivatives can be reduced. To find derivatives on or near the edge of the integration region, special treatment of boundary integrals is required. Values of normal derivative at points on the edge are obtainable by the method described. Example results are given for derivatives of analytically known functions, as well as results from finite element analysis.

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).

To study and to analyze a second order finite‐element boundary‐flux approximation using isoparametric numerical integration.

The numerical finite‐element integration is the main method used in this research. Since a domain with curved boundary is considered we apply an isoparametric approach. The lumped flux formulation is another method of approach in this paper.

This research study presents a careful analysis of the combined effect of the numerical integration and isoparametric FEM on the boundary‐flux error. Some L₂‐norm estimates are proved for the approximate solutions of the problem under consideration.

The authors offer a general study within the framework of the boundary‐flux approximation theory, which completes the results of published works in this scientific field of research.

A useful application is to employ appropriate quadrature formulae without violating the precision of the boundary‐flux FEM. The lumped mass approximation is also an important practical approach to the problem in question.

The paper presents an entire investigation in FE boundary‐flux approximation theory, in particular, elements of arbitrary degree and domains with curved boundaries. The work is addressed to the possible related fields of interest of postgraduate students and specialists in fluid mechanics and numerical analysis.

Most studies of power‐law fluids are carried out using stress‐based system of Navier‐Stokes equations; and least‐squares finite element models for vorticity‐based equations of power‐law fluids have not been explored yet. Also, there has been no study of the weak‐form Galerkin formulation using the reduced integration penalty method (RIP) for power‐law fluids. Based on these observations, the purpose of this paper is to fulfill the two‐fold objective of formulating the least‐squares finite element model for power‐law fluids, and the weak‐form RIP Galerkin model of power‐law fluids, and compare it with the least‐squares finite element model.

For least‐squares finite element model, the original governing partial differential equations are transformed into an equivalent first‐order system by introducing additional independent variables, and then formulating the least‐squares model based on the lower‐order system. For RIP Galerkin model, the penalty function method is used to reformulate the original problem as a variational problem subjected to a constraint that is satisfied in a least‐squares (i.e. approximate) sense. The advantage of the constrained problem is that the pressure variable does not appear in the formulation.

The non‐Newtonian fluids require higher‐order polynomial approximation functions and higher‐order Gaussian quadrature compared to Newtonian fluids. There is some tangible effect of linearization before and after minimization on the accuracy of the solution, which is more pronounced for lower power‐law indices compared to higher power‐law indices. The case of linearization before minimization converges at a faster rate compared to the case of linearization after minimization. There is slight locking that causes the matrices to be ill‐conditioned especially for lower values of power‐law indices. Also, the results obtained with RIP penalty model are equally good at higher values of penalty parameters.

Vorticity‐based least‐squares finite element models are developed for power‐law fluids and effects of linearizations are explored. Also, the weak‐form RIP Galerkin model is developed.

The purpose of this paper is to discuss about evaluating the integrals involving B-spline wavelet on the interval (BSWI), in wavelet finite element formulations, using Gauss Quadrature.

In the proposed scheme, background cells are placed over each BSWI element and Gauss quadrature rule is defined for each of these cells. The nodal discretization used for BSWI WFEM element is independent to the selection of number of background cells used for the integration process. During the analysis, background cells of various lengths are used for evaluating the integrals for various combination of order and resolution of BSWI scaling functions. Numerical examples based on one-dimensional (1D) and two-dimensional (2D) plane elasto-statics are solved. Problems on beams based on Euler Bernoulli and Timoshenko beam theory under different boundary conditions are also examined. The condition number and sparseness of the formulated stiffness matrices are analyzed.

It is found that to form a well-conditioned stiffness matrix, the support domain of every wavelet scaling function should possess sufficient number of integration points. The results are analyzed and validated against the existing analytical solutions. Numerical examples demonstrate that the accuracy of displacements and stresses is dependent on the size of the background cell and number of Gauss points considered per background cell during the analysis.

The current paper gives the details on implementation of Gauss Quadrature scheme, using a background cell-based approach, for evaluating the integrals involved in BSWI-based wavelet finite element method, which is missing in the existing literature.