Файл: Васильев, С. П. Приближенные методы расчета сооружений на устойчивость и динамику с применением ЭВМ Наири (учебное пособие).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—^212 - — ; Г23—г320;

Взят из книги [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

 

60,1 a

 

 

 

2

 

8

 

 

0

 

 

 

6

=

0 .

( б )

0,53a

 

0

 

204,8—4,85a

0

 

 

 

 

 

60,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