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

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

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

Узкополосные процессы с переменным средним значением. Рассмотрим гауссовский узкополосный процесс с переменным средним о (/) (рис. 15.1). Все характеристики этого процесса, за исключением его среднего значения, со временем не изме­ няются. В этом случае за среднее значение цикла нагружения принимается соответствующее значение функции ат (t) — д ((). Соотношение (15.2) в этом случае принимает вид

t ОО

V (0 = яС

j

/ (<Та) d o а d t

(15.3)

N[<Ja, от(/)] ‘

о

о

 

 

Пусть число циклов до разрушения и плотность распределения амплитуд определяются следующими соотношениями:

В этом случае

 

 

/

оо

2

 

v W=

I

S

[а*+ (/)1т

ехр

do•dL

 

 

0 - 1

о

о а - ф т

 

 

 

 

 

 

 

 

 

 

 

 

(15.4)

 

При малых значениях параметра ф можно принять

 

 

 

(Оа+ Ф<7т)

 

,т—1

 

 

 

 

« СТа + ФтОт Оа

 

 

Внутренний интеграл в соотношении (15.4) можно тогда пред­

ставить так:

 

 

 

 

 

 

 

J =

2m/2smr

_(_ i f z)

+

rn\pam2<m- 1)/2sm~ir

y i- , z) ,

где

Г

(v, z) —

неполная

гамма-функция;

z =

( P - l — ф р т )

2s*

^

0-1

fO-lOm

 

 

 

 

 

 

 

 

 

 

 

~2s* s*

Усталостное повреждение за время t

v {t) = тфй12m/2smJ^ [T- +

0 - &<T« Л+

 

. e

*

m + 1

 

+

/n\|:2<m 1)/2sm_1 J amtr [

a - bom (/)] dt , (15.5)

где a =

<r_j/(2s2); b =

фсг.х/в*.

 

 

Рассмотрим частный случай, когда среднее напряжение про­ цесса нагружения изменяется со временем по линейному закону ат (t) = ct, где с — заданный параметр. Вводя новую перемен-

162

Рис. 15.1. Процесс с переменным сред-

Рис. 15.2. Процесс с переменной дис-

ним значением a (t)

Персией

ную г = а bet, приходим при вычислении v (t) к табличным интегралам от неполной гамма-функции [48]:

Ji = j Г (v, г) dz = (z — v) Г (v, г) — zve~z;

J2 = ^ гГ (v, z) dz = (z2— v2 — v) Г (v,

z) — (z + v +

1) zve~z.

Теперь расчет долговечности может

быть доведен до числен­

ного результата.

дисперсией.

Рассмотрим

Случайные процессы с переменной

в качестве примера узкополосные гауссовские процессы с изменя­ ющейся во времени дисперсией и со средним значением о = О (рис. 15.2). В этом случае усталостное повреждение за один цикл нагружения можно вычислить по формуле (14.8), которую за­ пишем в следующем виде:,

v(1) = N o 'x -m,2r

+ 1 , х ),

(15.6)

где х = oii/(2s2).

Пусть дисперсия процесса нагружения является заданной

детерминированной функцией времени s2 =

s®(t). Тогда устало­

стное повреждение за время t

 

 

 

t

 

 

v (*) =

I х~тп W Г [“Г + 1’

х <*>]

(15J)

 

о

 

 

где ? — средний период нагружения.

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

i

Jf,-vr (v, х) dx,

(15.8)

ДГф

 

 

где v = - j 4- 1.

11*

163

Интеграл в соотношении (15.8) табличный. Используя полу­

ченный в работе [48] результат, получаем

 

v (0 = axJN .V -v) [Х^ Г (V* Х) - (Х+ 1) в~Х] Й (1+0°-

(15‘9)

Рассмотрим случай, когда х — х0 ехр (—at). Тогда

 

I

x_Vr (v> x)d * =

 

х0 exp (—a/)

 

 

 

* ) J f c

(15.10)

Аналогично могут быть рассмотрены и другие виды нестационарности.

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

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

+ п2 +

... +

пк смен режимов нагружения, обозначим через Р (пъ

«2* •••» nh/n),

а функцию распределения вероятности события,

что за

все время эксплуатации конструкции произойдет ровно

п смен режимов нагружения, обозначим через Р (п). Тогда без­ условная функция распределения вероятности события того, что

