Динамическая модель для расчета пружинного пакета волнового шагового двигателя Текст научной статьи по специальности «Физика»

УДК 621.8

С. С. К л е н и к о в, Т. А. Фомина

ДИНАМИЧЕСКАЯ МОДЕЛЬ ДЛЯ РАСЧЕТА ПРУЖИННОГО ПАКЕТА ВОЛНОВОГО ШАГОВОГО ДВИГАТЕЛЯ

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

E-mail: [email protected]

Ключевые слова: волновой шаговый двигатель, импульсы сил, волна деформирования, метод конечных элементов.

Основные понятия. Гибкое колесо волнового шагового двигателя (ВШД), описание которого приведено в работе [1], представляется в виде совокупности плоских концентрических круговых колец, соединенных между собой жесткими перемычками (рис. 1). Далее описан-

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

Рис. 1. Пружинный пакет:

1 — жесткая корпусная деталь; 2 — гибкое тонкостенное кольцо; 3 — гибкий зубчатый венец; 4 — жесткие перемычки

ная совокупность _______________ как пружинный пакет.

4

4

пакета использовать теорию плоских круговых колец с учетом условии их сопряжения.

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

/ d6 v d4 v d2 v

d2

m

= — r

d^2

+ m — r

' dqr N

т~ + qt

(1)

где V — окружные перемещения; ^ — угловая координата; — из-гибная жесткость соответствующего кольца; г — радиус его средней

