Метод прогнозирования усталостной долговечности в условиях квазимонотонной деградации физико-механических свойств объекта

ных которого следует, что максимум скорости выщелачивания оксидов марганца (MnO2 и Mn2O3) находится при рН = 1.6±0,2.

Таким образом, оптимальный режим растворения диоксида марганца приходится на рН 1,6±2, соотношение концентраций ионов марганца к оксалат-ионам 1:5, температура 353°К.

Выводы

Лимонная и щавелевая кислоты могут быть использованы для растворения диоксида марганца как в чистом виде, так и в смеси с серной кислотой. Скорость растворения возрастает с увеличением концентрации цитрат- и оксалат ионов. Порядок по оксалат - ионам выше, поэтому использование растворов щавелевой кислоты для растворения MnO2 предпочтительнее. Оптимальное значение рН для проведения процессов растворения в щавелевой кислоте 1,6±0,2 при соотношении концентраций ионов марганца и оксалат ионов не менее 1:5.

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

Работа выполнена при поддержке договора № П 205 Программы: «Научные и научно-педагогические кадры инновационной России на 2009-2013 годы».

Литература

1. Куприков Ю.П. Производство тяжелых цветных металлов из лома и отходов. Харьков, изд. «Основа» при Харьковском гос. университете, 1992, 269c.

2. Пат. N 2164955 Российская Федерация, МПК7 С 22 В 7/00, 19/00, 47/00 Способ утилизации отработанных химических источников тока / Птицын А.Н., Галкова Л.И., Ледвий В.В., Скопов С.В., заявитель и патентообладатель ОАО «Елизаветинский опытный завод». - № 99115669/02, заявл. 14.07.1999; опубл. 10.04.2001, Бюл. № 10. - 14с.

3. Blesa M.A., Morando P.J., and Regazzoni A.E. Chemical Dissolution of Metal Oxides. CRC Press, Boca Raton, FL., 1994. 310p.

4. Н.А. Кишкина, И.Г. Горичев, Б.Е. Зайцев, А.Д. Изотов, В.С. Духанин Влияние 1-гидроксиэтилидендифосфоновой кислоты на кинетику и механизм взаимодействия Mn3O4 с растворами серной кислоты. Журнал физической химии, 1999, Т. 73, №4, с. 653660.

5. Марченко З., Бельцежак М. Методы спектрофотометрии в УФ и видимой областях в неорганическом анализе/ Пер. с польск. - М.: БИНОМ. Лаборатория знаний, 2007, 712с.

6. Дельмон Б. Кинетика гетерогенных реакций. М.: Мир, 1972. 554c.

7. Браун М., Доллимор Д., Галвей А. Реакции твердых тел. М.: Мир, 1983. 360c.

8. Барре П. Кинетика гетерогенных реакций. М.: Мир, 1976. 399c.

Метод прогнозирования усталостной долговечности в условиях квазимонотонной деградации физико-механических свойств объекта

д. т.н. проф. Боровских В. Е., Подвойский А. О. Саратовский государственный технический университет + 7-937-221-18-75, podvoiskijalex@mail.ru

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

Ключевые слова: стохастический процесс; деградационный процесс; вычет предела выносливости; энергия выброса.

Проблемам теории управления ресурсом посвящено большое количество монографий и обзоров, но, несмотря на это, задача прогнозирования усталостной долговечности объекта в

поле стохастических напряжении остается весьма актуальной и сегодня.

Современная наука о прочности насчитывает несколько десятков методов оценки ресурса при стохастических напряжениях. Аналитический аппарат этих методов чрезвычайно разнообразен: одни основываются на многомерных пуассоновских потоках, другие - на марковских процессах и уравнении Фоккера-Планка-Колмогорова, третьи - на понятии кривой усталости и линейной гипотезе накопления усталостных повреждений (гипотезе Пальмгрена-Майнера).

