Projects per year
Abstract
Rotational motion is a frequently occurring aspect in many engineering applications such as automobiles, rotating machinery or wind turbine rotors. The present thesis is organized in four parts  all concerning development of efficient computational methods for modeling the dynamic behavior of rigid bodies and flexible beam structures with emphasis on the rotational motion.
The first part deals with motion in a rotating frame of reference. A novel approach where the equations of motion are formulated in a hybrid statespace in terms of local displacements and global velocities is presented. It is shown that particular simplifications are obtained when the same interpolation functions are used for both statespace variables, whereby all inertia effects can be represented via the classic constant mass matrix. The hybrid statespace constitutes the basis for developing a conservative time integration scheme and an associated algorithmic dissipation, that affects only the local motion.
In the second part a conservative integration scheme for rigid body motion in a global frame of reference is presented. A fully algebraic representation of the rotational motion is obtained by using either four quaternion parameters or nine convected base vector components. In both cases, the equations of motion are obtained via Hamilton’s equations by including the kinematic constraints associated with the redundant rotation description by means of Lagrange multipliers. A special feature of the formulation is that these can be replaced by a projection operator applied to the external potential gradient and possible external constraint gradients. The third part presents a novel twonode beam element capable of undergoing arbitrary large displacements and finite rotations. The element is expressed explicitly in terms of the global components of the position vectors and associated convected base vectors for the element nodes. The kinematics is expressed in a homogeneous quadratic form and the constitutive stiffness is derived from complementary energy of a set of equilibrium modes, each representing a state of constant internal force or moment. This approach avoids local interpolation of kinematic variables, which makes the formulation inherently lockingfree and frameinvariant. In the final part of the thesis a multilevel optimization procedure for wind turbine blades is presented. Computational procedures similar to the ones presented in this thesis constitute the basis for a preliminary optimization using a multibody beam model of the full wind turbine. The novel aspects considered here consist of automatic generation of a 3D FEmodel of a single blade for detailed stress and buckling analysis which cannot be performed on the beam model, and utilization of the results for updating the beam model for a subsequent iteration step. It is demonstrated that convergence between the two models is obtained in very few iterations.
The first part deals with motion in a rotating frame of reference. A novel approach where the equations of motion are formulated in a hybrid statespace in terms of local displacements and global velocities is presented. It is shown that particular simplifications are obtained when the same interpolation functions are used for both statespace variables, whereby all inertia effects can be represented via the classic constant mass matrix. The hybrid statespace constitutes the basis for developing a conservative time integration scheme and an associated algorithmic dissipation, that affects only the local motion.
In the second part a conservative integration scheme for rigid body motion in a global frame of reference is presented. A fully algebraic representation of the rotational motion is obtained by using either four quaternion parameters or nine convected base vector components. In both cases, the equations of motion are obtained via Hamilton’s equations by including the kinematic constraints associated with the redundant rotation description by means of Lagrange multipliers. A special feature of the formulation is that these can be replaced by a projection operator applied to the external potential gradient and possible external constraint gradients. The third part presents a novel twonode beam element capable of undergoing arbitrary large displacements and finite rotations. The element is expressed explicitly in terms of the global components of the position vectors and associated convected base vectors for the element nodes. The kinematics is expressed in a homogeneous quadratic form and the constitutive stiffness is derived from complementary energy of a set of equilibrium modes, each representing a state of constant internal force or moment. This approach avoids local interpolation of kinematic variables, which makes the formulation inherently lockingfree and frameinvariant. In the final part of the thesis a multilevel optimization procedure for wind turbine blades is presented. Computational procedures similar to the ones presented in this thesis constitute the basis for a preliminary optimization using a multibody beam model of the full wind turbine. The novel aspects considered here consist of automatic generation of a 3D FEmodel of a single blade for detailed stress and buckling analysis which cannot be performed on the beam model, and utilization of the results for updating the beam model for a subsequent iteration step. It is demonstrated that convergence between the two models is obtained in very few iterations.
Original language  English 

Publisher  DTU Mechanical Engineering 

Number of pages  200 
ISBN (Print)  9788774753735 
Publication status  Published  2013 
Series  DCAMM Special Report 

Number  S156 
ISSN  09031685 
Keywords
 Finite rotations
 Conservative time integration
 Structural dynamics
 Constrained mechanics
 Beam elements
 Quaternions
 Cosine directors
 Wind turbines
Fingerprint Dive into the research topics of 'Dynamics of Rigid Bodies and Flexible Beam Structures'. Together they form a unique fingerprint.
Projects
 1 Finished

Modeling of Rotating Structures
Nielsen, M. B., Høgsberg, J. B., Cardona, A., Romero, I. & Krenk, S.
Technical University of Denmark
01/09/2010 → 28/04/2014
Project: PhD