Search results

In this work, the authors aim to employ the so-called linked-interpolation concept already tested on beam and quadrilateral plate finite elements in the design of displacement-based higher-order triangular plate finite elements and test their performance.

Starting from the analogy between the Timoshenko beam theory and the Mindlin plate theory, a family of triangular linked-interpolation plate finite elements of arbitrary order are designed. The elements are tested on the standard set of examples.

The derived elements pass the standard patch tests and also the higher-order patch tests of an order directly related to the order of the element. The lowest-order member of the family of developed elements still suffers from shear locking for very coarse meshes, but the higher-order elements turn out to be successful when compared to the elements from literature for the problems with the same total number of the degrees of freedom.

The elements designed perform well for a number of standard benchmark tests, but the well-known Morley's skewed plate example turns out to be rather demanding, i.e. the proposed design principle cannot compete with the mixed-type approach for this test. Work is under way to improve the proposed displacement-based elements by adding a number of internal bubble functions in the displacement and rotation fields, specifically chosen to satisfy the basic patch test and enable a softer response in the bench-mark test examples.

A new family of displacement-based higher-order triangular Mindlin plate finite elements has been derived. The higher-order elements perform very well, whereas the lowest-order element requires improvement.

A simple shape-free high-order hybrid displacement function element method is presented for precise bending analyses of Mindlin–Reissner plates. Three distortion-resistant and locking-free eight-node plate elements are proposed by utilizing this method.

This method is based on the principle of minimum complementary energy, in which the trial functions for resultant fields are derived from two displacement functions, F and f, and satisfy all governing equations. Meanwhile, the element boundary displacements are determined by the locking-free arbitrary order Timoshenko’s beam functions. Then, three locking-free eight-node, 24-DOF quadrilateral plate-bending elements are formulated: HDF-P8-23β for general cases, HDF-P8-SS1 for edge effects along soft simply supported (SS1) boundary and HDF-P8-FREE for edge effects along free boundary.

The proposed elements can pass all patch tests, exhibit excellent convergence and possess superior precision when compared to all other existing eight-node models, and can still provide good and stable results even when extremely coarse and distorted meshes are used. They can also effectively solve the edge effect by accurately capturing the peak value and the dramatical variations of resultants near the SS1 and free boundaries. The proposed eight-node models possess potential in engineering applications and can be easily integrated into commercial software.

This work presents a new scheme, which can take the advantages of both analytical and discrete methods, to develop high-order mesh distortion-resistant Mindlin–Reissner plate-bending elements.

As one of the earliest high-level programming languages, Fortran with easy accessibility and computational efficiency is widely used in the engineering field. The purpose of this paper is to present a Fortran implementation of isogeometric analysis (IGA) for thin plate problems.

IGA based on non-uniform rational B-splines (NURBS) offers exact geometries and is more accurate than finite element analysis (FEA). Unlike the basis functions in FEA, NURBS basis functions are non-interpolated. Hence, the penalty method is used to enforce boundary conditions.

Several thin plate examples based on the Kirchhoff-Love theory were illustrated to demonstrate the accuracy of the implementation in contrast with analytical solutions, and the efficiency was validated in comparison with another open method.

A Fortran implementation of NURBS-based IGA was developed to solve Kirchhoff-Love plate problems. It easily obtained high-continuity basis functions, which are necessary for Kirchhoff formulation. In comparison with theoretical solutions, the numerical examples demonstrated higher accuracy and faster convergence of the Fortran implementation. The Fortran implementation can well solve the time-consuming problem, and it was validated by the time-consumption comparison with the Matlab implementation. Due to the non-interpolation of NURBS, the penalty method was used to impose boundary conditions. A suggestion of the selection of penalty coefficients was given.

Gives a bibliographical review of the finite element methods (FEMs) applied for the linear and nonlinear, static and dynamic analyses of basic structural elements from the theoretical as well as practical points of view. The range of applications of FEMs in this area is wide and cannot be presented in a single paper; therefore aims to give the reader an encyclopaedic view on the subject. The bibliography at the end of the paper contains 2,025 references to papers, conference proceedings and theses/dissertations dealing with the analysis of beams, columns, rods, bars, cables, discs, blades, shafts, membranes, plates and shells that were published in 1992‐1995.

This paper aims to describe the formulation of a displacement-based triangular membrane finite element with true drilling rotational degree of freedom (DOF).

The presented formulation incorporates the true drilling rotation provided by continuum mechanics into the displacement field by way of using the polynomial interpolation. Unlike the linked interpolation, that uses a geometric transformation between displacement and vertex rotations, in this work, the interpolation of the displacement field in terms of nodal drilling rotations is obtained following an unusual approach that does not imply any presumed geometric transformation.

New relationship linking the mid-side normal displacement to corner node drilling rotations is derived. The resulting new element with true drilling rotation is compatible and does not include any problem-dependent parameter that may influence the results. The spurious zero-energy mode is stabilized in a careful way that preserves the true drilling rotational degrees of freedom (DOFs).

Several works dealing with membrane elements with vertex rotational DOFs have been published with improved convergence rate, however, owing to the need for incorporating rotations in the finite element meshes involving solids, shells and beam elements, having finite elements with true drilling rotational DOFs is more appreciated.

