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




Скачать 276.9 Kb.
НазваниеПовышение качества стабилизации выходного напряжения преобразователей частоты на основе инверторов тока за счет применения адаптивных регуляторов
страница1/3
Дата конвертации12.11.2012
Размер276.9 Kb.
ТипАвтореферат
  1   2   3


На правах рукописи


Миргородская Екатерина Евгеньевна




ПОВЫШЕНИЕ КАЧЕСТВА СТАБИЛИЗАЦИИ

ВЫХОДНОГО НАПРЯЖЕНИЯ ПРЕОБРАЗОВАТЕЛЕЙ ЧАСТОТЫ

НА ОСНОВЕ ИНВЕРТОРОВ ТОКА ЗА СЧЕТ ПРИМЕНЕНИЯ АДАПТИВНЫХ РЕГУЛЯТОРОВ


Специальность: 05.09.12 – Силовая электроника


Автореферат


диссертации на соискание ученой степени

кандидата технических наук


Саратов – 2010

Работа выполнена в ГОУ ВПО «Саратовский

государственный технический университет»


Научный руководитель: доктор технических наук, профессор

Митяшин Никита Петрович


Официальные оппоненты: доктор технических наук, профессор

Рогинская Любовь Эммануиловна


кандидат технических наук, доцент

Беспалов Николай Николаевич


Ведущая организация: Институт проблем точной механики и управления РАН, г. Саратов


Защита состоится «23» декабря 2010 г. в 1300 часов на заседании диссертационного совета Д 212.242.10 при ГОУ ВПО «Саратовский государственный технический университет» по адресу: 410054, Саратов, ул. Политехническая, 77, Саратовский государственный технический университет, корп. 1, ауд. 414.


С диссертацией можно ознакомиться в научно-технической библиотеке ГОУ ВПО «Саратовский государственный технический университет».


Автореферат разослан «19» ноября 2010 г.

Автореферат размещен на сайте www.sstu.ru «19» ноября 2010 г.


Ученый секретарь

диссертационного совета Томашевский Ю.Б.


Общая характеристика работы

Актуальность проблемы. Преобразование частоты требуется для потребителей, нуждающихся в нестандартной частоте. В частности, таким потребителем является групповая двигательная нагрузка (ГДН) цеха внутришлифовальных станков. Другой важной областью применения мощных преобразователей частоты (ПЧ) являются локальные системы электроснабжения. В указанных областях в качестве инверторов целесообразно применение автономных инверторов тока (АИТ) ввиду требований большой мощности при высоком качестве кривой питающего напряжения. Функциональное назначение и условия эксплуатации нагрузки предъявляют к источнику электропитания ряд требований, тесно связанных с его динамическими свойствами. Одним из основных является требование стабилизации напряжения на нагрузке с заданной точностью при изменениях нагрузки и напряжения питающей сети. Допустимые значения отклонения напряжения частоты от номинальных значений в цепях различных приемников электроэнергии устанавливаются требованиями ряда ГОСТов. Так, ГОСТ 13109-97 устанавливает допустимые отклонения на выводах приемников электрической энергии в пределах от –5 до +5 % от номинального напряжения.

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

Определяющими для всего ПЧ являются динамические свойства автономного инвертора. Связь динамических свойств АИТ с эффективностью системы стабилизации выходного напряжения очевидна: АИТ наряду с фильтром звена постоянного тока (ЗПТ) и компенсатором реактивной мощности являются основными инерционными звеньями этой системы. Это является основанием для глубокого исследования динамических свойств АИТ, работающего совместно с входным фильтром и компенсатором. Важность таких исследований определяется также тем, что динамика АИТ взаимосвязана с другими показателями ПЧ. В частности, влияние динамики на КПД преобразователя проявляется при резкопеременной нагрузке, поскольку при этом возникают дополнительные потери энергии тем большие, чем больше колебательность переходного процесса и его длительность. Повышение входного коэффициента мощности ПЧ часто достигается применением дополнительных фильтрокомпенсирующих устройств, что также усложняет динамику преобразователя.

