Полученное решение на каждом шаге нагружения можно уточнить,
используя алгоритм итераций (уточняя приращение вектора кривизны
при каждой итерации). При итерациях на каждом шаге нагружения
получаем уравнения равновесия
d
Y
(
j
)(
k
)
dη
+ A
(
j
)(
k
)
Y
(
j
)(
k
)
= b
(
j
)
,
(6)
где индекс
k
— номер итерации;
Y
(
j
)(
k
)
= ΔQ
(
j
)(
k
)
ΔM
(
j
)(
k
)
θ
(
j
)(
k
)
u
(
j
)(
k
)
T
;
A
(
j
)(
k
)
=
A
(
j
)(
k
)
æ
A
−
1
A
(
j
)(
k
)
Q
A
(
j
)(
k
)
L
0
A
1
A
(
j
)(
k
)
æ
+
A
−
1
A
(
j
)(
k
)
M
0
0
0
−
A
−
1
A
(
j
)(
k
)
æ
0
0
0
A
1
A
(
j
)(
k
)
æ
;
A
(
j
)(
k
)
æ
=
=
0
−
æ
(
j
−
1)(
k
j
−
1
)
3
−
Δæ
(
j
)(
k
−
1)
3
æ
(
j
−
1)(
k
j
−
1
)
2
+Δæ
(
j
)(
k
−
1)
2
æ
(
j
−
1)(
k
j
−
1
)
3
+Δæ
(
j
)(
k
−
1)
3
0
−
æ
(
j
−
1)(
k
j
−
1
)
1
−
Δæ
(
j
)(
k
−
1)
1
−
æ
(
j
−
1)(
k
j
−
1
)
2
−
Δæ
(
j
)(
k
−
1)
2
æ
(
j
−
1)(
k
j
−
1
)
1
+Δæ
(
j
)(
k
−
1)
1
0
;
A
(
j
)(
k
)
Q
=
0
Q
(
j
−
1)(
k
j
−
1
)
3
−
Q
(
j
−
1)(
k
j
−
1
)
2
−
Q
(
j
−
1)(
k
j
−
1
)
3
0
Q
(
j
−
1)(
k
j
−
1
)
1
Q
(
j
−
1)(
k
j
−
1
)
2
−
Q
(
j
−
1)(
k
j
−
1
)
1
0
;
A
(
j
)(
k
)
L
=
0
l
(
j
−
1)(
k
j
−
1
)
32
−
l
(
j
−
1)(
k
j
−
1
)
22
−
l
(
j
−
1)(
k
j
−
1
)
32
0
l
(
j
−
1)(
k
j
−
1
)
12
l
(
j
−
1)(
k
j
−
1
)
22
−
l
(
j
−
1)(
k
j
−
1
)
12
0
;
A
(
j
)(
k
)
M
=
=
0
M
(
j
−
1)(
k
j
−
1
)
3
+Δ
M
(
j
)(
k
−
1)
3
−
M
(
j
−
1)(
k
j
−
1
)
2
−
Δ
M
(
j
)(
k
−
1)
2
−
M
(
j
−
1)(
k
j
−
1
)
3
−
Δ
M
(
j
)(
k
−
1)
3
0
M
(
j
−
1)(
k
j
−
1
)
1
+Δ
M
(
j
)(
k
−
1)
1
M
(
j
−
1)(
k
j
−
1
)
2
+Δ
M
(
j
)(
k
−
1)
2
−
M
(
j
−
1)(
k
j
−
1
)
1
−
Δ
M
(
j
)(
k
−
1)
1
0
.
16 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2008. № 2
