Mathematical Problems in Engineering

Volume 2012 (2012), Article ID 949834, 22 pages

http://dx.doi.org/10.1155/2012/949834

## Parallel-Distributed Model Deformation in the Fingertips for Stable Grasping and Object Manipulation

^{1}Facultad de Ingeniería y Ciencias Aplicadas, Universidad de los Andes, Av. San Carlos de Apoquindo 2200, Las Condes, Santiago, Chile^{2}Departamento de Ingeniería Eléctrica, Universidad de Chile, Av. Tupper 2007, Santiago, Chile

Received 27 April 2012; Revised 13 July 2012; Accepted 25 July 2012

Academic Editor: J. Rodellar

Copyright © 2012 R. García-Rodríguez and G. Díaz-Rodríguez. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The study on the human grip has inspired to the robotics over the past decades, which has resulted in performance improvements of robotic hands. However, current robotic hands do not have the enough dexterity to execute complex tasks. Recognizing this fact, the soft fingertips with hemispherical shape and deformation models have renewed attention of roboticists. A high-friction contact to prevent slipping and the rolling contribution between the object and fingers are some characteristics of the soft fingertips which are useful to improve the grasping stability. In this paper, the parallel distributed deformation model is used to present the dynamical model of the soft tip fingers with *n*-degrees of freedom. Based on the joint angular positions of the fingers, a control scheme that fuses a stable grasping and the object manipulation into a unique control signal is proposed. The force-closure conditions are defined to guarantee a stable grasping and the boundedness of the closed-loop signals is proved. Furthermore, the convergence of the contact force to its desired value is guaranteed, without any information about the radius of the fingertip. Simulation results are provided to visualize the stable grasping and the object manipulation, avoiding the gravity effect.

#### 1. Introduction

From physiological point of view, the human hands are considered as a powerful tool whereby the human brain interacts with the world, that is, how it perceives and acts with the environment [1]. In order to increase the dexterity in robotic hands, some intelligent human-like functions have been imitated.

In general, the dexterous manipulation in robotics, to emulate pinching motions, have been formulated in terms of the object, that is, the forces/torques exerted on it to produce the desired movements and how they behave [2]. The grasping and the object manipulation are based on the assumption that the contact between the object and fingers is frictionless, so the finger can only exert a force along the common normal axis at the contact point [3]. Then, to grasp the object without slipping the standard friction cone is used, generating a complex motion control since their evolution is governed by the laws of Coulomb friction, which is nonlinear and imposes constraints on the system. On the other hand, some authors have considered fingers with a very sharp curvature assuming that the contact point between the fingers and object does not change significantly. Although, some manipulation tasks are executed by robotic fingers successfully, this assumption is not valid for several manipulation tasks because the rolling between object and fingers is essential in the human manipulation tasks. Moreover*,* some attempts to determine the best grasp configuration and manipulation tasks are presented in 2D and 3D [4–11]. Unfortunately, a lot of them require an exact knowledge of the system parameters and the object localization [12–15]. Some authors, to reproduce more characteristics of the human fingertips, have considered the use of deformation models with hemispherical soft tips. Many hemispherical soft tip fingers have been designed and constructed to execute several manipulation tasks [16–19] where a high-contact friction to prevent slipping and the rolling of the finger tip on the object surface are some characteristics [3, 20, 21]. Nevertheless, in these approaches the contribution of the fingertip deformation on the manipulation tasks in a dynamic sense is not evident, considering that the rolling constraint is defined in kinematic or semidynamics sense [21–24].

On the basis that human hands have fingers with soft tips, in this paper the grasping and the object manipulation is presented, using a pair of robotic fingers with *n*-degree of freedoms and deformable tips. The parallel deformation model is based on a virtual spring with infinitesimal section, where the normal and tangential deformations are taking into account [25]. So, a tangential movement of the object without slipping and a dependency of the relative orientation between the object and the finger are considered. Inclusion of the normal and tangential deformations in the deformation model, contribute to reproduce some intrinsic characteristics of the deformable material and a better grip. The key of our approach is to introduce the parallel deformation model to grasp an object using a pair of robotic fingers with *n*-degrees of freedom. Moreover, the grasping controller guaranteed that the contact force converges to the desired value, avoiding a direct dependence with the radius of the tip [26, 27]. An approximation of the object angle based on the joint angular position of the fingers is used to control orientation of the object. Finally, a control signal for translation of the object is defined. To carry out the grasping and the object manipulation, the superposition principle is used [28, 29], which allows us to separate a complex task into a set of basic tasks, where each task has a unique stationary point which represents the desired action [27]. Boundedness of all closed-loop signals is proved, while the asymptotic stability is guaranteed using the stability on the manifold [27, 28, 30]. It is important to notice that forces-closure conditions, to grasp firmly an object, are satisfied dynamically during manipulation task execution, rather than a static equilibrium. The proposed approach is validated by numerical simulations through a pair of robotic fingers with soft tips in the horizontal plane.

