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


Скачать 1.34 Mb.
Название Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет»
страница 5/10
Тип Реферат
rykovodstvo.ru > Руководство эксплуатация > Реферат
1   2   3   4   5   6   7   8   9   10
ГЛАВА 4 РАЗРАБОТКА АНАЛОГО-ЦИФРОВОГО КОМПЛЕКСА ДЛЯ ПАРАЛЛЕЛЬНОГО МОДЕЛИРОВАНИЯ ФИЗИЧЕСКИХ И МАТЕМАТИЧЕСКИХ ПРОЦЕССОВ

4.1 Назначение аналого-цифрового комплекса

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

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

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

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

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

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

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

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

Основными преимуществами математического моделирования перед физическим являются следующие:

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

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

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

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

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

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

Для реализации приведённых рекомендаций предлагается использование аналого-цифрового комплекса [34]

Задачей, на решение которой направлен АЦК, является создании эффективного способа моделирования процесса газификации остатков жидкого ракетного топлива в баках отделяющейся части ступени РКН и устройства на основе физической и математической моделей, путем повышения точности, достоверности, верификации и валидации получаемых результатов [25,35]

Основной проблемой при разработке АЦК являются:

- определение структуры АЦК;

- выбор сложности математической и физической моделей;

- организация функционирования (внутренних связей) внутри АЦК, в том числе и синхронизация процессов математического и физического моделирования;

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

Возможность использования моделей различной сложности:

- различные физические модели, например, элементы емкости бака (ЭММУ) и топливный отсек реального бака, модельный компонент (вода) и реальный компонент (керосин);

- различные математические модели, например, уравнение Навье - Стокса и метод теплового баланса;

- оценка на ранних этапах исследования возможных направлений конструирования бортовых систем, например, системы газификации остатков топлива в баках [36]

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

Конкретные виды физической и математической моделей определяют из условия совместимости по исследуемым параметрам, скоростям протекания процессов, например, использование уравнения Навье - Стокса, решаемого в пакете ANSYS требует на 1 секунду реального процесса в среднем 300-600 секунд вычислительного времени. Таким образом, использование таких сложных математических моделей в составе АЦК нереально.




Рисунок 4.1- принципиальная схема АЦК: 1- компрессор, 2- теплоэлектрического нагревателя, 3- дроссель, 4,5,6,8,9,10 – датчики, 7- ЭММУ, 11- контроллер, 12,13 – АЦП с интегрированной системой верификации данных, 14 – управляющий комплекс, 15 – исходные данные, 16 – математическая модель, 17 – блок обработки данных, 18 сервер совместного расчета параметров и критериев.
Последовательность работы АЦК (рисунок 4.1):

-запуск стенда (запуск компрессора 1 или забор теплоносителя из ресивера, теплоэлектрического нагревателя 2, открытие дросселя 3 максимальное, для наиболее быстрого прогрева системы);

- до выхода на рабочий режим (достижения заданной температуры на выходе из теплоэлектрического нагревателя 2 и выходе из экспериментальной емкости) осуществляется измерение только температуры внутри экспериментальной емкости 7;

- после выхода стенда на рабочий режим,датчики давления 4,10, температуры 5, 8 (вход и выход экспериментальной емкости 7), расхода 6 на входе в экспериментальную емкость 7 и процентного состава газа 9, подают данные на контроллер11.

Контроллер 11 обрабатывает данные и посылает их на аналого-цифровой преобразователь с интегрированной системой верификации данных 12 для перевода показаний датчиков 4-6, 8-10в цифровой.

Управляющий комплекс13, совмещенный с хранилищем базы данных (показание всех датчиков), накапливает получаемые данные и задает различные варианты рабочего давления подачи теплоносителя через дроссель3 в магистрали давления, температуру теплоэлектрического нагревателя 2, согласно заранее заданной программе проведения эксперимента 14 (циклограммы секундного расхода теплоносителя и его температура).

Управляющий комплекс13 контролирует поступающие измерения с датчиков 4-6, 8-10.

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

Таблицы полученных параметров поступают в блок совместного расчета параметров и критериев для обработки сервером совместного расчета параметров и критериев17 [37]

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

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

