1882
.pdf
|
|
Окончание табл. 2 . 3 3 |
|
1 |
2 |
3 |
4 |
Значения lnA |
|
|
|
С=84,5 % |
22,64089284 |
21,2204301 |
17,8329533 |
С=65 % |
22,9789164 |
19,7204139 |
17,6477951 |
С=39 % |
24,1673998 |
22,4104682 |
19,5893143 |
Значения А |
|
1,6441E+09 |
|
С=84,5 % |
6,8048E+09 |
5,5559E+07 |
|
С=65 % |
9,5415E+09 |
3,6683E+08 |
4,6168E+07 |
С=39 % |
3,1316E+10 |
5,4043E+09 |
3,2176E+08 |
Графические зависимости lnA(Т,С) и b(Т,С), построенные по данным таблицы приведены на рис. 2.35,а и 2.35,б.
а
б
Рис. 2.35. а – график lnA(Т,С); б – график b(Т,С)
Для аппроксимации поверхности lnA(Т,С) использовалась функция вида:
lnA(Т,С) = x1 C +x2 T + x3 C2 +x4 T2+x5 TC + x6,
для определения коэффициентов x1,x2,x3,x4,x5,x6 которой применялся метод наименьших квадратов. Результаты расчета приведены в табл. 2.34.
191
Для аппроксимации поверхности b(Т,С) использовалась функция вида: b(Т,С) = у1 C +у2 T + у3 C2 +у4 T2+у5 TC + у6,
значения коэффициентов которой, определенные с помощью метода наименьших квадратов, приведены в табл. 2.35.
Таблица 2 . 3 4 Значения коэффициентов зависимости lnA(Т,С)
x1 |
-28,35286386 |
x2 |
0,025724277 |
x3 |
21,46029289 |
x4 |
-0,000768897 |
x5 |
-0,015860578 |
x6 |
34,52385881 |
Таблица 2 . 3 5 Значения коэффициентов
зависимости b(Т,С)
y1 |
-2,444115069 |
y2 |
0,065908239 |
y3 |
1,731809485 |
y4 |
-0,000409049 |
y5 |
-3,69119E-05 |
y6 |
1,089807767 |
С использованием модели коррозионного растрескивания (2.98) при найденных зависимостях коэффициентов этой модели от температуры и концентрации агрессивной среды проведено исследование влияния программы изменения напряжения, температуры и концентрации среды на кинетику коррозионного растрескивания арматуры. Предварительно, с помощью численного эксперимента была определена величина шага интегрирования по времени, необходимая для достижения требуемой точности.
Шаг по времени, час |
Погрешность при определении |
|
долговечности, % |
||
|
||
0,1 |
3,61 % |
|
0,01 |
0,48 % |
|
0,001 |
0,06 % |
Дальнейшие расчеты выполнялись с шагом по времени 0,01 ч. Результаты исследования влияния различных программ изменения
напряжения, температуры и концентрации среды на кинетику коррозионного растрескивания арматуры приведены на рис. 2.36.
192
Рис. 2.36а. Кинетика коррозионного растрескивания преднапряженной арматуры (модель Работнова) при различных программах внешних воздействий:
1 – =865 МПа = const; С = 39 % = const; T = 110 C = const; 2 – =865 МПа = const; С = 39 % = const; T =110 -5* t, C; 3 – =865 МПа = const; С = 39 + 10*t, %; T = 110 C = const
Рис. 2.36б. Кинетика коррозионного растрескивания преднапряженной арматуры (модель Работнова) при различных программах внешних воздействий:
1 – =865 МПа = const; С = 39 % = const; T = 110 C = const; 2 – = 865 – 50*t, МПа; С = 39 % = const; T = 110 C = const; 3 – = 300 +50*t, МПа; С = 39 % = const; T = 110 C = const
193
2.6.3.Моделирование деформирования
икоррозионного растрескивания арматурного стержня
Рассмотрим задачу построения модели деформирования и разрушения арматурного стержня, рассматриваемого в общем случае как объемный цилиндрический стержень, находящийся под действием механических нагрузок в условиях коррозионного растрескивания, приводящего к изменению механических свойств материала. Так как коррозионное растрескивание происходит при практически полном отсутствии пластической деформации, то полагаем материал арматурного стержня нелинейно-упру- гим, а его свойства принимаются зависящими от концентрации водорода в точках стержня. Считаем, что и нагрузка, и водородсодержащая среда действуют на цилиндрический арматурный стержень осесимметрично, причем и закон распределения нагрузки и характер водородного воздействия не меняются по длине арматурного стержня.
Применим для моделирования поведения арматурного стержня в условиях совместного действия нагрузки и коррозионной среды теорию деформирования толстостенных труб, полагая, что внутренний радиус трубы бесконечно мал. Отнесем арматурный стержень к цилиндрической системе координат r, , z и обозначим его внутренний радиус R1, а внешний радиус R2. На стержень действует внешнее давление – P2 и продольная сила N. Схема загружения стержня показана на рис. 2.37.
Рис.2.37. Схема загружения стержня
194
В силу принятых допущений распределение концентрации водорода по сечению стержня описывается уравнением диффузии:
C |
|
|
C |
|
D C |
|
|
|
t |
|
|
D |
|
|
r r |
, |
(2.101) |
|
||||||||
|
r |
r |
|
|
|
|||
где C – концентрация водорода в точке стержня; t – время; D – переменный коэффициент диффузии водорода, учитывающий вид напряженного состояния.
Граничные условия для уравнения диффузии (2.101) записываются в следующем виде:
а) среда действует снаружи
C r R2 C,
C |
0 |
(2.102) |
r |
|
|
r R1 |
|
где С – равновесная концентрация.
Начальное условие принимается в виде C |
t 0 |
0. |
Геометрические соотношения, связывающие |
перемещения и дефор- |
|
мации в каждой точке стержня, записываются в виде |
|
||||
r dw |
, |
w |
, |
z e const, |
(2.103) |
dr |
|
r |
|
|
|
где r, , z – компоненты тензора деформаций; w – радиальное перемещение.
Уравнение равновесия и уравнение неразрывности деформаций принимаются в виде:
d r |
r |
0, |
(2.104) |
|
dr |
|
r |
|
|
d |
|
(2.105) |
||
dr |
r , |
|||
|
|
|
||
где r, , z – компоненты тензора напряжений.
Уравнения, связывающие компоненты тензоров напряжений и деформаций, записываются следующим образом:
|
|
|
1 |
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
|
|||||||
|
r |
|
u , ,C |
r |
|
|
|
z |
|
|
|
|||
|
1 |
|
|
r |
z , |
(2.106) |
||||||||
|
|
|
||||||||||||
|
|
|
u , ,C |
|
|
|
|
|
|
|
|
|||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
|
|||||||
|
z |
|
u , ,C |
z |
|
|
|
r |
|
|
|
|||
|
|
|
|
195 |
|
|
|
|
|
|
|
|
||
где – коэффициент Пуассона: |
1 |
|
1 2 0 u ,,C ; |
|
C |
|
(2.107) |
||
|
2 |
|
2E0 |
|
( u, ,C) – функция интенсивности напряжений и других |
параметров; |
|||
0, Е0 – коэффициент поперечной деформации и модуль упругости материала оболочки в начальный момент времени соответственно; u, u – интенсивности напряжений и деформаций, причем
|
|
|
1 |
|
|
|
|
2 |
|
|
e 2 |
|
|
e 2 |
1/2 . |
(2.108) |
|
|
|
|
|
|
|||||||||||
|
u |
|
2 1 |
r |
|
|
|
|
r |
|
|
|
|
|
|
|
Из третьего уравнения (2.106) находим z:z r e
подставляя которое в два остальных, получим:
|
1 2 |
|
|
|
e, |
|
r |
|
r |
|
|
||
|
1 |
|||||
|
|
|
|
|||
|
1 2 |
|
|
|
e. |
|
|
|
|
|
r |
||
|
1 |
|||||
|
|
|
|
(2.109)
(2.110)
Подставив (2.100) в уравнение неразрывности деформаций (2.105),
получим: |
|
|
|
|
|
|
r d r |
|
|
||||||
|
|
|
|
r |
, |
|
|||||||||
|
|
|
|
|
|
|
|
dr |
|
(2.111) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|||||
d |
|
|
|
|
d |
|
|
|
d 2 |
||||||
|
2 |
r |
r |
r , |
|||||||||||
|
|
|
|
|
dr2 |
||||||||||
dr |
|
|
dr |
|
|
|
|||||||||
аналогично: |
|
|
|
|
1 2 0 |
|
|
|
|
||||||
|
d |
|
. |
|
(2.112) |
||||||||||
|
|
|
|
|
|
||||||||||
|
dr |
|
|
|
2E0 |
|
|
|
|||||||
С достаточной общностью можно принять зависимость между интен-
сивностями напряжений и деформаций в следующем виде |
|
u A C um B C un , |
(2.113) |
где A(C), B(C) – функции, учитывающие изменение диаграммы деформирования при изменении поля концентрации водорода по сечению стержня.
Поэтому функция (u,,C) имеет вид |
|
|
||||
u ,,C |
u |
A C um 1 |
B C un 1 1 m C, , |
(2.114) |
||
|
u |
|
|
|
|
|
C , |
|
|
1, 0 |
|
(2.115) |
|
|
|
0 a , 0 , |
||||
|
|
|
1 / 1 C |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
196 |
|
|
– параметр, характеризующий жесткость схемы напряженного состояния,0 – среднее напряжение, и – интенсивность напряжений.
Выражение (2.115) предложено на основе анализа экспериментальных данных по влиянию концентрации водорода и схемы напряженного состояния на механические свойства материалов.
Основное разрешающее уравнение деформирования арматурного стержня окончательно может быть записано в следующем виде:
d 2 r S u |
,C d r |
U u ,C r F u ,C , |
||||||||||||||||||
dr2 |
|
|
|
|
|
|
|
dr |
|
|
|
|
|
|
|
|
||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
S |
3 |
|
|
2 |
|
|
|
|
|
, |
|
|
|
|
|
|||||
r |
1 2 |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
2 2 |
2 / |
|
|||||||
U |
|
1 4 |
|
|
|
|
|
1 |
|
, |
||||||||||
|
|
|
|
|
|
|
|
|
|
r |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 2 |
|
|
|
|
||
F |
|
|
|
e |
|
|
|
|
. |
|
|
|
|
|
|
|
|
|||
|
r |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Штрих в (2.117) означает дифференцирование по r. Граничные условия
r 01 при r R1,
r P2 при r R2.
(2.116)
(2.117)
(2.118)
Продольная деформация e находится из уравнения равновесия стержня в осевом направлении
|
|
R2 |
|
|
|
||
|
|
N 2 |
zrdr |
|
|
(2.119) |
|
|
|
R1 |
|
|
|
|
|
и определяется выражением |
|
|
|
|
|||
|
N |
R2 |
|
|
R2 |
|
|
|
|
|
|
||||
e |
|
r rdr |
rdr. |
(2.120) |
|||
2 |
|||||||
|
R1 |
|
|
R1 |
|
||
|
|
|
|
|
|||
Коэффициент диффузии D в общем случае является функцией от интенсивности напряжений, концентрации водорода и параметра повреж-
денности П: |
|
D D u ,C,П . |
(2.121) |
197
Кинетическое уравнение, описывающее процесс накопления повреждений, должно соответствовать реальной картине разрушения арматуры.
Анализ экспериментальных данных по коррозионному разрушению арматуры показывает, что этот процесс может описываться кинетическим уравнением накопления повреждений вида:
dП dt А(С) э k (С) / 1 П k (С) , |
(2.122) |
где A(С), k(C) – коэффициенты, определяемые из экспериментальных данных, позволяющие учесть влияние концентрации водорода в точке стержня на кинетику процесса коррозионного растрескивания; э – некоторое эквивалентное напряжение.
2.6.4. Применение теории длительной прочности А.Р. Ржаницына для моделирования коррозионного растрескивания арматуры
Коррозионное растрескивание будем трактовать как постепенное снижение предельного сопротивления материала, вызванное влиянием коррозионной среды и напряженного состояния. По аналогии с [78, 79] для описания кинетики коррозионного растрескивания введем параметр мгновенной прочности r, изменяющийся во времени и равный напряжению, при котором происходит разрушение арматуры в данный момент времени.
Вначальный момент времени r=rо= пр0, а в момент разрушения r(tр) = .
Впромежуточных состояниях мгновенная прочность уменьшается по закону:
r F r, ,a , |
(2.123) |
где а – вектор параметров, отражающих влияние агрессивной среды и других факторов на кинетику коррозионного растрескивания. Вместо точной функции F(r, ) для приближенного решения задачи достаточно задаться каким-либо упрощенным ее выражением, правильно отражающим качественные свойства этой функции. Например, уравнению (2.123) можно придать следующий вид [78]:
r A n e r . |
(2.124) |
Вэтом уравнении переменные разделяются, и его интеграл получается
взамкнутом виде:
|
t |
|
e r |
A ne dt e r0 . |
|
|
0 |
|
|
|
|
Для того чтобы отсюда найти время разрушения, следует положить r(tр) = (tр),
198
что дает:
e |
(t |
p |
) |
r |
tp |
|
|
|
|
e 0 |
A n (t) e (t)dt . |
(2.125) |
|||
|
|
|
|
||||
|
|
|
0 |
|
|
||
|
|
|
|
|
|
|
|
В частном случае при =const получим: |
|
||||||
|
|
|
|
tp |
e r0 e r |
. |
(2.126) |
|
|
|
|
|
|||
|
|
|
|
|
A n e |
|
|
Недостатком уравнения (2.124) является то, что получаемая из него кривая длительного сопротивления (2.126) асимптотически приближается к линии = 0, а не к = дл – аналогу предела длительной прочности.
Закон изменения r может быть принят в виде:
dr А em (r |
r) (r |
)1 ;r(t 0) r . |
(2.127) |
|
dt |
0 |
0 |
0 |
|
|
|
|
|
|
Так как кривая коррозионного растрескивания получена при постоянном напряжении, то, интегрируя (2.127) при = const и учитывая начальное условие, получим:
|
e |
m 1/ 1 |
. |
(2.128) |
r ,t r0 r0 A 1 |
t |
Подставив условие разрушения r(tp) = (tp) в это выражение, получим формулу, связывающую уровень действующего напряжения с временем до разрушения:
|
tp |
B e m , |
(2.129) |
где |
B 1/ |
A(1 ) . |
(2.130) |
Здесь А, , m – коэффициенты, определяемые по экспериментальным данным.
Используя зависимость (2.129) по известным коэффициентам B и m можно определить долговечность арматуры при заданном напряжении. Однако применимость зависимости (2.129) ограничена требованием постоянства уровня напряжений, поэтому для решения задачи прогнозирования долговечности арматуры при меняющихся напряжениях необходимо применять уравнение (2.127), в котором для определения значения коэффициента следует использовать результаты эксперимента по коррозионному растрескиванию при какой-либо программе изменения напряжения.
Для идентификации уравнения кривой коррозионного растрескивания и нахождения значений коэффициентов B и m зависимости (2.129) используем экспериментальные данные из табл. 2.32.
199
Применяя логарифмическое выравнивание и метод наименьших квадратов получим значения B и m при разных температурах Т и концентрациях агрессивной среды С (табл. 2.36).
|
Значения коэффициентов В и m |
Таблица 2 . 3 6 |
||
|
|
|
||
|
|
|
|
|
Для значений a |
Т=70 C |
Т=90 C |
|
Т=110 C |
С=84,5 % |
0,004331297 |
0,00432554 |
|
0,00384487 |
С=65 % |
0,00438284 |
0,00428728 |
|
0,00372849 |
С=39 % |
0,00443965 |
0,00447273 |
|
0,00405401 |
Для значений B |
|
|
|
|
С=84,5 % |
1171,758612 |
232,600691 |
|
40,316679 |
С=65 % |
1533,20903 |
277,663012 |
|
55,1153316 |
С=39 % |
2113,63531 |
391,020021 |
|
103,594668 |
Для значений |
|
|
|
|
lnB |
|
|
|
|
С=84,5 % |
7,066260986 |
5,44932321 |
|
3,69676525 |
С=65 % |
7,33511823 |
5,62640819 |
|
4,00942793 |
С=39 % |
7,65616464 |
5,96875876 |
|
4,64048586 |
Графические зависимости lnВ(Т,С) и m(Т,С), построенные по данным таблицы приведены на рис. 2.38а и 2.38б.
а
б
Рис. 2.38. а – график lnB; б – график m(Т,С)
200
