Физика \ Электрогидродинамика

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

nsel,s,loc,x,-ax,ax

nsel,r,loc,y,-0.0005,-dy

Определяется область в которой задается объемный заряд.

*get,N1,node,0,num,min

*get,N2,node,0,num,max

Определяется максимальный и минимальный номер узлов принадлежащих области. Формат команд: параметру N1 присваивается значение равное максимальному (ключ max) (минимальному) номеру (ключ num) узла (ключ node) из выделенного набора узлов.

*do,N,N1,N2,1

Создается цикл по номерам узлов от N1 до N2.                                                                        N – переменная цикла.                                 

*get,X1,node,N,LOC,X

*get,Y1,node,N,loc,y

Определяются Х и Y координаты узла с номером N. Формат команды: параметру X1(Y1) присваивается значение равное X(Y) координате (ключ LOC) узла (ключ NODE) с номером N.

*if,Х1,Lt,-ax,or,Х1,gt,ax,cycle

*if,y1,gt,-0.0005,or,y1,lt,(-0.99*dy),cycle

По координатам определяется, принадлежит ли узел   нужной области. Если нет – переход к следующему                                                                       узлу. (Если параметр Х1 меньше параметра –ах или Х1 больше ах то следующий исполняемый оператор *enddo)

CHG(N)=z

Вычисляется значение плотности заряда в узле как                                                                        функция координат (в данном случае заряд                                                                        однородный) и записывается в массив CHG в                                                  элемент с номером N.

bf,N,chrgd,CHG(N)     

В узле с номером N задается плотность объемного заряда (ключ chrgd) равная значению N-го элемента в массиве CHG.

*enddo

Конец цикла, переход к следующему узлу.

allsel,all

finish

/SOLUTION

SOLVE

FINISH

Переход к решению и решение электростатической задачи.

3. Вычисление распределения плотности электростатических сил. Третья программа называется Fpost1 и осуществляет постпроцессорную обработку результатов полученных предыдущей программой. Она вычисляет значения компонент плотности электростатических сил в узлах и записывает их в файлы. Для визуализации этих результатов можно воспользоваться методикой описанной в работе «Программа вычисления электростатических cил в Ansys».

/post1

Вход в постпроцессор.

*dim,FieldX,array,nodemx

*dim,FieldY,array,nodemx

Создаются массивы (FieldX и FieldY), в которые будут записываться значения X и Y компонент плотности электростатических сил.

*do,N,1,nodemx,1

Открывается цикл по всем узлам, N – переменная цикла.

*get,EX,node,N,ef,x              

*get,ey,node,N,ef,y

Для узла с номером N определяются значения X и Y  компонент электрического поля.

FieldX(N)=CHG(N)*ex

Fieldy(N)=CHG(N)*ey

Вычисляются компоненты плотности электростатических сил для данного узла и записываются как N-е элементы соответствующих массивов.

*enddo

Конец цикла

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

*cfopen,foX

Команда *cfopen открывает «командный» файл, создает файл foX и открывает его для записи данных. По умолчанию файл создается в рабочей папке, но при желании, можно указать путь к другой папке (см. описание команды *cfopen).

*vwrite,fieldx(1)

Данные из массива fieldx (начиная с первого элемента) записываются в файл foX.

(F12.6)

Формат записи чисел: F – число с плавающей точкой, 12 знаков отводится на число, 6 знаков после точки. Эта команда должна следовать непосредственно за  командой *vwrite

*cfclos

Закрывается «командный» файл
Скачать файл