Определение ускорения центра масс тела

Размещено на http://www.allbest.ru/

Южно-Российский государственный политехнический университет имени М.И. Платова, г. Новочеркасск

Определение ускорения центра масс тела

С.А. Алексеев

О.А. Бурцева

А.С. Алексеева

Н.Р. Абуладзе

На основе информации о полных ускорениях трех произвольных точек абсолютно твердого тела, решается задача о нахождении мгновенных угловых скоростей и ускорений тела, а также определяется полное ускорение его центра масс, выбираемого в качестве полюса. Задача актуальна ввиду ее практического использования в машиностроении, строительстве и прочих производственных областях.

Рис. 1. Расположение акселерометров и радиус-векторов расстояний до них

Имеется некоторое количество акселерометров, измеряющих вектор линейных ускорений в трех произвольных точках движущегося твердого тела, см. рис. 1. Требуется определить мгновенную угловую скорость, мгновенное угловое ускорение и ускорение центра масс этого тела. Положение центра масс тела известно. Разобьем задачу на три этапа. Первый заключается в определении мгновенного углового ускорения тела. Второй - в определении его мгновенной угловой скорости. Третий - в нахождении полного ускорения его центра масс.

Определение углового ускорения тела

Введем неподвижную и подвижную систему отсчета . Неподвижная система отсчета обычно связана с Землей. Начало подвижной системы отсчета совпадает с мгновенным центром ускорений (МЦУ) - Q. Примем за полюс Полюсом называют точку, относительно которой происходит вращение во вращательном движении. точку [1], ускорение которой известно. Тогда ускорение точки можно записать в виде:

(1)

В выражении (1) использовано свойство двойного векторного произведения и введены следующие обозначения:

, i=2,3 - линейные ускорения точек и ;

, i=2,3 - радиус-векторы точек и от МЦУ;

- мгновенное Мгновенное значение - в данный момент времени. угловое ускорение и мгновенная угловая скорость тела; * - знак транспонирования.

Поскольку ускорения точек и известны, то система уравнений (1) позволяет определить необходимые нам параметры. Домножим выражение (1) скалярно справа на вектор , получим модуль мгновенной угловой скорости:

(2)

Домножив выражение (1) векторно слева на , получим модуль мгновенного углового ускорения

(3)

При использовании формул (2) и (3) знать положения МЦУ не нужно, см. рис. 1. Пусть точка - начало отсчета подвижной системы координат. Тогда в силу разности радиус-векторов имеем:

.

Нахождение угловой скорости тела

Начало подвижной системы отсчета выберем в одной из точек установки акселерометра. Запишем ускорения точек и , выбирая в качестве полюса точку :

(4)

Векторные уравнения (4) в проекциях на оси подвижной системы координат дают шесть уравнений относительно шести неизвестных. Распишем их более подробно:

(5)

Система уравнений (5) в векторно-матричной форме имеет вид

, (6)

где введены следующие обозначения:

В вычислительном отношении при решении системы уравнений (6) есть ряд «подводных камней». Например, значение может получиться отрицательным. В этом случае для определения значений угловых скоростей лучше использовать уравнения:

где - полученные числовые значения.

Другая «неприятность» системы (6) это сингулярность матрицы А, что не позволяет обратить ее [2]. В этом случае можно использовать метод сопряженных градиентов [3-8] или регуляризацию матрицы А [9]. Тем более, что матрица А подвергается обращению один раз для данного расположения акселерометров ввиду абсолютной твердости тела.

Итак, в результате решения системы уравнений (6), получаем значения угловых скоростей тела.

Определение ускорения центра масс тела

Поскольку положение центра масс тела известно, тогда выбрав в качестве полюса любую из точек установки акселерометров, получим:

. (7)

При известных значения мгновенной угловой скорости и мгновенного углового ускорения, вычисление выражения (7) не представляет никаких трудностей.

Авторы статьи решали поставленную задачу независимо друг от друга. Поэтому далее представлен другой алгоритм ее решения. Предложенный ниже алгоритм лишен указанных недостатков.