Исследованиям динамики преобразователей частоты посвящены работы многих советских и российских ученых. Особенно большой вклад в решение проблем динамики преобразователей частоты внесли Г. А. Белов, А. А. Булгаков, А. С. Васильев, Г. С. Зиновьев, А. Д. Поздеев, Ю. К. Розанов, Ю. Г. Толстов, В. А. Чванов, В. П. Шипилло. Несмотря на это, ряд вопросов, связанных с исследованием влияния отдельных факторов на динамические характеристики ПЧ на основе АИТ не получил достаточного освещения в литературе. Сказанное, в частности, касается использования этой информации при синтезе адаптивных систем управления преобразователями для указанных областей применения. Таким образом, исследования по повышению качества стабилизации выходного напряжения ПЧ за счет адаптации к изменению величины и характера нагрузки являются актуальными.

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

Для достижения поставленной цели в диссертации решаются следующие задачи:

  1. Разработка методики построения дискретных динамических моделей инверторов тока.

  2. Разработка методики исследования динамических свойств инверторов тока по спектрам матриц их дискретных динамических моделей.

  3. Исследование динамических характеристик основных схем АИТ.

  4. Разработка и исследование адаптивных нечетких систем стабилизации ПЧ на основе АИТ, построенных по методам Такаги-Сугено и Мамдани.

Методы исследования. Математическое моделирование, аппарат функций от матриц, методы нахождения собственных значений и собственных векторов матриц.

Научная новизна

  1. Разработана оригинальная методика построения дискретных динамических моделей АИТ, основанная на учете структурной симметрии трехфазных мостовых схем в сочетании с использованием аппарата функций от матриц.

  2. Проведенные исследования спектров матрицы дискретных динамических моделей АИТ позволили установить зависимость параметров переходных процессов от величины и характера нагрузки и параметров схемы.

  3. На основании проведенных исследований дискретных динамических моделей АИТ доказано, что для повышения качества стабилизации выходного напряжения ПЧ, работающего на изменяющуюся в широких пределах нагрузку, необходимо применение адаптивных регуляторов.

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

Практическая ценность и реализация результатов работы

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

Основные результаты диссертационной работы внедрены в ЗАО «Промышленная электроника» (г. Саратов) и в научно-производственной фирме ЭлЭс при Саратовском государственном техническом университете при синтезе систем управления преобразователей частоты для технологических целей.

Апробация работы

Основные положения диссертационной работы докладывались и обсуждались на IV Международной конференции «Методы и средства управления технологическими процессами» (Саранск, 2007), XX Международной научной конференции «Математические методы в технике и технологиях» (Ростов-на-Дону, 2007), XXI Международной научной конференции «Математические методы в технике и технологиях» (Саратов, 2008), Международной научно-практической конференции «Интернет и инновации» (Саратов, 2008), XXII Международной научной конференции «Математические методы в технике и технологиях» (Псков, 2009), V Международной конференции «Методы и средства управления технологическими процессами» (Саранск, 2009), XXIII Международной научной конференции «Математические методы в технике и технологиях» (Саратов, 2010).

Положения, выносимые на защиту

  1. Учет структурной симметрии трехфазных мостовых схем АИТ в сочетании с аппаратом функции от матрицы позволяет увеличить частоту дискретизации их динамических моделей в шесть раз, что способствует повышению статической и динамической точности систем стабилизации выходного напряжения ПЧ.

  2. Анализ спектра матриц дискретных динамических моделей АИТ позволяет установить характер влияния параметров силовой схемы преобразователя и величины и характера нагрузки на колебательность и время переходных процессов в ПЧ.

  3. Параметры дискретных динамических моделей и субоптимальных регуляторов системы стабилизации выходного напряжения ПЧ существенно меняются при изменении величины и характера нагрузки. Так, величина субоптимального значения демпфирующего коэффициента усиления регулятора k1 при изменении параметров нагрузки в диапазонах, характерных для ГДН, изменяется в 1.8-2 раза, что требует применения адаптивных регуляторов.

  4. Адаптивные нечеткие регуляторы, построенные по методам Мамдани и Такаги-Сугено, позволяют достичь высокого качества стабилизации выходного напряжения ПЧ, характеризующегося временем затухания переходных процессов, не превосходящим 3-7 периодов выходной частоты при изменении нагрузки в диапазоне, характерном для ГДН.

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