первый режим нагружения будет реализован

раз, второй режим

нагружения будет реализован я2 раза и т. д.,

оо

 

р (П1, я2, ... , пк) = 2J р (til, ti2, ...,

nk/n) Р (п). (15.11)

Л=0

 

Если для каждого режима нагружения вычислена мера уста­ лостного повреждения в единицу времени v* (i = 1,2,..., k), то, вводя в рассмотрение для каждой конкретной конструкции единый расчетный режим нагружения, эквивалентный всему процессу нагружения конструкции за все время ее эксплуатации,

164

получаем расчетное значение усталостного повреждения в еди­ ницу времени

 

<=*

 

 

v = 4

' 2 "

iVi"

(15.12)

 

i= 1

 

 

Так как величины nt (i =

1, 2,

k)

являются случайными,

то и v также будет величиной случайной. Распределение этой величины может быть найдено как распределение функции с за­ данными случайными аргументами. Статистические моменты этой величины

<vm> = 2J

•••

Е

[v(/ix,

п2, .... nh)]m Р (ях, п..........nh),

 

0 (к)

п1=0

 

 

 

 

 

 

(15.13)

где т = 1, 2,

... номер

момента.

На практике в

качестве

распределений Р (пи п2, .... n jn )

и Р (я) удобно использовать полиномиальное распределение и распределение "Пуассона. Эти распределения задаются соответ­ ственно соотношениями

Р (пи па, .... nk/n) =

• • •

а **; (15.14)

 

 

05.15)

где а 1( а 2, ..., а* — вероятности

воспроизведения

первого, вто­

рого и т. д. режимов нагружения,

причем а г + а2 + ... 4- а к —

= 1; я = п0Т — ожидаемое число смен режимов нагружения за все время эксплуатации конструкции; я0 — интенсивность смены режимов нагружения (число смен режимов в единицу времени).

Среднее значение усталостного повреждения при общем числе я

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

 

t=k

 

v - S a , v „

(15.16)

где <xj = «i/л — ожидаемое число реализаций t-го режима нагружения).

Аналогично можно вычислить моменты величины v более высо­ ких порядков. В частности, при большом значении я момент второго порядка величины v

i=k

<v2) « 2] afv2. ‘ г=1

165

§16. Расчеты при случайном нагружении

исложном напряженном состоянии

Постановка задачи. Расчеты статической прочности элементов конструкций при сложном напряженном состоянии сводятся к определению по заданным постоянным значениям компонент напряженного состояния в точке расчетного напряжения по той или иной теории прочности и сопоставлению его с опасным для конструкции напряжением. Если компоненты напряженного со­ стояния со временем изменяются и представляют, например, случайные процессы, возникают дополнительно две новые задачи: 1) определение расчетного напряжения, предопределяющего на­ копление в материале конструкции усталостных повреждений; 2) определение расположения площадки с максимальным значе­ нием накопленного усталостного повреждения.

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

жение принимается октаэдрическое касательное напряжение, то

l/'jT

можно принять x_! = -!у - о_! « 0,7а.!. При использовании

гипотезы наибольших касательных напряжений принимается т_х = = О.бст,!.

Рассмотрим план решения поставленной задачи для случая, когда компоненты тензора напряжений изменяются во времени случайно (несинхронно и несинфазно). Ограничимся случаем плоского напряженного состояния, характеризуемого гауссов­ скими стационарными и стационарно связанными процессами изменения напряжений ах (/), ау (t) и т (t) (рис. 16.1). Эти про­ цессы различаются как по интенсивности воздействия (по диспер­ сиям), так и по частотному составу.

Выберем вначале расчетное напряжение в некоторой пло­ щадке, расположение которой характеризуется направляющими косинусами I = cos (х, v), т — cos (у, v), п = cos (г, v), где v — ось, перпендикулярная к рассматриваемой площадке (рис. 16.2).

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

av = axl2 -f оут2 + 2m/т;

(16.1)

х%= (о*/ + тт)2 + (т/ + оут)2— о2-

166

бу
X
t
Рис. 16.1. Плоское напряженное состоя­ ние с компонентами напряжений, изме­ няющимися по случайным законам

б,

z

X

Рис. 16.2. Напряжения в площад­ ке общего вида, определяемой нор­ малью v

