CC BY
63
Серiя: Техшчш науки
УДК 629.423
Воропай В.С. *
ПРОГНОЗИРОВАНИЕ СРОКА СЛУЖБЫ ВАГОНОВ-ЦИСТЕРН
В статье представлены обобщенный показатель величины расхода остаточного ресурса вагона-цистерны и специализированные критерии, учитывающие действие на элементы базовых частей вагона коррозии, изменения механических характеристик металла, много- и малоцикловой усталости, возможности хрупкого разрушения и условия ползучести. Предложено использование метода экстраполяции при прогнозировании допускаемого срока службы вагонов-цистерн по обобщенному показателю величины расхода остаточного ресурса.
Ключевые слова: параметр технического состояния, аппроксимация, остаточный ресурс, метод экстраполяции, процедура сглаживания, прогнозная величина.
Воропай В.С. Прогнозування строку служшня вагошв-цистерн. У статт1 представленi узагальнений показник величини витрати залишкового ресурсу вагону-цистерни i спецiалiзованi критерп, що враховують дiю на елементи базових час-тин вагону корозИ, змти мехатчних характеристик металу, багато- i малоцикло-вт втомi, можливостi крихкого руйнування i умови повзучостi. Запропоновано ви-користання методу екстраполяцИ' при прогнозувант термту служби ваготв-цистерн, що допускаеться, по узагальненому показнику величини витрати залиш-кового ресурсу.
Ключовi слова: параметр техтчного стану, апроксимащя, залишковий ресурс, метод екстраполяцИ', процедура згладжування, прогнозна величина.
V.S. Voropaj Forecasting of service lives of tank cars. Described in the article were generalized indices of residual resources of tank cars and special criteria, taking into account the action of corrosion upon car elements, as well as changes in metal mechanical characteristic, poly and small cycle fatigue, the possibility of destruction and creeping conditions, The use of the method of extrapolation for forecasting of a possible service lives of tank cars, according to generalized index of the value of th flow of residual resource was proposed.
Keywords: parameters of technical conditions, approximation, a residual resource, an extrapolation method, smoothing procedure, forecasted value look.
Постановка проблемы. С точки зрения диагностирования вагонов-цистерн, наибольший интерес представляет математическое описание характера зависимости параметра технического состояния от наработки. От выбора математической функции зависит качество, достоверность и простота постановки диагноза и прогнозирования остаточного ресурса. Эта функция должна быть универсальной - характеризующей линейную, степенную и другие зависимости изменения параметра от наработки, и простой - содержащей небольшое число коэффициентов и обеспечивающей простоту построения номограмм, таблиц и других справочных материалов. В статье предложено применение метода экстраполяции при прогнозировании допускаемого срока службы вагонов.[1,2,3]
Анализ последних исследований и публикаций. Созданию методического обеспечения проектирования и испытания котлов железнодорожных цистерн, элементов ударно-тяговых устройств, деталей ходовых частей и тормозов посвящены работы В.М.Бубнова, В.Н.Котуранова, В.А.Камаева, М.Б.Кельриха, А.Н.Савоськина, В.Д.Хусидова и др. В этих работах рассматриваются конкретные задачи выбора параметров отдельных узлов вагона на различных стадиях проектирования и модернизации.
Исследованиями, направленными на изучение и обоснование ресурса ответственных элементов металлоконструкции грузового подвижного состава, занимались следующие ученые:
ассистент, ГВУЗ «Приазовский государственный технический университет», г.Мариуполь
Серiя: Техшчш науки
A.А. Битюцкий, В.В. Болотин, Г.М. Волохов, Н.Н. Воронин, Г.И. Герасименко, А.Д. Кочнов, Ю.П. Портнов, А.М. Соколов, М.Г. Сыровцев, В.П. Сычев, А.В. Третьяков, С.В. Урушев,
B.Н, Цюренко, Ю.М. Черкашин и ряд других исследователей. В данных работах проведен достаточно полный структурный анализ состояния универсальных грузовых вагонов. На основании длительного опыта эксплуатации произведена оценка основных повреждений несущих элементов вагонов. Предложен алгоритм оценки ресурса железнодорожного подвижного состава и рассмотрены основные этапы его реализации. Представлены результаты испытаний после ремонта увеличенного объема, которые свидетельствуют о наличие ресурса у вагонов с истекшим назначенным сроком эксплуатации.
Цель статьи - выделение наилучшего описания тренда и определение прогнозных значений путем их экстраполяции .
Изложение основного материала. Использование метода экстраполяции при прогнозировании срока службы вагонов-цистерн основано на том, что рассматриваемый процесс изменения параметров технического состояния представляет собой сочетание двух составляющих: детерминированной и случайной
y(x) = f (а, x) + т( x) (1)
Считается, что детерминированная составляющая f (а, x) представляет собой гладкую функцию от аргумента, описываемую конечномерным вектором параметров |а|, которые сохраняют свои значения на периоде упреждения прогноза. Эта составляющая называется трендом или тенденцией.
Случайная величина составляющей г (x)является величиной некоррелированного случайного процесса с нулевым математическим ожидаемым. Ее оценки необходимы для дальнейшего определения точностных характеристик прогноза
Специфической чертой прогнозной экстраполяции является предварительная обработка числового ряда с целью преобразования его к виду, удобному для прогнозирования, с учетом анализа логики и физики прогнозируемого процесса.
Для предварительной обработки числового ряда прогноза нового срока службы вагона используется процедура сглаживания и выравнивания.
Процедура сглаживания направлена на минимизацию случайных отклонений точек ряда от некоторой гладкой кривой предполагаемого тренда процесса.
В работе сглаживание проводилось с помощью многочленов. В этой процедуре наилучшее сглаживание достигалось для средних точек группы. Так по первым точкам Yj, Y2, Y3, Y4, Y5 сглаживалась средняя точка Y3, затем по следующей пятерке Y2, Y3, Y4Y5Y6 сглаживалась Y4 и т.д. Оставшиеся крайние точки сглаживались по специальным формулам. Применялась наиболее распространенная линейная форма сглаживания, т.е. использовался многочлен первой степени.
Для сглаживания по трем точкам применялись следующие формулы:
~0 =1 (y-1 + Уо + y+i); =1 (5 y-1 +2 Уо + y+i); ~+1 =1 (-y-i +2 Уо +5 у+i ), (2)
3 6 6
где y0, ~0 - значения исходной и сглаженной функции в средней точке;
у-1уУ-1 - значения исходной и сглаженной функции в левой от средней точке; У+1~+1 - значения исходной и сглаженной функции в правой от средней точке.
Формулы и J+ применялись по краям интервала.
Аналогичные формулы применялись для сглаживания ряда по пяти точкам:
~о = -1(у-2 + у-1 + Уо + У1 + у 2); ~-1 = 10(4 у-2 +3 У-1 +2 Уо + У1); ~+1 = 1о(У-1 + 2Уо + 3У1 + 4У2);
~-2 = IT (3 У-2 + 2 У-Ч + Уо - У 2); 2 = 5 (-У-2 + Уо + 2 У1 + 3 У 2 ). (3)
Сглаживание в линейном варианте, как показала практика эксплуатации вагонов, позволяет достаточно точно выявлять величину тренда и определять срок полезного использования.
При формировании прогноза технического состояния вагона учитывались следующие оценки соотношений между временем основания прогноза, временем упреждения и заданной
Серiя: Технiчнi науки
точностью.
Для линейной функции экстраполяции стандартное отклонение прогнозной величины ресурса вагона определялась по формуле:
^р - Бу
Л+1 + (4)
п п - 2 '
п Е (t-t)2
7-1
где Бу - стандартное отклонение исходного ряда от расчетных значений у;
п - число наблюдений ретроспективного периода; ^ - время упреждения плюс время ретроспекции; I - середина ретроспективного периода. Заменив в (4) значения
7 - ПИ- е (I -7)2 - п(п2 -1 и I -Г = п + 2Ь-1 2 ; Е ( ) 12 Ь 2 '
получали
5, - ^ - (5)
Р У У\п п(п2-1)
Тогда доверительный интервал прогноза определялся как
€ + Ь ± 5уК7а, (6)
где 7а - значение I - критерия Стьюдента.
Обозначив К'- К7а, получали §у - БуК'. Таким образом, К' рассматривалось как отношение стандартных отклонений прогнозных значений к отклонениям фактических значений на ретроспективном участке тренда:
К . (7)
5у
Прогнозирование остаточного ресурса вагона-цистерны специализированными критериями выполнялось следующим образом.[3,4] Остаточный ресурс базового элемента вагона, подвергающегося действию коррозии, определялся по формуле:
\(ГЭ) -^^^ (8)
а
где 5ф - фактическая минимальная толщина стенки элемента, мм; 8р - расчетная толщина стенки элемента, мм;
а - скорость равномерной коррозии (эрозийного износа), мм/год; Формула (8) используется, если число замеров N толщины стенок за время эксплуатации вагона не превышает 3. Скорость равномерной коррозии «а» определялась следующим образом. Если после проведения очередного обследования имелось только одно измерение контролируемого параметра Бф (^), полученное при рассматриваемом обследовании, то скорость коррозии определялась по формуле
а - + с - Бф , (9)
где 8и - исполнительная толщина стенки элемента, мм; С0 - плюсовой допуск на толщину стенки, мм;
71 - время от момента начала эксплуатации до момента обследования, лет. Если после проведения очередного обследования имелось два измерения контролируемого параметра Бф 0-2), Бф (^), то скорость коррозии определялась по формуле:
а-ЗДЬ^, (10)
0-2 - 71) К К 2
где Бф (-2), Бф - фактическая толщина стенки элемента при первом и втором обсле-
Серiя: Технiчнi науки
дованиях соответственно, мм; ^ - время от момента начала эксплуатации до момента первого и второго обследования, соответственно, лет; К1 - коэффициент, учитывающий отличие средней ожидаемой скорости коррозии (эрозии) от гарантированной скорости коррозии (эрозии) с доверительной вероятностью у = 0,7 + 0,95 ; К2 - коэффициент, учитывающий погрешность определения скорости коррозии по линейному закону, от скорости коррозии, рассчитанной по более точным (нелинейным) законам изменения контролируемого параметра.
Для большинства грузовых вагонов значение К1 и К 2 можно принимать в пределах
К1 = 0,5-0,75; К 2 =0,75-1,0.
Если после проведения очередного диагностирования имелось три значения контролируемого параметраБф Бф (^)^ф(Ц), полученные при обследованиях в моменты времени /1, ^, /3 , то для определения скорости коррозии а проводились следующие вычисления.
Определялись величины:
3 3 3 3
Б = Е Бф (/,); Б 2 = Е Бф /)/,; ^ = Е ; Х2 = Е (11)
г=1 1 =1 1 =1 1 =1
После чего «а» определялась по формуле:
а = -. (12)
(3X2 - X12)К К2
Во всех трех случаях скорость коррозии - величина случайная, и ее оценку необходимо производить в доверительном интервале прогноза - по зависимости (2), выбирая наиболее неблагоприятный ход развития событий - максимально возможную скорость коррозии.
Как было показано выше, в процессе эксплуатации вагонов, в ряде случаев, происходит снижение механических свойств материалов отдельных элементов. Такое снижение механических свойств выявляется в процессе испытания образцов, изготовленных из контрольных вырезок, либо путем замера твердости металла и определения механических характеристик элементов вагона по существующим корреляционным зависимостям. В случае, если снижение механических свойств оказалось более 5% нормативных, то все расчеты проводились по фактическим механическим свойствам материала, а скорость снижения свойств определялась аналогично скорости коррозии и путем экстраполяции определялись механические свойства материала к концу ожидаемого остаточного периода эксплуатации вагона.
Оценка ресурса вагона-цистерны по критерию мало- и многоцикловой усталости проводилась следующим образом.
Величина эквивалентной приведенной амплитуды динамических напряжений при расчете многоцикловой усталости определялась по дискретной функции плотности распределения амплитуд напряжений по формуле[5]:
Тк
N
^Е (а*)т • Р1 + Nс2 Е (*£)т • Р1 + ^ е )т • Р1
(13)
где а аэ - расчетная величина амплитуды динамического напряжения условного симметрического цикла, приведенная к базовому числу циклов N, эквивалентная по повреждающему действию реальному режиму эксплуатационных случайных напряжений за расчетный срок эксплуатации; т - показатель степени в уравнении кривой усталости в амплитудах; Тк - расчетный срок эксплуатации в годах; N - базовое число циклов; ^123 - число циклов динамических напряжений, действующих на вагон соответственной через автосцепку, от колебаний на рессорах, ремонтных и др.; ^ - амплитуды динамических напряжений в , диапазонах
ударных продольных сил; - амплитуды динамических напряжений от колебаний на рес-
сорном подвешивании (в к диапазонах); - апмлитуды динамических напряжений от испытательного внутреннего давления, ремонтных нагрузок и т.д. (в п диапазонах); Р.| кп - частость возникновения амплитуд при соответствующих напряжениях.
Величина расчетного срока эксплуатации в годах по критерию многоцикловой усталости
а а,э = т
к
п
Серiя: Технiчнi науки
определялась по формуле:
°а,N
(-ТТ)т • N0
Тк =-П-, (14)
N1Е К)т • р] + N С2 Е (*0к)т • Р1 + ^з е (^1)т • РЩ1
] к п
где <7аN - предел выносливости по амплитуде для контрольной зоны при симметричном цикле и установившемся режиме нагружения при базовом числе циклов N ; [ п] - допускаемый коэффициент запаса сопротивления усталости.
При расчете конструкции на малоцикловую усталость , согласно [5] принималось:
N■
Тк—- = 1, (15)
к Г [ N ],■ ' ' '
где N1 - число циклов 1 - го нагружения за 1 год эксплуатации; Тк - расчетный срок эксплуатации, в годах; [К] - допускаемое число циклов ь го нагружения.
Из условия циклической прочности расчетный срок эксплуатации определялся по формуле:
Тк • (16) "[ N ],
Определение остаточного ресурса по критерию хрупкого разрушения (трещиностойко-сти) проводился в следующих случаях. Во-первых, если при проведении дефектоскопии элемента были обнаружены дефекты, близкие к установленным предельным значениям, и ремонт их связан с большими техническими трудностями. Во-вторых, если при проведении дефектоскопии были выявлены отдельные трещины, которые при ремонте были заварены и при повторном контроле было установлено отсутствие этих дефектов. Условие сопротивления хрупкому разрушению проверялось на основании выполнения соотношения:
Кст1 < [Кет1],
где Ка1 - коэффициент интенсивности напряжений;
[ Ка1 ] - допускаемый коэффициент интенсивности напряжений.
При отсутствии информации о дефекте, что является наиболее распространенным случаем для базовых частей вагонов-цистерн, при определении Ка1 принималась условная поверхностная трещина глубиной а = 0,25S и полудлиной С=1,5а [116].
Остаточный ресурс с учетом ползучести материалов (длительной прочности) определялся для вагонов-цистерн, работающих при повышенных температурах, таких как вагон-цистерна для перевозки пека модели 15-1532, когда в расчетах на прочность допускаемое напряжение определяется по пределу длительной прочности или 1% предела ползучести для заданного срока эксплуатации. Если нет таких данных, то температура, когда учитывается ползучесть, принималась равной и выше 380 °С для углеродистых сталей, 420 °С - для низколегированных сталей, 525 0С - для аустенитных сталей.
На основании ретроспективной информации об изменении фактического размера Lф(t)
элемента базовых частей вагона, в местах действия наиболее высокой температуры, где была выявлена остаточная деформация ползучести, ресурс элемента Тп определялся по следующей зависимости:
Тп = —, (17)
ап
где ап - скорость установившейся ползучести, %/год.
Остаточный ресурс ТП элемента в этом случае определялся по формуле:
Т0 = Тп - Тэ (18)
где Тэ - продолжительность эксплуатации от начала до последнего обследования.
Скорость установившейся ползучести определялась по формуле:
Серiя: Техшчш науки
о _ юосЬф^!)-^(12)] (19)
" Lф (11)ДК1 К2 ' ^ }
где Lф Lф (12), - фактический размер контролируемого элемента, мм; Дt - время между первым и вторым обследованием, лет; К1 -коэффициент, учитывающий отличие средней ожидаемой скорости ползучести с доверительной вероятностью у _ 0,7 + 0,95 ; К2 - коэффициент, учитывающий погрешность определения скорости ползучести, рассчитанной по более точным нелинейным законам измерения контролируемого параметра.
Значения коэффициентов К1 и К 2 для грузовых вагонов-цистерн принимались в пределах: К = 0,5-0,75; К 2 =0,75-1,0.
Если после проведения очередного диагностирования имелись три значения контролируемого параметра Lф Lф (12)^ф(13), полученные в моменты времени t1,^.2, t3, то для определения скорости ползучести ап проводились следующие вычисления:
3 3 3 3
Ll _ £ Lф (^); L2 _ £ Lф (^ %; ^ _ £ ^; Х2 _ £ . (20)
г_1 1_1 1 _1 г_1
Скорость ползучести определялась по формуле:
100(LjXJ -3L2)
Ln (X2 -3X2)K,K2
(21)
На основании прогнозной оценки остаточного ресурса вагона и с учетом оперативных комплексных показателей и экономических показателей работы вагона за ретроспективный период времени можно сделать вывод об экономически обоснованных управляющих воздействиях, которые должны быть применены к вагону. Возможными являются следующие управляющие воздействия: пролонгация (продление) срока службы вагона, его модернизация, реконструкция, модификация и, наконец - исключение из инвентаря. Кроме того, проведение прогностических расчетов позволяет определять оптимальный дополнительно назначаемый срок службы вагона-цистерны, после истечения его нормативного срока использования.
Выводы
1. Предложено применение метода экстраполяции при прогнозировании срока службы как метода, при котором происходит выделение наилучшего описания тренда, определение прогнозных значений путем их экстраполяции. Процедура сглаживания в линейном варианте позволяет достаточно точно выявлять величину тренда и определять срок полезного использования.
2. Результаты расчетов коэффициента запаса усталостной прочности n, расчетной эквивалентной величины амплитуды динамических напряжений оа э исследуемых вагонов-цистерн показали, что при п = 1,6 подтверждается наличие у испытуемой цистерны остаточного ресурса, срок ее эксплуатации может быть продлен до 10 лет (т.к. коэффициент запаса усталостной прочности находится в пределах, допускаемый « Нормами...» ), при наименьшем , полученном расчетно-экспериментальным методом п= 2,2, подтверждается наличие у испытуемой цистерны остаточного ресурса свыше 5 лет и срок ее эксплуатации может быть продлен до 5 лет.
Список использованных источников:
1 .Вальд А. Последовательный анализ. - М.:Физмат, 1960. - 328 с.;
2. Галлагер. Методы получения матриц жесткости элементов - М.: В сб.: Ракетная техника и космонавтика. 1963, вып.6. - с. 187-189;
3. Длин А.М. Математическая статистика в технике. - М.:Советская наука. 1958 - 466 с.;
4. Дробов Г.М. и др. Экспертные оценки в научно-техническом прогнозировании. - Киев.: Наукова думка, 1974. - 159 с.;
5.Нормы для расчета на прочность и проектирования вагонов железных дорог МПС колеи 1520 мм (несамоходных). - М.: ГосНИИВ-ВНИИЖТ, 1996. - 315 с.
Рецензент: В.К. Губенко
д-р техн. наук, проф., ГВУЗ «ПГТУ»
Статья поступила 29.11.2010
