Научная статья на тему 'Применение квадратичных функций Ляпунова при решении прикладных динамических задач'

Применение квадратичных функций Ляпунова при решении прикладных динамических задач Текст научной статьи по специальности «Математика»

CC BY
170
10
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
НЕПРЕРЫВНАЯ ДИНАМИЧЕСКАЯ СИСТЕМА / CONTINUOUS DYNAMIC SYSTEM / МАКРОСТРУКТУРА ПРОСТРАНСТВА СОСТОЯНИЙ / MACROSTRUCTURE OF THE STATE SPACE / МЕТОД ФУНКЦИЙ ЛЯПУНОВА / LYAPUNOV FUNCTION METHOD

Аннотация научной статьи по математике, автор научной работы — Антоновская О.Г.

В настоящей статье предлагается методика определения момента окончания переходного процесса при изменении параметров непрерывной динамической системы на основе применения прямого метода Ляпунова, основанная на оценивании областей притяжения состояний равновесия. За оценку области, которую траектория системы не покинет с течением времени, принимается окрестность состояния равновесия, в которой первая производная функции Ляпунова отрицательна. Функция Ляпунова выбирается в виде положительно определенной квадратичной формы, для которой знакоотрицательность ее первой производной в силу линеаризованной системы обеспечивается с заданным запасом.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Антоновская О.Г.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

APPLICATION OF LYAPUNOV'S QUADRATIC FUNCTIONS IN SOLVING APPLIED DYNAMIC PROBLEMS

In the following paper, we propose a technique for determining the moment of the end of the transient process when the parameters of a continuous dynamic system are changed on the basis of the application of the direct Lyapunov method, based on the estimation of the regions of attraction of equilibrium states. For an estimate of the region the trajectory of the system does not leave with time, we take a neighborhood of the equilibrium state where the first derivative of the Lyapunov function is negative. The Lyapunov function is chosen in the form of a positive definite quadratic form, where the negative sign of its first derivative is provided with a given margin due to a linear system.

Текст научной работы на тему «Применение квадратичных функций Ляпунова при решении прикладных динамических задач»

2. Goryunov V.I. K analizu bystrodeystviya sintezatora pri pereklucheniyah po diapazonu [On synthesizer speed analysis when switching in a range] / V .I. Goryunov, V. N. Eruslanov, M. N. Zaytseva // Tehnika sredstv svyasi. Ser. Tehnika radiosvyasi [Techniques of communication. Series Techniques of radio-communication]. - 1986. Vyp. 3. P. 44-49. [in Russian]

3. Antonovskaya O. G. K analizu formy I dlitelnosti perehodnyh processov pri pereklyucheniyah sintezatora s delitelem chastity I proportsionalno-integrtiruyushim filtrom po diapazonu [On the analysis of form and duration of transient processes when switching over a range in the synthesizer with counter and proportional-integrating filter] / O. G. Antonovskaya, V. I. Goryunov, N. I. Lobashov // Dinamika system: Mezhvuz. sbornik [System dynamics: Inter-university Collection]. - Gorky: Izd-vo GGU, 1989. P. 59-72. [in Russian]

4. Lyapunov A. M. Obshchaya zadacha ob ustoychivosti dvizheniya [General problem of the stability of motion]. / A. M. Lyapunov. - Moscow-Leningrad : Izdat. Nthn.-teor. Lit., 1950. 472 p. [in Russian]

5. Barbashin E.A. Funktsii Lyapunova [Lyapunov functions]. /E. A. Barbashin. - Moscow: Nauka, 1970. 240 p. [in Russian]

6. Chetaev N. G. Ustoychivost' dvizheniya [Stability of motion]. / N. G. Chetaev. - Moscow: Nauka, 1966. 207 p. [in Russian]

7. Antonovskaya O. G. On the construction of quadratic Lyapunov function with given properties / O. G. Antonovskaya // Differential equations. - 2013. - V. 49. - № 9. - P. 1187-1191.

8. Antonovskaya O. G. Determination of coefficients of a quadratic Lyapunov function with given properties / O. G. Antonovskaya // Differential equations. - 2016. - V. 52. - № 3. - P. 275-261.

