Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

3340

.pdf
Скачиваний:
0
Добавлен:
15.11.2022
Размер:
4.32 Mб
Скачать

Расчет индивидуального рейтинга

Ранжирование рейтингов по профилям кафедр и категориям преподавателей

Оценка результатов работы преподавателей

Разработка мероприятий по стимулированию работы ППС, кафедр, факультетов

Рис. 2

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

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

Литература

1.Макаров Г.В. Положение о рейтинговой оценке деятельности факультетов и кафедр ВГТУ/ Г.В. Макаров, Э.Д. Поликарпов, Т.А. Бурковская и др. Воронеж: Изд-во ВГТУ, 2000. - 17 с.

2.Макаров Г.В. Положение о рейтинговой оценке деятельности профессорско-преподавательского состава ВГТУ/ Г.В. Макаров, Э.Д. По-ликарпов, Т.А. Бурковская и др. Воронеж: Изд-во ВГТУ. 2000. - 26 с.

3.Поликарпов Э.Д., Халявина А.В. Оценка личного рейтинга ППС ВГТУ в 1999 году. /Сборник научных трудов «Высокие технологии в

технике, медицине, экономики и образования». Часть 1. Воронеж, 2001. – 200 с.

Получено 6.05.01

Воронежский государственный

 

технический университет

153

УДК 621.313

С. А. Горемыкин, Ю. В. Писаревский ИСПОЛЬЗОВАНИЕ ЧИСЛЕННОГО МЕТОДА И МЕТОДА ЭКВИВАЛЕНТНЫХ ТЕПЛОВЫХ СХЕМ ЗАМЕЩЕНИЯ ДЛЯ РАСЧЕТА ТЕПЛОВЫХ ПРОЦЕССОВ В ЭЛЕКТРИЧЕСКИХ МАШИНАХ МАЛОЙ МОЩНОСТИ

В работе рассмотрены вопросы расчета температуры в электродвигателе постоянного тока с полым якорем с

использованием метода эквивалентных тепловых схем замещения и программы Quick Field.

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

метода ЭТСЗ, тепловая система с непрерывно распределенными параметрами замещения эквивалентной системой однородных тел (материал, условия выделения потерь, характер контакта с соседними элементами или средой). К недостаткам данного метода следует отнести сложности, связанные с разбиением электрической машины на однородные тела (т.е. число узлов схемы замещения, количество конструктивных элементов, входящих в каждый узел).

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

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

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

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

Выполним расчет в программе Quick Field для стационарного теплового процесса. Картина распределения температурного поля в электродвигателе ДПР 52 для номинального режима представлена на

рис.1.

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

154

Рис. 1. Распределение температурного поля в электродвигателе ДПР – 52

ниях учитываются анизотропные в тепловом отношении свойства материалов. Так в таблице 1 приводятся результаты расчѐта радиальной p

и аксиальной

a

составляющих коэффициента теплопроводности обмотки

 

 

якоря, в зависимости от коэффициента заполнения обмотки по меди.

 

 

 

 

 

Таблица 1

 

 

 

 

 

 

p , Вт м K

0,401

0,435

0,499

0,582

0,754

 

 

 

 

 

 

a , Вт м K

60,4

108,4

144,4

216,4

276,4

 

 

 

 

 

 

кз

0,14

0,2

0,29

0,4

0,5

 

 

 

 

 

 

Таким образом, в общем случае Quick Field позволяет учесть:

а) нелинейную зависимость термического сопротивления от температуры;

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

В итоге, на основе анализа теплового поля ДПР, рассчитанного в программе Quick Field для стационарного режима, была уточнена упрощенная тепловая схема замещения, указанная в /3/. При детальном анализе распределения теплового поля были конкретизированы узлы тепловой схемы, коэффициенты теплопроводности неоднородных материалов, а также геометрические размеры, связывающие отдельные блоки в тепловую сеть. Схема замещения для расчѐта стационарного теплового процесса изображена на рисунке 2.

155

 

 

 

 

 

 

 

R

а

Р

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я

 

 

к

Р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

R

4

R

2

R

 

 

R

 

 

 

 

 

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

Рис. 2

 

 

 

