Книга S.Gran A Course in Ocean Engineering. Глава Усталость (стр. 5 из 10)
Наконец, распределение Рэлея для размахов напряжений, т.е. d=1 и k=2, дает
где F (_) - нормированный нормальный интеграл, определенный с помощью (2.3A.1).
Усталость вызванная неустановившейся нагрузкой.(xxx) До сих пор мы рассматривали только стационарные (т.е. с постоянной амплитудой), случайные напряжения. До того как отойти от формул усталости замкнутого вида, будет уместно обратить внимание на конструкции, которые испытывают неустановившиеся (т.е. с переменной амплитудой) колебания после импульсной нагрузки. Прибрежный кран, нижняя запись на рис. 4.7.1b, может послужить примером этого явления. Более схематичное изображение дано на рис. 4.7.7. Когда часть груза поднимается краном, конструктивный элемент в кране испытает изменение в статическом уровне напряжений Z. Если статическое напряжение возвращается на начальный уровень, когда нагрузка снята, то элемент испытал один усталостный цикл напряжений с размахом напряжений Z. Использование основной логарифмической S-N кривой (4.7.9) показывает, что эта единичная операция подъема увеличила коэффициент использования h на
Рис. 4.7.7 Последовательность размахов напряжений неустановившейся реакции, полученная с помощью метода дождевого потока для подсчета циклов (the rain-flow cycle counting method).Жирная линия показывает квазистатический цикл напряжений с размахом Z.
Однако, это заниженная оценка, т.к. не учтены динамические явления. Напряженное состояние свидетельствует о том, что за отклонением от начального значения, описанным коэффициентом динамичности Y, следует последовательность неустановившихся циклов, размах напряжений которых последовательно уменьшается на величину определяемую показателем e-aT=e-pl. Здесь, a- коэффициент затухания, T – период колебаний и l=aT/2p- это относительное демпфирование (доля критического демпфирования).
До сих пор не было сказано о том, как подсчитать циклы напряжений. Предполагалось, что цикл может иметь симметричный период, полученный либо из спектральной функции (2.5.78), либо по методу порогового пересечения (threshold crossing procedure), как в главе 2.3.3(i). В случае неустановившихся колебаний это не работает. Однако мы можем использовать более общий метод, известный как метод дождевого потока для подсчета циклов (rain-flow counting method). Более подробное описание посмотрите, например, в работе /2/. Этот метод дает размахи напряжений обозначенные на рис. 4.7.7 как No 1, 2, 3 и т.д.
Подставляя эти циклы напряжений в формулу Палмгрена-Майнера (4.7.10) и проводя суммирование (без предела усталости), мы получим уточненное усталостное значение
Сравнивая с (4.7.34), очевидно, что выражение в фигурных скобках это коэффициент усиления для процесса усталости, вызванного неустановившимися колебаниями. В качестве типичных значений, мы можем подставить m=1 и l=0,025, что дает
Скажем, коэффициент динамичности Y=1, из-за неустановившихся процессов, дает коэффициент усиления усталости равный 7. Это соответствует снижению ресурса, полученного только из изменений статической нагрузки, до 1/7. Подробности даны в работе /9/.
Глава 4.7.4 Естественная дисперсия.
В части 4.3, мы придерживались того, что экстремальное значение в последовательности случайных амплитуд имеет некоторую дисперсию или погрешность, а именно (4.3.19), даже если параметры распределения амплитуд известны точно. Мы назвали это явление естественной дисперсией экстремального значения.
Мы имеем подобный эффект и в усталости. Поскольку коэффициент использования h это сумма вкладов в усталость отдельных циклов напряжений, то эта сумма также обязательно будет случайной, и будет иметь некоторую дисперсию. Если все переменные данного материала и параметры распределения напряжений рассматриваются как заданные, то коэффициент использования h, а также прогнозируемый ресурс, все равно будут иметь погрешность вызванную естественной дисперсией.
Распределение отдельных скачков. Мы можем рассмотреть коэффициент использования как точку движущуюся скачкообразно вдоль координатной оси h. Первоначально, когда конструкция новая, эта точка расположена на h=0. Предполагается, что конструкция изношена или требует ремонта, когда точка проходит через h=1. Коэффициент использования h, на этом отрезке, перемещаетсяскачками
движимый вперед отдельными циклами напряжений. Вообще, Dhj– это возрастание, вызванное j-м циклом напряжений.
Средний период напряжений можно обозначить через T. По истечении времени t=nT конструкция испытывает n циклов напряжений и значение коэффициента использования можно записать как
Длины отдельных скачков Dhпринадлежат одному и тому же статистическому ансамблю и можно предположить, что они имеют одно и то же распределение вероятностей. Поэтому, для удобства, мы опишем длину скачка Dhс помощью случайной величины Dh=xi. Эта переменная связана с размахом напряжений S действительных циклов напряжений по всей S-N кривой. Учитывая, для удобства, основную кривую (4.7.9) это дает
Теперь, S – размах напряжений вызванный действием волн на конструкцию, который подчиняется гамма распределению с плотностью вероятности для больших интервалов времени (4.7.7), т.е. g(d, k, D; S).
Т.к. соотношение (4.7.39) согласовывается с преобразованием энергии (2.6.31), то длина скачка xi также подчиняется гамма распределению с функцией плотности вероятности f(xi)
как было получено в (2.6.33). Математическое ожидание
длины скачка xi получено из (2.6.17) как момент первого порядка 
Соответствующие статистические моменты порядка 2 и 3 около нуля
и
соответственно. Для последующего использования, мы подставили обобщенные скорости U, Vи W, определенные из выражений
В частности, U может интерпретироваться как средний рост коэффициента использования h за один цикл.
Дисперсия
длины скачка xi может быть получена как центральный момент второго порядка (2.4.3), что дает 
Параметр n- это среднеквадратическое отклонение 
относительно математического ожидания длины скачка
, этот параметр можно найти и в (2.4.3). Есть сходство с (2.8.34) по ширине диапазона.Аналогично, центральный момент m3(xi) длины скачка xi может быть получен из (2.4.4) и его можно записать
l-это коэффициент асимметрии длины отдельного скачка. Например, для экспоненциального распределения он равен 2, а для нормального распределения 0.
Распределение вероятностей длины скачка Dh=xi имеет характеристическую функцию f(s), определенную в общем виде в (2.4.8) как
где s, в общем, может быть комплексным параметром. Разложение экспоненты в интеграле в ряд даст
Почленно интегрируя по xi и учитывая выражения (4.7.41) – (4.7.44) получим
Говоря физическим языком, отдельные вклады xi в коэффициент использования, которые вносятся напряжениями вызванными волнами, обычно довольно-таки нерегулярны. Если размах напряжений распределен экспоненциально, что часто бывает в морских конструкциях, то отдельный вклад xi для m=1 имеет среднеквадратическое относительное отклонение n=4,36 и коэффициент асимметрии l=19,6. Следовательно, функция плотности вероятности для отдельных приростов коэффициента использования очень широкая и в значительной степени асимметричная.
Уравнение движения для коэффициента использования.Коэффициент использования h в момент времени t описывают при помощи функции плотности вероятности r(h,t). Соответствующую характеристическую функцию обозначим через F(s,t), где s – такая же переменная, как и в (4.7.48). Ее получают с помощью преобразования плотности вероятности r(h,t)
