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

УДК 623.8 Дата подачи статьи: 08.04.2014

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

В.В. Курятников, к.т.н.

(НИИ «Центрпрограммсистем», просп. 50 лет Октября, 3а, г. Тверь, 1 70024, Россия);

Г.В. Лушин, к.т.н., доцент; М.С. Марков; Б.Г. Иванов

(Военно-морской политехнический институт ВУНЦ ВМФ «ВМА», Кадетский бульвар, 1, г. Санкт-Петербург, г. Пушкин, 196604, Россия, kqfedra_vks@mail.ru).

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

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

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

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

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

ITERATIVE METHOD OF SHIP LANDING AND STABILITY CALCULATION IN SYSTEMS OF THE AUTOMATED SURVIVABILITY ASSESSMENT Kuryatnikov V. V., Ph.D. (Engineering) (Research Institute «Centerprogramsystem», 50 let Oktyabrya Av. 3a, Tver, 170024, Russian Federation);

Lushin G. V., Ph.D. (Engineering), Associate Professor;

Markov M. S.; Ivanov B. G.

Received 08.04.2014

Abstract. The article describes the method of calculating landing and stability of an intact and damaged ship in systems of the automated survivability assessment. The iterative method for a given load is based on V.G. Vlasov's method of successive approximations. It allows solving the problem with given and controlled accuracy degree. It does not solve some spatial problem of in trim landing characterization. It calculates the diagram of crosscut statistic stability of a freely trimmed ship. Calculated parameters of crosscut diagram are used to assess the ship state.

The method uses a mathematical model of the ship. It is specified by the point coordinates on the theoretical surface and watertight compartments. It allows calculating the parameters of landing, stability and others that characterize the safety of the ship for a given load (flooded compartments) with controlled accuracy. The article contains load calculation principles for a damaged ship at various ways of flooding watertight compartments, including doubtful flooded compartments. There are certain formula dependances that allow assessing stability of an intact and damaged ship.The authors conclude that using this iterative method can improve the efficiency and reasonability of decision-making and assignment the measures to restore stability.

Keywords: damaged ship, intact ship, calculation, landing parameters, safety parameters, chart of crosscut static stability, ship loading, waterproof compartments.

Итерационный метод расчета посадки и остойчивости корабля для заданного варианта нагрузки (варианта затопления) основывается на методе последовательных приближений В.Г. Власова [1], позволяющем решать задачу с заданной и контролируемой степенью точности. Основное его отличие от метода В.Г. Власова в том, что в методе расчета, используемом в системе автоматизированной оценки живучести корабля, не решается пространственная задача определения всех трех параметров равновесной посадки корабля Tp, ©p путем системы из трех уравнений равновесия, а рассчитывается диаграмма поперечной статической остойчивости свободно дифферентующего-ся корабля (ДСО).

Итерационный метод расчета посадки и остойчивости корабля

В процессе расчета поперечной ДСО для каждого фиксированного угла крена 0, в диапазоне от -65° до 65° с шагом Д0=О,5° с заданной точностью

|ST\ < 0,01м,

|ST|< 0,0005 рад.

решаются первое и второе уравнения равновесия

— (Т, Т,при ©; = const) =

Р

= V (T, Т, при ©,. = const) = VXc (Т, Т, при ©,. = const) -

—Xg (Т, Т, при ©,. = const) =

Р

= g (Т, Т, при ©,. = const) -

Р

-VZc (Т, Т, при ©; = const))tgy

где M, xg, zg - масса и координаты центра масс заданной нагрузки; V, xc, zc - объемное водоизмещение и координаты его центра величины.

Уравнения равновесия (2) решаются методом последовательных приближений. Причем параметры посадки каждого i-го приближения определяются по формулам

T = T-1+ST

где приращения параметров посадки 5 Т и 5Т, вычисляются по упрощенным формулам перехода:

M

- V

V-1

MXgi-

- Vi-1XOi_

ST =-

, sy = -

vH

p

p

S

В формулах перехода (4) SCT, VCT, Нст - площадь ватерлинии, погруженный объем и продольная метацентрическая высота (МЦВ) при стандартном водоизмещении; индекс i-1 соответствует предыдущему i-1-му приближенно.

Если приращения параметров посадки 5 Ti и 5Y„ вычисленные по формулам (4), удовлетворяют условиям точности (1), то параметры Ti и вычисленные по формулам (3), принимаются за решение продольной плоской задачи (2). Затем вычисляются поперечный восстанавливающий момент (рис. 1) по формуле