Данные вносятся в базу данных для сравнительного анализа и моделирование изменений текущих условий протекания физического процесса (изменение параметров во время эксперимента, например температура теплоносителя и его массовый расход) [38]
4.2 Математические модели АЦК

Математическая модель процесса испарения в общем виде одинакова как для ЭМУ №1, ЭМУ №2, так и для топливного бака, однако наличие различных устройств, их количества приводит к необходимости конкретного описания каждого объекта.

На рисунке 4.2 приведена общая схема взаимодействия тепловых потоков физических объектов (ТН, стенки бака, газ наддува, испаряемая жидкость), участвующих в процессе теплообмена внутри ёмкости (бака), которая иллюстрирует достаточно сложный механизм теплообмена.

Рисунок 4.2 – Тепловые потоки и их составляющие (конвективные, кондуктивные, излучение) между участниками взаимодействия

На рисунке 4.3 приведено фото ЭМУ№1, а на рисунке 4.4 приведено фото ЭМУ №2

c:\users\валерий\yandexdisk\рффи_2012\отчёты _10-08-00064\1.jpg

Рисунок 4.3 - Экспериментальная модельной ёмкости №1 для исследования процессов испарения модельных жидкостей (вода, спиртовая смесь) и керосина при параметрах ТН: избыточное давление в ёмкости 0 – 1,5 атм., массовый секундный расход 100 – 500 л/мин.; углы ввода теплоносителя 00 , 300 , 450; температура теплоносителя 1000 С – 1500 С; скорость входа ТН ….м/с; состав теплоносителя – воздух; состав газовой подушки – воздух.
c:\users\валерий\yandexdisk\рффи_2012\отчёты _10-08-00064\2.jpg
Рисунок 4.4 - Экспериментальная модельной ёмкости №2 в вакуумной камере для исследования процессов испарения модельных жидкостей (вода, спиртовая смесь) и керосина в вакуумную камеру. Диапазоны давления в вакуумной камере 0,1 – 1,0 атм. Производительность форвакуумного насоса для откачки газифицированной смеси до 500 л/мин. Абсолютное давление в ёмкости 0,1 – 3,5 атм., массовый секундный расход 100 – 500 л/мин.; углы ввода теплоносителя 00 , 300 ; температура теплоносителя 1000 С – 1500 С; скорость входа ТН ….м/с. Состав теплоносителя – воздух, азот; состав газовой подушки – воздух, гелий.

Исходная система уравнений для исследования процесса испарения жидкости взята в соответствии с работой [14], при этом она преобразована с учётом использования не химического механизма получения тепла, а за счёт подачи горячих газов, т.е. термодинамического источника тепла.

Изменение массы продуктов испарения в объеме бака можно записать в следующем виде:

, (4.1)

где: - массовый расход ТН;

- массовая скорость испарения жидкости;

где - коэффициент массоотдачи определяется экспериментально [24];

- парциальное давление испарившейся жидкости у поверхности жидкости, принимаемое равным давлению насыщенных паров при температуре поверхности (имеются табличные значения[24]);

- текущее парциальное давление испарившейся жидкости в свободном объёме бака;

- площадь испарения жидкости, является граничным условием, например, «зеркало»;

- массовый расход жидкости уходящей через дренажное отверстие,

где - коэффициент расхода;

- площадь выходного сечения дренажного отверстия;

- давление в смеси газов в баке (газ наддува испарившаяся жидкость и ТН);

- температура газовой смеси в баке;

- газодинамическая функция.

Начальные условия в τ=0 , где - масса воздуха в баке.

Уравнение изменения плотности газа в емкости имеет вид:

, (4.2)

где - объем бака.

Начальное условие - плотность воздуха.

Уравнение изменения давления в газовой подушке в соответствии с [14]:

, (4.3)

где - суммарное поступление тепловой энергии в емкость;

- энтальпия ТН;

- энтальпия смеси газов (ТН, воздух, испарившаяся жидкость) на выходе из бака;

- адиабатическая постоянная.

Температура газовой подушки в емкости определяется по формуле Менделеева –Клапейрона [39]:

(4.4)

где R – газовая постоянная смеси газов.

Распределении тепловой энергии (уравнение теплового баланса) имеет вид:

, (4.5)

где

- энергия расходующаяся на нагрев газа,

- конвективная составляющая энергии,

