книги / Прямые и обратные задачи механики упругих композитных пластин и оболочек вращения
..pdf82 |
Гл. 3. Методы решения краевых задач |
Отметим, что подпространство, базисом которого являются векторыстроки G p должно содержаться в подпространстве натянутом на вектор-столбцы Y. Противное означает либо некорректность краевой задачи, либо вырождение набора Y в результате численного интегриро вания. Итак, пусть существует такое Л, что выполняется (3.38), тогда для любого вектора с имеем
Y • с = (Y • Л -1 ) • (Л • с) = |
G jc " + Y 'c', |
(3.39) |
где Л • с = ||с"т |с/т ||т . Подставляя (3.39) |
в (3.37), получаем |
|
Gp • ( G j с" + Y 'c' + уо) = gp. |
(3.40) |
Так как по построению Y '_LG j и так как G p — ортонормированный набор векторов, выражение принимает вид
с" = gр - G py0. |
(3.41) |
При этом из (3.36) получаем у'0 = G j • gp. Таким образом определяют ся Уд и связь между с и с' и остается описать алгоритм построения Y ' и Л.
Построение Y ' будем вести способом, аналогичным применявшему ся для нахождения Y/ в случае двухточечной краевой задачи. Отличие будет заключаться в том, что вместо векторов е* будут использоваться вектор-столбцы у* матрицы Y и для того, чтобы построить матрицу преобразования Л, будут определяться коэффициенты в разложениях векторов на проекции и ортогональные дополнения.
Опишем соответствующий алгоритм по шагам: |
|
|
|||
1) |
задаются начальные значения г = 1, J = О |
|
|
||
2) |
определяются проекция y i и ортогональное дополнение |
у j |
век |
||
|
тора Уг |
на пространство, |
натянутое на вектор-столбцы |
G j, |
j = |
|
= 1, |
, SP матрицы G j и |
построенные векторы у'-, j = |
1, |
,J: |
|
|
Sp |
j |
|
|
|
Уг = |
5^(G j,yi) • Gj + |
^ ( y ',y i ) -у', Уг = Уг Уг’ |
(3-42) |
|
|
|
3 = 1 |
3 = 1 |
|
|
при этом элементы матрицы Л, находящиеся на пересечении г-го столбца и j -ой (j = 1, ... , Sp + J) строки равны соответствующим коэффициентам в разложении:
М, = |
(G j,y i), j |
Sp, |
AJ |
=(y'j-sp>yi)> |
j > S P |
3)если норма ||yi|| не равна нулю, то переходим к пункту 5,
4)полагаем элементы матрицы Л равными
=0 , j > Sp + J,
|
3.4. Алгоритм решения краевых задач для систем ОДУ |
83 |
|||||||
5) |
полагаем вектор y j +1 и элементы матрицы Л равными |
|
|||||||
|
|
|
|
|
y S + i |
= Уу <( 1 Г■1‘ . |
|
||
|
|
Л/+1 = 1, |
Ц = |
0. |
j > S p + J + 1, |
|
|||
|
увеличиваем г на единицу и переходим к пункту 6 |
|
|||||||
6) |
если г = |
S'p + |
1 и |
J |
= |
Sp — Sp + |
1, то процедура |
завершена |
|
|
успешно, |
если |
г = |
5 ' |
+ |
1 либо J |
= S'p — Sp + 1, то |
процедура |
завершена неуспешно, иначе переходим к пункту 2.
Здесь все суммирования ведутся от меньшего значения индекса к боль шему, т.е. если верхнее значение индекса суммирования меньше ниж него, то операция не выполняется. Отметим также, что при численной реализации, как и ранее, равенство нулю нормы вектора нужно заме нить приближенным.
В узле х р\р=р данный обобщенный |
алгоритм восстанавливает ре |
шение задачи, которое будет равно |
вектору уо, и если х р < х Г, |
то дальнейшее решение задачи эквивалентно решению задачи Коши. В процессе обратного хода при прохождении через точки хр восстанав ливается вектор с(с') и по соответствующим формулам вычисляется решение, аналогично алгоритмам методов дискретной ортогонализации и начальных параметров.
3.4.2. Обеспечение устойчивости расчетов. Для решени методом дискретной ортогонализации краевых задач, возникающих при определении напряженно-деформированного состояния оболочек вращения в рамках теории оболочек Кирхгофа-Лява либо теории типа Тимошенко, достаточно использовать равномерные сетки узлов ортогонализации и интервалов интегрирования. При этом типичное количество узлов ортогонализации, достаточное для устойчивого сче та, относительно невелико — до 20, а количество интервалов ин тегрирования, необходимое для достижение приемлемой относитель ной точности (порядка 10-5) составляет от 200 до 400. Однако при переходе к использованию уточненных теорий пластин и оболочек, учитывающих поперечные сдвиги, таких значений оказывается явно недостаточно и появляется необходимость в дополнительном контроле точности и устойчивости расчетов и в их обеспечении.
В методе дискретной ортогонализации, в дополнение к шагу инте грирования, появляется еще один управляющий параметр — расстояние между узлами ортогонализации. Совместно, эти параметры позволяют контролировать и управлять устойчивостью численного расчета, его точностью и объемом вычислений. Дополняя алгоритм метода дис кретной ортогонализации процедурами выбора шага интегрирования и расстояния между узлами ортогонализации, можно построить алго ритм решения краевых задач для систем ОДУ с малыми и большими параметрами, наподобие тех, что возникают при определении НДС пластин и оболочек в рамках уточненных теорий [9, 129].
84 |
Гл. 3. Методы решения краевых задач |
Проведем экспериментальную оценку значений расстояния между узлами ортогонализации и шага интегрирования, необходимых для обеспечения устойчивости расчета на соответствующих равномерных сетках. Для этого рассмотрим задачу изгиба слоистой длинной пря моугольной пластины, нагруженной равномерно распределенным дав лением. Обезразмеренная разрешающая система ОДУ, описывающая поведение такой пластины в уточненной постановке, имеет вид [9]
|
|
|
|
a ^ U " - ' y 2a2W ' " + а3П" = О, |
|
|
|
|
|||||
|
|
|
а2и"' - 7 2a 4W r" " + а5П'" = |
7 “ 2, |
|
|
(3.43) |
||||||
|
|
|
аъи " - |
72abW"' + а6П" - 7 ” 2а7П = |
0. |
|
|
||||||
Здесь |
|
|
|
к |
|
|
к |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
о\ |
= £ ^ 1 . ь |
а,2 = |
|
|
|
|
|
|
||
|
|
|
|
|
fc=i |
|
|
k=i |
|
|
|
|
|
|
|
|
К |
|
|
|
|
|
|
|
|
|
|
|
|
аз = |
|
|
+ 2(1 + |
(б4.к |
- |
|<S3.fc) ] |
, |
|
|||
|
|
|
к=1 |
|
|
|
|
|
|
|
|
|
|
|
|
К |
|
|
К |
|
|
|
|
|
|
|
|
|
а4 = ^ > Л |
, ь |
а5 = |
[тк6 2,к + 2(1 + Vi) |
(%,* - |
|
. |
||||||
|
|
|
|
|
к=1 |
|
|
|
|
|
|
|
|
06 = |
Ё - * |
[т|^1 ,к + |
|
|
(&4.к - |
+ |
|
|
|
|
|
||
|
к= 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
4(1 + |
Uk) |
( £ |
Q£ |
|
| 9 Г |
\ |
, (3.44) |
|
|
|
|
|
н---------2----- |
^7,fc ~~ Мб.к + ^o5ifcJ |
|||||||
|
|
|
a7 = V |
M i + i i ) № |
_ 2 St ,k + f „ ) , |
|
|
||||||
|
|
|
|
* |
|
.QL |
|
|
|
|
|
|
|
|
|
|
|
k= l |
|
|
|
|
|
|
|
|
|
Ti |
= 0 , T j |
= T j - t + |
2[(1 + n |
- l)/»j_ i - |
(1 + ^ ) / s 3] |
( $ _ , |
- |
| $ _ , ) . |
|||||
|
h k = ( 0 l - P L l)/i |
(i e N; к = |
1.......if; j |
= |
2........ if), |
||||||||
|
|
|
|
|
|
J |
|
|
|
|
|
|
|
|
|
|
7 = |
/i/Z , 0 j |
= ^ 2 h k /h {j = 0 , . . . , K ) , |
|
|
||||||
|
|
|
|
|
|
k= l |
|
|
|
|
|
|
|
|
|
sfc = |
E fc(l |
- |
i/£)“ |
|
= E k/ E {{k = |
1,..., ЛГ), |
|
где Л,, / — полная толщина и длина пластины; hk,Ek,Vk — толщина, модуль Юнга и коэффициент Пуассона /с-го слоя пластины; /Г — количество слоев; 1У, П, U — безразмерные обобщенные перемещения.
Выбирая в качестве искомой вектор-функции
у(ж) = \\W, W ' , W ”, W m |
, U , U' ||т , |
3.4. Алгоритм решения краевых задач для систем ОДУ |
85 |
представляем систему (3.43) в виде (3.11). Ненулевые коэффициенты матрицы А и вектора f будут при этом равны
/4 |
- 7 - 4 a i ( c3 + c j c - 1) - 1, |
А \2 - А 2з - А 34 - А 56 = А 78 = 1, |
||
|
А 4Ъ = a7ci с2 164, А м = |
—72CI |
1, |
Л65 = 7 _2ai a7c^1, |
|
As4 = 7 2a, '(аг + a3Cic2 '), |
A8s = 7 _2аза7с2 \ |
||
где ci |
= <2 2 (1 3 —aids, c2 = aids — a2, C3 = a2 |
—a ja 4. Разрешающая си |
стема уравнений является системой с постоянными коэффициентами. Шесть собственных значений матрицы системы Ai = А2 = A3 = А4 = = А5 = Аб = 0, оставшиеся два А7 = —Ав = А, где А = у/А&Г-^А^Жм е GM [9]. При этом, величина А пропорциональна отношению l/h.
С учетом введенных обозначений общее решение системы (3.43) имеет вид
W{x) = Ci + С2х + Сзх2 + С4х 3 + С5еХх + С6е~Хх + С х \
Щ х) = - 6 ^ С |
4 + А ^ Х \ С ъеХх - С%е~Хх) + Ьх, |
|
Щх) = CS + С7х + 3(А84 - |
A 85j ^ ) C 4x 2+ |
|
|
^65 |
|
+ ± (24А84С + А 8ЬЪ)хъ + А(А84 + ^ ) { С ъеХх - |
С6е~Хх), (3.45) |
|
О |
А 4 6 |
|
где b = (24С - f 4)/ A 46, С |
= f 4A 88/(24X2), С \ , . . . , С 8 |
— произвольные |
постоянные. Полученные аналитические решения, наряду с полиноми альными, содержат экспоненциальные компоненты вида ехр(А(ж — 1)) (кривая 1 на рис. 3.2, а), ехр(—Хх) (кривая 2 на рис. 3.2, а) (в случае, приведенном на рис. 3.2, а, А = 155,66).
Рис. 3.2
Жесткое защемление пластины соответствует граничным условиям вида
[W, W', П, U]T {0) = [W, W', П, С/]т (1) = 0. |
(3.46) |
86 |
Гл. 3. Методы решения краевых задач |
|||
В этом случае произвольные постоянные будут равны |
||||
|
Ci = -(С 5 + С6), |
|
С2 = -А (С 5 - С6), |
|
Сз = (1 —еА+ А —В)С$ + (1 —е~х |
- А + В)Се — С, С4 = В{С5 - С е ), |
|||
C5 = - ( 2 C + D6C6) / D 5, |
С6 _ ЪАыРб/{ЬАм)-2В(ех - 1 ) С |
|||
|
|
|
|
" B [ ( e x - l ) D 6 + ( e - x - l ) D 5] ’ |
С7 = - |
Лб4\ |
, |
24AsiC + Asbb |
|
с 8 + з с 4 ( Ч 4 - а 85^ |
) |
+ |
+ |
|
|
+ |
А (л 84 + ^ ) ( С 5еЛ - С 6е - Л) , (3.47) |
с 8 = —А(А84 + А 85/А46)(С5 — Се),
где D5 = 2 + А - (2 - А)ел + В, D6 = 2 - А - (2 + А)е"А - В, В =
= А3А б 5 /(6 -А 4б44б4).
Рассмотрим трехслойную пластину, внешние слои которой имеют толщину Л,] и выполнены из материала с модулем Юнга Е\ и коэффи циентом Пуассона v\. Внутренний слой толщины h2 выполнен из мате риала с модулем Юнга Е 2 и коэффициентом Пуассона v2. Из рис. 3.2, б видно, что спектральный радиус матрицы системы А существенно зави сит от механически характеристик материалов и структуры пластины и эта зависимость нелинейная.
о
О 0,25 0,5 0,75 х
Рис. 3.3
На рис. 3.3 приведен вид решения задачи изгиба жестко защем ленной на обоих краях трехслойной пластины в разрешающих функ циях, нормированных в равномерной метрике (геометрические и ме
ханические параметры пластины: |
l/h |
= 15, h2/hi = 8 , Е \ / Е 2 = 10, |
v\ = v2 = 0,3). Сплошные линии |
на |
рис. 3.3 соответствуют W*(x), |
U*{x); штриховые — Wl(x), С/'(ж); пунктирные — W"(x), П*(ж); штрихпунктирные — W " ', П*. Видно, что получаемые решения имеют ярко выраженные краевые эффекты.
Таким образом, данная задача является очень удобным инструмен том для численного исследования проблемы обеспечения устойчивости
3.4. Алгоритм решения краевых задач для систем ОДУ |
87 |
иточности вычислений, так как она описывается системой уравнений с постоянными коэффициентами, спектральный радиус матрицы системы которой много больше единицы и легко регулируется с помощью пара метров задачи. При этом полученное аналитическое решение позволяет вычислять как глобальные, так и локальные погрешности вычислений.
Исследовалось влияние расстояния между узлами ортогонализации
иколичества интервалов интегрирования между ними на точность расчетов, которая использовалась как критерий их устойчивости. Про веденные численные эксперименты показали, что с увеличением жест кости задачи начинает проявляться проблема размазанности области перехода от неустойчивых расчетов к устойчивым при уменьшении расстояния между узлами ортогонализации. При этом данная проблема проявляется сильнее при большем числе интервалов интегрирования между соседними узлами ортогонализации.
а |
|
б |
|
|
1 |
|
|
1 |
щ |
\ |
Ч |
|
1ЦЦ.1 * |
ч |
|
ч ; |
|
|
|
•i |
70 |
85 |
100 |
115 |
130 |
810 |
100 |
160 |
520 |
|
|
|
|
Рис. 3.4 |
|
|
|
|
Так, на рис. 3.4 представлены зависимости погрешности вычислений (логарифмированной по основанию 10) от числа узлов ортогонализации J для различного числа интервалов интегрирующей процедуры: штри ховые кривые соответствуют трем интервалам, пунктирные — четырем, сплошные — пяти. Зависимости получены для описанной выше задачи
при значении l / h = 100 |
(рис. 3.4, а) и |
l/h = |
500 (рис. 3.4, б), спек |
тральные радиусы матриц |
системы при |
этом |
равняются Л = 1037,73 |
и 5188,63 соответственно. Видно, что при меньшем спектральном ради усе прослеживается четкая граница перехода от неустойчивых решений к устанавливающимся, т. е. таким, погрешность которых стабильно уменьшается при уменьшении шага сетки. При больших спектраль ных радиусах очень существенной становится проблема сдерживания экспоненциального роста погрешности численного интегрирования, что проявляется в возникновении осциляций погрешности на графиках, соответствующих большему числу интервалов интегрирования между ортогонализациями. Отметим, что при значении l / h = 1000 осциляции возникают уже при использовании трех интервалов интегрирования.
88 |
Гл. 3. Методы решения краевых задач |
Таким образом, чем больше Л уравнений, тем меньшее количество интегрирований допускается между узлами ортогонализации.
Рассмотрим теперь, какое расстояние между узлами ортогонализа ции необходимо для устойчивого расчета. Для этого построим зави симость числа узлов ортогонализации J, при котором расчет является устойчивым (по признаку стабильно небольшой погрешности) от спек трального радиуса системы Л.
6
Рис. 3.5
На рис. 3.5, а представлены два варианта таких зависимостей: при нефиксированном (больше трех, подбираемом в зависимости от ре зультатов) количестве интегрирований между узлами ортогонализа ции (пунктирная кривая) и при фиксированном (три интегрирова ния, сплошная кривая). Во втором случае фиксировался и критерий устойчивости — достижение погрешности вычислений порядка 10_3. Видно, что использование числа интервалов интегрирования больше трех целесообразно при относительно небольших Л, иначе это приводит к нестабильности в поведении метода, о чем сказано выше, а выиг рыш в количестве ортогонализаций при этом составляет < 3 0 % для А = 5000 и уменьшается с ростом А. В обоих случаях зависимость приближенно стабилизируется на некоторой линейной функции при увеличении А. Однако для фиксированного количества интегрирований между узлами ортогонализации эта зависимость более стабильная: J « 0,095 • А + 25.
При описании метода дискретной ортогонализации [77] С. К. Году нов указывает в качестве признака устойчивости счета, кроме есте ственного условия невырожденности набора векторов Y во всех точках интервала решения, близость к единице определителей матриц W j. Эти матрицы имеют особую структуру — они верхнетреугольные; та ким образом, их определители равны произведению их диагональных элементов, которые в свою очередь являются нормирующими множите лями для векторов-решений. Последнее и накладывает ограничения на
3.4. Алгоритм решения краевых задач для систем ОДУ |
89 |
их величины — нормирование не должно приводить к существенным потерям в значащих цифрах компонент векторов-решений.
Рассмотрим поведение двух различных величин, характеризующих
матрицы W j\ Wm = max{—1пю \Wu j|} |
и Wd = max{—1пю | det W j|}. |
|
ij |
' |
j |
На рис. 3.5,6 приведены значения |
Wm , полученные при построении |
представленных на рис. 3.5, а зависимостей. Видно, что близость к еди нице нормирующих множителей весьма относительна — они могут отличаться от нее на 2-8 порядков, хотя расчет остается устойчивым. При увеличении спектрального радиуса матрицы системы рассмат риваемые величины стабилизируются, гиперболически приближаясь к значениям порядка 4-5. Отметим, что для задач с относительно небольшим спектральным радиусом матрицы Якоби системы уравне ний, устойчивость счета можно обеспечить при меньшем числе узлов ортогонализации, увеличивая число интервалов интегрирования между ними.
Посмотрим теперь, как меняются матрицы W j при изменении числа узлов ортогонализации. В табл. 3.1 представлены данные для задач с
J-1/h |
1 |
2 |
4 |
l/h |
|
|
Wm |
1 0 |
4,5 |
2,25 |
1,13 |
50 |
4,5 |
2,25 |
1,58 |
250 |
4,5 |
2,5 |
2,28 |
1 0 0 0 |
4,28 |
3,1 |
2 . 8 8 |
Т а б л и ц а 3.1
|
1 0 |
1 |
2 |
4 |
|
1 0 |
|
|
|
|
w d |
|
|
5,33 |
■10- 1 |
4,64 |
2,32 |
1,16 |
4,63 |
■10“ ' |
8,83 |
■10“ ' |
4,53 |
2,27 |
1,13 |
4,53 |
■10“ ' |
0 |
1 |
4,51 |
2,26 |
1,13 |
4,51 |
• 10_ 1 |
сл0 |
о |
|||||
|
— |
4,28 |
2,24 |
1,13 |
|
— |
разным соотношением l/h. В качестве второго параметра рассматрива ется величина J - l/h, что позволяет в определенной степени оставаться независимыми от жесткости задачи при исследовании величин Wm и Wd. Действительно, если посмотреть на правую половину табл. 3.1, то видно, что Wd мало изменяется по вертикали и при этом обратно пропорционален предложенному параметру. Для Wm наблюдаются схо жие тенденции. Это позволяет при реализации в алгоритме выбирать расстояние до следующего узла ортогонализации исходя из предполо жения линейной зависимости от него величин Wd, W m.
Резюмируя проведенные исследования можно сказать, что для задач с небольшими спектральными радиусами матрицы системы диффе ренциальных уравнений большой вклад в обеспечении устойчивости расчета имеет вопрос устойчивого численного интегрирования, обеспе чиваемого выбором шага. С ростом спектрального радиуса матрицы си стемы на первый план выходит необходимость сдерживания взрывного роста погрешности интегрирования, что обеспечивается уменьшением расстояния между узлами ортогонализации.
90 |
Гл. 3. Методы решения краевых задач |
|
|
Таким образом, в качестве критерия устойчивости счета на этапе |
|
ортогонализации выбираем |
|
|
|
W m < const |
(3.48) |
либо |
(3.49) |
|
|
Wd < const, |
где константа зависит от таких факторов, как разрядность математи ческих вычислений. Далее в расчетах const « 4,0. При выборе рассто яния до следующего узла ортогонализации будем руководствоваться полученной величиной W m или Wd и предположением о линейной зависимости между выбранной для контроля величиной (Wm или Wd) и расстоянием между узлами ортогонализации.
Другой фактор обеспечения устойчивости счета на прямом ходе метода дискретной ортогонализации — устойчивость численного инте грирования задач Коши. Напомним, что одной из главных особенно стей системы ОДУ при решении краевых задач является возможность существования и больших положительных, и больших отрицательных собственных значений в ее матрице Якоби, как, например, для рас сматриваемой тестовой задачи. Также нужно отметить, что, прежде чем провести ортогонализацию и оценить ее успешность, необходимо численно проинтегрировать задачи Коши, не допустив при этом взрыва погрешности и обеспечив устойчивость расчета.
В работе [256] для контроля и обеспечения устойчивого числен ного интегрирования задач Коши для жестких систем ОДУ методом Рунге-Кутты-Мерсона предлагается при выборе шага интегрирования руководствоваться следующим критерием
12 |
( k 3 - к |
2 ) г |
|
(3.50) |
— max |
1, |
|||
7 » |
(k2- k,)i |
^ |
|
где kj — стадии метода. Однако прямое применение этого критерия приводит к неоправданному измельчению сетки интегрирующей проце дуры вслед за узлом, в котором производится ортогонализация. Экс перименты на тестовых задачах при использовании равномерных сеток показали, что на интервале, следующем за узлом ортогонализации,
происходит скачок величины 1СТ = 12/7 • max |(кз —кг^Д кг —ki)*|
г
на несколько порядков. На последующих интервалах интегрирования рассматриваемая величина нормализуется, при этом для устойчивых расчетов критерий (3.50) начинает выполняться, а в случаях, когда он не выполняется на нескольких интервалах интегрирования подряд, результаты расчета неверны. Таким образом в качестве условия устой чивости на этапе решения задач Коши выбираем выполнение критерия (3.50) начиная со второго либо с третьего от узла ортогонализации интервала интегрирования.
На основе сформулированных критериев были реализованы и до бавлены в алгоритм метода дискретной ортогонализации процедуры
3.4. Алгоритм решения краевых задач для систем ОДУ |
91 |
контроля и управления устойчивостью расчета выбором шага интегри рующей процедуры и расстояния между узлами ортогонализации.
3.4.3. Схема алгоритма. Принципиально схему алгоритма можно разбить на пять блоков:
1)обработка входных данных и вычисление начальных векторов на левом крае интервала решения на основе граничных условий;
2)прямой ход метода: решение задач Коши и ортогонализации на бора векторов решений с автоматическим подбором шагов соот ветствующих сеток;
3)разрешение граничных условий на правом крае и вычисление век тора свободных коэффициентов для последнего интервала между узлами ортогонализации;
4)обратный ход метода: вычисление векторов свободных коэффици ентов на всех интервалах решения и восстановление там прибли женного решения краевой задачи;
5)подготовка выходных данных.
Вкачестве исходных данных алгоритму требуются:
—порядок системы уравнений;
—левая и правая граница интервала решения;
—количество граничных условий на левом и правом крае;
—функция, вычисляющая матрицу Якоби и вектор свободных эле ментов системы;
—функция, вычисляющая матрицы и векторы граничных условий;
—ряд вещественных параметров:
—константа, определяющая жесткость критерия (3.48),
—максимально и минимально допустимые расстояния между узлами ортогонализации,
—максимально и минимально допустимый шаг интегрирова ния.
На выходе алгоритм выдает следующую информацию о решении:
—узлы сетки интегрирующей процедуры,
—узлы ортогонализации,
—набор векторов — приближенное решение краевой задачи во всех узлах сетки интегрирующей процедуры.
На основе полученных данных можно также восстановить в тех же узлах векторы производных решения, которые можно использовать для построения интерполяционных сплайнов. Кроме того, в результате работы алгоритма на прямом ходе возникает служебная информация о структуре решения: наборы векторов-решений задач Коши и линейные операторы ортогонализации в соответствующих узлах, которые могут быть использованы для проведения уточняющих расчетов.
На первом этапе алгоритма на основе граничных условий на ле вом крае вычисляем начальные векторы по алгоритму, описанному в п. 3.4.1. Определяем исходные шаг интегрирования и расстояние до