Задачи по теме: Расчёт параметров оптических кабелей
Стр 1 из 2Следующая ⇒ ОСНОВНЫЕ ВЕЛИЧИНЫ, ФОРМУЛЫ И ТАБЛИЧНЫЕ ДАННЫЕ Показатель преломления среды: , (1) где e – диэлектрическая проницаемость среды; m – магнитная проницаемость среды.
Относительная разность показателей преломления «сердцевина-оболочка» оптического волокна (ОВ): , (2) где nс – показатель преломления сердцевины ОВ; nо – показатель преломления оболочки оптического волокна, на практике D = 0,003 – 0,01. Показатель преломления сердцевины градиентного ОВ меняется от сердцевины к оболочке: , (3) где r – текущий радиус; nc – максимальное значение показателя преломления на оси ОВ; D - относительная разность показателей преломления «сердцевина-оболочка»; a – радиус сердцевины; u – показатель степени изменения профиля показателя преломления, равный 1 для треугольного профиля, 2 – для параболического, ¥ - для прямоугольного. Показатель преломления сердцевины градиентного оптического волокна (при параболическом распределении показателя преломления): , (3а) где nс – показатель преломления в центре сердцевины; r – текущий радиус; а – радиус сердцевины волокна. Для направляемых мод в ОВ выполняется условие полного внутреннего отражения: . (4) В кварцевых ОВ для систем связи значение jпред полного внутреннего отражения составляет на практике 85-86о. Числовая апертура. Апертура - это угол между оптической осью и одной из образующих светового конуса, попадающего в торец волоконного световода, при котором в световоде в ОВ выполняется условие полного внутреннего отражения. Этот телесный угол qА определяет числовую апертуру (5) где noc – показатель преломления окружающей среды у торца волокна. Числовая апертура для ступенчатых световодов (при noc = 1):
. (6) Для градиентных световодов локальная числовая апертура: . (6а) Нормированная (характеристическая) частота – обобщённый параметр ОВ: , (7) где d – диаметр сердцевины волокна; l – длина волны.
Таблица 1 – Соотношения нормированной частоты, длины волны и радиуса сердцевины для nc=1,49 и при различных значениях коэффициента преломления оболочки
Одномодовый режим реализуется, если нормированная частота V < 2,405. Чем меньше разность Dn = nс – nоб, тем при большем радиусе световода и меньшей длине волны обеспечивается одномодовый режим. Число распространяемых мод: (8)
где u –показатель степени изменения профиля показателя преломления (ППП), равный 1 для треугольного профиля, 2 – для параболического, ¥ - для прямоугольного. Для ступенчатого профиля показателя преломления u = ∞ и число распространяемых мод: . (9) Для градиентных световодов с параболическим профилем u = 2 и число распространяемых мод: . (9а) Затухание в оптических волокнах. Потери излучения в световодах измеряются отношением интенсивности выходного (прошедшего) и входного (падающего) световых потоков: , где Pвых – мощность прошедшего через ОВ светового потока; Pвх – мощность входящего (падающего) на ОВ светового потока. Экспериментально потери излучения в волоконных световодах принято выражать в децибелах на единицу длины волокна 1 км (линейные потери): (10) Затухание является важнейшим параметром оптических ка́белей. Затухание обусловлено собственными потерями в волоконном световоде aс и дополнительными потерями, так называемыми кабельными, aк, обусловленными скруткой, а также деформацией и изгибами оптических волокон при наложении покрытий и защитных оболочек в процессе изготовления ОК:
. (11) Собственные потери волоконного световода состоят из потерь поглощения молекулами кварца и потерь рассеяния вследствие изменения направления распространения лучей на нерегулярностях и их высвечивания в окружающее пространство: . (11а) Часть мощности поглощается также и посторонними примесями (aприм), выделяясь в виде джоулева тепла. Примесями могут являться ионы металлов (никель, железо, кобальт и др.) и гидроксильные группы (ОН), приводящие к появлению резонансных всплесков затухания. В результате суммарные потери в ОК определяются из выражения: (11б) Затухание за счёт поглощения связано с потерями на диэлектрическую поляризацию, линейно растет с частотой, существенно зависит от свойств материала оптического волокна (tgd) и рассчитывается по формуле: , дВ/км, (12) где nс – показатель преломления сердцевины ОВ; tgd – тангенс угла диэлектрических потерь материала сердцевины ОВ (для кварцевого материала tgd» 10-11); l – длина волны, мкм. Потери на поглощение растут линейно с уменьшением частоты. Потери релеевского рассеяния обусловлены неоднородностями материала волоконного световода, расстояния между которыми меньше длины волны, и тепловой флуктуацией показателя преломления. Потери на рассеяние, называемое рэлеевским, определяются: (13) или подробнее , (13а) где Кр – коэффициент рассеяния (для кварца равный 0,8 (мкм4×дБ)/км), l – длина волны, мкм, k=1,38 ×10-23 Дж/К – постоянная Больцмана, b=8,1×10-11 м2/Н – коэффициент сжимаемости для кварца, Т»1500 К – температура затвердевания стекла при вытяжке. Потери на рэлеевское рассеяние определяют нижний предел потерь, присущих волоконным световодам. Потери на рассеяние (в отличие от потерь на поглощение) растут с изменением частоты существенно быстрее - по закону F4. Потери на рэлеевское рассеяние сильно проявляются в коротковолновой области спектра. Поглощение в ИК области уменьшаются с ростом частоты и определяются показательной функцией: , (14) где С и k – постоянные коэффициенты, характерные для поглощающей среды. Для кварца С = 0,9; k = (0,7 … 0,9) ×10-6 м.
В кварцевом световоде существует три окна прозрачности: первое окно прозрачности на длине волны 0,85 мкм, второе окно прозрачности на длине волны 1,3 мкм и третье окно прозрачности на длине волны 1,55 мкм. Наименьшее затухание имеется в третьем окне прозрачности. Второе и третье окна прозрачности разделены между собой так называемым пиком водяного поглощения с максимумом на длине волны 1383 нм. Затухание в оптических кабелях. Дополнительное затухание, обусловленное кабельными потерями (aк), состоит из суммы, по крайней мере, семи видов парциальных коэффициентов затухания: (15) где a¢1 – возникает вследствие приложения к ОВ термомеханических воздействий в процессе изготовления кабеля; a¢2 – вследствие температурной зависимости коэффициента преломления материала ОВ; a¢3 – вызывается микроизгибами ОВ; a¢4 – возникает вследствие нарушения прямолинейности ОВ (скрутка); a¢5 – возникает вследствие кручения ОВ относительно его оси (осевые напряжения скручивания); a¢6 – возникает вследствие неравномерности покрытия ОВ; a¢7 – возникает вследствие потерь в защитной оболочке ОВ. Дополнительные потери определяются в основном процессами рассеяния энергии на неоднородностях, возникающих вследствие перечисленных влияний, и частично увеличением потерь на поглощение энергии. Причинами увеличения потерь на поглощение являются остаточные осевые и поперечные напряжения в ОВ, могущие возникнуть при изготовлении кабеля. В ряде случаев микроизгибы могут существенно влиять на прирост ak. Значение потерь на одном микроизгибе может изменяться в пределах 0,01 – 0,1 дБ. Приращение затухания от микроизгибов a¢3 зависит от мелких локальных нарушений прямолинейности ОВ, характеризуемых смещением оси ОВ в поперечных направлениях на участке микроизгиба. Основными причинами появления микроизгибов являются локальные неосесимметричные механические усилия различного происхождения, приложенные к очень малым участкам ОВ. К микроизгибам следует отнести такие поперечные деформации ОВ, для которых максимальное смещение оси ОВ соизмеримо с диаметром сердцевины волокна. Особенностями микро-изгибов является то, что они, как правило, многочисленны, а расстояние между соседними микроизгибами существенно больше их размера. Общий вклад потерь, создаваемых микроизгибами, может быть значителен. Вследствие микроизгиба происходит ограничение апертурного угла излучения, распространяющегося по ОВ, и часть энергии излучается из ОВ. Зависимость приращения затухания от микроизгибов a¢3 можно определить из выражения:
(16) где k3 = 0,9 – 1,0; Nн - число неоднородностей в виде выпуклостей со средней высотой yн на единицу длины; а - радиус сердцевины; b - диаметр оптической оболочки; D – относительное значение показателя преломления; n1 и n2 – показатели преломления сердцевины и оболочки; Ео и Ес - модули Юнга оболочки и сердцевины ОВ. Затухание в местах соединений оптических волокон. Внешние потери. В отличие от традиционных кабелей затухание в сростках ОВ может достигать больших величин, соизмеримых с километрическим затуханием. На затухание сростка большое значение оказывает поперечное смещение и смещение осей. На рис. 1 приведены возможные варианты дефектов сопряжения ОВ (радиальное, угловое и осевое смещение).
Рис.1 Радиальное d (а), угловое q (б) и осевое S (в) смещение ОВ
При радиальном смещении дополнительные потери определяются из формулы: (17) где d – диаметр сердцевины ОВ; d – радиальное смещение ОВ. При угловом смещении дополнительные потери определяются из формулы: (18) где qА - апертурный угол волокна; q – угловое смещение ОВ. При осевом смещении дополнительные потери определяются из формулы: (19) где qА - апертурный угол волокна; d - диаметр сердцевины ОВ; S - осевое смещение ОВ. Приведенные выражения достаточно хорошо согласуются с экспериментальными, если значения d, q и S достаточно малы (d / d £ 1). Наиболее жесткие требования предъявляются при радиальном и угловом смещении. Наличие зазора между торцами волокон существенно меньше влияет на потери. Потери, вызванные различием показателей преломления оптических сред (френелевские потери): (20) Внутренние потери. Производство оптических волокон оставляет некоторые допуски на воспроизводимость их параметров. Потери в ОВ обусловлены различием: числовых апертур, диаметров сердцевины, диаметров модового поля, некруглостью, неконцентричностью сердцевины/оболочки. Если передающее ОВ имеет бо́льшую числовую апертуру, чем приёмное, то возникают потери. Для ступенчатых ОВ: (21) Когда диаметр передающего ОВ больше, чем диаметр приёмного, то возникают потери:
(22) Допуски на числовую апертуру и диаметры различны для разных ОВ. Суммарные потери в соединении многомодовых ОВ определяются: (23) При соединении одномодовых ОВ необходимо учитывать различие диаметров модовых полей. При гауссовом распределении поля диаметр модового поля определяется соотношением: (24) где а – радиус сердцевины ОВ, V – нормированная частота. Дисперсия сигналов в ОВ
Уширение импульсов в оптическом волокне: , (25) где tвх – ширина импульса на входе; tвых – ширина импульса на выходе (значения берутся на уровне половины амплитуды импульсов). Рис. 2.
Существует целая группа причин возникновения дисперсии (рис.2). В многомодовых ОВ преобладает модовая дисперсия, которая зависит от типа профиля показателя преломления в используемом волокне. Модовая (межмодовая) дисперсия возникает в многомодовых волокнах и обусловлена различием времени распространения различных мод в световоде. Модовая дисперсия определяется в основном профилем показателя преломления. Модовая дисперсия зависит от длины ОВ. Длина связи мод - это длина оптического волокна, после прохождения которой, в результате взаимного преобразования мод на нерегулярностях (обмен энергии между модами и их высвечивание) соотношение между мощностями различных мод становится практически постоянным. Модовая дисперсия в этом случае возрастает уже не по линейному, а по корень-квадратичному закону. Рис. 3.
Модовая дисперсия в многомодовых волокнах со ступенчатым профилем показателя преломления: при L £ Lc, (26) при L ³ Lc, (26а) где D – относительная разность показателей преломления; nс – показатель преломления сердцевины; с – скорость света (с = 3×105 км/с); L - длина оптического волокна; Lc – длина связи мод (5–7 км для ступенчатого многомодового волокна). В градиентных волокнах различные лучи также проходит различный путь. Однако их скорость v=c/n различна. Околоосевые лучи распространяются по короткой траектории, однако в среде со сравнительно высоким n, т. е. с малой скоростью, а периферийные лучи - по длинной траектории, но в основном в среде с низким n, т. е. с большой скоростью. В результате модовая дисперсия в градиентных ОВ уменьшается. Модовая дисперсия в градиентных волокнах: при L £ Lc, (27) при L ³ Lc, (27а) где D – относительная разность показателей преломления; nс – максимальное значение показателя преломления сердцевины; с – скорость света (с = 3×105 км/с); L – длина оптического волокна; L c – длина связи мод (10 – 15 км для градиентных МОВ).
Таблица 2. Значения модовой дисперсии tмод в ступенчатых и градиентных волокнах для различных длин линии и при различных соотношениях показателей преломления сердцевины и оболочки.
Хроматическая (частотная) дисперсия возникает из-за некогерентности источников излучения и делится на материальную и волноводную (иногда выделяют ещё так называемую профильную дисперсию). Материальная дисперсия определяется частотной зависимостью показателя преломления материала ОВ n = f (l). Уширение импульсов из-за материальной дисперсии: , (28) где Dl – ширина спектральной линии источника излучения (равная 0,1 – 4 нм для ЛД и 15 – 80 нм для СИД); с – скорость света; L – длина линии. Для инженерных расчётов в первом приближении можно использовать упрощённую формулу, не учитывающую форму профиля показателя преломления (для идеального ступенчатого профиля показателя преломления): , (28а) где М(l) – удельная материальная дисперсия.
Таблица 3. Типичные значения удельной материальной дисперсии для стандартных кварцевых ОВ.
С увеличением длины волны значение tмат увеличивается, а затем проходит через нуль и приобретает положительное значение. Для расчёта хроматической дисперсии используется эмпирическая формула Селмейера (для ОВ с несмещённой дисперсией): (29) где So – наклон дисперсионной кривой, пс/(нм2×км), lо – длина волны нулевой дисперсии.
Волноводная (внутримодовая) дисперсия обусловлена процессами внутри моды. Она характеризуется зависимостью коэффициента распространения моды от длины волны. Являясь составной частью хроматической дисперсии (так же как и материальная дисперсия), волноводная дисперсия зависит от ширины передаваемого спектра частот.
Уширение импульсов из-за волноводной дисперсии tвв находится из выражения: (30) где Dl – ширина спектральной линии источника излучения; nc – показатель преломления сердцевины ОВ; с – скорость света в вакууме; D – относительная разность показателей преломления; L –длина линии. Для инженерных расчётов можно использовать упрощённую формулу: (30а) где В(l) - удельная волноводная дисперсия, значения которой приведены в таблице.
Таблица 4. Типичные значения волноводной дисперсии в кварцевых ОВ.
Профильная дисперсия проявляется в реальных оптических волокнах, которые могут быть регулярными (например, с регулярной, геликоидальной скруткой), нерегулярными (например, нерегулярное изменение границы раздела ППП), неоднородными (например, наличие инородных частиц). К причинам возникновения профильной дисперсии относятся поперечные и малые продольные отклонения (флуктуация) геометрических размеров и формы волокна, например: небольшой эллиптичности поперечного сечения волокна; изменения границы профиля показателя преломления (ППП); осевые и внеосевые провалы ППП, вызванные особенностями технологии изготовления ОВ. Продольные флуктуации могут возникать в процессе изготовления ОВ и ОК, строительства и эксплуатации ВОЛС. В ряду случаев профильная дисперсия может оказать существенное влияние на общую дисперсию. Профильная дисперсия может проявляться как в многомодовых, так и в одномодовых ОВ. Физически происходит перекачка энергии между направляемыми, оболочечными и излучаемыми волнами.
Уширение импульсов из-за профильной дисперсии: (31) где nэф – эффективный показатель преломления; b – нормированная постоянная распространения; nгр – групповой показатель преломления сердцевины; Г – коэффициент локализации по мощности; V – нормированная частота; с – скорость света; D – относительная разность показателей преломления; L – длина линии; l – длина передаваемой световой волны. Это выражение справедливо для одномодовых волокон при реальной флуктуации границы раздела ППП. При наличии других внешних влияющих факторов её значение может значительно увеличиваться. Для инженерных расчётов можно использовать упрощённую формулу: (31а) где П(l) – удельная профильная дисперсия, типичные значения которой представлены в табл.5; Dl – ширина спектра излучения источника; L – длина линии.
Таблица 5. Типичные значения профильной дисперсии в ОВ.
Результирующее уширение импульсов за счёт модовой, материальной, волноводной и профильной дисперсий: (32) Уширение импульсов возникает не только исключительно при прохождении сигнала через оптическое волокно, но и за счёт прохождения сигнала через соединители, модулирующие, демодулирующие и другие устройства. В случае гауссовой формы импульсов все эти приращения длительности сигнала суммируются по квадратичному закону: (33) где tвых – длительность импульса на выходе фотоприемника; tвх – длительность импульса на входе излучателя; ti – уширение импульса в i-м элементе тракта. Как правило, наибольшие искажения в сигнал вносит оптический кабель. Поэтому, согласно (27), зная длительность импульсов на входе и выходе тракта, можно рассчитать величину уширения импульсов в оптическом кабеле по формуле (20), причем значения tвых и tвх берутся на уровне половины амплитуды импульсов. Критическая частота и критическая длина волны. По оптическим волокнам возможна передача электромагнитных волн длиной меньше критической длины волны. Значения критической частоты волны и критической длины можно определить из выражений: , , (34) где Pnm – значения корней функции Бесселя для различных типов волн; с – скорость света; d – диаметр сердцевины оптического волокна; nc и noб – показатели преломления сердцевины и оболочки. В этих выражениях параметр Pnm характеризует тип волны (моду). Для критических значений нормированной частоты Vnm = Pnm. Для оценки Pnm следует воспользоваться таблицей 6. Таблица 6. Значения корней бесселевых функций в зависимости от порядка n и номера корня m.
Как видно из таблицы, только одна волна НЕ11 не имеет ограничения по критической частоте. Расчёт элементов конструкций оптических ка́белей. Из большого количества разнообразных типов конструкций ОК можно выделить четыре группы: ка́бели концентрической повивной скрутки, ка́бели с профилированным сердечником, плоские ка́бели ленточного типа и ка́бели пучковой скрутки. При механическом расчёте ОК элементы ка́беля рассматривают как систему независимо деформирующих цилиндров, пренебрегая их поперечным взаимодействием при растяжении вдоль оси кабеля. Допустимое растяжение конструкции ка́беля, состоящей из n элементов, можно оценить по формуле:
, (35)
где d – коэффициент допустимого продольного растяжения конструкции кабеля; ki – коэффициент, учитывающий расположение i-гo элемента относительно оси конструкции кабеля; Ei – модуль продольной упругости материала i-гo элемента конструкции кабеля, Па; Si – поперечное сечение i-го элемента конструкции кабеля, м2.
Коэффициент, учитывающий расположение элементов конструкций ка́беля, имеет следующие значения: • ki = 1 для элементов, параллельных оси кабеля (испытывающих механическое напряжение); • ki = 0 для элементов, не испытывающих механическое напряжение (обычно оптических волокон); • ki = cosyi для элементов, испытывающих напряжения и расположенных по спирали с углом подъема y1 (обычно для оптических модулей можно считать угол подъема 30 – 35°).
Таблица 7. Значения некоторых параметров элементов конструкций ОК.
Коэффициент допустимого продольного растяжения d для разных типов конструкций ОК меняется в пределах 0,01 – 0,025. В случае продольного армирования кабеля одной центральной жилой из высокопрочного материала, ее сечение может быть определено из выражения: (36) где d – коэффициент допустимого продольного растяжения конструкции кабеля; ki – коэффициент, учитывающий расположение i-го элемента относительно оси конструкции кабеля; Еi – модуль продольной упругости материала i-го элемента конструкции кабеля; Si – поперечное сечение i-гo элемента конструкции кабеля, м2. Длину регенерационного участка ограничивает один из двух факторов: затухание или дисперсия. При определении длины регенерационного участка необходимо на первом этапе найти максимально допустимое расстояние (ограниченное затуханием световодного тракта), на которое можно передать сигнал, а затем его восстановить. Вторым этапом определяют пропускную способность оптического кабеля и находят длину трассы, на которую ещё возможно передавать оптические сигналы с заданной скоростью. В многомодовых ОВ длина регенерационного участка обычно лимитируется дисперсией, а в одномодовых ОВ лимитируется затуханием. При определении длины регенерационного участка, лимитированного затуханием, следует пользоваться выражением: (37) где Э – энергетический потенциал системы передачи; С – энергетический запас системы; aпас – дополнительные потери в пассивных компонентах ВОЛС (на вводе-выводе), дБ; aкаб – коэффициент затухания оптического кабеля, дБ/км; aсоед – потери в неразъемном соединении, дБ; Lстр – строительная длина кабеля, км. Энергетический потенциал системы передачи (Э) определяет максимально допустимое затухание оптического сигнала в оптическом кабеле, в разъемных и неразъемных соединениях на участке регенерации, а также другие потери в узлах аппаратуры. Энергетический потенциал определяется как разность между уровнем мощности оптического сигнала, введённого в волокно, и уровнем мощности на входе приемного устройства, при котором коэффициент ошибок регенератора не превышает заданного значения, установленного для данной системы передачи. Энергетический потенциал зависит от скорости передачи, технического уровня элементов электрооптических и оптоэлектронных преобразователей, длины волны используемого источника излучения и других факторов и задаётся для каждого вида аппаратуры ВОСП. В таблице 8 приведены справочные данные энергетического потенциала различных отечественных систем.
Таблица 8. Справочные данные различных отечественных систем ВОЛС.
Энергетический запас системы (С) обычно составляет 3 – 6 дБ, он необходим для компенсации эффекта старения элементов аппаратуры и оптического ка́беля, компенсации дополнительных потерь при ремонте оптического кабеля (потери на стыках кабельных вставок) и других отклонений параметров участка в процессе эксплуатации. Дополнительные потери в пассивных компонентах ВОЛС (Аа) составляют 3–5 дБ и возникают за счет разъемных соединителей, устройств соединения линейного кабеля со станционным и т. д. Ограничение длины регенерационного участка дисперсией. Длина регенерационного участка ограничивается также пропускной способностью оптического кабеля. Пропускная способность DF является одним из основных параметров ВОЛС, так как она определяет полосу частот, пропускаемую оптическим волокном, и соответственно объем информации, который можно передать по оптическому кабелю на длину регенерационного участка. Пропускная способность оптического кабеля существенно зависит от используемых в них типов оптических волокон (одномодовые, многомодовые ступенчатые, градиентные), которые могут иметь различные дисперсионные параметры. Связь между величиной уширения импульсов и полосой частот выражается соотношением: ; (38) где k – коэффициент, учитывающий форму оптического импульса (от 0,44 при гауссовой форме импульса до 0,6 при прямоугольных импульсах); t – уширение импульса.
Дисперсионные искажения существенно зависят от длины оптического волокна, поэтому DF нормируется на один километр оптического кабеля. Так, если на километровой длине оптического волокна происходит уширение импульса на t =10 нc, то его пропускная способность DF ограничена 44 МГц (при гауссовой форме импульса). Для того, чтобы оценить способность какого-либо участка ВОЛС (Lx) передавать информацию с определенной шириной полосы частот, при известной нормированной полосе пропускания ОК на один километр (DF1), для коротких линий, меньших, чем длина установившегося режима (Lx <Lc), следует использовать выражение: (при Lx < Lc). (39) Для линий, больших, чем длина установившегося режима (Lx > Lc), следует использовать выражение: (при Lx > Lc). (40) Длина установившегося режима передачи для ступенчатого многомодового волокна составляет 5 – 7 км, для градиентного волокна 10 – 15 км. Для одномодовых волокон, в которых распространяется один тип волны, следует считать Lс = 25 – 30 км. Для одномодовых волокон Lc является длиной установившегося режима (либо длиной модовой связи) ортогонально-поляризованных двух мод (направляемая мода представляется двумя ортогонально-поляризованными модами одного типа). Если известно не нормированное значение полосы пропускания ОК на один километр, а известна пропускная способность DF оптического кабеля длиной L, то можно воспользоваться выражениями: для коротких линий, у которых Lx<Lc; (41) для длинных линий, у которых Lx > Lc. (42)
Как видно из выражений, с увеличением длины линии уменьшается пропускная способность световода DF. На рис. 4 показана зависимость дисперсии и пропускной способности от длины линии.
Рис. 4. Рис. 5.
Для определения длины регенерационного участка, исходя из формул (35) и (36), строят график изменения DFx от длины трассы ВОЛС Lx. На основании этого графика для требуемой системы передачи определяют длину регенерационного участка. Скорость передачи системы соответствует необходимой ширине полосы пропускания регенерационного участка. На рисунке 5 показан пример определения длины регенерационного участка третичной цифровой системы. После нахождения длины регенерационного участка, лимитированного дисперсией, проводят её сравнение с найденной ранее длиной регенерационного участка, лимитированного затуханием, и выбирают наименьшее значение.
ПРИМЕРЫ РЕШЕНИЯ ЗАДАЧ
Воспользуйтесь поиском по сайту: ©2015 - 2024 megalektsii.ru Все авторские права принадлежат авторам лекционных материалов. Обратная связь с нами...
|