Альтернативный алгоритм решения задачи

Пусть в каждый момент времени известны линейные суммарные ускорения a_i в i-тых точках твердого тела [9]:

Необходимо определить вектор мгновенных угловых е ускорений и скоростей щ твердого тела, а также ускорение центра его масс a_С.

Проекции ускорений i-той точки по оси декартовой системы координат имеют вид [10]:

Учитывая, что:

Для определения 6 компонент ускорений и 3 компонент угловой скорости определенной точки твердого тела достаточно показаний акселерометров с трех точек (i = 1, 2, 3), таким образом, получим систему (8) из 9 нелинейных уравнений относительно 9 неизвестных:

;

; (8)

Решение нелинейной системы (8) вызывает определенные трудности [11-13] и требует существенных затрат машинного времени. Для линеаризации решения сделаем замену переменных (9):

(9)

После подстановки получим систему линейных уравнений:

(10)

Для трех точек твердого тела получим систему 9 уравнений с 12 неизвестными. Недостающие 3 уравнения получим, разместив 4-ый акселерометр в любой другой точке твердого тела и сняв показания с него. В итоге получим систему из 12-ти линейных уравнений, решив которую определим вектор неизвестных:

.

Фактические значения угловых скоростей щ и ускорений е получим из выражений замены переменных (9).

В матричном виде систем уравнений из 12-ти линейных уравнений выглядит следующим образом:

(11)

где введены следующие обозначения:

M - матрица коэффициентов перед неизвестными, определитель которой не равен нулю;

- вектор линейных ускорений четырех точек твердого тела, полученные из показаний акселерометров.

Для решения поставленной задачи использованы два различных подхода. Оба метода привели к одинаковым результатам при решении тестовых задач. Таким образом, удалось значительно сократить расчётное время каждой итерации по определению ускорений заданной точки твёрдого тела путём линеаризации разрешающей системы уравнений.

Литература

ускорение акселерометр угловой скорость

1. Голубев Ю.Ф., Основы теоретической механики. - М.: Изд-во МГУ, 2000. - 719 с.

2. Гантмахер Ф.Р., Теория матриц. - М.: Физматлит, 2010. - 560 с.

3. Акулич И.Л., Математическое программирование в примерах и задачах. - М.: Высш. шк., 1986. - 319 c.

4. Гилл Ф., Мюррей У., Райт М., Практическая оптимизация. Пер. с англ. - М.: Мир, 1985. - 509 c.

5. Коршунов Ю.М., Коршунов Ю.М., Математические основы кибернетики. - М.: Энергоатомиздат, 1987. - 496 c.

6. Максимов Ю.А., Филлиповская Е.А., Алгоритмы решения задач нелинейного программирования. - М.: МИФИ, 1982. - 52 c.

7. Корн Г., Корн Т., Справочник по математике для научных работников и инженеров. - М.: Наука, 1970. - 832 c.

8. Henk A. van der Vorst. Iterative Krylov Methods for Large Linear System. - Cambridge University Press, 2003. - 221 p.

9. Тихонов А.Н., Арсенин В.Я., Методы решения некорректных задач. - М.: Наука, 1986. - 288 с.

10. Лачуга Ю.Ф., Теоретическая механика. - М.: Изд-во КолосС, 2004. - 575 с.

11. Herbert S. Wilf. Algorithms and Complexity. A.K. Peters, second edition, 2002. - 135 p.

12. Косенко, Е.Е., Бурцева, О.А., Нефедов, В.В., Косенко, В.В., Черпаков, А.В. Моделирование напряженного состояния арматурных стержней, применяемых при производстве преднапряженных железобетонных конструкций // Инженерный вестник Дона, 2011, №4 URL: ivdon.ru/ru/magazine/archive/n4y2011/549.

13. Косенко Е.Е., Косенко В.В., Черпаков А.В. Исследование колебаний полнотелой стержневой модели кантилевера с дефектом // Инженерный вестник Дона, 2013, № 4 URL: ivdon.ru/ru/magazine/archive/n4y2013/2153.

Размещено на Allbest.ru