This paper aims to focus on the temporal and spatial fourth‐order finite volume discretization of the incompressible form of the Navier‐Stokes equations on structured uniform grids. The main purpose of the paper is to assess the accuracy enhancement with the inclusion of a high‐order reconstruction of the point‐wise velocity field on a fourth‐order accurate numerical scheme for the solution of the unsteady incompressible Navier‐Stokes equations.

The present finite volume method uses a fractional time‐step for decoupling velocity and pressure. A Runge‐Kutta integration scheme is implemented for integrating the momentum equation along with a polynomial interpolation and Simpson formula for space‐integration. The formulation is based on step‐by‐step de‐averaging process applied to the velocity field.

The reconstruction of the point‐wise velocity field on a higher‐order basis is essential to obtain solutions that effectively stand for a fourth‐order approximation of the point‐wise one. Results are provided for the Taylor vortex decay problem and for co‐ and counter‐rotating vortices to assess the increase in accuracy promoted by the inclusion of the high‐order de‐averaging procedure.

High‐order reconstruction of the point‐wise velocity field should be considered in high‐order finite volume methods for the solution of the unsteady incompressible form of the Navier‐Stokes equations on structured grids.

The inclusion of a high‐order reconstruction of the point‐wise velocity field is a simple and effective method of enhancing the accuracy of a finite volume code for the computational fluid dynamics analysis.

The paper develops an improved version of a fourth‐order accurate finite volume projection method with the inclusion of a high‐order reconstruction step.

Gives introductory remarks about chapter 1 of this group of 31 papers, from ISEF 1999 Proceedings, in the methodologies for field analysis, in the electromagnetic community. Observes that computer package implementation theory contributes to clarification. Discusses the areas covered by some of the papers ‐ such as artificial intelligence using fuzzy logic. Includes applications such as permanent magnets and looks at eddy current problems. States the finite element method is currently the most popular method used for field computation. Closes by pointing out the amalgam of topics.

The purpose of this paper is to investigate the use of the finite element (FE) method in calculating accurate field gradients to achieve good estimates of aberration coefficients in accelerator tubes.

Gradient fields up to third order are calculated using FE calculations with higher order basis functions. A commercial code with the ability to use high order basis functions was embedded within the MathCad system which provides a convenient interface for evaluating the aberration coefficients. Analytic and realistic models are used to test the methodology.

It is shown that the ability to use higher order FEs achieves sufficient accuracy and smoothness to allow the third order aberration coefficients to be determined with confidence.

The results demonstrate the importance of using higher order basis FEs in determining optical properties of accelerator tubes and other particle optical devices.

Describes a new finite element scheme for the large‐scale analysis of compressible and incompressible viscous flows. The scheme is based on a combined compressible‐ incompressible Galerkin least‐squares (GLS) space‐time variational formulation. Three‐ dimensional unstructured meshes are employed, with piecewise‐constant temporal interpolation, local time‐stepping for steady flows, and linear continuous spatial interpolation in all the variables. The scheme incorporates automatic adaptive mesh refinement, with a choice of various error indicators. It is implemented on a distributed‐memory parallel computer, and includes an automatic load‐balancing procedure. Demonstrates the ability to solve both compressible and incompressible viscous flow problems using the parallel adaptive framework via numerical examples. These include Mach 3 flow over a flat plate, and a divergence‐free buoyancy‐driven flow in a cavity. The latter is a model for the steady melt flow in a Czochralski crystal growth process.

Modelling accurately the transient behaviour of natural convection flow in enclosures been a challenging task because of a variety of numerical errors which have limited achieving the higher order temporal accuracy. A fourth-order accurate finite difference method in both space and time is proposed to overcome these numerical errors and accurately model the transient behaviour of natural convection flow in enclosures using vorticity–streamfunction formulation.

Fourth-order wide stencil formula with appropriate one-sided difference extrapolation technique near the boundary is used for spatial discretisation, and classical fourth-order Runge–Kutta scheme is applied for transient term discretisation. The proposed method is applied on two transient case studies, i.e. convection–diffusion of a Gaussian Pulse and Taylor Vortex flow having analytical solution.

Error magnitude comparison and rate of convergence analysis of the proposed method with these analytical solutions establish fourth-order accuracy and prove the ability of the proposed method to truly capture the transient behaviour of incompressible flow. Also, to test the transient natural convection flow behaviour, the algorithm is tested on differentially heated square cavity at high Rayleigh number in the range of 10³-10⁸, followed by studying the transient periodic behaviour in a differentially heated vertical cavity of aspect ratio 8:1. An excellent comparison is obtained with standard benchmark results.

The developed method is applied on 2D enclosures; however, the present methodology can be extended to 3D enclosures using velocity–vorticity formulations which shall be explored in future.

The proposed methodology to achieve fourth-order accurate transient simulation of natural convection flows is novel, to the best of the authors’ knowledge. Stable fourth-order vorticity boundary conditions are derived for boundary and external boundary regions. The selected case studies for comparison demonstrate not only the fourth-order accuracy but also the considerable reduction in error magnitude by increasing the temporal accuracy. Also, this study provides novel benchmark results at five different locations within the differentially heated vertical cavity of aspect ratio 8:1 for future comparison studies.