This paper is organized as follows. Section 2 presents dynamical equations of the fingers-object system. The blind control law is proposed in Section 3. Simulation results to confirm the validity of our approach are presented in Section 4. Finally, the conclusions are presented in Section 5.

#### 2. Dynamical Equations

Consider a pair of soft tip fingers, with three degree of freedom each one, grasping an object in the horizontal plane, as shown in Figure 1. In the fingers-object system is the origin for the left finger and it is considered as the reference frame, is the origin for the right finger, and is the distance between the origins of each finger. In addition, is the joint angular positions of the finger , is the radius of the soft tip finger , is the length of the object, is the length of the link for the finger , is the center position of the deformable fingertips with , is the center of mass of the object, and is the orientation of the object. Unlike deformation model proposed by [22, 28], where the force applied to the object produce a distribute pressure and it is parametrized as a normal force with respect to the object surface. In this paper we use the deformation model proposed in [25], where a virtual spring inside the soft tip finger allows us to known the normal and tangential movements in the deformable material. In such a way, the elastic energy induced by the soft fingertip is given as where is the maximum radial deformation, is the object relative orientation angle, is the contact tangential displacement of the object, and is the Young's modulus of the finger tip material, as shown in Figure 2. Thus, the total potential energy of the deformable fingertips is expressed as The constraint between the radial deformation of the fingertip and the object is given as which guarantees that exists a distance that limits the grasping on the object in normal direction. A particular case of (2.3) is considered when which represents the normal constraint for a rigid fingertip as reported in [30]. Accordingly, the normal constrained for the fingers-object system is defined as where is the Lagrange multiplier and represents the contact force .

On the other hand, taking into account the curvature effects of the fingertip*, *the rolling of the object on the soft tip finger will be defined as movement of the contact area, where the relative velocity of the contact area between the finger tips and the object is zero. Assuming that the normal deformation is smaller than the radius of the tip [25], the angular displacement of the contact area and the projection displacement of the center of mass of the object is defined as
where
and of the same size as the vector , for . Moreover, according to the deformation model, a tangential displacement on the deformable fingertip arises when the object is rolling on the fingertip. Then, the rolling constraint between finger tip and the object surface is given as
where is the integration constant with respect to the initial conditions of contact. To avoid initial conditions in (2.7) a velocity constraint is defined as [25]
The Lagrangian of the system with holonomic constraint is described by
where and are defined in (2.2) and (2.4), respectively, and is the kinetic energy of the system defined as
with is the inertia matrix of the finger , , and , are the normal and tangential mass deformations, respectively. Applying the Lagrangian variational principle, incorporating the rolling velocity constraint, the equations of motion are expressed for each component of the vector as follows:
where is the vector of generalized coordinates and is the Lagrange multiplier which represents the tangential force exerted for the finger on the object surface. Note that treatment of the velocity constraint should not be done in the Lagrangian (2.9), but rather in the equations of motion [31].

Thus, the equations of motion for the fingers are given as where , is the Jacobian of the point with respect to the joint variables , , represent the matrix of Coriolis and centripetal forces, and stands for the torque input. It is important to notice that the moments induced by the tangential and normal forces contribute to grasp and manipulate an object more securely. In addition, the terms of deformation model give us information about the behavior of the deformable material. Furthermore, the movement equations of the object are given as The last two terms of (2.15) represent the contribution of the deformable fingertips to assure a stable grasping through induced forces/moments. Finally, dynamical equations related to the normal and tangential movements on the fingertips are defined as