9. Sarybekov R.A. Ekstremal'nye kvadratichnye funktsii Lyapunova system uravneniy vtorogo poryadka [Extremal quadratic Lyapounov functions for systems of second order] / R.A. Sarybekov // Sib. Matem. Zhyrn. [Siberian Mathematical Journal] - 1977. - V. 18. - № 5. - P. 1159-1167.

10. Antonovskaya O. G. O postroenii I primenenii uslovno-extremalnoy funktsii Lyapunova pri izuchenii ddinamiki sistemy s pomoshyu tochechnogo otobrazheniya ploskosti v ploskost' [On a construction and use of extreme-conditional Lyapunov function when studying system dynamics by poinj mapping plane to plane] / O. G. Antonovskaya // Matematicheskoe modelirovanie I optimal'noe upravlenie [Mathematical modeling and optimal control] Vestnik NNGU. N. Novgorod: Izd-vo NNGU. - 2006. - № 3(32). - P. 110-117. [in Russian]

DOI: https://doi.org/10.23670/IRJ.2017.62.066 Антоновская О.Г.

ORCID: 0000-0002-5688-7996, Кандидат физ.-мат. наук, Нижегородский государственный архитектурно-строительный университет ПРИМЕНЕНИЕ КВАДРАТИЧНЫХ ФУНКЦИЙ ЛЯПУНОВА ПРИ РЕШЕНИИ ПРИКЛАДНЫХ

ДИНАМИЧЕСКИХ ЗАДАЧ

Аннотация

В настоящей статье предлагается методика определения момента окончания переходного процесса при изменении параметров непрерывной динамической системы на основе применения прямого метода Ляпунова, основанная на оценивании областей притяжения состояний равновесия. За оценку области, которую траектория системы не покинет с течением времени, принимается окрестность состояния равновесия, в которой первая производная функции Ляпунова отрицательна. Функция Ляпунова выбирается в виде положительно определенной квадратичной формы, для которой знакоотрицательность ее первой производной в силу линеаризованной системы обеспечивается с заданным запасом.

Ключевые слова: непрерывная динамическая система, макроструктура пространства состояний, метод функций Ляпунова.

Antonovskaya O.G.

ORCID: 0000-0002-5688-7996, PhD in Physics and Mathematics, Nizhny Novgorod State University of Architecture and Civil Engineering APPLICATION OF LYAPUNOV'S QUADRATIC FUNCTIONS IN SOLVING APPLIED DYNAMIC

PROBLEMS

Abstract

In the following paper, we propose a technique for determining the moment of the end of the transient process when the parameters of a continuous dynamic system are changed on the basis of the application of the direct Lyapunov method, based on the estimation of the regions of attraction of equilibrium states. For an estimate of the region the trajectory of the system does not leave with time, we take a neighborhood of the equilibrium state where the first derivative of the Lyapunov function is negative. The Lyapunov function is chosen in the form of a positive definite quadratic form, where the negative sign of its first derivative is provided with a given margin due to a linear system.

Keywords: continuous dynamic system, macrostructure of the state space, the Lyapunov function method.

При исследовании динамики систем одной из наиболее важных является проблема определения длительности переходных процессов в системе. В ходе решения задачи определения длительности переходных процессов в существенно нелинейных системах при переключении параметров в заданном диапазоне возникает задача точного определения момента окончания переходного процесса в системе [1], [2], [3]. Дело в том, что критерием его окончания может оказаться попадание траектории на заданное, быть может, неограниченное множество, содержащее

состояние равновесия, соответствующее рабочему режиму системы [1], [3]. И тот факт, что траектория попала на заданное множество, не гарантирует того, что она в дальнейшем этого множества не покинет. Таким образом, возникает задача получения оценок областей притяжения, т.е. таких областей, которые траектория с течением времени не покинет. Для решения этой задачи может быть использован второй (прямой) метод Ляпунова или метод функций Ляпунова, когда изучается поведение некоторой вспомогательной функции вдоль траекторий системы [4 ], [5], [6].

