способ управления процессом изготовления оптического волокна и способ изготовления оптического волокна

Формула изобретения

1. Способ управления процессом изготовления оптического волокна, включающий измерение диаметра волокна при его вытягивании и стабилизацию натяжения волокна, отличающийся тем, что в процессе вытягивания измеряют поперечные колебания волокна, анализируют множество составляющих частот поперечных колебаний и определяют основную частоту, для чего из множества составляющих частот выбирают одну частоту с максимальной амплитудой, подтверждают, что основная частота является частотой с максимальной амплитудой, рассчитывают величину текущего натяжения волокна по измеренному диаметру и основной частоте поперечных колебаний, при отклонении рассчитанного натяжения от заданного стабилизацию ведут в соответствии с этим отклонением.

2. Способ по п.1, отличающийся тем, что для подтверждения основной частоты удваивают частоту с максимальной амплитудой, сравнивают удвоенную частоту с частотой других составляющих и при наличии составляющей с частотой, близкой к удвоенной частоте, выбранную частоту с максимальной амплитудой принимают в качестве основной.

3. Способ по п.1, отличающийся тем, что для подтверждения основной частоты из множества составляющих частот дополнительно выбирают последующую, вторую частоту с максимальной амплитудой, удваивают частоты с максимальной амплитудой, каждую удвоенную частоту с максимальной амплитудой сравнивают с соответствующей близкой частотой, а за основную частоту принимают ту частоту с максимальной амплитудой, результат сравнения удвоенной частоты которой с соответствующей близкой частотой наименьший.

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

5. Способ по п.4, отличающийся тем, что для определения основной частоты поперечных колебаний из множества составляющих частот выбирают и удваивают частоту колебаний с максимальной амплитудой, сравнивают удвоенную частоту с частотами других составляющих и при наличии составляющей с частотой, близкой к удвоенной частоте, выбранную частоту с максимальной амплитудой принимают в качестве основной.

6. Способ по п.4, отличающийся тем, что для определения основной частоты поперечных колебаний из множества составляющих частот выбирают и удваивают частоты по крайней мере двух составляющих с максимальными амплитудами, каждую удвоенную частоту с максимальной амплитудой сравнивают с соответствующей близкой частотой, а за основную частоту принимают ту частоту с максимальной амплитудой, результат сравнения удвоенной частоты которой с соответствующей близкой частотой наименьший.

Описание изобретения к патенту

Изобретение относится к изготовлению оптических световодных волокон и к способам контроля и управления натяжением волокон при их вытягивании.

При изготовлении стеклянных оптических световодных волокон волокно вытягивается из конца стеклянной заготовки, которая была нагрета до достаточно высокой температуры. Одним из более важных параметров процесса при производстве оптических световодных волокон является натяжение в волокне при его вытягивании и в особенности натяжение в волокне в области между горячей зоной и первым устройством для нанесения покрытия. Величина этого натяжения влияет на конечные свойства волокна, включая его диаметр, прочность на разрыв и, за счет явления влияния напряжения на оптические свойства, также на оптические свойства волокна.

С точки зрения технологии производства, натяжение волокна также влияет на общую устойчивость и производительность процесса вытягивания. Чрезмерное натяжение ведет к быстрому образованию сужений и окончательному разрыву волокна в области горячей воды. Если не приводить тщательного контроля, увеличение температуры в горячей зоне для снижения натяжения в волокне может привести к резонансным явлениям при вытягивании и колебаниям корня, которые, в свою очередь, могут дать колебательные изменения диаметра волокна, которые трудно контролировать с помощью известных приборов для контроля диаметра волокна. Колебания волокна при вытягивании могут также неблагоприятно повлиять на процесс нанесения покрытия.

Натяжение волокна связано с вязкостью стекла в корневой части заготовки, из которой волокно вытягивают, и со скоростью вытягивания волокна. Так как вязкость стекла является функцией температуры, натяжением при вытягивании можно управлять, регулируя температуру печи.

Температура самой печи может измеряться с использованием известных методов и приборов, таких как пирометры или термопары. Однако вследствие термодинамики процесса вытягивания, это управление не обеспечивает управления температурой корня заготовки. Часть заготовки выше корня служит поглотителем тепла, который снижает температуру корня. Если печь поддерживается при постоянной температуре, корень становится горячее с уменьшением размера заготовки. Печь с постоянной температурой поэтому даст уменьшение натяжения при вытяжке волокна с уменьшением размера заготовки в процессе вытягивания из нее волокна.

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