Summing the products between with (2.12), with (2.13), with (2.14), with (2.15), with (2.16), and with (2.17) yields where corresponds to the total energy of the system.

#### 3. Controller Design

##### 3.1. Immobilization on the Object

As first step before to execute manipulation tasks, a stable grasping must be established. To grasp stably an object, the force-closure is used to guarantee that the object should be held securely by the fingers. This mean that maintaining the contact between the fingers and the object, that is, and for any , the forces and torques applied on the object should immobilize it. Let the forces and torques applied on the object be defined as Then, if we choose that the first two equations are equal to zero, while the third equation is given as Thus, a force-closure can be established, if the forces acting on the object are defined as, where is the desired normal force. Once fingers grasp an object and hold it securely, we are in conditions to execute manipulation tasks on the object as orientation and move it at coordinates.

Inspired that humans can execute some manipulation tasks without any object information, in this paper a stable grasping and the object manipulation based only on the center position of the soft tip finger are presented, so that the orientation and the object parameters are avoided. Using the superposition principle*,* the joint torque applied to each finger can be defined as
where is a diagonal symmetric positive definite matrix, , is the desired contact force, , is the desired angle of rotation, and
Notice that refers to the stable grasp exerted on the object, while indicates the orientation control of the object. In the former case, a stable grasp is achieved minimizing the distance between the centers of the fingertips through the normal and tangential forces, which guarantee the control objectives , . In the latter case, to avoid the measurement of , an approximation of is proposed by the trigonometric tangent function . This approximation has the feature that the convergence of to is guaranteed once is satisfied. This implies that the conditions of the stable grasping , , and should be satisfied.

To start the analysis and to ensure a stable grasping, it is considered that the joint torque is defined as . Substituting in (2.12), the closed-loop system equations are defined as where , , and .

Expressing the closed-loop system equations (3.7) in a vector-matrix equation we have that where with , for .

Taking the inner product between and (3.8) yields where is called artificial potential energy. Then, the system satisfys the passivity condition in closed loop.

Although is negative definite along the solution trajectories of fingers-object system, cannot play a role of a Lyapunov function for the closed-loop system. Because it is neither define positive in the 26-dimensional state space nor in 18-dimensional constrained manifold defined by

In order to define a constrained manifold where the trajectories of the system guarantees a stable grasping; a first step is to show the boundedness of solutions in the closed loop. Considering that and we have that and , respectively. In fact, if the matrix from previous expressions we have that . Then,

Now, if we multiply (3.8) by we obtain that where from (3.12). Taking into account that , , , , and are bounded due to ; we have that is bounded under assumption that the matrix in nondegenerate. Then, from (3.8) we have that is uniformly bounded which implies that is uniformly continuous. Now, given that from (3.10) by Barbalat lemma [32] we have that as which implies that as . Due to and are uniformly continuous, we have that is uniformly continuous too. This in turn implies that at and (3.8) will be defined as Assuming that the matrix is nondegenerate, a point that minimizes the right side of the (3.14) on exists. The unique critical point that minimizes (3.14) can be defined as This means that the closed-loop trajectories converge to critical point in the equilibrium point manifold (EP) defined as Hence, the total potential energy is minimizing and the object is held securely. It is important to notice that initial conditions of the system and the length of the object are very important parameters to guarantee the convergence on . Furthermore, the matrix is non-degenerate if and are full rank matrices, that is, the matrix is degenerate if or for . These joint values represent a special configuration of the fingers which can be excluded as a possible configuration in the manipulation tasks. At the same time, it is possible to show that the matrix is nondegenerate if the and are full rank matrices.

Now, we are in conditions to define the following.(i)Let the neighborhood singularity-free be with a radius around the critical point defined as
where , . That is, singular configurations of the soft finger tips are avoided inside the neighborhood provide that is chosen adequately.(ii)Let the neighborhood around the critical point be where the closed-loop trajectories remain with , that is,
(iii)Definition** **(see [26]) If for any given there exists and another constant being less than and independent of such that the solution tracking trajectories starting from any initial condition lying on remains on and converges asymptotically to the set as , then it is called that the state is stable.

Finally, we have the following result.

Theorem 3.1. *Considering that the desired reference state and initial state lying on . The trajectories in the closed-loop system remains on and converges asymptotically to the set as under assumption that the matrix is nondegenerate in a neighborhood . Thus, is assured that , , as .*

