Research Papers

Integration of Nonlinear Models of Flexible Body Deformation in Multibody System Dynamics

[+] Author and Article Information
Martin Schulze

e-mail: Martin.Schulze@simpack.de

Bernhard Burgermeister

Friedrichshafener Straße 1,
Gilching 82205, Germany

Andrey Tuganov

7-4270 Rue Saint-Dominique,
Montreal, QC H2W 2B1, Canada
e-mail: Andrey.Tuganov@gmail.com

Holger Lang

Chair of Applied Dynamics,
University of Erlangen-Nuremberg,
Konrad-Zuse-Str. 3/5,
Erlangen 91052, Germany
e-mail: Holger.Lang@ltd.uni-erlangen.de

Joachim Linn

Fraunhofer Institut für Techno-und
Wirtschaftsmathematik (ITWM),
Fraunhofer-Platz 1,
Kaiserslautern 67663, Germany
e-mail: Joachim.Linn@itwm.fraunhofer.de

Martin Arnold

Institute of Mathematics,
Martin Luther University Halle-Wittenberg,
Halle (Saale) 06099, Germany
e-mail: Martin.Arnold@mathematik.uni-halle.de

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received November 5, 2012; final manuscript received April 26, 2013; published online October 14, 2013. Assoc. Editor: Hiroyuki Sugiyama.

J. Comput. Nonlinear Dynam 9(1), 011012 (Oct 14, 2013) (10 pages) Paper No: CND-12-1191; doi: 10.1115/1.4025279 History: Received November 05, 2012; Revised April 26, 2013

Current challenges in industrial multibody system simulation are often beyond the classical range of application of existing industrial simulation tools. The present paper describes an extension of a recursive order-n multibody system (MBS) formulation to nonlinear models of flexible deformation that are of particular interest in the dynamical simulation of wind turbines. The floating frame of reference representation of flexible bodies is generalized to nonlinear structural models by a straightforward transformation of the equations of motion (EoM). The approach is discussed in detail for the integration of a recently developed discrete Cosserat rod model representing beamlike flexible structures into a general purpose MBS software package. For an efficient static and dynamic simulation, the solvers of the MBS software are adapted to the resulting class of MBS models that are characterized by a large number of degrees of freedom, stiffness, and high frequency components. As a practical example, the run-up of a simplified three-bladed wind turbine is studied where the dynamic deformations of the three blades are calculated by the Cosserat rod model.

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


