Research Papers

Projective Constraint Stabilization for a Power Series Forward Dynamics Solver

[+] Author and Article Information
Paul Milenkovic

Department of Electrical and Computer Engineering,
University of Wisconsin—Madison,
1415 Engineering Drive,
Madison, WI 53706
e-mail: phmilenk@wisc.edu

The top panel of Fig. 4 from the earlier paper [5] is inverted in relation to Table 1 in that paper. The corrected orientation is given in Fig. 4 of the current paper.

Contributed by the Dynamic Systems Division of ASME for publication in the JOURNAL OF DYNAMIC SYSTEMS, MEASUREMENT, AND CONTROL. Manuscript received March 8, 2012; final manuscript received November 30, 2012; published online February 21, 2013. Assoc. Editor: Nariman Sepehri.

J. Dyn. Sys., Meas., Control 135(3), 031004 (Feb 21, 2013) (9 pages) Paper No: DS-12-1079; doi: 10.1115/1.4023212 History: Received March 08, 2012; Revised November 30, 2012

A power series expression for the forward dynamics of a closed kinematic chain provides an explicit time-step update of the system state. The resulting numerical differential equation solver applies kinematic constraints to the power series terms for acceleration and higher derivatives of motion. Integrating acceleration determines velocity and position time histories that approximate the constraints to a high degree of precision when using a high order of the expansion. When high precision is not required, a lower order achieves shorter computation times, but that condition results in violation of the constraints in the absence of any correction. Projecting the velocities and positions onto the constraint manifold after each time step produces step changes. This paper determines which choices of linear subspace for this projection give step changes that are equal to the residues of truncating the power series solution for the kinematic portion of the problem. The limit of that power series gives position and velocity time histories that approximate the dynamics while giving an exact kinematic solution. Thus projection onto the constraints in this procedure determines sample values of an underlying solution for the motion trajectories, where that underlying solution is continuous in both velocity and position and also satisfies the kinematic constraints at all times. This property is confirmed by numerical simulation of a Clemens constant-velocity coupling.

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



Grahic Jump Location
Fig. 1

Clemens CV coupling at initial posture: direction vectors ω1,…,ω9 are placed along the joint axis lines

Grahic Jump Location
Fig. 2

Planar projection of the screw motion of a point following an arc of radius R and angle θ

Grahic Jump Location
Fig. 3

Dynamic response of Clemens CV coupling input shaft rotation angle

Grahic Jump Location
Fig. 4

Dynamic response of Clemens CV coupling output shaft pitch and yaw deflection angles

Grahic Jump Location
Fig. 5

Projective correction to the order 2 solution of the input shaft rotation angle

Grahic Jump Location
Fig. 6

Projective correction to the order 2 solution of the input shaft rotation rate

Grahic Jump Location
Fig. 7

Projective correction to the order 4 solution of the input shaft rotation rate

Grahic Jump Location
Fig. 8

Relative computation times for the order 8 forward-dynamics solution




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