Наиболее близким по технической сущности является способ управления процессом изготовления оптических волокон, в котором натяжение контролируют механически посредством устройства с тремя колесами, два из которых прилегают к волокну с одной стороны, а третье колесо прилегает к волокну с другой стороны. Положение третьего колеса относительно двух других используется как мера натяжения волокна. Сигнал измерителя натяжения используется для управления температурой печи, из которой вытягивается волокно, только в процессе начальной настройки, когда на волокно еще не наносится покрытие.

Трехколесный подход имеет многочисленные недостатки. Трудно точно согласовать устройство с волокном так, чтобы не изменялся первоначальный путь волокна. Контакт трехколесного устройства с волокном действует на петлю регулирования диаметра волокна так, что уменьшается скорость вытягивания волокна. Также подвижное волокно может рваться при контактировании с трехколесным устройством. Измеритель натяжения предпочтительно устанавливать прямо под печью при вытягивании волокна с покрытием. Разрыв волокна в этом месте снижает производительность, так как требуется проведение всех операций, как в начале процесса вытягивания волокна.

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

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

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

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

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

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

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

На фиг. 1 представлена блок-схема устройства изготовления волокна; на фиг. 2 схематичное изображение устройства контроля положения волокна; на фиг. 3 типичная запись зависимости положения от времени с выхода полосового фильтра по фиг. 1; на фиг.4 преобразование Фурье сигнала положения волокна фиг. 3; на фиг. 5 единичная импульсная функция, используемая для цифрового сглаживания формы волны по фиг.4; на фиг. 6 сглаженный спектр Фурье; на фиг. 7 модифицированная единичная ступенчатая функция, используемая для получения производной спектра фиг.6.

На фиг.1 представлена система изготовления оптических волокон, в которой стеклянная заготовка 10 расположена вертикально в муфеле 11 печи для вытягивания волокон. Нагревательный элемент 12 подает тепло к по меньшей мере нижней части заготовки 10. После известного процесса запуска волочения волокно 14 вытягивают из нижней части 13 заготовки 10 тянущими валками 20. После выхода из муфеля 11 волокно проходит контроль диаметра в устройстве 15, которое выдает сигнал, используемый в петле обратной связи, которая регулирует скорость тянущих валков 20 для поддержания постоянного диаметра волокна. Волокно 14 затем проходит устройство контроля положения 16, охладительную трубку 17 и устройство нанесения покрытия 18. Покрытое волокно может также проходить через устройство полимеризации покрытия и, если желательно, через добавочные устройства нанесения покрытий (не показаны).

В соответствии с данным изобретением выполняют следующие этапы для избежания выдачи ошибочных сигналов натяжения при определении натяжения движущегося волокна:

а) с помощью устройства контроля положения 16 воспринимают движение в поперечном к направлению вытягивания направлении волокна 14. Может использоваться любое известное устройство определения положения;

в) воспринятое перемещение анализируют для определения составляющих частот, а также амплитуды каждой составляющей;

c) выбирают составляющую частоту fm₁ c наибольшей амплитудой;

d) анализируют остальные составляющие для уточнения, имеется ли вторая гармоника частоты fm₁. Этот этап проводят для подтверждения того, что составляющая с максимальной амплитудой действительно является основным тоном гармонических колебаний волокна. Вторая гармоника определяется как составляющая в определенном интервале частот выше или ниже 2fm₁, т.е. удвоенной частоты основной составляющей 2/fm₁/;

е) затем определяют натяжение волокна по подтвержденной основной составляющей гармонических колебаний волокна.

Этапы от в) до е) могут выполняться отдельными процессорами, цифровым компьютером, многопроцессорным компьютером и т.п. Этап в), например, может осуществляться анализатором спектра.

Соответственно предпочтительному исполнению, вторая верхняя составляющая частоты fm₂ также рассматривается как возможный кандидат на основную частоту гармонического колебания волокна. Остальные составляющие анализируют для определения, имеется ли вторая гармоника вблизи удвоенной частоты 2/fm₂/, т. е. удвоенной частоты fm₂. Компоненту вблизи 2/fm₁/ называют первой второй гармоникой, а составляющую вблизи 2/fm₂/ называют второй второй гармоникой. Если вторая вторая гармоника ближе по частоте к 2/fm₂/, чем первая вторая гармоника к 2/fm₁/, то тогда основной частотой гармонического колебания принимают fm₂. Подобным же образом можно было определить, имеется ли частотная составляющая вблизи 2fm₃, удвоенной частоты третьей более высокой частотной составляющей.