Существующие методы прогнозирования учитывают довольно широкий спектр факторов, влияющих на способность объекта (по ГОСТ 13377-75 понятия «система» и «элемент» можно объединить единым термином объект) сопротивляться усталостному разрушению, однако прогностический алгоритм этих методов основывается, главным образом, на стационарной гипотезе (гипотезе о неизменности механических свойств), а между тем известно, что реальные системы представляют собой сложные стохастические системы, свойства которых изменяются с течением времени, а случайные процессы, совершающиеся в этих системах, являются «существенно нестационарными».

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

Ю

Л m = 0

У\\\ — > 0

> OOj \ \ Y\~

od3 > аз2 Ig^c

lg AT

Рисунок 1 - Изменение кривых выносливости по мере накопления в материале усталостных повреждений (при = const): - априорное значение предела

выносливости; Ю - мера усталостного повреждения; N - число циклов; No - абсцисса

точки перегиба кривой усталости

Замечание: в настоящей работе параметры кривой выносливости определялись не по кривой усталости Велера, а по модифицированной кривой усталости Д.И. Беренова, учитывающей зависимость предела выносливости от числа циклов до разрушения.

Из многочисленных экспериментальных исследований известно, что по мере накопления усталостных повреждений предел выносливости постепенно понижается (рисунок 2) [1].

Для учета этого изменения рекомендуется заменять истинное его значение расчетной величиной, равной 0,5-0,7 первоначального значения предела выносливости. Вместе с тем уровень снижения предела выносливости зависит от спектра нагрузок: чем большая часть этого спектра находиться выше предела выносливости, тем больше следует снижать его расчетное значение.

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

Рисунок 2 - Изменение кривых выносливости по мере накопления в материале усталостных повреждений

Однако прямое экспериментальное выявление закономерности усталостного повреждения для конкретных материалов и конструкций связано с почти непреодолимыми техническими трудностями [1, с.140].

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

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

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

Для решения указанной задачи рассмотрим элементарную площадку в условиях сложного напряженного состояния в некоторой окрестности опасной точки конструкции. Под опасной точкой конструкции понимается точка, в которой поле напряжений реализует максимум, а поле энергии прочности - минимум.

Предположим, что на основе априорной информации о статистических закономерностях нагруженности конструкций-аналогов определены вероятностные характеристики логарифмически нормального закона распределения, аппроксимирующего реальный закон распределения амплитуд деформаций 8х (/), 8у (/) , 8у (/) [2, 3]. Другими словами, имеется пакет осциллограмм, зарегистрированных в эксперименте, ординаты которых подчиняются логарифмически нормальному закону распределения (распределению-аппроксиманту).

Зная составляющие трехмерного 8 -процесса (т.е. располагая оцифрованными осциллограммами, снятыми с розеток тензорезисторов), можно с помощью известных соотношений линейной теории упругости определить так называемые псевдонапряжения [1]

Ч У) = Е -8 х «)

Ч

