тэмп метода. Практикум 1 Задание 1. Электростатическое поле коаксиального
 Скачать 1.35 Mb.
|
|
ТЭМП/ELCUT/Лабораторный практикум 1 Задание №1. Электростатическое поле коаксиального кабеля Целью работы является изучение влияния конструктивных параметров коаксиального кабеля на картину электростатического поля во внутренней области. 1. Сведения из теории Рис.1 Конструкция коаксиального кабеля Одножильный кабель на рабочее напряжение U 0 имеет радиус жилы r 1 и внутренний радиус оболочки r 2 (рис.1). В идеале оба проводника должны иметь общую ось, однако в процессе изготовления или эксплуатации возможно смещение проводников на расстояние d, что может привести к пробою изоляции кабеля, если максимальная напряженность электростатического поля превысит допустимое значение. Обычно при проектировании коаксиального кабеля задается радиус жилы r 1 , который выбирается по величине токовой нагрузки. Внешний радиус изоляции (внутренний радиус оболочки r 2 ) определяется исходя из допустимого значения максимальной напряженности электростатического поля и рабочего напряжения кабеля E max0 : ) ln( 1 2 1 0 0 max r r r U E = , откуда r 2 = r 1 ⋅ exp (U 0 /r 1 E max0 ) (1) Для уменьшения внешнего радиуса кабеля часто применяется изоляция из кабельной бумаги с различной диэлектрической проницаемостью слоев. При этом стараются обеспечить наименьшую неравномерность распределения напряженности электростатического поля вдоль радиуса. Существенного выравнивания поля можно добиться при выполнении определенных соотношений между диэлектрическими проницаемостями слоев и их радиусами. ТЭМП/ELCUT/Лабораторный практикум 2 Рис.2 Кабель с градированной изоляцией В частности, для кабеля с двухслойной изоляцией (рис.2) оптимальное распределение напряженности поля достигается при соблюдении равенства: ε 1 r 1 = ε 2 r 3 (2) Радиус оболочки r 2 в этом случае может быть найден из равенства: ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ ε ε + = ) ln( ) ln( 3 2 2 1 1 3 1 0 0 max r r r r r U E (3) При расчете электростатического поля, создаваемое такой системой проводников, принимают, что потенциал на поверхности внутреннего проводника (жилы) равен U 0 , а потенциал на поверхности наружного проводника равен нулю. При этом поле сосредоточено в области диэлектрика, заполняющего кабель, и во внешнее пространство не проникает. Здесь учитывается и тот факт, что электростатическое поле в проводник не проникает, поэтому поверхности проводников становятся границами области существования поля. Учитывая, что длина кабеля велика по сравнению с его диаметром, поле можно считать плоскопараллельным, т.е. не изменяющимся вдоль оси проводников. В однородном изотропном диэлектрике ( ε = const) такое поле описывается двухмерным уравнением Лапласа для электростатического потенциала U: 0 2 2 2 2 = ∂ ∂ + ∂ ∂ y U x U (4) В общем случае электростатическое поле для плоскопараллельной задачи описывается уравнением Пуассона: ρ − = ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ ∂ ∂ ε ∂ ∂ + ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ ∂ ∂ ε ∂ ∂ y U y x U x y x , (5) где ρ - объемная плотность заряда в рассматриваемой области, а диэлектрические свойства различны вдоль осей координат. Следовательно, для того, чтобы уравнение (4) и (5) были эквивалентны, в (5) необходимо положить ρ = 0 и указать изотропность свойств диэлектрика ( ε x = ε y = ε ). Решение уравнения Лапласа может быть получено при учете условий на границах рассматриваемой области. Применительно к поставленной задаче следует задать потенциалы на поверхности проводящих цилиндров. ТЭМП/ELCUT/Лабораторный практикум 3 2. Расчетная модель задачи Перед началом работы в пакете ELCUT следует разработать расчетную модель решаемой задачи. При моделировании обычно стараются выделить часть конструкции, используя наличие симметрии в расположении элементов конструкции и приложенных нагрузок. Это позволяет повысить точность расчетов. В частности для анализируемой задачи в качестве оси симметрии может быть рассмотрена горизонтальная ось. Условие симметрии означает равенство нулю производной от электростатического потенциала по направлению нормали к поверхности ( 0 = ∂ ∂ n U ). Рис.3 Расчетная модель В результате расчетная модель приобретет вид, изображенный на рис.3. Такая модель соответствует начальному расположению элементов. Если по условию задачи задан однослойный диэлектрик, то полуокружность радиуса r 3 прорисовывать не следует. Требуемые изменения модели после проведенного базового расчета будут проводиться путем смещения соответствующих элементов или изменения их радиусов. Смещение жилы предполагается осуществлять вдоль оси Ox , а изменение радиуса обеспечивать смещением узловых точек. 3. Этапы решения задачи 1. Создание новой, пустой задачи ELCUT и ввод параметров задачи 2. Создание геометрической модели (рисунка расчетной области) и меток объектов 3. Описание физических свойств материалов и ввод граничных условий 4. Построение сетки конечных элементов 5. Решение задачи; 6. Обработка результатов решения Внимание!! Предварительно на студенческом диске (диск H ) должна быть создана личная папка студента, в которой следует сохранять все файлы ELCUT 4. Инструкция к выполнению работы 1. При создании новой задачи ELCUT ( «Файл» → «Создать» → «Задача ELCUT» ) в соответствующих диалоговых окнах указать: имя задачи – task1 ; тип задачи – электростатика ; класс модели – плоская ; единицы длины – миллиметры ; система координат – декартова . В результате на рабочей панели слева появится «Окно задачи» 2. Построение модели: Двойной клик мыши по разделу «Геометрия» в «Окне задачи» откроет окно геометрического редактора – появится координатная сетка, на которой красным цветом нанесены координатные оси (горизонтальная – Ох , вертикальная – Оу ). Следует установить шаг сетки равным 1 мм ( «Правка» → «Сетка привязки» ). ТЭМП/ELCUT/Лабораторный практикум 4 2.1. Прорисовка расчетной области: Перейти в режим «Вставка объектов» и используя инструменты «Половина круга» и «Прямая линия» изобразить расчетную модель заданных начальных размеров. Выбрав нужный инструмент, курсор помещается в точку с начальными координатами линии (координаты курсора отображаются справа в нижней строке «Окна задачи» ) и далее при нажатой левой кнопкой мыши перемещается в точку с координатами конца линии. 2.2. Ввод меток объектов: Перейти в режим «Выделение объектов» и присвоить имена: блокам – замкнутым геометрическим областям, обладающим материальными свойствами (для данной задачи диэлектрической проницаемостью); ребрам – линиям, ограничивающим расчетную область, на которых будут заданы граничные условия. Клик правой кнопки мыши по объекту выделяет его красным цветом и вызывает контекстное меню, в котором следует выбрать пункт «Свойства» . В появившемся диалоговом окне задать имя объекта. 3. Описание физических свойств материалов: Двойной клик мыши по имени блока в «Окне задачи» вызывает диалоговое окно для ввода свойств. Ввести заданное значение для относительной диэлектрической проницаемости. 4. Задание граничных условий Двойной клик мыши по метке ребра в «Окне задачи» вызывает диалоговое окно для ввода свойств. Соответствующие разным границам условия показаны на рис.4. Рис.4 Граничные условия Ребро b 1 b – U = 0; ребро a 1 a – U = U 0 (дана в условии задачи); ребра b 1 a 1 и a b – σ = 0. 5. Построение сетки конечных элементов Для построения сетки конечных элементов предварительно необходимо задать шаг дискретизации (разбиения) расчетной области. Шаг дискретизации определит густоту сетки конечных элементов в различных областях модели. Для рассматриваемой задачи можно предположить существенную неоднородность поля вблизи жилы, поэтому здесь сетка должна быть гуще. Чтобы задать значение шага дискретизации следует: • Двойным кликом мыши выделить ребро и выбрать команду «Свойства» в контекстном меню; • Установить переключатель «Шаг дискретизации» в положение «ручной» и ввести нужное число (для ребра a 1 a – 0.1 мм (в зад.№3 – 0.02 мм), для ребра b 1 b – 1 мм). Для построения сетки воспользоваться файловым меню ( «Правка» → «Построить сетку» → «Во всех блоках» ) или соответствующей кнопкой на панели инструментов. 6. Решение задачи: «Задача» → «Решить задачу» ТЭМП/ELCUT/Лабораторный практикум 5 После окончания расчета в правой части основного окна будет выведена рассчитанная картина поля текущей задачи 7. Анализ результатов: 7.1 Построить зависимость напряженности поля вдоль контура a b : 1.1.1. «Контур» → «Добавить (Линия/Ребро/Блок)» → клик левой кнопки мыши по ребру. При необходимости сменить направление контура, чтобы он шел от узла a к узлу b ( «Контур» → «Сменить направление» ). 1.1.2. «Вид» → «График» и из предложенного набора величин выбрать «Напряженность» 1.1.3. График распечатать: «Файл» → «Печать» 7.2 Определить максимальное значение напряженности (соответствует узлу a ) и занести в таблицу результатов 1.2.1. «Вид» → «Локальные значения» . В результате появляется новое окно. 1.2.2. Кликом правой кнопки мыши по строке «Укажите точку» вызвать контекстное меню и выбрать пункт «Координаты точки» 1.2.3. В предложенном поле ввести координаты узла a 1.2.4. Полученное значение записать в таблицу. 7.3 Закрыть окно «Анализ результатов» 8. В соответствии с программой работы произвести необходимые изменения в геометрии задачи, используя возможности геометрического редактора (смещение объектов): 8.1 Левой кнопкой мыши выделить объект – ребро или узел (ребро a 1 a при изучении влияния смещения жилы, узел a ( a 1 ) при изменении радиуса жилы, узел b ( b 1 ) при изменении радиуса обмотки). 8.2 «Правка» → «Передвинуть выделенное» . В появившемся диалоговом окне указать координаты вектора переноса (смещение производится только по оси Ох !). 8.3 При изменении радиуса жилы или оболочки п.8.2 повторяется дважды (для каждого из узлов), при этом вектора переноса антисимметричны, например (2, 0) и ( −2, 0). 9. Построить сетку (при этом не надо менять шаг дискретизации, а просто воспользоваться файловым меню – см.п.5), решить задачу и выполнить пункт 7.2. 10. Пункты 8 и 9 повторить необходимое число раз. Для найденного окончательного варианта построить зависимость напряженности поля вдоль контура a b в соответствии с п.7.1 и распечатать картину эквипотенциальных линий. 5. Задание и содержание отчета Задача №1. Определить, на какую величину может быть допущено смещение оси жилы кабеля по отношению к оси оболочки, чтобы максимальная напряженность поля в диэлектрике не превысила двойной по сравнению с максимальной напряженностью при совпадении осей. Изоляцию кабеля считать однородной. Построить зависимость E max (d) (результаты расчета E max для различных d должны быть занесены в таблицу). Для построения использовать не менее 8 ÷10 точек. Привести зависимости напряженности поля вдоль контура a b для исходной и конечной конфигурации. Привести картину эквипотенциальных линий для конечной модели. ТЭМП/ELCUT/Лабораторный практикум 6 № U 0 , кВ ε r r 1 , мм r 2 , мм 1 3 2 5 20 2 3 2.2 6 22 3 6 2.2 7 25 4 6 2.4 8 25 5 10 2.4 9 30 6 10 3 10 30 Задача №2. При заданном размере оболочки кабеля, рассчитанного на U 0 = 3 кВ, определить оптимальный радиус жилы (соответствует минимальному значению максимальной напряженности поля в конструкции). Результирующее значение найти с точностью до десятых долей миллиметра. Исходное значение радиуса жилы приведено в таблице. Построить зависимость E max (r 1 ) (результаты расчета E max для различных r 1 должны быть занесены в таблицу). Для построения использовать не менее 8 ÷10 точек. Привести зависимости напряженности поля вдоль контура a b для исходной и оптимальной конструкции и картину эквипотенциальных линий для оптимальной модели. № ε r r 1 , мм r 2 , мм 1 2 10 15 2 2 11 17 3 2.2 12 20 4 2.2 13 22 5 2.4 14 23 6 2.4 15 25 Задача №3. Исследовать, как изменится радиус оболочки кабеля при использовании двухслойного диэлектрика. Вначале по формуле (2) рассчитать радиус внутреннего слоя изоляции (r 3 ), а потом изменять радиус оболочки так, чтобы напряженность не превысила E max0 = 2.5 ⋅10 6 В/м. Базовый вариант – кабель с однородной изоляцией ε r2 Построить зависимость E max (r 2 ) (результаты расчета E max для различных r 2 должны быть занесены в таблицу). Для построения использовать не менее 6 ÷8 точек. Привести зависимости напряженности поля вдоль контура a b для базового и оптимального варианта и картину эквипотенциальных линий для оптимальной модели. № ε r1 ε r2 r 1 , мм r 2 , мм 1 3.4 2 2 14 2 4.2 2 2.5 12.5 3 3.6 2 3 11.5 4 4.1 2 3.5 11.5 5 3.2 2 4 11 В этой задаче при прорисовке геометрии модели вводятся три полуокружности радиусами r 1 , r 2 и r 3 (r 3 находят по формуле (2)). В результате получают 2 замкнутые области (между r 1 и r 3 , и между r 3 и r 2 ), они обе должны быть проименованы. В предварительном расчете их свойства одинаковы и определены ε r2 . В дальнейших расчетах свойства ближайшего к жиле блока изменяют на ε r1 ТЭМП/ELCUT/Лабораторный практикум 1 Задание №2. Проводники в электростатическом поле Целью работы является изучение влияния тонких проводников на выравнивание на- пряженности электростатического поля на примере проходного изолятора. 1. Сведения из теории Рис.1. Эскиз проходного изолятора: 1 – токоведущий стержень; 2 – заземленный фланец; 3 – изоляция (бумажно-масляная или бакелитовая) Для ввода высокого напряжения внутрь металлических баков высоковольт- ных трансформаторов, масляных выключа- телей, силовых конденсаторов и других ви- дов оборудования высокого напряжения ис- пользуются проходные изоляторы. Конст- рукцию проходных изоляторов (рис.1) отли- чает весьма неблагоприятное расположение электродов, приводящее к крайне неравно- мерному распределению напряженности электростатического поля (как радиальной, так и аксиальной компоненты). В первом приближении можно считать, что заземлен- ный фланец и токоведущий стержень обра- зуют цилиндрический конденсатор, для ко- торого свойственна неравномерность рас- пределения напряженности поля по радиусу. В соответствии с этим радиальная напря- женность поля вблизи стержня значительно превышает напряженность поля вблизи фланца. Если не учитывать влияние краев, то радиальную компоненту напряженности можно определить по формуле: c r r r r U E 0 ln = Кроме того, к неравномерности поля (аксиальной компоненты напряженности, направ- ленной вдоль поверхности изоляции) приводит отличие продольных размеров обкладок кон- денсатора. В результате, наибольшая напряженность поля имеет место у края фланца, вблизи точки А. Для обеспечения хорошей электрической прочности изоляции необходимо сгладить не- равномерность распределения напряженности. Для достижения более равномерного распределения радиальной компоненты напря- женности изоляцию делают расслоенной: между слоями диэлектрика равной толщины про- кладывают тонкие металлические цилиндры из алюминиевой фольги толщиной 10 ÷20 мкм. В результате получается как бы несколько вложенных друг в друга цилиндрических конденсато- ров. Чтобы получить приблизительно равномерное распределение радиальной компоненты напряженности поля во всех слоях изоляции, размеры цилиндров (радиус и длина) должны обеспечивать равные заряды на каждом конденсаторе. ТЭМП/ELCUT/Лабораторный практикум 2 Рис.2. Эскиз проходного изолятора конденса- торного типа с проводящими обкладками Элементарный конденсатор, обра- зуемый цилиндрическим слоем длиной h и радиусами r и r+ ∆r, имеет емкость (рис.2): r rh C ∆ πε = 2 Заряд этого конденсатора будет r E rh r U rh U C q ⋅ πε − ≈ ∆ ∆ ⋅ πε = ∆ ⋅ = 2 2 , где r U E r ∆ ∆ − ≈ Для выполнения условий E r = const и q = const радиусы и длины конденсатор- ных обкладок должны подчиняться равен- ству r ⋅h = const (1) Для улучшения распределения на- пряженности электростатического поля вблизи краев обкладок, например в области |