Компьютер управления вытяжкой 47, который может быть типа Digital 11/73, содержит алгоритмы управления скоростью вытяжных валков, скоростью подачи вниз заготовки и температурой печи. Диаметр волокна 14 измеряется контрольным устройством 15, выдающим на управляющий вытяжкой компьютер сигнал, который позволяет динамически рассчитать линейную плотность волокна. Характеризующая длина вибрирующего волокна может быть измерена по фактическим размерам оборудования, но предпочтительно определяется экспериментально. Характеризующая длина может быть определена вместе с калибровкой системы посредством измерения основной гармонической составляющей и использованием трехколесного контактного измерителя натяжения для периодического измерения истинного натяжения. Этот калибровочный этап требуется лишь при начальном вводе устройства в эксплуатацию, если не меняется конфигурация оборудования в эксплуатации. Для верности калибровка может проводиться периодически, например ежегодно.

Сигнал, пропорциональный основной составляющей частоты на компьютер управления вытяжкой, который затем рассчитывает натяжение волокна Т посредством комбинации этой частоты с текущим диаметром волокна, рассчитывая по формуле

T = +d²² (2),

d диаметр волокна; - основная гармоническая частота, а коэффициенты и - экспериментально определяются и являются объединенными параметрами из уравнения (1).

Как показано на фиг.1, сигнал с алгоритма управления 48 подается на схему управления температурой 49 (которая может быть частью компьютера вытяжки 47), которая, в свою очередь, управляет температурой нагревательного элемента 12. Температурой можно, например, управлять так, чтобы поддерживать натяжение волокна постоянным.

Для реализации данного изобретения использовалась следующая система.

Измерение вибрации волокна осуществлялось оптическим устройством, показанным на фиг. 2. Свет лазера 24 рассеивался по горизонтальной оси с помощью цилиндрической линзы 26, для создания сравнительно широкой области, где осуществляется подсветка вибрирующего волокна. Пунктирные линии 27 показывают пределы расширения света. Свет попадает на волокно 14, отражается и преломляется во всех направлениях; преломляется в основном в передних углах. Так как волокно относительно мало по сравнению с размерами остальной части системы, волокно практически является точечным источником рассеянного света. Линзы 30 и 31 расположены примерно на расстоянии 1,2 дюйма (30 мм) от волокна для рефокусировки подсвеченного волокна в точки на поверхности детекторов 32 и 33 соответственно. Расстояние волокна от линз, расстояние детекторов от линз и фокусное расстояние линз определяют увеличение системы. Увеличение это отношение движения волокна перпендикулярно центральной оси линзы к движению рефокусированного пятна по соответствующему детектору. Каждый из этих детекторов является одноэлементным кремниевым детектором, имеющим два анода и один катод. Относительная разность между потенциалами на двух анодах и общем катоде является мерой положения света, сфокусированного на детекторе. С тем, как световое пятно движется вдоль детектора, разность потенциалов между анодами и общим катодом также меняется. Если пятно находится в центре детектора, потенциалы будут равны. Так как линза рефокусирует на детекторе изображения подсвеченного волокна, разность потенциалов может быть связана с положением волокна. Используются два детектора для одновременного измерения положения волокна в двух перпендикулярных плоскостях. Сигналы с обоих детекторов используются в петле обратной связи для позиционирования заготовки в горизонтальной плоскости (по осям x и y).

Так как для получения вибрационной информации достаточно сигнала от одного детектора, показана выходная схема только одного детектора 33. В конкретном исполнении вибрации измерялись в плоскости, перпендикулярной к поверхности волочильной башни. Два выходных сигнала с детектора 33 усиливаются усилителями 34 и 35, и выходы усилителей вычитаются друг из друга в схеме 36 для получения разностного сигнала на выходе 37. Для осуществления усиления и вычитания использовался усилитель, имеющийся в продаже. Результирующий сигнал после усиления и вычитания на выходе 37 представляет мгновенное положение волокна.

Подобный усилительно-вычитающий блок использовался с детектором 32 для получения добавочного сигнала, требуемого для позиционирования заготовки 10 в плоскости x-y.