Центральное место во втором (прямом) методе Ляпунова занимает проблема построения функции Ляпунова. При этом функцию Ляпунова нелинейной динамической системы часто ищут в классе положительно определенных квадратичных форм, исходя из того условия, что построенная квадратичная форма является функцией Ляпунова для соответствующей линеаризованной системы [5, С. 120-132]. Кроме того, может ставиться вопрос о построении квадратичных функций Ляпунова с некоторыми заданными свойствами, которые определяются особенностями задачи. В частности, в динамических задачах, когда интерес представляют собой ухе не только качественные, но и количественные характеристики системы, возникает необходимость ограничений на первую производную функции Ляпунова вдоль траекторий системы [7,8].

В настоящей работе предлагается способ получения аналитических оценок окрестности асимптотически устойчивого состояния равновесия непрерывной динамической системы, которую траектория системы с течением времени не покинет, основанный на применении прямого метода Ляпунова и использующий методику работ [9,10].

Пусть нелинейная непрерывная динамическая система задана системой дифференциальных уравнений

xi = fi(x1,x2,...,xn), (i = 1,2,...,п) (1)

имеющей асимптотически устойчивое состояние равновесие Xy = х2 = ... = хп = 0. И пусть функции f■ имеют непрерывные частные производные до второго порядка включительно. Запишем (1) в виде [4, С. 74-75]

х, =

I

]=1

alJxJ + Д1(х1,х2,...,хп), (I = 1,2,...,п)

(2)

где

а] =

X1

■■ х-

дх,-V ] У

: х„

- постоянные, равные значениям частных производных функций f^ при - остаточные члены второго порядка в разложении по формуле Тейлора

1

п п

Щх1,х2.....хп) = тХХ

д 21,

2 ]=1к=1дх] дхк

хх . (I = 1,2,...,п)

(3)

В силу непрерывности вторых производных для остаточных членов (3) на любом ограниченном множестве, содержащем состояние равновесия, можно указать оценку

Д (х1,х2,...,хп)\ < С1 х2, С = Ьп/2, (I = 1,2,...,п)

]=1

где

д2/,

дх]дхк

< Ь (i,],k = 1,2,...,п).

Выберем функцию Ляпунова в классе положительно определенных квадратичных форм

п п

У(х1,х2,...,хп) =11 К^х] (К] = Кр,1,] = 1,2,...,п)

I=1]=1

так, чтобы она была функцией Ляпунова соответствующей (2) линеаризованной системы

(4)

(5)

х1 = I а]х] ]=1

(I = 1,2,...,п)

(6)

и удовлетворяла при этом условию [6,7] ШОх(У /V) = —§ (2тах{^еЯ,} < < 0), где

0

I = 1,п

]У(х1,х2,...,хп) = I I К](хгх] + х1х])

I=1]=1

(7)

первая производная (5) в силу линеаризованной системы, а Я1,Я2,..,Яп, Л- < 0 (I = 1,2,...,п) -собственные числа матрицы системы А = (О] ) . В [8] показано, что выбор коэффициентов (5), удовлетворяющей этому условию, можно осуществить с помощью простых явных соотношений.

п

а

п

Оценим первую производную такой функции Ляпунова для нелинейной системы (2) подобно [4]. Следует отметить, что, согласно [10], в точках поверхности уровня V(x^,^2,■■■,^n) = V имеют место неравенства

^ < с1л[¥0,C1 = л/Aii /detK . (I = 1,2,...,п), (8)

где Aii - алгебраическое дополнение элемента Kii в матрице К = (Ку ), а det K - определитель матрицы K. Но тогда, согласно (4),

п

\aI(Xl,x2,...,xn)\ <му0, м = с£ с; (г = 1,2,..,п), (9)

з=1

а первая производная (5) в силу нелинейной системы (2) будет иметь вид

V п п

V (Х1,Х2,...,Хп) = V (Х1, Х2,..., Хп ) + ££ К3(х1П] +Х]П1) ' (10)

i=13=1

где V (Х1, Х2 ,■■■, Хп ) есть первая производная (5) в силу уравнений линейной системы. Так как на сечении V(Х1 ,Х2 ,■■ .,Хп ) < —SVg, то, учитывая ограничения (8), (9), получаем, что

~ пп _

V <—Ж0 + ££ Ку (С1 + С.

у ^г ' ^3/ г о^г о 1=1 з=1

(11)

Значит, неравенство V < 0 будет обязательно выполняться, если

п п _

— ^ + £ £1К у \(С1 + Су )М^¥0 < 0, (12)

i=1 3=1

поскольку Vо > 0 . То есть при

п п

0 <s (£ £ |Ку \(С + с} )М)-1 (13)

