ОБОСНОВАНИЕ НОРМАТИВНОГО ИЗНОСА ЭЛЕМЕНТОВ КОНСТРУКЦИЙ ПЛАВУЧЕГО ДОКА В ТРЕБОВАНИЯХ К ОБЩЕЙ ПРОЧНОСТИ И УСТОЙЧИВОСТИ ПРАВИЛ РМРС НА ОСНОВЕ ВЕРОЯТНОСТНОГО МОДЕЛИРОВАНИЯ

Полный текст

Введение

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

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

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

Требования Правил РМРС к общей прочности и устойчивости продольных связей

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

Добавка на коррозионный износ  – го элемента конструкции ПД при этом определяется по формуле

                                                                      (1)

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

 – нормативный расчетный срок службы дока, принимаемый равным 50 лет, если нет других условий.

Требования к общей продольной прочности корпуса плавучих доков Правил РМРС [1] подробно были рассмотрены в [2, 3].

В общем виде, согласно [1], требуемый момент сопротивления поперечного сечения корпуса дока на начало срока службы, , см³, определяется по формуле

                                                                          (2)

где         – требуемый момент сопротивления поперечного сечения к концу срока службы дока, см³;

 – коэффициент, определяющий добавку к моменту сопротивления на износ элементов конструкций, находится по формуле, учитывающей износ одновременно всех элементов эквивалентного бруса;

 – наибольший расчетный изгибающий момент, кН·м;

 – изгибающий момент, вызывающий прогиб дока, кН·м;

 – изгибающий момент, вызывающий перегиб дока, кН·м;

 – коэффициент допускаемых напряжений от общего продольного изгиба, который принимается равным 1,0;

 – нормативный предел текучести по нормальным напряжениям, МПа.

Условие устойчивости согласно [1] имеет вид

                                                                         (3)

где         – коэффициент запаса устойчивости;

 – расчетные сжимающие напряжения  – го элемента, определяемые для поперечного сечения дока к концу расчетного срока службы, МПа;

 – критические напряжения для  – го элемента ЭБ определяемые с учетом 100% износа, МПа;

 – расчетный изгибающий момент, вызывающий сжатие  – го элемента ЭБ, кН·м;

 – фактический момент инерции поперечного сечения ЭБ, определяемый к концу срока службы дока, см⁴;

 – отстояние  – го элемента ЭБ от нейтральной оси, м.

Расчетные сжимающие и критические напряжения в условии устойчивости (3) также определяются с учетом 100%-го износа элементов эквивалентного бруса.

Из условий общей прочности (2) и устойчивости (3) следует, что они могут приводить к необоснованному увеличению размеров продольных связей корпуса дока. Это подтверждается опытом практического использования требований Правил.

Методика и алгоритм определения среднестатистического износа корпуса дока по состоянию на конец срока службы

Для корректировки требований к конструкциям плавучих доков Правил РМРС в плане смягчения требований по общей прочности и устойчивости при общем продольном изгибе предлагается методика учета неравномерности износа элементов конструкции дока в процессе эксплуатации. Такой же подход принят в Правилах МАКО, регламентирующих требования к конструкции корпуса нефтеналивных и навалочных судов. Характеристики поперечного сечения корпуса из условия общей прочности и расчетные сжимающие напряжения в условиях устойчивости листовых и балочных элементов конструкций таких судов регламентируются по состоянию на середину срока службы, т.е. при 50%-ном износе. Учитывая, что расчетный срок службы судна в Правилах МАКО принимается 25 лет, а также существенно большие значения скоростей среднегодового уменьшения элементов конструкции нефтеналивных и навалочных судов по сравнению с доковыми конструкциями, нормативное состояние конструкций дока по износу может быть иным.

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

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

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

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

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

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

1. Построение расчетных зависимостей для определения добавок на коррозионный износ к требуемым характеристикам поперечного сечения корпуса дока.
2. Моделирование износа конструкций с использованием датчиков случайных чисел.
3. Построение статистических законов (гистограмм) распределения добавок на коррозионный износ.
4. Аппроксимация статистического распределения теоретическим законом. Проверка согласования теоретического и статистического распределений с использованием методов проверки правдоподобия статистических гипотез.
5. Обоснование нормативного износа конструкций плавучего дока в требованиях к общей прочности и устойчивости.

