DOI: 10.18698/0236-3941-2016-2-93-112
УДК 531.36
ДВИЖЕНИЕ АБСОЛЮТНО ЖЕСТКОЙ ДВУХОСНОЙ ТЕЛЕЖКИ С УПРУГИМ ТЕЛОМ ПО ДВУХРЕЛЬСОВОМУ КРИВОЛИНЕЙНОМУ ПОЛОТНУ
В.А. Гнездилов, С.В. Русских
МАИ (НИУ), Москва, Российская Федерация e-mail: gnezdilov_07@mail.ru; sergey.russkih@rambler.ru
Рассмотрено нестационарное движение упругого тела, расположенного на абсолютно жесткой тележке, скользящей на упругих роликах по произвольному криволинейному полотну с двумя эквидистантными направляющими рельсами. Полотно задано уравнениями пространственного изгиба и кручения его средней линии. Контакт тележки с рельсами осуществляется через нелинейно-упругие ролики задней ведущей оси, жестко связанной с тележкой, и передней оси, соединенной с ней сферическим шарниром. Разработана новая нелинейная математическая модель для расчета переносного и относительного движений тележки с упругим телом по полотну с определением контактных реакций и перегрузок. Проведено геометрическое моделирование кинематики движения двухосной тележки, в котором считается расстояние между осями, малыми по сравнению с радиусом кривизны средней линии полотна. Получены соотношения, позволяющие дополнительно учесть кривизну средней линии между осями тележки. Динамика относительного движения упругого тела на тележке описана конечными деформациями в квадратичном приближении. Уравнения движения системы составлены с учетом нелинейной упругости и демпфирования роликов, на которых тележка безотрывно скользит по полотну с заданной скоростью или под действием заданной тяги.
Ключевые слова: упругое тело, тележка, движение по криволинейной направляющей, нелинейно-упругий ролик, математическая модель, реакция, перегрузка.
THE MOTION OF AN ABSOLUTELY RIGID BIAXIAL TROLLEY
WITH AN ELASTIC BODY ALONG A DOUBLE RAIL CURVILINEAR
TRACK
V.A. Gnezdilov, S.V. Russkikh
Moscow Aviation Institute (National Research University), Moscow, Russian Federation
e-mail: gnezdilov_07@mail.ru; sergey.russkih@rambler.ru
The research deals with the unsteady motion of an elastic body located on an absolutely rigid trolley moving on the elastic rollers along an arbitrary curvilinear track with two equidistant guide rails. The track is defined by the equations of bending and torsion of its midline. The contact between the trolley and the rails is carried out through the nonlinear elastic rollers ofrear drive axle, rigidly connected with the trolley, and those of the front axle, coupled with it in the middle by the spherical hinge. We developed a new nonlinear mathematical model for calculating the relative and portable motions of a trolley with an elastic body along the surface and defined contact reactions and overloads. Initially, we conducted a geometric modeling of kinematic motion of a biaxial trolley, where the distance between the axes is considered to be small in comparison with the curvature radius of the trolley centerline. Next, we obtained the relations allowing us to further consider the curvature of the midline between the axes of the trolley. Then, we described
dynamics of the relative motion of an elastic body on the trolley by finite deformations in the quadratic approximation. Finally, we worked out equations of system motion taking into account the nonlinear elasticity and damping of the rollers, on which the trolley slides along the surface at a predetermined speed or under a given traction.
Keywords: elastic body, trolley, motion along a curvilinear track, nonlinear elastic rollers, a mathematical model, reaction, overloads.
Во многих областях техники и жизнедеятельности человека таких, как транспортные средства, перевозящие деформируемые (упругие, жидкие, сыпучие) материалы и людей, управляемые летательные и космические аппараты, движущиеся по заданной программной траектории, быстродействующие роботы-манипуляторы, возникают задачи механики, связанные с движением (скольжением или качением) твердых и деформируемых тел по криволинейным поверхностям и направляющим. К таким задачам относятся: пуск беспилотных самолетов и реактивных снарядов по направляющей балке; движение автомобилей и мотоциклов по неровной дороге; движение железнодорожного транспорта по длиннопролетным мостам, совершающим колебания при сильных порывах ветра и землетрясении; взлет самолета с палубы авианосца с трамплином; взлет и посадка самолета на тонком льду; сложное пространственное движение тележек с пассажирами по изогнутым и закрученным рельсовым направляющим, расположенным на эстакадах ("катальные горы"), динамика некоторых спортивных дисциплин и цирковых аттракционов и пр. [1-5].
Контактное взаимодействие между движущимся телом и поверхностью в общем случае осуществляется в нескольких точках посредством податливых колесных шин или роликов с амортизацией и описывается нелинейными упруговязкими двухсторонними, односторонними или устранимыми связями. При решении таких задач ставятся различные требования, основными из которых являются функциональная безопасность, прочностная надежность, физиологическая переносимость перегрузок и вибраций, комфорт экипажа и пассажиров. Некоторые из этих задач могут быть сформулированы и решены в линейной постановке, а некоторые — в геометрически и/или физически нелинейной постановке. При этом в некоторых случаях для упрощения задачи тело, скользящее или катящееся по поверхности, может считаться абсолютно твердым или заменяться сосредоточенной массой, а поверхность, по которой движется тело, считается недефор-мируемой и слабоискривленной.
Некоторые задачи в настоящее время решены в той или иной постановке. Тем не менее, имеются нерешенные задачи и задачи, требующие решения в уточненной нелинейной постановке, в частности, при движении абсолютно твердых и упругих тел с большими углами поворота по сильно искривленным поверхностям (направляющим кривым). Цель настоящей работы — разработка общей математической
модели с обоснованиями для расчета движения упругого тела на абсолютно жесткой тележке, скользящей без отрыва на упруговязких роликах по произвольному двухрельсовому изогнутому и закрученному пространственному полотну. Настоящая статья является продолжением исследований плоского движения твердого тела на упругих колесах (роликах) по криволинейной направляющей [6].
Постановка задачи. Кинематические соотношения. Рассмотрим пространственную изогнутую и закрученную полосу с прямолинейными поперечными сечениями, с двумя направляющими линиями (рельсами) и постоянным расстоянием между ними. Геометрия этой полосы характеризуется параметрами изогнутой в двух плоскостях и нерастяжимой средней линии Б.0(в) = [Хм(в) Х02(в) Х03(в)]т и углом закручивания в\ (в) поперечного сечения относительно касательной к этой линии, где 5 — длина дуги, отсчитываемая вдоль этой линии; Х01, Х02, Х03 — координаты точки во средней линии в неподвижной системе 0*Х\Х2Х3, ось Х2 которой направлена вверх противоположно вектору ускорения свободного падения д0 (рис. 1).
По полотну на упругих роликах, размеры которых малы, безотрывно движется абсолютно жесткая двухосная тележка с расположенным на ней деформируемым телом. Связи между полотном и тележкой считаются двухсторонними, безотрывное качение роликов осуществляется двумя роликами, расположенными в вертикальной плоскости по разную сторону от рельса, однако в рассматриваемой модели движения два ролика заменяются одним с единственной точкой контакта. Трение между полотном и роликами отсутствует, т.е. в рассматриваемой модели полотно считается абсолютно гладкой поверхностью.
Рис. 1. Геометрическая постановка задачи
Задняя ось жестко связана с тележкой, а передняя ось в середине связана с ней сферическим шарниром. При этом в расчетной модели два одинаковых ролика передней оси заменяются одним роликом, эквивалентным им по податливостям и движущимся по средней линии. Таким образом, по существу рассматривается трехроликовая тележка.
Подвижная система координат Ох1х2х3 с началом О на средней линии связана с тележкой: ее плоскость Ох\х3 проходит через контактные точки недеформированных роликов передней оси (точка А) и задней оси (точки В и С); ось Ох1 проходит через точку А (ОА = I), а ось Ох3 — через точки В и С (ОВ = ОС = а) (см. рис. 1). Точки В и С движутся по кривым, являющимся эквидистантами к И0(в). Сначала для упрощения задачи рассмотрим случай, когда кривизной И0(в) в пределах расстояния между осями можно пренебречь. Тогда подвижную систему координат Ох1Х2Х3 можно заменить на систему Ох1х2х3 (см. рис. 1), в которой Ох1 является касательной, Ох2 — нормалью к кривой К0(в) в точке 8 = во, а ось Ох3 дополняет их до правой системы координат. При этом можно считать в а — во ~ I. Далее для учета кривизны И0(в) на длине вА — во систему Ох1 х2х3 довернем на малые углы для совмещения ее с системой Ох1Х2Х3.
Скорость тележки V0(t) в точке О при условии, что она создается вращающимися роликами ведущей задней оси, направлена по касательной в точке О к кривой И0(в); при этом ^0(Ь)| = У0(Ь) = ¿в/ёЬ. Считается, что задана скорость У0(Ь) или приложенная в точке О по касательной к кривой И0(в) сила тяги Р0(Ь). На движущуюся тележку со стороны полотна действуют реакции (см. рис. 1): Р1 и Р2 — реакции на заднюю ось, параллельные оси Ох2, приложенные соответственно в точках В и С; реакция Р3, приложенная в точке А на переднюю ось параллельно Ох2; Р4 — реакция на заднюю ось вдоль Ох3 и приложенная в точке О; реакция Р5 на переднюю ось вдоль Ох3 и приложенная в точке А.
Ролики из физически нелинейного вязкоупругого материала за счет контактного смятия допускают малые перемещения в направлении нормалей к плоскостям Ох1х2 и Ох1х3, за счет этого тележка как твердое тело перемещается в направлении осей Ох3 и Ох2 и поворачивается относительно Ох1, Ох2, Ох3.
Установленное на тележке деформируемое тело считается упругим. Точка М с координатами х1, х2, х3 совершает относительные перемещения в направлениях этих координат ик(х1, х2, х3, Ь), к = 1, 2, 3, в состав которых входят перемещения тележки как твердого тела за счет деформаций роликов и перемещения, обусловленные деформациями тела. Конечные деформации (относительные удлинения е11, е22, е33 и углы сдвига 712, 723, 731) определяются через проекции перемещения
на оси подвижной системы координат «1(^1, х2, ж3, ¿), м2(ж1, х2, х3, ¿), м3(ж1, х2, х3, ¿) через квадратично нелинейные соотношения [7, 8]:
en =
dui dxi
1
+ 2
du2
dx1
+
du3 dx1
ди1 ди2 ди3 ди3
Т12 гч I гч I гч гч ; . . . ,
дх2 дх1 дх1 дх2
которые справедливы при ди1/дх1 ^ 1, (дм1/дж2)2 ^ 1, ...; (1 ^ 2 ^ 3 ^ 1).
Введем векторные обозначения: г = [х1 х2 ж3]т, и = [и1 и2 и3]т, г = г + и, И = И0 + Лтг, где Л — матрица перехода от системы О^Х1Х2Х3 к системе Ож1ж2х3 (Л-1 = Лт). Вращение системы Ож1ж2ж3 относительно 0*Х1 Х2Х3 будем характеризовать заданными самолетными углами вь в2 и в3, представляющими собой соответственно углы крена, рыскания и тангажа, с последовательностью поворотов в2 ^ в3 ^ вь Вектор углов поворота будем обозначать через в = [01 в2 в3]т. Матрица Л записывается в виде [1, 7, 8]
Л =
c2 Сз
SiS2 - C1C2S3
S3 С1С3
-S2C3 S1C2 + C1S2S3
c1S2 + S1c2S3 — s1c3 C1C2 — S1S2S3
где sk = sin 9k, ck = cos 9k (k = 1, 2, 3). Поскольку ось Ox1 подвижной системы координат в рассматриваемом приближении является касательной к Ro(s) в точке sO, углы поворота 92 и 93 выражаются через проекции вектора R0(s) с помощью дифференциальных соотношений [9-11]:
dX01 dX02 dX03
= С2 С3; "— = s3; :- = -s2C3,
ds
ds
ds
а угол 01(в) задается независимо.
Вектор мгновенной угловой скорости ш = [ш1 ш2 ш3]т подвижной системы координат Ож1ж2ж3 выражается через вектор углов поворота в как [1, 7, 8]:
ш = Л0; A =
1 0
0 С1С3 0 -Я1С3 С1
При решении геометрически нелинейных задач по методу Ритца или по методу конечных элементов перемещение упругого тела в общем случае может быть представлено в виде [7, 8] и = ^^ +
г
I"! ^^^, где — обобщенные координаты; г, ] = 1, 2,..., п;
2
2
^¿(х1, х2, х3), ф ^(х1, х2, х3) — заданные вектор-функции, которые выбираются в классе возможных перемещений тела согласно принятой расчетной модели (стержень, пластина и др.) с учетом определенных допущений. Функции ф ^ выражаются через ^ и ^ и являются симметричными, т.е. ф^ = ф^. Обобщенные координаты при г = 1, 2,..., п — упругие деформации тела.
Скорость и ускорение произвольной точки М упругого тела, связанного с подвижной системой координат Ох1 х2х3 записываются в виде
V = Vо + ш х г + и; а0 = V/0 + ш х Vо; а = а0 + ш х г + ш х (ш х г) + 2ш х и + и,
где а0 — вектор ускорения точки О. Вектор угловой скорости ш и его производная по времени ш находятся из следующих соотношений:
ёв
ёв
- = A f |Vo|
ds
dt
ds ds ds2
V0 + A^
V 0
Для удобства записи векторное произведение заменяется матричным:
axb = a b;
a =
а1 " bi " v ; a = 0 -аз a2
02 ; b = b2 a3 0 —а1
a3 b3 —a2 ai 0
Если движущееся по направляющему полотну тело, связанное с подвижной системой координат Ох1х2х3, является абсолютно жестким, а ролики — недеформируемыми, то задача в заданный момент времени Ь является задачей кинематики. Тогда после определения ускорения а(х1, х2, х3, Ь) из уравнений динамического равновесия твердого тела с учетом силы тяжести находят главные вектор Р и момент М реакций и затем — сами реакции. После этого можно выполнить расчеты на прочность конструкции тележки и всего сооружения вместе с полотном и опорными устройствами.
Если учитывается относительное движение, представленное вектором перемещений и(х1, х2, х3,Ь), обусловленное упругостью роликов и амортизации, а также упругим деформированием присоединенных к тележке масс, то задача существенно усложняется. При заданных векторных функциях И0(в), в(в), скорости V0(t) и найденной векторной функции ш(Ь) задача сводится к задаче динамики относительного движения для обобщенных координат, характеризующих неизвестную векторную функцию и(х1, х2, х3, Ь).
Нелинейные уравнения относительного движения. Разрешающие уравнения динамики тела, движущегося по произвольной поверхности, можно получить по аналогии с уравнениями возмущенного движения упругих летательных аппаратов по заданной программной траектории [1, 7, 8, 12].
Получим нелинейные уравнения динамики рассматриваемой системы для наиболее общего случая, учитывая относительное движение, обусловленное упруговязким деформированием роликов и конечными деформациями присоединенного упругого тела.
Уравнения движения системы для заданных векторов У0, ш и обобщенных координат дг составляются на основе принципа Даламбера-Лагранжа как для свободной системы с освобожденными связями, включая неизвестные реакции в число внешних сил:
5П = 5 Ар + 5Аин, (1)
где П[и] — потенциальная энергия деформации системы; 5АР — вариация работы гравитационных сил, силы тяги Р0 и реакций Р1... Р5 со стороны полотна, приложенных в точках с координатами , ,
; 5Аин — вариация работы инерционных сил. Выражения для 5АР и 5Аин имеют следующий вид:
5Ар = 5и^ + )Р; (2)
V
5Аи-=-/5иТа (3)
V
где 5и* — вариация перемещения любой точки тела с учетом вариаций смещения 5и0 = А5И0 и малого поворота 5в подвижной системы координат
5и* = 5и0 + 5в х г + 5и; (4)
g — вектор массовых сил тяготения (будет определен позднее); ^ш = р^ — масса элемента тела, сохраняющая свою величину в процессе его деформирования; р(ж , ж2,ж3) — плотность тела; ¿V — элемент объема тела. Из выражения для вектора и с учетом ф^ = получим
5и = ^2и = ^2 ^;
г _ г (5)
и = ^2 ^ + ^2 ^ ; ^ = ^ + ^2 ъ .
г г ] ]
Поскольку вариации 5и0, 5в и 5дг произвольны и независимы, из принципа Даламбера-Лагранжа (1) с учетом (2)-(5) следуют уравнения движения:
J а ^ш = Р; J га ^ш = М;
Г т , дП ^ (6)
^га ¿ш + ^ = Уг, г = 1, 2, ..., п, У ддг
V
где Р = Рё + Рд; М = Мё + Мл; Рд и Мд — векторы равнодействующих сил и моментов от реакций; Ql = Векторы Рё (Рё = — Рё ), Мё и Qgll вычисляются как
Pg = ^ g ^ш; Mg = ^ г g ёш; Qg¡l = ^ g ёт. (7)
V V V
После преобразований система уравнений (6) с добавлением сил вязкого демпфирования по обобщенным координатам ql записывается в виде (г = 1, 2, ..., п):
г
ша0 — Ь ш + X Ь д —
V V 3
(ш Ь0 +2 X Ь3 ъ)ш + X X Ьзкъ^ = Р;
3 3 к
v
- L Tao + ÍqqC^ + X S0j <ij +
j
+ (cv Íoo + 2 X Íoj ójj )ш + X X Soj ójj qk = M;
j j k
Lt ao + S0i" + X] mij- (^TÍtí + 2 X Sj^j)W + jj
+ XX m-ijk q?j 4h + X dij 4j +1" = Qi,
j k j
где
(8)
m = J dm; L0 = J rdm; L% = J ^%dm;
V V V
/ryv
^¿j dm; Ioo = - rr dm;
V V
/у y fy_
r ^ dm; So¿ = r ^%dm; (9)
V V
/у _ r y
r -фгзdm; S¿j = ip% ^dm;
V V
m%j = J фтф3 dm; m%ijfe = J фтг-фзк dm; i, j, k = 1, 2,..., n;
VV
d%j = dj% — коэффициенты вязкого демпфирования, которые обычно определяются на основании экспериментальных данных. Первые два уравнения (8) используются для определения реакций P1... P5 и силы тяги P0, а третье уравнение представляет собой уравнение относительного движения. В выражениях (9) инерционные характеристики,
отмеченные сверху волной, вычисляются для деформированного тела и зависят от обобщенных координат в силу того, что от них зависят функции г = г + и и фг.
Потенциальная энергия деформации упругого тела записывается в виде [13]
П = 1 У (е - вт)= 1 ^ (е - ет)ТС(е - е^У,
V V
где е = [еп £22 £33 712 723 7з1]т, & = [сп ^22 ^зз Т12 Г23 Т31]Т — векторы, составленные из компонент деформаций и напряжений; вектор ет представляет начальные нестесненные температурные или технологические деформации; С — симметричная матрица коэффициентов упругости, входящих в уравнения закона Гука.
Вектор массовых сил тяготения вычисляется следующим образом:
g = -gоV; V = Л [0 1 0]Т , (10)
где g0 — ускорение свободного падения (модуль вектора g0); V — единичный вектор, направленный вдоль оси Х2 неподвижной системы координат О^Х1Х2Х3 и записанный в проекциях на оси подвижной системы. Гравитационные силы и моменты (7) с учетом (10) описываются следующими выражениями:
v
т
р = ; Щ = -£0 £ V; = -§0^гV.
В результате получается замкнутая система нелинейных дифференциальных уравнений: шесть уравнений второго порядка для шести компонент векторов Рд, Мд и п уравнений для обобщенных координат дьд2 ...Яп.
Если движущееся по направляющему полотну тело и ролики считаются абсолютно жесткими (я1 = д2 = ... = дп = 0), то уравнения (8) упрощаются
v v v v т v
ша0 — Ь0ш — ш Ь0ш = Р; — Ь0Та0 + 100шШ + ш 100ш = М. (11)
При известных векторах У0, ш, полученных при решении кинематической задачи, из уравнений (11) определяются главные вектор Р и момент М сил реакции и затем — сами реакции.
Учет влияния кривизны полотна между осями. В задаче, сформулированной ранее, предполагалось, что расстояние между осями тележки мало по сравнению с кривизной центральной линии, ролики передней и задней осей в данный момент времени находятся на прямой, касательной к кривой И0 (в) в точке 8 = во .В случае двухосной тележки со свободно поворачивающейся относительно центра
передней осью при движении по сильно искривленному пространственному полотну ролики задней и передней осей будут находиться в точках в = во и в = в а кривой И0(в) с различными углами наклона касательной. Локальная кривизна полотна между осями окажет некоторое влияние на угловое положение тележки в данный момент времени и на кинематику и динамику ее движения.
Учтем влияние локальной кривизны направляющей И0(в) между осями тележки. Введем систему координат Ож1Х2Ж3, ось ОХ1 которой проходит через точку А, (см. рис. 1). Матрицу перехода от системы координат Ож1ж2ж3 к ОХ1Х2Х3 будем обозначать через А*, а матрицу перехода от системы координат О*Х1Х2Х3 к АХ1Х2Х3 — через Л. Очевидно, что выполняются соотношения Л = Л*Л, Л* = ЛЛт.
Естественная координата 8а точки А находится из следующего
г
трансцендентного уравнения при известной координате во =J|(Ь) | ёЬ:
/2= (xoi(sa)-xOI(sO))2 + (xO2(sA)-xO2(sO))2 + (Xc3(sA)-Xo3(so))2 .
Если кривизна кривой R0(s) небольшая, то можно приближенно полагать, что sA ~ + l.
Введем вектор AIR, который характеризует разность соответствующих координат точек O и A:
xOI(sA) - Xol(so) " AX "
АН = X02 ( s A ) - Xo2(so) = AXX2
_ xo3(sa) - Xo3(so) _ AX3
Углы поворота системы координат Axix2x3 будем обозначать как
01 = 01 + A01; 02 = 02 + Д02; 0з = 0з + А0з, (12)
где А01, А02 и А03 — малые углы (sinA0i w A0¿, cosA0i w 1, i = 1, 2, 3). При этом с учетом (12) будем иметь
Si = sin (0i + A0i) w Si + A0iCi; 0 = cos (0i + A0i) w Ci - A0iSi.
Заметим, что если угол A02 немалый, то необходимо, чтобы задняя ось соединялась с тележкой шарнирно, как и передняя ось. Матрица перехода Л(01,02,03) определяется так же, как матрица Л(01,02,03) с заменой si ^ 0i, ci ^ 0i. С точностью до линейных членов c A01, A02 и A03 получим
- 4 д Л . п д Л . п д Л А Л
Л = Л + eel ^ + a"'2 +
(13)
дЛ
дЛ
0 0 0
cls2 + slc2s3 — slc3 С1С2 — sls2s3 — sls2 + clc2s3 —clc3 —slc2 — Cis2s3
-S2C3 0 -C2C3
SlC2 + ClS2S3 0 —SlS2 + ClC2 83 clc2 — sl s2s3 0 —cls2 — slc2s3
дЛ
— c2s3 c3 s2s3
—ClC2C3 —Cl83 cls2c3 SlC2C3 SlS3 —SlS2C3
" дхх1 " ■ l '
ДХ2 = 0
ДХ3 0
Уравнение для вектора г в точке А гА = ЛДИ при х\а = I, х2А = 0, Х3А = 0 с учетом (13) запишем в виде
дЛ дЛ дЛ
Л + eil ^ + ^ + й- ^
Это уравнение эквивалентно трем скалярным уравнениям, из которых определяются углы Д$ь Д$2 и Д$3:
lio + /ПД$1 + /12 Д $2 + /^Д$3 = /;
/20 + /21Д $1 + /22 Д $2 + 123Д$3 = 0; (14)
1зо + /31Д$1 + /32 Д $2 + /33Д$3 = 0.
Здесь
/10 = С2С3ДХ1 + S3ДХг_- 82в3ДхХ3;
/20 = (S1S2 - C1C2S3) Дх^1 + С1С3ДХ2 + (S1C2 + С1S2S3) ДХ3; /30 = (C1S2 + S1C2S3) ДХ - S1С3Д+ (С1С2 - S1S2S3) ДХ3; /11 = 0;
/12 = -^2С3ДХ1 - С2С3ДХ3;
/13 = C2S3 ДХ1 + С3ДХ2 + «2в3ДхХ3;
/21 = (C1S2 + S1C2S3) ДХ1 - S1С3Д+ (С1С2 - S1S2S3) ДХ3;
/22 = (S1C2 + C1S2S3) Дх^1 -JS1S2 - C1C2S3) ДХ3;
/23 = -С1С2С3ДХ1 - C1S3ДХ2 + C1S2C3ДХ3;
/31 = (-S1S2 + C1C2S3) ДХ - С1С3ДХ2 - (S1C2_ + C1S2S3) ДХ3
/32 = (С1С2 - S1S2S3) ДХ1 - (С1S2 + S1C2S3) ДХ3;
/33 = S1C2C3+ S1S3ДХ^2 - 8182в3ДхХ3.
В результате с учетом найденных значений Д$1, Д$2 и Д$3 определяется переходная матрица Л для подвижной системы координат Ox1X2X3, которая по сравнению с матрицей Л учитывает влияние кривизны полотна направляющей.
Исходные линейные и угловые скорости движения V0 и ш, а также их первые производные по времени V0 и Ш, задаются в системе
координат ОЖ1Ж2 Х3, а затем пересчитываются в системе координат
ОХ1Х2Х3:
V0 = A*V0; ш* = Л*ш;
V 0 = f Vo + aV о;
* dA* ш = —т— ш + Л ш;
dA* dA*
dt
dt
ds
|Vo(t)|.
Учет упругости и демпфирования роликов. За счет упругости и демпфирования роликов тележка будет совершать дополнительные движения как твердое тело с пятью степенями свободы относительно подвижной системы координат Ож1ж2ж3 или ОХ1Х2Х3. В качестве обобщенных координат, представляющих эти движения, будем рассматривать: перемещения д0з = и02, 9о3 = и03 полюса (точки О) вдоль координатных осей х2, х3 (или Х2, Х3) и углы поворота тележки д04 = Дб1, д0б = Дб2, д0б = Дб3 в квадратичном приближении (втДбг « Дбг, совДбг « 1 — Дб2/2, г = 1, 2, 3) относительно осей Х1, Х2, Х3 (или относительно осей Х1, Х2, Х3) (рис. 2).
В результате точки О, А, В и С, которые по-прежнему связаны с траекторией движения, смещаются в точки О', А', В' и С', которые в данном случае характеризуют положение тележки после введения упругостей роликов, т.е. О'А' = I и В 'С' = 2а (смещения точек В и С по оси х3 или Х3 равны).
В этом случае вектор дополнительных перемещений твердого тела, обусловленных упругими деформациями роликов, будет иметь вид
6 1 6 6 й = ^ ^ ^ + 2 ^ ^ ^ % Ф00,
г=2 г=4 _у=4
Рис. 2. Схема учета упругости и демпфирования роликов
0 0 0 X3
^02 — 1 ; ^0з— 0 ; ^04— -X3 ; ^0Б — 0
0 i X2 —X1
^06 =
-Х2
Xi 0
Ф0404 =
0
-X3 -X2
ф0д0к —
X3 —X1 0 —X1
ф0л0б — 0 ; ф0в0в— 0 ; ф0в0б— 0 ; ф0б0б— — X2
0 —X3 X2 0
X2 0 0
При использовании подвижной системы координат ОХ1Х2Х3 вместо 0ж1ж2ж3 необходимо выполнить замену ж» ^ ж», г = 1, 2,3.
В общем случае при физически нелинейных деформациях, включая
возможные зазоры, потенциальная энергия роликов тележки Пр как составная часть потенциальной энергии П записывается в обобщенных координатах как Пр = Пр (д02, , ?о4, ?оБ, 5о6). При малых упругих деформациях роликов (в пределах закона Гука) и малых углах поворота тела, обусловленных этими деформациями, будем иметь
1
1
nP = q0iq0j или nP = 2 qTKpЧс где Kp = [kp] пРи
i=2 j=2
i, j = 2... 6 — матрица коэффициентов жесткости, q0 =
= [q02 q03 q04 q05 q06]T-
В общем случае вязкоупругих роликов вариация работы сил демпфирования записывается в следующем виде: ÍAp = — £qTDpq0, где Dp = [dp] при i, j = 2 ... 6 — матрица коэффициентов демпфирования.
Вектор перемещений тележки ü, обусловленный податливостью роликов, складывается с вектором относительных упругих перемещений присоединенных масс, характеризуемых обобщенными координатами qi, q2 ... qn, и суммарный вектор обозначается через ü. В результате с учетом податливостей роликов вектор u будет выражаться через обобщенные координаты q02, q03... q06, q1, q2 ... qn.
Пример расчета. Приведем пример расчета параметров движения абсолютно жесткой двухосной тележки с абсолютно жестким телом по пространственному полотну, образованному из центральной винтовой линии R0(s). Будем считать, что оси подвижной системы координат Ox1x2x3 без учета локальной кривизны R0(s) совпадают с осями трехгранника Френе: ось Ox1 — касательная t, ось Ox2 — главная нормаль n, ось Ox3 — бинормаль b.
Уравнение винтовой линии R0(s) в неподвижной системе координат O*X1X2X3 имеет вид [10]: X01(s) = R cos <£, X02(s) = R sin X03(s) = A^, где ^ — угол закрутки линии относительно оси O*X3;
Рис. 3. Общий вид траектории движения
R — радиус винтовой линии (радиус цилиндрической образующей); Л — шаг подъема винтовой линии за угол поворота (, равный 1 рад (2пЛ — подъем за один виток) (рис. 3). В любой точке винтовой линии касательная t, нормаль n и бинормаль b имеют следующие координаты (здесь и далее R = R¡л/R2 + Л2, Л =
= Л//R2 + Л2): t = [—Rsin (Rcos ( Л]T; n = [— cos ( — sin (О] T;
b = [Л sin ( — Л cos (R ]T.
Матрица Л, как матрица направляющих косинусов единичных векторов t, n и b (без использования углов 9\, в2 и 93), записывается в следующем виде:
(—R sin ( R cos ( Л — cos ( — sin ( О Лsin ( — !\cos ( R
Функциональная связь между длиной дуги s (она отсчитывается от точки M(R, О, О) при ( = О) и углом закрутки ( выражается следующим образом: s(f) = (\/R2 + Л2. Кривизна k(s) и кручение ж(з) винтовой линии в любой ее точке постоянны и равны [10]
Поскольку подвижная система координат Ox\x2x3 связана с трехгранником Френе, то вектор мгновенной угловой скорости ш может быть определен через векторное соотношение ш = + кЬ. Поскольку матрица Л и вектор ш в данном примере могут быть записаны
без использования углов , #2 и #3, то необходимо разработать иной алгоритм определения матриц Л и Л*, чем приведенный ранее. Матрицу поворота Л* получим как матрицу поворота единичного вектора t = [t1 t2 t3]T (ось Ox1) на угол Ф = Z (Ox1, Ox1) относительно некоторого единичного вектора $ = £2 £3]T, чтобы ось Ox1 совпала с Ox1 и вектор t совпал с единичным направляющим вектором Р = [p1 Р2 P3]T оси Ox1. Координаты вектора р определяются с помощью уравнения прямой Ox1 как прямой, проходящей через две заданные точки O и A, после нормирования:
xcH(sa) — X01(s0) _ X02(sA) — X02(so) P1 = l ; P2 = l ;
xO3(sA) — Xo3(so) P3 = l . Угол поворота Ф найдем из скалярного произведения единичных векторов t и Р:
cos Ф = ¿1P1 + Í2P2 + ¿3P3.
Координаты вектора $, вокруг которого осуществляется поворот, определяются как координаты вектора нормали плоскости Ox1X1 или через векторное произведение $ = t х р:
6 = ¿2P3 — ¿3 P2; £2 = ¿3P1 — ¿1P3; £3 = ¿ 1P2 — ¿2P1.
Матрица поворота Ox1x2x3 ^ Ox1X2Х3, учитывающая локальную кривизну направляющей поверхности, имеет следующий вид:
(cos Ф + (1 — cos Ф) £2 (1 — cos Ф) £^2 + (sin Ф) £3 (1 — cos Ф) £26 — (sin Ф) £3 cos Ф + (1 — cos Ф) £2 (1 — cos Ф) £3£1 + (sin ф) £2 (1 — cos Ф) £3£2 — (sin Ф) £1
(1 — cos Ф) £1£3 — (sin Ф) £2 \ (1 — cosФ) £2£3 + (sinФ) £1 I . cos Ф + (1 — cos Ф) £2 /
Отсюда определяется матрица Л = Л*Л.
Исходные данные для расчета: R = 7м, A = 5/2п = 0,796 м, V0(t) = 3t + 7sin(t/2) м/с, к = 0,141 1/м, ж = 0,016 1/м, R = 0,994, Л = 0,113, а = 1 м, I = 2 м, m = 100 кг — масса тележки. Координаты центра тяжести тела (в вертикальной плоскости симметрии) — r =[10, 5 0]T. Рассматриваемый интервал движения t £ [0, 20] с. Число разбиений рассматриваемого отрезка времени — 5000. Задача решалась в двух постановках: без учета локальной кривизны направляющей (расчетный случай 1 — РС № 1), с учетом локальной кривизны направляющей (расчетный случай 2 — РС № 2). Функция V0(t) была
подобрана таким образом, чтобы на начальном этапе движения скорости были относительно небольшие (0... 10 м/с), а в конце интервала немонотонно возрастали до величин 40... 50 м/с. Это позволяет более детально рассмотреть необходимость учета локальной кривизны. Все расчеты выполнены в программном комплексе PTC MathCad 15 M030.
В табл. 1 приведено сравнение результатов решения задачи по компонентам вектора скорости V = [V V2 V3]T и компонентам ускорения а = [а\ а2 а3] для различных моментов времени. Из табл. 1 следует, что в наибольшей степени учет локальной кривизны влияет на кинематические параметры, записанные в проекции на бинормаль b. Эти параметры могут не только изменить свой знак, но и изменить свой функциональный вид — в РС № 1 функции V3 (t) и а3 (t) являются прямыми горизонтальными линиями (const), а вид этих же функций для РС № 2 показан на рис. 4. Этот факт следует учитывать при определении перегрузок и биомеханическом анализе [14].
Таблица 1
Кинематические параметры Расчетный случай
№ 1 № 2 № 1 № 2 № 1 № 2
при t = 5 c при t = 10 c при t =15
У, м/с 19,119 18,928 23,217 22,985 51,495 50,981
У2, м/с 0,141 -0,239 0,141 -0,321 0,141 -0,879
У3, м/с 8,017-3 0,048 8,017-10-3 0,044 8,017-10-3 -0,058
ах, м/с2 0,176 0,227 3,973 3,998 4,193 4,293
а2, м/с2 2,696 2,639 3,274 3,131 7,263 7,036
а3, м/с2 2,261 •ÍO-3 0,049 2,261 •ÍO-3 -0,112 2,261 •ÍO-3 -0,693
В табл. 2 приведено сравнение результатов решения задачи по реакциям со стороны тележки на кривую Р1-Р5 и Р0. Отрицательные значения реакций Рх-Р3 не означают отрыв тележки от полотна. Это говорит о том, что при учете обхвата рельса двумя роликами в вертикальной плоскости деформироваться начнет нижний ролик, а не верхний. В наибольшей степени учет локальной кривизны влияет на реакции Р4 и Р5, т.е. реакции, векторы которых направлены по бинормали Ь, причем различие проявляется на конечном этапе движения при больших скоростях. Следует отметить, что: 1) вид кривых Р1-Р5 и Р0 качественно для РС № 1 и РС № 2 не изменяется; 2) значения Р4 и Р5 равны; 3) значения Рх и Р2 близки. На рис. 5 приведены графики изменения реакции Р4 (Р5) по времени в РС № 1 и РС № 2.
Рис. 4. Проекции скорости (я) и ускорения (б) на бинормаль
Таблица 2
Реакция тележки, Н Расчетный случай
№ 1 №2 № 1 № 2 № 1 №2
при t = 5 c при t = 10 c при t = 15 c
Pi -134,335 -129,267 223,449 221,546 494,055 469,680
P2 -182,455 -178,614 181,821 186,080 520,960 530,626
Рз 96,693 103,757 571,804 562,907 567,684 541,930
Р4 48,120 49,348 41,628 35,466 -26,905 -60,946
Р5 48,120 49,348 41,628 35,466 -26,905 -60,946
Ро - 826,967 -826,967 -333,068 -310,563 894,662 916,751
Рис. 5. Реакции со стороны тела на направляющую по бинормали:
1 - РС № 1; 2 - РС № 2
Таким образом, учет локальной кривизны направляющей позволяет уточнить значения проекций скорости и ускорения по бинормали, а также влияет на значения реакций со стороны тела в направлении бинормали.
Выводы. Разработана новая общая математическая модель для расчета движения тела на упруговязких роликах по произвольной двухрельсовой изогнутой и закрученной пространственной направляющей с учетом относительного движения упругих присоединенных масс. С учетом разработанной методики расчета динамических нагрузок (реакций движущегося тела) в качестве дальнейшего исследования темы можно построить алгоритм, позволяющий проектировать и моделировать криволинейную направляющую поверхность, которая удовлетворяет определенному набору кинематических, прочностных и биомеханических критериев. Методики и алгоритмы, разработанные в настоящей работе, могут быть использованы для расчета различных механических систем в проектной практике, в частности, для проектирования систем запуска беспилотных летательных аппаратов, экстремальных катальных гор и т.д.
ЛИТЕРАТУРА
1. Горшков А.Г., Морозов В.И., Пономарев А.Т., Шклярчук Ф.Н.Аэрогидроупру-гость конструкций. М.: Физматлит, 2000. 592 с.
2. Гнездилов В.А. Проектирование и изготовление металлоконструкций для сложных механизированных аттракционов // Монтажные и специальные работы в строительстве. 2000. № 6. С. 20-24.
3. Wayne T. Roller Coaster Physics, Wayne, 1998. 154 p.
URL: http://vip.vast.org/BOOK/HOME.HTM (дата обращения 24.04.2015).
4. Кулаков Н.А. Воздействие динамической нагрузки на наземные транспортные средства. Избранные проблемы прочности современного машиностроения. М.: Физматлит, 2008. 204 с.
5. Селифонов В.В. Теория автомобиля. М.: Гринлайт, 2009. 206 с.
6. Русских С.В. Движение твердого тела на двух колесах по плоской кривой // Известия высших учебных заведений. Машиностроение. 2014. № 2 (647). С. 5258.
7. Шклярчук Ф.Н.Нелинейные и линеаризованные уравнения движения упругих космических конструкций // Изв. РАН. МТТ. 1996. № 1. С. 161-175.
8. Гришанина Т.В., Шклярчук Ф.Н.Динамика упругих управляемых конструкций. М.: Изд-во МАИ, 2007. 328 с.
9. Бюшгенс С.С. Дифференциальная геометрия. М.: Комкнига, 2006. 302 с.
10. Рашевский П.К. Курс дифференциальной геометрии. М.: Изд-во ЛКИ, 2013. 432 с.
11. Banchoff T., Lovett S.T. Differential Geometry of Curves and Surfaces. AK Peters, Taylor & Francis, 2010. 352 p.
12. Космодемьянский А.А. Динамика космического полета. М.: Либроком, 2011. 248 с.
13. Работнов Ю.Н.Механика деформируемого твердого тела. М.: Наука, 1988. 712 с.
14. Рабинович Б.А. Безопасность человека при ускорениях (Биомеханический анализ). М.: Книга и бизнес, 2007. 208 с.
REFERENCES
[1] Gorshkov A.G., Morozov V.I., Ponomarev A.T., Shkliarchuk F.N. Aerogidrouprugost' konstruktsiy [Construction aerohydroelasticity]. Moscow, Fizmatlit Publ., 2000. 592 p.
[2] Gnezdilov V.A. Design and fabrication of steel structures for complex mechanized attractions. Montazhnye i spetsial'nye raboty v stroitel'stve [Erecting and special works in construction], 2000, no. 6, pp. 20-24 (in Russ.).
[3] Wayne T. Roller Coaster Physics, Wayne, 1998. 154 p. Available at: http://vip.vast.org/BOOK/HOME.HTM (accessed 24 April 2015).
[4] Kulakov N.A. Vozdeistvie dinamicheskoi nagruzki na nazemnye transportnye sredstva. Izbrannye problemy prochnosti sovremennogo mashinostroeniia [The impact of the dynamic loads on ground vehicles. Selected problems of modern engineering strength]. Moscow, Fizmatlit Publ., 2008. 204 p.
[5] Selifonov V.V. Teoriya avtomobilya [The theory of car]. Moscow, Grinlait Publ., 2009. 206 p.
[6] Russkih S.V. Motion of a rigid body on two wheels on a plane curve. Izv. Vyssh. Uchebn. Zaved., Mashinostr. [Proc. of Higher Educational Institutions. Machine Building], 2014, no. 2 (647), pp. 52-58 (in Russ.).
[7] Shklyarchuk F.N. Non-linear and linearized equations of motion of the elastic space structures. Izv. RAN. MTT [Proceedings of the Russian Academy of Sciences. Mechanics of rigid body], 1996, no. 1, pp. 161-175 (in Russ.).
[8] Grishanina T.V., Shklyarchuk F.N. Dinamika uprugikh upravlyaemykh konstruktsii [Dynamics of elastic controlled constructions]. Moscow, MAI Publ., 2007. 328 p.
[9] Biushgens S.S. Differentsial'naia geometriia [Differential geometry]. Moscow, Komkniga Publ., 2006. 302 p.
[10] Rashevskiy P.K. Kurs differentsial'noi geometrii [A course of differential geometry]. Moscow, LKI Publ., 2013. 432 p.
[11] Banchoff T., Lovett S.T. Differential Geometry of Curves and Surfaces. AK Peters, Taylor & Francis, 2010. 352 p.
[12] Kosmodem'yanskiy A.A. Dinamika kosmicheskogo poleta [The dynamics of space flight]. Moscow, Librokom Publ., 2011. 248 p.
[13] Rabotnov Yu.N. Mekhanika deformiruemogo tverdogo tela [Mechanics of a deformable solid]. Moscow, Nauka Publ., 1988. 712 p.
[14] Rabinovich B.A. Bezopasnost' cheloveka pri uskoreniyakh (Biomekhanicheskiy analiz) [Human safety when accelerating (Biomechanical analysis)]. Moscow, Kniga i biznes Publ., 2007. 208 p.
Статья поступила в редакцию 17.06.2015 Гнездилов Владимир Алексеевич — канд. техн. наук, заведующий кафедрой "Прочность авиационных и ракетно-космических конструкций" Московского авиационного института (национального исследовательского университета) (Российская Федерация, 125993, Москва, Волоколамское ш., д. 4).
Gnezdilov V.A. — Cand. Sci. (Eng.), Head of Strength of Aviation and Aerospace Structures Department, Moscow Aviation Institute (National Research University) (Volokolamskoye shosse 4, Moscow, 125993 Russian Federation).
Русских Сергей Владимирович — канд. физ.-мат. наук, доцент кафедры "Прочность авиационных и ракетно-космических конструкций" Московского авиационного института (национального исследовательского университета) (Российская Федерация, 125993, Москва, Волоколамское ш., д. 4).
Russkikh S.V. — Cand. Sci. (Phys.-Math.), Assoc. Professor of Strength of Aviation and Aerospace Structures Department, Moscow Aviation Institute (National Research University) (Volokolamskoye shosse 4, Moscow, 125993 Russian Federation).
Просьба ссылаться на эту статью следующим образом:
Гнездилов В.А., Русских С.В. Движение абсолютно жесткой двухосной тележки с упругим телом по двухрельсовому криволинейному полотну // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2016. №2. C. 93-112. DOI: 10.18698/0236-3941-2016-2-93-112
Please cite this article in English as:
Gnezdilov V.A., Russkikh S.V. The motion of an absolutely rigid biaxial trolley with an elastic body along a double rail curvilinear track. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr. [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 2016, no. 2, pp. 93-112. DOI: 10.18698/0236-3941-2016-2-93-112