*Remark 3.2. *Once the forces applied to the object have been compensated to hold it stably, the object can be rotated to the desired angle using the superposition principle, that is, the control law now is defined as
Premultiplying by we have that
where . This means that there exists a constant such that
where .

Now, taking the inner product between and the closed-loop system equations with defined in (3.19) we have that
where and . As in stable grasping, the closed-loop trajectories converge to critical point on a constrained manifold where , , and as .

*Remark 3.3. *Now, when the stable grasping and object orientation tasks has been established, the objective will be to move the object to desired coordinates by the following control law
where , , is the Cartesian coordinates and is the estimated position of the object which is defined as an average distance between center positions of the deformable fingertips, that is,
Using the superposition principle, the control law for the stable grasping and object manipulation is defined as
As in the previous case, if we multiplying by we have that
where is the shifting energy to move the object in -coordinates. Consequently, exists a constant such that
where .

Now, taking the inner product between and the closed-loop system equations with defined in (3.25) we have where and . Then, the passivity condition is satisfy in the closed-loop. Hence, the closed-loop system trajectories converge to critical point on a constrained manifold where , , , , and as .

#### 4. Simulation Results

In order to demonstrate usefulness of our scheme for stable grasping and object orientation, numerical simulations were carried out on a pair of deformable fingertips in the horizontal plane, see Figure 1. The simulations were implemented on stiff numerical solver on Matlab R2007b, under 1 ms sampling time. Additionally, to approximate the holonomic constraints was used the Constrained Stabilization Method (CSM) [33].

The physical parameters of the fingers and object are shown in Table 1 where , and are the length, mass and moment of inertia of the link for the finger , and , , are the mass, moment of inertia and the length of the object, respectively (Table 2). Moreover, is the distance between fingers, is the Young's modulus of the fingertips and is the radius of the hemispherical finger tip for .

The simulation study is divided in three steps. As first step before any manipulation task is necessary to guarantee that the stable grasp is achieved through the control law defined as . The initial conditions used in the simulations are , , , , and which establish the following conditions: , . For reference, these initial conditions are called normal initial conditions .

In Figure 3 we observe that the forces on the object and converge rapidly to the desired values using the normal initial conditions defined previously. At the same time, the fast convergence of to zero is shown in Figure 3. On the other hand, the Figures 4, 5, and 6 show the performance of , , and for different values of damping gains . Notice that the damping gains are related with the convergence velocity towards equilibrium point that minimizes the potential energy of the system. However, the system takes longer to converge and several oscillations arise when the damping gains are small. The control parameters used in this step are and .

Finally, the Figures 7, 8, and 9 show the convergence of and to zero under more extreme initial conditions. The initial conditions used in these simulations are and which are defined as , , , and , , , , respectively.

Notice that the convergence to zero of and present a transient responses for different initial conditions, but still converge to the desired values in few seconds. Additionally, the control objectives converge in different times. Specially, the convergence of takes more time for establishing a stable grasping when the initial conditions are more extreme, . Thus, it is possible to describe the stable grasping in two phases. In the first phase, the stabilization of the normal and tangential forces on object is performed while the stabilization of the rotational moments to stop the angular motion of the object is carried out in a second phase.

Once the all forces applied to the object have been compensated, it is possible to execute any manipulation task. In this case the rotation of the object, to a desired angle , will be realized. Using the superposition principle the control law is given as . The convergence to zero of , , , and using the initial conditions are shown in Figures 10 and 11. It is important to notice the special role of as a parameter to increase the dexterity. In this case, the convergence of is closely associated to object orientation through . The control parameters using in this second step are , , , and .

As the final step of this study, a stable grasping and object manipulation which include the orientation and translation of the object to a desired reference is presented. Using a superposition principle, the control law for this step is defined as .

The Figures 12, 13, and 14 show the convergence of , , , and to zero using initial conditions. The control parameters used in this step are , , , , , and .

#### 5. Conclusions

A scheme to grasp and manipulate an object using soft tip fingers is presented. In order to include more characteristics of the human fingertips a parallel deformation model is used. The control law proposed ensures stability on a constrained manifold. Furthermore, the control law avoids information of the radius of the tips and the convergence to the desired force value is guaranteed.