Ч(0 = Е - 8у( О

где Е - модуль упругости.

Тогда зависимость между напряжениями и деформациями примет вид [1]

1

а

а y =

1—ц2 1

1—ц2 1

Оа + ц  гау)

2 •

где Ц - коэффициент Пуассона.

Расчетное напряжение можно определить следующим образом [1]

а р = А •а х + В •а у + С • т,

где: А, В, С - коэффициенты, зависящие от углов наклона элементарной площадки (рисунок 3).

У У У

Рисунок 3 - Частные случаи ориентации элементарной площадки

Эти коэффициенты определяются по формулам [1]: для площадки 1

. 2 — X 2 — X ~ X — 1 -

А =-+--cos 2а +--sin 2а

2 2 2

D 2 — X 2 — X _ X — 1

B =----cos 2а---sin 2а,

2 2 2

C = (X — 1)  cos 2а — (2 — X)  sin 2а

для площадки 2

B

для площадки 3

—  cos 2p — 1  sin 2p; A = C = 0,

, 2 — X 2 — X А =-:-+-:--cos 2y +

2

X — 1 .

sin 2y; B = C = 0,

2 2 2

где: X - коэффициент, определяющий чувствительность материала к нормальным и каса-

С—1

тельным напряжениям, X

т—1

Если, например, за расчетное напряжение принимается октаэдрическое касательное напряжение, то т-1 ~ 0,7 • а-1.

Расположение опасной площадки и расчетное усталостное повреждение определяются из условия максимума накопленного в ней усталостного повреждения [1] 164 Известия МГТУ «МАМИ» № 2(10), 2010.

ю = шах[ю(а), ю(в), ю(у)], аДу

а углы определяются из уравнений

( ю

( ю ( ю

0.

( а ( в ( у

В первом приближении углы наклона площадки можно определить из условия максимума дисперсии расчетного процесса нагружения [1]

((Б а

(б г

( а

0,

( в ( у

где Б а - среднее квадратическое отклонение расчетного а р -процесса.

Трехмерный случайный процесс {*ах(/),'оу(/)/ау(^)} полностью характеризуется матрицей автокорреляционных и взаимных корреляционных функций [1]

кх

х у'

Г ух

где: кх (т), ку (т), ку (т) - автокорреляционные функции стационарных случайных процес-

сов;

Гух(т) ,...,Гуу(т) - взаимные корреляционные функции стационарно связанных случайных процессов.

Корреляционная функция расчетного а р - процесса определяется следующим образом

[1]

ка (т) = А2 - кх (т) + В2 - ку (т) + С2 - ку (т) + 2 - А - С1 - Гху(т) + ...

... + 2 - А - В - Гху(т) + 2 - В - С

Г у у

где: А1 =

А + В - ц

С

1 -ц

метры.

2-

, Вх =

С

А - ц + В 1 -ц2 - 2 -

С1

С

1 + ц

(10) пара-

Положив в уравнении (10) т = 0 , найдем дисперсию ур. Ба2 = ка(0) = А? - + В2 - Бу2 + С? - Б2 + 2 - А - С - Г™ + 2 - А1 - В - Г™ + 2 - В, - С1 - Гуу,

_>у I иу I ^ ^ »ху I ^ »ху I -Ь- 'ур

где: Гху, Гху, Гуу - корреляционные моменты.

Подставив (11) в (8), получим соотношение для определения угла наклона площадки

[1]

А1 - В[А1] - + В1 - ОД] - б2 + С1 - ДС^] - б2 +

где: - дифференциальный оператор первого порядка по углу наклона а ^ = {а, в, у} .

Ниже приведены реализации составляющих расчетного а р (/) -процесса с логарифмически нормально распределенными компонентами.

Рисунок 4 - Реализация 'ах (/) -процесса

Рисунок 5 - Реализация 'ау (/) -процесса

Рисунок 6 - Реализация 'а (/) -процесса

Рисунок 7 - Реализации: 1 - а х

Рисунок 8 - Реализация а p(t) -процесса

Замечание: в данном случае удобнее работать не с ансамблем реализаций расчетного а p(t) -процесса, а с одной, но репрезентативной реализацией, поэтому предположим, что

а p(t) -процесс обладает свойством эргодичности; для справедливости эргодической теоремы необходимо и достаточно, чтобы автокорреляционная функция случайного а р (t) -процесса удовлетворяла условию [4]

1 г ( И

lim — -11l--- ка

т T 0 ^ T) aW v ;

На практике необходимое и достаточное условие (13) заменяют более простым достаточным [4] (автокорреляционная функция должна стремиться к нулю при т ^ ^ ).

Рисунок 9 - Автокорреляционные функции: 1 - кY

Рисунок 10 - Взаимные корреляционные функции: 1 - ГуY

Рисунок 11 - Автокорреляционная функция а р

Компоненты корреляционной матрицы определялись по формуле [4]

г шТ ^ 1 »-»

к

ху

П

П - ш

I [х+ш)-шхЦу)-шу], ш < 1/4 • п,

г = 1

где: ш , ш - оценка математического ожидания; п - число разбиений.

Явление деградации механических свойств объекта как результат исчерпания предела выносливости можно описать в первом приближении с помощью соотношения, предложенного А.С. Гусевым.

Согласно этому предположению, скорость исчерпания предела выносливости зависит от скорости накопления усталостных повреждений и априорного значения предела выносливости [1]:

а-1,ю = а-1,0 • (1 - ю)^

где: а-1 ю - текущее значение предела выносливости; а-ю - априорное значение предела выносливости; ю - мера усталостных повреждений; ^ - постоянная материала.

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

При больших длительностях разрушение происходит при незначительных пластических деформациях. Металлографические исследования показывают, что разрушению предшествует образование микротрещин перпендикулярно оси действия первого главного напряжения. С течением времени трещины растут и происходит разрыв [5].

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

Скорость распространения трещины есть функция среднего напряжения, температуры и структурного состояния материала.

Следовательно, можно записать [5]

( ю „ ч

-— = и,

ш

где: и - скорость роста трещины; t - условное время нагружения.

Если трещина имеет удлиненную форму и распространяется в одном направлении, то величина и равна скорости роста фронта трещины с точностью до численного множителя. В случае округлой трещины [5]

— = ехр[1пи + 'в • 1п ю],

где 'в - постоянная материала.

Относительно функции и можно делать различные предположения, опираясь на теоретические и опытные факты. В случае расчета повреждаемости при длительном усталостном нагружении и (а) логично представить функцией, зависящей от величины циклических напряжений ау с учетом уменьшения действующего сечения элемента за счет роста повреждаемости. В этом случае уравнение (17) запишем в виде [5]

( ю

= А -

а

а в-

где: N - число циклов; А, П - эмпирические постоянные; а в - временное сопротивление материала.

Решая уравнение (18) при Ю(N=0) = 0 , получим [5]

' А

ю(]') = 1 -

1 -

{ [а(])]сdN

1/

,(0 ' < N

(п + 1) -аВ 0 .

Зная историю нагружения, можно вычислить повреждаемость в любой момент времени. На стадии проектирования можно принять ю = 1 [6]; в этом случае параметры уравнения (8) определяются следующим образом

п = т, А = с -1 - [ав -

где: т - котангенс угла наклона левой ветви кривой усталости, описываемой уравнением

а2 - N = с .

Возвратимся теперь к выборочной функции ар -процесса: согласно Когаеву В.П. [7], амплитуды а а < ка • а-1 0 (по опытным данным можно принять ка = 0, 5) не оказывают

влияния на процесс накопления усталостных повреждений.

Учитывая это обстоятельство можно предположить, что вычет предела выносливости

Ла-1 (ВПВ) (разность, составленная из значения предела выносливости, вычисленного на шаге к и значения предела выносливости, вычисленного на шаге к + 1) связан с энергией выброса ^^ (ЭВ) (энергией напряжений, превышающих текущее значение предела выносливости).

Подставив формулу (19) в соотношение (15) с учетом ка, получим (рисунок 12):

а-1,юN') = к -а-1,0

А

(п + 1) - а В

у /

| [а(])]П(]

Рисунок 12 - Деградационный процесс (рассчитан по формуле (21) для суммы

гармоник)

Тогда предел выносливости будет исчерпываться по правилу (рисунок 13):

к-1 к-1

.

Да = а -0) • П (< 1) У • [1 - (< 1) к ] а П (< 1) j к 1 = 1,

(к=1,...,4) у=1 у=1

а

-. - а7 •П (< 1)У,

У=1

(< 1),

ш+1

а

12

Л

а ш1,0 • Nс

, w > wth,

1

а

12

а ш1,0 • N0

, w < Wth,

I ш+1 |а х 0 )>а-1-1)

где: Да-1) - вычет предела выносливости; а

а -1 - текущее значение предела выносливости;

а-1) - априорное значение предела выносливости, а-1) = ка - а_ю; (< 1) к - выражение меньшее единицы;

11,12 - точки пересечения процессом а р текущего значения предела выносливости; к - текущий номер шага; # - число шагов; а - параметр; Wth- пороговое значение ЭВ.

Таким образом, если текущее значение ЭВ превышает пороговое Wth, то предел вы-

носливости а -1 мгновенно снижается на величину, определяемую по формуле (22), если же условие W > Wth не выполняется, то предел выносливости снижается только по окончании выброса (обозначим его через Wb ), энергия которого превышает пороговую Wth (важно заметить, что учитываются и выбросы предшествующие Wb ).

Рисунок 13 - К пояснению правила снижения предела выносливости: 'ту - относительное время пребывания объекта в у -ом состоянии

В итоге получается траектория деградационного процесса П ('ту), описывающего квазимонотонное ухудшение параметров качества системы в ходе ее эксплуатации (рисунок 14).

Рисунок 14 - Реализация дискретного процесса исчерпания предела выносливости

Замечание: для того чтобы определить средний ресурс объекта с помощью метода условных функций надежности В .В. Болотина, необходимо знать плотность распределения ординат деградационного (обратного кумулятивного) процесса. К сожалению, даже в случае простого кумулятивного процесса не всегда удается определить плотность распределения, поэтому в настоящей статье предлагается определять ресурс объекта как расстояние от начала координат до точки пересечения деградационного процесса с осью абсцисс. Для сокращения машинного времени можно ограничиться рассмотрением нескольких сотен или тысяч ступеней процесса, а затем с помощью метода Монте-Карло достроить деградационный процесс, т. е. аппроксимировать ступенчатый процесс прямой линией с углом наклона

1 - аг^ (Ла ^/Чк).

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

Кроме всего прочего, в данном случае отпадает необходимость составлять кинетическое уравнение (в настоящее время не существует четких рекомендаций относительно правой части кинетического уравнения!), стохастическое дифференциальное уравнение Ито, определять коэффициенты сноса, диффузии, решать уравнение Фоккера-Планка-Колмогорова, которое, как известно, имеет относительно простое решение лишь в ряде частных случаев, не представляющих особого интереса для практики.

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

Выводы

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

Литература

1. Гусев А.С. Сопротивление усталости и живучесть конструкций при случайных воздействиях. - М.: Машиностроение, 1989. - 248 с.: ил. - (Б-ка расчетчика/ Ред. кол.: Н.Н. Мали-нин (пред.) и др.).

2. Дмитриченко С.С., Борисов Ю.С., Русанов О.А. Накопление повреждений и характеристики сопротивления усталости узлов и деталей мобильных машин// Тракторы и сельскохозяйственные машины. - 2003, №8.

3. Ефремов Л.В. Практика инженерного анализа надежности судовой тех-ники. — Л.: Судостроение, 1980. — 178 с.

4. Пугачев В. С. Теория случайных функций и ее применение к задачам автоматического управления. М.: Физматгиз, 1960. - 883 с.

5. Питухин А.В. Вероятностно-статистические методы механики разрушения и теории катастроф в инженерном проектировании. Петрозаводск: изд-во ПетрГу, 1998. - 304 с.

6. Сильвестров И.Н. Расчет ресурса и длительной прочности с использованием критерия повреждаемости// Проблемы машиностроения и надежности машин, 2006, №6, С. 116-118.

7. Когаев В.П. Расчеты на прочность при напряжениях, переменных во времени. М., 1977. 232с.

Численный анализ возбуждаемых трением релаксационных автоколебаний рулевого электропривода при переменной скорости входного воздействия

к. т. н. Михайлова В. Л. МГТУ "МАМИ" (495) 223-05-23, доб. 1318

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

Ключевые слова: рулевой электропривод; численное моделирование релаксационных автоколебаний

В работе [1] представлены результаты численного моделирования фрикционных автоколебаний релаксационного типа в следящем (рулевом) электроприводе (с положительной обратной связью по моменту на валу электродвигателя) в случае постоянной скорости входного воздействия. Настоящая статья имеет целью на примере того же привода проанализировать влияние непостоянства скорости входного воздействия на картину подобных автоколебаний.

Рассматриваемый привод, функциональная схема которого представлена на рисунке 1,

должен обеспечивать поворот руля на угол в соответствии с заданной величиной управляющего воздействия.

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