Введение Построение систем поддержки принятия решений (СППР) для танкеров, перевозящих сжиженные газы, является сложной многофакторной научно-прикладной задачей, поскольку возникновение и развитие опасных ситуаций определяется большим количеством различных факторов, зависящих от ряда внешних, технических и технологических условий, значительная часть которых не поддается прогнозированию. Выработка СППР корректных рекомендаций для персонала, управляющего комплексом судовых технических средств (СТС), осуществима лишь при возможности прогнозирования будущих состояний СТС и всех контролируемых параметров судна, что невыполнимо без наличия в составе системы соответствующих математических имитационных моделей, позволяющих выполнять прогнозирование с учетом всех возможных статических и динамических воздействий. Категории параметров, контролируемых системой поддержки принятия решений В контексте проблематики исследования следует выделить четыре категории контролируемых параметров: - входные; - технологические; - функциональные; - диагностические. К входным параметрам СТС относятся исходные данные, характеризующие условия работы СТС в текущий момент времени. Например, для главного двигателя это мгновенные значения частоты вращения, нагрузки и пр. К группе технологических параметров относятся параметры, непосредственно характеризующие процессы, протекающие в рассматриваемом узле или механизме. Для двигателя это давление в масляной системе, рабочая температура, частота вращения и пр. К группе функциональных параметров следует относить так называемые «внешние характеристики»; тип функциональных параметров определяется характером режима функционирования конкретного СТС, например, для детерминированных систем это значение их передаточной функции. Все три вышеперечисленные группы параметров достаточно тесно связаны между собой, при этом во многих случаях эта связь носит комплексный многофакторный характер, часто нелинейный и вариативный. Более того, в некоторых частных случаях при различных взаимных сочетаниях некоторых параметров существенно изменяется и вид их зависимостей. К группе диагностических параметров следует относить те параметры, которые напрямую не вносят вклад в непосредственное функционирование СТС или их взаимодействие, но при обеспечении безопасности судна-газовоза имеют высокую информационную ценность. К этой группе параметров в первую очередь следует относить данные систем ранней углубленной диагностики узлов и агрегатов. В наших предыдущих исследованиях были достаточно подробно рассмотрены методики первичной параметрической, функциональной и глубокой ранней диагностики, ранее не применявшиеся на судах, реализация которых позволила бы существенно повысить отказоустойчивость СТС, значительно снизить вероятность возникновения внезапных неисправностей, а также повысить экономическую и технологическую эффективность технического обслуживания и ремонта применительно к конкретным СТС. Основные факторы опасности для судов-газовозов В контексте обеспечения безопасности эксплуатации танкеров, предназначенных для перевозки сжиженных газов, необходимо подробно рассмотреть особенности контроля параметров грузовых систем (в первую очередь, непосредственно грузовых танков), для которых можно выявить ряд критичных факторов: - давление груза; - давление внешней среды; - динамические нагрузки, вызываемые движением судна; - термические нагрузки; - нагрузки, возникающие в результате всплескивания сжиженного газа; - нагрузки от изгибов и прогибов корпуса судна; - вес груза, емкостей, изоляции, сопряженного технологического оборудования, ветровые нагрузки на надпалубные части танков и соответствующие реакции опор и креплений. Представляется очевидным, что наибольшей степенью опасности характеризуются динамические составляющие нагрузок, при этом следует отметить, что вышеприведенный перечень нагрузок можно классифицировать по важнейшему признаку - наличию возможности регулирования воздействующего фактора. Для нерегулируемых факторов (гравитационная составляющая нагрузок, ветровые нагрузки, температура окружающей среды) необходимо реализовать систему мониторинга и сигнализации, позволяющую в случае достижения контролируемыми параметрами граничных значений реализовать подготовку нештатных мер реагирования, таких как смена курса и пр. Для регулируемых факторов необходимо реализовать возможность соответствующих управляющих воздействий. Динамические воздействия на грузовую систему Вне зависимости от типов используемых танков, в соответствии с требованиями [1], для грузовой системы газовоза необходимо обеспечить соблюдение ряда технологических условий, связанных с параметрами перевозимого груза. Одним из наиболее важных параметров является внутреннее давление жидкости (сжиженного газа), которое возникает в результате действия ускорения силы тяжести груза, приложенного в его центре, в процессе движения или маневрирования судна, а также при продольно-горизонтальной, поперечно-горизонтальной, вертикальной, бортовой и килевой качках и при рыскании на нерегулярном волнении. В связи с невозможностью воздействия на причины возникновения таких нагрузок единственным эффективным методом ослабления вышеперечисленных факторов становится снижение скорости и изменение курса судна, при этом СППР должна сформировать соответствующие рекомендации для оператора. Величина внутреннего давления в танках, определяющаяся комбинированным воздействием гравитационного и динамического ускорений, рассчитывается согласно соотношению , где αβ - безразмерное ускорение, выражаемое в долях от ускорения силы тяжести, которое возникает в результате воздействия гравитационных и динамических нагрузок в произвольном направлении β; ρt - максимальная плотность груза при заданной расчетной температуре; Zβ - наибольшая высота жидкости над точкой, в которой определяется внутреннее давление, измеряемая от обшивки грузового танка в направлении вектора β, при этом объем куполов танка, рассматриваемых как часть принятого объема всего танка, должен приниматься в расчет при определении Zβ только в тех случаях, когда объем всех куполов танка Vd превышает значение, определяемое по формуле , где Vt - объем танка без куполов, FL - предел заполнения танка. Следует дополнительно отметить, что в расчет принимается то направление действия сил, в котором наблюдается максимальное значение Pgd. При разработке СППР необходимо предусмотреть контроль величины накапливаемой усталостной нагрузки, которая может быть определена в соответствии с расчетным соотношением , где ni - число циклов напряжений при каждом уровне напряжений в течение всего срока службы танка; Ni - число циклов до излома для соответствующего уровня напряжений по кривой Велера (S-N); Nj - число циклов до излома для усталостных нагрузок, возникающих при погрузочно-разгрузочных операциях; CW - величина, определяемая в зависимости от применяемого метода испытаний конструкций танков, используемого для построения кривой Велера, при этом стандартное значение CW ≤ 0,5, в обоснованных случаях [1] допускается отклонение значения 0,5 ≤ CW ≤ 1. Кроме того, необходимо осуществлять непрерывный контроль суммарных напряжений: , , , где индексы *.st соответствуют статическим, а *.dyn - динамическим составляющим напряжений корпуса судна и танков при прогибах и кручении. Для расчета величин ускорений, возникающих в результате воздействия качки судна, используется группа расчетных соотношений, справедливых для судов длиной более 50 м: - вертикальное ускорение: - поперечное ускорение: - продольное ускорение: , где L0 - длина судна, значение которой используется для определения размеров поперечных сечений связей корпуса в соответствии с принятыми стандартами; СB - коэффициент общей полноты; В - наибольшая теоретическая ширина судна; х - продольное расстояние от миделя до центра тяжести танка с его содержимым (при этом величина х характеризуется положительными значениями в нос от миделя и отрицательными - в корму от миделя); z - расстояние по вертикали от фактической ватерлинии судна до центра тяжести танка с его содержимым (при этом величина z характеризуется положительными значениями выше ватерлинии и отрицательными - ниже). Значение аz не включает составляющую ускорения статического веса; аy включает составляющую ускорения статического веса, действующую в поперечном направлении в результате бортовой качки; аx включает составляющую ускорения статического веса, действующую в продольном направлении в результате воздействия килевой качки. При этом ; ; , где V - эксплуатационная скорость; K - коэффициент, обычно принимаемый равным 1 (при этом его величина находится в зависимости от конкретных условий загрузки и характеристик обводов корпуса); GM - метацентрическая высота. Динамические нагрузки на силовую установку Одним из наиболее неблагоприятных внешних факторов, оказывающих существенное влияние на рабочие характеристики, техническое состояние и, как следствие, уровень безопасности эксплуатации, для всех типов двигателей является внезапные набросы и сбросы нагрузок. В результате таких воздействий возникают знакопеременные динамические и инерционные моменты, величина которых может соответствовать и даже в некоторых случаях превышать аналогичные рабочие значения, что оказывает крайне негативное влияние на двигатель внутреннего сгорания (ДВС). Оценку влияния наброса нагрузки можно выполнить с использованием уравнения вращающихся масс ДВС [2]: , где - постоянная времени разгона (J - момент инерции вращающихся масс; TtqH - крутящий момент на номинальном режиме; ωH - угловая скорость на номинальном режиме); - относительное приращение величины ω; (ηМЕХ.H - механический КПД на номинальном режиме); - величина относительного индикаторного КПД; - относительные механические потери; - коэффициент, учитывающий изменение механических потерь; - давление наддува; PeH - среднее эффективное давление; - величина, находящаяся в зависимости от закона регулирования нагрузки, как правило, принимающаяся в пределах 0,995-1,003. Для контроля и выявления ситуаций, связанных с набросом значительных нагрузок, при которых в главном двигателе возникают значительные механические напряжения, необходимо в составе проектируемой СППР создать адаптивный модуль, отвечающий за повышение оборотов двигателя перед набросом нагрузок. Математическое описание любого ДВС характеризуется, в первую очередь, основным целесообразным критерием идентификации [3]: , где L - функционал от четной функции q(e); e - погрешность идентификации: e = y - w; w = G [u; α], где u - входное управляющее воздействие, x и w - выходные процессы объекта исследования и модели, G - оператор связи. Настройка G осуществляется посредством изменения параметров αT = (α1, α2, …, αm) в соответствии со значением градиента E: , где γ = α(0) - начальное условие. Таким образом, ; , при этом представляет собой функцию чувствительности параметра αj. Градиент Е определяется расчетным соотношением, где . Множество K[u, α] позволяет получить все функции чувствительности параметров α. Использование данных расчетных соотношений позволяет определить все функции чувствительности изменяемых параметров виртуальной модели, характеризующей динамику исследуемого ДВС, позволяющей идентифицировать все мгновенные состояния двигателя. Детерминирование функции чувствительности позволяет впоследствии при использовании математической модели исследуемого двигателя производить анализ его частотных и временных характеристик, а также выявлять основные зависимости, определяющиеся варьированием значений исходных параметров и прочих внешних факторов. Для такого анализа представляется целесообразным применение функции чувствительности в логарифмической форме - логарифмической функции чувствительности (ЛФЧ), поскольку она позволяет характеризовать величину относительного изменения выходного параметра R (или его оценки) для исследуемой модели, которое вызвано соответствующим изменением i-го показателя или модели bi. В первом приближении линейная функция чувствительности характеризуется следующей зависимостью [3]: где представляет собой искомую ФЧ, при этом нулевое значение индекса означает, что частные производные располагаются рядом с начальным значением bi(0). В случае одновременного изменения m параметров исследуемой математической модели ДВС величины абсолютного и относительного изменения величины R будут соответствовать соотношениям . Для дальнейшей реализации модели ДВС необходимо рассмотреть следующие зависимости: где ω(t) - динамика угловой скорости, ε(t) - динамика углового ускорения коленчатого вала в режимах, близких к квазистатическому режиму (ω = ω*). В этом случае уравнение чувствительности модели принимает следующий вид: , или , где и . Заключение Результаты исследования позволяют с высокой степенью точности синтезировать математические модели грузовой системы и ДВС с учетом внешних и технологических динамических воздействий. Продолжением работы является разработка и апробирование моделей остальных СТС, так или иначе задействованных в обеспечении безопасности эксплуатации танкеров, перевозящих сжиженные газы.