CC BY
36
№2
2006
8. Кузнецов В. А. Ползучесть стохастически неоднородных сред в условиях плоского напряженного состояния // Математическая физика. — Куйбышев: КуАИ, 1976. — С. 107—110.
9. С в е ш н и к о в А. А. Прикладные методы теории случайных функций, — М.: Наука, 1968. — 464 с.
10. Закономерности ползучести и длительной прочности: Справочник/ Под ред. С.А, Шестерикова. — М.: Машиностроение, 1983.— 101 с.
Исследуется процесс вытялски топких стержней из разогретой стекломассы.Решается стационарная задача поведения стекломассы в зоне формирования оптических стержней. Задача сводится к интегрированию нелинейного дiiфференциапьного уравнения 2-го порядка относительно скорости. Уравнение не может быть проинтегрировано в замкнутом виде, поэтому, для нахождения решения использовалась система Рун-ге-Кугга.Решение позволяет вводить в полученное дифференциальное соотношение воз-лiyи(ения различных управляющих факторов. Таким образом, лtажио проследить> как меняется с течением времени отклик радиуса сечения волокла в определенном сечении на возмущение вязкости и найти распределение функции отклика радиуса сеченгш по длине волокна в произвольный фиксированный момент времени.
Process of extraction of pin-type rods from fluid glass is examined. Stationary problem of glass melt behavior in a zone of optical rods for mat ion is solved\ It's reduced to integration of a nonlinear differential second-degree equation against velocity The equation cannot be integrated in a self-contained aspect; hence, for finding a solution the Runge-Kugg system was used. Perturbations of various operating factors could be introduced into the obtained differential relation owing to the solution. Thiis, it is possible to trace the response of beam cross-section radius in certain circumstances on viscous resistance and find a distribution of the response function of cut radius throughout the fiber length in the arbitrary fixed instant.
Рассматривается процесс вытяжки оптических стержней из разогретой стекломассы. Качество стержней зависит от множества физических и технологических факторов, оптимальное значение которых обеспечивается системой управления процессом вытяжки [1]. Для решения задач управления в выбранной области параметров необходимо знать статические и динамические свойства объекта управления. Статические свойства определяются чувствительностью процесса к различного рода возмущениям, в том числе и технологических параметров, в установившемся состоянии. Автором ищутся решения для установившегося движения стекломассы в зоне формирования оптического стержня.
В работе [2] показано, что основные характеристики в зоне формирования оптического стержня при установившемся движении могут быть определены для модели одномерного потока из следующей системы дифференциальных уравнений:
666.1.4: 681.7.068.4
ЧИСЛЕННОЕ И АНАЛИТИЧЕСКОЕ РЕШЕНИЯ ЗАДАЧИ УСТАНОВИВШЕГОСЯ ДВИЖЕНИЯ СТЕКЛОМАССЫ
Канд. техн. наук, доц. В.А. ИЛЬИЧЕВ
dz 3|i-v2A,p
dv р-Хрgv
у
d v р-\рgV
№2
2006
где V,/?
формирования
динамическая вязкость стекломассы; р — плотность стекломассы; Х = \х/С — время максвелловской релаксации; С — модуль сдвига; г — аксиальная координата. Рассмотрим интегрирование системы (1) при заданных граничных условиях: г = О,
V - уп , г = Ь, V = где Ь массы в зону формирования
формирования скорость вытяжки.
скорость подачи стекло-
таких граничных условиях известной величиной можно полагать длину зоны фор
мирования, а искомой
функцию
В этом случае систему (1) целесообразнее представить в виде
с1\>
с1г
с1 v с1 р
Зц- V2 Хр
Р-Ь>Р8У
Р(ур-Зця)
(2)
В систему уравнений (2) входит неизвестная функция р = р (г). Чтобы решить задачу с заданными граничными условиями, можно систему (2) интегрировать при граничных условиях вида
г = У
1 Р= Ро
где р0 —значения напряжения в конечном сечении зоны формирования.
Задав р0 и решив систему (2), найдем vn (р0)-
Зададим необходимую погрешность А выполнения граничного условия г = Ь, V
v
(1
Если
V
п (А>)
V
п
< А , будем считать задачу решенной. Если уп (р{)) - уп > А , то итератив-
ным образом будем подбирать р0 до выполнения заданного неравенства уп (р()) - уп
< А.
Ясно, что при таком численном методе встает вопрос о сходимости решения и трудоем-
кости метода. Этих проблем можно избежать, если исключить функцию р = р (г) из (2). Из первого уравнения системы (2) следует:
Р = \>'(А + ВУ2)-С\>,
(3)
где А = В
Хр; С
?Xg\ V
с1\> с1 г
После дифференцирования по £ выражения (3), имеем
р = /(А + В V2) + 2 Въ (у')2 - С у'.
(4)
Из системы (2) находим
где И = р; Е- -Зцр g .
с1 р _ О V р л- Е с1 г
А + Ву
1 *
(5)
Правую часть уравнения (5) можно выразить через функцию V = используя
(4). Тогда (5) примет вид
Р
А + В V
Используя (4), (б), получим
2006
№2
v
//
1
А+ВУ-
А + В\>
+
Е
\
А + ВУ
(7)
У
Возвращаясь к исходным физическим величинам, после преобразований уравнение
(7) примет вид
с1 I1
р
(
Зц-^ру2
2 А,
\
/
\
сЬ_
Л
\г
V
и)
у
¿V
й г
Граничные условия
z = 0, у = v , v = у „ .
в
Заметим, что из уравнения неразрывности следует:
\
8
(8)
У
(9)
5
V = V
Г!
5
14 ^
Í^L = V — ,
II 1 '
(10)
Г
где 5„, г
П 1 Г1
соответственно площадь и радиус поперечного сечения зоны формирования
при г = 0; г — текущий радиус поперечного сечения зоны формирования.
Введем безразмерную функцию
5
х
(1
5
(11)
Тогда при г — 0, д: = 1, г = Ь, х-к, где к — коэффициент перетяжки.
Из соотношений у = V х , у
ух' > у* = у / в (8) можно перейти к искомой функции х
X
н
р
1 1
3 ц- А. р у
/
8
\
у
\
п
(12)
/
Для возможной оценки влияния в дифференциальном уравнении (12) отдельных параметров приведем диапазоны и отдельные значения физических и геометрических факторов: V, = 1м/мин... 10м/с, уе(0,017...10)[м/с]; £> = 10...20мм, Ое [0,01 ...0,02] [м]; Яе [0,005...0,01][м]; г = 0,25мм...0,5мм; ге[0.25-10"\..0,5-Ю"3][м];
у е К. V, ] = V
в
/
Г
Я
\ /
л
; V,, <0,01ув; 8= 9,8м/с:; р = 2,5-105[кг/м3];
й = 2500кг/мм2 = 2,5 -Ю10 Па ;ие
106...10
8
Па-с]; Х = [ю-4...10"2][с].
С 5
1
Введем нормированные величины
цо-106; С = СоЛ010', р = ро-103.
Представим (12) через эти нормированные параметры
х
н
Ро
ЗЦо-^Р()(упЛ-)2-10
X
-7
С
о
X
г
9
Г)
в
/
О
V ю-з о7
\
С
\
V
о
х
8
\
у
II
• ю3
У
(13)
№2
2006
Обращает на себя внимание присутствие в (13) малых слагаемых в правой части урав-
нения: 2 —х(упх')~ х-10 ,
/
в
о
\
С
V
о
* по сравнению с
8
)
уп-10
т , а также
Но
а,
Ро (упхУ ' 1по сравнению с 3 ц0.
Конечно, роль этих слагаемых повышается с увеличением коэффициента перетяжки. Можно выявить диапазон управляющих параметров, при котором (13) можно заменить уравнением
х
Р 8
Зц0
при г = 0, х = 1, г = Ь, х = к
(14)
Решение (14) может быть представлено в виде
х
+
V
к-1
I
Р£
\
L
г + 1.
У
На рис. 1 представлена относительная погрешность 5 аналитического решения уравнения (14) по отношению к численному решению уравнения (13). Расчет велся для параметров ц0 = 1; р0 = 2,5; й0 = 2,5; уп =4,2 -КГ4м-с""1; £=400; L = 0,3м . В рассмотренном случае погрешность не превышала 0,3 %.
Таким образом, в зависимости от параметров процесса можно выбрать метод решения задачи.
Полученные решения стационарной задачи движения стекломассы в зоне формирования оптического стержня соответствует стабильным условиям протекания процесса
вытяжки. Если к параметрам, соответствующим невозмущенному движению стекломассы г, V , [I, добавить отклонения первого по-
г, м
Рис. 1. Погрешность приближенного решения
рядка этих величин /\ V, ц, то можно получить уравнения, связывающие эти отклонения. Считая одни из отклонений возмущениями, а другие откликами на эти возмущения, можно оценить чувствительность зоны формирования к изменению технологических параметров, а следовательно, опреде-
лить точность поддержания этих параметров на заданном уровне.
СПИСОК ЛИТЕРАТУРЫ
1.
Уваров В, П., И л ь и ч е в В. А. Математические модели процесса вытяжки оптических стержней. СПб.: ХИМИЗДАТ, 2003. — 136 с.
Ильичев В. Ам У в а р о в В. П. Определение формы зоны вытягивания вязкой стекломассы в стационарном режиме // Проблемы машиноведения и машиностроения: Межвуз. сб. — СПб.: СЗТУ, 2002. — Вып. 25. — С. 170—173.