Публикации. По результатам диссертационной работы опубликовано 16 печатных работ, из них 1 в издании, рекомендованном ВАК РФ.

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

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

К ПЧ для централизованного электропитания технологической нагрузки предъявляется ряд требований: стабилизация напряжения на общих шинах с заданной точностью при изменениях нагрузки и питающего напряжения, высокое качество кривой выходного напряжения, высокие КПД и коэффициент мощности во всем диапазоне изменения параметров нагрузки и параметров питающего напряжения и ряд других. В работе показано, что выполнение большинства из перечисленных требований тесно связано с динамическими свойствами основных функциональных модулей ПЧ: регулируемого звена постоянного тока, фильтров на стороне постоянного и переменного тока, автономного инвертора, источников и компенсаторов реактивной мощности. Однако определяющими для характеристики всего ПЧ являются свойства автономного инвертора. Это тем более верно, поскольку исследование таких элементов ПЧ, как фильтров и источников и компенсаторов реактивной мощности, целесообразно проводить совместно с автономным инвертором. Таким образом, при исследовании динамических свойств ПЧ основное внимание должно быть уделено исследованию автономного инвертора.

В работе рассматривается динамика АИТ. Помимо высокого качества выходной кривой в пользу применения АИТ в преобразователях большой мощности говорят следующие обстоятельства. Во-первых, они строятся на основе дешевых однооперационных тиристоров. Во-вторых, в АИТ легко достигается значительное превышение напряжения выходной сети над напряжением входной сети без применения трансформаторов, что невозможно в АИН. В-третьих, в схемах ПЧ на основе АИТ легко достигается агрегирование нескольких преобразователей для увеличения выходной мощности с одновременным улучшением качества кривой выходного напряжения также без применения трансформаторов. Эта задача для инверторов напряжения пока не имеет эффективных решений. Кроме того, из-за наличия разветвленной непрерывной части схемы – батарея коммутирующих конденсаторов (БКК) на стороне переменного тока и реактор на входе – динамика АИТ значительно сложнее динамики инверторов напряжения и, следовательно, представляет больший теоретический интерес.

В качестве базовой схемы в работе принята схема трехфазного мостового параллельного инвертора, как наиболее часто применяемая на практике. Кроме того, рассматриваются схема АИТ с «расщепленной» БКК и ее двухмостовой вариант, и схемы АИТ с вентильно-реакторными компенсаторами. Динамика всех указанных инверторов исследуется отдельно, а также совместно с включенными в их состав входными LC-фильтрами. Это определено тем, что фильтры существенно влияют на исследуемые динамические свойства преобразователя, а их уравнения естественным образом стыкуются с уравнениями динамики самого инвертора.

При выборе метода исследований и типа динамической модели анализируются существующие подходы к анализу динамики инверторов: метод переключающих функций, z-преобразований, метод функций от матриц. В работе используется последний метод в сочетании с преобразованием уравнений динамики инвертора к интервалу симметрии.

При выборе вида модели в сравнении с моделями в виде уравнения «вход-выход» предпочтение в работе дается модели преобразователей в пространстве состояний (ПС). Это обусловлено, в частности, тем, что все координаты ПС рассматриваемых преобразователей могут быть измерены, то есть они являются полностью наблюдаемыми объектами. Этот факт может быть использован при построении регуляторов систем управления преобразователями. Поскольку модели в ПС наиболее приспособлены для аналитического синтеза регуляторов, то целесообразно строить дискретные динамические модели в ПС.

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

Поскольку единственной переменной управления в большом числе важных случаев является величина напряжения на выходе звена постоянного тока, то в работе рассматривается скалярное управление u(i), представляющее среднее за интервал инвертирования напряжения на выходе звена постоянного тока преобразователя. В результате дискретная динамическая модель (ДДМ) в данной работе для всех инверторов строится в виде