i=1 3=1

первая производная в силу нелинейной системы будет отрицательной. Таким образом, доказана следующая теорема.

Теорема. Область

^ =кХ1>Х2>...'Хп)^(Х1>Х2>...>Хп) < Ку } (14)

где

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

пп

Vоy = S2(££ Ку(Сг + С,)М)-2, (15)

0у = о2( ££ Ку(с + сз"' -2

,=1 3=1

является оценкой снизу области притяжения асимптотически устойчивого состояния равновесия Х1 — Х2 — ■■■ — Хп

= 0 системы дифференциальных уравнений (2).

Замечание 1. Оценка (14) может быть использования при определении точного момента окончания длительности переходного процесса в динамической системе, математическая модель которой представляет собой нелинейную систему дифференциальных уравнений вида (2), при численном исследовании поведения ее траекторий. Пусть

критерием окончания переходного процесса является попадание траектории системы (2) на заданное множество , содержащее устойчивое состояние равновесия системы. В случае, если множество (14)-(15) $ С Оа, его можно

принять за окрестность состояния равновесия, которую траектория системы не покинет с течение времени. Если это не так, рассмотрим множество

$ ={(xl,x2,■■■,xn):V(xl,x2,■■■,xn) < Кв } (16)

для поверхности уровня V(= , вписанной в . Например, в случае, когда есть «7 -полоса»

п

\ £Р^к\<7 , (П)

к=1

значение

а2 detK

К0 в = —-. (18)

II А}ргр} 1=1 ]=1

Скажем, при п = 2

тг а2(К11К22 — Ю

V в = ~-2 " 22-^-2 .. (19)

К22р2 — 2к12р1р2 + К22р1

Так как в этом случае К0в < V'0у , то V < 0 на множестве (16), и его можно принять за искомую окрестность

состояния равновесия. То есть в качестве окрестности состояния равновесии (2), которую траектория системы не покинет с течением времени, можно взять множество

={(х1,х2,...,х„) ■^(х1,х2,...,хп) < К0 ^ (20)

где V = тт{К0вУ0у}, а момент окончания переходного процесса в системе определять следующим образом.

За начальную точку траектории выберем точку, координаты которой соответствуют координатам состояния равновесия при значении параметра, с которого происходит переключение, по отношению к координатам состояния

равновесия х^ = х2 = ... = хп = 0, соответствующего значению параметра, на которое переключение

происходит. Вычисления будем проводить до тех пор, пока траектория системы не попадет на множество

содержащее х^ = х2 = ... = хп = 0 . Запомнив момент попадания траектории на , продолжим вычисления до

тех пор, пока либо траектория не попадет в область (19), либо не покинет множество . В первом случае будем

считать моментом окончания переходного процесса запомненный ранее момент попадания в , во втором -

продолжим вычисления до следующего момента попадания траектории на , с последующим повторением всей процедуры.

Замечание 2. Доказанная теорема может быть использована для получения оценок области притяжения состояния равновесия нелинейной системы (2) [9]. Пусть мы ищем оценку области притяжения состояния равновесия

х1 = х2 = ... = хп = 0 системы (2), принадлежащую ограниченному множеству |х;-1 < Ь^ (I = 1,2,...,п). Поскольку поверхности уровня К(ху ,х2,...,хп ) = V , вписанной в |х; | < Ь1 (I = 1,2,...,п), соответствует

значение К0т = тш{Ь^ detK/ Ац } [9], то искомой аналитической оценкой области притяжения является 1

множество

Ят ={(х1,х2,...,хп):К(х1'х2'...'хп) < т1п{К0тУву}\. (21)

Замечание 3. 1. Пусть все корни характеристического уравнения, соответствующего состоянию равновесия х1 = х2 = ...хп = 0, действительны и различны, причем < К2 < ... < Кп < 0. В этом случае всегда существует линейное невырожденное преобразование координат, приводящее систему (10) к каноническому виду (столбцы матрицы преобразования координат являются собственными векторами матрицы А , соответствующими собственным значениям Х,Х2,...,К соответственно). Рассмотрим систему дифференциальных уравнений вида

х1 =Л1х1 + Д1(х1,х2,...,хп) , (I = 1,2,...,п) (22)

Следует отметить, что согласно [8], квадратичную функцию Ляпунова для системы (22) можно выбрать в виде

п—2

К(х,,х7,...,х ) =!К х2 + К , ,х2 ,+ 2К , х ,х + К х2 (23)

\ 1 1 2 ' ' п / II I п-1,п-1 п—1 п—1,п п—1 п пп п к '

1=1

где

К2—1,п = (1 — Щ8))Кп—щ—Кпп, ЩЗ) = (Хп—; — \)2 (К—1+ Х — ¿Г2 (24)

2

(2Л3п < З < 0). Здесь Кц > 0,1 = 1,2,...,п,Кп — 1,п — 1Кпп — Кп—1,п > 0 . Тогда, в силу доказанной теоремы, получим, что в этом случае

п

К0у = З2 (2(IIК |С +4(1 — Щ8))КЯ_,п—}Кпп (Сп + С^))Ы)-2. (25)

1=1

2. Пусть все корни характеристического уравнения, соответствующего состоянию равновесия Xj = Х2 = ...Xn = 0, таковы, что A1 < A2 < ... < Ann_2 < Ц < 0, т.е. действительны и различны. И пусть

^n-l п = М — . В этом случае всегда существует линейное невырожденное преобразование координат,

приводящее линеаризованную систему к каноническому виду [4]. Рассмотрим систему дифференциальных уравнений

Xi = Aixi + Qi(x1,x2,...,xn), (i = 1,2,...,n - 2)

Xn-1 = Mxn-1 -Vxn +^n-1(xi,x2,...,xn) , (26)

xn =™n-1 +Mxn +OJxi,x2,...,x„) ,.

