Пиротехническая химия
Главная Начинающим пиротехникам Статьи Добавить статью Добавить материалы на сайт Поиск по сайту Карта книг Карта сайта
Книги в помощь
Военная история Изготовление и применение ВВ Пиротехника в военном деле Разное по пиротехнике Физика в пиротехнике Химия ВВ и составов
Новые книги
Суворов С. "Бронированная машина пехоты БМП -3 часть 1" (Военное дело)

Яковлев Г.П. "122 мм самоходная пушка образца 1944 г." (Военное дело)

Суарес Г. "Тактическое преимущество " (Военное дело)

Стодеревский И.Ю. "Автобиография записки офицера спецназа ГРУ " (Военное дело)

Семиколенков Н.П. "стрельба из танковых пулеметов " (Военное дело)
Применение ЭВМ для термодинамических расчетов металлургических процессов - Синярев Г.Б.
Синярев Г.Б., Ватолин Н.А., Трусов Б.Г., Моисеев Г.К. Применение ЭВМ для термодинамических расчетов металлургических процессов . Под редакцией Щепкин А.А. — М.: Наука, 1982. — 267 c.
Скачать (прямая ссылка): primenenevm1982.djvu
Предыдущая << 1 .. 28 29 30 31 32 33 < 34 > 35 36 37 38 39 40 .. 103 >> Следующая

где, пх, It2 — номера алфавитно-цифровых символов согласно списку, расположенному в массиве AZ. Полученное значение сравнивается с шаблонами, размещенными в массивах NR (простые переменные), NM (имена числовых последовательностей) и ND (дублирующие имена).
Символы, помещенные между знаком равенства и запятой или между двумя запятыми, обрабатываются как число и присваиваются той переменной из общего блока /DA/, которая была установлена сравнением предыдущего ключевого слова с шаблонами.
В подпрограмме производится синтаксический контроль обрабатываемого списка символов и при обнаружении ошибок печатается предупреждение. В этом случае оставшиеся перфокарты с исходными данными прогоняются через устройство ввода и производится возврат в вызывающую программу с присвоением переменной Kl значения 777.
CALAPP. Модуль реализует алгоритмы вычисления и аппроксимации термодинамических функций индивидуальных веществ.
Основная исходная информация к началу работы подпрограммы содержится в общем блоке /DA/. Массивы переменных, входящие в него, предназначены для размещения следующей информации: E(500) - энергия электронных состояний атомов или молекул (е / ) ; G (500) — статистические веса энергетических уровней (g,- ); OM(IOO) - частоты колебаний молекулы (о&); F(IOO) — приведенный изобарно-изотермический потенциал (Фр; T(IOO) - температура в узлах таблиц (7)); CP(IOO) -
удельная теплоемкость при постоянном давлении (CDf ); H(IOO) - энтальпия вещества [Hj (T) - Щ (О) ]; S (100) - энтропия (S^); R (30) - значения молекулярных и термохимических констант.
Ветвление алгоритма (выбор пути расчета) производится на основе анализа введенных исходных данных. Те из переменных, значения которых не заданы в текущем варианте расчета, к моменту начала работы модуля оказываются равными 0,77777. Это позволяет в одной подпрограмме и при единообразном задании исходных данных производить вычисления и аппроксимацию свойств по разным алгоритмам. В подпрограмму заложены блоки вычисления термодинамических функций газообразных атомов и атомарных ионов; определения свойств двухатомных газообразных молекул методом Майера и Гепперт-Майер; вычисления термодинамических функций газообразных молекул в приближении "жесткий ротатор—гармонический осциллятор". 79
N= [...[(W1 -1)60+ (и2-1)] 60 + ...,
Аппроксимация выполняется методом наименьших квадратор для пяти полиномов разной сложности. Из них выбирается самый простой, удовлетворяющий заданной точности. При невозможности гарантировать требуемую погрешность аппроксимации выполняется разбиение температурного интервала и по шедовательное нахождение полиномов щ каждом из участков. Ошибки, выявленные в результате обработки исходных данных или в процессе счета, выводятся на печать. /
Из подпрограммы производится вызов модулей N^MASS и POLIN.
SYSTEM. Эта подпрограмма осуществляет рещ^ние системы линейных уравнений, коэффициенты которой готовятся в модуле POLIN. Решение производится с двойной точностью общим методом Гаусса. Предварительно выполняется перестановка строк и столбцов матрицы с тем, чтобы ее главные элементы имели наибольшее значение.
Модуль вызывается из подпрограммы POLIN. При этом порядок системы уравнений N задается как формальный параметр, а матрица коэффициентов А предварительно-заносится в неименованный общий блок. Результаты решения располагаются в массиве TR общего блока /REZ/.
FLAGS. Заданная в алфавитно-цифровом виде химическая формула рассматриваемого индивидуального вещества в этой подпрограмме преобразуется в условное число — две целочисленные константы JA и JB. Исходная строка символов, имеющих формат Al, располагается в массиве CS(IS).
Для преобразования используются алфавитные константы, а также коды символов химических элементов, расположенные в общем Ьлоке /SYMB/.
Условное число, так же как и правила загчеи исходной химической формулы, остались теми же, что и для программы определения параметров равновесного состояния. При обнаружении ошибок в записи или кодировке исходной информации программа присваивает переменной JA значение 777 и возвращается в вызывающий модуль READDA, где организуется печать предупредительного сообщения, а также выполняются действия, приводящие к переходу на новый вариант исходных данных или к завершению задания.
POLIN. Модуль с этим именем реализует алгоритм нахождения полинома, аппроксимирующего методом наименьших квадратов исходные термодинамические данные, заданные в виде таблиц.
Значения энтальпии, энтропии, удельной теплоемкости, изобарно-изотермического потенциала и температуры размещаются в общем блоке /DA/ в массивах Н, SS, CP, FhT. Аппроксимации подлежат термодинамические функции, находящиеся в узлах таблиц и номерами от Nl до N2. Те переменные, значения которых не определены, предварительно заполняются константами 0,77777, что позволяет игнорировать их в процессе обработки.
В подпрограмме производится нахождение коэффициентов для пяти полиномов, а также максимальных погрешностей по каждому из них, вычисленных относительно аппроксимируемых функций. Коэффициенты помещаются в массив R (5,9), а наибольшие отклонения и соответствующие им температуры - в массив С (5,8). Решение системы линейных уравнений, требуемое по методу наименьншх квадратов, производится с помощью подпрограммы SYSTEM.
Предыдущая << 1 .. 28 29 30 31 32 33 < 34 > 35 36 37 38 39 40 .. 103 >> Следующая
Реклама
 
 
Авторские права © 2010 PiroChem. Все права защищены.