Будем считать, что сопро­ тивление материала кон­ струкции усталости обуслов­

ливается наличием в площадках как касательных, так и нормаль­ ных напряжений. Тогда расчетное напряжение в некоторой рас­ сматриваемой площадке

ар = (2 — X) <yv + (А, — 1) ту,

( 16.2)

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

Коэффициент

 

X = OL

 

где o_i — предел

выносливости по нормальным напряжениям;

т_х — предел выносливости по касательным напряжениям.

Для идеально пластичных материалов X =

2, а для идеально

хрупких X = 1.

В первом случае прочность

материала опреде­

ляется только касательными напряжениями, а во втором — только нормальными.

Вероятностные характеристики расчетного случайного про­ цесса нагружения, определяемого соотношением (16.2), пол­ ностью определяются по заданным вероятностным характери­ стикам исходных процессов ах (t), о„ (t) и т (/). В частности, для расчетного процесса <тр (t) могут быть определены плотность распределения амплитуд напряжений f (о, /, т, п) и частота со (/, т, п) как функции направляющих косинусов, определяемых расположением рассматриваемой площадки. По формулам, при­ веденным в §§ 13—15, может быть вычислено накопленное к не­ которому заданному моменту времени усталостное повреждение v (/, т, п) как функция параметров I, т и п . Исследуя эту функцию как максимум, определяем расположение опасной площадки. Долговечность конструкции теперь может быть вычислена, на­

167

Рис. 16.3. Частные случаи площадок бесконечно малого эле­ мента

пример, по формуле типа (14.11), где в качестве расчетного пре­ дела выносливости принимается

о_1)Р — (2 — X) а_х + (Я — 1) т_х.

(16.3)

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

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

Нормальные и касательные напряжения в этих площадках можно вычислить по следующим формулам:

для

площадки

1

 

 

оа — °х °и

+ .?* ^ av . cos 2а — т sin 2а;

 

 

 

(16.4)

 

 

тв =

sin 2а -f 1 cos 2а;

для

площадки

2

 

<*3= (1 — cos 2р);

(16.5)

тр = — Ц- sin 20;

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

<*v = - у - (1 + cos 2у);

(16.6)

т, = - у - sin 2у.

168

Подставив

(16.4),

(16.5) и (16.6) в фор-

у

 

мулу (16.2), получим

соответствующие вы­

 

 

ражения для определения расчетного напря­

 

 

жения.

Во всех рассмотренных

площадках

А

 

это напряжение

представляется

в виде сле­

 

 

дующего линейного соотношения:

 

 

\ z

 

 

Op = Аах -|- By -f- Ст,

(16.7)

 

 

 

 

где коэффициенты А,

В а С определяются

формаций

р^ Т а де­

по формулам:

 

 

 

 

 

 

 

 

для

площадки 1

 

 

 

 

 

 

 

А =

2Т Х + ——ip—cos 2а +

 

1 sin 2а;

 

 

 

 

2

2

 

 

 

 

 

 

 

 

2

X

2 К

cos 2 а -----—

1 sin 2а;

(16.8)

 

 

 

 

2

2

--------

2

 

 

 

С =

— 1) cos 2а — (2 — к) sin 2а;

 

для

площадки 2

А = С = 0;

 

 

 

 

 

 

2 —к

 

 

(16.9)

 

В =

 

- ^ c o s 2Р---- ^Z-L_sin 20;

 

 

 

для

площадки 3

+ 2 — к

 

 

 

 

 

 

А = -

 

^

cos -j— - я—- sin 2у;

(16.10)

 

 

 

 

 

В =

С = 0.

 

 

 

 

 

 

 

 

 

 

Исходной информацией для расчёта обычно являются осциллографические записи деформаций, соответствующие трем выбран­ ным направлениям. При этом, например, используется прямо­ угольная розетка тензодатчиков (рис. 16.4). Тогда в результате испытаний регистрируются три случайных процесса вх (/), еу (/) и ву (t). В этом случае удобно перейти к псевдонапряжениям

дх (0 = Евх (*),

Ьу (0 = Е&у (Г),

by (t) = Еву (t).

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

<3х = Т=7ТГ (°Х + рду);

= i l p - К

+ ^* ):

(16.11)

Т = 2(1 + Ц)

6v)'

 

169

