Обоснование нормативного износа элементов конструкций плавучего дока в требованиях к общей прочности и устойчивости Правил РМРС на основе вероятностного моделирования
- Авторы: Иванова Д.Т.1, Тряскин В.Н.1
-
Учреждения:
- Санкт-Петербургский государственный морской технический университет
- Страницы: 181-188
- Раздел: Машиностроение
- Статья получена: 13.02.2025
- Статья одобрена: 03.04.2025
- Статья опубликована: 30.06.2025
- URL: https://vestnik.nvsu.ru/2414-1437/article/view/655836
- DOI: https://doi.org/10.52899/24141437_2025_02_181
- EDN: https://elibrary.ru/VOKWZV
- ID: 655836
Цитировать
Полный текст
Аннотация
Обоснование. Согласно действующим Правилам РМРС, Часть II «Корпус», п. 3.12 Плавучие доки, допускаемые напряжения при общем продольном изгибе дока и расчетные сжимающие напряжения в требованиях к устойчивости продольных связей нормируются по состоянию на конец расчетного срока службы дока и определяются в предположении, что все элементы эквивалентного бруса (ЭБ) имеют 100% нормативный износ. Между тем, такая ситуация в реальных условиях маловероятна.
Цель. В рамках совершенствования требований Правил РМРС к общей продольной прочности и устойчивости при общем изгибе дока необходимо оценить допускаемый уровень нормативного износа элементов корпусных конструкций плавучих доков в расчетных зависимостях, регламентирующих общую прочность и устойчивость, с целью обоснованного смягчения этих требований.
Методы. Для достижения поставленной задачи предполагается применить специальную методику оценки среднего износа продольных связей корпуса дока по состоянию на конец срока службы, основанную на случайном моделировании уменьшения размеров каждого элемента ЭБ в связи с коррозионным износом с использованием датчиков случайных чисел. Указанная методика позволит построить гистограммы распределения значений добавок на износ для характеристик эквивалентного бруса (площади сечения, минимального момента сопротивления, центрального момента инерции), аппроксимировать их неким теоретическим законом, найти математическое ожидание и значения этих характеристик с заданной обеспеченностью.
Выводы. Представленные в настоящей работе оценки показывают, что добавки на коррозионный износ следует принимать на уровне: 65-70% от величин, регламентируемых Правилами на конец нормативного срока службы дока. Это приведет к существенному смягчению требований к общей продольной прочности и устойчивости при общем продольном изгибе корпуса дока.
Полный текст
ВВЕДЕНИЕ
В процессе эксплуатации плавучего дока (ПД) элементы корпусных конструкций (листы, балки набора) подвергаются коррозионному износу. Как известно из практического опыта эксплуатации плавучих доков и других морских сооружений, износ элементов конструкций, работающих в одинаковых условиях, довольно редко принимает одинаковые значения.
Возникает проблема определения фактических характеристик поперечного сечения эквивалентного бруса на конец срока службы с учетом неравномерного износа его элементов.
С целью решения такой задачи предлагается произвести моделирование износа элементов конструкций, участвующих в общем продольном изгибе дока, с применением датчика случайных чисел.
ТРЕБОВАНИЯ ПРАВИЛ РМРС К ОБЩЕЙ ПРОЧНОСТИ И УСТОЙЧИВОСТИ ПРОДОЛЬНЫХ СВЯЗЕЙ
При разработке действующих требований Правил к общей продольной прочности плавучего дока нормативные допускаемые напряжения определялись в предположении, что все элементы эквивалентного бруса имеют 100% износ.
Добавка на коррозионный износ i-го элемента конструкции ПД при этом определяется по формуле
, (1)
где uzon_i = kzonui — нормативная скорость среднегодового уменьшения толщины i-ой связи эквивалентного бруса, определяемая с учетом зонального коэффициента kzon; T — нормативный расчетный срок службы дока, принимаемый равным 50 лет, если нет других условий.
Требования к общей продольной прочности корпуса плавучих доков Правил РМРС [1] подробно были рассмотрены в [2, 3].
В общем виде, согласно [1], требуемый момент сопротивления поперечного сечения корпуса дока на начало срока службы, W, см3, определяется по формуле
, (2)
где — требуемый момент сопротивления поперечного сечения к концу срока службы дока, см3; ωk — коэффициент, определяющий добавку к моменту сопротивления на износ элементов конструкций, находится по формуле, учитывающей износ одновременно всех элементов эквивалентного бруса; — наибольший расчетный изгибающий момент, кН·м; — изгибающий момент, вызывающий прогиб дока, кН·м; — изгибающий момент, вызывающий перегиб дока, кН·м; kσ — коэффициент допускаемых напряжений от общего продольного изгиба, который принимается равным 1,0; σn — нормативный предел текучести по нормальным напряжениям, МПа.
Условие устойчивости согласно [1] имеет вид
kσсi ≤ σcri , (3)
где k — коэффициент запаса устойчивости; — расчетные сжимающие напряжения i-го элемента, определяемые для поперечного сечения дока к концу расчетного срока службы, МПа; — критические напряжения для i-го элемента ЭБ определяемые с учетом 100% износа, МПа; M — расчетный изгибающий момент, вызывающий сжатие i-го элемента ЭБ, кН·м; — фактический момент инерции поперечного сечения ЭБ, определяемый к концу срока службы дока, см4; zi — отстояние i-го элемента ЭБ от нейтральной оси, м.
Расчетные сжимающие и критические напряжения в условии устойчивости (3) также определяются с учетом 100%-го износа элементов эквивалентного бруса.
Из условий общей прочности (2) и устойчивости (3) следует, что они могут приводить к необоснованному увеличению размеров продольных связей корпуса дока. Это подтверждается опытом практического использования требований Правил РМРС.
МЕТОДИКА И АЛГОРИТМ ОПРЕДЕЛЕНИЯ СРЕДНЕСТАТИСТИЧЕСКОГО ИЗНОСА КОРПУСА ДОКА ПО СОСТОЯНИЮ НА КОНЕЦ СРОКА СЛУЖБЫ
Для корректировки требований к конструкциям плавучих доков Правил РМРС в плане смягчения требований по общей прочности и устойчивости при общем продольном изгибе предлагается методика учета неравномерности износа элементов конструкции дока в процессе эксплуатации. Такой же подход принят в Правилах МАКО, регламентирующих требования к конструкции корпуса нефтеналивных и навалочных судов. Характеристики поперечного сечения корпуса из условия общей прочности и расчетные сжимающие напряжения в условиях устойчивости листовых и балочных элементов конструкций таких судов регламентируются по состоянию на середину срока службы, т.е. при 50%-ном износе. Учитывая, что расчетный срок службы судна в Правилах МАКО принимается 25 лет, а также существенно большие значения скоростей среднегодового уменьшения элементов конструкции нефтеналивных и навалочных судов по сравнению с доковыми конструкциями, нормативное состояние конструкций дока по износу может быть иным.
Предлагается методика оценки износа конструкций дока к концу нормативного срока службы с заданной обеспеченностью, основанная на вероятностном моделировании. Предполагается, что в процессе эксплуатации дока происходит коррозионный износ конструкций, обеспечивающих его общую прочность. Скорости износа для каждого конструктивного элемента могут изменяться от нуля до максимального значения, регламентируемого Правилами. Для каждого i-го элемента ЭБ среднее значение скорости износа в течение расчетного срока службы T и, соответственно, добавки на износ , определяемые по формуле (1), генерируется с использованием датчика случайных чисел для равномерного распределения.
По результатам моделирования строятся статистические законы распределения добавок на износ или их относительных величин в долях или в % от максимального значения.
Подбирается теоретическое распределение наилучшим образом аппроксимирующее статистическое распределение (гистограмму). Для этого используются известные методы и критерии проверки статистических гипотез (критерии χ2 — Пирсона, критерий Колмогорова).
Далее определяется нормативный износ конструкций по состоянию на конец срока службы с заданной обеспеченностью (например, с использованием правила 3-х стандартов).
Похожая методика рассмотрена в работе [4] при обосновании так называемого эквивалентного срока службы судна. В этой работе с использованием вероятностного подхода моделировались результаты замеров остаточных толщин конструкций судов в эксплуатации.
Ниже на основе предложенной методики рассмотрен алгоритм решения задачи обоснования нормативного износа конструкций корпуса дока по состоянию на конец срока службы, который предполагает решение следующих задач:
- Построение расчетных зависимостей для определения добавок на коррозионный износ к требуемым характеристикам поперечного сечения корпуса дока.
- Моделирование износа конструкций с использованием датчиков случайных чисел.
- Построение статистических законов (гистограмм) распределения добавок на коррозионный износ.
- Аппроксимация статистического распределения теоретическим законом. Проверка согласования теоретического и статистического распределений с использованием методов проверки правдоподобия статистических гипотез.
- Обоснование нормативного износа конструкций плавучего дока в требованиях к общей прочности и устойчивости.
ПОДГОТОВКА РАСЧЕТНЫХ ЗАВИСИМОСТЕЙ ДЛЯ ОПРЕДЕЛЕНИЯ ДОБАВОК НА КОРРОЗИОННЫЙ ИЗНОС К ХАРАКТЕРИСТИКАМ ПОПЕРЕЧНОГО СЕЧЕНИЯ КОРПУСА ДОКА
В качестве примера рассмотрим поперечное сечение монолитного плавучего дока грузоподъемностью 12 000 т, схема эквивалентного бруса которого соответствует рис. 1.
Рис. 1. Схема поперечного сечения эквивалентного бруса монолитного дока грузоподъемностью 12 000 т.
Fig. 1. Cross-sectional girder diagram of a cast-in-place dock with a carrying capacity of 12,000 tons.
Добавки, учитывающие износ элементов для каждого j-го эксперимента, можно определить по известным аналитическим зависимостям:
- добавка к площади — , см2;
- добавка к моменту инерции —
, м2см2;
- добавка к моменту сопротивления на уровне верхнего пояска (топ-палубы) ЭБ —
, см3;
- добавка к моменту сопротивления на уровне нижнего пояска (днища понтона) ЭБ —
, см3.
При этом:
j = 1...n — количество рассматриваемых вариантов износа ЭБ (общее число экспериментов в выборке или объем выборки);
— добавка к площади поперечного сечения i-го элемента ЭБ, учитывающая его износ, см2;
— добавка к собственному моменту инерции i-го элемента ЭБ, учитывающая его износ, м2см2;
zi — отстояние центра тяжести от оси сравнения (ось сравнения принята на уровне основной плоскости ПД), м;
D — расчетная высота плавучего дока, м.
Согласно положениям, описанным в [5], при малом объеме выборки n распределение существенно отличается от функции распределения Колмогорова, но уже при n ≥ 15...20 ошибка при вычислении вероятности «согласия» оказывается малой и при дальнейшем увеличении объема выборки практически не изменяется. По критерию согласия Пирсона, по рекомендациям [6] объем выборки должен быть не менее 100, согласно этому, в настоящей работе принято n = 100.
МОДЕЛИРОВАНИЕ ИЗНОСА КОНСТРУКЦИЙ С ИСПОЛЬЗОВАНИЕМ ДАТЧИКОВ СЛУЧАЙНЫХ ЧИСЕЛ
В работе моделирование износа выполняется в процессе реализации в ячейке вычислительной процедуры следующего вида
ui = СЛУЧМЕЖДУ (4)
где — максимально допустимая Правилами РМРС скорость коррозионного износа, которая соответствует табличным значениям, представленным в п. 3.12.1.4.2 [1], мм/год.
Функция СЛУЧМЕЖДУ(минимальный порог; максимальный порог), реализуемая в Microsoft Excel, генерирует случайное целое число между значениями минимального и максимального порога включительно.
Используя сгенерированную скорость коррозионного износа согласно (4), можно определить добавку на износ i-го элемента, , мм, по (1).
Гистограммы распределения значений добавок на износ характеристик эквивалентного бруса построены в результате анализа n различных вариантов износа ЭБ. Для каждого j-го эксперимента определяются характеристики: ; ; ; и износ в процентах, с учетом, что — соответствует полному износу всех элементов ЭБ дока по состоянию на конец срока службы (т.е. при , мм/год), 100% износ; — соответствует началу срока службы дока (т.е. при ui = 0, мм/год), износ составляет 0%, см. рисунок 2. Уровень износа для характеристик: ; ; определяется аналогично .
Рис. 2. Изменение добавки к площади сечения ЭБ в зависимости от степени износа дока.
Fig. 2. Change in the cross-sectional area allowance of the girder depending on the degree of wear of the dock.
ПОСТРОЕНИЕ СТАТИСТИЧЕСКИХ ЗАКОНОВ (ГИСТОГРАММ) РАСПРЕДЕЛЕНИЯ ДОБАВОК НА КОРРОЗИОННЫЙ ИЗНОС
Получена выборка, состоящая из n вариантов различных износов поперечного сечения ПД (объем выборки), для которых определены фактические геометрические характеристики ЭБ и соответствующий им уровень износа в процентах. Для полученной выборки выполняется анализ частот попадания в тот или иной разряд гистограммы (с шагом 5% от T или года), с использованием функции ЧАСТОТА(массив_данных; массив_интервалов), реализуемая в Microsoft Excel. В результате принимается наибольшее полученное значение износа, по которому строится результирующая гистограмма.
Гистограмма распределения фактического износа поперечного сечения плавучего дока будет построена с принятым в настоящей работе шагом равным , результаты представлены на рис. 3.
Рис. 3. Гистограмма распределения износа дока в %.
Fig. 3. Wear distribution bar graph (%).
Важно, что функция СЛУЧМЕЖДУ( ) обновляет сгенерированное случайное значение при каждом пересчете (клавиша F9), таким образом можно рассмотреть результаты анализа для нескольких выборок, что позволит более достоверно оценить полученный результат.
АППРОКСИМАЦИЯ СТАТИСТИЧЕСКОГО РАСПРЕДЕЛЕНИЯ ТЕОРЕТИЧЕСКИМ ЗАКОНОМ
Аппроксимация статистического закона распределения некоторым теоретическим законом позволяет получить оценки нормативного износа дока с заданной обеспеченностью (вероятностью превышения).
Построение теоретических законов распределения износа конструкций корпуса ПД заключается в том, чтобы подобрать теоретическую плавную кривую распределения, описывающую данное статистическое распределение. Вопрос о согласовании теоретического и статистического распределения решается методами проверки правдоподобия гипотез с использованием двух критериев согласия:
- критерия χ2 (критерия Пирсона);
- критерия Колмогорова.
В настоящей работе анализируется соответствие статистического распределения нормальному закону (закон Гаусса), гамма-распределению, распределению Вейбулла. Показано, что эти теоретические распределения приводят практически к одинаковым результатам, поэтому принято решение использовать нормальный закон.
Ниже представлены основные математические зависимости, которые использовались для построения нормального закона распределения, и расчетные формулы для определения характеристик критериев согласия [7–9].
Нормальный закон распределения
Плотность вероятности
. (5)
Математическое ожидание
, (6)
. (7)
Дисперсия случайной величины
, (8)
. (9)
Среднее квадратическое отклонение
. (10)
Функция распределения
. (11)
Критерий согласия χ2.
Находим значение критерия по формуле
χ2 (12)
где nz — эмпирические частоты, определенные с помощью функции ЧАСТОТА(массив_данных; массив_интервалов) для полученной выборки;
— теоретические частоты, которые определяются как произведение плотности вероятности на участке от xk до xk+1 на общее число значений в рассматриваемой выборке n (n = 100), используя функции Excel можно записать для нормального распределения
= (НОРМРАСП(; m; σ; ИСТИНА) – – НОРМРАСП(; ; σ; ИСТИНА))*n;
k — число интервалов настоящей выборки.
Определяем критическое значение критерия на уровне значимости α и с степенью свободы, в среде Excel для определения может быть использована функция ХИ2ОБР(вероятность; степени_свободы), которую можно записать в виде
= ХИ2ОБР(α = 0,05; ).
Если , гипотезу о распределение можно принять, иначе стоит отвергнуть.
Критерий А.Н. Колмогорова
Имея эмпирические частоты для каждого рассматриваемого интервала nj ( года), можем определить накопленные эмпирические частоты , например, для второго разряда накопленная частота будет равна сумме частот для первого и второго интервалов и т.д.
Определяются теоретические частоты (как для критерия Пирсона) и соответствующие им накопленные частоты .
Вычисляется разность накопленных частот для всех интервалов k, максимальное по модулю отклонение теоретических от эмпирических частот обозначено, как .
При гипотеза о распределении может быть принята, иначе — отклоняется, при этом определяется по формуле
, при α = 0,05. (13)
Математическое ожидание значений нормативного износа конструкций, m, может быть его осредненной оценкой. Более жесткой оценкой нормативного износа конструкций будет значение, смещенное относительно математического ожидания (m) на (1÷3)σ в сторону больших значений.
Рис. 4. Гистограмма распределения эмпирических и теоретических частот износа монолитного дока грузоподъёмностью 12 000 т.
Fig. 4. Distribution bar graph of empirical and theoretical wear rates of a cast-in-place dock with a carrying capacity of 12,000 tons.
Рис. 5. Гистограмма распределения накопленных эмпирических и теоретических частот износа монолитного дока грузоподъёмностью 12 000 т.
Fig. 5. Distribution bar graph of accumulated empirical and theoretical wear rates of a cast-in-place dock with a carrying capacity of 12,000 tons.
Рис. 6. Нормативный износ поперечного сечения монолитного дока грузоподъёмностью 12 000 т для общей продольной прочности.
Fig. 6. Specified cross-sectional wear of a cast-in-place dock with a carrying capacity of 12,000 tons for overall longitudinal strength.
ОБОСНОВАНИЕ НОРМАТИВНОГО ИЗНОСА КОНСТРУКЦИЙ ПЛАВУЧЕГО ДОКА В ТРЕБОВАНИЯХ К ОБЩЕЙ ПРОЧНОСТИ И УСТОЙЧИВОСТИ
Таким образом, в (2) требуемый момент сопротивления поперечного сечения к концу срока службы дока () будет соответствовать не 100%-ному износу ЭБ, а экспериментально оцененному в настоящей работе, и с ошибкой в безопасную сторону можно принять равным 70%. Такой подход повлечет за собой пересмотр уровня допускаемых напряжений.
В условии устойчивости (3), предлагается, что расчетные сжимающие напряжения (σсi) должны быть определены с учетом износа поперечного сечения дока, соответствующего 70%, при этом (σe) эйлеровы напряжения i-го элемента ЭБ должны определяться при условии 100%-го износа рассматриваемого элемента.
Настоящие изменения позволят существенно смягчить действующие требования Правил РМРС к общей прочности и устойчивости при продольном изгибе дока.
ЗАКЛЮЧЕНИЕ
По результатам анализа нормативного износа поперечного сечения корпуса плавучего дока с применением специальной методики, основанной на случайном моделировании скоростей износов с использованием датчиков случайных чисел, были получены результаты, которые показали:
- математическое ожидание принимает значения в диапазоне: от 50 до 54%;
- стандартное отклонение принимает значения в диапазоне: 4–5%.
Тогда с вероятностью 99,7% нормативный износ эквивалентного бруса дока, участвующего в общем продольном изгибе, на конец срока службы может быть принят 65–70%.
ДОПОЛНИТЕЛЬНАЯ ИНФОРМАЦИЯ
Вклад авторов. Все авторы внесли существенный вклад в разработку концепции и методологии организации расчетов, реализацию расчетной части и подготовку статьи, прочли и одобрили финальную версию перед публикацией. Личный вклад каждого автора: Д.Т. Иванова — разработка методологии, реализация расчетов, анализ результатов; В.Н. Тряскин — разработка концепции и методологии.
Источники финансирования. Авторы заявляют об отсутствии внешнего финансирования при проведении исследования.
Конфликт интересов. Авторы декларируют отсутствие явных и потенциальных конфликтов интересов, связанных с публикацией настоящей статьи.
ADDITIONAL INFO
Author contributions. All authors made substantial contributions to the conceptualization, investigation, and manuscript preparation, and reviewed and approved the final version prior to publication. Personal contribution: D.T. Ivanova: methodology, investigation, formal analysis; V.N. Tryaskin: conceptualization, methodology.
Funding sources. The authors declare no external funding was received for conducting the study.
Disclosure of interests. The authors have no explicit or potential conflicts of interest associated with the publication of this article.
Об авторах
Дианна Тимофеевна Иванова
Санкт-Петербургский государственный морской технический университет
Автор, ответственный за переписку.
Email: D.T.Ivanova@smtu.ru
SPIN-код: 1314-7682
старший преподаватель кафедры конструкции и технической эксплуатации судов
Россия, Санкт-ПетербургВладимир Николаевич Тряскин
Санкт-Петербургский государственный морской технический университет
Email: vladimir.tryaskin@smtu.ru
SPIN-код: 7663-8210
д-р техн. наук, профессор, заведующий кафедры конструкции и технической эксплуатации судов
Россия, Санкт-ПетербургСписок литературы
- Правила классификации и постройки морских судов. Ч. II. Корпус. Российский Морской Регистр Судоходства. М.: Минтранс, 2024.
- Тряскин В.Н. Проектирование конструкций корпуса плавучего дока по требованиям к общей прочности. В кн.: Регистр СССР: науч.-техн. сб. М.: Транспорт, 1991. Вып. 18. С. 3–18.
- Отчет по научно-исследовательской работе. Ч. I. Обоснование требований Правил Регистра СССР по корпусу металлических плавучих доков. Л.: ЛКИ, 1964.
- Лам Ван Хунг. Методические основы, алгоритмы и программное обеспечение задач проектирования и контроля технического состояния конструкций корпуса судна. дисс. … канд. техн. наук. Санкт-Петербург, 2007.
- Р 50.1.037-2002 Рекомендации по стандартизации. Прикладная статистика, Правила проверки согласия опытного распределения с теоретическим. Ч. II Непараметрические распределения. М.: Госстандарт России, 2002.
- Мицель А.А. Прикладная математическая статистика. Томск: ТУСУР, 2019. EDN: TSMVGT
- Бородина А.В., Некрасова Р.С. Статистические критерии в анализе данных. Петрозаводск: ПетрГУ, 2023.
- Гмурман В.Е. Теория вероятностей и математическая статистика. М.: Высшая школа, 2003. ISBN: 5-06-004214-6 EDN: QJLKXP
- Шилова З.В., Шилов О.И. Теория вероятностей и математическая статистика. Киров: ВГГУ, 2015. ISBN: 978-5-906172-62-4 EDN: UGPGCP
- Екимов В.В. Вероятностные методы в строительной механике корабля. Л.: Судостроение, 1966.
Дополнительные файлы