Заметим [7], что в этом случае квадратичную функцию Ляпунова, удовлетворяющую условию

max(V /V) = -8 (2max{Re A^j <-8 < 0), можно выбрать в виде (23), где V =V0 i=1,n

(К„—1п—1 + Кпп)2 = С(8)(Кп_пКп — К2—1п), (27)

С(8) = (—8 — 2ц)2У 2 + 4. В частности, в случае Кп_1п_ 1 = Кпп,

ГГЦ + 8/ + 4у2)К11п = (м + 8)2Кп_1п—1Кпп. (28)

Тогда, в силу доказанной теоремы, в этом случае

V0y = 82 (2(K^ W(2М - 8)2 /((2ц - 8)2 + 4v2)Knn(Cn + Cn^1))M)-.

i=1

Список литературы / References

1. Горюнов В. И. Техническая полоса захвата одноконтурного синтезатора частоты / В. И. Горюнов, В. Н. Ерусланов, Н. И. Лобашов // Техника средств связи. Сер. Техника радиосвязи. - 1990. - Вып. 2. - С. 88-94.

2. Горюнов В. И. К анализу быстродействия синтезатора при переключениях по диапазону / В. И. Горюнов, В. Н. Ерусланов, М. Н. Зайцева // Техника средств связи. Сер. Техника радиосвязи. - 1986. - Вып. 3. - С. 44-49.

3. Антоновская О. Г. К анализу формы и длительности переходных процессов при переключениях синтезатора с делителем частоты и пропорционально-интегрирующим фильтром по диапазону / О. Г. Антоновская, В. И. Горюнов, Н. И. Лобашов // Динамика систем: Межвуз. сб. Горький: Изд-во ГГУ, 1989. - С.59-72.

4. Ляпунов А. М. Общая задача об устойчивости движения. / А. М. Ляпунов.. - М.-Л.: Изд-во техн.-теор. лит., 1950. - 472 с.

5. Барбашин Е. А. Функции Ляпунова. / Е. А. Барбашин. - М. Наука, 1970. - 240 с.

6. Четаев Н. Г. Устойчивость движения. / Н. Г. Четаев. - М.: Наука, 1965. - 207 с.

7. Антоновская О. Г. О построении квадратичной функции Ляпунова с заданными свойствами / О. Г. Антоновская // Дифференциальные уравнения. - 2013. - Т.49. - № 9. - С. 1220-1224.

