Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Теоретические основы динамики кольца при его обработке по мобильной технологии. Монография

.pdf
Скачиваний:
0
Добавлен:
07.09.2026
Размер:
2 Мб
Скачать
☆
41
центробежная сила становится больше силы упругости и возникает расходящийся процесс прецессионный резонанс.
Подобный эффект описан в работе [1] для конической свободной оболочки, совершающей сложное вращение без учета внутреннего трения. Для вращающихся на опорах кольца, оболочки описания этого явления и его исследование в научной литературе автор не нашел.
42
ij i j
S G
− + = λ +
3. Метод получения численного значения элементов матриц системы дифференциальных уравнений
При получении решений для технических задач, процессы в которых описываются дифференциальными уравнениями, необходимо выполнить два этапа. Первый – составление дифференциальных уравнений, описывающих поведение системы, и второй – разработка программы для их решения и полу­чения интересующих результатов. При выполнении первого этапа в случае сложных систем возникает возможность ошибок при аналитических выкладках, что приводит к неверным уравнениям и большим потерям времени на поиск ошибок.
Имеется метод получения численного решения системы обыкновенных дифференциальных уравнений без знания аналитических формул для этих уравнений [6], основанный на использовании уравнения Лагранжа второго рода. При его использовании необходимо иметь аналитическое выражение для кинетической и потенциальной энергий системы, для которой необходимо получить решение. В методе получают сразу численное значения правой части системы уравнений, необходимое для численного интегрирования. Его определяют с помощью численного дифференцирования и интегрирования выражений для кинетической и потенциальной энергии.
Предлагаемый в данной работе метод применим для линейных систем, и отличается от рассмотренного тем, что позволяет получить численные значения матриц при обобщенных координатах и их производных без знания формул для них. Это может использоваться для анализа свойств исследуемой системы или сравнения значений матриц с полученными аналитическим способом, с целью контроля правильности аналитических выкладок.
Рассмотрим метод получения решения для технических систем, в частности оболочек, колец, описываемых линейными дифференциальными уравнениями, не требующий первого этапа - получения формул для коэффициентов матриц системы уравнений. Метод основан на использовании уравнения Лагранжа второго рода, а также на том, что известна структура окончательного уравнения.
При использовании формализма Лагранжа требуемую систему диффе­ренциальных уравнений получают по зависимости
d T T
∂ ∂ ∂∏
&
dt q q q
∂ ∂ ∂
j j j
K
∑
i
1
=
, (3.1)
D1(q) = 0, D2(q) = 0,…, DK(q) = 0, (3.2)
где , – соответственно кинетическая и потенциальная энергия; q
j -
j – я
координата, (j=1,2,…,N); Sij – производная i – го уравнения связи Di по j – й
координате, (i = 1,2,…,K); v
- неопределенный множитель Лагранжа; Gj –
внешняя сила, действующая по j – й координате.
Рассмотрим метод при отсутствии правой части. Вычисляя производные из уравнения (3.1) и группируя одинаковые слагаемые, получим в общем случае окончательную систему уравнений
A
2 tÅ.
+ A
+ A0 q=0, (3.3)
1 t:
43
0.
+ + =
( , , , )
T F q q V t dV
( , , )
F q V t dV
( , , , )
dV
dVq
= +
&
( , , , ) ( , , , )
dVq dV
( , , , ) ,
F q q V t dV
( , , ) ,
F q V t dV
,
dV
∂
∂
( , , )
dV
∂
∂
где A2, A1, A0 – матрицы постоянных или переменных коэффициентов, которые полностью определяют характер решения и формулы для которых ищут при аналитических выкладках. Для j – й строки системы (3.3) можем записать
N N N
ji ji ji
&& &
A q A q A q
∑ ∑ ∑
2 1 0
i i i
1 1 1
i i i
= = =
(3.4)
Формулы для кинетической и потенциальной энергии тела, в общем случае, имеют вид
&
T
, (3.5)
, (3.6)
=
=
∫
V
∫
V
где V – координаты, задающие объем тела, по которым вычисляем полную энергию; t – время.
Наибольшую проблему по трудоемкости и возможности ошибиться представляет аналитическое вычисление интегралов в (3.5), (3.6), производных от них и нахождение результирующих выражений при получении формул для элементов матриц системы дифференциальных уравнений. Преобразуем первое
слагаемое в (3.1) с целью выделения tçÅ. Получим
&
T
,
2 2
& &
F q q V t F q q V t
T T
& & &
q q q t
∂ ∂ ∂ ∂
j i j
&&
+
i
. (3.7)
2
 
