Development of an Educational Code of Deriving Equations of ‎Motion and Analyzing Dynamic Characteristics of Multibody ‎Closed Chain Systems using GNU Octave for a Beginner

Document Type : Research Paper


Department of Mechanical Engineering, Yuhan University, 590, Gyeongin-ro, Bucheon-si, Gyeonggi-do, Republic of Korea


In this study, an automatic GNU Octave code, a free high-level language, for the educational purposes was developed to derive equations of motion and constrain equations of a multibody closed chain system and to calculate the response of the system. The code for calculating the dynamic response was developed by formulating several equations in symbolic expression and extracting differential-algebraic equations in matrix form. The code has a similar structure to the previous code for the open chained system, but it deals with the constraint equation and different numerical integration. The examples of closed chain systems provide an additional procedure to derive the constraint equations by using Lagrangian multiplication theory and to solve the differential-algebraic equations using the Runge-Kutta method. The code was made to understand the theory of analysis and the structure of calculation easily. In addition, the code has an automatic process of the derivation of the Lagrange equation and the constraint equations in matrix form after inputting the number of symbolic information such as position and velocity coordinates and design variables of the system that the user wants to review. The code was validated by comparing the dynamic response of the four-bar linkage with the same design variables and initial conditions of the previous work. By using the code, the reader's ability to exchange information such as symbols and matrices will be expected to be improved.


Main Subjects

Publisher’s Note Shahid Chamran University of Ahvaz remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

[1] Yoo, W. S., Park, S. J., Dmitrochenko, O., Pogorelov, D., Verification of Absolute Nodal Coordinate Formulation in Flexible Multibody Dynamics via Physical Experiments of Large Deformation Problems, Journal of Computational and Nonlinear Dynamics, 1(1), 2006, 81-93.
[2] Yoo, W. S., Kim, T. Y., Jung, S., History and Future of Multi-body Dynamics, Transactions of the Korean Society of Mechanical Engineers - A, 43(7), 2019, 483-491.
[3] Munro, N., Symbolic methods in control system analysis and design, Institution of Engineering and Technology, Stevenage, 1999.
[4] Van Khang, N., Kronecker product and a new matrix form of Lagrangian equations with multipliers for constrained multi-body systems, Mechanics Research Communications, 38(4), 2011, 294-299.
[5] Negrut, D., Dyer, A., ADAMS/Solver Primer, MSC Software Documentation, Ann Arbor, 2004.
[6] Bonev, I. A., Geometric analysis of parallel mechanisms, Université Laval, Qubec City, 2002.
[7] vander Linde, R. Q., Schwab, A. L., Lecture Notes Multi-body Dynamics B, Delft University, Delft, 1998.
[8] Sadati, S. M. H., Naghibi, S. E., Naraghi, M., An automatic algorithm to derive linear vector form of Lagrangian equation of motion with collision and constraint, Procedia Computer Science, 76, 2015, 217-222.
[9] Sadati, S. M. H., Naghibi, S. E., Shiva, A., Zschaler, S., Hauser, H., Walker, I., Althoefer, K., Nanayakkara, T., AutoTMTDyn: A Matlab software package to drive TMT Lagrange dynamics of series rigid-and continuum-link mechanisms, IROS 2018 Workshop on Soft Robotic Modeling and Control: Bringing Together Articulated Soft Robots and Soft-Bodied Robots, Madrid, Spain, 2018.
[10] de Jalón, J. G., Callejo, A., A straight methodology to include multibody dynamics in graduate and undergraduate subjects, Mechanism and Machine Theory, 46(2), 2011, 168-182.
[11] Park, Y. H., An Automatic Program of Generation of Equation of Motion and Dynamic Analysis for Multi-body Mechanical System Using GNU Octave, Journal of Applied and Computational Mechanics, 7(3), 2021, 1687-1697.
[12] Yoo, H. H., Kane's methodology for derivation of equation of motion, Journal of KSNVE, 14(4), 2004, 40-47.
[13] Yang, K. D., Bae, D. S., Yang, S. M., Choi, C. K., An Index 2 Differential-Algebraic Equation Formulation for Multibody System Dynamics, Transactions of the Korean Society of Mechanical Engineers, 19(11), 1995, 2769-2775.
[14] Lee, D. C., Lee, S. H., Han, C. S., A Formulation of the Differential Equation on the Equations of motion and Dynamic Analysis for the Constrained Multibody Systems, Transaction of the Korean Society of Automotive Engineers, 5(1), 1997, 154-161.
[15] Eaton, J. W., Bateman, D., Hauberg, S., Gnu octave, London: Network theory, London, 1997.
[16] Ateş, A., Yeroğlu, C. Two degrees of freedom FOPID control loop design via SMDO algorithm, Pamukkale University Journal of Engineering Sciences, 22(8), 2016, 671-676.
[17] Arda, M., Dynamic analysis of a four-bar linkage mechanism, Machines. Technologies. Materials, 14(5), 2020, 186-190.