На рисунке 2 приняты следующие обозначения:

Кл - узел, учитывающий температуру коллектора; ющий активную и лобовую часть якоря;

К - узел, включающий корпус, индуктор, подшипниковые щиты;

Pa 0,83Iа2 R a - электрические потери в данной части обмотки

якоря;

Рk , Вт - потери в коллекторе

P

0,17 I2

R

a

I

a

U

щ

k

тр

P

S

щ

кол

m ,

k

a

 

 

 

 

щ

 

 

где Ia , А — ток в обмотке якоря;

R a , Ом — активное сопротивление обмотки якоря;

U щ , В — переходное падение напряжения на пару щѐток; kтр = 0,2 — коэффициент трения щѐток о коллектор;

Pщ (6 10) , Нсм2 — удельное давление щѐтки на коллектор; Sщ , см2 — сечение щѐтки;

m — число щѐток;

кол , мс — окружная скорость коллектора.

156

Таблица 2

Название и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

обозначение

 

Расчётная формула

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

термического

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сопротивления

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1. R1 — между якорем

R1 R 2

 

 

 

 

1

 

 

 

 

 

hоб

 

 

 

 

 

 

 

 

 

 

 

и корпусом

 

 

 

2

 

 

p Dа,с

la

 

lл

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в,з

 

 

 

Dа,с

la

lл

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.

R 3

между

R

3

 

 

 

 

 

 

 

 

4

la

 

 

 

 

 

 

 

 

0,2 lоб

 

 

 

коллектором и якорем

 

 

 

 

 

 

 

 

 

Dа2

Dа2

 

 

 

 

 

 

 

qпр

S

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

м

 

 

3.

 

 

R 4

R

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

теплопередачи

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в,з

 

 

 

 

 

Dи,н

Dи,в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

воздушного зазора между

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

обмоткодержателем

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

индуктором

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.

R 5 — учитывающее

R

 

 

 

1

 

hщ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

теплоотвод через щетки

5

 

 

 

 

Sщ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

щ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.

 

 

R 6

R 6

1

 

 

 

 

 

 

0,75 dк,н

 

dв,ц

 

 

 

 

 

 

сопротивление,

 

 

 

пл

 

 

 

 

0,75 dк,н

 

 

 

dв,ц

 

 

lk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

учитывающее теплоотвод

 

 

r1

 

r3

 

 

 

 

r2

 

r3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

через

два

выходных

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r1

 

 

r2

2

r3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

конца вала

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

2

 

 

 

 

 

 

 

 

 

 

 

d2

 

 

 

 

 

 

где r

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

в

,

 

r

 

l

 

 

 

 

в

,

 

 

 

 

 

 

 

в,и

 

 

 

в

 

 

 

 

 

в,к в

 

 

 

 

 

 

 

1

 

 

 

 

 

 

4

 

 

 

 

 

2

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r3

 

1 ( т,в

 

 

 

dв lв )

 

 

 

 

 

 

 

 

 

 

 

6.

R 7 — теплоотдачи

R 7

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

наружной

поверхности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т,нп

 

 

 

Dк,н Lп

 

 

 

Dк,н

 

 

двигателя

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В таблице 2 геометрические размеры отдельных частей машины и

тепловые коэффициенты обозначены следующим образом:

 

 

 

в,з

136 Вт м2 K - коэффициент теплопередачи воздушного зазора;

м390 Втм K - коэффициент теплопроводности меди;

щ100 Втм K - коэффициент теплопроводности щѐток;

пл

в

0,2 Втм K - коэффициент теплопроводности пластмассы; 16 Втм K - коэффициент теплопроводности вала;

157

т,в

79 Вт м2

K

-

коэффициент

теплоотдачи

наружной

 

 

 

 

 

 

поверхности валов;

 

 

 

 

 

т,нп

13 Вт м2

K

-

коэффициент

теплоотдачи

наружной

 

 

 

 

 

 

поверхности машины;

 

 

 

 

hоб , м - толщина обмотки якоря;

 

 

Dа,н , м - наружный диаметр обмотки якоря;

 

Dа,в , м - внутренний диаметр обмотки якоря;

 

