ВУЗ:
Составители:
58
ственную функцию системы, нам нужны только относительные значения
компонент вектора Z. Поэтому можно считать, что одна из неизвестных
(скажем, z
1
) всегда равна единице.
2) Вычислить АZ
g1
.
3) Произведение АZ
g1
представляет собой вектор, который можно за-
писать в виде λ
g2
Z
g2
, где λ
g2
– множитель, такой, что компонента z
1
вектора
Z
g2
опять равна единице, а остальные переменные z
2
, z
3
, …, z
n
принимают
соответствующие значения.
4) Сравнить Z
g2
с Z
g1
или в общем случае Z
gr
с Z
g(r+1).
Если они не отли-
чаются (в пределах заданной точности) друг от друга, то полученное мно-
жёство значений образует собственный вектор, а множитель представляет
собой наибольшее собственное значение. В противном случае снова вер-
нуться к пункту 4.
В результате наибольшим собственным значением становится после-
дующее значение λ. После этого процесс итераций повторяется.
Условием окончания итерационного процесса являются следующие
данные: максимальное число итераций и точность определения собствен-
ного числа (малая положительная величина). Варьируя в исходных дан-
ных значениями максимального числа итераций и точностью определения
собственного числа можно получить желаемый результат.
Содержание основных блоков алгоритма МКЭ последовательность их
выполнения в задаче устойчивости состоит в следующем.
1. Образование расчётной схемы.
2. Вычисление матриц [K
r
] жёсткости КЭ в местной системе коорди-
нат XYZ.
3. Вычисление матриц [G
r
] потенциала нагрузки КЭ в местной систе-
ме координат.
4. Вычисление матриц [T
r
] ортогонального преобразования КЭ.
5. Вычисление матриц [K
0
r
] жёсткости КЭ в общей системе координат
X
0
Y
0
Z
0
:
[K
0
r
] = [T
r
]
Т
[K
r
] [T
r
].
Страницы
- « первая
- ‹ предыдущая
- …
- 56
- 57
- 58
- 59
- 60
- …
- следующая ›
- последняя »
