Search results

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.

The numerical quadrature of the stiffness matrices and force vectors is a major factor in the accuracy and efficiency of the finite element methods. Even in the analysis of linear problems, the use of too many quadrature points results in a phenomenon called locking whereas the use of insufficient quadrature points results in a phenomenon called spurious singular mode. Therefore, efficient numerical quadrature schemes for structural dynamics are developed. It is expected that these improved finite elements can be used more reliably in many structural applications. The proposed developed quadrature schemes for the continuum and shell elements are straightforward and are modularized so that many existing finite element computer codes can be easily modified to accommodate the proposed capabilities. This should prove of great benefit to many computer codes which are currently used in structural engineering applications.

Based on a lumped mass model and an incremental iteration method, an efficient simultaneous iteration procedure is developed for the finite element solution of the enthalpy model. This procedure uses Gauss elimination to solve the resulting algebraic equation system. A one‐point quadrature program based on the isoparametric quadrilateral element is incorporated for the calculation of the heat conductance matrix, leading to a significant reduction of the computation time.

A 4‐node flat shell quadrilateral finite element with 6 degrees of freedom per node, denoted QC5D‐SA, is presented. The element is an assembly of a modification of the drilling degree of freedom membrane presented by Ibrahimbegovic et al., and the assumed strain plate element presented by Bathe and Dvorkin. The part of the stiffness matrix associated with in—plane displacements and rotations is integrated over the element domain by a modified 5‐point reduced integration scheme, resulting in greater efficiency without the sacrifice of rank sufficiency. The scheme produces a soft higher order deformation mode which increases numerical accuracy. A large number of standard benchmark problems are analyzed. Some examples show that the effectiveness of a previously proposed “membrane locking correction” technique is significantly reduced when employing distorted elements. However, the element is shown to be generally accurate and in many cases superior to existing elements.

This paper is concerned with new formulations of local meshfree and finite element numerical methods, for the solution of two-dimensional problems in linear elasticity.

In the local domain, assigned to each node of a discretization, the work theorem establishes an energy relationship between a statically admissible stress field and an independent kinematically admissible strain field. This relationship, derived as a weighted residual weak form, is expressed as an integral local form. Based on the independence of the stress and strain fields, this local form of the work theorem is kinematically formulated with a simple rigid-body displacement to be applied by local meshfree and finite element numerical methods. The main feature of this paper is the use of a linearly integrated local form that implements a quite simple algorithm with no further integration required.

The reduced integration, performed by this linearly integrated formulation, plays a key role in the behavior of local numerical methods, since it implies a reduction of the nodal stiffness which, in turn, leads to an increase of the solution accuracy and, which is most important, presents no instabilities, unlike nodal integration methods without stabilization. As a consequence of using such a convenient linearly integrated local form, the derived meshfree and finite element numerical methods become fast and accurate, which is a feature of paramount importance, as far as computational efficiency of numerical methods is concerned. Three benchmark problems were analyzed with these techniques, in order to assess the accuracy and efficiency of the new integrated local formulations of meshfree and finite element numerical methods. The results obtained in this work are in perfect agreement with those of the available analytical solutions and, furthermore, outperform the computational efficiency of other methods. Thus, the accuracy and efficiency of the local numerical methods presented in this paper make this a very reliable and robust formulation.

Presentation of a new local mesh-free numerical method. The method, linearly integrated along the boundary of the local domain, implements an algorithm with no further integration required. The method is absolutely reliable, with remarkably-accurate results. The method is quite robust, with extremely-fast computations.

A simple beam element is developed for the solution of large deflection problems. The total Lagrangian formulation is based on the kinematic relations proposed by Reissner for finite rotations and stretching as well as shearing of plane beams. The motion is discretized by linear expansions of the global displacement components and the cross‐sectional rotation in two‐dimensional Euclidean space yielding a simple beam element with three degrees of freedom at the two nodes. The shear locking is reduced by selective integration in order to eliminate the spurious shear constraint similar to interdependent variable interpolation. The large rotation formulation is compared with two forms of moderate rotation theories which have been used in the past to develop the geometric stiffness properties for linear stability analysis of the so‐called Mindlin plate elements. The predictive value of different geometric stiffness approximations is assessed with several examples which range from the static and kinetic stability analysis of the classical Euler‐column to the large deflection problem of a clamped beam.

The first purpose of this paper is to design more accurate, efficient and robust gradient smoothing methods (GSMs) for spatial derivative approximations for computational fluid dynamics (CFD) application. The second purpose is to design an adaptive GSM-CFD solver for the compressible turbulent flows, with special focus on the shock-wave boundary layer interactions.

A new integration scheme is proposed for the node-associated GSM to improve the accuracy and robustness of the previous versions. A matrix-based algorithm and corresponding data structures are devised to improve the computational efficiency of GSM. The GSM-CFD solver is coupled with a mixed solution-based adaptive mesher to form a functional adaptive GSM-CFD solver.

The improved GSMs are insensitive to mesh qualities, and can achieve high accuracy on all kinds of hybrid meshes. The adaptive GSM-CFD solver can better capture the shock wave.

The matrix-based GSM and its corresponding data structure for improved GSM, and the development of the adaptive GSM-CFD solver for compressible turbulent flows are newly presented in this paper.

Inside the finite element framework of LAGAMINE code, the contact conditions are introduced with specific two‐node interface elements and four‐node quadrangular elements or four‐node one point quadrature elements. A non‐associated flow rule is involved for sliding unilateral contact modelling. Two methods of penalty factor computations in the penalty contact algorithms are presented. These methods are then used for contact modelling of two isothermal examples: axisymmetric tube expansion and asymmetric slab bending, the material bulk constitutive equation being isotropic and elasto‐plastic.

The constitutive relationship of a four‐node flat shell finite element with six degrees of freedom per node and a modified five‐point quadrature, previously presented by the authors, is extended to include symmetric and unsymmetric orthotropy. Through manipulation of the kinematic assumptions, provision is made for out‐of‐plane warp. A wide range of membrane and thin to moderately thick plate and shell examples are used to demonstrate the accuracy and robustness of the resulting element.

In empirical research, panel (and multinomial) probit models are leading examples for the use of maximum simulated likelihood estimators. The Geweke–Hajivassiliou–Keane (GHK) simulator is the most widely used technique for this type of problem. This chapter suggests an algorithm that is based on GHK but uses an adaptive version of sparse-grids integration (SGI) instead of simulation. It is adaptive in the sense that it uses an automated change-of-variables to make the integration problem numerically better behaved along the lines of efficient importance sampling (EIS) and adaptive univariate quadrature. The resulting integral is approximated using SGI that generalizes Gaussian quadrature in a way such that the computational costs do not grow exponentially with the number of dimensions. Monte Carlo experiments show an impressive performance compared to the original GHK algorithm, especially in difficult cases such as models with high intertemporal correlations.