Dа,c , м - средний диаметр обмотки якоря; Dи,н , м - наружный диаметр индуктора; Dи,в , м - внутренний диаметр индуктора;

lа , м - длина активной части обмотки якоря; lл , м - длина лобовой части обмотки якоря;

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

S - число секций обмотки якоря;

qпр , м 2 - сечение обмоточного провода; h щ , м - высота щѐтки;

Sщ , м 2 - сечение щѐтки;

dк,н , м - наружный диаметр коллектора; dв,ц , м - диаметр центральной части вала; lк , м - аксиальная длина коллектора;

lв,к , м - длина вала от центра коллектора до центра шпонки со

стороны коллектора;

lв,и , м - длина вала от центра коллектора до центра шпонки со

стороны индуктора;

lв , м - длина вылета выходного конца вала; dв , м - диаметр выходного конца вала;

Dк,н , м - наружный диаметр корпуса; Lп , м - полная длина корпуса.

Таким образом, использование программы Quick Field, позволяет более точно рассчитать плоскую тепловую задачу. А также является перспективным использование этой программы для уточнения существующих методов теплового расчѐта (тепловых схем замещения).

Литература

1. Борисенко А. И. и др. Охлаждение промышленных электрических машин // А. И. Борисенко, О. Н. Костиков, А. И. Яковлев. — М.:

Энергоатомиздат, 1983. — 296 с., ил.

158

2. Эквивалентные тепловые схемы замещения двигателя постоянного тока с полым ротором / Горемыкин С. А., Кононенко Е. В., Писаревский Ю. В., Фурсов В. Б.; ВГТУ. - Воронеж, 2000. Деп. в ВИНИТИ.

Получено 17.09.01

Воронежский государственный

 

технический университет

УДК 621.313.333

159

Г. В. Коробов

ОПРЕДЕЛЕНИЕ ОПТИМАЛЬНЫХ ПАРАМЕТРОВ ОДНОФАЗНОГО АСИНХРОННОГО ЭЛЕКТРОДВИГАТЕЛЯ С ЭКРАНИРОВАННЫМИ ПОЛЮСАМИ

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

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

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

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

Воснову одного из них [2] положено использование аппарата математического программирования, где вводится понятие функции цели (критерия) и фактора (параметра). Для АДЭ целесообразно совместное применение метода планирования эксперимента и штрафных функций [3].

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

пусковой ток Iп и коэффициент полезного действия при номинальном скольжении sн 0,16.

В

качестве переменных факторов были выбраны: величина воздушного зазора

, длина

мостика насыщения l, ширина мостика насыщения g, угол экранирования

и угол смещения оси

шунтов

(рис. 1).

 

 

160

g

Рис. 1. Варьируемые параметры статора

Для исследуемого электродвигателя ДА–60–05–3 были получены следующие зависимости.

M n

2,15 10

 

3

5,18 10

5

 

 

6,2 10 5

 

 

1,07

5,31 10

4 l

 

 

 

1,32 10

3

g

2,26 10

 

9

 

 

 

8,34 10 2

 

 

3,98 10

6

l

 

 

 

2,25 10

5

 

 

g

1,61 10

3

 

 

3 10

8

 

 

l

1,13 10 7

 

g

 

 

 

1,28 10 1

 

 

l

3,65 10 1

 

 

g

5,51 10

5

l

g

2,51 10

7

2

 

 

 

1,32 10

6

 

 

2

4,96 103

 

2

1,09 10

4

l 2

 

8,34 10

4 g 2 ;

 

 

 

 

I п

1,1 10 1

 

1,56 10 4

 

 

 

4,27 10 4

 

2,5 102

5,08 10 3

l

 

 

5,4 10 2

g

 

2,41 10

2

 

 

 

4,17 10 6

 

l

2,8 10 5

l g

 

 

(1)

 

 

 

 

 

 

 

 

3 10 3

 

 

 

 

2

 

6,28 10 6

2

 

 

 

2,75

l

4,17

g

 

l

g

1,11

 

 

 

 

 

5,65 105

 

2

7,94 10 5

l 2

2,17 10

2

g 2 ;

 

 

 

 

 

 

 

 

 

6,74

10

2

1,26

10

