0
Research Papers

Fast and Robust Full-Quadrature Triangular Elements for Thin Plates/Shells With Large Deformations and Large Rotations

[+] Author and Article Information
Hui Ren

MSC Software Corporation,
201 Depot Street,
Suite 100,
Ann Arbor, MI 48104
e-mail: hui.ren@mscsoftware.com

Manuscript received April 30, 2014; final manuscript received March 24, 2015; published online June 10, 2015. Assoc. Editor: Javier Garcia de Jalon.

J. Comput. Nonlinear Dynam 10(5), 051018 (Sep 01, 2015) (13 pages) Paper No: CND-14-1112; doi: 10.1115/1.4030212 History: Received April 30, 2014; Revised March 24, 2015; Online June 10, 2015

A new formulation for plates/shells with large deformations and large rotations is derived from the principles of continuum mechanics and calculated using the absolute nodal coordinate formulation (ANCF) techniques. A class of triangular elements is proposed to discretize the plate/shell formulation, which does not suffer from shear locking or membrane locking issue, and full quadrature can be performed to evaluate the integrals of each element. The adaptability of triangular elements enables the current approach to be applied to plates and shells with complicated shapes and variable thicknesses. The discretized mass matrix is constant, and the elastic force and stiffness matrix are polynomials of the generalized coordinates with constant coefficients. All the coefficients can be evaluated accurately beforehand, and numerical quadrature is not required in each time step of the simulation, which makes the current approach superior in numerical efficiency to most other approaches. The accuracy, robustness, and adaptability of the current approach are validated using both finite element benchmarks and multibody system standard tests.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.

References

MacNeal, R., 1978, “A Simple Quadrilateral Shell Element,” Comput. Struct., 8(2), pp. 175–183. [CrossRef]
Hughes, T. J. R., and Liu, W. K., 1981, “Nonlinear Finite Element Analysis of Shells: Part I. Three Dimensional Shells,” Comput. Methods Appl. Mech. Eng., 26(3), pp. 331–362. [CrossRef]
Belytschko, T., and Tsay, C. S., 1983, “A Stabilization Procedure for the Quadrilateral Plate Element With One-Point Quadrature,” Int. J. Numer. Methods Eng., 19(3), pp. 405–419. [CrossRef]
Belytschko, T., Wong, B. L., and Chiang, H. Y., 1992, “Advances in One-Point Quadrature Shell Elements,” Comput. Methods Appl. Mech. Eng., 96(1), pp. 93–107. [CrossRef]
Belytschko, T., Liu, W. K., and Moran, B., 2000, Non-Linear Finite Elements for Continua and Structures, John Wiley and Sons, New York.
Crisfield, M. A., 1991, Non-Linear Finite Element Analysis of Solids and Structures, Vol. 2, Wiley, New York, pp. 260–307.
Dmitrochenko, O. N., and Pogorelov, D. Y., 2003, “Generalization of Plate Finite Elements for Absolute Nodal Coordinate Formulation,” Multibody Syst. Dyn., 10(1), pp. 17–43. [CrossRef]
Dmitrochenko, O., 2008, “Finite Elements Using Absolute Nodal Coordinates for Large-Deformation Flexible Multibody Dynamics,” J. Comput. Appl. Math., 215(2), pp. 368–377. [CrossRef]
Shabana, A. A., 2008, Computational Continuum Mechanics, Cambridge University, Cambridge, pp. 272–275. [CrossRef]
Sedov, L. I., 1996, Mechanics of Continuous Media, 3rd ed., World Scientific, Singapore.
Barber, J. R., 2002, Elasticity, 2nd ed., Kluwer Academic, Dordrecht, p. 36.
Edelsbrunner, H., 2000, “Triangulations and Meshes in Computational Geometry,” Acta Numer., 9, pp. 1–81. [CrossRef]
Delaunay, B., 1934, “Sur la Sphère Vide,” Izv. Akad. Nauk SSSR, Otd. Mat. Estestv. Nauk, 7, pp. 793–800.
Zienkiewicz, O. C., Taylor, R. L., and Zhu, J. Z., 2005, The Finite Element Method: Its Basis and Fundamentals, 6th ed., Elsevier Butterworth-Heinemann, MA, p. 119.
Sugiyama, H., and Yamashita, H., 2011, “Spatial Joint Constraints for the Absolute Nodal Coordinate Formulation Using the Non-Generalized Immediate Coordinates,” Multibody Syst. Dyn., 26(1), pp. 15–36. [CrossRef]
MacNeal, R., and Harder, R. L., 1985, “A Proposed Standard Set of Problems to Test Finite Element Accuracy,” Finite Elem. Anal. Des., 1(1), pp. 3–20. [CrossRef]
Blevins, R. D., 1979, Formulas for Natural Frequency and Mode Shape, R. E. Krieger Publishing, Malabar, p. 240.
Den Hartog, J. P., 1952, Advanced Strength of Materials, Dover, New York, p. 128.
Taber, L. A., 1982, “Large Deflection of a Fluid-Filled Spherical Shell Under a Point Load,” ASME J. Appl. Mech., 49(1), pp. 121–128. [CrossRef]
Simo, J. C., Fox, D. D., and Rifai, M. S., 1990, “On a Stress Resultant Geometrically Exact Shell Model. Part III: Computational Aspects of the Nonlinear Theory,” Comput. Methods Appl. Mech. Eng., 79(1), pp. 21–70. [CrossRef]
Simo, J. C., Rifai, M. S., and Fox, D., 1990, “On a Stress Resultant Geometrically Exact Shell Model. Part IV: Variable Thickness Shells With Through-the-Thickness Stretching,” Comput. Methods Appl. Mech. Eng., 81(1), pp. 91–126. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