Подставив соотношения (16.11) в (16.7), получим расчетное напряжение

о = Агдх + Вхду + Сгду,

(16.12)

где

лA -f- Вр________ С

 

1 - р ’

2(1 + |i)

 

D

Ар + В_________С

 

(16.13)

 

1 — I**

2(1 +Ю

 

 

В формуле (16.12) и далее индекс р расчетного напряжения опущен.

Трехмерный случайный процесс (d* (t), by (t), b7 (()}, реги­ стрируемый в эксперименте, полностью характеризуется следу­ ющей матрицей корреляционных и взаимных корреляционных функций:

Г К х (т)

К ху (т)

К Ху (т) Л

 

м = I КуХ(т)

Ку (т)

Куу (т)

I ,

(16.14)

L КуХ(т)

Куу(х)

к у (т)

]

 

где Кх (т), К у (т), Ку

(т) — корреляционные функции

процессов

Ьх (t), Ьу (t) и by (/);

КХу (т), .... КуУ (т) — взаимные

корреля­

ционные функции этих процессов.

 

Корреляционную функцию расчетного процесса нагружения, определяемого соотношением (16.12), можно записать в следу­

ющем виде:

 

 

Ко (т) = А2гКх (т) + В\Кд (т) +

С\Ку (т) + 2Л,С, к ф )

+

+ 2А1В1Кху (т) +

2B1C1Kvy (т).

(16.15)

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

So — A\sx -f- B\sy -(- C\sy -J- 2A\BiKXy -I- 2A\C\Kxy “b 2B\C\Kyy\

So A \S x

-(-

B{S$ -f- C JS^I -f- 2 A i B i K i $

2 A i C i K t y -f-

2 B \ C i K g ^ '

 

 

 

 

(16.16)

So A 2sl

 

B\s2y -f- CjS^ -f- 2A\B\Kij -)- 2A\C\Kii ~b 2B\C\Kg

где si, si,

si,

Kxy, Kxy, •••. Kgy— дисперсии и взаимные корреля­

ционные моменты процессов bx (t), by (t) и ач (т) и

их первых

двух производных.

 

нагружения

Средние

значения частот расчетного процесса

по нулям и экстремумам соответственно равны

 

©О— Sd/Sa\

(16.17)

©э ~~~ S{j/S<3r.'

170

Подставив соотношения (16.16) и (16.17) в формулу типа (14.8), получим усталостное повреждение v как функцию углов а, р и у в соответствии с тем, какая из трех площадок рассматривается (см. рис. 16.3). Расположение опасной площадки и расчетное усталостное повреждение определяются из условия максимума накопленного в ней усталостного повреждения

v =

max {v(a),

v(P),

v(y)},

(16.18)

 

а. Р. V

 

 

 

 

 

где углы а, р и у определяются

из уравнений

 

 

dv

_

dv _

dv

__Л

(16.19)

 

da

~~

dfi

~~

dy

~

 

 

В первом приближении расположение опасной площадки

можно определить из условия максимума дисперсии:

 

ds

_

dsa

_

dsa

 

(16.20)

da

~

dp

~

dy

 

 

 

Подставив в уравнения (16.20) первое из выражений (16.16), придем к следующему соотношению для определения расположе­ ния опасной площадки:

A\A\sx +

B\B\s2y 4* CiCiSy -(- (А\В\

+ А\В\) Кху +

+ (А &

+ ХгСг) К*v + (В & +

В & ) К „ = 0 , (16.21)

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

Выпишем некоторые выражения для определения коэффи­ циентов, входящих в соотношение (16.21). Так, для площадок, расположение которых определяется углом а,

A l = 2(1 ц ) + 2 (1 + (1) C0S 2 а + 2(1 + ц ) S ln 2rt;

B l = 2 (1 —ц)

2 (1 + ц)

cos2oc+ 2 ( 1 + ц ) sm 2а: (!6.22)

~

X — 1

о

2 А, , о

С, = —7-:— cos 2 а ----i—;— sin 2а.

1

1+!*■

 

1+И.

Аналогично можно выписать и все другие коэффициенты. Рассмотрим более подробно'случай, когда X — 2, т. е. случай,

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

А \ -

V 2

sin

2(1 + р)

V T

(16.23)

2 ( 1 + Ц )

 

т-п -с°$2а .

 

171