3

 

 

5,62

10

4

 

 

3,83 10

3

 

3,32

10

3

l

2,28 10

2

 

g

3.53 10

7

 

 

 

1,58

 

 

4,78 10

5

l

 

 

 

 

2,68 10

4

 

 

g 1.08 10 1

 

 

 

3,78 10

6

 

l 2,13 10

5

g

 

 

4,71 10 1

 

 

l

8,92

 

g

2,43 10 3

l

g

 

8,7 10

6

2

 

 

 

 

1,02 10

5

 

 

2

1,28 106

 

2

1,78 10

3

l 2

6,63 10 3

g 2 .

 

 

 

Анализ зависимостей (1) показал, что улучшение показателей АДЭ возможно за счѐт изменения

,.

Зависимости (1) использовались для решения следующих задач квадратичного

программирования: определить

в

заданном

пятимерном

факторном

пространстве

координаты

, , , l, g , максимизирующие

функции

M п F1 (

, , , l, g) ,

F3 ( ,

, , l, g) ;

минимизирующие - I п F2 ( ,

,

, l, g) при линейных ограничениях на факторы ,

, ,l, g .

161

При решении задач оптимизации использовался метод штрафных функций. Это позволило перейти от оптимизации функций с ограничениями (F1,F2,F3) к оптимизации функции без ограничений ( 1, 2, 3), что значительно упростило поиск решения.

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

 

 

u

1

 

2 ( , , ,l, g,r) F2 ( , , ,l, g) r0

;

 

 

 

 

j 1 c j ( , , ,l, g)

(2)

min

2 (

, , ,l, g,r0 ).

 

 

( , , ,l,g ) D

 

 

где c j (

,

, ,l, g) - ограничения, накладываемые на факторы или функцию цели; u- число

ограничений, r0- целое положительное число; D- допустимая область определения факторов.

При r0

0

минимум

функции

2 ( , , ,l, g,r0 ) без ограничений совпадает с

минимумом функции I п

F2 ( ,

, , l, g)

с ограничениями.

На основе этого метода был разработан алгоритм оптимизации зоны шунтов АДЭ, который реализован в виде программы на алгоритмическом языке Turbo Pascal 7. Результаты оптимизации приведены в таблице.

 

 

 

 

 

 

Таблица

Показатели

M п

 

I п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Н м

 

A

 

 

%

Базовый АДЭ

2,45 10

3

54,3 10

3

7,6

 

 

 

 

 

 

 

 

 

 

 

АДЭ с максимальным

2,59 10

3

57,8 10

3

7,61

пусковым моментом

 

 

 

 

 

 

 

 

 

АДЭ с минимальным

2,03 10

3

41 10

3

 

7,1

пусковым током

 

 

 

 

 

 

 

 

 

 

АДЭ с максимальным

2,57 10

3

56,2 10

3

8,7

КПД

 

 

 

 

 

 

 

 

 

Улучшение показателей асинхронных двигателей с экранированными полюсами после оптимизации по сравнению с базовым вариантом составило для пускового момента 5,7%, для пускового тока 24,5%, для коэффициента полезного действия 14,5%.

Результаты оптимизации подтверждают обоснованность применения разработанной методики для формирования конфигурации шунтов статора АДЭ. Это позволяет проектировщику улучшить один из интересующих его показателей (Mп, Iп, ), проведя дополнительную конструктивную проработку микродвигателя.

Литература 1. Кононенко Е.В., Пархоменко Г.А., Коробов Г.В. К анализу асинхронных двигателей с

экранированными полюсами

(статья) / Научно-практический

вестник ―Энергия‖, №3 (37), 1999,

г. Воронеж.

 

 

2. Ивботенко Б.А.,

Ильинский Н.Ф., Копылов И.П.

Планирование эксперимента в

электромеханике. – М.: Энергия, 1975. –185 с.

3. Коробов Г.В., Титова Л.Н. Определение оптимальных параметров электродвигателей с диагональной обмоткой полого якоря (статья) / Научно-практический вестник ―Энергия‖, №4 (42), 2000,

г. Воронеж.

 

Получено 25.10.01

Воронежский государственный

 

агроуниверситет

162

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]