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

мобилей с использованием перспективных электрических машин российского производства". Этап № 2 "Р^работка технического проекта комбинированной энергетической установки с комплексной системой энергообеспечения". ГК №16.516.11.6062.

1. Отчёт о НИР по теме: "Р^работка и изготовление комбинированной энергетической установки (КЭУ) в стендовом исполнении для среднетоннажных грузовых городских автомобилей с использованием перспективных электрических машин российского производства". Этап №3 "Изготовление комбинированной энергетической установки в стендовом исполнении". ГК № 16.516.11.6062.

2. Куликов И.А., Селифонов В В., Филонов А.И. Управление комбинированной энергоустановкой автомобиля: экология или энергоэффективность? //^томобильная промышленность - М.: Научно-те^ическое издательство Машиностроение - 2011 - № 1.

3. Карунин А.Л., Бахмутов С В., Селифонов В.В., Круташов А.В., Баулина Е Е., Авруцкий Е.В., Филонов А.И., Карпухин КЕ. Гибридные автомобили - столбовая дорога к экономичному и экологически чистому транспорту. // Журнал ассоциации автомобильных инженеров - Москва - 2007- № 3. - с. 38 - 45.

4. Карунин А.Л., Бахмутов СВ., Селифонов В В., Круташов А.В., Баулина ЕЕ., Авруцкий Е.В., Филонов А.И., Карпухин КЕ. Гибридные автомобили - столбовая дорога к экономичному и экологически чистому транспорту (продолжение). // Журнал ассоциации автомобильных инженеров - Москва - 2007- № 4. - с. 52 - 56.

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

д.т.н. проф. Боровских В.Е., Подвойский А О.

Саратовский государственный технический университет + 7-937-221-18-75, podvoiskijalex@mail.ru

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

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

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

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

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

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

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

lg СТ.

■1,0

со = 0

— COi > 0

GQ2 > OCij \

\ > Сй2 Ig^C k

lgtf

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Ч V ) = Я-Ех Ц )

'оу ^ ) = Е-ъу (?), (1)

) = Е-б ^ )

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

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

1 л ^ л

Ох + ц • ау)

1-Ц2 1

1-Ц3 1

т =

■Сау+Ц -Ох) 5

(2)

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

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

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

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

У

а

и ж {

у

X

1 1 3

_

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

у

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

' 2-А, 2-А, _ Х-1 . 0 А =--1---cos2a +--sm2a

2

2

_ 2-Я, 2-Х _ B =----cos2a-

2 Х-1

2 2 2 С = (X -1) • cos2a - (2 - А,) • sin2a

sin 2а,

(3)

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

B

2-Х 2-Х

2

2

cos

(2-Р)

Х-1

2

sin

(2-р); А = С = 0:

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

2 + —srn(2-Y); В = С = 0,

(4)

(5)

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

—. Если, например, за расчетное напряжение прини-

мается октаэдрическое касательное напряжение, то Т^ « 0,7 • а^ . Расположение опасной площадки и расчетное усталостное повреждение определяются из условия максимума накопленного в ней усталостного повреждения [1]:

ю= тах[со(а),ю(Р),ю(у)], а,Р,У

(6)

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

dra da

do

= 0.

(7)

^а ^Р dy

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

ds„ ds„ ds„

— = 0,

(8)

^ ^ dy

где: sст - среднее квадратическое отклонение расчетного ур -^оцесса.

Трехмерный случайный процесс {(7х ),лоу )} полностью характеризуется

матрицей автокорреляционных и взаимных корреляционных функций [1]:

кх (т) гху (т) г^ (т)

Г ух ф к у ф г„ ф , (9)

/ух СО г .у ф ку

где: кх(т), ку (т), ку(т) - автокорреляционные функции стационарных случайных процессов; гух (т),..., гуу (т) - взаимные корреляционные функции стационарно связанных случайных процессов. Корреляционная функция расчетного ар -^оцесса определяется по формуле [1]:

К (т) = A 2 • кх (т) + B2 • ку (т) + С2 ■ ку (т) + 2 • A-С- ^ (т) +... ... + 2- Ar Br ГхуBi Ci (т),

(10)

где: А =

Л + В-\л

С

Вг =

Л-р + В

С

С

1 + ц

параметры.

