Scientific journal
Modern high technologies
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

REFINED METHODS OF BOUNDARY CONDITION TRANSFER VINOGRADOV IN AN ARBITRARY POINT OF THE INTERVAL OF INTEGRATION FOR SOLVING STIFF BOUNDARY VALUE PROBLEMS

Vinogradov Yu.I. 1 Vinogradov A.Yu. 1
1 South Innovation and Technology University
The article defines the use of mathematical modeling techniques for the development of development of information systems in the aerospace industry. Compared with the above previously described method initially are known formula of the theory of matrices for systems of ordinary differential equations (ODE), the entrance is considered a variant of the method for non-rigid boundary-written formula incremental transfer of boundary conditions at the point of the interval of integration for non-rigid and rigid boundary value problems are other progressive formula orthonormality portable matrix equations of the boundary conditions of rigid boundary value problems, is another version of the calculation of the vector of a particular solution of inhomogeneous systems of hard ODE boundary value problems. The results of test calculations coincided with the results of another method of Vinogradov – method for solving boundary value problems without hard orthonormality. We describe the method used in spacecraft flight paths Conditional Split on specific areas, determining the best management programs and the pairing of the results. Application of the method avoids complex computational procedures used in the calculation of optimal trajectories of classical methods. The results of testing of the developed method in calculating trajectories for removing vertical and horizontal starts spacecraft.
rigid boundary value problems
the transfer of boundary conditions
method of calculation of the flight

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

vin01.wmf,

где Y(x) – искомая вектор-функция задачи размерности 8х1, vin02.wmf – производная искомой вектор-функции размерности 8х1, A – квадратная матрица коэффициентов дифференциального уравнения размерности 8х8, F(x) – вектор-функция внешнего воздействия на систему размерности 8х1. Здесь и далее векторы обозначаем жирным шрифтом вместо черточек над буквами. Краевые условия имеют вид

vin03.wmf

где Y(0) – значение искомой вектор-функции на левом крае х = 0 размерности 8х1, U – прямоугольная горизонтальная матрица коэффициентов краевых условий левого края размерности 4х8, u – вектор внешних воздействий на левый край размерности 4х1, Y(1) – значение искомой вектор-функции на правом крае х = 1 размерности 8х1, V – прямоугольная горизонтальная матрица коэффициентов краевых условий правого края размерности 4х8, v – вектор внешних воздействий на правый край размерности 4х1. В случае, когда система дифференциальных уравнений имеет матрицу с постоянными коэффициентами A = const, решение задачи Коши имеет вид [10]:

vin04.wmf,

где

vin05a.wmf

vin05b.wmf,

где E – это единичная матрица.

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

vin06.wmf.

Тогда решение задачи Коши может быть записано в виде:

vin07.wmf,

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

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

vin09.wmf.

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

vin10.wmf,

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

vin11.wmf,

где vin12.wmf.

Рассмотрим метод «переноса краевых условий» (прямой вариант метода) для решения краевых задач с нежесткими обыкновенными дифференциальными уравнениями. Предлагается выполнять интегрирование по формулам теории матриц [2] сразу от некоторой внутренней точки интервала интегрирования к краям:

vin13.wmf,

vin14.wmf.

Подставим формулу для Y(0) в краевые условия левого края и получим

vin15.wmf,

vin16.wmf,

vin17.wmf.

Аналогично для правых краевых условий получаем

vin18.wmf,

vin19.wmf,

vin20.wmf.

То есть получаем два матричных уравнения краевых условий, перенесенные в рассматриваемую точку x:

vin21.wmf,

vin22.wmf.

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

vin23.wmf.

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

vin24.wmf.

Или можно записать

vin25.wmf.

Подставляем это выражение для Y(0) в краевые условия левого края и получаем

vin26.wmf,

vin27.wmf,

vin28.wmf.

Или получаем краевые условия, перенесенные в точку x1:

vin29.wmf,

где vin30.wmf и vin31.wmf.

Далее запишем аналогично

vin32.wmf.

И подставим это выражение для Y(x1) в перенесенные краевые условия точки x1:

vin33.wmf,

vin34.wmf,

vin35.wmf.

Или получаем краевые условия, перенесенные в точку x2:

vin36.wmf,

где

vin37.wmf

и

vin38.wmf.

И так в точку x* переносим матричное краевое условие с левого края и таким же образом переносим матричное краевое условие с правого края. Покажем шаги переноса краевых условий правого края. Можем записать

vin39.wmf.

Подставляем это выражение для Y(1) в краевые условия правого края и получаем

vin40.wmf,

vin41.wmf,

vin42.wmf.

Или получаем краевые условия правого края, перенесенные в точку xn –1:

vin43.wmf,

где

vin44.wmf

и

vin45.wmf.

Далее запишем аналогично

vin46.wmf.

И подставим это выражение для Y(xn –1) в перенесенные краевые условия точки xn –1:

vin47.wmf,

vin48.wmf,

vin49.wmf.

Или получаем краевые условия, перенесенные в точку xn –2:

vin50.wmf,

