Теория движения составных стержней

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

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

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

Рис. 2.3.1

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

, ,

, , (2.3.1)

, .

Аэродинамические нагрузки будем считать зависящими не только от перемещений стабилизатора x 0, y0и φ 0, перемещений руля x, y и φ (2.2.2), но и от их производных. Запишем их в решении в соответствии с [16], понимая, что лучше было бы для этого использовать метод дискретных вихрей С.М. Белоцерковского. Следует заметить, что в 70-е годы, когда строилась приведенная теория, еще не было соответствующей вычислительной техники, а в последнее время в России предельно ослаб интерес не только к новым теориям, но и к авиационной технике вообще.

Для построения инерционных характеристик в произвольном сечении стержня (рис.2.3.2) закоординируем положение центра масс (ц.м.) в подвижных осях ξ η с началом в центре жесткости сечения вектором .

Перемещения ц.м. связаны с перемещениями ц.ж.

(2.3.2)

и определяют его ускорение (ц.м. руля)

(2.3.3)

Для стабилизатора можно записать

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

Рис. 2.3.2

Если и - погонные массы стержней (руля и стабилизатора), то погонная массовая нагрузка

, ,

, (2.3.4)

Здесь , ,

- масса единицы площади плоского стержня (руля, крыла)

и погонный массовый момент инерции стержня соответственно

, .

Запишем приближенные моменты и перерезывающие силы деформированных стержней на примере рулевой поверхности в соответствии с (2.2.3)

,

,

,

(2.3.5)

к которым ещё добавятся константы.

Аналогичные уравнения записываются для стабилизатора.

Выражения (2.3.5), представляющие собой при нижнем пределе интегрирования уравнения равновесия, вместе с условиями совместности(2.2.5), физическими зависимостями Кирхгофа-Клебша (2.1.3), уравнениями связи угловых и линейных перемещений(2.1.7) и связи перемещений x, y, z с деформациями æ x, æ y, τ z(2.1.5), записанными для руля и для стабилизатора, составляют разрешающие уравнения движения. Записанные через интегрирующие матрицы и приведенные к системе алгебраических эти уравнения могут быть представлены в виде вектор- функции

(2.3.6)

где .

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

Тогда, .

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

.

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

;

.

Решение уточняется итерационным методом Ньютона-Канторовича. На рис. 2.3.3 представлен фрагмент решения. Приводится сравнение нормальной и лобовой кривизны руля во времени.

Рис 2.3.3

Детальное обсуждение рис. 2.3.3. приведено ниже вместе с результатами эксперимента.