Shabana, A. A., 1997, “Flexible Multibody Dynamics: Review of Past and Recent Developments,” Multibody Syst. Dyn., 1(2), pp. 189–222. [CrossRef]
Wallrapp, O., 1994, “Standardization of Flexible Body Modeling in Multibody System Codes, Part I: Definition of Standard Input Data,” Mech. Struct. Mach., 22(3), pp. 283–304. [CrossRef]
Rulka, W., 1998, “Effiziente Simulation der Dynamik Mmechatronischer Systeme für Industrielle Anwendungen,” Ph.D. thesis, Vienna University of Technology, Department of Mechanical Engineering.
Huber, F., 2006, “Application of SIMPACK for Simulating Valve Trains and Timing Mechanisms at Mercedes-Benz Automotive Engines,” SIMPACK User Meeting 2006, Available at: http://www.simpack.com/fileadmin/simpack/doc/usermeeting06/um06_dc-huber.pdf
Meyer, M., and Matthies, H. G., 2003, “Efficient Model Reduction in Non-Linear Dynamics Using the Karhunen-Loève Expansion and Dual-Weighted-Residual Methods,” Comput. Mech., 31(1–2), pp. 179–191. [CrossRef]
Dietz, S., Wallrapp, O., and Wiedemann, S., 2003, “Nodal vs. Modal Representation in Flexible Multibody System Dynamics,” Proceedings of Multibody Dynamics 2003, July 1–4 2003, IDMEC/IST, Lisbon, Portugal. Available at: http://www.simpack.com/uploads/media/fb_01_nodal_modal-SD_OW_SW_01.pdf
Géradin, M., and Cardona, A., 2001, Flexible Multibody Dynamics: A Finite Element Approach, John Wiley and Sons, Chichester.
Lang, H., Linn, J., and Arnold, M., 2011, “Multibody Dynamics Simulation of Geometrically Exact Cosserat Rods,” Multibody Syst. Dyn., 25(3), pp. 285–312. [CrossRef]
Linn, J., Lang, H., and Tuganov, A., 2012, “Geometrically Exact Cosserat Rods With Kelvin–Voigt Type Viscous Damping,” Proceedings of the IMSD2012, The 2nd Joint International Conference on Multibody System Dynamics, May 29–June 1, Stuttgart, Germany, P.Eberhard and P.Ziegler, Eds.
Featherstone, R., 2007, Rigid Body Dynamics Algorithms, Springer, New York.
Gruttmann, F., Sauer, R., and Wagner, W., 1998, “A Geometric Nonlinear Eccentric 3D-Beam Element With Arbitrary Cross-Sections,” Comput. Methods Appl. Mech. Eng., 160, pp. 383–400. [CrossRef]
Brandl, H., Johanni, R., and Otter, M., 1988, “A Very Efficient Algorithm for the Simulation of Robots and Similar Multibody Systems Without Inversion of the Mass Matrix,” Theory of Robots, P.Kopacek, I.Troch, and K.Desoyer, Eds.,. Pergamon, Oxford, pp. 95–100.
Otter, M., Elmqvist, H., and Cellier, F. E., 1996, “Modeling of Multibody Systems With the Object-Oriented Modeling Language DYMOLA,” Nonlinear Dyn., 9(1–2), pp. 91–112. [CrossRef]
Schiehlen, W., 1997, “Multibody System Dynamics: Roots and Perspectives,” Multibody Syst. Dyn., 1(2), pp. 149–188. [CrossRef]
Eich-Soellner, E., and Führer, C., 1998, Numerical Methods in Multibody Dynamics, Teubner, Stuttgart, Germany.
Arnold, M., Burgermeister, B., Führer, C., Hippmann, G., and Rill, G., 2011, “Numerical Methods in Vehicle System Dynamics: State of the Art and Current Developments,” Vehicle Syst. Dyn., 49(7), pp. 1159–1207. [CrossRef]
Brenan, K. E., Campbell, S. L., and Petzold, L. R., 1996, Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations, 2nd ed., SIAM, Philadelphia, PA.
Hairer, E., and Wanner, G., 1996, Solving Ordinary Differential Equations. II. Stiff and Differential-Algebraic Problems, 2nd ed., Springer, Berlin.
Negrut, D., and Ortiz, J., 2006, “A Practical Approach for the Linearization of the Constrained Multibody Dynamics Equations,” ASME J. Comput. Nonlinear Dyn., 1(3), pp. 230–239. [CrossRef]
Coleman, T. F., Garbow, B. S., and Moré, J. J., 1984, “Algorithm 618: FORTRAN Subroutines for Estimating Sparse Jacobian Matrices,” ACM Trans. Math Software, 10(3), pp. 346–347. [CrossRef]
Arnold, M., Fuchs, A., and Führer, C., 2006, “Efficient Corrector Iteration for DAE Time Integration in Multibody Dynamics,” Comput. Methods Appl. Mech. Eng., 195, pp. 6958–6973. [CrossRef]
Tuganov, A., 2009, “Elaborated Ball Screw Simulation in SIMPACK,” Proceedings of Multibody Dynamics 2009, K.Arczewski, J.Fraczek, and M.Wojtyra, Eds., ECCOMAS Thematic Conference, Warsaw, Poland.
Davis, T. A., 2004, “Algorithm 832: UMFPACK v4.3—An Unsymmetric-Pattern Multifrontal Method,” ACM Trans. Math Software, 30(2), pp. 196–199. [CrossRef]
Brown, P. N., Hindmarsh, A. C., and Petzold, L. R., 1994, “Using Krylov Methods in the Solution of Large-Scale Differential-Algebraic Systems,” SIAM J. Sci. Comput., 15(6), pp. 1467–1488. [CrossRef]
Jonkman, J., Butterfield, S., Musial, W., and Scott, G., 2008, “Definition of a 5-MW Reference Wind Turbine for Offshore System Development,” TP 500-38060, National Renewable Energy Laboratory.
Gasch, R., and Twele, J., 2011, Windkraftanlagen, Grundlagen, Entwurf, Planung und Betrieb, Vieweg and Teubner, Wiesbaden.
David, and Laino, J., 2013, NWTC Design Codes. http://wind.nrel.gov/designcodes/simulators/aerodyn/
Moriarty, P. J., and Hansen, A. C., 2005, AeroDyn Theory Manual, National Renewable Energy Laboratory, CO.
Leishman, J. G., 2000, Principles of Helicopter Aerodynamics, Cambridge University Press, Cambridge, UK, pp. 78–127.
Bauchau, O. A., and Hong, C. H., 1988, “Nonlinear Composite. Beam Theory,” ASME J. Appl. Mech., 55(1), pp. 156–163. [CrossRef]
Hodges, D. H., 2006, Nonlinear Composite Beam Theory, The American Institute of Aeronautics and Astronautics Reston, VA.
Hodges, D. H., and Yu, W., 2007, “A Rigorous, Engineer-Friendly Approach for Modelling Realistic, Composite Rotor Blades,” Wind Energy, 10, pp. 179–193. [CrossRef]
Yu, W., Volovoi, V. V., Hodges, D. H., and Hong, X., 2002, “Validation of the Variational Asymptotic Beam Sectional Analysis (VABS),” AIAA J., 40, pp. 2105–2113. [CrossRef]


Grahic Jump Location
Fig. 1

Interaction of SIMPACK and the Cosserat rod model

Grahic Jump Location
Fig. 2

Difference approximation of a band structured Jacobian matrix with the minimum number of rhs evaluations, N = 6

Grahic Jump Location
Fig. 3

Time evolution of a swinging rubber beam

Grahic Jump Location
Fig. 4

Wind turbine rotor blade being deformed by a constant wind field

Grahic Jump Location
Fig. 5

Static deformation of the NREL 5 MW wind turbine rotor blade in the constant force distribution of a wind field: Comparison of the Abaqus B31, the Cosserat rod, and the segmented SIMBEAM models

Grahic Jump Location
Fig. 6

Comparison of the first and the 12th eigenmodes of the NREL 5 MW wind turbine rotor blade obtained with the Abaqus B31 and with the Cosserat rod models, respectively (normalized units)

Grahic Jump Location
Fig. 7

Simplified model of a three-bladed wind turbine

Grahic Jump Location
Fig. 8

Steady state of the simplified wind turbine. The left diagram shows the blade deformation perpendicular to the movement plane, the right diagram the main shaft velocity for one rotation.




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