где - коэффициент теплоотдачи ТН, является функцией от температуры и давления

();

- площадь стенок емкости, в которую заключен газ;

- температура ТН

- температура газа;

- лучистая составляющая энергии,

где - постоянная Стефана-Больцмана

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

Температура газа находится из уравнения[39]:

, (4.6)

где - теплоёмкость воздуха;

- масса воздуха в емкости.



где

- конвективная составляющая теплообмена между жидкостью и теплоносителем;

где - площадь поверхности жидкости;

- температура теплоносителя;

- температура жидкости;

- лучистая составляющая теплообмена между жидкостью и теплоносителем;

где - степень черноты жидкости;

- конвективная составляющая теплообмена между жидкостью и газом;

где - коэффициент теплоотдачи газа, находится экспериментально;

- лучистая составляющая теплообмена между жидкостью и газом;

- конвективная составляющая теплообмена между жидкостью и пластиной;

где - температура пластины;

- коэффициент теплоотдачи жидкости, находится экспериментально;

- лучистая составляющая теплообмена между жидкостью и пластиной;

Температуру жидкости находим из уравнения:

, (4.7)

где - теплоёмкость жидкости;

- масса жидкости.



где

- конвективная составляющая теплообмена между пластиной и газом;

где - площадь поверхности пластины;

- лучистая составляющая теплообмена между пластиной и газом;

где - степень черноты пластины;

- конвективная составляющая теплообмена между жидкостью и пластиной;

- лучистая составляющая теплообмена между жидкостью и пластиной;

- конвективная составляющая теплообмена между пластиной и теплоносителем;

- лучистая составляющая теплообмена между пластиной и теплоносителем;

Температуру пластины находим из уравнения:

, (4.8)

где - теплоёмкость пластины;

- масса пластины.



где

- конвективная составляющая теплообмена между стеклянной поверхностью стенки и газом;

где - площадь стеклянной поверхности стенки;

- температура стеклянной поверхности стенки;

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

где - степень черноты стеклянной поверхности стенки;

- конвективная составляющая теплообмена между стеклянной поверхностью стенки и теплоносителя;

- лучистая составляющая теплообмена между стальной поверхностью стенки и теплоносителем;

- конвективная составляющая теплообмена между стальной поверхностью стенки и газом;

где - площадь стальной поверхности стенки;

- температура стальной поверхности стенки;

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

где - степень черноты стальной поверхности стенки;

- конвективная составляющая теплообмена между стальной поверхностью стенки и теплоносителя;

- лучистая составляющая теплообмена между стальной поверхностью стенки и теплоносителем;

- лучистая составляющая теплообмена с стальной поверхности стенки с внешней средой.

- лучистая составляющая теплообмена с стеклянной стенки с внешней средой.

Температура стеклянной стенки находим из уравнения:

, (4.9)

где

- теплоёмкость стеклянной стенки;

- масса стеклянной стенки.

Температура стальной стенки находим из уравнения:

, (4.10)

где - теплоёмкость стельной стенки;

- масса стальной пластины.

Начальные условия для жидкости принимаются в виде «зеркала».

Значения теплофизических констант для ЭМУ принимаются в следующими:

- температура теплоносителя 373-423 К;

- массовый секундный расход 100-300 л/мин;

- температура в ЭММУ 293°К;

- масса газифицируемой жидкости 50 гр.

На рисунке 4.5 приведена блок-схема алгоритма расчёта параметров процесса испарения жидкости.
z:\kazakov\схема.jpg

Рисунок 4.5. Блок-схема алгоритма расчета параметров процесса испарения жидкости.

В приведённой системе уравнений проблемными вопросами является определение коэффициентов теплоотдачи, которые являются функциями конкретных условий взаимодействия и аналитическое определение их невозможно.
4.3 Физические модели АЦК

Физическая составляющая АЦК представляет собой набор измерительно-контрольной аппаратуры.

Таблица 4.1 Датчики контроля и регистрации

1. Датчик давления АИР-20М2

описание: безымянный

Максимальный верхний

предел измерений 600 кПА

Предел перегрузки

(предельно допускаемое

рабочее давление) 2,5 МПА

Основная погрешность от 0,075 %



2. Датчик температуры ТС1088



