Описание математической модели

Математическая модель - система уравнений, описывающих объект, и алгоритм (набор правил), определяющих последовательность решения уравнений модели и включающий набор значений параметров технологических производств.

Были приняты следующие допущения:

а) рассматриваемый объем гальванической ванны является реактором идеального смешения, т.к. обеспечивается непрерывное перемешивание электролита;

б) электролит на протяжении всего процесса анодирования алюминия не теряет своих свойств;

в) толщина и пористость покрытия постоянны на любом участке рабочей поверхности.

Существует три метода построения мат. моделей технологических объектах: экспериментальный, аналитический и комбинированный.

При экспериментальном методе построения формальных мат. моделей параметры определяются по опытным данным, полученным на действующем объекте.

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

Комбинированный метод заключается в объединении двух первых методов.

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

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

Классические теории проводимости АОП /1/ рассматривают зависимость ионного тока I_i, от напряженности электрического поля E:

 

, (1)

 

где А, В - постоянные.

Рост пленки происходит за счет переноса вещества ионами Me и O^2-.

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

 

, (2)

 

где a - ширина энергетического барьера в объеме оксида;

ν - частота колебания частиц в кристаллической решетке;

z - заряд катионов металлов;

W₀ - высота энергетического барьера;

δ - толщина оксида.

При малых E отсутствуют диффузионные затруднения в АОП, выполняется соотношение:

 

, (3)

 

где a - ширина энергетического барьера в объеме оксида;

ν - частота колебания частиц в кристаллической решетке;

z - заряд катионов металлов;

W₀ - высота энергетического барьера;

δ - толщина оксида.

Уравнения (3) и (4) достаточно хорошо согласуются с экспериментом, но при их выводе предполагается кристаллическое строение АОП. Однако, практически вовсех случаях сформированные AOП являются аморфными.

Эксперименты Дэвиса и Брауна /3/ c метками инертных газов показали, что новые слои АОП на алюминии в том числе, образуются как на внутренней, так и на внешней границах АОП. Было рассчитано число переноса металла (Al) t_M как отношение толщины АОП, образовавшейся на внешней поверхности, к общей толщине АОП. Для А1 по разным источникам составляет от 0,24 до 0,58 (для оксидирования А1 в этиленгликолевом растворе - 0,68).

Видно, что сейчас, при теоретической слабой базе, которая не до конца позволяет понять принцип процесса анодирования алюминия, аналитически невозможно вывести формулу нахождения коррозийной стойкости оксидного слоя. Поэтому создадим формулу, аппроксимировав табличные данные. Формулу представим в виде:

 

, (4)

 

где F - множитель полученный при аппроксимации данных взятых при изменении концентрации соли, плотности электролита и плотности тока;

T - множитель полученный при аппроксимации данных взятых при изменении во времени.

Данные для нахождения Т найдем из рисунка 6.1, где представлена зависимость коррозийной стойкости от времени протекания процесса, максимум находится в 1, потому что данные для аппроксимации функции F считаются при оптимальном значении времени анодирования.

 

Рисунок 6.1 - Зависимость коррозийной стойкости от времени протекания процесса.

 

Аппроксимировав данные получим:

 

, (5)

 

где t - время анодирования алюминия.

Проведем эксперимент, где будем менять плотность серной кислоты от 1,2 до 1,29 с шагом 0,03 и плотность тока от 2 до 8 А/дм² с шагом 1,5 А/дм², концентрацию соли будем изменять от 0,2 % до 2,6% с шагом 0,06. Результаты представлены в виде таблиц 6.1,6.2,6.3,6.4.

Таблица 6.1 - Зависимость коррозийной стойкости от концентрации соли и плотности тока при плотности электролита равной 1.2 кг/м³

	концентрация
Плотность Тока	6	7	7	7	7
	7	8	8	8	8
	8	9	9	9	8
	8	9	9	9	8
	7	8	8	8	7

 

Таблица 6.2 - Зависимость коррозийной стойкости от концентрации соли и плотности тока при плотности электролита равной 1.23 кг/м³

	концентрация
Плотность тока	7	8	8	8	7
	8	9	9	9	8
	9	10	10	10	9
	9	10	10	10	9
	8	9	9	9	8

 

Таблица 6.3 - Зависимость коррозийной стойкости от концентрации соли и плотности тока при плотности электролита равной 1.26 кг/м³

	концентрация
