Файл: Севостьянов, А. Г. Основы математического моделирования механико-технологических процессов текстильной промышленности из цикла лекций заочного факультета по технической кибернетике в текстильной промышленности.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 31.10.2024
Просмотров: 69
Скачиваний: 0
Рассмотрим задачу при следующих численных значениях па раметров:
т*=г8 с; т0 = 15 с; 70 = 20 с; /), = 1600 [г/м2]2.
Значения параметров, имеющих размерность времени, таковы, что моделирование можно вести в натуральном масштабе времени. Однако имеет смысл изменить масштаб времени, с тем чтобы ус корить процесс вычисления R2(t i ) . Поэтому машинную постоян ную времени Гм чесальной машины примем равной 1 с; тогда Mt =
— — . Соответственно изменится и |
т.: |
т ки — —-------- = 0,4 с. |
|
Тм |
к |
км |
20 |
Величина запаздывания в объекте т0 тоже изменится, но пос кольку это запаздывание, как говорилось выше, можно исключить при расчете R2(ri), значение том не представляет интереса.
Для величины Di можно принять любой |
удобный |
масштаб. |
|
Учитывая |
характер изменения зависимости |
(133) и |
вид схемы |
(см. рис. |
18, в), удобно величину Dx моделировать возможно боль |
шим допустимым для АВМ напряжением [напряжение в любой точке схемы будет ниже, чем выбранное значение «о]- Примем, к примеру, «о=80 В. Масштаб в этом случае равен:
М, _ А . |
1600 |
|
г 2 |
80 |
= |
20 |
|
щ |
|
|
|
Значения R2(ti) будут получаться |
на выходе схемы в машин |
ных переменных. Для фиксированного т, величину ^?2(тг) имеем на выходе интегратора, строго говоря, при т -*■счз. Однако практи чески постоянное значение, напряжения на выходе интегратора ус тановится уже при t > (5 ч -1 0 ) (тйм + ТЬи) . Найденные значе ния R2{Ti) пересчитываются в реальный масштаб величин с по мощью масштабов Mt и jM d . При использовании последнего сле дует иметь в виду, что имеющийся в схеме блок перемножения ‘осуществляет операцию умножения двух напряжений vx и ь2 по закону м = 0,01 v xv2- Отсюда, с учетом коэффициента £„ = 0,01,
М п = М п. 100=2000 — — . |
■ |
||
° |
п |
м2В |
|
Легко видеть также, что дисперсию чесальной ленты можно |
|||
рассчитать с помощью модели, |
показанной на |
рис. 18, в, приняв |
тг= 0, т. ,е. исключив из схемы блок запаздывания БПЗ.
Схема набора модели на ABM МН-7М приведена на рис. 19, б. Здесь операционные усилители 1, 2 и 3 образуют апериодическое звено с постоянной времени т;(М= 0,4 с; начальное напряжение на выходе интегратора составляет 80 В. При включении схемы на вы
ходе интегратора 3 напряжение будет,меняться по закону |
/?1 „(т)= |
= 8 О ~ т/0'4;тем самым моделируется корреляционная функция RДт). |
|
Операционный усилитель 4 моделирует апериодическое |
звено с |
постоянной времени 70м=1 с, т. е. моделирует объект. Блоки БП и Б П З — это стандартные блоки перемножения и запаздывания.
50
Усилитель 5 реализует операцию интегрирования. На выходе это го ОУПТ в машинных переменных получается значение Яг(хг) при t > (5 -г Ю) (iftM+ Г0м) для любого заданного на блоке-
запаздывания хи Применение масштабирования по времени позволило в данном
примере не только ускорить вычисление значений Яг{^г)), но и уп ростить подбор сопротивлений в схеме благодаря удобному вы бору машинной постоянной времени 7ом=1 с.
На рис. 20 точками обозначены значения ^г(т), полученные на модели. Плавная кривая соответствует теоретически рассчи танной кривой ^?г(т). Из рисунка видно хорошее совпадение ре зультатов, найденных с помощью модели, с расчетными значе ниями.
§ 4. М О Д Е Л И Р О В А Н И Е ПРОЦЕССА В Ы ТЯ ГИ В А Н И Я В В Ы Т Я Ж Н Ы Х П РИ БО РАХ
Вытяжной прибор (ВП) является весьма распространенным элементом машин прядильного производства, и поэтому исследо вание его работы имеет большое значение.
Основными функциями ВП являются утонение продукта (лен ты, ровницы) и распрямление волокон в продукте. Процесс уто нения лучше изучен. Исследование же второй функции представ ляет огромные трудности. Из-за вероятностного характера изме нения толщины продукта при рассмотрении работы ВП требует-' ся обработка статистического материала большого объема. Необ ходимость создания систем автовыравнив-ания поставила вопрос об упрощенном, но адекватном в существенных чертах описании динамики ВП.
51
С этой целью было разработано значительное число более или менее сходных математических моделей, описывающих динамику ВП. Как правило, изучались два канала динамической системы ВП: 1) от развеса продукта на входе ВП к развесу продукта на выходе ВП (передаточную функцию обозначим здесь через Н70); 2) от скорости питающей (или вытяжной) пары к развесу продук та на выходе ВП.
Ниже рассматривается первый канал.
Предположим, что все волокна имеют одинаковую длину L и полностью распрямлены. Пусть, кроме того, плотность передних концов волокон в ленте описывается функцией ni(x). Тогда мас са волокон, передние концы которых находятся на некотором участке 0—х, определится по формуле
s ( x ) = а Ц и ( х — X;) —- и ( х — х / — /.)] ~ |
|
^ п ( х ' ) W(x — x')dx'. |
(134) |
b |
|
Здесь и(х) — так называемая единичная ступенька, или функция Хевисайда [2]. Переходя в этом выражении к изображениям по Лапласу, получим:
sip) = П^р) Щр) = пг(р) 1 е —. |
(135) |
Р
Изображение развеса как функции, производной по х от мас сы s(x), имеет вид
g {p)^pW (p)nl(p). |
(136) |
Остановимся теперь на простейшей схеме работы ВП. Многочисленными экспериментами установлено, что волокна в
ВП можно разделить на две группы: контролируемые, меняющие свою скорость с V\ на v2 в одной и той же точке неконтролируе мые, короткие, меняющие скорость в различных точках ВП. Нали чие второй группы волокон приводит к появлению дополнительной неровноты в процессе вытягивания. Полное число волокон в ВП в момент времени t
n(t) = nK(t) + nH(t). |
(137) |
Предположим, что число неконтролируемых волокон пропор ционально развесу продукта на выходе ВП и длине волокон:
. . . kl. ds(t)
л„( О = —
dt
Коэффициент k выбирается исходя из дополнительных сообра жений о доле неконтролируемых волокон в ВП или из условия близости экспериментальных и теоретических результатов.
52
Обозначим через /(|) вероятность для одного волокна изме нить свою скорость с о на 02 в области точки Тогда полное чис ло неконтролируемых волокон в ВП в момент t
, , ч |
f kL |
ds(v,t) и ,, |
n a( t ) = |
— • — ^ f ( v 2t ) d t ; |
|
|
Оv2 |
dt |
здесь использовалось соотношение x = v 2t.
После преобразования этой формулы по Лапласу получаем:
ы
А'и(p) = - p S {p )F (p ) . (138)
о2
Наличие неконтролируемых волокон приводит и к изменению числа волокон, имеющих скорость щ. Чем большая доля неконт ролируемых волокон в данный момент относится к волокнам, име ющим скорость Ог, тем меньше доля этих волокон среди волокон, имеющих скорость о,. Учет этого фактора и уравнений (135), (137) и (138) дает возможность построить блок-схему процесса вытягивания, которая изображена на рис. 21, а. В отрицательную
обратную связь введена вытяжка |
Е — |
предусматривающая, |
что волокна на входе в ВП движутся со скоростью vt -до точки перехода
sa
Предполагаем, |
что |
f( |) |
|
имеет экспоненциальную |
форму: |
|||
/(?) = е хр | — у |
J. |
а = |
10 |
(по условию, что 7(E) |
существенно |
|||
отлично от 0 на отрезке |
0 < |
х |
< -^ )• Тогда |
|
|
|||
|
|
|
№ = |
|
1 |
|
(139) |
|
|
|
|
1 + р |
L |
|
|||
|
|
|
|
|
10 v 2 |
|
|
|
|
|
|
|
|
|
|
|
|
Преобразуя |
выражение |
|
W(p) = а 1— ехр(— p L /v s) |
с |
использо- |
|||
|
|
|
|
|
|
Р |
|
|
ванием разложения в ряд Паде, получаем:
W(p) = а |
Шг |
(140) |
|
l-\-p L /2 v a |
|
^ Здесь линейный размер L переведен во временной интервал Т——. |
||
|
|
Ъ'3 |
На основании выражений (139) и (140) строим структурную |
||
схему модели процесса вытягивания (рис. 21, б). |
Модель эта, |
разумеется, весьма приближенная. В ней не учитывается ни один из параметров ВП, за исключением скоростей Щ и v2 (именно поэтому исследование второго канала динамической системы на указанной модели затруднено).
Показанная на рис. 21,6 модель позволяет выявить роль не контролируемых волокон в процессе вытягивания, их влияние на неровноту продукта. Модель обладает свойством линейности, поэ тому на ней можно изучать зависимость неровноты выходящего продукта от неровноты входящего и от доли неконтролируемых волокон отдельно. Кроме того, поскольку в постоянные времени Т{ элементов модели входит средняя длина волокна L, можно исследо вать ее влияние на генерирование’дополнительной неровноты ВП.
§ 5. М АТЕ М А ТИ Ч Е С К О Е М О Д Е Л И Р О В А Н И Е ПРОЦЕССА С М Е Ш И В А Н И Я ВОЛОКОН
ПРИ В Ы Т Я ГИ В А Н И И П РО Д УКТА В В Ы Т Я Ж Н О М ПРИ БО РЕ
Материал данного параграфа является как бы развитием ма териала предыдущего: на основе модели процесса вытягивания рас сматривается процесс смешивания двух видов продукта, например двух лент, в ВП.
Две ленты А я В подаются на вход ВП и вытягиваются вместе. Во время вытягивания лент волокна взаимодействуют между со бой. На выходе ВП образуется одна лента, неровнота которой за висит от неровноты лент на входе, от работы ВП, от доли неконт ролируемых волокон и от степени взаимодействия волокон обеих лент.
54
Построим модель этого процесса.
Число передних концов волокон на входе ВП n ^ t) равно сум ме числа передних концов волокон в лентах А и В:
«i(0 = n lA(t) Н- niB(t).
При сложении возможны два случая:
1) когда взаимодействием волокон разных компонентов можно пренебречь. Тогда каждый компонент вытягивается независимо от другого, а толщина продукта на выходе ВП равна сумме тол щин каждого из компонентов на выходе из ВП. Процесс проис ходит так, как если бы ленты вытягивались на разных ВП, а по том складывались;
2) когда волокна обоих компонентов взаимодействуют в ВП, что проявляется в проникании волокон каждой из лент в другую ленту. При этом, однако, вытягивание нельзя рассматривать как вытягивание однородного продукта, так как ленты А и В вытя гиваются по-разному, в зависимости <?т средней длины волокон в компоненте и доли неконтролируемых волокон.
Блок-схема работы ВП для первого случая представлена на рис. 22, а.
Чтобы построить модель для второго случая, необходимо учесть взаимодействие волокон разных компонентов. Воспользуем ся здесь простейшей моделью: будем считать, что доля неконтро лируемых волокон каждого из компонентов пропорциональна раз весу всего продукта на выходе ВП (в действительности имеет место и вариант, когда доля неконтролирумых волокон компонен та пропорциональна развесу последнего на выходе ВП, что соот ветствует вытягиванию по описанной выше схеме).
55‘
Обозначим через а А и а в |
доли волокон каждого из компонен |
тов, вытягиваемых по схеме, |
представленной на рис. 22,а, а через |
Ра и Рв — доли волокон, вытягиваемых по схеме, когда коли чество неконтролируемых волокон зависит от развеса всего выхо дящего продукта. Очевидно, должны выполняться два условия:
ал-Ь«в = Рл+Рв— 1-
Используя блок-схему процесса вытягивания (см. рис. 21, а), построим блок-схему модели вытягивания двухкомпонентного продукта (рис. 22, б). Соответствующая структурная схема набо ра модели на АВМ представлена на рис. 23. Эта модель позволя-
Р йс. 23
ет провести широкий круг исследований по выявлению роли раз личных факторов на неровноту выходящего продукта. Влияние каждого из вариантов смешивания при вытягивании можно изу чить,* варьируя коэффициенты а и р . Варьируя же коэффициенты КА и Кв, можно выявить влияние неконтролируемых волокон на процесс смешивания. Наконец, варьируя средние длины волокон каждого из компонентов (ЬА и LB), можно проанализировать эф фективность смешивания в зависимости от длины волокон компо нентов.
56