8. Антоновская О. Г. О построении квадратичной функции Ляпунова с заданными свойствами / О. Г. Антоновская // Дифференциальные уравнения. - 2016. - Т.52. - № 3. - С. 276-281.

9. Антоновская О. Г. Об одном способе оценки области притяжения устойчивого решения системы дифференциальных уравнений / О. Г. Антоновская // Международный научный журнал «Инновационная наука». -2015. - №9. - С. 11-15.

10. Антоновская О. Г. Об одном способе оценки размеров области притяжения неподвижной точки нелинейного точечного отображения произвольной размерности / О. Г. Антоновская, В. И. Горюнов // Известия вузов. Математика. - 2016. - № 12. - С. 12-18.

Список литературы на английском языке / References in English

1. Goryunov V. I. Tehnicheskaya polosa zahvata odnokonturnogo sintezatora chastoty [Technical capture strip of one-contour frequency synthesizer] / V .I. Goryunov, V. N. Eruslanov, N. I. Lobashov // Tehnika sredstv svyasi. Ser. Tehnika radiosvyasi [Techniques of communication. Series Techniques of radio-communication]. - 1990. Vyp. 2. P. 88-94. [in Russian]

2. Goryunov V.I. K analizu bystrodeystviya sintezatora pri pereklucheniyah po diapazonu [On synthesizer speed analysis when switching in a range] / V .I. Goryunov, V. N. Eruslanov, M. N. Zaytseva // Tehnika sredstv svyasi. Ser. Tehnika radiosvyasi [Techniques of communication. Series Techniques of radio-communication]. - 1986. Vyp. 3. P. 44-49. [in Russian]

3. Antonovskaya O. G. K analizu formy I dlitelnosti perehodnyh processov pri pereklyucheniyah sintezatora s delitelem chastity I proportsionalno-integrtiruyushim filtrom po diapazonu [On the analysis of form and duration of transient processes when switching over a range in the synthesizer with counter and proportional-integrating filter] / O. G. Antonovskaya, V. I. Goryunov, N. I. Lobashov // Dinamika system: Mezhvuz. sbornik [System dynamics: Inter-university Collection]. - Gorky: Izd-vo GGU, 1989. P. 59-72. [in Russian]

4. Lyapunov A. M. Obshchaya zadacha ob ustoychivosti dvizheniya [General problem of the stability of motion]. / A. M. Lyapunov. - Moscow-Leningrad : Izdat. Nthn.-teor. Lit., 1950. 472 p.

5. Barbashin E.A. Funktsii Lyapunova [Lyapunov functions]. /E. A. Barbashin. - Moscow: Nauka, 1970. 240 p.

6. Chataev N. G. Ustoychivost' dvizheniya [Stability of motion]. / N. G. Chetaev. - Moscow: Nauka, 1966. 207 p.

7. Antonovskaya O. G. On the construction of quadratic Lyapunov function with given properties / O. G. Antonovskaya // Differential equations. - 2013. - V. 49. - № 9. - P. 1187-1191.

8. Antonovskaya O. G. Determination of coefficients of a quadratic Lyapunov function with given properties / O. G. Antonovskaya // Differential equations. - 2016. - V. 52. - № 3. - P. 275-261.

9. Antonovskaya O. G. Ob odnom sposobe otsenki oblasti prityazheniya ustoychivogo resheniya sistemy differentsialnyh uravneniy [On a method of evaluation of attraction domain for fixed point of the stable solution of differential system] / O. G. Antonovskaya // Mezhdunarodnyy zhurnal "Innovatsionnaya nauka" [International magazine "Innovation Science"].- 2015. -№ 9. - P. 11-15. [in Russian]

10. Antonovskaya O. G. On a method of evaluation of attraction domain for fixed point of nonlinear point mapping of arbitrary dimension / O. G. Antonovskaya, V. I. Goryunov // Russian Mathematics.- 2016. - V. 60. - № 12. - P. 10-14.

DOI: https://doi.org/10.23670/IRJ.2017.62.087 Долапчи С.М.1, Денисова О.А.2

1 ORCID: 0000-0001-9609-5937, Южно-Уральский институт управления и экономики

