Search results

Presents a procedure for obtaining an improved finite element solution of boundary problems by estimating the principle of exact displacement method in the finite element technique. The displacement field is approximated by two types of functions: the shape functions satisfying the homogeneous differential equilibrium equation and the full clamping element functions as a particular solution of the differential equation between the nodes. The full clamping functions represent the solution of the full clamping state on finite elements. An improved numerical solution of displacements, strains, stresses and internal forces, not only at nodes but over the whole finite element, is obtained without an increase of the global basis, because the shape functions are orthogonal with the full clamping functions. This principle is generally applicable to different finite elements. The contribution of introducing two types of functions based on the principle of the exact displacement method is demonstrated in the solution procedure of frame structures and thin plates.

A finite element formulation for flexure of isotropic plates based on a recent refined theory is developed. The refined theory incorporates effects of transverse shear, transverse normal stress and transverse normal strain. The Galerkin finite element method was used to develop the finite element equations for both plate bending and inplane problems. The performance of the proposed finite element model was evaluated by solving problems of uniformly loaded thick plates with different support conditions. The results of the present formulation are compared with Mindlin/Reissner and elasticity solutions.

Error indicators are introduced as part of a simple computational procedure for improving the accuracy of the finite element solutions for plate and shell problems. The procedure is based on using an initial (coarse) grid and a refined (enriched) grid, and approximating the solution for the refined grid by a linear combination of a few global approximation vectors (or modes) which are generated by solving two uncoupled sets of equations in the coarse grid unknowns and the additional degrees of freedom of the refined grid. The global approximation vectors serve as error indicators since they provide quantitative pointwise information about the sensitivity of the different response quantities to the approximation used. The three key elements of the computational procedure are: (a) use of mixed finite element models with discontinuous stress resultants at the element interfaces; (b) operator splitting, or additive decomposition of the finite element arrays for the refined grid into the sum of the coarse grid arrays and correction terms (representing the refined grid contributions); and (c) application of a reduction method through successive use of the finite element method and the classical Bubnov—Galerkin technique. The finite element method is first used to generate a few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Bubnov—Galerkin technique. The similarities between the proposed computational procedure and a preconditioned conjugate gradient (PCG) technique are identified and are exploited to generate from the PCG technique pointwise error indicators. The effectiveness of the proposed procedure is demonstrated by means of two numerical examples of an isotropic toroidal shell and a laminated anisotropic cylindrical panel.

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.

In this paper we derive a simple finite element formulation for geometrical nonlinear shell structures. The formulation bases on a direct introduction of the isoparametric finite element formulation into the shell equations. The element allows the occurrence of finite rotations which are described by two independent angles. A layerwise linear elastic material model for composites has been chosen. A consistent linearization of all equations has been derived for the application of a pure Newton method in the nonlinear solution process. Thus a quadratic convergence behaviour can be achieved in the vicinity of the solution point. Examples show the applicability and effectivity of the developed element.

Gives a bibliographical review of the error estimates and adaptive finite element methods from the theoretical as well as the application point of view. The bibliography at the end contains 2,177 references to papers, conference proceedings and theses/dissertations dealing with the subjects that were published in 1990‐2000.

This bibliography is offered as a practical guide to published papers, conference proceedings papers and theses/dissertations on the finite element (FE) and boundary element (BE) applications in different fields of biomechanics between 1976 and 1991. The aim of this paper is to help the users of FE and BE techniques to get better value from a large collection of papers on the subjects. Categories in biomechanics included in this survey are: orthopaedic mechanics, dental mechanics, cardiovascular mechanics, soft tissue mechanics, biological flow, impact injury, and other fields of applications. More than 900 references are listed.

This paper aims to develop a hybrid finite difference‐finite element method and apply it to solve the three‐dimensional energy equation in non‐isothermal fluid flow past over a tube.

To implement the hybrid scheme, the tube length is partitioned into uniform segments by choosing grid points along its length, and a plane perpendicular to the tube axis is drawn at each of the points. Subsequently, the Taylor‐Galerkin finite element technique is employed to discretize the energy equation in the planes; while the derivatives along the tube are discretized using the finite difference method.

