книги / Микропроцессорное управление технологическими процессами в радиоэлектронике
..pdfСледует отметить, что В — нелинейный оператор и ему соответствует нелинейная полугруппа, тем не менее метод Г-произведений можно применять в этом случае [30]. Точное решение задачи Коши (3.50) имеет вид
и (0 = |
lim п |
(ДОS® (ДОи „ |
|
N-+°° k=\ |
|
где числа /V, At |
и t связаны соотношением t = nAt. |
Приближенное решение задачи Коши (3.50) (а сле довательно, и начально-краевой задачи (3.34), (3.36), (3.49)) получается в виде конечного произведения при фиксации N
S ( 0 - S S*(A Q S® (A Q «.- |
(3.51) |
В уравнении (3.51) 5 (l) (ДО — полугруппа (в дан ном случае группа) операторов, отвечающая задаче
Коши |
|
|
|
|
|
|
|
|
~ |
= Ай, и (0) = |
|
|
(3.52) |
||
S® (Д/) — нелинейная группа операторов в |
отве |
||||||
чающая |
задаче Коши |
|
|
|
|
|
|
|
|
- = Ви, и (0) = |
ы0, |
(3.53) |
|||
где At = |
t/N (N — натуральное число), |
|
|||||
|
|
/ оо |
оо |
|
|
|
|
5 (1) (т) (/ (х, у), g (х , </)) = ( £ |
|
Е (5™ cos У Х й от + |
|||||
|
|
\т=1 п=1 |
|
|
|
||
“Ь Т$пт Sin ]/*Хт п0Я) Vnm (%* У)* |
|
X X |
|
( ПV ^пт В п щ Х |
|||
|
|
|
т=1 Пэа1 |
|
|
||
X sill V hnm ПТ “}” П |
knm Bntn COS |
Хпт ПТ) 1)^ (jf, у))| |
91
s m (t, r) ( { (x, y), g (x, у)) = |
(/(*, y), g (x + |
2tiT, x) + |
' t 1 |
2o (< — !),(/) d|), |
(3 55) |
+ £ F(6, * + |
t
здесь коэффициенты £ „m, £nm, Ялт и функции цлт (*, у) вычисляют по формулам
В пт — |
|
|
|
|
|
J / ) s i n - ^ X |
|
^ |
» |
Л |
|
j p * ’ ’ |
|
||
|
|
|
|
тя . . |
|
|
|
|
|
X sin — |
ydxdy; |
|
|||
В пт — |
|
/ |
|
|
7 |
|
1Уо |
|
|
|
|
|
|
|
|
|
X |
. |
ля |
. |
/ля |
, , |
|
|
sin |
- j — A: sm — |
ydxdy; |
||||
Vmn (X, У) ■- |
I / |
|
V |
\ a J } . пп |
. mn |
||
V |
-------jjj-------sm — |
* s m — y. |
Неоднородные дифференциальные уравнения с пе ременными коэффициентами. Описанный метод мож но применять для приближенного решения уравне ния (3.34), когда коэффициенты а и v зависят от t
1«г( 0 - « а (0) I F |
+ |
* 2 W - 0 - + 2* ( 0 |
т г - |
- ^ |
W |
= F ( x , y , t ) - |
(3.56) |
92
Решение задачи (3.35) с граничными условиями (3.49) получим в виде
N—1
и (/, х ,у ) = П S <1) ((i + 1) Д(, Ш ) S°> ((i + 1) X /=0
хШ )(и 0(х, у), WQ xt у)),
где 5 (0 — эволюционные операторы (Д/ = tJN):
S m (Г, 0 (/ (х, у), g (х, у)) - |
( £ |
£ |
4 « (Ох |
|||||
X cos О / й Ж |
а (О (Г - |
0) + |
Sm„(0 х |
|||||
X sin ( / M |
f l |
а (0 ( Т - |
0» Олт(0 (*. У). |
|||||
£ |
( - а ( 0 К ^пт (0 Впт (t) Sin (VKm(t) X |
|||||||
«=1 |
|
|
|
|
|
|
|
|
X а (/) (7 - |
0) + а (0 К |
й |
й |
5 » |
х |
|||
х cos (V Кт (0 |
в (0 (7 — /))) о™. (0 (*, у)); |
|||||||
S'2' (7, |
о V (*. У), 8 (X, у)) = |
(/ (*, у), g (х + |
||||||
|
|
|
|
Т |
|
|
|
|
+ 2о(0 (Г - 0. У) + |
^ |
(6, * + 2» (0 (Г - |
||||||
|
|
|
|
t |
|
|
|
|
— * — Е). |
|
|
|
|
|
|
|
|
где коэффициенты |
|
(/), |
Bnm it), зависящие от / |
|||||
и g, коэффициент A,rtmи функции |
(0 (*. У) вычисля |
|||||||
ют по следующим формулам: |
|
|
||||||
В пт (0 — |
|
2 |
|
|
1Уо |
|
||
4 / |
1 т \» ( f Н *' d X |
|||||||
X sin |
» |
sin |
»0 |
|
|
|
|
|
93
~$nm(0 — ' |
|
|
|
|
'« w |
K w o / » . / |
|
l |
Уо |
n n |
(3.58) |
f |
f |
ntn |
|
x ) ) |
g { x > y ) |
sin — |
X sin — y d x d y ; |
|
’V |
- ( |
. |
i f)‘ |
ЯП |
. ntn |
|
|
Unm (0 (*, £/) = |
- |
> = ----- - sin —r~ JCsin ---y. |
||||||
|
|
Viy„ |
|
1 |
«о |
|
||
Формулы |
(3.57) |
и (3.58) |
дают решение |
поставлен |
||||
ной задачи, т. е. приближенное решение |
начально- |
|||||||
краевой задачи (3.36), (3.49) для |
уравнения (3.56). |
|||||||
Рассмотрим уравнение, |
описывающее |
колебания |
||||||
переменного |
вдоль |
ширины |
движущегося |
материала |
||||
натяжения а |
{у). Поскольку |
при моделировании про |
цесса функция а (//), описывающая это переменное натяжение, получается интерполяцией, обычно по двум точкам вдоль ширины материала, то ограничимся
семейством функций вида |
а {у) = |
с!(ау |
-+- Ь). |
Урав |
|
нение колебаний имеет вид |
|
|
|
|
|
( а Ч у ) - о г) ъ ? + а * ( у ) ^ - |
+ 2v |
д2и |
|
|
|
|
|
|
d x d l |
|
|
д ги |
г , . |
„ |
|
|
(3.59) |
|
|
|
|
|
|
Как и в предыдущем пункте, для решения началь |
|||||
но-краевой задачи с условиями и (0, х, у) = |
и0 (л:, у); |
||||
х, у) |/=о = |
(х, у); |
|
|
|
|
и (х, 0, () = ы (х, уа, 0 — и (0, y,t) = |
u (/, y ,t) |
= |
0 (3,6°* |
94
примени i метод Г-произведений. Для этого запишем условия (3.60) в виде эволюционного уравнения
= Аи + Ви + Си\
(3.61)
и (0) =
где операторы Л, В и С описываются правилами
А{и0, г , ) = |
(0,(ва0 / ) - ^ ) - ^ - ; |
|||
|
В (« 0, y o) = |
(o ,2 0 -^S- + F ); |
||
|
С (Ыо,1Г 0) = |
( г , а % |
) | ^ ) . |
|
Точное |
решение |
уравнения |
(3.61) и = (н0, W0) |
|
имеет вид |
|
|
|
|
и (0 = |
Пт П S® (ДОS® (ДО S® (ДОи (0) |
|||
|
Л'-*оо Л=1 |
|
|
|
и приближенное |
|
|
|
|
|
и (0 = |
(«о (0> ^о(0] = |
||
|
N |
|
|
|
= |
П S® (ДОS® (ДОS01 (ДО («0. W'o). |
|||
|
*=1 |
|
|
|
при этом функция и0 (t) задает решение задачи (3.56) —
(3.57). |
Во всех формулах |
S &) (т), |
S (2) (т) и 5 (1) (т) — |
группы |
преобразований фазового |
пространства |
|
|
Wl(lO, /) х fO, у0]) х |
В2([0, /] х [О, 0О]), |
отвечающие операторам С, В, А.
Группа нелинейных преобразований S {2) (т), от вечающая оператору В, описана формулами (3.57) и
95
(3.59). Группа |
S w (т) |
соответствует |
решениям |
на |
||||||
чально-краевой |
задачи |
|
|
|
|
|
||||
ди |
л. |
№ |
|
« 2 / .л |
ач |
д*и |
|
|
|
|
и (0, х, |
у) = |
и0; |
№ (0, |
г/) = |
UP0; |
|
|
|
||
,и (х, О, 0 = |
м (х, у0, () = |
и (О, y ,t) = u (/, |
*/,0 = |
0 |
||||||
и имеет вид |
|
|
|
|
|
|
|
|
|
|
5“’ М («.. ^о) = |
( “о. |
+ f |
(<*а (У) ~ |
о*) -здР-) • |
||||||
Наконец, |
группа |
S (3) (т) |
получается |
решением |
||||||
уравнения |
|
I H w f - P - |
|
(3-62) |
||||||
|
|
|
|
6 граничными условиями из уравнения (3.60) и началь ными условиями
и{о, X, у) = у0(х> у);
-§-(<. х,у) |<_о= «70(л:, у).
Вид переменного коэффициента позволяет сделать замену
|
1 |
ay + Ь \ |
|
(t,X — \n----1--- |
|
при которой соответствующее уравнение имеет вид |
||
dhi |
dhx |
(3.63) |
дР ~ |
дуа |
|
с начально-краевыми условиями |
|
|
и (0, х, у) ■=и0 (х, у) — (ау + Ь) и0 {х, |
||
ди |
|
Ь) W0(х, — х |
~ зг & X, у) |/„о = w 0(х, у) =* (ay + |
X In ау + Ь
Ь)
Q6
и на границе прямоугольника [О, I) х [0, (Ыа) ( ^ — 1)],
и — 0.
Обозначим (Ыа) (еа^ — 1) через у0. Тогда решение уравнения (3.63) можно получить следующим обра
зом: |
продолжая функции и0 и №0 с прямоугольника |
||
[0, |
Л X [0Г у$\ сначала на прямоугольник [0, /] X |
||
X [—у0, у0] нечетным образом по у, а затем перио |
|||
дически на [0, Л X |
[—оо, |
оо] (по у), получаем функ |
|
ции и0 на W [0, Л |
X I- о о , |
оо]. Формула д’Аламбера |
дает решение и уравнения (3.63) с требуемыми началь но-краевыми условиями
и (t, X, у) = -g- (и„ (X, у — 0 + и„(х, у + 0) +
Заметим, что поскольку нас интересует примене ние этой формулы при малых значениях t (ввиду ее использования в методе Г-произведения) достаточно взять_периодическое и нечетное продолжение функций
и0 и W0 лишь на прямоугольник [О, Л X [— е,г/0 + в], где е зависит от шага А/ в методе Т-произведения.
Возвращаясь к уравнению (3.62), запишем его решение в виде
97
|
ф/а)(е°У-1 )+ * |
|
|
|
2 (ay + b) |
I |
^aZ + |
|
ф/а)(еаУ— l)—t |
|
|
+ |
b)W0(x, - i- ln (- 2- Z + l ) d z ) |
||
Группа 5 (3) (г) действует по правилу |
|||
S 0)W |
( « „ , r o)(A:,<,) = |
( 4 ( |
l - ^ ) x |
|
х « о ( * 4 1п (e°!' - - f T)) + |
|||
+ 'F |
( 1+ S у Т ь ) и»{х' 4 ' 1п(в°!' + ~Т т)) + |
|||
|
ф/а){еаУ— 1)+т |
|
|
|
+ -ТЬ5ТГ) |
I |
{‘* + ьЩ х,±1п[+г + |
||
|
Ф/а)(еаУ—Х) — 1 |
|
|
|
+ 04S W |
T Т ) ^ ( х’ 1 Г ы (еаУ- 1 Г * ) ) х |
|||
|
|
+ 2 (ау + |
Ь) ф {еау— 1 ) 4- |
|
+ |
a t + b ) W, (х, ± |
In [еау + |
-5-т)) - |
+ b) W 0 ( т , + - In [еау |
F т ) ) • |
98
3.4. М О Д Е Л И Т Е Х Н О Л О Г И Ч Е С К И Х О Б Ъ Е К Т О В С « П О Д В И Ж Н Ы М » В О З М У Щ А Ю Щ И М В О З Д Е Й С Т В И Е М
Большой класс технологических объектов представ ляют собой станки с двухкоординатными каретками (сверлильные, для сборки печатных плат и др.), а также планшетные координатографы и графопострои тели. Для оптимизации управления этими объектами задаются обычно управляющие воздействия на испол нительные устройства каждой из координат, причем это управление обычно взаимосвязанное, т. е. выбор одного управляющего воздействия, как правило, накладывает ограничение на второе.
Положение платформы в пространстве двух коорди нат (х, у) определяется в любой момент времени соот ношением управляющих воздействий, заданных на каждой координате. Самым распространенным преоб разователем вращательного движения в поступатель ное является шарико-винтовая пара (ШВП). Посколь ку системы управления каретками станков требуют высокого быстродействия, которое обеспечивается оп тимальной структурой управления и допустимыми форсировками, в ШВП возникают большие ускорения. Приводы при этом постоянно работают в динамиче ских режимах с разрывной функцией управления.
Еще более остро стоит вопрос о разрывных функ циях управления, приложенных к ШВП, в коорди натографах или копировальных станках, предназна ченных для работы в динамических режимах. Когда копировальный станок воспроизводит профиль в виде прямого или острого угла, максимальную динамиче скую ошибку, возникающую в вершине угла, можно уменьшить лишь за счет резкого снижения скорости винта, т. е. за счет резкого торможения электродви гателя. Эти форсировки приводят к возникновению в винте различных колебаний: крутильных, продоль ных, поперечных.
99
На ошибку позиционирования каретки непосредст венно влияют продольные колебания, частоту, и ам плитуду которых необходимо оценить. Сформулируем задачу описания продольных колебаний для идеали зированного случая при следующих общепринятых допущениях: 1) считаем платформу абсолютно твер дым телом, жестко связанным с гайкой ШВП, т. е. платформа, перемещаясь вдоль направляющих по основанию станка, повторяет продольные движения гайки; 2) пренебрегаем трением во всех сочленениях станка; 3) пренебрегаем погрешностями изготовления механических узлов; 4) считаем, что опоры винта не допускают продольных перемещений.
Упругий стержень, жестко закрепленный на кон цах, в некотором сечении с координатой х нагружен массой Л4, совершающей вместе с ним продольные крлебания. Таким образом, подвижная платформа с массой М представляет собой для винта «подвиж ное» возмущающее воздействие. Если пренебречь рассеянием энергии, то продольные колебания стерж
ня |
описываются |
уравнением |
|
|
||
|
|
д2у |
Е |
д*у |
|
(3.64) |
|
|
дР ~ ~ ~ |
дх* |
’ |
||
|
|
|
||||
где Е — модуль |
упругости |
материала |
стержня; р — |
|||
его |
плотность. |
|
|
|
|
|
|
Краевые условия |
|
|
|
|
|
|
у (0, /) = |
</(/, 0 = |
0. |
(3.65) |
Выведем условия сопряжения в точке х0. По закону Гука силы упругости, приложенные к платформе (точка х0) со стороны стержня справа и слева, при ближенно равны: Fnpau = ^ле» = SEAylAx, где — площадь поперечного сечения стержня.
100