Введение При перевозке тарно-штучных грузов нередко необходимо производить их крепление, учитывая массу груза и действующие на него силы инерции. Для расчета сил инерции следует знать линейные ускорения, воздействующие на груз, которые зависят от закономерностей изменения линейных перемещений груза вместе с палубой судна относительно опорной системы координат. В общем случае движение судна характеризуется шестью степенями свободы и описывается шестью дифференциальными уравнениями. Все виды качки судна взаимосвязаны [1]. В [2] показано, что для расчета линейных ускорений в первом приближении можно применить линейные модели бортовой, килевой и вертикальной качки, используя соответствующие изолированные уравнения. Для того чтобы найти силы инерции, действующие на единицу груза в инерциальной системе отсчета, достаточно применить второй закон Ньютона, при условии, что масса этой единицы груза и соответствующего ускорения известны. Тогда ускорения будут определены как вторая производная по времени функций угла крена, угла дифферента и амплитуды вертикальной качки. Выражения для углов и амплитуды могут быть выведены из решения уравнений бортовой, килевой и вертикальной качки. Это возможно на основании утверждения о том, что уравнения поперечной, продольной и вертикальной качки в ситуации поиска сил инерции, воздействующих на груз, можно рассматривать независимо [3]. Параметры качки Так как бортовая качка порождает доминирующие силы инерции, вначале рассмотрим уравнение поперечной качки для определения выражения угла крена θ. Для этого воспользуемся исходным дифференциальным линейным уравнением, описывающим угол крена θ судна, которое приведено в [4]: (1) где Jx – момент инерции судна относительно продольной оси Х-Х; mx – обобщенные присоединенные массы воды относительно оси Х-Х; μx – коэффициент демпфирования относительно оси Х–Х; D – водоизмещение судна (сила тяжести); h0 – поперечная начальная метацентрическая высота; Χθ – редукционный коэффициент при бортовой качке; ωk – кажущаяся частота волнения. Разделив все члены уравнения (1) на коэффициент при старшей производной, получим , (2) где h – коэффициент затухания. . ω0 – собственная частота судна при бортовой качке. . Уравнение (2) является линейным неоднородным дифференциальным уравнением с постоянными коэффициентами, и его решение является суммой частного решения θr, описывающего вынужденные колебания судна относительно оси Х-Х под воздействием регулярного волнения, и решения соответствующего однородного уравнения, которое характеризует собственные затухающие колебания судна. Так как амплитуда собственных затухающих колебаний судна достаточно быстро обращается в нуль, то уравнение бортовой качки, как стационарный процесс, согласно [4], можно характеризовать только вынужденными колебаниями, т. е. (3) В [3, 4] показано, что, в случае продольной качки, судно совершает вынужденные колебания с частотой качки ωk, так же как и при поперечной. Изолированное уравнение продольной качки, а равно и его решение, по своей структуре аналогично уравнению поперечной качки, т. е. описывает как собственные затухающие колебания судна, так и вынужденные гармонические колебания с частотой качки. Вследствие этого выражение для текущего значения угла дифферента β также описывается вынужденными гармоническими колебаниями с частотой качки ωk: , (4) где Χβ – редукционный коэффициент при килевой качке; ω0β – собственная частота судна при килевой качке; hβ – коэффициент затухания при килевой качке. Отметим, что вертикальная поступательная качка, согласно [3, 5], вызвана орбитальным движением судна, которое считается круговым с радиусом, равным половине высоты волны. Вертикальное поступательное движение ζ носит гармонический характер с частотой качки ωk и описывается выражением (5) где ζ0 = 0,5 hw – амплитудное значение вертикального перемещения, причем hw – высота волны. Выражения (3)–(5) позволяют рассчитать угловые ускорения от килевой и бортовой качки, а также найти линейные ускорения и силы инерции, воздействующие на груз. Формализация сил инерции Найдем силы инерции, возникающие от бортовой, килевой и вертикальной качки, которые действуют на груз с массой mc. Наиболее значительной является поперечная сила инерции от бортовой качки Fθ. Очевидно, что где ay – линейное ускорение при бортовой качке. В свою очередь, линейное ускорение ay является произведением углового ускорения на радиус кривизны ry относительно продольной оси, проходящей через центр тяжести судна G, т. е. . Угловое ускорениенайдем как вторую производную угла крена, для чего дважды дифференцируем выражение (3) и получим , где , . Следовательно, сила инерции Fθ равна: . Продольную силу инерции Fβ рассчитываем аналогично, т. е. где ax – линейное ускорение при килевой качке. Линейное ускорение , причем rx – радиус кривизны относительно поперечной оси. Угловое ускорение получим, дважды дифференцируя выражение (4): , где , . Поэтому сила инерции Fβ описывается выражением . Сила инерции от вертикальной поступательной качки Fζ определяется из выражения , причем линейное ускорение получаем с помощью двойного дифференцирования выражения (5): , поэтому с учетом того, что ζ0 принимается равным половине высоты волны hw, т. е. , окончательно получим На рисунке приводится взаимосвязь опорной (невозмущенной) системы координат OXYZ, относительно которой получены выражения сил инерции, и судовой OX1Y1Z1, в которой будет производиться расчет реакций в связях крепления груза. Связь систем координат Вследствие этого для расчета реакций в связях крепления груза силы инерции Fθ, Fβ, Fζ и силу тяжести Pc необходимо спроектировать на судовую систему координат, которая отклонена от опорной на углы крена θ и дифферента β. Заключение Таким образом, нами описан метод определения сил инерции, действующих на груз, для последующего учета их значений в расчетах максимальной рабочей нагрузки для креплений тарно-штучного груза. Метод основан на предположении, что в рамках указанной задачи линейные модели бортовой, килевой и вертикальной качки применимы и могут считаться независимыми. Полученные выражения угла крена, угла дифферента и амплитуды вертикальной качки на волне позволяют рассчитать соответствующие угловые и линейные ускорения. Силы инерции определены в невозмущенной системе отсчета и могут быть легко перенесены в судовую систему координат, т. к. соотношение между двумя системами координат известно. Определенные зависимости в дальнейшем будут применяться нами в математической модели процесса безопасной укладки и крепления тарно-штучных грузов на судне.