& &
dt q q q
∂ ∂ ∂
N
∑
∫
1
i
=
j j i
V
F q q V td T
∂∂
( , , , )
T
F q q V t
T
∂
& &
q q
∂ ∂
&
i
i
∂
=
∫
j j
V
N
∂ ∂
∑
∫ ∫
1
i
=
V V
Величины остальных производных в (3.1) можно вычислить по зависимостям
T
∂ ∂
q q
∂ ∂
либо
Выбор конкретных выражений из этих зависимостей определяется простотой вычислений, которая зависит от вида FT, F.
При вычислении вторых производных от F использовать аналитические или численные методы. Интегралы вычисляем численным методом.
Таким образом, левая часть уравнения (3.1) представлена в виде, когда в нем фигурирует вектор – функция q – решение системы, и ее производные
=
j j
=
q q
∂ ∂
j j
T
∫
V
( , , , )
F q q V tT
T
∫
V
&
&
∂ ∂
q q
∂ ∂
q q
∂ ∂
j j
=
j j
F q V t
=
∫
V
и первых от FT, F можно
T
∫
V
. (3.8)
44
j j j
dt q q q
∂ ∂ ∂ ∏
∂ ∂ ∂
i i i j
A q A q A q Ñ
ij i j
A q Aq A q S G
+ + = λ +
, ( 1,2,..., ).
j N
t: tÅ. Они нам неизвестны. Зададим текущий момент времени t
, для
к
которого надо вычислить матрицы коэффициентов системы дифференциальных уравнений, а также зададим произвольными числами вектор – функцию решения q
и ее производные t
€
€
:
Å
, t
. Подставим их в (3.7), (3.8) и вычислим
зависимости. Результат подставим в левую часть (3.1). Получим число
C
не равное нулю. Подставляя t€ ,t
d T T
= − +
j
&
€
€
:
, t
Å
, в левую часть (3.3) и приравнивая j
, (3.9)
получим линейное алгебраическое уравнение
N N N
ji Ï ji Ï ji Ï
&& &
∑ ∑ ∑
1 1 1
i i i
= = =
в котором неизвестными являются элементы матриц W разных значений переменных t€ ,t
+ + =
2 1 0
€
€
:
Å
, t
получим систему из 3N алгебраических
, (3.10)



 W
. Задавая 3N
, W
уравнений с 3N неизвестными. Решая ее, найдем численное значение элементов j – й строки матриц 2, 1, 0. Осуществляя подобные вычисления для каждой строки этих матриц получим их числовые значения, которые ис­пользуем для получения требуемого решения системы дифференциальных уравнений (3.3).
Рассмотрим применение этого метода в случае уравнений со связями (3.1). Основная трудоемкость в получении матриц 2, 1, 0 выражения
A
2tÅ.
+ A
+ A0 q,
1t:
соответствующего левой части (3.1). Численные значения элементов матриц можно вычислить таким же способом, как и в случае уравнения без правой части. Значения элементов Sij в (3.1) легко могут быть вычислены по урав­нениям связи (3.2). Вычисление сил Gj не представляет труда. В результате получаем систему
k
&& &
2 . 1 0
∑
i
=
1
(3.11)
с уравнениями связи (3.2). Используя матричные преобразования систему (3.11) совместно с уравнениями связи можно привести к системе с N-K уравнений с N-K неизвестными и получить требуемое решение.
Рассмотрим применение этого метода для случая, когда в системе возникают возмущающие силы. Это будет при наличии неоднородностей в оболочке, кольце. Формулы для вычисления этих сил получим, вычисляя выражения
В этом случае левая часть (3.1) после преобразования примет вид
tÅ + A1 t:
A
2
+ A0q + R, где R – вектор возмущающих сил. Уравнение (3.10)
для определения элементов матриц тогда будет
d T T
∂ ∂
 