1 — М-2 2-(1 + Ю 1-ц2 2-(1 + ц)

Положив в уравнении (10) т = 0, найдем дисперсию СТр .

¿2 = К(0) =Л12 .¿2+ В2-¿2 + С2-s2+2^ А -Сг+ 2ЛгВ1 + 2В1 С-Тл, (11)

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

Подставив (11) в (8), определим угол наклона площадки [1]: Л- Д Л,]- ¿2+ Вг ДВ^ ¿2 + Сг ДС^ *;+( Лг ЩВХ] + Д Лl^ В1). Гху+... ... + (Лг Д[С ] + Д л] • С1). г^ (В г Д[С ] + Д [ В1 ]. С1). туу = 0, где: Д[-] - дифференциальный оператор первого порядка по углу наклона ау- = {а,|3,у}.

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

(12)

шШй

к*

12 16 20 24 28 32 36 40 44 43 52 56 60 64

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

72 76 30

г,с

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

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

Рисунок 7 - Реализации: 1 - стх (?), 2 - ау (?), 3 - )

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

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

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

-ч

з Т J

Нт

/у Т\ 1--

ту

о V

^ (т)<^т = 0.

(13)

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

~ I тТ\ 1 п—т

- =--Е [хЦг+т)-тх]•[)-ту],т< 1/4-п, (14)

V п ) п-т г=1

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

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

КЮ

,1 | 1 - 1 _ ----2

V >у Л, , Л ¿к 1, Л / ИЙг!

о* ■р V/ г '.у- 'ЧУ

0 2 4 6

10 12 14 16 13 20 22 24 26 28 30 32 34 36 38 40

Рисунок 9 - Автокорреляционные функции: 1 - ky.(т) ; 2 - kx (т); 3 - ky (т)

Рисунок 10 - Взаимные корреляционные функции: 1 - гу„ (т); 2 - г (т); 3 - г^ (т)

12 16 20 24 28 32 36 40 44 48 52 56 60 64 63 72 16 80

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

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

а-и^о-а-ю)*, (15)

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

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

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

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

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

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

dm

— = и, (16)

dt

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

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

— = ехр[1п и + Р-Ью], (0 < 'Р < 0,5), (17)

dt

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

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

d(£^

dN

= А-

(18)

ав<1-со)

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

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

( А N1 У/(И+1)

о(N') = 1- 1---Г[ст(N)]" dN , (0< N < N) (19)

I (п + 1)-а£ 00 )

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

п = т, А = С~Х-[а£-(и+1)],

(20)

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

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

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

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

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

гармоник)

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

а-1,И( N ) = ka

'-1,0

1

А

N

V

(п +1) • ст

п

в 0

Дсг( N)]

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

Да(? = < П(< 1)у • [1 -(< 1)k] л П(< 1),

k=1

= 1,

}—1

}— 1

а(?=а!М1(<1) у ,

]—1

(21)

(22) (23)

(<1) k =

т+1 т

Е

а

ГТт

1 (

т+ 1 1

-1) у V

'2

\<УХ ^ )dt

(24)

а

N 0

2

, w < w

А?

-1,0 t1 х

где: Ао,(1) - вычет предела вьшосливости; ст^ - текущее значение предела выносливости; С7-0) - априорное значение предела выносливости, СУ^ = ka^ 0,-10; (<1) k - выражение меньшее единицы; t1, ? 2 - точки пересечения процессом ср текущего значения предела вьшосливости; k - текущий номер шага; ^ - число шагов; а - параметр; wth - пороговое значение ЭВ.

1

0

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

Таким образом, если текущее значение ЭВ превышает пороговое wth, то предел выносливости СУ^ мгновенно снижается на величину, определяемую по формуле (22), если же условие w>wth не выполняется, то предел выносливости снижается только по окончании выброса (обозначим его через wb), энергия которого превышает пороговую wth (в^сно заметить, что учитываются и выбросы предшествующие wb).

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

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

-1-ст^/\).

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

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

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

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

Выводы

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

Литература

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

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

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

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

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

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

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

Моделирование устройства балансировки Li-ton аккумуляторной батареи с коммутируемыми конденсаторами для электромобиля

Варламов ДО. Университет машиностроения 8(916)513-72-16, varlamovd@mail.ru

Аннотация. В статье рассмотрены принцип работы и способ построения схемы моделирования устройства балансировки с коммутируемым конденсатором для Li-Ion аккумуляторной батареи в программе Matlab Simulink. Сделаны выводы о возможности применения устройства балансировки для аккумуляторной батареи