CC BY
60
УДК 621.316.9:620.193.7
РАСЧЕТ КОРРОЗИОННЫХ ТОКОВ ДЛЯ МОДЕЛИ ЗАЗЕМЛЯЮЩЕЙ СИСТЕМЫ
Л.Д. Сафрошкина, Ю.Р. Гунгер, В.Е. Дмитриев, Ю.В. Демин
Новосибирская государственная академия водного транспорта. г. Новосибирск E-mail: [email protected]
Рассматриваются расчетные модели коррозионных токов и потенциалов (напряжений), возникающих при грунтовой (электрохимической) коррозии подземного сооружения, состоящего из металлически связанных между собой горизонтальных и вертикальных линейных элементов с учетом сопротивлений поляризации элементов, их сопротивления растеканию и взаимных сопротивлений.
В процессе эксплуатации в агрессивных условиях заземляющие системы (ЗС) электроустановок подвергаются воздействию электрохимической коррозии и разрушаются [1]. Для разработки методов защиты от коррозии ЗС и оптимизации применения защит необходимо знать распределение коррозионных токов в указанных системах. С этой целью предлагается метод расчета коррозионных токов в ЗС, физической основой которого является теория многоэлектродных электрохимических систем. В основе математической модели лежит система нелинейных уравнений, связывающих значения электродных электрохимических потенциалов и коррозионных токов.
Для этого типа коррозии рассматривается эквипотенциальный вариант решаемой задачи. При этом пренебрегают продольным сопротивлением элементов для протекающих по ним токам по сравнению с их (элементов) поперечным сопротивлением для стекающих в грунт токов. В этом случае все элементы рассматриваемой системы соединены "звездой" и имеют единый компромиссный потенциал.
В рассматриваемом решении поляризационные характеристики элементов задаются в виде электродных потенциалов (напряжений) иэ и сопротивлений поляризации ЙА (анодное) и Й (катодное). При более строгом решении учитывается зависимость сопротивления поляризации каждого элемента от величины коррозионного тока. Последнее достигается заданием по каждому элементу исходной информации о значениях анодного и катодного напряжений, отвечающих фиксированным величинам коррозионных токов.
Решение рассматриваемой задачи осложняется отсутствием исходной информации о характере поляризации (анодная или катодная), соединенных в систему элементов. Существующие в настоящее время решения по определению типа поляризации элементов не учитывают их взаимных гальванических связей [2]. Использование же итерационных решений без задания исходной информации о характере поляризации элементов зачастую ведет к расхождению итераций. В предлагаемом нами решении исходный тип поляризации для каждого элемента определяется путем предварительного анализа решаемой системы уравнений с учетом гальванических связей всех элементов. И только в
случае несовпадения знаков коррозионного тока, полученных по предварительному анализу и при окончательном решении для отдельных элементов, для них производится смена типа поляризации с последующим запуском итерационного процесса. При этом, как правило, случаи использования итераций очень ограничены и обычно связаны с некорректной постановкой задачи.
В случае необходимости уточнения полученного решения с учетом нелинейной зависимости сопротивления поляризации элементов от величины коррозионных токов, это выполняется с помощью итерационных циклов, использующих вышеупомянутые зависимости напряжений поляризации (анодных и катодных) от величины коррозионных токов. Но при этом, как правило, уточняются только величина сопротивления поляризации, а их тип не меняется.
Таблица 1. Поляризационные характеристики для контрольной расчетной модели
Номера элементов модели Сопротивление поляризации, Ом Электродный потенциал Ц, В
Ансдное RA Катодное RK
1 82,22 36,30 -0,12
2 183,48 27,12 -0,10
3 400,27 36,91 -0,24
4 221,87 21,87 -0,28
5 57,14 60,00 -0,35
6 51,51 40,40 -0,21
7 18,69 20,56 -0,33
8 33,73 19,68 -0,28
9 61,92 16,54 -0,30
10 6,30 6,30 -0,53
11 40,00 53,53 -0,51
12 40,40 3737 -0,47
13 14,04 137,37 -0,37
Рассматриваемая расчетная модель предполагает однородное строение грунта. В основу этого положены следующие факторы:
- существенное (как правило, в несколько раз) преобладание собственных значений поперечных сопротивлений над взаимными;
- в отличие от обычного заземлителя (потенциалы на котором задаются относительно бесконечности коррозионные напряжения относятся только к ближней области [1, 2], поскольку свя-
заны с возникновением на поверхности металла двойного электрического слоя;
- коррозионные токи замыкаются между элементами внутри сооружения, а не вытекают (втекают) из него, как в случае заземлителя.
На основании вышеизложенного, при расчетах собственных значений поперечных сопротивлений для каждого элемента предусмотрена возможность использования своего удельного сопротивления грунта. При расчетах же взаимных поперечных сопротивлений задается общее для всего сооружения удельное сопротивление грунта. Значения сопротивления поляризации приведены в табл. 1.
Исходные системы уравнений и способы их решения для эквипотенциального варианта грунтовой коррозии
Исходная система уравнений для рассматриваемого вида коррозии может быть записана в виде следующих матриц:
Дь— Щ р ,— Щп,- 1 11 и,
Лр1,- щрр ,— Щрп ,— 1 X I, = Цр
Д,ь- Япр ,— Япп, — 1 Тп и„
1 1 1 0 и0 0
(1)
Щ11,— , Щ р ,— , Щ,п I, 1 - 0
Rp1,- , ЩРР,- , Щрп X Тр = Р 1 : 0
Щп1,- , Щпр, - — , Щпп Тп п 1 : 0
ладающей по отношению к сопротивлению Щ, что является решающим при отсутствии изоляционного покрытия элементов (Л“).
Для ориентировочного определения характера поляризации рассматриваемых элементов проведем предварительный анализ системы (1), записав их в следующем виде:
(2)
Форма записи системы (2) предполагает, что значение компромиссного потенциала и известно. При этом, поскольку собственные значения Щ поперечных сопротивлений существенно (как правило, в несколько раз) превосходя-т их взаимные значения Щ, знак поперечного тока I будет определяться зна--ом правой части Щ-Щ ур. (2). Общая сумма токов Р в ур. (2) должна равняться нулю, согласно последнему уравнению системы (1).
Если исключить из системы (2) р-ое уравнение путем его вычитания из остальных уравнений, а затем убрать из квадрат-ой матрицы столбец коэффициентов для тока I, то получим дополнительную (остаточную) систему уравнений без р-го элемента и без компромиссного потенциала и0;
где Щ - собственные и взаимные значения попе-речны—сопротивлений между /-ым и р-ым элементами; I - поперечный ток, отекающий (втекающей) с р-го элемента; и0 - компромиссный потенциал системы рассматриваемых элементов, соединенных "звездой"; и - исходные (задаваемые) электродные потенциалы одиночных элементов (до их соединения в систему).
В состав собственного (/=р) значения поперечного сопротивления Кц=Щ+Щ+Щк входят: Щ -сопротивление растеканию металл-грунт; Щ - поперечное сопротивление изоляционного покрытия элемента (при его наличии); Щк - сопротивление поляризации (анодное или катодное) рассматриваемого элемента.
Последнее уравнение системы (1) фиксирует равенство нулю суммы стекающих с элементов (катодных) и втекающих в них (анодных) токов. Последняя отрока квадратной матрицы, также представленный единицами (в конце - нуль), отвечает компромиссному -потенциалу Ц0, который как и поперечные токи I требует определения.
Основная трудность решения рассматриваемой системы (1) заключается в определении типа (анодное или катодное) сопротивления поляризации соединенных в систему э-лементов. Действительно, знак поперечного тока I (плюс - катодный, минус - анодный) определяет и величину сопротивления поляризации, поскольку значения Щ и Щ могут существенно отличаться друг от друга. К тому же, величина Щк является, как правило, преоб-
(Д,-
(Л
(Л
(Л.1-
■ад,—, - ад,- ЗД,-
■зд,—
(Д( р+1
(Л Р-(Л Р+
Р-
-лр(рМ)),(Д(р+1) -
- лр( р-1)),(Д р-1)( р+1
- Щр( р-1)),(Д р+1)( р-1
-Rp(p-1)),{■Щ(Р+1) -
Р+1)),-
) - ЛрР р+1) ) - ЛрРр+1 ЛР( Р+1)),-
(Дп-
(Щ Р-1)п (Щ Р + 1)п (Щп -
Лрп )
- Лрп )
- Лрп ) Лрп )
I, и -и,
Т( Р-1) и(Р- )-ир
I( Р+1) и(Р+ )-ир
Iп ип -и,
(3)
Тип сопротивления поляризации Щк в диагональных коэффициентах (Д-Щ) системы (3) будем задавать в соответствии со знаком правой части уравнений, т.е. при (Ц-Ц,>0) - используется катодное сопротивление Щ, при (Ц-Ц,<0) берется анодное сопротивление Щ, а при (Ц-^=0) -рекомендуется использовать среднее значение сопротивлений поляризации (ЛА+ЛК)/2.
Вспомогательная система (3) используется для ориент-ировочного определения суммарного значения (Щ) коррозионных токов без учета р-ого эле-мента.Р Для этого достаточно обратить квадратную матрицу (3) и провести суммирование обращенных коэффициентов по столбцам. Умножив "постолбцовые" суммы на соответствующие значения правых частей системы уравнений (3) и сложив полученные результаты, -найдем искомую сумму токов. Очевидно, что при Е/>0 - поляризация рассматриваемого р-ого элемента анодная (т.к. для баланса коррозион-ных токов- не хватает отрицательного
значения I), а при и,<0 - поляризация р-ого эле-1 Мр -1
мента катодная.
Вышеописанному анализу (с целью предварительного определения типа поляризации) подверга-
ются все элементы рассматриваемого сооружения. Исключение делается только для элементов с минимальным и максимальным значениями (с учетом знака) электродного потенциала, поскольку первый из них всегда будет анодом, а второй - катодом.
Определив ориентировочный характер поляризации всех элементов сооружения и, задействовав в соответствии с этим требуемые сопротивления поляризации Щ,к, решим исходную систему ур. (1) относительно коррозионных токов и компромиссного потенциала. Сопоставляя знаки полученных коррозионных токов с результатами предварительного определения типа поляризации, убедимся в их совпадении. В случае несовпадения знаков по отдельным элементам их тип поляризации меняется в соответствии с результатом последнего расчета, и этот расчет повторяется до полного совпадения знаков коррозионных токов по данным последнего и предыдущего расчетов, т.е. запускается итерационный процесс.
Аналогичные итерационные циклы используются и при уточнении (в случае необходимости) полученных коррозионных токов за счет учета нелинейной зависимости сопротивления по-ляризации (ЩА,к) от величины коррозионного тока (I). Но при этом, как правило, уточняется только величина сопротивления поляризации, а ее тип не меняется.
В табл. 2 приведены результаты контрольных расчетов коррозионных токов.
Расчетное значение компромиссного потенциала составило ио=-0,335 В. Отметим также, что знаки коррозионных токов по предварительному определению и при окончательных расчетах совпали и дополнительных итераций не потребовалось.
Таблица 2. Значения коррозионных токов для случая грунтовой коррозии
Номер элементов мА Режим работы элементов
1 5,10 катод
2 7,02 катод
3 2,30 катод
4 2,45 катсд
5 +0,01 катод
6 3,06 катод
7 0,15 катод
8 2,10 катод
9 1,46 катод
10 -16,23 анод
11 -3,96 анод
12 -2,84 анод
13 -0,61 анод
Как видим из приведенных результатов расчетов, компромиссный потенциал системы ио оказывается достаточно близким к электродным потенциалам элементов №№ 5 и 7 (табл. 1), чем и объясняется малая величина их коррозионных токов. Вместе с -ем, поскольку и5< ио и и7>—0, коррозионный ток 15 оказывается анодным, а 17 - катодным. При этом, хотя модули разности потенциалов
|и5-ио|=|и7-ио| равны, коррозионный ток |/7| практически в 1о раз превосходит по величине ток 151 что определяется существенно большим (почти в 3 раза) значением сопротивления поляризации Щ по отношению к сопротивлению ЩА (табл. 1).
5 - номера горизонтальных элементов (Ьэ = 25 м);
13 - номер вертикального элемента (Ьэ = 5 м).
Рис. 1. Расчетная модель заземлителя
Аналогичная корреляция между знаком и величиной коррозионного тока, с одной стороны, и разностью потенциалов (Ц-Цо), а также сопротивлением поляризации ЩАк, с другой, прослеживается и по остальным элементам рассматриваемого сооружения. Эта корреляция еще более усиливается при учете сопротивления растеканию металл-грунт Щ, которое для горизонтальных элементов составило 6,11 Ом, а для вертикального - 17,6 Ом (сопротивление изоляции Щи =о). Взаимные сопротивления растеканию Щ имеют существенно меньшую величину и даже для соприкасающихся элементов составили для:
- параллельных горизонтальных о,87 Ом,
- ортогональных горизонтальных 1,1о Ом,
- горизонтального и вертикального 1,93 Ом.
Для более удаленных элементов величина взаимного сопротивления Щ составляет десятые доли Ом. Таким образом, влияние сопротивления растеканию на величину коррозионных токов проявляется наиболее сильно на элементах № 1о (горизонтальный) и № 13 (вертикальный), где они оказываются сопоставимыми с действующими сопротивлениями поляризации (ЩАо=6,3о и Щ/3=14,о Ом).
Малая величина поперечного сопротивления Щ1о=Щи)+Щ1о=6,3о+6,11=12,41 Ом и большая разность электродного и компромиссного потенциалов Ц1о-Цо=-о,525+о,335=-о,19о В определила максимальную величину коррозионного анодного тока (-16,2 мА) у элемента № 1о. С другой стороны, несмотря на сравнительно большую разность потенциалов Ц2-Цо=-о,1о4+о,335=о,231 В у элемента № 2, относительно повышенная величина его поперечного сопротивления Щ2=ЩК+Щ2=27,12+6,11=33,23 Ом привела к с-нижению стекающего с него коррозионного тока до /2=7,о238 мА (максимальный катодный ток).
--------суммарные кривые по двум электродам;
...... — кривые поляризации Си<А), Си<к), Ее<А), Ее<к), Ре<А), Ре<к);
--------сопротивления растеканию элемент-грунт Щ.
Рис. 2. Поляризационные кривые экспериментальных измерений для трех электродов — Ре (сплав), Ее, Си <диаграмма Эванса)
Определение значений сопротивления поляризации. Будем определять сопротивление поляризации (В**) как отношение подаваемого на поляризуемый элемент напряжения Щ’к (за вычетом величины электродного потенциала Щк) и создав—.емо-му этим напряжением коррозионному току (7эАК):
яА,К = (иА,К —и АоК)/ 1эА,К
(4)
С этой целью измерим с помощью электрода сравнения электродный потенциал ра—сматривае-мого одиночного элемента (при этом 7эАД=0), а затем будем увеличивать или уменьшать (алгебраически) этот потенциал путем подачи на элемент дополнительных напряжений от внешнего источника. При этом будем одновременно измерять и значение искусственно создаваемого потенциала (с помощью электрода сравнения) и величину возникающего в результате поляризации элемента коррозионного тока. Определение тока выполняется с помощью подключения между элементом и землей измерительной цепи, содержащей стандартизированное сопротивление и милливольтметр.
Проводя вышеописанные измерения при алгебраическом уменьшении (Щ) искусственно создаваемых на элементе потенциалов (по отношению к электродному) или их увеличении (иА), получим со-
ответственно анодные или катодные поляризационные кривые, определяющие значения анодного (ВА) или катодного (В/) сопротивления поляризации для фиксированной величины анодного (7/=^^^ или катодного (7/=^^^ коррозионного тока (рис. 2).
Сопротивление поляризации нелинейно зависит от величины коррозионного тока. При этом, как правило, характер анодных и катодных поляризационных кривых (а, следовательно, и величины отвечающих им сопротивлений поляризации) оказывается существенно различным. Более того, электродные потенциалы (на—чальные значения поляризационных кривых при 7эАК=0) при анодной и катодной поляризации элемента могут иметь разные значения. Происходящие при этом изменения в правой части исходной системы уравнений (1) могут приводить к существенному изменению уровня компромиссного потенциала, а, следовательно, и характера поляризации ряда элементов, что резко усложняет ход итерационного процесса.
В этом плане предлагаемое определение сопротивления поляризации способствует более устойчивому ходу итерационного процесса, поскольку уточняется только один параметр — сопротивление поляризации, а значения электродных потенциалов остаются неизменными.
Из расчетов следует, что для снижения коррозии и повышения долговечности заземляющей системы необходимо повышать сопротивление растеканию горизонтальных шин, так как они объединяют всю систему в единое целое, а по условиям их работы в качестве искусственных заземлителей его наоборот надо понижать. Это противоречие можно разрешить поместив горизонтальные шины в засыпку, имеющую нелинейную электропроводность. Это в нормальном режиме работы означает повышение сопротивления растеканию, а в переходном резкое его уменьшение.
На основании расчетов коррозионных токов для модели заземляющей системы можно сделать следующие выводы:
СПИСОК ЛИТЕРАТУРЫ
1. Демин Ю.В., Демина Р.Ю., Горелов П.В. Обеспечение долговечности электросетевых материалов и конструкций в агрессивных средах. Книга 1. Теоретические основы / Под ред. д.т.н., проф. В.П. Горелова. - Новосибирск: НГАВТ, 1998. -209 с.
1. На коррозию искусственных заземлителей максимальное влияние оказывают сопротивления поляризации и сопротивления растеканию. Взаимные сопротивления оказывают на коррозионные токи минимальное влияние и ими можно пренебречь.
2. Для обеспечения растекания токов коротких замыканий (снижения сопротивления растеканию) и снижения токов коррозии (повышения сопротивления растеканию) необходимо помещать горизонтальные шины в засыпку, имеющую нелинейную электрическую проводимость.
2. Иоссель Ю.Я., Кленов ГЭ. Математические методы расчета электрохимической коррозии и защиты металлов. Справочник. -М.: Металлургия, 1984. -271 с.
УДК 621.374
СИНТЕЗ АВТОКОЛЕБАТЕЛЬНЫХ СИСТЕМ ГЕНЕРАТОРОВ ЭЛЕКТРИЧЕСКИХ СИГНАЛОВ
Ю.К. Рыбин
Томский политехнический университет E-mail: [email protected]
Периодические сигналы широко используются в измерительной технике в качестве тестовых сигналов, в медицине в качестве стимулов, в строительной индустрии при акустических испытаниях строительных конструкций, в промышленности при испытаниях механизмов и машин и т.д. Известны и широко применяются приборы и устройства для их воспроизведения — генераторы периодических сигналов. Развиты и методы построения автоколебательных систем генераторов для воспроизведения периодических колебаний простейших форм. Тем не менее, задача синтеза математических моделей автоколебательных систем до конца ещё не решена. В статье рассматривается метод стационарных автоколебаний, на основе которого несложно построить систему уравнений, решением которой являются колебания заданной формы.
Введение
Периодические электрические сигналы с предписанной (заданной) формой используются в измерительной технике для питания различных датчиков физических величин, измерения передаточных характеристик четырёхполюсников. Они применяются в метрологической практике при поверке средств измерений параметров переменных напряжений, в медицине для кардиостимуляции, в автоматике для создания программных движений, в строительстве при виброиспытаниях различных конструкций и т.д. Поэтому проблема получения сигналов с заданной формой или с заданным спектром продолжает оставаться актуальной. Для решения этой проблемы развиваются два способа: аналоговый и цифровой. Аналоговый способ основан на формировании сигналов непрерывной динамической системой, а цифровой базируется на воспроизведении сигналов дискретной динамической системой. К достоинствам
первого способа относят практически неограниченный диапазон частот, а к преимуществам второго -простоту воспроизведения сигналов сложных форм. Проблема создания динамической системы с заданным поведением имеет не только прикладное, но и самостоятельное теоретическое значение в теории колебаний. Для решения указанной проблемы можно использовать метод локализации [1], метод управления по производной [2] и др. В работе [3] автором предложен метод синтеза автоколебательных систем с заданной синусоидальной формой колебаний в стационарном режиме. Этот метод назван методом стационарных автоколебаний. Далее этот метод применяется для синтеза систем, воспроизводящих периодические колебания других форм.
Постановка задачи исследований
Метод стационарных автоколебаний основан на том, что в динамической системе второго порядка