2ORCID: 0000-0001-6374-3109, доктор физико-математических наук, доцент, Уфимский государственный нефтяной технический университет ВЛИЯНИЕ ПРОЦЕССА ТЕРМООБРАБОТКИ НА КАЧЕСТВО ПОВЕРХНОСТИ ИЗДЕЛИЙ ИЗ

КВАРЦЕВОГО СТЕКЛА

Аннотация

Основу прочной поверхности в силикатных стёклах составляют кремнекислородные тетраэдры, являющиеся наноструктурами для макрообразований, которые называются глобулами. Чем меньше размер глобул, тем прочнее поверхность, так как между ними образуется большее количество силоксановых связей. Данный фактор влияет на качество поверхности изделия, что улучшает индуктивность контура резонатора, следовательно, повышает коэффициент добротности. Размер глобул обеспечивает технология отжига. Процесс отпуска происходит при градиенте температур 300оС в час. Дальнейшее изменение температуры охлаждения ведёт к образованию трещин. Опыты, проведённые с помощью метода Виккерса и метода лежачей капли на дистиллированной воде и глицерине показали, что наибольшая прочность поверхности достигается в резонаторах, отожжённых в диапазоне температур 950 - 1080оС, а коэффициент добротности принимает максимальное значение.

Ключевые слова: энергия Гиббса, объёмная энергия, поверхностная энергия, отжиг, силоксановая связь, кремнекислородные тетраэдры, натрий, прочность, упругость, глобула, энтальпия, энтропия.

Dolapchi S.M.1, Denisova O.A.2

1 ORCID: 0000-0001-9609-5937, South Ural Institute of Management and Economics, 2ORCID: 0000-0001-6374-3109, PhD in Physics and Mathematics, Associate Professor, Ufa State Petroleum Technological University INFLUENCE OF THERMAL PROCESSING ON QUALITY OF PRODUCTS SURFACE MADE OF QUARTZ

GLASS

Abstract

Silicon-oxygen tetrahedra is the basis of a strong surface in silicate glasses. They are nanostructures for macroformations, which are called globules. The smaller the size of globules is, the stronger is the surface, because more siloxane bonds form between them. This factor affects the quality of the product surface, which improves the inductance of the resonator circuit, and therefore, increases the Q factor. The size of globules is provided by means of annealing technology. The tempering process takes place at a temperature gradient of300 °C per hour. A further change of the cooling tempera ture leads to the formation of cracks. Experiments carried out with the help of the Vickers method and the method of a recumbent drop on distilled water and glycerin showed that the greatest strength of the surface is achieved in resonators annealed in the temperature range 950 - 1080 °C, and the Q factor takes the maximum value.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Keywords: Gibbs energy, volume energy, surface energy, annealing, siloxane bond, silicon-oxygen tetrahedra, sodium, strength, elasticity, globule, enthalpy, entropy.

Большую роль в упрочнении поверхности изделии из кварцевого стекла играет такой вид термообработки как отжиг [1], [2]. Подбирая различные температурные режимы, можно добиваться различного размера глобул от

10 до 300 мкм. Перед нами стояла задача получить наиболее прочную силоксановую поверхность механических резонаторов, используемых в авиационной и космической отрасли. Для этой цели была разработана индукционная печь, управляемая микроконтроллером, и обеспечивающая равномерный градиент температур по всей поверхности изделия. Стекло кварцевое ультрафиолетовое КУ-1 (содержание натрия менее 1%) получают из кварцевого видимого стекла (КВ). Стёкла КУ-1 пропускают ультрафиолетовое излучение, КВ - нет. Резонатор представлял собой полусферу с внешним диаметром 20 мм и толщиной 2 мм, изготовленную из стекла КУ-1. Восемь резонаторов, выполненных по нашей технологии, показали коэффициент добротности (10 - 22) миллионов единиц.

Образцы для исследования изготавливали из кварцевого стекла марки КУ-1 первой категории производства ОАО «ТехноКварц». При термической обработке одну из пластин не подвергали отжигу, остальные прокаливали при температуре 900, 950, 1000, 1025, 1050, 1080, 1100, 1125, 1200, 1300оС в течение 5 часов, 3 часов прогрева и 3 часов отпуска.

i Надоели баннеры? Вы всегда можете отключить рекламу.