Диапазон измерений -50...+200 °С

Длина монтажной части L 60 мм

Условное давление 6,3 МПа

Показатель тепловой инерции 30 с



3. Измеритель-регулятор технологический ИРТ 5920Н





Диапазон измерений -50...+200 °С
Климатическое

исполнение -30...+65 °С

Время установления

рабочего режима не более 30 мин.

Пределы допускаемой основной приведенной погрешности
относительно НСХ ±0,2 %


4. Датчик расхода FESTO SFAB – 600U - 2SA

описание: c:\documents and settings\admin\мои документы\загрузки\sfab_6738u_1_240px.jpg


Диапазон измерений 6…600 л/м
Максимальное давление 6 атм
Предельная погрешность 0,3 %
Рабочая температура -20…+80 °С




5.Аналогово-цифровой преобразователь B – 480G

описание: c:\documents and settings\admin\мои документы\загрузки\b480.jpg

Число измерительных каналов - 8

Разрядность - 16 бит

Быстродействие –

250 Квыб/с (R ист. сигнала не более 1 кОм)

100 Квыб/с (R ист. сигнала не более 1 MОм)
Диапазон входного сигнала ± 10 В


7. Персональный компьютер

описание: computer_014

CPU Pentium 4

ОЗУ 1024

HDD 40Gb

ОС Windows XP

Установленный пакет программ Statistica 8.0

Система обработки данных АЦП ADC 2.3.8




Рисунок 4.6 - Схема подключения датчиков и управления пневмосистемой
Программное обеспечение для АЦК используется как стороннее так уникальные программные продукты разработанные специально для данной системы.

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

4.4 Выводы по главе 4

На основании проведённых в главе исследований получены следующие основные результаты

1. Рассмотрен состав и назначение аналого-цифрового комплекса для параллельного математического моделирования процессов, исследуемых на стенде. Показаны возможности АЦК для решения задач валидации и верификации исследуемых на стенде процессов.

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

Введены критерии, расчёт которых подтверждает (или опровергает) достоверность получаемых результатов.

3. Представлена физическая модель АЦК в составе экспериментального стенда. Сформированы основные требования к функционированию АЦК, в частности, разработки алгоритма синхронизация математического и физического времени.

4. Показана возможность повышения достоверности и точности получаемых результатов математического моделирования за счёт использования в математическом моделировании физических измерений, учитывающих нелинейные эффекты, которые затруднено описать в упрощённой математической модели, например, углы ввода теплоносителя, введение ультразвукового акустического взаимодействия с различной ориентацией зон воздействия и т.д.
1   2   3   4   5   6   7   8   9   10

Похожие:

Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Федеральное государственное бюджетное образовательное учреждение...
Министерство образования и науки российской федерации федеральное государственное бюджетное образовательное учреждение высшего профессионального...
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Федеральное государственное бюджетное образовательное учреждение...
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Федеральное государственное бюджетное образовательное учреждение...
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Федеральное государственное бюджетное образовательное учреждение...
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Федеральное государственное бюджетное образовательное учреждение...
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Отчет о самообследовании деятельности федерального государственного...
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Отчет о самообследовании деятельности федерального государственного...
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Федеральное агентство воздушного транспорта
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Московский государственный...
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon «донской государственный технический университет» (дгту) Кафедра «Иностранные языки»
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Программ а стратегического развития Федерального государственного...
Полное наименование вуза – федеральное государственное бюджетное образовательное учреждение высшего профессионального образования...
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Федеральное государственное образовательное учреждение высшего профессионального...
Федеральное государственное образовательное учреждение высшего профессионального образования
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Рекомендации по приведению в соответствие мероприятий и планов реализации
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Нижегородский государственный...
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Федеральное государственное бюджетное образовательное учреждение...
Сущность и этапы организационно управленческого анализа (на примере ОАО "вымпелком")
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Самарский государственный технический университет
Федеральное государственное бюджетное образовательное учреждение высшего образования
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Федеральное государственное бюджетное образовательное учреждение...
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования
Рф федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Омский государственный технический университет» icon Федеральное государственное бюджетное образовательное учреждение...
Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования

Руководство, инструкция по применению




При копировании материала укажите ссылку © 2024
контакты
rykovodstvo.ru
Поиск