m%¡ = у (VYc cos © + VZc sin ©) -

- PZg sin ©- PYg cos ©, тангенс угла наклона касательной к ДСО

dm„

= tga;, угол крена, при котором в во-

d ©

ду входит кромка наружной палубы ©к

f Z

W 1 1

1 G yv

С

1 s zc 1 ,P /

А-У^Г/ У

Рис. 1. Силы, действующие в поперечной плоскости

Fig. 1. Forces that influence in rolling plane

В формуле (5) Р=М^ и yg - сила тяжести и ордината центра масс корабля.

Когда поперечная ДСО пересекает ось углов крена (выполняется условие sign(m@i-1)=-sign(m@i)) и при этом tgaI>0 (остойчивость корабля положительная), методом линейной интерполяции вычисляются параметры равновесной посадки: угол крена ©р, осадка на миделе Тр, угол дифферента

(равновесный дифферент Др) (рис. 2).

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

- запас плавучести Ар с учетом герметичности надводного борта;

- поперечная МЦВ ^ численным дифференцированием поперечной ДСО;

- объем влившейся забортной воды к,вр.;

- минимальная высота надводного борта

min p?

- осадка носом Тнр;

- осадка кормой Ткр.

m0

©i- , ©, ©i m.

m .i+1 M ©

Рис. 2. Определение равновесного угла крена р Fig. 2. 0р equilibrium lurch angle detection

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

При ©,=0 также численным дифференцированием поперечной ДСО вычисляется начальная поперечная МЦВ hо.

При одновременном выполнении условий sign(tgai-1)=-sign(tgai) и tgaI-1<0 (рис. 3) имеет место минимум в ДСО третьего, четвертого и пятого типов [2]. Для этого случая вычисляются ©т1П и т0тп по формулам нелинейной интерполяции.

При одновременном выполнении условий я/^я^тт^^я/^я^тш,) и ^шт,<0 по формулам линейной интерполяции вычисляется угол входа в воду кромки наружной палубы ©крит (критический угол крена поврежденного корабля).

При одновременном выполнении условий sign(tgaI-l)=-sign(tgai) и tgO(i-1>0 (рис. 3) имеет место максимум поперечной ДСО, и для этого случая по формулам нелинейной интерполяции вычисляется предельный угол крена корабля ©пред (предельный угол наклонения корабля статически приложенным кренящим моментом).

При одновременном выполнении условий sign(m@i-1)=-sign(m@i) и tgai-1<0 (рис. 3) поперечная ДСО имеет угол заката (©зак), величина которого вычисляется по формулам линейной интерполяции.

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

Продольная остойчивость поврежденного и неповрежденного корабля рассчитывается следующим образом. В положении равновесия при ©=©p=const вычисляются две ординаты продольной ДСО при у=уР±Ду, причем Ду=0,1°, и численным дифференцированием - продольная МЦВ

Н = I Р

dM

У

d у

У = Ур

Алгоритм решения задачи строится следующим образом. При ©=©Р и у=уР±Ду методом по-

Рис. 3. Определение экстремумов ДСО Fig. 3. Static Stability Diagram extremum detection

следовательных приближений решается первое уравнение равновесия:

— (Т, у = уР ± Ду, ©=©Р) = Р

= V(Т, у = уР ±Ду, © = ©Р).

При этом для приращения осадки стГ, используются формулы перехода (4) и требования точности (1). После вычисления двух значений продольного восстанавливающего момента по трехточечной формуле численного дифференцирования рассчитывается продольная МЦВ.

Элементы нагрузки корабля в уравнениях (2) и (5) вычисляются по формулам:

М = М0 + + ХМнгу + ХРжгЛ,

Xg =

Yg = -

М0Xgo + ХМшг;Хшг + ХМнгуХнг/ + ХРжгк К )к

I ] к

м '

М0^?0 + ^Емшг?ш« + ЕМЛ + ХРжгк Ы )к

I ] к

м '

MoZgo + £Мшп+ 2Мнгу2нг/ + £Ржгк ()к

^=-!-ь-к-'

м

