Research Papers

Curvature Expressions for the Large Displacement Analysis of Planar Beam Motions

[+] Author and Article Information
Yinhuan Zheng

School of Mechanical
and Electronic Engineering,
Wuhan University of Technology,
Wuhan 430070, Hubei, China
e-mail: zhengyinhuan@whut.edu.cn

Ahmed A. Shabana

Department of Mechanical
and Industrial Engineering,
University of Illinois at Chicago,
Chicago, IL 60607
e-mail: shabana@uic.edu

Dayu Zhang

National Key Laboratory of
Aerospace Flight Dynamics,
School of Astronautics,
Northwestern Polytechnical University,
Xi'an 710072, Shaanxi, China
e-mail: dyzhang@mail.nwpu.edu.cn

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received April 3, 2017; final manuscript received June 19, 2017; published online October 31, 2017. Assoc. Editor: Dan Negrut.

J. Comput. Nonlinear Dynam 13(1), 011013 (Oct 31, 2017) (12 pages) Paper No: CND-17-1149; doi: 10.1115/1.4037226 History: Received April 03, 2017; Revised June 19, 2017

While several curvature expressions have been used in the literature, some of these expressions differ from basic geometry definitions and lead to kinematic coupling between bending and shear deformations. This paper uses three different elastic force formulations in order to examine the effect of the curvature definition in the large displacement analysis of beams. In the first elastic force formulation, a general continuum mechanics approach (method 1) based on the nonlinear strain–displacement relationship is used. The second approach (method 2) is based on a classical nonlinear beam theory, in which a curvature expression consistent with differential geometry and independent of the shear deformation is used. The third elastic force formulation (method 3) employs a curvature expression that depends on the shear angle. In order to examine numerically the effect of using different curvature definitions, three different planar beam elements are used. The first element (element I) is the fully parameterized absolute nodal coordinate formulation (ANCF) shear deformable beam element. The second element (element II) is an ANCF consistent rotation-based formulation (CRBF) shear deformable beam element obtained from element I by consistently replacing the position gradient vectors by rotation parameters. The third element (element III) is a low-order bilinear ANCF/CRBF finite element in which nonzero differential geometry-based curvature definition cannot be obtained because of the low order of interpolation. Numerical results are obtained using the three elastic force formulations and the three finite elements in order to shed light on the definition of bending and shear in the large displacement analysis of beams. The results obtained in this investigation show that the use of method 2, with a penalty formulation that restricts the excessive cross section deformation, can improve significantly the convergence of the ANCF finite element.

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


