Программа «Алмаз» по расчету ЭОС методом анализа

2.2. Программа «Алмаз» по расчету ЭОС методом анализа.

Для расчета ЭОС методом Анализа изложенном в параграфе 1.3.2 использована программа «Алмаз». Эта программа состоит из двух загрузочных модулей: aupr.exe – расчетный модуль, grafl.exe – графический модуль.

Для выполнения расчетов по программе aupr.exe необходимо предварительно подготовить файл исходных данный «fd». Затем выполнить расчеты с помощью программы aupr.exe. При этом по запросу ЭВМ указать файл вывода результатов расчета «frl». В процессе расчета программа сама создает следующие файлы для построения результатов расчета в графической форме:

geom - для построения геометрии,

traek - для построения траекторий электронов,

tok - для построения распределения плотности тока.

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

При этом в программе grafl.exe работают пункты меню в соответствии со следующей таблицей:

1.   Геометрия – работает.

2.   Траектории – работает.

3.   Эквипотенциали эл. поля – не работает.

4.   График распределения BZ – не работает.

5.   График распределения JZ – не работает.

6.   График токооседания – не работает.

При повторном выполнении расчетов необходимо предварительно ввести в файл «fd» новые данные и удалить файлы «frl», «geom», «traek», «tok», старого варианта.

При подготовке исходных данных ввод исходных данных осуществляется с дисплея. Данные вводятся в виде последовательности строк, содержащих наименования массивов и переменных и их числовые значения.

Числовые значения вводятся в виде строки чисел действительного типа форматом 8F8.0, причем последовательность ввода массивов и переменных должна строго выполнятся.

При подготовке исходных данных необходимо придерживаться следующих правил:

1. Сначала вводится строка, содержащая наименования переменных или массива и символы – разграничители «I», определяющие начало и конец поля, в котором располагаются числовые значения в соответствии с указанным форматом (8F8.0).

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

3. Числа могут располагаться в любом месте в пределах отведенного для них поля.

4. Количество символов в числе, включая знак и кодовую точку для чисел действительного типа не должно превышать ширины поля.

5. Знак «+» у чисел можно опускать.

Далее приведены наименования переменных, их назначение, последовательность массивов, которые объединяют эти переменные:

Массив 1 описывает общие данные о приборе, данные определяющие точность решения:

RU – максимальные размер области для поля по оси R.

RF – максимальный размер области для пространственного заряда по оси R.

ZU– максимальные размер области для поля по оси Z.

TTT – конец прибора по оси Z и последняя плоскость симметрии.

FH – число узлов, приходящихся на область перекрытия. FH = S /H, где S – протяженность участка прибора, общего для соседних областей (область перекрытия) S = 1.5dk, где dk – диаметр канала.

H – шаг разностной сетки, выбираемой из условий:

H>=RU/147,

H>=RF/147,

H>=ZU/297.

VQ – шаг интегрирования. VQ = (2-3)*H.

U – анодное напряжение в вольтах.

FK – расстояние от катода до плоскости. В которой определяется ток и первеанс прибора. FK = (2-3)*H.

RK – радиус кривизны катода.

HK – высота катода.

ZO – координата по оси Z центра окружности катода.

Y1, Z1 Y2, Z2 – координаты конического катода (временно не используются) и полагаются равными Y1 = Z1 = Y2 = Z2 0.

FE – число слоев электронного потока FE<=30.

GE – критерий сходимости потенциала при расчете поля. GE = 0.001 – 0.0001.

RM– критерий неоднородности магнитного поля.

NP – число плоскостей симметрии, которыми прибора разбивается на отдельные области NP<=30.

IWN – число линий второго порядка, описывающих контур прибора.

IWP – число линий первого порядка, описывающих контур прибора.

NPR – максимальное число последовательных приближений при расчете каждой области.

NS – параметр, задающий частоту вывода цифровой информации на печать. При NS = 0 информация выводится на каждом шаге интегрирования.

NPL – число вводимых в массив Х15 плоскостей.

TK – температура катода в градусах Кельвина. Для расчета без учета влияния тепловых скоростей положить TK = 0.

NEG – число энергетических групп для многоскоростного потока NEG = 1-3.

Массив 2 описывает плоскости, в которых определяется токопрохождение.

Х15 (1) – Z1

Х15 (2) – Z2

Х15 (30) – Z30

Примечание:

1. Первые плоскости должны находиться в местах геометрической симметрии прибора. Число их не должно превышать 30.

2. Х15 (1) равно начальной границе прибора.

3. Х15 (NP) = TTT.

4. При расстановке плоскостей симметрии, на которые разбивается прибор, должно выполняться условие:

max (Х15 (i + 1) - Х15 (i)) <= ZU – FH*H

Массив 3 описывает контур прибора и потенциалы на электродах. Геометрические размеры прибора должны быть заданы в (мм).

Примечание:

1. Все числа, описывающие данную линию набираются в одной строке.

2. Все линии второго порядка должны быть описаны в начале массива восемью числами (Z_H, R_H, Z_K, R_K, U_H, Z_O, R_O, R). Линии первого порядка описываются после линий второго порядка пятью числами (Z_H, R_H, Z_K, R_K, U_H).