. (1)

Здесь i – номер коммутации; – вектор дискретных компонент ПС инвертора; – вектор управления; D и – матрица размерности и вектор модели, определяемые, как и размерность модели n, схемой рассматриваемого преобразователя.

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

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

Система дифференциальных уравнений (ДУ) АИТ после преобразования вектора непрерывных компонент ПС инвертора к компонентам вектора первого интервала с помощью формулы



примет инвариантный вид

,

где u – величина напряжения звена постоянного тока; K – блочная матрица преобразования.

Применив теперь для интегрирования последней системы ДУ аппарат матричной экспоненты, получим искомую модель в форме (1), где

. (2)

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

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

Для этого вводится непрерывная вектор-функция y(t), связанная с дискретным вектором следующим образом

.

Систему ДУ для этого огибающего вектора ищем в виде

,

где матрица S находится с помощью матричного логарифма:

.

Связь между дискретной и непрерывной моделями АИТ используется для переноса методики оценки характера переходного процесса в непрерывной модели на соответствующую оценку дискретной модели.

Оценка динамических свойств АИТ в работе производится по спектру матрицы D дискретной модели. Для этого решается задача о собственных векторах и собственных значениях матрицы D при различных значениях параметров преобразователя, величины и характера нагрузки, то есть находятся решения уравнения

.

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

В качестве основных показателей характера переходного процесса в ПЧ в работе рассматриваются колебательность и число интервалов коммутаций , после которых все составляющие переходного процесса уменьшаются в раз. Эти величины оцениваются по собственным значениям матрицы D .

В третьей главе приводятся результаты исследований динамики инверторов по спектрам матрицы их дискретных моделей. Исследования динамики АИТ проводились по следующей схеме.

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

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

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

4. Для типичных значений B и строятся годографы СЗ матрицы в зависимости от параметров выделенных элементов схемы. Анализ этих годографов может служить основанием для рекомендаций по выбору этих параметров с точки зрения улучшения динамических свойств АИТ.

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

На рис. 1 приведены годографы собственных значений (СЗ) матрицы динамической модели базовой схемы в зависимости от коэффициента загрузки B от 0.1 (отмечено цифрой 1) до 1 (отмечено цифрой 10) для двух значений коэффициента мощности нагрузки =0.3 (слева) и =0.7 (справа).

Аналогичные годографы компенсированного инвертора приведены на рис. 2. Как видно из этих рисунков, годографы содержат пять ветвей для базового АИТ и шесть ветвей для компенсированного АИТ, что соответствует размерностям пространств состояний указанных инверторов. При учете в базовом инверторе влияния входного Г-образного LC-фильтра в годографе появляются две дополнительных ветви (рис. 3). Расчеты проводились для следующих базовых параметров = 500 мкГн, = 2 мФ, = 2 мГн.

Проведенные исследования характера влияния отдельных параметров инвертора и нагрузки на положение и движение СЗ на комплексной плоскости позволяют идентифицировать ветви годографов, связав каждую из них с определенной частью преобразователя. Такая привязка является приближенной, но отражающей преобладающие связи между обсуждаемыми математическими объектами и физическими процессами. На рис. 4 приведен вид годографа базового инвертора с входным фильтром и компенсатором на выходе. Две пары ветвей комплексных СЗ, выделенные в квадратах 1 и 1’ (в левой полуплоскости) и 2 и 2’ (в правой полуплоскости), наблюдаются на всех приведенных выше годографах. Установлено, что пара ветвей 1 – 1’ наиболее сильно определяется параметрами нагрузки и БКК, то есть B и . В значительно меньшей степени СЗ этих ветвей зависят от параметров входного реактора инвертора . Пара ветвей 2 – 2’ более равномерно связана со всеми частями преобразователя, что является препятствием для декомпозиции динамической модели инверторов тока на модели входной и выходной цепей, то есть цепей постоянного и переменного тока.