dv 1 ( d2 V \

линии; т = —---радиальные перемещения; в = - V + -— —

г \ d(p2)

угол поворота нормали; — радиальная распределенная нагрузка; д — тангенциальная, или окружная распределенная нагрузка; т — распределенный изгибающий момент (рис. 2).

Матрицы жесткости, масс и демпфирования конечного элемента. Решение динамической задачи системы из к колец (к = 1, 2, 3 ... — число колец в пружинном пакете) предполагает решение матричного дифференциального уравнения

[МШ + [С]{д} + [К]{д} = {Р}, (2)

где [М] — матрица масс; [С] — матрица сопротивления (демпфирования); [К] — матрица жесткости; {Р} — вектор внешней нагрузки на систему; {(¡}, {д}, {д} — векторы обобщенных ускорений, скоростей и перемещений узлов конечных элементов.

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

Рис. 2. Кольцо пружинного пакета

Глобальные матрицы жесткости и масс формируются из соответствующих локальных матриц для элементов колец и жестких перемычек.

Для решения стационарной задачи (1) в работе [2] использована конечно-элементная модель. Конечным элементом в модели является криволинейный стержень — дуга радиуса r с углом, изменяющимся от до <^2. Функции формы имеют вид

Nk = Nk (<р, r, Ck) = Cki + Ck2

+ (Ck3 + Ck4^) cos ^ + (Ck5 + Ck6^) sin k = 1,..., 6, (3)

где Ck = Ck1,..., Ck6 — константы, определяемые из граничных условий для каждого конечного элемента.

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

(Si)

¿2 S3

w i ..... ¿4

¿5 ¿6

= [N(p)] i

(4)

где уг = ¿1; Wi = ¿2; вг = ¿3; vj = ¿4; Wj = ¿5; = ¿6; v¿, Wi, v, — узловые перемещения конечного элемента;

[N (*>)] =

Wj

Ni(p) dNi(p)

dp

N2(p) dN2(p) dp

Ns(p)

dNa(p) dp

N4(p)

dNj(p)

dp

N5 (p)

dN5(p) dp

N6(p) dN6(p) dp

В работе [3] показано, что элементы матрицы жесткости для элементов каждого кольца можно записать как

EJ

V2

kij = — / {(Ni" + Ni) (Nf + Nj)} dp, i = 1,..., 6, j = 1,..., 6.

vi

(5)

При построении зависимости для вычисления элементов матрицы масс в настоящей работе используется формула

mij = / [Ni]Tp[Nj]dV,

(6)

Ve

где Ve — объем конечного элемента; dV = F dl = Fr dp (F — площадь поперечного сечения того или иного кольца).

В предположении, что плотность колец постоянна (р = const), выражение для элемента матрицы масс запишется в виде

'-Р2 '-Р2

тг] = pFr J[Щт[р]dp = pFr J (nM + ^dP ^

Аналогичную формулу для элемента матрицы масс можно получить, используя выражение для кинетической энергии

E = 1 (S }т[М](5}. (8)

Для отдельного элемента кольца формулу кинетической энергии можно представить как

'-Р2

E = pF J(w2 + v2)rdp.

<pi

Из формулы (4) следует, что

6 6 dN.

v = J2 NiSi, w = Si,

i=l i=l

v2 = £ £ Ni M SSj, w2 = £ £ dp dpU,-

i= 1 j=l i=l j = l Y Y

Таким образом,

E = j{t t »N^ + tt f fl«P =

<f>1 % 3

pFr dN% dNj

-— (??( N% N3 8% S3 + , % , 3

2 3=i ^ % 3 % 3 d, dp

<pi

8%Sj I dp.

Сравнивая полученное выражение для кинетической энергии с формулой (8), получаем то же выражение для элементов матрицы масс (7).

Для однородного материала матрицу внутреннего демпфирования можно вычислить через матрицу масс, умножая ее на р, где ^ — коэффициент демпфирования. Таким образом, элементы матрицы демп-

фирования вычисляются по формуле

'-Р2

Cij = ^Fr / NiNj + ) dp;

¿я» Щ ^р ^р

VI

элементы матрицы жесткости вычисляются через функции формы следующим образом:

V2

Е 7 С

кг3 = Е-^ {(N1" + + ¿р, г = 1,..., 6, ^ = 1,..., 6.

VI

(9)

Элементы матрицы масс через функции формы вычисляются с помощью выражения

V2

т„ = + ^ (ю)

VI

элементы матрицы внутреннего демпфирования вычисляются через функции формы по формуле

V2

Cij = ^Fr / NiNj + -^-fj ) dp. (11)

^р ^р

VI

Матрицы жесткости, масс и демпфирования пружинного пакета. С учетом изложенного, пружинный пакет можно представить в виде совокупности колец, соединенных между собой перемычками. Для каждого кольца в отдельности можно составить матрицы жесткости, масс и демпфирования, используя соотношения (9)—(11) для отдельных конечных элементов. Далее необходимо объединить матрицы для отдельных колец в единую систему для всей конструкции. Полученные соотношения должны учитывать влияние колец друг на друга за счет перемычек, которые могут быть расположены неравномерно.

Перемычка при деформировании пружинного пакета движется как единое целое, поэтому, зная значения перемещений крайних точек перемычки, можно определить положения всех остальных точек.

Иными словами, в месте стыковки колец пружинного пакета должны выполняться следующие условия (рис. 3):

VI + £01 = г>2, = ^2, 01 = 02.

При модифицировании матрицы жесткости используется способ, предложенный в работе [2]. Он основан на введении в матрицу жесткости стыковочных условий с помощью больших чисел.

Матрица масс для перемычки вычисляется по формуле (6). Для перемычки ^^ = Б (Б — площадь поперечного сечения перемычки),

Рис. 3. Перемычка

тогда

0 0 x 0 0 ' Vi '

L

wi

1 1 6i

0 0 0 <

2 2 V2

1 1 W2

0 0 0 Л

2 2_

Г 2

I [Ni ]т [Nj ]d

ri

x.

Из условий стыковки для смежных колец в зоне перемычки следует, что

/ ^2 - VI \ (л х\ X

V = VI + х01 = VI + XI —£— I = 11 - ^ VI + ^2;

1 1

™ = = ^2 = ^ + 2 ^2;

0 = 01 = 02 = 101 + 102;

' v 4

w б

= [N (x)]

Vi Wi 6i V2 W2

- L)

0 0-00

0 0-0

0-00 2

Vi W1 6i V2 W2

v62,

Матрица масс для перемычки будет иметь вид

0

0

[m] = pS x

x

1 r3 — Г3 3 Ar2

0 0

1 r3_r3

0 2 <r- + r' >— 3-Ar^

00

-Ar 0 4

0 -Ar 4

1 1 r3 r3 2 (r- + r >— 3 rA5i 0

0 1 Ar

4

0

0

0

1

-Ar 4

1 r3 — r3 3 Ar2

0

-Ar 0 4

0 -Ar 4

00

1

-Ar 4

0

0

1

-Ar 4

где Тг и г^ — радиусы большего и малого стыкующихся колец,

Аг = Тг — Тз .

Выражение для элементов матрицы демпфирования перемычки аналогично выражению для матрицы масс.

Матрицы масс и демпфирования для каждого из колец и всех перемычек объединяются в глобальные матрицы системы. Таким образом, строится система дифференциальных уравнений (2).

Численные методы решения задачи динамики. Численное решение системы дифференциальных уравнений (2) возможно несколькими способами, например методом линейного ускорения Вилсона или методом трапеций Кранка-Николсона. В разработанном программном комплексе предпочтение отдается методу обобщенного ускорения Ньюмарка [3].

На каждом шаге определяются значения перемещений, скоростей и ускорений сначала узловых точек, а затем и внутренних. Для получения гладкого решения достаточно п конечных элементов (п должно быть равно числу точек приложения сил, в данном случае п = 8) [4].

Значения переменных по времени и направлению импульсов сил в правой части системы (2) являются расчетными. Их значения, вычисляются с помощью модели, описанной в работе [4]. Возможные формы импульсов подаваемых сил (давлений) приведены на рис. 4.

При расчете предполагается, что силы прикладываются к внешнему кольцу (Ы3) пружинного пакета. Внутреннее кольцо (Ы1) является толстостенным, поэтому при расчете оно считается жестко закрепленным (рис. 5). Среднее кольцо (Ы2) получает упругие перемещения через перемычки от внешнего кольца.

0

0

0

0

0

Рис. 4. Гистограмма импульсов сил, подаваемых в полости пружинного пакета (п = 50 Гц):

1, 2, 3 и 4 — отверстия № 1, № 2, № 3 и № 4 соответственно

Как показывают расчеты, деформированная форма кольца не только поворачивается, но и видоизменяется, что обусловливается дискретностью самой системы и дискретностью воздействующих на нее импульсов сил. Кроме того, форма кольца видоизменяется в зависимости от угловой скорости подачи импульсов сил (т.е. частоты вращения внутренней части распределителя) [4].

Рис.5. Деформированная фор- Без учета демпфирования при угло-ма колец пружинного пакета вых скоростях подачи импульсов сил,

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

На рис. 6, а приведены значения радиальных перемещений w колец пружинного пакета в момент времени £ = 0,001 с. Частота вращения внутренней части распределителя в рассматриваемом случае принята равной 50 Гц. Там же приведен график зависимости радиальных перемещений, полученный с помощью динамического метода, описанного в работе [2]. Как следует из рис. 6 а, при заданной частоте амплитуды различаются на величину порядка 2... 2,5 %.

Графики радиальных перемещений колец в момент времени £ = 0,001 с при увеличении частоты вращения распределителя и соответствующих импульсов сил приведены на рис. 6, б,... ,д. На рис. 6, б,в

W, мм

О 1 2 3 4 5 6 <р, рад

а

_I_I_I_I_I_I_I_I_I_

О 0,001 0,001 0,002 0,002 0,003 0,003 0,004 0,004 0,005 t, с

W, мм

в

Рис. 6 (начало)

г

w, мм 0,3

0,2

ОД 0

0 1 2 3 4 5 6 ср, рад

Рис. 6 (продолжение)

частота вращения внутренней части распределителя принята равной 100 об/с, на рис. 6, г,д — 167 об/с.

Графики радиальных перемещений внешнего и второго колец пружинного пакета ВШД при частотах вращения внутренней части распределителя, равных 50, 100 и 167 Гц приведены на рис. 6, е,ж. В момент времени £ = 0,0025 с на кольца действуют примерно одинаковые по значению импульсы сил. Из графиков следует, что частота вращения внутренней части распределителя сильно влияет на значения радиальных перемещений колец и, соответственно, на их деформирование.

W, мм

0,24 0,20 0,16 0,12 0,08 0,04 0

-0,04 -0,08 -0,12 -0,16 -0,20

О 1 2 3 4 5 (р, рад

е

Рис.6 (окончание). Радиальные перемещения внешнего (N3) и среднего (N2) колец пружинного пакета и гистограммы импульсов подаваемых сил (б, г) при t = 0,001 с (а-д); 0,0025 с (ж); п = 50 (а); 100 (б, в); 167 Гц (г, д):

1,2,3,4 — силы N1, N2, N3, N4 соответственно

Заключение. Приведена последовательность формирования глобальных матриц жесткости, масс и демпфирования для решения методом конечных элементов динамической задачи о деформировании пружинного пакета. При этом в задаче учтено, что импульсы деформирующих сил переменны и их значения в каждый момент времени вычисляются с помощью приложения, приведенного в работе [4]. Ре-

зультаты вычислений приведены в виде графиков радиальных перемещений точек деформируемых колец при разной частоте вращения внутренней части распределителя [4].

СПИСОК ЛИТЕРАТУРЫ

1. Клеников С. С., Майков А. И. Разработка математической модели волнового редуктора с круговой формой зубьев и упругими звеньями в виде кольцевых пружинных пакетов // Тез. докл. и выступлений студентов, аспирантов и преподавателей на VII Межвуз. науч.-практ. конф. - М.: МГИУ, 2006. - С. 327335.

2. Клеников С. С., Майков А. И. Разработка математической модели расчета пружинного пакета волнового редуктора методом конечных элементов (МКЭ) // Изв. вузов. Машиностроение. - 2008. - № 11. - С. 25-30.

3. Голованов А. И., Тюленева О. Н., Шигабутдинов А. Ф. Метод конечных элементов в статике и динамике тонкостенных конструкций. -М.:Физматлит, 2006. - 392 с.

4. Фомина Т. А. Математическая модель рабочего вещества волнового шагового двигателя // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. - 2010. -№ 4.-С.98-107.

Статья поступила в редакцию 6.10.2010

Сергей Сергеевич Клеников родился в 1942 г., окончил МВТУ им. Н.Э. Баумана в 1971 г. Д-р техн. наук, профессор, заведующий кафедрой "Техническая физика и прикладная механика" в филиале Московского государственного университета приборостроения и информатики в г. Сергиев Посад. Автор около

É100 научных работ в области динамики и прочности машин, приборов и аппаратуры.

S.S. Klenikov (b. 1942) graduated from the Bauman Moscow Higher Technical School in 1971. D. Sc. (Eng.), professor, head of "Technical Physics and Applied Mechanics" department of Sergiyev Posad Branch of the Moscow State University. Author of about 100 publications in the field of dynamics and strength of machines, instruments, and apparatus.

Татьяна Александровна Фомина родилась в 1985 г., окончила Московский государственный индустриальный университет в 2007 г. Аспирант, преподаватель кафедры "Прикладная математика и информатика" в филиале Московского государственного университета приборостроения и информатики в г. Сергиев Посад. Автор двух научных работ в области динамики и прочности машин.

T.A. Fomina (b. 1985) graduated from the Moscow State Industrial University in 2007. Post-graduate, teacher of of "Applied Mathematics and Information Technologies" department of Sergiyev Posad Branch of the Moscow State University. Author of 2 publications in the field of dynamics and strength of machines.