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.

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.

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.

Based on the asymptotic solution for predicted natural frequencies of a two‐dimensional elastodynamic problem from the finite element analysis, presents the concept of the asymptotic error, which is an approximate error but tends to the exact error when the characteristic length of elements approaches zero, and a practical error estimator. The present practical error estimator contains two criteria: one is the error estimator criterion, the other the finite element mesh design criterion. Using this practical error estimator, not only can the accuracy of a finite element solution for natural frequencies of a two‐dimensional elastodynamic problem be directly evaluated without any further finite element calculation, but also a new target finite element mesh for the desired accuracy of solution can be immediately designed from the relevant information of an original finite element solution. Generally, for the purpose of designing a new target finite element mesh, this original finite element solution is obtainable from a very coarse mesh of a few elements and usually does not satisfy the accuracy requirement. Since the new target finite element mesh could result in a finite element solution with a desire accuracy, the finite element solution so obtained can be used for a structural design in engineering practice. The related numerical results from vibration problems of three representative plates of different shapes under plane stress conditions have demonstrated the correctness and applicability of the present practical error estimator.

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.

We describe a visual method—stress band plots—for displaying the stress solution within a two‐dimensional finite element mesh. The stress band plots differ from conventional stress contour plots because stress band plots display unaveraged stresses (the stresses are computed directly from the solution variables) and stress discontinuities in the finite element solution are directly displayed. Stress band plots are useful in judging the accuracy of a finite element solution, in the comparison of different finite element solutions and during mesh refinement. These uses are demonstrated in an axisymmetric pressure vessel analysis.

This paper gives a bibliographical review of the finite element and boundary element parallel processing techniques from the theoretical and application points of view. Topics include: theory – domain decomposition/partitioning, load balancing, parallel solvers/algorithms, parallel mesh generation, adaptive methods, and visualization/graphics; applications – structural mechanics problems, dynamic problems, material/geometrical non‐linear problems, contact problems, fracture mechanics, field problems, coupled problems, sensitivity and optimization, and other problems; hardware and software environments – hardware environments, programming techniques, and software development and presentations. The bibliography at the end of this paper contains 850 references to papers, conference proceedings and theses/dissertations dealing with presented subjects that were published between 1996 and 2002.

Contact problems are among the most difficult ones in mechanics. Due to its practical importance, the problem has been receiving extensive research work over the years. The finite element method has been widely used to solve contact problems with various grades of complexity. Great progress has been made on both theoretical studies and engineering applications. This paper reviews some of the main developments in contact theories and finite element solution techniques for static contact problems. Classical and variational formulations of the problem are first given and then finite element solution techniques are reviewed. Available constraint methods, friction laws and contact searching algorithms are also briefly described. At the end of the paper, a bibliography is included, listing about seven hundred papers which are related to static contact problems and have been published in various journals and conference proceedings from 1976.

Presents a review on implementing finite element methods on supercomputers, workstations and PCs and gives main trends in hardware and software developments. An appendix included at the end of the paper presents a bibliography on the subjects retrospectively to 1985 and approximately 1,100 references are listed.

This study aimed to overcome the challenging issues involved in providing high-precision eigensolutions. The accurate prediction of the buckling load bearing capacity under different crack damage locations, sizes and numbers, and analysing the influence mechanism of crack damage on buckling instability have become the needs of theoretical research and engineering practice. Accordingly, a finite element method was developed and applied to solve the elastic buckling load and buckling mode of curved beams with crack damage. However, the accuracy of the solution depends on the quality of mesh, and the solution inevitably introduces errors due to mesh. Therefore, the adaptive mesh refinement method can effectively optimise the mesh distribution and obtain high-precision solutions.

For the elastic buckling of circular curved beams with cracks, the section damage defect analogy scheme of a circular arc curved beam crack was established to simulate the crack size (depth), position and number. The h-version finite element mesh adaptive analysis method of the variable section Euler–Bernoulli beam was introduced to solve the elastic buckling problem of circular arc curved beams with crack damage. The optimised mesh and high-precision buckling load and buckling mode solutions satisfying the preset error tolerance were obtained.

The results of testing typical examples show that (1) the established section damage defect analogy scheme of circular arc curved beam crack can effectively realise the simulation of crack size (depth), position and number. The solution strictly satisfies the preset error tolerance; (2) the non-uniform mesh refinement in the algorithm can be adapted to solve the arbitrary order frequencies and modes of cracked cylindrical shells under the conditions of different ring wave numbers, crack positions and crack depths; and (3) the change in the buckling mode caused by crack damage is applicable to the study of elastic buckling under various curved beam angles and crack damage distribution conditions.

This study can provide a novel strategy for the adaptive mesh refinement for finite element analysis of elastic buckling of circular arc curved beams with crack damage. The adaptive mesh refinement method established in this study is fundamentally different from the conventional finite element method which employs the user experience to densify the meshes near the crack. It can automatically and flexibly generate a set of optimised local meshes by iteratively dividing the fine mesh near the crack, which can ensure the high accuracy of the buckling loads and modes. The micro-crack in curved beams is also characterised by weakening the cross-sectional stiffness to realise the characterisation of locations, depths and distributions of multiple crack damage, which can effectively analyse the disturbance behaviour of different forms of micro-cracks on the dynamic behaviour of beams.