To demonstrate the validity of the proposed numerical scheme, three‐dimensional test cases have been solved using the method. The variation of L²‐norm of the error with mesh refinement shows that the numerical solution converges to the exact solution with mesh refinement. Moreover, comparison of the computational time duration shows that the proposed method is approximately three times faster than the 3D finite element method. In the non‐isothermal fluid flow around a tube for Re=250 and Pr=0.7, the results show that the Nusselt number decreases with the increase in the tube length and, for the tube length greater than six times the tube diameter, the average Nusselt number converges to the value for the two‐dimensional case.

A hybrid finite difference‐finite element method has been developed and applied to solve the 3D transient energy equation for different test cases. The proposed method is faster, and computationally more efficient, compared with the 3D finite element method.

This study aims to provide a reliable and effective algorithm that is suitable for addressing the problems of continuous orders of frequencies and modes under different boundary conditions, circumferential wave numbers and thickness-to-length ratios of moderately thick circular cylindrical shells. The theory of free vibration of rotating cylindrical shells is of utmost importance in fields such as structural engineering, rock engineering and aerospace engineering. The finite element method is commonly used to study the theory of free vibration of rotating cylindrical shells. The proposed adaptive finite element method can achieve a considerably more reliable high-precision solution than the conventional finite element method.

On a given finite element mesh, the solutions of the frequency mode of the moderately thick circular cylindrical shell were obtained using the conventional finite element method. Subsequently, the superconvergent patch recovery displacement method and high-order shape function interpolation techniques were introduced to obtain the superconvergent solution of the mode (displacement), while the superconvergent solution of the frequency was obtained using the Rayleigh quotient computation. Finally, the superconvergent solution of the mode was used to estimate the errors of the finite element solutions in the energy norm, and the mesh was subdivided to generate a new mesh in accordance with the errors.

In this study, a high-precision and reliable superconvergent patch recovery solution for the vibration modes of variable geometrical rotating cylindrical shells was developed. Compared with conventional finite element method, under the challenging varying geometrical circumferential wave numbers, and thickness–length ratios, the optimised finite element meshes and high-precision solutions satisfying the preset error limits were obtained successfully to solve the frequency and mode of continuous orders of rotating cylindrical shells with multiple boundary conditions such as simple and fixed supports, demonstrating good solution efficiency. The existing problem on the difficulty of adapting a set of meshes to the changes in vibration modes of different orders is finally overcome by applying the adaptive optimisation.

The approach developed in this study can accurately obtain the superconvergent patch recovery solution of the vibration mode of rotating cylindrical shells. It can potentially be extended to fine numerical models and high-precision computations of vibration modes (displacement field) and solid stress (displacement derivative field) for general structural special value problems, which can be extensively applied in the field of engineering computations in the future. Furthermore, the proposed method has the potential for adaptive analyses of shell structures and three-dimensional structures with crack damage. Compared with conventional finite element methods, significant advantages can be achieved by solving the eigenvalues of structures with high precision and stability.

The paper starts with a review of constitutive equations for rubber‐like materials, formulated in the invariants of the right Cauchy—Green deformation tensor. A general framework for the derivation of the stress tensor and the tangent moduli for invariant‐based models, for both the reference and the current configuration, is presented. The free energy of incompressible rubber‐like materials is extended to a compressible formulation by adding the volumetric part of the free energy. In order to overcome numerical problems encountered with displacement‐based finite element formulations for nearly incompressible materials, three‐dimensional finite elements, based on a penalty‐type formulation, are proposed. They are characterized by applying reduced integration to the volumetric parts of the tangent stiffness matrix and the pressure‐related parts of the internal force vector only. Moreover, hybrid finite elements are proposed. They are based on a three‐field variational principle, characterized by treating the displacements, the dilatation and the hydrostatic pressure as independent variables. Subsequently, this formulation is reduced to a generalized displacement formulation. In the numerical study these formulations are evaluated. The results obtained are compared with numerical results available in the literature. In addition, the proposed formulations are applied to 3D finite element analysis of an automobile tyre. The computed results are compared with experimental data.