Discrete and Continuous Models and Applied Computational ScienceDiscrete and Continuous Models and Applied Computational ScienceDiscrete and Continuous Models and Applied Computational Science2658-46702658-7149Peoples' Friendship University of Russia named after Patrice Lumumba (RUDN University)8417ArticlesСтатьиResearch ArticleSolving Differential Equations of Motion for Constrained Mechanical SystemsРешение дифференциальных уравнений движения для механических систем со связямиMukharlyamovR GМухарлямовРоберт ГарабшевичDepartment of Theoretical MechanicsКафедра теоретической механикиrobgar@mail.ruBeshawA WБешауА ВDepartment of MathematicsКафедра математикиassayewalelgn@gmail.comPeoples’ Friendship University of RussiaРоссийский университет дружбы народовBahr Dar UniversityБахрдарский университет150320133NO3 (2013)№3 (2013)819108092016Copyright ©; 2013, Мухарлямов Р.Г., Бешау А.В.2013Мухарлямов Р.Г., Бешау А.В.http://creativecommons.org/licenses/by/4.0https://journals.rudn.ru/miph/article/view/8417This paper presents an investigation of modeling and solving system of differential equations in the study of mechanical systems with holonomic constraints. A method is developed for constracting equation of motion for mechanical system with constraints. A technique is developed how to approximate the solution of the problem that is obtained from modeling of kinematic constraint equation which is stable. A perturbation analysis shows that velocity stabilization is the most efficient projection with regard to improvement of the numerical integration. How frequently the numerical solution of the ordinary differential equation should be stabilized is discussed. A procedure is indicated to get approximate solution when the systems of differential equations can’t be solved analytically. A new approach is applied for constructing and stabilyzing Runge-Kutta numerical methods. The Runge-Kutta numerical methods are reformulated in a new approach. Not only the technique of formulation but also the test developed for its stability is new.Finally an example is presented not only to demonstrate how the stability of the solution depends on the variation of the factor but also how to find an approximate solution of the problem using numerical integration.В работе рассматривается задача построения систем дифференциальных уравнений по известным частным интегралам. Приводится метод определения правых частей систем дифференциальных уравнений, основанный на определении общего решения системы линейных алгебраических уравнений с прямоугольной матрицей коэффициентов. Предлагается использовать для численного решения построенной системы дифференциальных уравнений метод Рунге-Кутта. Для рассматриваемой задачи ранее были использованы простейшие разностные схемы первого порядка и метод Рунге Кутта для случая линейных дифференциальных уравнений возмущений связей с постоянными коэффициентами. В статье получены ограничения на коэффициенты уравнений возмущений связей, зависящие от фазовых координат системы, при решении дифференциальных уравнений методом Рунге-Кутта. Подробно рассмотрены случаи разностных уравнений первого порядка, состоящих из нескольких стадий. Получена общая форма условий стабилизации уравнений связей. Метод иллюстрируется на примере решения кинематической задачи кривошипно-шатунного механизма.numerical integrationkinematic constraintstable solutionTaylor expansionrow decompositionдифференциальные уравнениячисленное интегрированиекинематические ограничениястабилизацияряд Тейлора1.Arabyan A., Wu F. An Improved Formulation for Constrained Mechanical Systems // Kluwer Academic Pub. — 1998. — Pp. 49–69.2.Muharlyamov R. Equations of Motion of Mechanical Systems. — PFUR, 2001. — (in russian).3.Gonze D. Numerical Methods to Solve Ordinary Differential Equations // P.J. Brief Bioinform. — 2009. — Pp. 53–64.4.Lakoba T. Runge-Kutta Methods. — University of Vermont, 2006. — Pp. 15–20.