dt q q
&
∂ ∂
j j
=
45
.
i i i j j
1
q
&&
2
q
&&
N
q
&&
1
q
&
2
q
&
q
&
1
q
2
q
q
1
W
2
W
3
N
W
3 1
+
N N N
ji ji ji
&& &
A q A q A q R
∑ ∑ ∑
1 1 1
i i i
= = =
+ + + =
2 1 0
(3.12)
В нем появилась дополнительная неизвестная Rj, соответствующая неизвестному значению возмущающей силы. В этом случае вместо
задаваемого вектора t€ ,t
€
:
Å
, t
из 3N компонент будет вектор t€ ,t
€
:
, t
€
Å
и еще
€
одна компонента , т.е. вектор из 3N + 1 элементов. Значение последнего элемента всегда равно единице. Задавая систему из 3N + 1 линейно
€
€
:
независимых векторов t€ ,t получим систему линейных уравнений для нахождения элементов матриц W

W

 W
и вектора возмущающих сил Rj.
При выборе числовых значений векторов t€ ,t
Å
, t
, , вычисляя соответствующие им 1
€
€
:
Å
, t
необходимо обеспечить
 
хорошую обусловленность матрицы коэффициентов системы линейных алгебраических уравнений. Получить это можно заданием ортогональной сис­темы векторов W. Однако, в нашем случае задать ортогональную систему нельзя, т.к. компонента всегда равна единице. Наиболее удобной будет система, представленная в таблице 3.1.
Таким образом, рассмотренный метод позволяет быстро получить решение, не затрачивая время на аналитические преобразования с целью получения формул для матриц коэффициентов. Он может применяться при исполь­зовании физических и обобщенных координат для получения уравнений.
Это целесообразно в том случае, когда важно получить численное значение элементов матриц, а знание формул для коэффициентов матриц не имеет значения. Метод может использоваться для контроля правильности полученных аналитических формул для матриц коэффициентов системы линейных дифференциальных уравнений.
Таблица 3.1
Задание системы векторов
W
…
…
N
…
N
1 0 0 0 0 0 0 0 0 0 0 0 1
,
0 1 0 0 0 0 0 0 0 0 0 0 1
. . . . . . . . . . . . .
0 0 0 0 0 0 0 0 0 0 0 1 1
W
0 0 0 0 0 0 0 0 0 0 0 0 1
N
46
4. Pезультаты физического и вычислительного экспериментов по динамике кольца
Для проверки адекватности полученных уравнений динамики вращающихся на опорах оболочки и кольца реальным физическим процессам были проведены физические и вычислительные эксперименты. Физический эксперимент проводился не на реальных оболочках и кольцах, а на их физических моделях вследствие проблематичности использования реальных объектов. Результаты, полученные на физических моделях, могут быть пересчитаны на реальные объекты с использованием теории подобия дина­мических процессов. Для проведения физических экспериментов были изготовлены металлическое кольцо и оболочка, а также устройство для их вращения на двух опорах.
Необходимо сказать, что используемые при расчетах геометрические и физические параметры указанных тел являются некоторой идеализацией. На самом деле эти параметры в реальных телах меняются в зависимости от угловой и линейной координаты точки тела, т.е. непостоянны. Объясняется это сложностью изготовления с выдерживанием точных параметров. Так толщина используемого в физическом эксперименте кольца меняется по периметру в интервале от 1,1 до 1,3 мм, а его диаметр средней линии от 371,5 до 372,5 мм. Диаметр срединной поверхности одной стороны используемой оболочки меняется в интервале от 372,5 до 373мм, а ее толщина от 0,7 до 0,8 мм. Для другой стороны диаметр меняется от 372 до 373 мм, а толщина от 1,1 до 1,4 мм. Поэтому при проведении расчетов брались усредненные параметры. Для кольца радиус средней линии r = 0,18619 м, толщина h = 0,001 м, ширина a = 0,0605 м, модуль упругости = 2 1011 н/м1. Радиус срединной поверхности оболочки r = 0,1874м, толщина h = 1,2 10-3м, ширина = 0,2м, модуль упругости = 2 1011н/м1. При проведении эксперимента по определению собственных частот использовались два подхода. Первый заключается в том, что задают начальные возмущения кольцу или оболочке по первой или второй частоте и с помощью измерительных средств определяют собственную частоту. Оказалось, что такой подход может использоваться только в случае невращающихся кольца или оболочки. Объясняется это тем, что изготовленные реально, они являются телами с непостоянными параметрами, т.е. радиус, толщина, модуль упругости их непостоянны, а зависят от координат точки. Это обстоятельство ведет к возникновению возмущений в динамике вращающихся кольца, оболочки. Кроме того, на динамику вращающихся кольца, оболочки влияет переменный радиус, погрешность в задании цилиндрической формы опорных роликов, погрешности в их параллельной выставке, нестабильность угловой скорости вращения, обусловленной нестабильностью питающего напряжения и величины внутреннего трения в приводе. Эти обстоятельства не позволяет выделить собственную частоту колебаний, возбужденную с помощью начального возмущения, на фоне действия возмущающих факторов от неидеальности физических и геометрических характеристик. Поэтому использовался резонансный способ возбуждения требуемой частоты.
47
Рис.4.1. Установка для физического моделирования
Резонансная частота возбуждалась с помощью электромагнита, на который подавался периодический сигнал, частота которого задавалась низкочастотным генератором периодических колебаний. Величина напряжения периодического сигнала, снимаемого с генератора, недостаточна для получения
48
-
2
-
2
-
4
-4
-4
-4
-
4
-4
-
4
-4
электромагнитных сил, вызывающих возникновение резонанса в кольце, оболочке. Поэтому периодический сигнал подавался на управляющую обмотку реле, контакты которого включены в цепь, содержащую источник регули­руемого постоянного напряжения и обмотку электромагнита для воздействия на кольцо или оболочку.
На ней были проведены эксперименты для определения динамических характеристик вращающегося на опорах кольца при действии различных факторов.
4.1. Определение коэффициента внутреннего трения кольца с опорами
В этом эксперименте вычислялся коэффициент внутреннего трения невращающегося кольца на двух опорах. На физической установке возбу­ждались колебания кольца по первой собственной частоте и с помощью измерительных средств записывался переходной процесс, представленный на рис. 4.2
Затем, по разработанным в данной работе математическим моделям дина­мики кольца с растяжимой и нерастяжимой средней линией вычислялись коэффициенты внутреннего трения, задающие такой же переходной процесс, как и в физическом эксперименте. Результаты представлены в табл. 4.1.
Рис.4.2. Затухающий процесс колебаний кольца
Таблица 4.1
Коэффициенты внутреннего трения
Угол между
11,96 25 50,62 64,96
опорами, град Коэф. затух. перех.
- 4,97 10
- 7,60 10
- 0,110 - 0,116
проц.в эксперим
Коэф. трен н с/м
Нерастяж. сред. лин.
Растяжим.
2
сред. лин.
0,98 10
0,66 10
1,38 10
0,95 10
1,52 10
1,12 10
1,32 10
1,04 10
Из таблицы видно, что отличие коэффициентов внутреннего трения для моделей с учетом растяжимости средней линии и без учета составляет, при-
49
-
4
τ
1
ω
τ
2
ω
τ
1
ω
τ
2
ω
-
11
-
10
10
-9
-
11
-12
-
12
-12
-
10
-12
-
11
-13
-
12
-11
-12
-14
-
11
-9
-10
-11
-
11
-11
-
11
мерно, 20 процентов. Различие коэффициентов трения для разных углов между опорами объясняется тем, что частота колебаний зависит от угла α и рассматривается частотнозависимая модель трения. Полученные значения коэффициентов использовались для расчетов динамики кольца и оболочки в зависимости от угловой скорости его вращения.
4.2. Определение первой и второй собственных частот колебаний кольца в зависимости от угловой скорости его вращения
Как уже указывалось, анализ полученных в данной работе дифференциальных уравнений поведения кольца позволил автору открыть новый физический эффект – уменьшение до нуля собственной частоты враща­ющегося на опорах кольца с увеличением угловой скорости вращения до неко­торого критического значения, и возникновение прецессионного резонанса при превышении скоростью этого значения.
Для подтверждения этих результатов автором были проведены расчеты, по полученной математической модели, влияния угловой скорости вращения кольца Ω на величину первой и второй собственных частот его колебаний. Эти результаты представлены в таблице 4.2 для случаев коэффициент внутреннего трения ТР = 0, и коэффициент внутреннего трения равен вычисленному в эксперименте, т.е. ТР = 1,52 10
-4
н с/м. В таблице обозначено: è,è
-
показатели экспоненты, характеризующие вид колебательного процесса по
первой и второй собственной частоте соответственно; é, é
- соответственно
первая и вторая собственные частоты.
Таблица 4.2
Прецессионный резонанс
ТР = 0 ТР = 1,52 10
Ω
0 9 10
2 10
20 1 10 30 1 10
40 1 10 50 1 10
60 9 10 70 1 10
75 3 10 79 4 10
80 2 10
1
39,78 6 10
39,33 3 10 37,92 8 10
35,50 5 10 31,88 4 10
26,84 1 10 19,97 6 10
10,66 5 10 4,72 2 10
0,92 4 10 2,50 8 10
2
1
2
82,64 -0,12 39,78 -0,52 82,64 82,06 -0,13 39,33 -0,52 82,06
80,32 -0,17 37,93 -0,51 80,32 77,42 -0,23 35,50 -0,51 77,42
73,38 -0,34 31,89 -0,50 73,38 68,21 -0,48 26,84 -0,48 68,21
61,87 -0,69 19,98 -0,46 61,89
54,20 -0,99 10,68 -0,41 54,26
49,69 -1,20 4,74 -0,34 49,77
0,92 -1,45 0,83 -0,24 45,68 2,50 -1,53 2,47 -0,16 44,52
90 10,85 26,58 9,54 26,58 -10,72 24,39 8,77 29,00 100 21,26 34,93 20,77 34,93 -20,20 32,33 21,60 36,66
50
0 2 4 6 8 10 12 14 16 18 20
8
Из таблицы для  = 0 видно, что при увеличении угловой скорости вращения от нуля до, примерно, 79 1/с, назовем ее критической, первая
собственная частота колебаний кольца é уменьшается до нуля, а величина è
,
характеризующая расходимость или сходимость периодического
процесса, равна нулю. Равенство è
нулю при наличии в дифферен-
циальном уравнении матрицы коэффициентов при первой производной происходит вследствие того, что матрица является кососимметричной и задает действие гироскопических сил. При угловой скорости вращения больше критической значение τ1 увеличивается, процесс делается расходящимся и начинает возрастать первая собственная частота. Поведение второй соб­ственной частоты при
= 0 при увеличении Ω аналогично поведению

