12

Решение жестких краевых задач строительной механики (расчет оболочек составных и со шпангоутами) методом виноградовых (без ортонормирования)

Алексей Юрьевич Виноградов, кандидат физико-математических наук, freelance

Юрий Иванович Виноградов, доктор физико-математических наук,

профессор МГТУ им. Н.Э.Баумана, Москва, AlexeiVinogradov@yandex.ru

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

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

SOLUTION OF STIFF BOUNDARY VALUE PROBLEMS OF BUILDIND MECHANICS (CALCULATION OF SHELLS COMPOUND AND WITH FRAMES) BY VINOGRADOVS’ METHOD (WITHOUT ORTHONORMALIZATION)

Alexei Yurievich Vinogradov, candidate of physical and mathematical sciences, freelance

Yuri Ivanovich Vinogradov, doctor of physical and mathematical sciences,

professor of MSTU n.a. Bauman, Moscow, AlexeiVinogradov@yandex.ru

The simplest method of solution of the stiff boundary value problems of building mechanics is offered for calculation of compound shells and shells with frames. The method does not need orthonormalization because of dividing an integration interval into matching sectors.

Keywords: stiff boundary value problems, shells, frames, orthonormalization.

1. Введение.

Введение дается на примере системы дифференциальных уравнений цилиндрической оболочки ракеты – системы обыкновенных дифференциальных уравнений 8-го порядка (после разделения частных производных методом Фурье).

Система линейных обыкновенных дифференциальных уравнений имеет вид:

,

где – искомая вектор-функция задачи размерности 8х1,– производная искомой вектор-функции размерности 8х1,– квадратная матрица коэффициентов дифференциального уравнения размерности 8х8,– вектор-функция внешнего воздействия на систему размерности 8х1.

Краевые условия имеют вид:

где – значение искомой вектор-функции на левом крае х=0 размерности 8х1,– прямоугольная горизонтальная матрица коэффициентов краевых условий левого края размерности 4х8, – вектор внешних воздействий на левый край размерности 4х1,

– значение искомой вектор-функции на правом крае х=1 размерности 8х1,– прямоугольная горизонтальная матрица коэффициентов краевых условий правого края размерности 4х8, – вектор внешних воздействий на правый край размерности 4х1.

В случае, когда система дифференциальных уравнений имеет матрицу с постоянными коэффициентами =const, решение задачи Коши имеет вид [1]:

,

где , где- единичная матрица.

Матричная экспонента ещё может называться матрицей Коши или матрициантом и может обозначаться в виде . Тогда решение задачи Коши может быть записано в виде,

где это вектор частного решения неоднородной системы обыкновенных дифференциальных уравнений.

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

Вычисление этого вектора частного решения неоднородной системы дифференциальных уравнений производиться при помощи представления матрицы Коши под знаком интеграла в виде ряда и интегрирования этого ряда поэлементно. Это справедливо для случая системы дифференциальных уравнений с постоянной матрицей коэффициентов =const. Для случая дифференциальных уравнений с переменными коэффициентами в приведенной выше формуле для каждого участка может использоваться осредненная матрицакоэффициентов системы дифференциальных уравнений.

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

Если участок не мал, то его можно поделить на подучастки и тогда можно предложить следующие рекуррентные (итерационные) формулы для вычисления частного вектора: имеем, также имеем формулу для отдельного подучастка. Тогда:,.

Подставим ви получим:

.

Сравним это выражение с формулой

и получим, очевидно, что и для частного вектора получаем формулу. То есть вектора подучастковне просто складываются друг с другом, а с участием матрицы Коши подучастка.

Аналогично можем записать и подставить сюда формулу дляи т.д.

Из теории матриц [1] известно свойство перемножаемости матричных экспонент (матриц Коши): .

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

,

где матрицы Коши приближенно вычисляются по формуле:

, где.

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

В тоже время можно производить вычисления и с заранее известной точностью. Для этого следует вычислять матрицы Коши не как матричные экспоненты от осредненных аргументов, а можно использовать для вычисления векторов, входящих в матрицы Коши, методы типа методов Рунге-Кутты. Для вычисления матриц Коши методами типа Рунге-Кутты используется стартовая (начальная) единичная матрица. А для вычисления вектора частного решения неоднородной системы ОДУ берется начальный нулевой вектор. В случае применения методов типа Рунге-Кутты оценки погрешностей хорошо известны.