№ 4 2007
СПИСОК ЛИТЕРАТУРЫ
1. Андреева Л. Е. Упругие элементы приборов. — М.: Машиностроение, 1981. — 392 с.
2. 11 о л и щ у к Д. Ф. Обобщенная теория цилиндрических пружин. — Ижевск: Ихц-во Удм. ун-та. — 216 с.
3. С в е т л и ц к и й В. А. Механика стержней: Учеб. для втузов. В 2-х ч. Ч. 1. Статика. — М.: Высшая школа, 1987. —320 с.
4. С в с т л и ц к и й В. А. Механика стержней: Учеб. для втузов. В 2-х ч. 4.2. Динамика. — М.: Высшая школа, 1987. — 304 с.
5. Расчеты па прочность и машиностроении / С. Д. Пономарев, В. Л. Ьидерман и др.......М.: Машпп, 1959.
Т. 3, —1120 с.
666.1.4:681.7.068.4
ОЦЕНКА СТАТИЧЕСКИХ СВОЙСТВ ОБЪЕКТА УПРАВЛЕНИЯ, СОДЕРЖАЩЕГО УПРУГО-ВЯЗКУЮ СРЕДУ
Канд. техн. наук, доц. НА. ИЛЬИЧЕВ. осп. Л.Б. АЛЕКСЕЕВА, д~р техн. наук, проф. В.П. УВАРОВ
IIр о и е с с в ы т яжк и с • т е к я о из д ел и и (с в е т о в о д о в) из разогр е т о и с т е к л ш i асе ы /; аее. матр ивает с я как объект управления, статические свойства которого определяются на основе интегрирования системы нелинейных дифференциальных уравнений в частных прог гз во д /1 ых. Полу 'чепные решения п озволя/о т о цен ит ь 1 ivecm в ит ел ьно cm ь п р о ц ее а i к во зм\ чцеп и ям в уст сто вившем с я р е. жиме д вижения cm екл ол ¡асе ы.
Process of glassware (optica! paths) extract from a liquid glass is considered as a plant of control which static properties are determined on the basis of integration in a system of nonlinear d iffei 'en t ial part ial ecj nations. Th e / ~ece ived dec is ions alio vv t о e va l и ate pert и r hat ion i 'esp ons ivity of the process in the established conditions of glass maftrftow.
К объектам управления, содержащим упруго-вязкую среду, можно отнести разнообразные технологические процессы изготовления одножильных и многожильных стеклянных стержней и волокон. Это могут быть световоды, полые трубки, тонкие штабики с диаметром поперечного сечения 0,5.,.2 мм. При их изготовлении применяется метод вытягивания из фильеры или заготовки, конец которой размягчен под действием температуры. Размягченное таким образом стекло может вытягиваться под действием собственной силы тяжести, внешних растягивающих усилий или выдавливаться из замкнутого объема.
К стеклянным стержням (световодам), используемым в оптической промышленности, предъявляются высокие требования к стабильности диаметра и формы поперечного сечения по их длине. Поэтому технологический процесс изготовления таких стержней требует постоянного контроля и управления.
Выбор структуры управления, определение требований к измерительным комплексам зависит от свойств объекта управления, который включает в себя как собственно процесс вытяжки, так и устройства, осуществляющие вытяжку.
Решение возникающих в связи с этим задач требует, прежде всего, разработки математических моделей, связывающих входные и выходные параметры процесса. Это позволит оценивать статические свойства объекта управления, определяющие чувствительность выхода процесса к изменению технологических параметров в установившемся режиме. На основе такой оценки можно определить допустимые пределы изменения этих параметров.
№ 4 2007
Оценка динамических свойств связана с изучением переходных режимов, возникающих при управлении процессом вытяжки.
Выделим зону формирования световода (рис. 1), которая представляет собой переход от заготовки (или фильеры) до готового изделия. Па рис. 1 обозначено: г\ г — радиальная и осевая координаты; Ь —длина зоны формирования.
Рис. 1. Зона формирования световода
Геометрия зоны формирования определяется после интегрирования системы дифференциальных уравнений, которая состоит из уравнения неразрывности, уравнений равновесия сил (кинстостатичсскис уравнения) и уравнений, определяющих физическое состояние среды (модель Максвелла) [1]. Выбор этой модели обусловлен тем, что стекломасса в зоне формирования обладает как вязкими, так и упругими свойствами. Это среда, способная к релаксации.
Если считать, что напряжения, скорость, вязкость стекломассы являются функциями только осевой координаты, то получим модель одномерного потока, состоящую из следующих уравнений [ 1 ]:
, 'Л
С1 0 2
(
ОI) оо — + и —
01 дI
\
Р£ +
ор
/
р +к
о р о р —— + и дг дг
3 ц —
(1)
где г — время; и, р —соответственно скорость и напряжение в сечении в
момент
№4
2007
времени /; 5 —площадь рассматриваемого сечения; р, X, ц —соответственно плотность, время релаксации и динамическая вязкость стекломассы в рассматриваемом сечении зоны формирования; £ —ускорение свободного падения.
Граничные условия задачи определяются скоростью вытяжки и и скоростью подачи оп заготовки в зону нагрева
О : и = и ; z~L: и = i>.
(2)
Система уравнений (Г) с граничными условиями (2) описывает движение упруго-вязкой среды с подвижными нижней границей. Поскольку найти аналитическое решение этой системы не представляется возможным, то для необходимо применить ту или иную разностную схему. Но использование разностных методов дает надежный результат только при фиксированных границах. Переход может быть осуществлен с помощью
8 1 д
дифференциального оператора —~- —. Такого рода задача решена в [2]. Однако
dz t ин Заполученное решение не приспособлено для анализа перехода процесса вытяжки к установившемуся режиму. Была осуществлена попытка получить приближенное решение. С этой целью определяли реализуемые диапазоны технологических и геометрических параметров
ив - 1 м/мин. ..10 м/с, о 6 (0,017... 10) [м/с];
D = 10.. .20 мм, D е [0,01... 0,02] [м/с], диаметр заготовки; R е [0,005 ...0,01] [м], радиус заготовки; г = 0,25...0,5 мм;
г е
о е
0,25 • 10 3... 0,5 • 10 3 [м], радиус световода;
7 г > 2 1
, 1
П
; un <0,01ин; g = 9,8 м/с2; р = 2,5 • 10' [кг/м3];
G = 2500 кг/мм2 = 2,5-1010 Па, модуль сдвига; ц е
X -
ю\..ю8
[Па • с]
G 2,5
10"*... 10"
[с].
Было выяснено, что распределение скорости сечения о = и (z) имеет квазилинейный характер. Исходя из этого, в соответствии с ( 1 ) было получено соотношение, определяющее длину зоны формирования в момент времени /,
/ 1 Л
1 _ Г 1 Ï
кК)
V /
о,
Sn
(3)
где К — коэффициент перетяжки, равный К = — = —, К » 1; 50— площадь сечения заготовки. ип ^
Из соотношения (3) можно найти рабочие режимы, соответствующие определенной длине зоны формирования.
Зная длину зоны формирования, можно определить ее геометрию при установившемся режиме. Для этого в системе (1) следует приравнять частные производные по времени нулю.
Второе и третье уравнения полученной таким образом системы могут быть проинтегрированы независимо от первого уравнения. В граничные условия (2) напряжение р
№4
2007
не входит. Поэтому можно свести систему из двух дифференциальных уравнений относительно функций р - /7(2) и о = и (г) к одному дифференциальному уравнению второго порядка вида
с/~\)
dz' Зц-лри~
2Х
dv_ dz
+(l )-Àg)
¿/и
¡77
Это уравнение допускает приближенное аналитическое решение, полученное в [3]. Как показывают расчеты для диапазона выбранных параметров относительная погрешность аналитического решения по отношению к численному не превышает 1%. Кроме того, напряжение оказывается практически постоянным по всей длине зоны формирования.
После выхода на ус тановившийся режим могут происходить отклонения технологических параметров. Эти отклонения вызовут возмущения в зоне формирования световода, что приведет к изменению его геометрических характеристик. Необходимо выяви ть степень влияния такого рода возмущений на геометрию вытягиваемого световода.
Откликом на возмущения, т. с. параметром, определяющим качество световода, является площадь сечения S - S (z, /), где / — момент времени наблюдения. Уравнения возмущений могут быть получены из уравнений движения стекломассы. Решения задачи ищутся в виде
u(z,0 = uU) + í)(z,/) M(ZJ) =m(Z) + M(ZJ) S(z,t)^S(z) + S(zJ)
где o (z), fi(z), S (z) — решения для скорости, вязкости, площади поперечного сечения зоны
формирования, соответствующие установившемуся движению стекломассы; ü(z,/), jl(z,/), S (z, /) — возмущения скорости, вязкости, площади поперечного сечения зоны формирования световода.
После подстановки соотношений (4) в уравнении движения стекломассы для установившегося движения можно получить уравнения возмущений
Ov , сю 1 сю _ сЗи — = -(о + и) -+----U -
dt dz о dz dz
др _ 8г( м) dt л
•uS Su
(4)
СЮ , _ч dp „ dp р 3 м ¿)\J
--(i) + I))-—--I)—-------V--
rte oz oz л л oz
dS - (Зи . _ч S и-— = -(5' + 6)—■-(и + и)----
г; / riz r;z и
oz
(5)
Система (5) допускает приближенное интегрирование при обоснованных практическими исследованиями допущениях. На первом этапе можно проинтегрировать первые два уравнения этой системы, поскольку функция 5 - 5 явным образом в них не входит. С учетом этого получено уравнение связи возмущений по вязкости и скорости
с - \
ai,=aM+(au й-°ОеХР
■ОС,
-ос,.
(6)
OLll Î
где аи = — ; а а, = —; т = -ин; а
и ц К L
времени.
возмущение скорости в начальный момент
№4
2007
После интегрирования третьего уравнения системы (5) можно получить уравнение связи возмущения с откликом
l + av - (1 + а1}) ',
(7)
где а, =318.
Используя соотношение (7), можно нормировать технологические параметры таким образом, чтобы обеспечить заданное качество световода.
Если ввести отклонение не площади ал., а радиуса аг, то (7) можно представить в
виде
ом =
1 + а,
По последнему выражению на рис. 2 представлена связь между безразмерным отклонением радиуса аг и возмущением по скорости и вязкости схи. Как следует из графика, требования к возмущению по скорости и вязкости носят довольно жесткий характер. Так, для того, чтобы отклонение радиуса световода не превышало 6 %, возмущение аи не должно превышать 0,4 %.
I а ■102
7 а „103
Рис. 2. Влияние возмущений па геометрию световода
СПИСОК ЛИТЕРАТУРЫ
1. Уваров В. П., И л ь имев В. А. Математические модели процесса вытяжки оптических стержней. — С.-Петербург, Химиздат, 2003. — 136 с.
2. И л i» и ч с в В. Д. Решение нестационарной задачи движения вязко-упругой ереды с подвижными границами
// Известия вузов. Машиностроение. — 2005. — № 9. — С. 15—19.
3. II л ь и ч с в В. А. Численное и аналитическое решение задачи установившегося движения стекломассы // Известия вузов. Машиностроение. —• 2006. — № 2. — С. 11—14.