Плотность тока	7	8	8	8	7
	8	9	9	9	9
	9	10	10	10	9
	9	10	10	10	10
	8	9	9	9	8

 

Таблица 6.4 - Зависимость коррозийной стойкости от концентрации соли и плотности тока при плотности электролита равной 1.29 кг/м³

	концентрация
Плотность Тока	7	8	8	8	7
	8	9	9	9	9
	9	10	10	10	9
	9	10	10	10	9
	8	9	9	9	8

 

Необходимо определиться каким образом будем получать значения функции не в узлах таблицы.

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

Первой трудностью является то, что если представить многочлен в виде:

 

P (x,y) =a₀₀+a₁₀x+a₀₁y+a₂₀x²+ a₁₁xy+a₀₂y² +…+a_omy^m, (6)

 

где а_ij - коэффициенты.

То подставляя данные координаты точек и приравнивая левую часть соответствующему значению z_i, получим систему n+1 линейных алгебраических уравнений относительно 1+2+…+ (m+1) = (m+1) (m+2) /2 неизвестных коэффициентов a_ij. Вообще говоря эти уравнения независимы. Следовательно, если не накладывать на P (x,y) никаких дополнительных условий, то n+1 должно быть равно (m+1) (m+2) /2. Поэтому мы не можем решить поставленную задачу при произвольном количестве узлов интерполирования.

Второе принципиальное затруднение это то, что узлы интерполирования не могут располагаться произвольно. Рассмотрим на примере n=2 и n=5: для первого случая определитель будет обращен в нуль, если три точки (x₀,y₀), (x₁,y₁) и (x₂,y₂) лежат на одной прямой; во втором случае определитель будет обращаться в нуль если 6 точек интерполирования лежат на одной кривой второго порядка. Аналогично, если взять 10 узлов интерполирования, то определитель системы обратиться в нуль, если все они лежат на одной прямой третьего порядка. Проверка того, что определители не обращается в нуль, чрезвычайно затруднительна.

Третье принципиальное затруднение возникает при оценке остаточных членов. Так как теорема Ролля, для данного случая действовать не будет.

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

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

При аппроксимации главным допущением является то, что значения аргумента x₀,x₁,.,x_n найдены значительно точнее, чем значения функции f (x_i). Также будем предполагать, что систематические погрешности, а также грубые ошибки в значениях функции f (x_i) исключены (если нет то необходимо сгладить таблицные значения).

Для получения значения функции в точке x расположенных между x_i и x_i+1, будем применять сглаживающую функцию основанную на методе наименьших квадратов, предполагая, что x₀,x₁,.,x_n равностоящие, а все значения f (x_i) имеют одинаковую точность.

Пусть φ₀ (x), φ₁ (x),…, φ_m (x) - какая-то система линейно независимых функций на интервале [a,b], m≤n. Будем разыскивать обобщенный многочлен, составленный из этих функций:

 

, (7)

 

где f (x_i) - значение функции взятой из таблицы;

Ф (x_i) - значение подбираемой функции;

р_{i -} вес точности.

Далее необходимо, чтобы этот многочлен имел наименьшее значение. В нашем случае значения f (x_i) имеют одинаковую точность, поэтому:

 

, (8)

 

где р_{i -} вес точности.

Тогда получаем:

 

, (9)

 

где f (x_i) - значение функции взятой из таблицы;

Ф (x_i) - значение подбираемой функции.

Аппроксимировав данные взятые из таблиц 1,2,3,4 алгебраическими многочленами 1,x,x²,…,x^m, т.к. они образуют систему Чебышева на любом отрезке, а следовательно линейно не зависимы, получим зависимость коррозийной стойкости от плотности тока, концентрации соли и плотности серной кислоты. Для удобства читаемости запишем формулу в виде:

 

, (10)

 

где Ch - числитель в формуле,

Zn - знаменатель.

 

), (11)

 

где p - плотность заданного электролита,

c - концентрация соли в растворе электролита,

i - плотность тока.

 

, (12)

 

где p - плотность электролита.

Далее с учетом формул (4) и (5), получим:

 

, (13)

 

где p - плотность заданного электролита, c - концентрация соли в растворе электролита, i - плотность тока; t - время анодирования алюминия.

На основании этих допущений и зависимостей полученных экспериментальным путем на какой-либо конкретной установке и зная свойства электролита можно судить о любом похожем процессе.