Research Papers

Singularities in Differential-Algebraic Boundary-Value Problems Governing the Excitation Response of Beam Structures

[+] Author and Article Information
Mehdi Saghafi

Department of Mechanical
Science and Engineering,
University of Illinois at Urbana-Champaign,
Urbana, IL 61801
e-mail: msagha2@illinois.edu

Harry Dankowicz

Department of Mechanical
Science and Engineering,
University of Illinois at Urbana-Champaign,
Urbana, IL 61801
e-mail: danko@illinois.edu

1Corresponding author.

2Elements of this paper have appeared in Saghafi & Dankowicz, “Nondegenerate Continuation Problems for the Excitation Response of Nonlinear Beam Structures,” Proceedings of ASME 2013 IDETC/CIE, Paper #13115, 2013.

Manuscript received September 10, 2013; final manuscript received March 14, 2014; published online October 14, 2014. Assoc. Editor: Arend L. Schwab.

J. Comput. Nonlinear Dynam 10(1), 011017 (Oct 14, 2014) (10 pages) Paper No: CND-13-1218; doi: 10.1115/1.4027208 History: Received September 10, 2013; Revised March 14, 2014

The objective of this paper2 is to identify and, where possible, resolve singularities that may arise in the discretization of spatiotemporal boundary-value problems governing the steady-state behavior of nonlinear beam structures. Of particular interest is the formulation of nondegenerate continuation problems of a geometrically-nonlinear model of a slender beam, subject to a uniform harmonic excitation, which may be analyzed numerically in order to explore the parameter-dependence of the excitation response. In the instances of degeneracy investigated here, the source is either found (i) directly in a differential-algebraic system of equations obtained from a finite-element-based spatial discretization of the governing partial differential boundary-value problem(s) together with constraints on the trial functions or (ii) in the further collocation-based discretization of the time-periodic boundary-value problem. It is shown that several candidate spatial finite-element discretizations of a mixed weak formulation of the governing boundary-value problem either result in (i) spatial group symmetries corresponding to equivariant vector fields and one-parameter families of periodic orbits along the group symmetry orbit or (ii) temporal group symmetries corresponding to ghost solutions and indeterminacy in a subset of the field variables. The paper demonstrates several methods for breaking the spatial equivariance, including projection onto a symmetry-reduced state space or the introduction of an artificial continuation parameter. Similarly, the temporal indeterminacy is resolved by an asymmetric discretization of the governing differential-algebraic equations. Finally, in the absence of theoretical bounds, computation is used to estimate convergence rates of the different discretization schemes, in the case of numerical calibration experiments performed on equilibrium and periodic responses for a linear beam, as well as for the full nonlinear models.

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


Arena, A., Formica, G., Lacarbonara, W., and Dankowicz, H., 2011, “Nonlinear Finite Element-Based Path Following of Periodic Solutions,” Proceedings of the ASME Design Engineering Technical Conference, Vol. 4, pp. 787–793.
Formica, G., Arena, A., Lacarbonara, W., and Dankowicz, H., 2013, “Coupling FEM With Parameter Continuation for Analysis of Bifurcations of Periodic Responses in Nonlinear Structures,” ASME J. Comput. Nonlinear Dyn., 8(2), p. 021013. [CrossRef]
Dankowicz, H., and Schilder, F., 2013, Recipes for Continuation, SIAM, Philadelphia, PA.
Lacarbonara, W., and Yabuno, H., 2006, “Refined Models of Elastic Beams Undergoing Large In-Plane Motions: Theory and Experiment,” Int. J. Solids Struct., 43(17), pp. 5066–5084. [CrossRef]
Hughes, T. J. R., 2000, The Finite Element Method: Linear Static and Dynamic Finite Element Analysis, Dover Publications, Mineola, NY.
Aguiar, R. M., Moleiro, F., and Soares, C. M. M., 2012, “Assessment of Mixed and Displacement-Based Models for Static Analysis of Composite Beams of Different Cross-Sections,” Compos. Struct., 94(2), pp. 601–616. [CrossRef]
Aguiar, R. M., 2010, “Numerical Assessment of Theories and Finite Element Models for Analysis of Laminated Composite Beams,” Master's thesis, Technical University of Lisbon, Lisbon, Portugal.
Krauskopf, B., Osinga, H. M., and Galán-Vioque, J., 2007, Numerical Continuation Methods for Dynamical Systems: Path Following and Boundary Value Problems, Springer, New York.
Bai, Y., 1992, “A Modified Lobatto Collocation for Linear Boundary Value Problems of Differential-Algebraic Equations,” Computing, 49(2), pp. 139–150. [CrossRef]
Wang, Z.-Q., and Guo, B.-Y., 2012, “Legendre-Gauss-Radau Collocation Method for Solving Initial Value Problems of First Order Ordinary Differential Equations,” J. Sci. Comput., 52(1), pp. 226–255. [CrossRef]
Siminos, E., and Cvitanović, P., 2011, “Continuous Symmetry Reduction and Return Maps for High-Dimensional Flows,” Physica D: Nonlinear Phenom., 240(2), pp. 187–198. [CrossRef]
Sheu, M.-G., 1978, “On Theories and Applications of Mixed Finite Element Methods for Linear Boundary-Value Problems,” Comput. Math. Appl., 4(4), pp. 333–347. [CrossRef]
Oden, J. T., and Reddy, J. N., 1976, “On Mixed Finite Element Approximations,” SIAM J. Numer. Anal., 13(3), pp. 393–404. [CrossRef]
Boffi, D., Brezzi, F., and Fortin, M., 2013, Mixed Finite Element Methods and Applications, Springer-Verlag, Berlin, Heidelberg.
Ergüven, M. E., and Gedikli, A., 2003, “A Mixed Finite Element Formulation for Timoshenko Beam on Winkler Foundation,” Comput. Mech., 31, pp. 229–237.
Kim, W., 2008, “Unconventional Finite Element Models for Nonlinear Analysis of Beams and Plates,” Ph.D. thesis, Texas A&M University, College Station, TX.


Grahic Jump Location
Fig. 1

Frequency-response curves for different values of the nonlinearity coefficient α (reprinted from Ref. [2])

Grahic Jump Location
Fig. 2

Schematic of the beam with clamped-clamped boundary condition

Grahic Jump Location
Fig. 3

Illustration of the effect of oddness and evenness of m (4 in left panels and 3 in right panels) and N (2 in upper panels and 3 in lower panels) on the satisfaction of periodicity condition of the ghost solution δχ(t)

Grahic Jump Location
Fig. 4

Equivariance results in a symmetry group orbit of periodic solutions shifted along v∧. A unique solution is obtained by appending a boundary condition at t = 0.

Grahic Jump Location
Fig. 5

Frequency-response curves obtained using the index-2 discretization (dashed) for M = 2, 4, and 6 and p = 4 compared to that obtained using the discretized scalar model (solid) for M = 8 and p = 4



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