Three types of triangular elements discussed in this work: (a) linear strain element; (b) Lagrange quadratic strain element; and (c) Hermite quadratic strain element

Grahic Jump Location
Fig. 2

The elements in the membrane and bending plate patch tests

Grahic Jump Location
Fig. 3

The (a) schematic and (b) tessellation of the twisted beam

Grahic Jump Location
Fig. 4

Time histories of the displacements at the tip under (a) the out of plane force and (b) the in-plane force, calculated using the linear strain element approach (dotted), the Hermite quadratic element approach (dashed), and the rectangular thin shell elements in Ref. [7] (dashed-dotted); and the static benchmark results (solid) are also presented for comparison

Grahic Jump Location
Fig. 5

(a) Configuration and (b) tessellation of the Scordelis–Lo roof

Grahic Jump Location
Fig. 6

(a) Nonlinear displacement distribution and (b) vertical displacement field of the Scordelis–Lo roof

Grahic Jump Location
Fig. 7

Tessellation of the circular plate using (a) 600 linear strain elements and (b) 216 quadratic strain elements

Grahic Jump Location
Fig. 8

The mode shapes of the first six modes of a free circular plate

Grahic Jump Location
Fig. 10

Equilibria of the clamped circular plate under (a) a uniform pressure and (b) a concentrated force at the center

Grahic Jump Location
Fig. 9

Deflections of a clamped circular plate under (a) a uniform pressure and (b) a concentrated force, calculated from the linear strain elements (solid), and the linear plate theory (dashed)

Grahic Jump Location
Fig. 11

(a) Schematic of a hemisphere under a concentrated load and (b) tessellation of the hemisphere

Grahic Jump Location
Fig. 12

(a) The deflection-load response curve of the rubber hemisphere and (b) the collapsed configuration of the hemisphere

Grahic Jump Location
Fig. 13

The (a) x, (b) y, and (c) z tip motions of the plate pendulum, calculated using 5 × 5 × 2 linear strain triangular elements (dashed), 3 × 3 × 2 Hermite quadratic strain triangular elements (dashed-dotted), and 4 × 4 rectangular thin plate elements in Ref. [7] (dotted), and the results are compared with those calculated using 16 × 16 rectangular thin plate elements (solid)

Grahic Jump Location
Fig. 14

The (a) x, (b) y, and (c) z tip motions of the plate pendulum, calculated by 16 × 16 × 2 linear strain triangular elements (dashed), 8 × 8 × 2 Hermite quadratic strain triangular elements (dashed-dotted), and 8 × 8 rectangular thin plate elements in Ref. [7] (dotted), and the results are compared with those calculated using 16 × 16 rectangular thin plate elements (solid)

Grahic Jump Location
Fig. 15

The configurations of the flexible plate pendulum at time = 0, 0.2, 0.3, 0.4, 0.5 s

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In