Сигнал положения волокна на выходе 37 подается на преобразователь аналог-цифра 41 через усилитель с полосовым фильтром 39. Дополнительно к усилению в 50 раз, схема 39 пропускает только частоты в диапазоне от 1 до 40 Гц на преобразователь аналог-цифра (АЦП) 41. Фильтрация необходима для удаления сигнала номинального положения волокна, чтобы оставить только вибрационную информацию, а также чтобы удалить высокочастотный шум, который может вызвать ошибки в последующих расчетах. Типичная форма сигнала положения в зависимости от времени со схемы 39 показана на фиг.3.

Усиленный сигнал положения волокна преобразуется в цифровой сигнал преобразователем аналог-цифра 41, который берет пробы сигнала положения 100 раз в секунду с разрешением 12 битов (разрядов двоичного числа). Период отбора проб точно выдерживается точными электронными часами, входящими в состав преобразователя 41, которые используются в качестве опорного генератора. За период времени 10, 24 с получают 1024 точек положения волокна в виде цифр. Время отбора проб и количество точек положения определяют разрешение и рабочий диапазон устройства.

Фактическое натяжение волокна непрерывно меняется на несколько граммов. Основной источник этих изменений схема управления диаметром волокна. Скорость вытягивающих валков периодически корректируется, например раз в секунду. С увеличением диаметра, или уменьшением его относительно заданной величины, валки протягивают волокно быстрее или медленнее среднего значения скорости. Натяжение волокна изменяется, так как оно является функцией скорости вытягивания. По этой причине спектры амплитуды устремляются, как показано ниже, в этапе 3 расчета. Движение волокна также содержит многие компоненты, которые не являются элементами системы стоячей волны. Были разработаны алгоритмы с целью получения точного, надежного и устойчивого показания натяжения в присутствии всех этих помех шумового характера. В девяти этапах расчета, обсужденных ниже подробно, этапы 1-5 выполняются многопроцессорным компьютером 45. Data Translator модель 7020, a этапы 6-9 были потом выполнены обзорным компьютером 46, основанным на микропроцессоре 80386. Чтобы выполнить этапы расчета 1-9, эти имеющиеся в продаже компьютеры были запрограммированы по обычной технологии.

Расчеты осуществляются следующим образом.

1. Быстрое преобразование Фурье (БПФ) является первым из ряда расчетов, производимых для определения основной составляющей частоты вибрации волокна, БПФ рассчитывается на последнем наборе оцифренных данных, содержащем 1024 точек данных. Расчет дает 512 однозначных комплексных коэффициентов ряда Фурье.

2. Величина комплексных коэффициентов ряда Фурье определяется извлечением квадратного корня из суммы квадратов действительной и мнимой частей комплексных коэффициентов Фурье. Это дает 512 действительных чисел (половину числа точек данных) и называется амплитудным спектром. Типичный результат этого расчета можно видеть на фиг.4.

3. Усредняют текущий амплитудный спектр вместе с последними восемью ранее рассчитанными спектрами (можно усреднять не более или не менее чем 8 спектров). Результат здесь по-прежнему 512 действительных чисел.

4. Результат операции усреднения сглаживается цифровым способом. Это осуществляется сверткой результата этапа 3 с единичной импульсной функцией. Операция подобна "усреднению проб с запоминанием". Единичная импульсная функция показана на фиг.5. Многопроцессорный компьютер, используемый для этого расчета, не в состоянии проводить операцию свертки, так что на практике БПФ рассчитывается как для единичного импульса, так и для среднего спектра, и два комплексных набора множатся друг на друга. Затем рассчитывается обратная функция быстрого преобразования Фурье результата. Это свойство сверток и быстрых преобразований Фурье очень распространено. Результатом этого этапа опять является набор из 512 действительных точек данных. Для увеличения сглаживания единичная импульсная функция делается более широкой. Если ширина единичной импульсной функции равна только единице, то сглаживания не будет. Типичный результат этой операции можно видеть на фиг.6, который получен цифровой операции сглаживания, где ширина единичной ступенчатой функции была равна 6.

5. Результат этапа 4 основа подвергается свертке, но на этот раз с измененной ступенчатой функцией, показанной на фиг.7, где ширина равна 6. Фактическая операция осуществляется, как по этапу 4 (перемножением БПФ). То есть, многопроцессорный компьютер рассчитывает БПФ как для модифицированной ступенчатой функции фиг. 7, так и для сглаженного сигнала, показанного на фиг. 6, и два набора комплексных чисел перемножаются. Затем рассчитывают обратной БПФ результата. В результате этого этапа получают набор из 512 действительных чисел, представляющих производную результата этапа 4. Производная является практически производной по ширине модифицированной единичной функции, использованной при свертке. Если ширину модифицированной ступенчатой функции увеличить, результат расчета будет более гладким, но теряется разрешение.

