Structural Engineering and Mechanics

Techno-Press (테크노프레스)
권호 표지
DOI QR코드

DOI QR Code

A function space approach to study rank deficiency and spurious modes in finite elements

  • Sangeeta, K. (CSIR Centre for Mathematical Modelling and Computer Simulation (C-MMACS)) ;
  • Mukherjee, Somenath (Structures Division, National Aerospace Laboratories) ;
  • Prathap, Gangan (CSIR Centre for Mathematical Modelling and Computer Simulation (C-MMACS))
  • Received : 2005.02.24
  • Accepted : 2005.08.01
  • Published : 2005.11.30
⟨ Previous Next ⟩

Abstract

Finite elements based on isoparametric formulation are known to suffer spurious stiffness properties and corresponding stress oscillations, even when care is taken to ensure that completeness and continuity requirements are enforced. This occurs frequently when the physics of the problem requires multiple strain components to be defined. This kind of error, commonly known as locking, can be circumvented by using reduced integration techniques to evaluate the element stiffness matrices instead of the full integration that is mathematically prescribed. However, the reduced integration technique itself can have a further drawback - rank deficiency, which physically implies that spurious energy modes (e.g., hourglass modes) are introduced because of reduced integration. Such instability in an existing stiffness matrix is generally detected by means of an eigenvalue test. In this paper we show that a knowledge of the dimension of the solution space spanned by the column vectors of the strain-displacement matrix can be used to identify the instabilities arising in an element due to reduced/selective integration techniques a priori, without having to complete the element stiffness matrix formulation and then test for zero eigenvalues.

Keywords

References

  1. Bathe, K.J. (2003), Finite Element Procedures, Prentice Hall
  2. Belytschko, T., Liu, W.K. and Moran, B. (2000), Nonlinear Finite Elements for Continua and Structures, Wiley
  3. Cook, R.D., Malkus, D.S. and Plesha, M.E. (1989), Concepts and Applications of Finite Element Analysis, Third Ed., John Wiley, New York
  4. Doherty, W.P., Wilson, E.L. and Taylor, R.L. (1969), 'Stress analysis of axi-symmetric solids utilizing higher order quadrilateral finite elements', Structural and Material Research SESM Report 63-9, Department of Civil Eng., University of California, Berkeley, California
  5. Mukherjee, S. and Prathap, G (2001), 'Analysis of shear locking in Timoshenko beam elements using the function space approach', Commun. Numer. Meth. Eng., 17, 385-393 https://doi.org/10.1002/cnm.413
  6. Mukherjee, S. and Prathap, G (2002), 'Analysis of delayed convergence in the three-noded Timoshenko beam element using the function space approach', Sadhana, 27(5), 507-526 https://doi.org/10.1007/BF02703292
  7. Prathap, G (1996), 'Finite element analysis and the stress correspondence paradigm', Sadhana, 21(5), 525-546 https://doi.org/10.1007/BF02744102
  8. Prathap, G and Mukherjee, S. (2003), 'The engineer grapples with Theorem 1.1 and Lemma 6.3 of Strang and Fix', Current Science, 85(7), 989-994
  9. Prathap, G (1993), The Finite Element Method in Structural Mechanics, Kluwer Academic Press, Dordrecht
  10. Sangeeta, K., Mukherjee, S. and Prathap, G (2003), 'Generalization of projection theorem for finite element analysis', Research report C-MMACS: CM 0305
  11. Strang, G and Fix, G.J. (1966), An Analysis of the Finite Element Method, Prentice-Hall Series in Automatic Computation, Prentice-Hall, Englewood Cliffs, N.J
  12. Strang, G (2003), Introduction to Linear Algebra, Wellesley-Cambridge Press
  13. Zienkiewicz, O.C., Taylor, R. and Too, J. (1971), 'Reduced integration techniques in the general analysis of plates and shells', Int. J. Num. Meth. Eng., 3, 275-320 https://doi.org/10.1002/nme.1620030211

Cited by

  1. A technique to avoid aspect-ratio locking in QUAD8 element for extremely large aspect-ratios vol.37, pp.6, 2005, https://doi.org/10.12989/sem.2011.37.6.633