Файл: Решение задач машиноведения на вычислительных машинах [сборник]..pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 21.10.2024
Просмотров: 92
Скачиваний: 0
Разработанные методы расчета летучих ножниц с кривошипнокоромысловым механизмом резания обеспечивают возможность создания справочного атласа номограмм и графиков, отражающих кинематические и динамические характеристики. Использование справочных данных позволит непосредственно довести методы синтеза до инженерной практики.
Некоторые образцы справочных карт приведены на рис. 3 и 4. На рис. 3 представлены кривые экстремальных значений i, V, к в зависимости от длины коромысла L4 для симметричного режу щего механизма при Z/3= 1300 мм, ср3|ш=190о, ср4й)1=90°. На рис. 4 даны графики изменения вертикального перекрытия в зависимости от длины коромысла L4 для различных величин Ьъ при ^p3l{M= 190°,
D = 0, Дм=0,2 мм.
Пользуясь приведенными данными, конструктор, в зависимости от технологического назначения ножниц, может выбрать параметры режущего механизма, обеспечивающие при заданном боковом зазоре между ножами нужную величину вертикального перекры тия, с учетом инерционных нагрузок, воздействующих как на звенья самого механизма, так и на привод ножниц.
Расчет кривошипно-коромыслового механизма резания ле тучих ножниц проводился на ЭЦВМ «Минск-32». Алгоритм расчета написан на алгоритмическом языке Алгол-60. Использование вычислительной техники позволило всесторонне изучить особен ности исследуемого механизма и составить справочные данные для конструкторов.
ЛИ Т Е Р А Т У Р А
1.Н. И. Крылов, С. Н. Сужений. Определение максимально возможного
перекрытия ножен летучих ножниц. — Труды ВНИИМЕТМАШ, 1967,
сб. 19.
2.С. Н. Сужений. Новая методика определения зазоров между ножами летучих ножниц. — НИИИНФОРМТЯЖМАШ. Металлургическое обо рудование. Оборудование для прокатного производства. 1-68-21. М., 1969.
3.Н. И. Левитсний. Применение электронных цифровых машин для не которых задач анализа и синтеза четырехзвенных шарнирных механиз мов. — Сб. «Анализ и синтез механизмов». М., ГНТИ машиностроитель ной литературы, 1963.
4.И. И. Артоболевский. Теория механизмов. М., «Наука», 1965,
84
АЛГОРИТМ НАБЛЮДЕНИЯ ПАРАМЕТРОВ И ФАЗОВОГО СОСТОЯНИЯ НЕЛИНЕЙНЫХ ОБЪЕКТОВ
А. И. Медник
Постановка и решение задачи о наблюдении линейного объекта по неполной информации известны [1—5]. Задача о наблюдении параметров и фазового состояния нелинейного объекта состоит в следующем. Известны его динамические свойства, выражающиеся уравнением
~/(а*, с, t), |
{/?”}; с£{Я*}; |
с = const; t > |
10. (1) |
Начальные данные |
(/0) и параметры |
с . неизвестны, |
однако |
на траекториях (1) измеряются величины, составляющие /-вектор информации]
y(t) = h(x, с, /), y £ { R 1}; /<> -(-& . |
(2) |
Уравнения (2) неразрешимы относительно х, с. Требуется опре делить параметры с . (/ = 1, . . ., г) и восстановить движение,
т. е. найти х (/), используя (1), (2) в некотором промежутке [а, (3] предыстории процесса. Трудность решения задачи в случае не линейности системы (1) состоит в том, что здесь невозможно по строить обратную разрешающую операцию, как это может быть сделано для линейных систем. Эта трудность преодолена [6, 7] путем построения операции-процедуры, включающей решение некоторой экстремальной задачи.
Интегральные операции наблюдения задаются построением наблюдаемого сигнала Y и производящей А^-вектор-функции Ф, определяющей процедуру наблюдения. Используя оператор инте грирования
x(t) = g(xu, с, г0, г), г > го |
(3) |
и предполагая обычные свойства интегрируемости и дифференци руемости по начальным данным и параметрам в (2), (3), можем построить следующую операцию наблюдения неизвестных пара метров и фазового состояния.
Сигнал |
11|)ОИЗНОДЯ1ЦЛЯ функция |
||
Y0(t) = |
y(t) |
Ф0(С, с, t) = h(Z, |
с, t) |
|
13 |
1„ |
(4) |
Y, (О = |
У5СО d-z ф ,С»> с, 0 = 5 |
с> *> 0>ХЫ Х |
|
/, = /+ <)«,, О > 0 , |
0 < а , < 1 ■ ( « = ! , . . . , Р), |
85
Объединив векторы с, с в единый тг+/с-вектор z, рассмотрим си стему N =1 (р + 1) уравнений с п-\-к неизвестными
Ф(*, t )=Y(t ), Y£{R*}; z £ { R n'k), N ^ n + h. |
(5) |
Система (5) всегда имеет решение, которым является совокупность истинных значений параметров с . и фактически реализовавшегося
в (1) состояния х (t). Однако функция Ф (z, t) существует лишь в виде процедуры, реализуемой на ЦВМ, а множество решений (5) может быть конечным, счетным или континуальным. Эти обстоя тельства приводят к необходимости построения и решения экстре мальной задачи вместо уравнения (5).
Выбрав подходящую норму вектора и задав замкнутое мно жество Z (t), предположительно содержащее искомые значения параметров с. и фазовых переменных x i (t), сформулируем задачу
математического программирования. |
|
F(z, *) = ЦФ(*, t ) - Y ( t ) I inin, z £ Z ( t ) Q { R nhk). |
(6) |
z |
|
Процедура наблюдения параметров с . и фазового состояния х (t)
состоит из следующих этапов: формирование наблюдаемого сигнала и производящей функции, решение задачи (6) и поиск среди «ну левых» {F (z, 0 —0} решений фактических параметров с . и фазо-
вых переменных xi (I), интегрирование (1) при х0=х (t).
В таком виде, с однократным решением экстремальной задачи, процедура представляет собой систему прогнозирования. Если переставить в (4) пределы интегрирования, задав 0 < 0, и вести непрерывное решение последовательности экстремальных задач, то получается следящая система, пригодная для использования
вуправляющем устройстве. Наличие в (1) управляющего воз действия существенных осложнений не вносит.
Достаточное условие наблюдаемости параметров и фазового состояния нелинейного объекта, по аналогии с [7], может быть сформулировано следующим образом.
Теорема. Производящая функция Ф (z, t) обеспечивает полную наблюдаемость параметров и фазового состояния системы (1)—(2)
взамкнутой окрестности 0 точки z0 в момент t, если ранг ее ма трицы Якоби №/dz в точке z0 равен п-\-к. Здесь п — порядок системы (1), к — число параметров Cj.
Локальное свойство, сформулированное в теореме, естествен ным образом распространяется на замкнутую область Z (t) в силу ограниченности множества нулевых решений в {Rn+k}Z)Z (t) (если множество нулевых решений задачи (6) бесконечно в Z Ц), система ненаблюдаема [7 ]) и гарантирует реализуемость вычисли тельной процедуры наблюдения в области Z (t).
Разработана АЛГОЛ-программа [8] решения задачи о наблю дении параметров и фазового состояния, основными блоками ко-
86
торой являются следующие (см. рисунок): формирование (1) наблюдаемого сигнала и (2) производящей функции (процедура наблюдение); решение (3) оптимизационной задачи (6); поиск (4) среди «нулевых» {F (z, t)=0} решений фактических параметров и фазовых переменных с•/., (t) (процедура проверка истинности
наблюдения); интегрирование (5) системы (1) при Х 0 «3/ (0 ДО за данного момента Т > t (процедура ДИФУР).
1_
Кроме этих блоков, программа содержит стандартную проце дуру (6) интегрирования уравнений (1) методом Рунге-Кутта, реализующую оператор (3); блок (7) моделирования измеритель ного устройства (2) для обеспечения возможности проверки системы наблюдения до ее фактического включения в объект; блок (8) формирования минимизируемого функционала из (6) (процедура функционал); оптимизатор (9), построенный на основе случайного поиска (процедура поиск); блок (10) проверки достаточных усло вий полной наблюдаемости (процедура Якоби).
Структура алгоритма предусматривает помимо решения задачи о наблюдении параметров и фазового состояния системы проведе ние некоторых исследований с помощью ЦВМ, необходимых для построения эффективной системы управления или наблюдения. В частности, с помощью блока (7), моделирующего измерительное устройство (2), может быть произведена оценка эффективности различных систем измерения и выбрана на основе машинного экспе
рт
римента такая конструкция измерительного устройства (и соот ветственно Z-вектор измерений y=h (х , t), при этом I — не фикси ровано), которая обеспечит наилучшую сходимость вычислитель
ного процесса. Эффективность процедуры |
наблюдения зависит |
от выбора параметров процедуры 0, as (s= l, |
. . .,/?). Как показано |
в [6], возможно даже вырождение системы наблюдения вследствие неудачного выбора моментов измерений. Поэтому важно, еще до на чала функционирования системы наблюдения объекта, выявить область допустимых значений параметров процедуры б, а8, обеспе чивающих хорошую сходимость оптимизационного процесса. Эта задача также решается путем проведения машинного экспери мента.
П р и м е р . Рассмотрим задачу о наблюдении гирокомпаса с нелиней ной восстанавливающей силой [5]. Дифференциальные уравнения прецес сионного движения гирокомпаса можно представить в виде
dxi А |
А |
dt = ~В + ~В (1 — р) |
|
dx> |
(7) |
~dl~ — —Axi — |
dx%
dt — P (x24* хз)■
Л = 1,5392; 5=41,1368; |
p = 0,38; |
F = l,5 - 1 0 - 3; i]=0,4 • 10“«. |
|
|||||||||
Измеряется |
фазовая координата хг, т. |
е. |
|
|
|
|
|
|
|
|||
|
|
|
У (t)=x1, |
(t |
f0); |
|
|
|
|
|
(8) |
|
|
|
£о= 0 ; 0=600; |
а3={0,25; |
0,5; |
0,75; |
1,0}. |
|
|
|
|||
Необходимо |
определить |
начальное состояние |
х |
(t0). |
|
|
|
|||||
В качестве измерительного устройства была использована процедура |
||||||||||||
интегрирования |
системы (7) |
с |
начальными |
условиями |
[5]: |
а:1(^0)= |
0,3; |
|||||
^2(^о)=а:з(^о)= ^ПО _3. Множество |
Z(t0) |
задавалось |
в |
виде: 0,1 |
< zx < |
0,6; |
||||||
10_3 < z2 < |
10~2; 10-3 ^ z3 < |
10 ~2. В результате работы программы получены |
||||||||||
следующие |
значения начальных |
условий: |
д:1(^0)= 0,299; |
х2 (£0)= 3,98-10“3; |
||||||||
х3(^0)= 4 ,0 4 -1 0 _3. |
|
|
|
|
в реальных системах управления |
|||||||
Программа может быть использована |
и наблюдения; для проверки и оценки эффективности различных измери тельных устройств, входящих в систему управления; для решения задач идентификации нелинейных систем. Предусмотрена возможность иденти фикации систем с неизвестными параметрами — функциями cj (t), допускаю щими разложение по некоторой конечной системе известных функций.
ЛИ Т Е Р А Т У Р А
1.Р. Е. Калмап. Об общей теории систем управления. — Труды I Конгресса ИФАК, 2. М., Изд. АН СССР, 1961.
2.Н. Н. Красовский. Проблемы управляемости, наблюдаемости и стабилизируемости динамических систем. — Труды II Всесоюзного съезда по
теоретической прикладной механике. М., «Наука», |
1965, |
вып. 7. |
3. Н. Н. Красовский. Теория управления движением. |
М., |
«Наука», 1968. |
88
4. Э. Г. Альбрехт, II. |
I I . Красовский, О наблюдении нелинейной |
управля |
емой системы в окрестности заданного движения. — Автоматика |
и телеме |
|
ханика. 1964, 25, № |
7. |
|
5.Я . Н. Ройтенберг, Некоторые задачи управления движением. М., Физматгиз, 1963.
6. Е. А. Гальперин. Процедуры наблюдения нелинейных систем. — Изв. АН СССР, Техническая кибернетика, 1972, № 1.
7.Е. А. Гальперин. О наблюдаемости нелинейных систем. — Изв. АН СССР, Техническая кибернетика, 1972, № 2.
8.С. С. Лавров. Универсальный язык программирования. М., «Наука», 1972.
МОДЕЛИРОВАНИЕ НА ЭЦВМ ВЛИЯНИЯ РАЗБРОСА ПАРАМЕТРОВ
НА УСТОЙЧИВОСТЬ ДИНАМИЧЕСКОЙ СИСТЕМЫ МЕТАЛЛОРЕЖУЩИХ СТАНКОВ
Г. И. Фирсов
Изучению динамического поведения металлорежущих станков при условии детерминированности значений параметров динами ческой системы и отсутствии возмущающих воздействий, нося щих случайный характер, посвящено большое количество работ
[1, 2, 3].
Опыт эксплуатации, а также выполненные в последние годы исследования [4, 5] динамики металлорежущих станков с учетом стохастического характера внешних воздействий и разброса пара метров убедительно доказывают, что наиболее достоверная оценка динамических процессов в станках может быть получена лишь с теоретико-вероятностных позиций. Исходные данные для расчетов динамического поведения станков (экспериментальные или расчетные) отличаются, как правило, от фактических. В част ности, жесткостные параметры металлорежущих станков значи тельно меняются в зависимости от сборки, пригонки сопрягае мых поверхностей, от зазоров и натягов в сочленениях. Напри мер, при обследовании жесткости 150 токарных станков с высотой центров 200 мм, в заводских условиях установлено, что жест кость находилась в пределах от 1000 до 5000 кГ1мм, в то время как средняя ее величина составляла 1640 кГ1мм [3].
Для 50 станков 1Е61М повышенной точности с диаметром об рабатываемого изделия 320 мм были определены законы рас пределения жесткости передних бабок. По результатам испыта нии построены гистограмма и выравнивающая ее кривая плот ности нормального распределения (рис. 1).
где ] — жесткость в кПмкм [6].
7 Р еш ен и е за д а ч |
89 |