6. Результат этапа 5 теперь просматривают для нахождения всех частотных "пиков", т. е. точек, где наклон сглаженного частотного спектра равен нулю. Это производится нахождением всех точек функции из этапа 6, нулевого значения. В большинстве случаев, вследствие конечного количества точек, точное прохождение производной через нуль находится между точками функции. Поэтому пик также определяется тем, что он находится между двумя точками с положительным наклоном ниже по частоте (левее), и отрицательным наклоном выше по частоте (правее). Добавочно, чтобы предотвратить попадание шумовых точек спектра, также требуется, чтобы величина сглаженного спектра в точке предполагаемого пика была бы выше порога. В результате этого этапа получают список всех точек функции этапа 5, где величина равна нулю, или соседние величины имеют противоположные знаки, а величина спектра по этапу 4 выше порога шумов.

7. Пикам частотного спектра из этапа 6 придается ранг по величине пика. Пику с максимальной амплитудой дается высший ранг и т.д.

8. Три наибольших пика анализируют на наличие вторых гармоник. Пик, который имеет составляющую, самую близкую к положению идеальной второй гармоники, признается как основная частота вибрации. Если ни один из максимальных пиков не имеет другого пика в пределах 1,5 Гц от положения идеальной второй гармоники, тогда для безопасности считается, что фундаментальной частоты нет.

9. Этапы 1-8 должны периодически повторяться. Однако после определения основной частоты вибрации предпочитают "прицепиться" к этой частоте в дальнейших расчетах. Для последующих расчетов по алгоритму, алгоритм рассматривает лишь частоты пиков в пределах "диапазона захвата" предыдущих расчетов для этапов 7 и 8. Если в этом диапазоне не оказывается пиков частоты, тогда работают по полному алгоритму. Это используется для предотвращения рассмотрения в расчетах неправильных сигналов. На практике основная частота не может очень быстро меняться, так что диапазон частот, который должен рассматриваться алгоритмом, сужается, что соответствует физике процесса. Полезным диапазоном частот, в котором следует искать основную частоту вибрации, является 2 Гц. Выполнение требования по наличию второй гармоники должно проверяться, если пик частоты находится в диапазоне захвата.

Следующий гипотетический пример типичен для работы по способу изобретения. Выполнение этапа 6 в компьютере 46 дает в результате следующие показатели наличия пиков частот:

Частота, Гц Амплитуда

8,1 75

11,0 90

16,5 35

23,1 20

30,0 17

Выполнение этапа 7 дает следующее расположение по рангам амплитуды:

Частота, Гц Амплитуда

11,0 90

8,1 75

16,5 35

23,1 20

30,0 17

В этапе 8 три верхних пика анализируются на присутствие второй гармоники (см. таблицу).

Решают, что 8,1 Гц является основной гармоникой, так как пик частоты удален лишь на 0,3 Гц от идеальной второй гармоники. Частотная составляющая, имеющая наибольшую амплитуду, исключается, так как разность между 22 Гц (удвоенная 11,0 Гц) и ближайший пик частоты составляет 1,1 Гц, что более 0,3 Гц, и поэтому исключается из рассмотрения. Третий по величине пик частоты 16,5 Гц исключается из рассмотрения, так как ближайший ко второй гармонике частотный пик отличается от нее на 3 Гц, что за пределами допустимой разности 1,5 Гц. Таким образом, пики на 10,5 Гц и 16,5 Гц являются шумовыми.

Надписи на чертежах:

фиг. 1: 15 устройство контроля диаметра; 16 устройство контроля положения; 18 устройство нанесения покрытия; 39 усилитель с полосовым фильтром; 41 преобразователь аналог-цифра; 45 многопроцессорный компьютер; 46 наблюдательный компьютер процессора натяжения; 47 компьютер управления вытяжкой; 28 алгоритм управления; 49 управление температурой;

Фиг.2: 24 лазер; 36 устройство вычитания;

Фиг.3: по оси абсцисс секунды;

Фиг.4,6: по оси абсцисс герцы.