где М0, Xg0, Yg0, Zg0 - масса и координаты центра тяжести корабля порожнем, считываются из БД корабля; Мшп, ХшгЬ Yшгi, Zшгi - масса и координаты центра тяжести /-го штатного твердого груза, счи-тываются из БД корабля для их фактического наличия на борту, вводимого оператором системы автоматизированной оценки живучести корабля; Мщ, Хнф Yнгj, Zнгj - масса и координаты центра тяжести '-го нештатного твердого груза, вводятся оператором системы автоматизированной оценки живучести корабля; ук, хук, zvk - объем жидкого груза (забортной воды) в отделении № k и коор-

динаты его центра величины для текущих значений параметров посадки; значения величин вычисляются для заданной оператором системы автоматизированной оценки живучести корабля; р^. - плотность жидкого груза (забортной воды) в отделении № ^

Особенности вычисления элементов

водонепроницаемых отсеков с жидкими грузами

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

1) неповрежденный водонепроницаемый отсек (цистерна) заполнен жидким грузом плотности ржг на заданное количество (объем);

2) поврежденный ВО затоплен забортной водой по ватерлинию;

3) поврежденный ВО затоплен полностью и связан с забортной водой;

4) поврежденный ВО затоплен на заданный объем;

5) поврежденный ВО затоплен на заданную долю (%) от его полного объема;

6) поврежденный ВО затоплен фильтрационной водой на заданный уровень t метров;

7) поврежденный ВО затоплен по комингс водонепроницаемого закрытия (люка, двери);

k

8) поврежденный ВО имеет надводную пробоину;

9) затопление ВО не установлено - сомнительно затопленный ВО;

10) ВО разгерметизирован (в ВО открыто водонепроницаемое закрытие - люк, дверь, иллюминатор и т.п.).

Второй и третий виды информации о повреждении подпадают под тип затопления ВО по III категории, то есть ВО затоплен по ватерлинию (отсеки 1 и 2 на рисунке 4). Первый, четвертый, пятый и шестой виды информации о заполнении (затоплении) подпадают под тип затопления ВО на заданный объем уоф, в частности:

