Special Issue Article
Concise method to the dynamic modeling of climbing robot
Advances in Mechanical Engineering 2017, Vol. 9(2) 1–10 Ó The Author(s) 2017 DOI: 10.1177/1687814017691670 journals.sagepub.com/home/ade
Yaru Xu and Rong Liu
Abstract With the aim of dynamic modeling of the climbing robot with dual-cavity structure and wheeled locomotion mechanism, a succinct and explicit equation of motion based on the Udwadia–Kalaba equation is established. The trajectory constraint of the climbing robot, which is regarded as the external constraint of the system, is integrated into the dynamic equation dexterously. A modified numerical method is considered to reduce the errors because the numerical results obtained by integrating the constrained dynamic equation yield the errors. The trajectories are almost coincident by comparing the modified numerical value and the theoretical value. The driving torques required to guarantee the climbing robot to move along the given trajectory are obtained explicitly, which overcomes the disadvantage of obtaining dynamical equation from traditional Lagrange equation by Lagrange multiplier. The simulations are performed to demonstrate that the dynamical equation established by this method with brevity and accuracy is in accordance with reality status. Keywords Dynamic modeling, climbing robot, Udwadia–Kalaba equation, trajectory constraints, errors reducing
Date received: 29 September 2016; accepted: 11 January 2017 Academic Editor: Chuanzeng Zhang
Introduction Climbing robots have been a very attractive research topic since there are various potential applications to increase operational efficiency and protect human health and safety in environments such as the exteriors of buildings, bridges or dams storage tanks, nuclear facilities, and reconnaissance within building. In general, climbing robots can be categorized based on the surface adhesion (e.g. magnetic type,1 gripping type,2 rail guided type,3 biomimetic type,4 and suction type5) and the locomotion mechanism (e.g. legged type,6 tracked type,7 translation type,3 cable-driven type,8 and wheeled type9), which are the two major issues in the design of climbing robots. Suction adhesion, which is the most commonly used adhesion method, can be widely adopted for less rough surfaces because it enables strong attachment to the surfaces regardless of materials such as glass, ceramics tiles, and cement. Using wheeled type as the locomotive mechanism, the
robots have high moving speed, simple structure, and feasible control system. For climbing robot, the accurate kinematics and dynamical models are the basis for completing all kinds of tasks. F Xu et al.10 derived the kinematics and dynamical models for the obstacle-climbing capabilities of the driving and driven wheels of the robot; WR Provancher et al.11 developed the dynamical model of the city-climber robot that has the capability to move on floors, climb walls, walk on ceilings, and transit between them when it travels on different surfaces. WH Ko et al.12 simplified and redeveloped the reducedorder two-arm dynamic model on climbing robot (i.e. not computer-aided design (CAD) model or kinematic Institute of Robotics, Beihang University, Beijing, China Corresponding author: Yaru Xu, Institute of Robotics, Beihang University, Beijing 100191, China. Email:
[email protected]
Creative Commons CC-BY: This article is distributed under the terms of the Creative Commons Attribution 3.0 License (http://www.creativecommons.org/licenses/by/3.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/ open-access-at-sage).
2 model) based on Lagrangian mechanics to explore the relation between the individual model parameters and the resultant dynamic behavior of the model. S Nam et al.13 derived the dynamics model of a compliant multibody climbing robot with magnetic adhesion using the Lagrangian formulation. However, the processes of above-mentioned dynamic modeling are complex and complicated. Obtaining equations of motion for constrained mechanical systems is one of the central issues in multibody dynamics.14 Constrained motion was initially described by Lagrange in 1787 using Lagrange multiplier method, which relies on problem-specific approaches to determine the multipliers that are often very difficult to find. Gauss’ principle was proposed by Gauss in 1829 which gives a clear description of the general nature of constrained motion through minimization of a function of the accelerations of the particles of a system. Gibbs15 and Appell16 had independently developed equations of constrained motion when the constraints satisfy d’Alembert’s principle in 1879 and 1899, respectively. However, Gibbs–Appell equations require specific quasi-coordinates. Dirac17 developed how to determine the Lagrange multipliers using Poisson brackets for singular Hamiltonian systems in 1964, in which the constraints do not exactly depend on time. There are some marvelous discoveries, called the Udwadia–Kalaba equation, proposed by Firdaus E. Udwadia and Robert E. Kalaba from the University of Southern California in 1992. The significant advantages are threefold. (1) The governing motion equations of multibody system subject to holonomic constraints are proposed based on d’Alembert’s principle and Gauss’ principle.18,19 (2) Given that the system with nonideal constraints may not satisfy d’Alembert’s principle, explicit set of equations under the case of nonholonomic constraints are added to perfect the previous motion equations of multibody system.20,21 (3) General and explicit equations of motion to handle systems whether or not their mass matrices are singular are developed in the case of singular mass matrices arise. The new theory opens up a new way of modeling complex multibody systems. An additional perspective has been proposed by Udwadia and Kalaba which is useful to help us understand nature’s law from new points of view. The explicit equations of constraint force without Lagrange multiplier can be obtained according to Udwadia–Kalaba theory, which is a simple, aesthetic, and thought-provoking description of the world at a very fundamental level, to realize a great breakthrough in the field of analytical mechanics.18,22–24 Udwadia and colleagues24–27 also successfully addressed the dynamics and control of nonlinear uncertain systems
Advances in Mechanical Engineering using the fundamental equation to provide considerable reference value for the basic research and the further application of this method. Many researchers have studied and established dynamical model in some application fields by employing the Udwadia–Kalaba approach, such as satellite systems,28 industry mechanical arm,29 parallel manipulator,30 flexible multibody systems,31 machine fish,14 heavenly bodies’ movements (especially on Kepler’s law and the inverse square law of gravitation),32 and falling cat’s movements.33 In view of its strong attachment to the surfaces, high moving speed and simple control system, climbing robot based on dual-cavity structure and wheeled locomotion mechanism is studied in this article. Like many other robots, climbing robots are often required to move along some specified trajectories (i.e. trajectory constraint). However, only a few scholars focus on the explicit dynamical equations for climbing robot subject to constrains. In this article, a succinct and exact dynamic modeling of constrained climbing robot is presented by applying the Udwadia–Kalaba approach. However, the numerical results obtained by integrating the constrained dynamic equation yield the errors. Therefore, a modified numerical method is considered to reduce the errors. Numerical simulations are performed to demonstrate the efficacy and accuracy of this method. Finally, conclusions of this work and the future tasks are given.
Udwadia–Kalaba equation This section shows how to obtain the explicit equations of motion for a constrained dynamical system using a simple straightforward three-step procedure:18,24 1.
2. 3.
In terms of the generalized coordinates, equations of motion of the unconstrained dynamical system written using Newtonian or Lagrangian mechanics are considered. Trajectory constraint required to model of the given constrained system is described. Additional generalized forces of constraint imposed on the system are expressed.
Step 1: unconstrained dynamics Consider an unconstrained dynamical system, the configuration of which is uniquely specified by the n generalized coordinates q :¼ ½q1 q2 qn T , in which the superscript T represents the transpose of a vector or a matrix. Using Newtonian or Lagrange equation, the equation of unconstrained motion of the system can be expressed in the following form24
Xu and Liu
3 _ tÞ Mðq, tÞ€q = Qðq, q,
ð1Þ
with the initial conditions qð0Þ = q0 , q_ ð0Þ = q_ 0
ð2Þ
in which M(q, t) 2 Rn 3 n is a mass matrix, which is a _ t) 2 Rn collects the normal function of q and t; Q(q, q, and Coriolis inertial terms and the applied forces and is _ and t; q 2 Rn is the generala known function of q, q, ized coordinate n-vector; q_ 2 Rn is the velocity; €q 2 Rn is the acceleration. From equation (1), the generalized acceleration of the unconstrained dynamical system is _ t), which is of the form denoted by a(q, q, € _ tÞ :¼ aðq, q, _ tÞ q = M1 ðq, tÞQðq, q,
ð3Þ
Step 2: constraint description The system is subjected to the p sufficiently smooth control requirements as constraints provided by _ tÞ = 0, ui ðq, q,
i = 1, 2, . . . , p
ð4Þ
These p constraints include all the usual varieties of holonomic and/or nonholonomic constraints and then some.24Equation (4) can be simplified as F = u1
u2
up
T
ð5Þ
Differentiating the usual constraint equation (5) in Lagrangian mechanics which are usually in Pfaffian form with respect to time t once (for nonholonomic constraints) or twice (for holonomic constraints) yields the following set of consistent constraint equations24 _ tÞ€q = bv ðq, q, _ tÞ Aðq, q,
ð6Þ
_ t) 2 R is referred to as constraint in which A(q, q, _ t) is a p-vector. matrix and bv (q, q, The second-order constraint is linear in €q so that it is believed to be the most appropriate form for further dynamic analysis and ensure the completeness of the information. p3n
Step 3: constraint force description Additional constrained forces arise due to the constraints applied to the unconstrained system. Accordingly, the actual explicit equation of motion of the constrained system could be assumed to take the form
recently, a simple and insightful approach was proposed by Udwadia and colleagues21–27 based on Gauss’ principle of least constraint, which gives the constraint _ t) explicitly as21,22,34,35 force Qc (q, q, _ tÞ Qc ðtÞ :¼ Qc ðq, q,
1 _ tÞ Aðq, q, _ tÞM1 ðqÞQðq, q, _ tÞ = M2 (q)B+ ðq, tÞ bv ðq, q,
ð8Þ 12
_ t)M (q), and the superscript in which B(q, t) = A(q, q, ‘‘ + ’’ represents the generalized Moore–Penrose _ t), which miniinverse. The exact control force Qc (q, q, _ t)T M1 ½Qc (q, q, _ t) at each mizes the quantity ½Qc (q, q, instant of time t, enables all the constraints to be exactly satisfied at every instant of time t and can be derived explicitly by equation (8). Accordingly, the explicit dynamical equation of the constrained system is shown in the following general form _ tÞ MðqÞ€q = Qðq, q, 1 + 2 _ tÞ Aðq, q, _ tÞM1 ðqÞQðq, q, _ tÞ + M ðqÞB ðq, tÞ bv ðq, q,
ð9Þ In this article, equation (9) is referred to as the Udwadia–Kalaba equation. Premultiplying both sides of equation (9) with M1 , the acceleration of the constrained system can be obtained by 1 1 + €q :¼ a + M1 Qc = M1 Q + M2 AM2 bv AM1 Q
Dynamic modeling of climbing robot This section develops the dynamic model of the climbing robot with dual-cavity structure and wheeled locomotion mechanism. Dual cavities are adsorbed on the wall with negative pressure absorption. As the driving wheels, the two rear wheels are driven by two DC servo motors, respectively, while the power is transmitted to the two front wheels by synchronous belt. Four-wheel move steering mechanisms are employed to enhance the driving forces on the wall. Several important assumptions are made to simplify the modeling: 1.
ð7Þ
2.
_ t) 2 Rn is the additional constraint in which Qc (q, q, force imposed on the system to ensure the control requirements in equation (6) are satisfied. There are _ t), and already many ways to determine Qc (q, q,
3. 4.
_ tÞ _ tÞ + Qc ðq, q, Mðq, tÞ€q = Qðq, q,
ð10Þ
The climbing robot is considered to be a single rigid body. The center of mass is in the geometric center of the climbing robot. There is no-slip motion on the wheels. The moment of inertia of the wheels and the rolling friction between the wheels and the wall are ignored.
4
Advances in Mechanical Engineering wall; P is the pressure difference between the atmospheric pressure and the pressure inside the seal ring; and S is the effective adsorption area of the seal ring. The velocity of the climbing robot in the generalized coordinates q = ½x, y, uT can be expressed as v = x_ cos u + y_ sin u
ð12Þ
and the acceleration of the climbing robot can be represented by v_ = €x cos u + €y sin u x_ sin uu_ + y_ cos uu_
in which x_ and y_ are the velocities of the mass center of the climbing robot in the X- and Y-directions, respectively; and €x and €y are the accelerations of the mass center of the climbing robot in the X- and Y-directions, respectively. Substituting equation (13) into equation (11), the following equation is derived
Figure 1. Climbing robot subject to trajectory constraint.
5.
6.
The friction between the seal ring and the wall distributes on the center of the robot with uniformity and symmetry. The climbing robot is symmetrical, so the force analysis model is simplified as a twodimensional (2D) model.
Consider an inertial frame of reference XOY in which a climbing robot of mass m (with center of mass is at o) moves on the 2D surface (see Figure 1). The generalized coordinates q = ½x, y, uT are denoted as the configuration of the mass center of the climbing robot. In which v is the velocity of the climbing robot; u is the orientation angle of the climbing robot; FR and FL are the forces produced by the right and left motors, respectively; f is the friction between the seal ring and the wall; G1 is the gravity component in the v-direction; J is the moment of inertia around the mass center of the climbing robot; and d is the distance between the right or left wheel and the mass center of the climbing robot.
_ x mr cos uu_ _y mr cos u€x + mr sin u€y = mr sin uu_ ð14Þ + t R + tL fr mgr sin u 2. The momentum balance can be written as
J €u = FR d FL d
ð15Þ
Jr € u = tR tL d
ð16Þ
or, equivalently
3. The structural constraint of the climbing robot is expressed as
x_ sin u y_ cos u = 0
1. According to Newton’s second law of motion, the following relation can be obtained
ð11Þ
The different forces in equation (11) can be given as FR = tR =r, FL = t L =r, G1 = mg sin u, f = (1 k)mM FS , and FS = PS, in which tR and t L are the driving torques of the right and left motors, respectively; r is the radius of the wheels; k is the ratio of the pressure on the wheels to the pressure on the whole robot; mM is the coefficient of rolling resistance between the tracked wheels and
ð17Þ
Differentiating equation (17) with respect to t yields _ x sin uu_ _y sin u€x cos u€y = cos uu_
Explicit dynamic modeling
m_v = FR + FL G1 f
ð13Þ
ð18Þ
Accordingly, equations (14), (16), and (18) can be expressed in matrix equation 2
mr cos u 6 0 4
mr sin u 0
32 3 €x Jr 76 7 y5 d 54 € € 0 u 0
sin u cos u 2 3 _ x mr cos uu_ _ y fr mgr sin u mr sin uu_ 6 7 =4 5 ð19Þ 0 _ _ cos uu_x sin uu_y 2 3 1 1 6 7 tR + 4 1 1 5 tL 0 0
Xu and Liu
5
Constraint force description
Equation (19) can be simplified as _ tÞ _ tÞ + Qc ðq, q, MðqÞ€q = Qðq, q,
ð20Þ
Using equation (8), the constraint torque Qc , which is equivalent to the constraint force, can be derived as
in which 2 3 2 3 2 3 €x x x_ 6 7 _ 6 7 € 6 7 q = 4 y 5, q = 4 y_ 5, q = 4 €y 5 €u u u_ 2 3 mr cos u mr sin u 0 6 Jr 7 7 0 0 MðqÞ = 6 4 d5 sin u cos u 0 2 3 _ _ y fr mgr sin u mr sin uu_x mr cos uu_ 6 7 Q=4 5 0 _ x sin uu_ _y cos uu_ 2 3 1 1 6 7 tR c Q = 4 1 1 5 tL 0 0 in which Qc is the constraint force to satisfy the trajectory constraint as described below.
Trajectory constraint description The climbing robot moves on the wall along the following trajectory to perform specific tasks. The constrained trajectory is expressed as
x = 10 sin t y = 6 cos t
ð21Þ
Differentiating equation (21) with respect to time t twice yields the second-order constrains
€x = 10 sin t €y = 6 cos t
ð22Þ
1 1 + bv AM1 Q Qc = M2 AM2 2 3 1 1 6 7 tR = 4 1 1 5 = dt tL 0 0
Therefore, the driving torques of the right and left motors t based on Udwadia–Kalaba equation can be obtained as 1 1 + bv AM1 Q t = d+ M2 AM2
Assuming the initial configurations are x(0) = 0, _ = 6=10. y(0) = 6, u(0) = 0; x_ (0) = 10, y_ (0) = 0, u(0) The parameters for simulation are m = 5 kg, J = 0:33 kg m2 , r = 0:027 m, d = 0:3 m, mM = 0:8, k = 0:53, P = 2:5 kPa, and S = 0:06 m2 . The solution of the equations can be obtained through ode45 algorithm in MATLAB 2010b. The error tolerance is set to 10212 and the simulation time is 20 s. The explicit analytic results given in Figures 2–6 are verified by numerical simulations: 1.
Figures 2 and 3 represent x-coordinate curve and y-coordinate curve as a function of time t when a certain level driving torque obtained from equation (25) is generated by the two rear wheels to realize the desired trajectory, in which the solid curves and the dashed curves represent the numerical value and the theoretical value,
ð23Þ
in which A=
1 0
0 1
0 0
2 3 €x €q = 4 €y 5 €u and bv =
10 sin t 6 cos t
ð25Þ
Result and simulation analysis
Equation (22) can be simplified as A€q = bv
ð24Þ
Figure 2. Comparison of x-coordinate curves between numerical value and theoretical value.
6
Advances in Mechanical Engineering
Figure 3. Comparison of y-coordinate curves between numerical value and theoretical value.
Figure 5. Error curve in y-coordinate between numerical value and theoretical value.
Figure 4. Error curve of x-coordinate between numerical value and theoretical value.
Figure 6. Comparison of trajectories between numerical value and theoretical value.
2.
respectively. x-num and y-num are numerical values of x-coordinate and y-coordinate, respectively. There are obvious errors as a function of time t for the trajectory constraint described by equation (21) between numerical value and theoretical value, even though the solid curve approximates the dashed curve in Figures 2 and 3, as shown in Figures 4 and 5. x-error and yerrors are errors between numerical value and theoretical value of x-coordinate and y-coordinate, respectively. Compared with x = 10 m and y = 6 m, however, the obvious displacement errors, which is about 1.8 3 1023 m in the xdirection and 0.6 m in the y-direction, still exist and increase with time.
3.
The numerical trajectories are slightly off the theoretical trajectory with time, as illustrated in Figure 6, in which the solid curve and the dashed curve represent the numerical trajectory and the theoretical trajectory, respectively. Therefore, it is necessary to consider the numerical method for reducing the errors.
Errors reducing Any numerical integration scheme on the second-order differential equation needs a set of two first-order differential equations.36 One can define the constrained acceleration, €q = v_ , and rewrite the equation of motion based on the second-order differential equation of equation (10) in the following first-order form
Xu and Liu
7 q_ = v 1
v_ = a + M2 B+ ðbv AaÞ
ð26Þ
It is known that the responses of the constrained dynamic system must satisfy the constraint equations at all times during numerical integration of the dynamic equation. However, it is found that the numerical results tend to deviate from the constrained path as shown in Figure 6. The basic idea for this numerical algorithm for reducing errors is to include the constraint effect into equation (26). It is considered that the errors appeared because the derived dynamic equation (see equation (6)) considers only the acceleration-based constraint and neglects the velocity-based constraint of given trajectory constraints that differentiate the constraint equations once with respect to time. It can be expressed as Aq_ = bq
Figure 7. Comparison of x-coordinate curves between modified numerical value and theoretical value.
ð27Þ
With the aim of incorporating this error source, one can add a new term considered as a kinematic correction term to equation (26) to compensate for numerical errors along the integration37 q_ = v + M1 AT m
ð28Þ
Substituting equation (28) into equation (27), one obtains 1 bq Av m = AM1 AT
ð29Þ
Then, the dynamic equation can be modified as q_ = v + M1=2 B+ bq Av v_ = a + M1=2 B+ ðbv AaÞ
ð30Þ
in which q_ is the constrained velocity, v is the constrained acceleration, and a = M1 Q is the unconstrained acceleration. Differentiating equation (21) with respect to time once, one obtains
1 0 0 A= 0 1 0 2 3 x_ q_ = 4 y_ 5 u_ 10 cos t bq = 6 sin t The dynamic responses and constraint forces can be obtained by substituting the dynamic equation of the unconstrained system and the constraint equations into equation (30). The explicit analytic results given in Figures 7–11 are verified by modified numerical
Figure 8. Comparison of y-coordinate curves between modified numerical value and theoretical value.
simulations. The required driving torques generated by motors enable the climbing robot to move along the ellipse trajectory presented above accurately: 1.
2.
Figures 7 and 8 represent the comparison of xcoordinate curves and y-coordinate curves as a function of time t between modified numerical value and theoretical value, respectively. x-numM and y-num-M are modified numerical values of x-coordinate and y-coordinate, respectively. Comparing with the x-error between numerical value and theoretical value, the amplitude of xerror between modified numerical value and numerical value increases but not continuously increases with time, as shown in Figure 9. It is observed that the y-errors’ growth trend
8
Advances in Mechanical Engineering
Figure 9. Comparison of x-error curves between modified numerical value and numerical value.
Figure 10. Comparison of y-error curves between modified numerical value and numerical value.
3.
obtained by the modified numerical integration is improved obviously and decreased by one order of magnitude, as shown in Figure 10. x-error-M and y-error-M are errors between modified numerical value and theoretical value of x-coordinate and y-coordinate, respectively. In general, the errors are within an acceptable range that the constraint relations in the x-direction and y-direction are being maintained well as desired. The trajectories are almost coincident by comparing modified numerical value and numerical value, as illustrated in Figure 11, in which the solid curves and the dashed curves represent the modified numerical value and the numerical value, respectively. The overall effect of constraint trajectory obtained by modified
Figure 11. Comparison of trajectories between modified numerical value and theoretical value.
Figure 12. Constraint torque curves.
4.
numerical integration gets clear improvement; still, the errors do not perfectly disappear. The constraint torques Qc1 and Qc2 as a function of time t which are required to satisfy the trajectory constraint of the climbing robot are shown in Figure 12. The driving torques TR and TL as a function of time t generated by the right and the left motors to realize the theoretical trajectory are shown in Figure 13.
Conclusion and future work With the aim of dynamic modeling of the climbing robot with dual-cavity structure and wheeled locomotion mechanism, this study establishes the dynamic equation based on the Udwadia–Kalaba equation. The research process results in several conclusions as follows:
Xu and Liu
9 Funding The author(s) received no financial support for the research, authorship, and/or publication of this article.
References
Figure 13. Driving torque curves generated by the motors.
1.
2.
3.
The trajectory constraint of the climbing robot is regarded as the external constraint of the system and integrated into the dynamic equation dexterously based on the idea of the Udwadia–Kalaba equation. The consideration of velocity-based constraint can reduce the error to guarantee that the trajectories obtained by modified numerical integration are close to the numerical trajectory as possible. The driving torques required to guarantee the climbing robot to move along the given trajectory are obtained explicitly by solving Udwadia–Kalaba equation, which overcomes the disadvantage of obtaining dynamical equation from traditional Lagrange equation by Lagrange multiplier effectively. The methodology described in this article can also be applied to many other kinds of climbing robots, such as climbing robots with magnetic adhesion or tracked type. In addition, the trajectory constraints are geometrical shapes of many kinds.
However, the simulation in this article is executed on the condition which the initial condition satisfies the constrained equation. Therefore, the future work is the dynamics modeling of climbing robot which the initial condition does not satisfy the constrained equation and the error reduction of the dynamic modeling by further study. Moreover, in numerous factors, the friction, no doubt, is one of the biggest nonlinear factors in the dynamic modeling of climbing robot. It requires in-depth discussion to obtain the precise model of the friction. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
1. Han SC, Kim J and Yi HC. A novel design of permanent magnet wheel with induction pin for mobile robot. Int J Precis Eng Man 2009; 10: 143–146. 2. Balaguer C, Gimenez A, Huete AJ, et al. The MATS robot: service climbing robot for personal assistance. IEEE Robot Autom Mag 2006; 13: 51–58. 3. Lee D, Lee S, Ku N, et al. Development of a mobile robotic system for working in the double-hulled structure of a ship. Robot Cim: Int Manuf 2010; 26: 13–23. 4. Kwak JS and Kim TW. A review of adhesion and friction models for gecko feet. Int J Precis Eng Man 2010; 11: 171–186. 5. Koo IM, Trong TD, Lee YH, et al. Development of wall climbing robot system by using impeller type adhesion mechanism. J Intell Robot Syst 2013; 72: 57–72. 6. Chen X, Wang L, Ye X, et al. Prototype development and gait planning of biologically inspired multi-legged crablike robot. Mechatronics 2013; 23: 429–444. 7. Gao X, Xu D, Wang Y, et al. Multifunctional robot to maintain boiler water-cooling tubes. Robotica 2009; 27: 941–948. 8. Akinfiev T, Armada M and Nabulsi S. Climbing cleaning robot for vertical surfaces. Ind Robot 2009; 36: 352–357. 9. Tavakoli M and Viegas C. Analysis and application of dual-row omnidirectional wheels for climbing robots. Mechatronics 2014; 24: 436–448. 10. Xu F, Shen J and Jiang GP. Kinematic and dynamic analysis of a cable-climbing robot. Int J Adv Robot Syst. Epub ahead of print 23 July 2015. DOI: 10.5772/60865. 11. Provancher WR, Jensen-Segal SI, Fehlberg MA, et al. ROCR: an energy-efficient dynamic wall-climbing robot. IEEE/ASME T Mech 2011; 16: 897–906. 12. Ko WH, Chiang WH, Hsu YH, et al. A modelbased two-arm robot with dynamic vertical and lateral climbing behaviors. J Mech Robot 2016; 8: 044503-1–044503-9. 13. Nam S, Oh J, Lee G, et al. Dynamic analysis during internal transition of a compliant multi-body climbing robot with magnetic adhesion. J Mech Sci Technol 2014; 28: 5175–5187. 14. Zhao H, Zhen SC and Chen YH. Dynamic modeling and simulation of multi-body systems using the Udwadia– Kalaba theory. Chin J Mech Eng 2013; 26: 839–850. 15. Gibbs JW. On the fundamental formulae of dynamics. Am J Math 1879; 2: 49–64. 16. Appell P. Sur une Forme Generale des Equations de la Dynamique. Cr Acad Sci I: Math 1899; 129: 459–460. 17. Dirac PAM. Lectures on quantum mechanics. New York: Yeshiva University, 1964. 18. Schutte AD and Udwadia FE. New approach to the modeling of complex multibody dynamical systems. J Appl Mech 2011; 78: 856–875.
10 19. Udwadia FE and Kalaba RE. Analytical dynamics: a new approach. Cambridge: Cambridge University Press, 1996. 20. Udwadia FE and Kalaba RE. On the foundations of analytical dynamics. Int J Nonlin Mech 2002; 37: 1079–1090. 21. Udwadia FE and Kalaba RE. What is the general form of the explicit equations of motion for constrained mechanical systems? J Appl Mech 2002; 69: 335–339. 22. Udwadia FE and Kalaba RE. A new perspective on constrained motion. P Roy Soc A: Math Phy 1992; 439: 407–410. 23. Udwadia FE and Schutte AD. Equations of motion for general constrained systems in Lagrangian mechanics. Acta Mech 2010; 213: 111–129. 24. Udwadia FE and Wanichanon T. Control of uncertain nonlinear multibody mechanical systems. J Appl Mech 2014; 81: 041020-1–041020-11. 25. Udwadia FE and Wanichanon T. A closed-form approach to tracking control of nonlinear uncertain systems using the fundamental equation. Earth Space 2012; 10: 1339–1348. 26. Udwadia FE, Wanichanon T and Cho H. Methodology for satellite formation-keeping in the presence of system uncertainties. J Guid Control Dynam 2014; 37: 1611–1624. 27. Udwadia FE and Koganti PB. Dynamics and control of a multi-body planar pendulum. Nonlinear Dynam 2015; 81: 845–866. 28. Cho H and Yu A. New approach to satellite formationkeeping: exact solution to the full nonlinear problem. J Aerospace Eng 2009; 22: 445–455.
Advances in Mechanical Engineering 29. Liu J and Liu R. Simple method to the dynamic modeling of industrial robot subject to constraint. Adv Mech Eng 2016; 8: 1–9. 30. Huang J, Chen YH and Zhong Z. Udwadia–Kalaba approach for parallel manipulator dynamics. J Dyn Syst: T ASME 2013; 135: 1012–1030. 31. Pennestri E, Valentini PP and de Falco D. An application of the Udwadia–Kalaba dynamic formulation to flexible multibody systems. J Frankl Inst 2010; 347: 173–194. 32. Zhang B, Zhen S, Zhao H, et al. A novel study on Kepler’s law and inverse square law of gravitation. Eur J Phys 2015; 36: 035018. 33. Zhen SC, Huang K, Zhao H, et al. Why can a free-falling cat always manage to land safely on its feet? Nonlinear Dynam 2015; 79: 2237–2250. 34. Udwadia FE. Equations of motion for constrained multibody systems and their control. J Optimiz Theory App 2005; 127: 627–638. 35. Kalaba RE and Udwadia FE. Equations of motion for nonholonomic, constrained dynamical systems via Gauss’s principle. J Appl Mech 1993; 60: 662–668. 36. Rahmatalla S, Lee ET and Eun HC. Numerical integration scheme to reduce the errors in the satisfaction of constrained dynamic equation. J Mech Sci Technol 2013; 27: 941–949. 37. Braun DJ and Goldfarb M. Eliminating constraint drift in the numerical simulation of constrained dynamical systems. Comput Method Appl M 2009; 198: 3151–3160.