Ding, J. , Wallin, M. , Wei, C. , Recuero, A. M. , and Shabana, A. A. , 2014, “ Use of Independent Rotation Field in the Large Displacement Analysis of Beams,” Nonlinear Dyn., 76(3), pp. 1829–1843. [CrossRef]
Shabana, A. A. , 2015, “ ANCF Consistent Rotation-Based Finite Element Formulation,” ASME J. Comput. Nonlinear Dyn., 11(1), p. 014502. [CrossRef]
Zheng, Y. , and Shabana, A. A. , 2017, “ A Two-Dimensional Shear Deformable ANCF Consistent Rotation-Based Formulation Beam Element,” Nonlinear Dyn., 87(2), pp. 1031–1043. [CrossRef]
Reissner, E. , 1972, “ On a One Dimensional Finite Strain Beam: The Plane Problem,” J. Appl. Math. Phys., 23(5), pp. 795–804. [CrossRef]
Reissner, E. , 1973, “ On a One Dimensional Large-Displacement, Finite Strain Beam Theory,” Stud. Appl. Math., 52(2), pp. 87–95. [CrossRef]
Reissner, E. , 1981, “ On Finite Deformations of Space-Curved Beams,” J. Appl. Math. Phys. (ZAMP), 32(6), pp. 734–744. [CrossRef]
Antman, S. S. , 1972, “ The Theory of Rods,” Handb. Phys., VIa/2, pp. 641–703.
Antman, S. S. , 1974, “ Kirchhoff's Problem for Nonlinearly Elastic Rods,” Q. Appl. Math., 32(3), pp. 221–240. [CrossRef]
Simo, J. C. , and Vu-Quoc, L. , 1986, “ On the Dynamics of Flexible Beams Under Large Overall Motions—The Plane Case: Parts I and II,” ASME J. Appl. Mech., 53(4), pp. 849–863. [CrossRef]
Omar, M. A. , and Shabana, A. A. , 2001, “ A Two-Dimensional Shear Deformable Beam for Large Rotation and Deformation Problems,” J. Sound Vib., 243(3), pp. 565–576. [CrossRef]
Dmitrochenko, O. , and Mikkola, A. , 2011, “ Digital Nomenclature Code for Topology and Kinematics of Finite Elements Based on the Absolute Nodal Co-Ordinate Formulation,” IMechE J. Multibody Dyn., 225(1), pp. 34–51.
Hu, W. , Tian, Q. , and Hu, H. Y. , 2014, “ Dynamics Simulation of the Liquid-Filled Flexible Multibody System Via the Absolute Nodal Coordinate Formulation and SPH Method,” Nonlinear Dyn., 75(4), pp. 653–671. [CrossRef]
Liu, C. , Tian, Q. , and Hu, H. Y. , 2011, “ Dynamics of Large Scale Rigid-Flexible Multibody System Composed of Composite Laminated Plates,” Multibody Syst. Dyn., 26(3), pp. 283–305. [CrossRef]
Nachbagauer, K. , 2013, “ Development of Shear and Cross-Section Deformable Beam Finite Elements Applied to Large Deformation and Dynamics Problems,” Ph.D. dissertation, Johannes Kepler University, Linz, Austria. https://hal.archives-ouvertes.fr/hal-00855840/document
Orzechowski, G. , 2012, “ Analysis of Beam Elements of Circular Cross Section Using the Absolute Nodal Coordinate Formulation,” Arch. Mech. Eng., 59, pp. 283–296. http://yadda.icm.edu.pl/baztech/element/bwmeta1.element.baztech-ca63dc64-4389-4e6b-8091-6bad7462c8ef
Orzechowski, G. , and Frączek, J. , 2012, “ Integration of the Equations of Motion of Multibody Systems Using Absolute Nodal Coordinate Formulation,” Acta Mech. Autom., 6(2), pp. 75–83. https://www.researchgate.net/publication/268060752_Integration_of_the_equations_of_motion_of_multibody_systems_using_absolute_nodal_coordinate_formulation
Orzechowski, G. , and Frączek, J. , 2015, “ Nearly Incompressible Nonlinear Material Models in the Large Deformation Analysis of Beams Using ANCF,” Nonlinear Dyn., 82(1–2), pp. 451–464. [CrossRef]
Tian, Q. , Chen, L. P. , Zhang, Y. Q. , and Yang, J. Z. , 2009, “ An Efficient Hybrid Method for Multibody Dynamics Simulation Based on Absolute Nodal Coordinate Formulation,” ASME J. Comput. Nonlinear Dyn., 4(2), p. 021009. [CrossRef]
Tian, Q. , Sun, Y. L. , Liu, C. , Hu, H. Y. , and Paulo, F. , 2013, “ Elasto-Hydro-Dynamic Lubricated Cylindrical Joints for Rigid-Flexible Multibody Dynamics,” Comput. Struct., 114–115, pp. 106–120. [CrossRef]
Shabana, A. A. , 2012, Computational Continuum Mechanics, 2nd ed., Cambridge University Press, New York.
Kreyszig, E. , 1991, Differential Geometry, Dover Publications, Mineola, NY.
Simo, J. C. , 1985, “ A Finite Strain Beam Formulation. The Three Dimensional Dynamics Problem—Part I,” Comp. Meth. Appl. Mech. Eng., 49(1), pp. 55–70. [CrossRef]
Simo, J. C. , and Vu-Quoc, L. , 1986, “ A Three-Dimensional Finite-Strain Rod Model—Part II: Computational Aspects,” Comput. Methods Appl. Mech. Eng., 58(1), pp. 79–116. [CrossRef]
Goldstein, H. , 1950, Classical Mechanics, 1st ed., Addison Wesley, Boston, MA.
Roberson, R. E. , and Schwertassek, R. , 1988, Dynamics of Multibody Systems, 1st ed., Springer-Verlag, Berlin. [CrossRef]
Shabana, A. A. , 2010, Computational Dynamics, 3rd ed., Wiley, Chichester, UK. [CrossRef]
Shampine, L. , and Gordon, M. , 1975, Computer Solution of ODE: The Initial Value Problem, Freeman, San Francisco, CA.


Grahic Jump Location
Fig. 2

Shear angle θ at a middle point (–□– M=1 N⋅m, –○– M=10 N⋅m)

Grahic Jump Location
Fig. 3

Spatial derivative of shear angle θ at a middle point (–□– M=1 N⋅m, –○– M=10 N⋅m)

Grahic Jump Location
Fig. 4

Free falling flexible pendulum

Grahic Jump Location
Fig. 5

Vertical displacement at tip point using 6 ANCF finite elements (–□– method 1, –○– method 2, and –Δ– method 3)

Grahic Jump Location
Fig. 6

Magnitude of the ry vector at tip point using method 2with 6 ANCF finite elements (–□– k=0, –○– k=103, –Δ– k=104)

Grahic Jump Location
Fig. 7

Different curvature definitions at tip point using method 3 with 40 ANCF finite elements (–□– |rss|, –○– ∂θ/∂s)

Grahic Jump Location
Fig. 8

Magnitude of the rss vector at the middle point using 40 ANCF finite elements (–□– method 1, –○– method 2, and –Δ– method 3)

Grahic Jump Location
Fig. 9

Shear angle at middle point using 40 ANCF elements (–□– method 1, –○– method 2, and –Δ– method 3)

Grahic Jump Location
Fig. 10

Energy balance obtained using method 2 with 6 ANCF finite elements (–□– kinetic energy, –○– elastic energy, –Δ– potential energy, and — total energy)

Grahic Jump Location
Fig. 11

Energy balance obtained using method 3 with 6 ANCF finite elements (–□– kinetic energy, –○– elastic energy, –Δ– potential energy, and — total energy)

Grahic Jump Location
Fig. 12

Derivative of θ with respect to the arc length at tip point using 40 bilinear CRBF elements (–□– method 1 and –Δ– method 3)

Grahic Jump Location
Fig. 13

Shear angle at tip point using 24 bilinear CRBF elements (–□– method 1 and –Δ– method 3)




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