Наличие компенсатора расщепляет ветвь действительных СЗ (см. рис. 1) на две ветви комплексно сопряженных СЗ 3 – 3’. Это означает, что по всем координатам исчезают апериодические составляющие переходных процессов. Это связано с образованием сложной колебательной системы, содержащей БКК и реакторы и . Наконец, учет влияния входного фильтра связан с парой собственных значений 4 – 4’. Из приведенных результатов следует, что при наличии фильтра именно она определяет



Рис. 1. Годографы СЗ матрицы базовой схемы АИТ при изменении B от 0.1 (отмечено цифрой 1) до 1 (отмечено цифрой 10) для двух значений коэффициента мощности нагрузки

=0.3 (слева) и =0.7 (справа)



Рис. 2. Годографы СЗ матрицы базовой схемы АИТ с компенсатором выпрямительного типа

при изменении B от 0.1 (отмечено цифрой 1) до 1 (отмечено цифрой 10) для двух значений коэффициента мощности нагрузки =0.3 (слева) и =0.7 (справа)



Рис. 3. Годограф базовой схемы

с входным Г-образным LC-фильтром

при изменении B от 0.1 до 1 и=0.3

Рис. 4. К идентификации ветвей

корневых годографов АИТ
  1   2   3

Добавить в свой блог или на сайт

Похожие:

Повышение качества стабилизации выходного напряжения преобразователей частоты на основе инверторов тока за счет применения адаптивных регуляторов iconСинтез схем автономных инверторов напряжения с улучшенным гармоническим составом выходного напряжения на основе эволюционного моделирования
Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Саратовский...

Повышение качества стабилизации выходного напряжения преобразователей частоты на основе инверторов тока за счет применения адаптивных регуляторов iconИнверторов напряжения и активных выпрямителей учебно
Шим. В приложении приведена программа для расчета выходного фильтра автономного инвертора напряжения

Повышение качества стабилизации выходного напряжения преобразователей частоты на основе инверторов тока за счет применения адаптивных регуляторов iconИзучение измерительнх преобразователей тока и напряжения цель работы
Цель работы: Изучение устройства и принципа действия измерительных преобразователей тока и напряжения

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

Повышение качества стабилизации выходного напряжения преобразователей частоты на основе инверторов тока за счет применения адаптивных регуляторов icon6 22013, рф, Свердловская обл
В частности это может быть достигнуто за счет применения бесконтактных пускателей и параметрических преобразователей напряжения -устройств...

Повышение качества стабилизации выходного напряжения преобразователей частоты на основе инверторов тока за счет применения адаптивных регуляторов iconВходной фильтр инвертора с шим напряжения, формирующий синусоидальное напряжение на нагрузке
Дов переменного тока, инверторов тока и напряжения на igbt-приборах в силу ряда присущих им преимуществ перед другими видами преобразовательных...

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

Повышение качества стабилизации выходного напряжения преобразователей частоты на основе инверторов тока за счет применения адаптивных регуляторов iconВнедрение преобразователей частоты компании Schneider Electric – эксперимент
Отчет комиссии о проведении эксперимента по внедрению преобразователей частоты пр-ва Schneider Electric (Франция) на повысительных...

Повышение качества стабилизации выходного напряжения преобразователей частоты на основе инверторов тока за счет применения адаптивных регуляторов iconСистема оперативного постоянного тока
Фск еэс россии и состоит из аккумуляторных батарей, блоков с предохранителями, зарядно-подзарядных агрегатов, щитов постоянного тока,...

Повышение качества стабилизации выходного напряжения преобразователей частоты на основе инверторов тока за счет применения адаптивных регуляторов iconИсследование электроприводов, электротранспорта и станций быстрой зарядки в случае применения различных типов аккумуляторов, каналов передачи данных систем и преобразователей питания и инфраструктуры «умных электросетей»
Преобразователи постоянного тока с инвертором на базе транзисторов igbt и промежуточным звеном высокой частоты для питания бортовых...


Разместите кнопку на своём сайте:
lib.convdocs.org


База данных защищена авторским правом ©lib.convdocs.org 2012
обратиться к администрации
lib.convdocs.org
Главная страница