Моделирование структуры ЭГД-течений в несимметричной системе электродов, страница 4

DX=0.02                       Ширина кюветы 

DY=0.01                       Межэлектродное расстояние

DY1=DY+DY                    

r=0.00035                     Радиус электрода

AX=0.001                     Размеры центральной струи

k,1,0,0                           Построение начинается с опорных точек, далее по ним строятся линии,

k,2,-dx,-dy                    по линиям строятся области

k,3,dx,-dy

k,4,dx,dy1

k,5,-dx,dy1                    

circle,1,r,,,,8                  Построение электрода, параметр 8 определяет из скольких линий

lcomb,1,8,0                   строится окружность. Это необходимо для моделирования

lcomb,2,3,0                   центральной струи

lcomb,4,5,0

lcomb,6,7,0

l,5,9

l,4,7

l,2,11

l,3,13

l,3,4

l,4,5

L,2,5

l,2,3

al,1,8,9,5

al,6,7,12,8

al,4,7,11,3

al,2,3,10,5

Построение сетки.

lsel,s,,,3,7,2                         Выбираются линии с 3 по 7 с шагом 2

lsel,a,,,8                               Дополнительно выделяется линия 8

lesize,all,,,35,0.1                 Выбранные линии разбиваются на 35 частей с соотношением 0.1

lsel,all                                  Далее разбиваются остальные линии

lesize,6,,,20,-0.3

lesize,12,,,20

lesize,4,,,10,-0.3

lesize,2,,,10,-0.3

lesize,1,,,10,-0.3

lsel,s,,,9,11,1

lesize,all,,,10

lsel,all

mshkey,1                            Строится упорядоченная сеть из прямоугольников.

amesh,all

После этого база данных сохраняется.

2. Электростатическая задача. В этой части определяются напряжение на электродах, распределение объемного заряда, и решается электростатическая задача. Ниже описана программа VCGunif, которая совершает данные действия.

et,2,plane121                           Выбор другого элемента (плоский электростатический Plane 121)

esel,all

EMODIF,all,TYPE,2               Замена элементов Mesh200 на Plane121

V1=0                                        Задается потенциал на плоском электроде

V2=8000                                  Задается потенциал на цилиндрическом электроде

mp,perx,1,2.2                           Значение относительной диэлектрической проницаемости

                                                 (perx) задается равным 2,2 для материала под номером 1.

mat,1                                        Всем элементам присваивается материал под номером 1.

*get,nodemx,node,0,count       Определяется максимальный номер узла в выделенном наборе и

                                                  и сохраняется в виде параметра с названием nodemx.

*dim,CHG,array,nodemx         Создается массив с названием CHG, тип массива – array,

                                                  размерность массива - nodemx´1´1                               

NSEL,S,loc,Y,-DY                   К плоскому электроду прикладывается значение потенциала

D,all,volt,V1                             равное V1.

LSEL,S,,,2,6,2                          

lsel,a,,,1                                     

nsll,s,1                                       К цилиндрическому электроду прикладывается значение

D,all,volt,V2                             потенциала равное V2.

z=v2*(1E-7)                             Начальное значение заряда определяется через напряжение на

                                                  цилиндрическом электроде.

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