где vin51.wmf и vin52.wmf.

И так во внутреннюю точку x* интервала интегрирования переносим матричное краевое условие, как показано, и с левого края и таким же образом переносим матричное краевое условие с правого края и получаем

vin53.wmf, vin54.wmf.

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

vin55.wmf.

Известны формулы ортонормирования систем линейных алгебраических уравнений [2]. В случае «жестких» дифференциальных уравнений предлагается применять построчное ортонормирование матричных краевых условий в процессе их переноса в рассматриваемую точку [3, 5]. То есть, получив vin56.wmf, применяем к этой группе линейных алгебраических уравнений построчное ортонормирование и получаем эквивалентное матричное краевое условие: U1ортоY(x1) = u1орто. И теперь уже в это проортонормированное построчно уравнение подставляем

vin57.wmf.

И получаем

vin58.wmf,

vin59.wmf.

Или получаем краевые условия, перенесенные в точку x2:

vin60.wmf,

где

vin61.wmf

и

vin62.wmf.

Теперь уже к этой группе линейных алгебраических уравнений применяем построчное ортонормирование и получаем эквивалентное матричное краевое условие:

vin63.wmf.

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

vin64.wmf.

Вместо формулы для вычисления вектора частного решения неоднородной системы дифференциальных уравнений в виде [5]

vin65.wmf

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

vin66.wmf.

Правильность приведенной формулы подтверждается следующим:

vin67.wmf,

vin68.wmf,

vin69.wmf.

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

vin70.wmf

Эта формула справедлива для случая системы дифференциальных уравнений с постоянной матрицей коэффициентов A = const.

Рассмотрим вариант, когда шаги интервала интегрирования выбираются достаточно малыми, что позволяет рассматривать вектор F(t) на участке (xj – xi) приближенно в виде постоянной величины F(xi) = constant, что позволяет вынести этот вектор из-под знаков интегралов:

vin71.wmf

Известно, что при T = (at + b) имеем vin72.wmf (при n ≠ –1).

В нашем случае имеем vin73.wmf (при n ≠ –1).

Тогда получаем vin74.wmf.

Тогда получаем ряд для вычисления вектора частного решения неоднородной системы дифференциальных уравнений на малом участке (xj – xi):

vin75.wmf

Для случая дифференциальных уравнений с переменными коэффициентами для каждого участка может использоваться осредненная матрица Ai =A(xi) коэффициентов системы дифференциальных уравнений. Если рассматриваемый участок интервала интегрирования не мал, то предлагаются следующие итерационные (рекуррентные) формулы. Приведем формулы вычисления вектора частного решения, например, vin76.wmf на рассматриваемом участке vin77.wmf через векторы частного решения vin78.wmf, vin79.wmf, vin80.wmf соответствующих подучастков vin81.wmf, vin82.wmf, vin83.wmf. Имеем vin84.wmf. Также имеем формулу для отдельного подучастка:

vin85.wmf.

Можем записать

vin86.wmf,

vin87.wmf.

Подставим Y(x1) в Y(x2) и получим

vin88.wmf

vin89.wmf.

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

vin90.wmf

и получим, очевидно, что

vin91.wmf

и для частного вектора получаем формулу

vin92.wmf.

То есть вектора подучастков vin93.wmf не просто складываются друг с другом, а с участием матрицы Коши подучастка. Аналогично запишем vin94.wmf и подставим сюда формулу для Y(x2) и т.д.:

vin95.wmf

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

vin96.wmf,

очевидно, получаем, что

vin97.wmf

и вместе с этим получаем формулу для частного вектора:

vin98.wmf

То есть именно так и вычисляется частный вектор – вектор частного решения неоднородной системы дифференциальных уравнений, то есть так вычисляется, например, частный вектор vin99.wmf на рассматриваемом участке vin100.wmf через вычисленные частные вектора vin101.wmf, vin102.wmf, vin103.wmf соответствующих подучастков vin104.wmf, vin105.wmf, vin106.wmf.

Взято из [2]. Пусть дана система линейных алгебраических уравнений порядка n: vin107.wmf. Здесь над векторами (как в первоисточнике) поставим черточки вместо их обозначения жирным шрифтом. Будем рассматривать строки матрицы A системы как векторы:

vin108.wmf

Ортонормируем эту систему векторов. Первое уравнение системы vin109.wmf делим на vin110.wmf. При этом получим

vin111.wmf

vin112.wmf,

где

vin113.wmf, vin114.wmf, vin115.wmf

Второе уравнение системы заменяется на

vin116.wmf

vin117.wmf,

где

vin118.wmf, vin119.wmf,

vin120.wmf vin121.wmf

Аналогично поступаем дальше. Уравнение с номером i примет вид

vin122.wmf vin123.wmf,

где

vin124.wmf, vin125.wmf,

vin126.wmf

vin127.wmf

Процесс будет осуществим, если система линейных алгебраических уравнений линейно независима. В результате мы придем к новой системе vin128.wmf, где матрица C будет с ортонормированными строками, то есть обладает свойством vin129.wmf, где E – это единичная матрица.

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