Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Сопротивление усталости и живучесть конструкций при случайных нагрузках

..pdf
Скачиваний:
8
Добавлен:
12.11.2023
Размер:
13.19 Mб
Скачать

Подставив соотношения (16.23) в (16.16), получим

где

s| = 60 + bi sin 4а +

62 cos 4а,

 

(16.24)

 

 

 

 

 

 

bo = 2(1

(0,5sl +

0,5s^ + s$ — КхУ— Куу),

 

=

2 (1 + |i)a <°’5s* - °-5s* +

-

*ет).

(16.25)

^2 =

2(1 -f- Ji)2 (SV+

Ary

Ayy).

 

Введем обозначение t = tg а. Тогда

 

 

 

s? = b0 + bt - j l f a - +

b3-1 ^ - .

 

(16.26)

Проведем исследование функции (16.26) на

экстремум. Из

условия dso/dt =

0 получаем, что дисперсия s® принимает экстре­

мальное значение

при

 

 

 

 

 

 

U . i ^ — ^ ±

V ЪЦЬ\ +

1.

 

(16.27)

Расположение

опасной площадки

будет

определяться

углом

а = arctg tu2.

 

 

 

 

 

 

После подстановки выражения (16.27) в формулу (16.16)

максимальная дисперсия расчетного процесса

 

 

 

max sl = b0 + \ f ь\ + bb

 

 

(16.28)

где b0, bx и bt — величины,

определяемые по

вероятностным

характеристикам составляющих процесса нагружения в соответ­ ствии с формулами (16.25).

Аналогично определяются расположение опасной площадки и максимальная дисперсия расчетного процесса нагружения для

случаев 2 и 3,

показанных на

рис. 16.3. Так, для

площадок,

определяемых углами

р,

 

 

 

 

Al =

2(1—fi») sln 2Р:

 

 

= -