- при затоплении ВО № ( на заданную долю полного объема п забортной водой плотности р получаем vоф=nxvоп[2] (м3), где уоп - полная вместимость «нетто» ВО № (;

- при заполнении ВО № ( на заданное количество жидким грузом плотности ржт получаем

усф [б ]=т1;

г жг

- при затоплении ВО № ( фильтрационной водой на заданный уровень t (м), который измеряется при 0=0 и ¥=0 над наиболее углубленной его точкой; сначала производится вычисление объема воды vоф (м3) в затопленном объеме ВО № ( по заданный уровень t при 0=0 и ¥=0, а затем расчет элементов по алгоритму вычисления элементов ВО, затопленного на заданный объем vоф.

Z© AZ

W

A / /\

У

Рис. 4. Затопление отсека по ватерлинию

Fig. 4. Flooding a compartment up to water line

Главной особенностью вычисления элементов отсека, затопленного по ватерлинию, является то, что в процессе расчета поперечной ДСО уровень воды в ВО всегда совпадает с уровнем ватерлинии (рис. 4).

Элементы заполненного объема ВО рассчитываются интегрированием по правилу трапеции по длине ВО элементов погруженных площадей расчетных шпангоутов ВО.

Главной особенностью алгоритма расчета элементов отсека №

на заданный объем vоф[2] является то, что в процессе расчета поперечной ДСО при каждом фиксированном угле крена 0, объем жидкого груза (забортной воды) в отсеке должен оставаться постоянным и с заданной точностью равным объему Vоф[2] (рис. 5).

Рис. 5. Затопление отсека на заданный объем Fig. 5. Flooding a compartment according to given volume

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

кф[2]- ^ I <0,01 ^[2],

где vоф[2] - заданный оператором объем жидкого груза (забортной воды) в отсеке № (; v0i - объем жидкого груза (забортной воды) при /-м приближении; vоп [(] - полный объем «нетто» отсека № 2.

Главной особенностью алгоритма расчета элементов водонепроницаемого отсека, затопленного по уровень ^ является то, что сначала при 0=0 и ¥=0 для заданного уровня затопления t в процессе ввода данных вычисляется объем воды в отсеке vоф[2]. Последующий расчет происходит при условии, что в ВО № 2 объем воды все время постоянный и равен voф, то есть последующий расчет идет в соответствии с алгоритмом расчета элементов отсека затопленного (заполненного) на заданный объем.

Главной особенностью расчета элементов ВО, затопленного по комингс водонепроницаемого закрытия, является то, что при заданной посадке корабля Т, ¥, 0 уровень воды в ВО проходит всегда по кромку комингса параллельно плоскости ватерлинии, и соответственно элементы забортной воды в ВО вычисляются по эту кромку.

Координаты кромки комингса

считываются из БД корабля.

Главной особенностью алгоритма расчета элементов отсека, имеющего одну или несколько надводных пробоин, является то, что отсек остает-

ся пустым, пока кромка одной из пробоин не вошла в воду, а когда она входит в воду (гпроб<Т0проб), отсек заполняется по ватерлинию и расчет продолжается по алгоритму вычисления элементов объема отсека, затопленного по ватерлинию.

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

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

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

Координаты разгерметизированного водонепроницаемого закрытия (Х„„г, 7„„г, Zvnz) для введенных оператором разгерметизированных ВО считываются из БД корабля.

Главная особенность алгоритма расчета элементов сомнительно затопленного ВО в том, что его элементы уга=0, тхуо1=0, туго,=0 рассчитываются по худшему, но реальному варианту затопления. В соответствии с этим принимаются во внимание следующие положения:

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

- при расчете величины кренящего момента, чтобы не недооценить величину удерживающего момента [2] при пятом типе поперечной ДСО, сомнительно затопленный ВО на вышедшем из воды борту считается затопленным по ватерлинию;

- при расчете величины кренящего момента, чтобы не переоценить величину спрямляющего момента при спрямлении поврежденного корабля, сомнительно затопленный ВО на вошедшем в воду борту считается пустым [3];

- при расчете величины начальной поперечной МЦВ сомнительно затопленный междудон-

ный ВО считается незаполненным, а ВО на палубах выше второго дна - затопленным фильтрационной водой на 10 % его высоты.

Главная идея расчета элементов отсека, затопленного по показаниям датчиков уровня, состоит в следующем:

- если в отсеке установлен датчик давления (датчик уровня) и он сработал, то элементы отсека при текущей равновесной посадке корабля вычисляются по тот уровень, который этот датчик регистрирует, и при последующем вычислении поперечной ДСО данный отсек учитывается как отсек, затопленный на заданный объем;

- если в отсеке установлен датчик нижнего уровня (сигнализатор) и он сработал, то элементы отсека при текущих равновесных параметрах посадки вычисляют по уровень воды, касательный к этому датчику, и в дальнейшем расчете поперечной ДСО этот отсек учитывается как затопленный на заданный объем;

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

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

Литература

1. Власов В.Г. Статика корабля. Л.: СУДПРОМ ГИЗ, 1961. Т. 6. С. 205-211.

2. Муру Н.П. Статика корабля. Л.: Изд-во ВВМИУ им. Ф.Э. Дзержинского, 1982. С. 169-173.

3. Шлемов А.Ф., Коротков Б.П., Горев Ю.Н. Теория корабля. СПб: Кораблестроение, 2007. С. 60.

4. Борисов Р.В., Луговский В.В., Мирохин Б.В., Рождественский В.В. Статика корабля: учеб. пособие. СПб: Судостроение, 2005.

5. Мальцев Н.Я., Дорогостайский Д.В., Прытков Ю.К. Теория непотопляемости судна. Л.: Судостроение, 1973. 320 с.

References

1. Vlasov V.G. Statika korablya [The ship statics]. SUDPROM Publ., vol. 6, Leningrad, 1961, pp. 205-211.

2. Muru N.P. Statika korablya [The ship statics]. Leningrad, F.E. Dzerzhinskiy Naval Engineer College Publ., 1982, pp. 169-173.

3. Shlemov A.F., Korotkov B.P., Gorev Yu.N. Teoriya ko-rablya [The ship theory]. St. Petersburg, Korablestroenie Publ., 2007, 60 p.

4. Borisov R.V., Lugovskiy V.V., Mirokhin B.V., Rozh-destvenskiy V.V. Statika korablya [The ship statics]. Study guide, St. Petersburg, Korablestroenie Publ., 2005.

5. Maltsev N.Ya., Dorogostayskiy D.V., Prytkov Yu.K. Teo-riya nepotoplyaemosti sudna [The theory of the ship resistance to flooding]. Leningrad, Sudostroenie Publ., 1973, 320 p.