Файл: Васильев, С. П. Приближенные методы расчета сооружений на устойчивость и динамику с применением ЭВМ Наири (учебное пособие).pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 02.11.2024
Просмотров: 31
Скачиваний: 1
III. ПРИМЕНЕНИЕ МЕТОДА «ОБОБЩЕННЫХ КООРДИНАТ»
ДЛЯ РАСЧЕТОВ НА УСТОЙЧИВОСТЬ И ДИНАМИКУ
Для расчетов сложных рам на устойчивость и динамику вышеизложенный метод, являющийся совокупностью метода «малых возмущений» с методом сил, может оказаться весьма трудоемким в смысле формирования исходных матриц. Для расчета сложных многоконтурных рам более приемлем метод «обобщенных координат», предложенный В. В. Болотиным [5], в основе которого лежит метод перемещений.
Метод «обобщенных координат» заключается в следующем. Деформированный вид рамы в момент потери устойчивости или колебаний описывается рядом «обобщенных координат», то есть выбираются такие аппроксимирующие функции, кото рые достаточно характеризуют формы изгиба элементов рамы. В качестве таких функций принимаются уравнения, описыва ющие формы изгиба элементов рамы в основной системе ме тода перемещений. Пусть п0— число неизвестных при расчете рамы по методу перемещений (степень кинематической неоп ределимости). Тогда в качестве первых щ аппроксимирующих функций можно принять формы статического изгиба стержней
рамы от соответствующих единичных воздействий |
(рис. 10, 6). |
Остальные функции выбираются при Z\ = Z2 — |
Z nQ=0, |
причем их число определяется необходимой точносшо расче-
ча. Так, например, в раме |
при расчете |
на |
устойчивость |
|
(рис. 10, а) для сжатых стержней 1—2 и 3—4 введение |
обоб |
|||
щенных координат Z6 крайне |
желательно |
(рис. |
10, в), |
в то |
время как уточняющее влияние координат Z4 и Z5 малосуще ственно и их можно не учитывать.
Количество обобщенных координат (Zb Z2, ..., Z„) дефор мированного состояния рамы и определяет число неизвестных в системе уравнений, составляемых при расчетах на устойчи вость и динамику по изложенной методике. Заметим, что чис ло этих координат п бывает обычно больше числа неизвест ных метода перемещений п0.
В работе Г5] произведен расчет элементарных стержней на воздействие единичных значений обобщенных координат и ре зультаты приводятся в виде таблиц. При расчетах на устойчи вость необходимо пользоваться табл. 2, а на динамику — табл. 3. Существенная особенность этих таблиц состоит в том, что параметры, определяющие критические нагрузки и собст венные частоты, входят в.реактивные усилия в явном виде,
40
Рис. 10
Т а б л и ц а 2
ю
Т а б л и ц а 3
Это приводит к тому, что отпадает надобность в решений трансцендентных уравнений, получающихся после раскрытия определителей. Данные табл. 2 и 3 можно получить и путем преобразования формул, реактивных усилий «классического» метода перемещений [1 ], разлагая их трансцендентные функ ции в ряды, удерживая при этом первые члены ряда. Это, ко нечно, снижает точность расчетов и делает невозможным по лучение полного спектра критических значений параметров. Однако, как было отмечено ранее, в инженерных расчетах не обходимо обычно знать лишь низшие члены спектра, которые по указанным методикам вычисляются достаточно точно.
Практическое применение метода «обобщенных коордиг мат» рассмотрим на примерах.
§ 6. Примеры расчета рам на устойчивость
Пример 1*. Подсчитать P Kp(mm) для рамы, показанной на рис. 11, а. Основная система и обобщенные координаты, кото рые приняты для описания деформированного вида этой рамы,
приведены на рис. |
11,6. На рис. 11, в |
показаны деформиро |
||||||
ванные виды основной системы |
от единичных Z u Z2, Z3 и Z4, |
|||||||
а на рис. |
11 , г — соответствующие эпюры изгибающих момен |
|||||||
тов, построенные |
по данным табл. 2. |
Реактивные моменты |
||||||
и реактивные |
сосредоточенные |
силы |
ri} определяются как в |
|||||
методе перемещений на основании данных табл. 2. |
||||||||
ЛЕГ |
t |
4 £ / |
4ЕГ |
2 |
Р1_ |
11-0,133 |
||
I |
Л~ |
Г |
|
I |
~ |
15 |
||
|
|
|
||||||
4EI |
|
ЛЕГ |
= 8- |
Е1_ . |
|
|
||
г,, = ---- ++-----— |
|
|
||||||
|
|
|
|
|
/ |
’■ |
|
|
24-1,2
^12—^21— 2 - — ; Г23—г32—0;
Взят из книги [51.
44
м . L |
s n |
N
|
6EI |
Р_ |
|
|
№_\ Ш_ _ |
Гн=г41= — - |
10 |
6- |
0,1 |
EI ) /2 ’ |
|
|
|
|
|
||
_ |
_л £/ |
|
|
|
|
' 24 —/ 42 —O -7J-; |
|
|
|
|
/■34==-/'43==0.
Необходимая для определения Ркр система уравнений
принимает вид |
|
|
|
|
|
|
|
|
|
|
|
||
(11-0,133а) - - |
Z,+2 у-,%+0,53а ^ Z 3+ (6-0,1 a) - ^ Z 4=0 |
||||||||||||
|
|
FI |
|
FI |
|
|
|
EI |
|
|
|
||
|
2 - Z j + 8 — |
Z2+ 0 + 6 — Z 4= 0 |
|
|
|
||||||||
|
|
/ |
|
/ |
|
|
|
|
/2 |
|
|
(а) |
|
0,53a^Z,+04-(204,8-4,85a) ^ - Z 3+0 = 0 |
|||||||||||||
| |
|||||||||||||
(6-0,la) ^ - Z 1+ 6 ^ Z 2+04-(24-l,2a) ^ |
Z4= 0. |
||||||||||||
|
|
|
p/2 |
|
|
|
|
|
|
|
/ |
||
Здесь принято a = |
|
|
|
|
|
|
имеет не нуле- |
||||||
— . Эта система уравнений |
|||||||||||||
вое решение при |
EI |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||||
11-0,133a |
2 |
|
0,53a |
|
6—0,1 a |
|
|
|
|||||
2 |
|
8 |
|
|
0 |
|
|
|
6 |
= |
0 . |
( б ) |
|
0,53a |
|
0 |
|
204,8—4,85a |
0 |
||||||||
|
|
|
|
|
|||||||||
6—0,1a 6 |
|
|
0 |
|
|
24—1,2a |
|
|
|
||||
Для этого |
уравнения |
am!n =15,33, откуда Якр =15,33 — . |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
У |
[2 |
|
Перепишем уравнения (а) |
в матричной форме |
|
|||||||||||
где |
|
|
|
|
(/? -aS )Z = 0 , |
|
|
(19) |
|||||
|
|
|
|
r |
* |
|
|
|
|
|
|
||
|
|
I |
12 ... |
ln |
|
S n |
S 12 |
|
S i n |
||||
/ |
11 |
r |
|
|
|
||||||||
**И* |
r *22 |
|
Г |
2n |
|
S 21 |
S 22 |
■■ |
S o n |
||||
Г |
21 |
|
, |
||||||||||
R = |
|
|
|
|
|
|
|
5 = |
|
|
|
||
|
|
r * |
|
. . |
Г * |
|
|
|
S n2 |
|
|
||
r |
n 1 |
' |
n2 |
nn\ |
|
S n i |
•. : < n |
||||||
**• |
|
|
|
||||||||||
|
|
|
|
|
|
|
I |
• |
|
|
|
|
|
|
|
|
|
Z = i |
|
|
|
|
z я
46
Здесь R — матрица единичных реакции от обобщенных пере мещений Zi= I, Z 2 = 1, Z„ = 1; S — матрица поправок к реакциям в связях по стержням, подвергаемым сжатию; Z матрица-столбец обобщенных координат.
В пашем случае
11 |
2 |
0 |
6 |
|
|
0,133 |
0 |
-0,53 |
0,1 |
2 |
8 |
0 |
6 |
|
с _ |
0 |
0 |
0 |
0 |
0 |
0 |
204,8 |
0 |
; |
Л — |
-0,53 |
0 |
4,85 |
0 |
6 |
6 |
0 |
24 |
|
|
0,1 |
0 |
0 |
1,2 |
-Умножив равенство (19) слева на матрицу R~l и разделив на
(—а),'получим |
(R~lS — —■Е) Z= 0, |
(20) |
|
У |
|
которое представляет собой вековое уравнение
\А -к Е \= 0 ,
где A ^ R - 'S , ),= “ ..Чтобы получить amin, необходимо найти
Унах для матрицы А. Таким образом, решение задачи на ус тойчивость рам приведено к формированию матриц R и S, а затем — к отысканию лшвх для матрицы A=R~'S по програм ме, изложенной в Приложении. Тогда
Я к к и и с - г ^ - . |
(21) |
Лт?х * |
|
Следует отметить, что для рам, имеющих линейную по движность узлов, нет надобности вводить обобщенные коор динаты, учитывающие изгибные деформации стоек (в нашем случае Z3) . Действительно, учитывая только три обобщенных координаты Zb Z2 и Z4, получим определитель из (б) путем вычеркивания третьего столбца и третьей строки, то есть
1 1 — 0 , 1 З З о с 2 |
6 — 0 , 1 а |
- 0 . |
|
|
||
2 |
8 |
0 |
|
|
|
|
6 — 0 , 1 а |
6 |
2 4 — 1 , 2 а |
|
|
|
|
Решив это уравнение, |
получим amin =15,38, что отличается от |
|||||
предыдущего значения на 0,5%. |
|
трехэтажной |
рамы |
|||
Пример 2*. Для |
пятипролетной |
|||||
(рис. 12, а) требуется |
определить |
параметр |
/ \ т(min). |
Причем |
||
у этой рамы жесткость стоек EI, |
а жесткость ригелей 2ЕЕ |
|||||
Так как рама и нагрузка симметричны, |
то |
необходимо рас |
смотреть симметричную и обратно симметричную форму по тери устойчивости.
* Данные для примера и результаты взяты из [81.
47