Numerical simulations, on a pair of deformable fingertips in horizontal plane, allow us to visualize the convergence of the closed-loop trajectories to the desired point. In addition, the special role of in the manipulation task and the effects of the superposition principle that have been observed.

#### Acknowledgments

This work was partially supported by CONICYT, Departamento de Relaciones Internacionales, Programa de Cooperación Científica Internacional, CONICYT/CONACYT 2011-380, and the Universidad de los Andes, Chile, FAI project ICI-002-11.

#### References

- J. R. Flanagan and R. S. Johansson,
*Hand Movements, Encyclopedia of the Human Brain*, vol. 2, Academic Press, San Diego, Calif, USA, 2002. - A. M. Okamura, N. Smaby, and M. R. Cutkosky, “Overview of dexterous manipulation,” in
*Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '00)*, pp. 255–262, April 2000. View at Scopus - J. K. Salisbury and B. Roth, “Kinematic and force analysis of articulated mechanical hands,”
*Journal of Mechanisms, Transmissions, and Automation in Design*, vol. 105, no. 1, pp. 35–41, 1983. View at Google Scholar · View at Scopus - E. A. Al-Gallaf, “Multi-fingered robot hand optimal task force distribution: neural inverse kinematics approach,”
*Robotics and Autonomous Systems*, vol. 54, no. 1, pp. 34–51, 2006. View at Publisher · View at Google Scholar · View at Scopus - A. Bicchi, “Hands for dexterous manipulation and robust grasping: a difficult road toward simplicity,”
*IEEE Transactions on Robotics and Automation*, vol. 16, no. 6, pp. 652–662, 2000. View at Publisher · View at Google Scholar · View at Scopus - W. Gueaieb, S. Al-Sharhan, and M. Bolic, “Robust computationally efficient control of cooperative closed-chain manipulators with uncertain dynamics,”
*Automatica*, vol. 43, no. 5, pp. 842–851, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J.-W. Li, H. Liu, and H.-G. Cai, “On computing three-finger force-closure grasps of 2-D and 3-D objects,”
*IEEE Transactions on Robotics and Automation*, vol. 19, no. 1, pp. 155–161, 2003. View at Publisher · View at Google Scholar · View at Scopus - H. Kawasaki, R. Bin Ramli, and S. Ueki, “Decentralized adaptive coordinated control of multiple robot arms for constrained tasks,”
*Journal of Robotics and Mechatronics*, vol. 18, no. 5, pp. 580–588, 2006. View at Google Scholar - T. Naniwa and K. Wada, “Experimental study of coordinated control of multifingered hands with the kinetostatic filtering method,”
*Advanced Robotics*, vol. 19, no. 2, pp. 191–206, 2005. View at Publisher · View at Google Scholar · View at Scopus - E. Rimon, “A curvature-based bound on the number of frictionless fingers required to immobilize three-dimensional objects,”
*IEEE Transactions on Robotics and Automation*, vol. 17, no. 5, pp. 679–697, 2001. View at Publisher · View at Google Scholar · View at Scopus - X.-Z. Zheng, R. Nakashima, and T. Yoshikawa, “On dynamic control of finger sliding and object motion in manipulation with multifingered hands,”
*IEEE Transactions on Robotics and Automation*, vol. 16, no. 5, pp. 469–481, 2000. View at Publisher · View at Google Scholar · View at Scopus - A. B. A. Cole, J. E. Hauser, and S. S. Sastry, “Kinematics and control of multifingered hands with rolling contact,”
*Institute of Electrical and Electronics Engineers*, vol. 34, no. 4, pp. 398–404, 1989. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Nakashima, K. Nagase, and Y. Hayakawa, “Control of contact points of two- fingered robot hand with manipulating an object,” in
*Proceedings of the ICASE/SICE Joint Workshop*, pp. 120–125, 2004. - A. I. Tuneski and G. Rafajlovski, “Adaptive control of multiple robots manipulation on dynamical environment,” in
*Proceedings of the 7th IEEE International Conference on Emerging Technologies and Factory Automation (ETFA '99)*, vol. 1, pp. 623–632, October 1999. View at Scopus - M. Zribi, J. Chen, and M. S. Mahmoud, “Coordination and control of multi-fingered robot hands with rolling and sliding contacts,”
*Journal of Intelligent and Robotic Systems*, vol. 24, no. 2, pp. 125–149, 1999. View at Google Scholar · View at Scopus - C. Borst, M. Fischer, and G. Hirzinger, “Grasp planning: how to choose a suitable task wrench space,” in
*Proceedings IEEE International Conference on Robotics and Automation*, pp. 319–325, May 2004. View at Scopus - I. Kao and F. Yang, “Stiffness and contact mechanics for soft fingers in grasping and manipulation,”
*IEEE Transactions on Robotics and Automation*, vol. 20, no. 1, pp. 132–135, 2004. View at Publisher · View at Google Scholar · View at Scopus - S. K. Song, J. B. Park, and Y. H. Choi, “Grasping control of 3-joint dual finger robot: lyapunov stability approach,” in
*Proceedings of the American Control Conference (ACC '09)*, pp. 2879–2884, St. Louis, Mo, USA, June 2009. View at Publisher · View at Google Scholar · View at Scopus - Y. Xue and I. Kao, “Dextrous sliding manipulating using soft fingertips,” in
*Proceedings of the IEEE International Conference on Robotics and Automation*, vol. 4, pp. 3397–3402, May 1994. View at Scopus - P. Akella and M. Cutkosky, “Manipulating with soft fingers,” in
*Proceedings of the IEEE International Conference on Robotics and Automation*, pp. 767–769, 1989. - K. B. Shimoga and A. A. Goldenberg, “Soft materials for robotic fingers,” in
*Proceedings of the IEEE International Conference on Robotics and Automation*, pp. 1300–1305, May 1992. View at Scopus - I. Kao and F. Yang, “Stiffness and contact mechanics for soft fingers in grasping andmanipulation,”
*IEEE Transactions on Robotics and Automation*, vol. 20, no. 1, pp. 132–135, 2004. View at Publisher · View at Google Scholar · View at Scopus - D. J. Montana, “Kinematics of contact and grasp,”
*International Journal of Robotics Research*, vol. 7, no. 3, pp. 17–32, 1988. View at Google Scholar · View at Scopus - Y. Yokokohji, M. Sakamoto, and T. Yoshikawa, “Vision-aided object manipulation by a multifingered hand with soft fingertips,” in
*Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '99)*, vol. 4, pp. 3201–3208, May 1999. View at Scopus - T. Inoue and S. Hirai,
*Mechanics and Control of Soft-Fingered Manipulation*, Springer, 2009. - S. Arimoto, “Intelligent control of multi-fingered hands,”
*Annual Reviews in Control*, vol. 28, no. 1, pp. 75–85, 2004. View at Google Scholar - M. Yoshida, S. Arimoto, and J. H. Bae, “Blind grasp and manipulation of a rigid object by a pair of robot fingers with soft tips,” in
*Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '07)*, pp. 4707–4714, Rome, Italy, April 2007. View at Publisher · View at Google Scholar · View at Scopus - S. Arimoto, P. T. A. Nguyen, H. Y. Han, and Z. Doulgeri, “Dynamics and control of a set of dual fingers with soft tips,”
*Robotica*, vol. 18, no. 1, pp. 71–80, 2000. View at Publisher · View at Google Scholar · View at Scopus - M. L. Latash and M. L. Zatsiorsky, “Principle of superposition in human prehension,” in
*Advances in Robot Control*, pp. 249–261, Springer, 2006. View at Google Scholar - S. Arimoto, “Reduction of complexity in learning dexterous multi-fingered motions: a theoretical exploration into a future problem C. E. Shannon raised,”
*Communications in Information and Systems*, vol. 1, no. 1, pp. 1–14, 2001. View at Google Scholar · View at Zentralblatt MATH - H. Goldstein, C. P. Poole, and J. L. Safko,
*Classical Mechanics*, Addison Wesley, 3rd edition, 2002. - F. L. Lewis and C. T. Abdallaah,
*Control of Robot Manipulators*, Macmillan, New York, NY, USA, 1994. - J. Baumgarte, “Stabilization of constraints and integrals of motion in dynamical systems,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 1, pp. 1–16, 1972. View at Publisher · View at Google Scholar · View at Zentralblatt MATH