Подготовка расчетных зависимостей для определения добавок на коррозионный износ к характеристикам поперечного сечения корпуса дока

В качестве примера рассмотрим поперечное сечение монолитного плавучего дока грузоподъемностью 12 000 т, схема эквивалентного бруса которого соответствует рисунку 1.

Рисунок 1 – Схема поперечного сечения эквивалентного бруса монолитного дока грузоподъемностью 12 000 т

Добавки, учитывающие износ элементов для каждого j-го эксперимента можно определить по известным аналитическим зависимостям:

- добавка к площади –                          , см²;

- добавка к моменту инерции –          , м²см²;

- добавка к моменту сопротивления на уровне верхнего пояска (топ-палубы) ЭБ –                                                                                     , см³;

- добавка к моменту сопротивления на уровне нижнего пояска (днища понтона) ЭБ –                                                                     , см³.

При этом             – количество рассматриваемых вариантов износа ЭБ (общее число экспериментов в выборке или объем выборки);

 – добавка к площади поперечного сечения  – го элемента ЭБ, учитывающая его износ, см²;

 – добавка к собственному моменту инерции  – го элемента ЭБ, учитывающая его износ, м²см²;

 – отстояние центра тяжести  от оси сравнения (ось сравнения принята на уровне основной плоскости ПД), м;

 – расчетная высота плавучего дока, м.

Согласно положениям, описанным в [5], при малом объеме выборки  распределение существенно отличается от функции распределения Колмогорова, но уже при 0 ошибка при вычислении вероятности «согласия» оказывается малой и при дальнейшем увеличении объема выборки практически не изменяется. По критерию согласия Пирсона, по рекомендациям [6] объем выборки должен быть не менее 100, согласно этому, в настоящей работе принято .

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

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

                        (4)

где         – максимально допустимая Правилами РМРС скорость коррозионного износа, которая соответствует табличным значениям, представленным в п. 3.12.1.4.2 [1], мм/год.

Функция СЛУЧМЕЖДУ(минимальный порог; максимальный порог), реализуемая в Microsoft Excel генерирует случайное целое число между значениями минимального и максимального порога включительно.

Используя сгенерированную скорость коррозионного износа согласно (4), можно определить добавку на износ  – го элемента, , мм, по (1).

Рисунок 2 – Изменение добавки к площади сечения ЭБ в зависимости от степени износа дока

Гистограммы распределения значений добавок на износ характеристик эквивалентного бруса построены в результате анализа  различных вариантов износа ЭБ. Для каждого  – го эксперимента определяются характеристики: ; ; ;  и износ в процентах, с учетом, что  – соответствует полному износу всех элементов ЭБ дока по состоянию на конец срока службы (т.е. при , мм/год), 100% износ;  – соответствует началу срока службы дока (т.е. при , мм/год), износ составляет 0%, см. рисунок 2. Уровень износа для характеристик: ; ;  определяется аналогично .

Построение статистических законов (гистограмм) распределения добавок на коррозионный износ

Получена выборка, состоящая из  вариантов различных износов поперечного сечения ПД (объем выборки), для которых определены фактические геометрические характеристики ЭБ и соответствующий им уровень износа в процентах. Для полученной выборки выполняется анализ частот попадания в тот или иной разряд гистограммы (с шагом 5% от  или  года), с использованием функции ЧАСТОТА(массив_данных; массив_интервалов), реализуемая в Microsoft Excel. В результате принимается наибольшее полученное значение износа, по которому строится результирующая гистограмма.

Гистограмма распределения фактического износа поперечного сечения плавучего дока будет построена с принятым в настоящей работе шагом равным , результаты представлены на рисунок 3.

Рисунок 3 – Гистограмма распределения износа дока в %

Важно, что функция СЛУЧМЕЖДУ( ) обновляет сгенерированное случайное значение при каждом пересчете (клавиша F9), таким образом можно рассмотреть результаты анализа для нескольких выборок, что позволит более достоверно оценить полученный результат.

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

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

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

- критерия  (критерия Пирсона);

- критерия Колмогорова.

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

Ниже представлены основные математические зависимости, которые использовались для построения нормального закона распределения, и расчетные формулы для определения характеристик критериев согласия [7, 8, 9]:

- Нормальный закон распределения:

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

                                                  (5)

математическое ожидание

                                                             (6)

                                                           (7)

дисперсия случайной величины

                                                (8)

                                                                      (9)

среднее квадратическое отклонение

                                                                               (10)

функция распределения

           (11)

Критерий согласия .

Находим значение критерия по формуле

                                                             (12)

где         – эмпирические частоты, определенные с помощью функция ЧАСТОТА( ) для полученной выборки;

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

= (НОРМРАСП( ; ; ; ИСТИНА) – НОРМРАСП( ; ; ; ИСТИНА))* ;

 – число интервалов настоящей выборки.

Определяем критическое значение критерия  на уровне значимости  и с  степенью свободы, в среде Excel для определения  может быть использована функция

= ХИ2ОБР( ; ).

Если , гипотезу о распределение можно принять, иначе стоит отвергнуть.

Критерий Колмогорова А.Н.

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

Определяются теоретические частоты  (как для критерия Пирсона) и соответствующие им накопленные частоты .

Вычисляется разность накопленных частот для всех интервалов , максимальное по модулю отклонение теоретических от эмпирических частот обозначено, как .

При  гипотеза о распределении может быть принята, иначе – отклоняется, при этом  определяется по формуле

,               при                                  (13)

Рисунок 4 – Гистограмма распределения эмпирических и теоретических частот износа монолитного дока грузоподъёмностью 12 000 т

Рисунок 5 – Гистограмма распределения накопленных эмпирических и теоретических частот износа монолитного дока грузоподъёмностью 12 000 т

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

Рисунок 6 – Нормативный износ поперечного сечения монолитного дока грузоподъёмностью 12 000 т для общей продольной прочности

Обоснование нормативного износа конструкций плавучего дока в требованиях к общей прочности и устойчивости

Таким образом, в (2) требуемый момент сопротивления поперечного сечения к концу срока службы дока ( ) будет соответствовать не 100%-ному износу ЭБ, а экспериментально оцененному в настоящей работе, с ошибкой в безопасную сторону можно принять 70%. Такой подход повлечет за собой пересмотр уровня допускаемых напряжений.

В условии устойчивости (3), предлагается, что расчетные сжимающие напряжения ( ) должны быть определены с учетом износа поперечного сечения дока соответствующего 70%, при этом ( ) эйлеровы напряжения  - го элемента ЭБ должны определяться при условии 100%-го износа рассматриваемого элемента.

Настоящие изменения позволят существенно смягчить действующие требования Правил РМРС к общей прочности и устойчивости при продольном изгибе дока.

Заключение

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

- математическое ожидание принимает значения в диапазоне: от 50 до 54%;

- стандартное отклонение принимает значения в диапазоне: 4 – 5%.

Тогда с вероятностью 99,7% нормативный износ эквивалентного бруса дока, участвующего в общем продольном изгибе, на конец срока службы может быть принят 65 – 70%.

Об авторах

Дианна Тимофеевна Иванова

Федеральное государственное бюджетное образовательное учреждение высшего образования «Санкт-Петербургский государственный морской технический университет»; Акционерное общество "Центральное морское конструкторское бюро "Алмаз"

Автор, ответственный за переписку.
Email: D.T.Ivanova@smtu.ru
SPIN-код: 1314-7682
Scopus Author ID: 1271952

старший преподаватель кафедры конструкции и технической эксплуатации судов; инженер-конструктор 3 кат.

Россия, 190121, Северо-Западный Федеральный округ, Санкт-Петербург, улица Лоцманская, дом 3; 196128, город Санкт-Петербург, Варшавская ул., д.50

Владимир Николаевич Тряскин

Федеральное государственное бюджетное образовательное учреждение высшего образования «Санкт-Петербургский государственный морской технический университет»

Email: vladimir.tryaskin@smtu.ru
SPIN-код: 7663-8210
Scopus Author ID: 542496

Профессор, доктор технических наук, заведующий кафедры конструкции и технической эксплуатации судов

Россия, 190121, Северо-Западный Федеральный округ, Санкт-Петербург, улица Лоцманская, дом 3

Список литературы

Дополнительные файлы