2 ( l - ^ r

Stn 2Р* Cl = 0:

(16>29)

St =

-4 (1 ^ ) »

(Р2^ +

2р/С*„ + 4) sin2 2р.

 

Расположение площадки с максимальной дисперсией расчет­ ного процесса нагружения определяется углом р = 45°. При этом

max s* = 4(i-ii»)»' + 2р К*у + 4)- (16.30)

172

Для площадок, расположение которых определяется углами у,

^ 1 = -2(Г- ^ )

sln2v;

(16.31)

 

 

 

Вг —

2(1 — ца)

sln

Сх — 0;

s? = 4(i

(si +

Р24 +

2p/G„) sin5 2у.

Расположение площадки с максимальной дисперсией расчет­ ного процесса нагружения определяется углом у = 45°. При этом

шах s| = 4(1^ .|t«)« (si + И24 + 2рКхУ)-

(16.32)

Частные случаи расчета. Рассмотрим процесс нагружения {d* (t), ду (/), ду (f)} с компонентами, которые изменяются во времени синхронно и синфазно. В этом случае Кху = sxs„, Кху = = sxsy\ из соотношений (16.16) следует, что среднее квадрати­ ческое отклонение расчетного процесса нагружения в площадках, определяемых углом а,

sa = AiSx -f- BxSy -f- CiSy,

(16.33)

где функции Alt и Cx определяются соотношениями (16.22). Проведя исследование функции (16.33) на экстремум, уста­ навливаем, что величина s0 принимает экстремальные значения в площадках, расположение которых определяется из уравнения

 

ia On =

^ (^ —

sy sx + (2Я — 3) su

(16.34)

 

8

Sy (2Я, — 3) s* +

2 (1 — Я) sv

 

 

При Я =

1

 

 

 

 

 

 

tg 2а =

^

 

(16.35)

При Я =

2

s

sv—sx

 

 

 

 

 

 

 

tg 2а =

sy

sx

(16.36)

 

 

sy + sx 2sv

При детерминированном нагружении расположение площадок с максимальными нормальными и максимальными касательными напряжениями определяется соответственно по следующим фор­ мулам:

tg 2а =

 

(16.37)

tg 2а =

езс

(16.38)

6х -(- ву 2&у

где ех, Е и, Еу деформации, зарегистрированные в некоторый момент времени в прямоугольной розетке тензодатчиков (см. рис. 16.4).

Формулы (16.35) и (16.36) являются вероятностными аналогами формул (16.37) и (16.38).

173

Подставив соотношения (16.35) и (16.36) в формулу (16.33), определим максимальное среднее квадратическое отклонение рас­ четного напряжения:

при X = 2

max зт =

V 2

1

- V (s y -

sv)a + (sx - sv)2;

(16.39)

a

(I +

p)

 

 

при X = 1

 

 

 

 

 

° 4 xs« = ^ i t %

~

К г о '- н о

- s»)' +

(,6-40)

Формулы (16.36) и (16.39) можно также получить из общих соотношений (16.27) и (16.28). Действительно, в рассматриваемом случае

Ьо — 4 (j

[(s* — sv)2 +

(sy ~ sv)2l = V

+

^2 !

=

4 ^j _|_

(Sy

sx) (Sy -f- sx 2sv);

 

(16.41)

 

^2 = 2(1-)- (l)2 (Sx

sv) (sy sv)-

 

 

Подставив соотношения (16.41) в (16.28), получим

 

 

 

maxst = 2&0-

 

(16.42)

 

 

а

 

 

 

Соотношения (16.39) и (16.42) совпадают. Аналогично из фор­ мул (16.41) и (16.27) следует равенство (16.36).

При синхронном и синфазном изменениях компонент процесса нагружения формулы (16.30) и (16.31) принимают следующий вид:

(16.43)

; S * ’ = 2(1 —ц*) f a + !“ »)•

Для дальнейшего расчета необходимо перейти от регистри­ руемого процесса {дх (t), ау (t), ov (t)} к анализу напряженного состояния, определяемого трехмерным процессом {ах (t), <г„ (/),

т (0 ) .

с помощью формул

(16.11),

Этот переход осуществляется

а обратный — с помощью следующие соотношений:

 

Ох — Ох

pOyi

 

■&у = оу — ро*;

(16.44)

dv = (1 + р) т +

(ах + ау).

 

Пример 1. Рассмотрим процесс синхронного и синфазного изме­ нений напряжений ах (t) и ои (t) при т (t) = 0 (рис. 16.5).

174

бх

* 1

Рис. 16.5.

Синхронное и

синфазное из­ Рис.

16.6. Синхронное н

асинфаз-

менения напряжений ох и оу

ное изменения напряжений ох и оу

Пусть

Sox = ksay,

где к — параметр. Обозначим Soy =

s. Из

соотношений (16.44) следует

 

 

 

 

 

Sx =

(k — Р) s;

 

 

 

 

s „ =

(1 —> £ )s;

(16.45)

 

 

=

— fi)(l +

k)s.

 

Подставиввыражения (16.45) в формулы (16.35), (16.36)» (16.39) и (16.40), установим, что расположение площадки с макси­ мальными касательными напряжениями определяется углом а = = л/4, а площадки с максимальными нормальными напряже­ ниями — углом а = 0. При этом

max sx =

| k — 1 | s;

а=Я/4

£

 

1

max sx = -s- s;

е=я/4

г

max sx = 4-fcs.

V-n/4

1

При 0 < k < 1 опасной площадкой по критерию максималь­ ных касательных напряжений является площадка 2, при k > 1 — площадка 3 (см. рис. 16.3).

Рассмотренный пример показывает необходимость учитывать

возможность появления

наибольших

усталостных

повреждений

в площадках 2 и 3.

 

 

 

 

Пример 2. Рассмотрим случай синхронного и асинфазного

изменений напряжений

ох (t)

и сг„ (t)

= kax (t)

при т (/) = 0

(рис. 16.6).

 

 

 

 

В этом случае

 

 

 

 

max sx = 4 - (k +

1) s.

 

а= я/4

*

 

 

Отсюда следует, что опасной всегда будет площадка 1 (см. рис. 16.3).

17Б

Рассмотрим случай, когда процессы ах (i), ау (t) и т (t) яв­ ляются статистически независимыми. Из соотношений (16.44) следует

 

si =

s|x +

р2^

;

 

 

 

 

sl = sby+

\i2slx;

 

 

 

sy = (1 + p)2 si + 4 (1 -

p)2 (sbx +

sby),

К*? = 4

- I*) SK -

 

X

* (1 “

^

SV ’

<1646)

^Kyv----2~(1

P) soy

 

2~ ^ ^

f*) S°x*

 

Kxy = —p (s|x + sby),

 

 

 

где s^, Soy; s| — дисперсии

процессов

ox (t),

ay (/),

x (t).

Подставив соотношения (16.46) в формулы (16.25), (16.27),

(16.28), (16.30) и (16.32), получим

 

 

 

 

 

 

bo =

(4s|

slx -f- s|^); b\ =

0;

 

62 = 4 - ( 4S* - S’* - 4 ) ;

 

 

 

max

s|

 

 

 

sx

 

 

(16.47)

 

 

 

 

 

 

a=0, я/2

 

 

 

 

 

 

 

inax S J —■ Q S(j

 

 

 

 

3 = 4 5 °

 

1

 

 

 

 

 

шах

*

Sty

 

 

 

 

V = 4 5 °

 

z

*

 

 

 

 

Из полученных данных следует, что опасной по критерию

максимальных

касательных

напряжений

площадкой

является

в рассматриваемом случае площадка типа

1 (см. рис.

16.3).

Пример 3. Рассмотрим случай, когда процессы ах (t) и ау (t)

являются статистически

независимыми, а

х = 0. Пусть, как и

в примере 1, Sax = ks„v =

ks,

где k — параметр. Из соотношений

(16.47) следует,

что

 

 

 

 

 

max s? =

-^-(&2 + 1) s2.

 

Пример 4. Рассмотрим случай чистого сдвига с дисперсией касательных напряжений s?. Из соотношений (16.47) следует

max s\ = s|. а=0

176

Рис. 16.7. Корреляционные (а) и взаимные корреляционные (б) функции про­ цессов изменения деформаций в элементе кабины трактора ДТ-75М при дви­ жении по полю с бороздами со скоростью 6,6 км/ч

Методика расчета усталостной долговечности при плоском напряженном состоянии была апробирована на металлоконструк­

циях

рам мобильных машин

типа автомобилей и тракто­

ров

[12,

34].

представлены корреляционные

Для

примера на рис. 16.7

и взаимные корреляционные функции процессов изменения де­ формаций в одном из элементов кабины трактора ДТ-75М при движении по полю с бороздами со скоростью 6,6 км/ч.

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

12 Гусев А. С.

177

Рис. 16.8. Схема возникновения и распростране­ ния усталостной трещины при -кручении

и в большинстве случаев уже не будет совпадать с первоначаль­ ным направлением. Долом элемента конструкции с -появившейся трещиной может происходить в- течение довольно длительного времени, которое определяет степень его живучести. Характер­ ным является разрушение элемента конструкции, находящегося в состоянии чистого кручения (рис. 16.8). На таком элементе сначала появляется направленная вдоль образующей трещина, возникновение которой обусловлено действием касательных на* пряжений. В дальнейшем трещина отклоняется от первоначаль­ ного направления на угол 45° и развивается уже по спирали, что обусловлено действием нормальных растягивающих напряжений.

§17. Использование датчиков повреждений для прогнозирования остаточного ресурса конструкций

Расчетные методы определения накопленного к некоторому моменту времени усталостного повреждения позволяют оценить долговечность конструкции и ее остаточный ресурс, если будут заданы вероятностные характеристики процесса нагружения. При этом усталостная долговечность конструкции определяется по величине усталостного повреждения v = 1, а остаточный ресурс — по остаточному усталостному повреждению v0 = 1 — \ t (v, — усталостное повреждение, накопленное к моменту вре­ мени t).

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

Значительное уточнение расчета индивидуального ресурса конструкций достигается применением специальных датчиков — свидетелей истории нагружения данного экземпляра конструкции, которые работают по принципу изменения электрического сопро­

178

тивления, происходящего с увеличением числа циклов нагру­ жения.

Простейший вариант использования датчиков для диагности­ рования остаточного ресурса состоит в сопоставлении происшед­ шего изменения их электрического сопротивления с критическим значением этого сопротивления, соответствующего моменту раз­ рушения конструкции. Критическое значение определяется для данной серии датчиков при тарировочных испытаниях на образцах в лабораторных условиях. Обычно используются фольговые и полупроводниковые датчики [5, 10). Критическое изменение их электрического сопротивления составляет 3—10 % от их началь­ ного сопротивления. Вид зависимостей AR = / (N , е) (где N — число циклов нагружения; е — амплитуда деформаций в месте установки датчика) показан на рис. 17.1. Крестиками обозначены моменты разрушения образцов. Критическое значение AR Kpможно рассматривать как константу датчиков, значения которой не зависят от истории нагружения. Эта константа имеет некоторый статистический разброс, вероятностные характеристики которого также определяются по результатам тарировочного эксперимента. В результате испытаний обычно обнаруживается порог чувстви­ тельности е* датчиков, который определяет нижний уровень деформаций, при котором еще происходит изменение их электри­ ческого сопротивления по мере увеличения числа циклов нагру­ жения.

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

Процесс изменения электрического сопротивления ф (N) датчи­ ков можно описать следующим кинетическим уравнением:

-aJ-= /(*>. *). (i7.i)

где х — амплитуда деформаций в циклах процесса нагруже­ ния х (N).

Функция / (ф, х) определяется по тарировочным графикам, показанным на рис. 17.1. В рассматриваемом случае доведение образцов до разрушения не требуется.

Непосредственно по тарировочным графикам определяется функция ф = ф (х, N). Дифференцируя ее по N, получаем ф = = ф (*, N). Подставив в полученное соотношение решенйе уравне­ ния ф = ф (х, N) относительно N, получим функцию / (ф, х).

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

процесс ф = ф (N) и вычисляется его производная ф = ф (N) (рис. 17.2.) Тогда для определения амплитуды xt в момент, опре­ деляемый номером цикла нагружения N it достаточно решить уравнение (17.1) относительно х = xt при заданных значениях

ф = ф| и ф = ф(. Повторяя такие вычисления на каждом цикле

12*

179

 

N

 

*1

N

 

 

 

 

Рис. 17.1.

Зависимость сопротивления

Рис. 17.2.

Анализ процесса

на-

датчиков от

числа циклов нагружения

копления

повреждений _

 

нагружения, можно выявить

всю историю нагружения

в виде

последовательности амплитуд

xt (t = 1, 2, 3, ...). Однако

такой

путь регистрации процесса нагружения малоэффективен.

 

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

Будем считать, что функцию / (ф, х) можно представить в виде произведения двух функций (ф) и /2 (х), первая из которых зависит только от ф, а вторая — только от х. Тогда уравнение (17.1) будет дифференциальным уравнением с разделяющимися переменными. Проинтегрировав его, получим

(17.2)

где |Т*1 — значение вычисленного интеграла.

Осредняя обе части равенства (17.2) и полагая, что его левая часть и число циклов N0 имеют малые флуктуации, получаем

со

(17.3)

О

где р (х, а) — плотность распределения амплитуд нагрузок, за­ данная с точностью до неизвестного параметра распределения а.

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

180

Рис. 17.3. Возможные виды процессов Рис. 17.4. Мультипликаторы накопления повреждений

процесс хг (t) = hx (t), где h — коэффициент усиления мульти­ пликатора. На рис. 17.4 показан один из возможных мультипли­ каторов, представляющий собой пластину с локальным сужением. Принцип его работы заключается в том, что датчик устанавли­ вается в зоне концентрации напряжений и поэтому показывает увеличенное значение изменения его электрического сопротивле­ ния. Коэффициент усиления Л устанавливается при тарировке. На рис. 17.5 показаны изменения сопротивлений первого и вто­ рого датчиков. По аналогии с соотношениями (17.2) и (17.3) можно записать для второго датчика

оо

 

g*2 = g (Ф) |£*» = No J f t (х2) р (х„ a)dx.s,

(17.4)

О

 

где Дф2 — изменение сопротивления второго датчика. Соотношения (17.3) и (17.4) представляют собой систему двух

уравнений с двумя неизвестными N0 и а. Чтобы решение этой системы было единственным, необходимо, чтобы эти уравнения были функционально независимыми. Это накладывает определен­ ные ограничения на вид функций /2 (х) и р (х, а). В частности, функция /2 (х) должна иметь порог чувствительности х*.

Рис. 17.5. Зависимость сопротивления датчиков, установленных на элементе конструкции (а) и на мультипликаторе (б), от числа циклов нагружения

181