первой собственной частоты. Критическая точка является точкой бифуркации. Из табл. 4,2 при  не равном нулю, видно, что внутреннее трение
существенно влияет на характер процесса. Не влияя на характер изменения первой собственной частоты, оно препятствует возникновению резонансного эффекта, делая колебательный процесс по первой собственной частоте затухающим.
Наличие трения ведет к изменению характера колебаний по второй
собственной частоте. При Ω = 0 коэффициент è
отрицательный, с увеличением
Ω он увеличивается и в окрестности точки бифуркации делается равным нулю, а затем резко растет, т.е. процесс по второй собственной частоте делается расходящимся.
На рис. 4.3, 4.4 представлены, для колец радиусов 5м, 4м, 3м, зависимости
коэффициента è собственной частоты é
расходимости колебательного процесса, а также первой
в рад/с от скорости вращения кольца Ω в рад/с.
6
4
Рис. 4.3. Зависимость показателя расходимости колебательного
процесса от скорости вращения кольца
Число гармоник ряда Фурье равно 6. Вычисление этих характеристик осуществлялось решением частотного уравнения численным методом. Коэффициент внутреннего трения равен нулю, толщина колец 0,05м, угол между опорами 40 градусов. Цифрой 1 обозначена зависимость для кольца радиуса 5 метров, цифрой 2 – радиуса 4 метра, цифрой 3 – радиуса 3 метра. Из
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]