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

книги / Микропроцессорное управление технологическими процессами в радиоэлектронике

..pdf
Скачиваний:
0
Добавлен:
12.11.2023
Размер:
10.3 Mб
Скачать

Следует отметить, что В — нелинейный оператор и ему соответствует нелинейная полугруппа, тем не менее метод Г-произведений можно применять в этом случае [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

Соседние файлы в папке книги