обозначить через Z, подвод материи или энергии в систему |
через |
Qi (при этом отвод обозначается как отрицательный |
подвод |
—Qi), то для текущего состояния справедливо: |
|
Накопление (аккумулирование) можно выразить в виде про изведения какой-то величины U на соответствующую емкость С. При условии, что C = const, можно придать уравнению следующую форму:
В конкретных случаях Q, Z, С, U могут |
означать: |
1. Q — подвод (или отвод) |
жидкости в резервуар; |
|
Z — объем жидкости в резервуаре; |
|
|
U — уровень жидкости в резервуаре; |
|
|
С — поперечное сечение |
резервуара. |
|
2. |
Q — подвод (или отвод) |
тепла; |
|
|
Z — аккумулированное |
тепло; |
|
|
U — температуру; |
|
|
|
|
С — теплоемкость |
(произведение |
удельной теплоемкости |
|
на массу) ; |
|
|
|
3. |
Q — подводимая (или отводимая) мощность; |
|
Z — кинетическая |
энергия; |
|
U= — у м 2 (со — угловая |
скорость); |
|
|
С •— момент инерции массы (по отношению к оси вращения). |
4. |
Q — пусковой (или тормозной) момент; |
|
Z — момент количества |
движения; |
|
|
U — угловая скорость; |
|
|
С — момент инерции массы по отношению к оси вращения. Б. Составление уравнения материального баланса при протека
нии |
химических |
реакций. Рассмотрим |
случай, когда |
п |
веществ |
(Ai, |
Ai |
Ап-і, Ап) |
участвует |
в |
m |
независимых |
реакций. |
Если |
стехиометрический |
коэффициент |
і-го вещества в /-й |
реакции |
обозначить |
<Xij, знаком |
« + » |
обозначить |
те вещества, |
которые |
являются продуктом реакции, |
а знаком |
« —» те вещества, |
которые |
вступают в реакцию, то тогда |
для |
каждой реакции будет спра |
ведливо уравнение: |
|
|
|
|
|
|
|
|
п |
|
|
|
|
|
|
|
|
|
|
|
2 * / Л = 0 |
M = |
1. 2. |
m). |
|
|
|
|
|
(185) |
Теперь введем понятие концентрации d вещества Ai, определяе мое количеством молей этого вещества в единице объема или массы реагирующей смеси, и понятие диапазона (степени превра щения) при /-й реакции с\ который можно выразить через изме-
нения |
концентраций |
веществ, обусловленных /-й реакцией: |
Аса |
Лег,- |
|
с* = —У- = —У-. |
(186) |
Через диапазон отдельных реакций можно выразить изменение концентрации вещества Ai, обусловленное реакциями т;
m |
m |
(187) |
Д^=2Д^=-2а^> |
где Acjj взято из уравнения (186). |
|
Дифференцированием уравнения |
(186) по времени получим |
скорость ri /-й реакции, которая является функцией местного тер модинамического состояния
rid |
1 dhCii |
1 |
dL\cki |
Л |
088) |
r ; ^ = ^ |
^ |
— |
= |
••• С « . Я . П . |
Дифференцируя уравнение (187), вновь получим скорость изме нения концентрации вещества Ai, обусловленного всеми реак циями т:
d Дс/ |
= 2««УО- |
|
|
|
<1 8 9 ) |
dt |
|
|
|
|
7 = 1 |
|
|
|
|
|
|
Это основное выражение определяет изменения веществ, обу |
словленные |
химическими |
реакциями. Эту |
зависимость |
включают |
в уравнение |
материального |
баланса, приведенное в п. А. |
|
Аналогично |
рассматривают, например, |
непрерывные |
реакции |
с постоянным |
объемом |
реагирующих смесей V. В этом случае |
сумма |
подводимых объемов |
(отводимые объемы — это подводимые |
объемы с обратным знаком) |
равна 0: |
|
|
р
2 Q* = °-
Для каждого вещества Ai можно тогда написать:
Т |
Г Ш |
= 2с ; Л + (4гг)реакЦия=4~ |
2cikQk-\- 2a U r J - О 9 0 ) |
В. Составление уравнения теплового баланса при эндотерми ческих и экзотермических реакциях. Изменением давления при этом пренебрегают, что для большинства практических случаев вполне допустимо. При постоянном давлении энергетическое равновесие (баланс) реакции является равновесием энтальпий. Для этого равновесия (баланса) необходимо знать количество генерирован ного тепла (—àHj) или количество тепла, израсходованного ( + AHj) в j-й реакции при единичном изменении ее диапазона.
Это количество тепла можно выразить на основе парциальных молярных энтальпий НІ при помощи выражения
я
Ші=^іаиНі. |
|
|
|
(191) |
« = і |
|
|
|
|
Для скорости генерирования тепла всех реакций на единицу |
объема смеси справедливо |
выражение |
|
m |
m |
п |
А ) - |
(192) |
-^-=2 (-ЛЯ;)гі/ = 2rJ |
2 |
Число приведенных зависимостей, естественно, не является |
исчерпывающим, |
кроме того, существуют различные |
модификации |
и комбинации этих зависимостей. Параметры, приведенные в от дельных уравнениях, определяются или по таблицам характерных свойств веществ (плотность, удельное тепло, парциальная моляр ная энтальпия и др.), или на основании конструктивных данных
(особенно размеры, |
объемы |
и |
др.) либо |
результатов |
измерений. |
Довольно часто, |
особенно |
у |
линейных |
моделей, коэффициенты |
подсчитываются на |
основании |
результатов измерений |
процессов |
в исследуемой системе. Подставляя их в структурные уравнения модели, можно рассчитать коэффициенты, например, методом наименьших квадратов, известным из статистики.
Э к с п е р и м е н т а л ь н ы е с т а т и с т и ч е с к и е м о д е л и [15]. В тех случаях, когда связь между входными и выходными величи нами является слишком сложной или когда она определена недо статочно ясно, математическую модель можно составить на основе экспериментальных данных. При этом из-за ошибок измерения, сложности и неясности зависимостей либо большой математиче ской сложности всего решения достичь правильной структуры мо дели не удается и действительные зависимости приходится заме нять предполагаемой оптимальной аппроксимацией. В качестве критерия точности модели чаще всего применяют минимум суммы квадратов отклонений. В этом разделе будут рассмотрены стати стические модели такого вида.
Составление статистической модели многопараметровой си стемы включает следующие операции:
а) выбор и соответствующее выражение входных и выходных величин;
б) ограничение области, в которой будут производиться изме рения и которая одновременно будет служить областью, для кото рой модель справедлива;
в) выбор аппроксимирующей зависимости; г) разработку программы измерений, т. е. плана проведения
экспериментов с учетом необходимости охвата всей исследуемой области;
д) обработку результатов измерений и вывод уравнения модели; е) проверку и уточнение полученной модели.
Первым шагом является выбор основных величин, которые характеризуют изучаемый процесс, и соответствующее их выраже ние. Обычно начинают с регулируемых величин, для которых модель построена. Это — качественные величины, которые описы вают основные продукты процесса (материальные продукты или производимая энергия). Затем рассматривают регулирующие воз действия и возмущения.
Целесообразно выражать входные и выходные величины в от носительных значениях, так, например, вместо количества продук
ции |
используют выход годного, т. е. отношение количества продук |
ции |
к количеству использованного сырья, и т. д. |
|
Область измерений, как правило, ограничивается |
определением |
границ для всех входных, то есть независимых |
переменных |
величин (в статистике их называют факторами). |
В противном |
случае ограничения выходных величин следует выразить, пользуясь имеющимся опытом, через ограничения входных величин. Оче видно, что выделение экспериментальных областей, так же как и выбор величин, предполагает наличие фундаментальных теоре тических знаний о рассматриваемом процессе и определенного практического опыта по ведению такого или аналогичного процесса.
Выбор аппроксимирующей функции чаще всего осуществляется разложением неизвестной истинной зависимости в ряд Тэйлора вблизи середины экспериментальной области, причем члены выс ших порядков в расчет не принимаются (ими пренебрегают). Между каждой выходной величиной Л* и входными величинами
существует |
неизвестная связь |
|
* , = / ( Г і . |
У2, |
Уп). |
(193) |
Если средние значения независимых переменных величин обо |
значить как |
У?, Y2, |
. • ., У°п, то функция (193) |
после разложения |
в ряд Тэйлора и после пренебрежения членами, имеющими поря док выше 2-го, примет вид:
Xl = |
b0 + bl(Y-YÏ) |
+ |
b2{Y2-Y$)+ |
. . . +Ьп{Уп-У°п) |
+ |
+ ô n ( r I - r I ) 2 + ô 2 2 ( r 2 - r 2 ° ) 2 = |
. . . |
+bm{Yn-YlJ+ |
+ |
Ь„ (Y, - Y ?) ( Г 2 - |
Yï) + bl3 (Y, - |
F?) (Y3 - |
Y°3) + .... |
+ |
+ |
bn_hJYn^-Y0n^)(Yn-Y°n), |
|
|
|
(194) |
где b0, b\, b2, . . ., bij — неизвестные коэффициенты.
Если обозначить значения разностей маленькими буквами: Уі—К? = г/і и т. д., то можно написать:
Xi = b0+blyl-Jrb2y2-\- |
. . . |
-\-b„yn + 6 ц у ? + . . . +bnny2n |
+ |
Ч-*і2УіУ2 + |
*ізУіУз+ |
••• |
+*я - іУ„ - іУя - |
(195) |
Последнее |
уравнение |
является полиномом второй степени. |
Это уравнение является |
наиболее часто применяемой |
аппроксими |
рующей функцией. Полиномы первой степени дают неудовле творительную аппроксимацию. Использование полиномов высших
степеней требует проведения экспериментов и вычислений для опре деления коэффициентов. Например, для четырех входных величин (п—4) и для второй степени показательного полинома число чле нов полинома, а тем самым и число неизвестных коэффициентов b
равно 15. Поэтому, чтобы можно было определить |
коэффициенты |
b экспериментальным путем, необходимо провести |
минимум 15 |
измерений при различных комбинациях входных величин. Подста новкой результатов измерения в уравнение (195) при п = 4 получим систему 15 линейных уравнений, из которых можно вычислить
коэффициенты |
Ь. Обычно, |
однако, |
проводят |
большее |
количе |
ство измерений, |
чем это |
необходимо |
(в |
данном случае |
необхо |
димо 15), и определяют коэффициенты |
b |
методом наименьших |
квадратов. |
|
|
|
|
|
|
|
При выборе |
аппроксимирующей |
зависимости |
большое |
влияние |
на точность аппроксимации оказывает число измерений и способ комбинации входных величин при отдельных измерениях.
Если каждую из п входных величин привести в соответствие с одной координатной осью, то каждое измерение будет опреде ляться определенной комбинацией значений входных величин (независимых) и его можно представить в виде точки в п-мерном (эвклидовом) пространстве, которое называют экспериментальным (иногда векторным) пространством. Путем ограничения входных величин в этом пространстве выделяют область эксперимента и разрабатывают план измерений, т. е. предлагают определенное число и определенное расположение экспериментальных значений (точек) в области эксперимента.
Так, например, рекомендуется проанализировать корреляцию между всеми входными и выходными величинами. Этот анализ покажет основные, а также менее значимые (которыми можно пре небречь) зависимости между этими величинами. Оценочные значе
ния коэффициентов корреляции |
вычисляются по выражению |
2 |
(Уік — Уі) (хщ — Хі) |
(196) |
k = l |
|
где tfjk, Xik — значения /-й входной_ и_і-й выходной величины, определенные при k-м измерении, a tjj и — это средние значения /-й входной величины и і-й выходной величины, N — число измере ний (число экспериментальных точек).
Аналогично определяют также корреляцию между отдельными входными и выходными величинами, коэффициент корреляции ко торых должен был бы быть теоретически равен или близок к нулю. Если это не так, то следует проверить, чем обусловлена корреля ция между входными и выходными величинами — способом их определения или случайными факторами. Выявленные ошибки необходимо устранить.