Файл: Васильев, С. П. Приближенные методы расчета сооружений на устойчивость и динамику с применением ЭВМ Наири (учебное пособие).pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 02.11.2024

Просмотров: 33

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Построим приближенную эпюру М р в данной раме (рис. 3, г), помня, что в сжатом стержне в момент потери устойчивости эпюра будет криволинейна. Охарактеризуем изображенную форму равновесия (рис. 3,в) двумя ординатами Ai и Аг, а эпюру Л1р в пределах стойки будем интерполировать квадратной параболой. Разбив раму на расчетные участки и сечения (рис. 3,<9), сформируем матрицы, необходимые для получения матрицы

C= L'BA.

Рис. 3

21

Для формирования матрицы податливости В запишем

В ,

В =

В 11

в.

\_

go

l_

 

 

 

 

2

1

6

4gi

6

 

 

 

6

1

2

 

g*

 

 

 

 

 

 

 

 

 

 

 

 

Делая сдвижку по диагонали,

ввиду

непрерывности эпюры

б сечении 2, имеем

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

8

4

Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

Запишем выражение

в сечениях

0,

1, 2, 3

через

реакции

(рис. 3, в ) .

 

М„ = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mi=Ai

 

 

 

 

 

 

 

 

М% —Д2

 

 

 

 

 

 

 

 

Af3=0.

 

 

 

 

 

Тогда матрица А примет вид

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0

 

 

 

 

 

 

 

 

0

1

 

 

 

 

 

 

 

 

0

0

 

 

 

 

Для формирования матрицы L

построим

эпюры М { и М2

(рис. 3, е)

от единичных сил, приложенных

по направлению

перемещений Ai и Д2Тогда.

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

0,5

0,5

 

 

 

 

 

 

 

0,5

1

 

 

 

 

 

 

 

0

 

0

 

 

 

 

Проделав операции над матрицами L, В и А по программе,

приведенной-в Приложении, получим

 

 

 

 

 

 

 

^тах '

6,83 ,

 

 

 

 

 

 

 

v =

1

 

 

 

 

 

 

 

1,41

 

 

 

 

 

 

 

 

 

 

 

 

22


соответствующая критическая сила с учетом ранее вынесемпых скаляров Р, / и E Iр

 

6£ /„

ЕТ„

■' к р ы п ' п )

/2)

0,878 2 -

1

1 ЛП1ЯХ

I

§ 4. Устойчивость статически неопределимых стержней и рам

При расчете на устойчивость статически неопределимых стержней можно, как и при расчете статически определимых систем, свести расчет к определению /-тах для матрицы

C* — L'*BA*,

где L*, А* — матрицы, формируемые для статически неоп­ ределимой системы. Однако формирование матриц L* и Л* для этих систем настолько трудоемко, что приходится искать иные пути решения задач устойчивости статически неопреде­ лимых систем.

При расчете подобных систем на устойчивость оказывает­ ся удобным метод малых возмущений, рассмотренный выше, в комбинации с методом сил. Это позволяет формировать ис­

ходные

матрицы

для

 

статически

определи­

и)

мой системы. Хотя при

 

этом

увеличивается

 

объем вычислений,

но

 

эта трудность устраня­

 

ется при помощи ЭВМ.

 

Рассмотрим

 

сущ­

 

ность этой методики на

 

простейшем

примере

 

(рис. 4,а). Деформиро­

 

ванная

ось

стержня в

г)

момент

потери

устой­

 

чивости

показана

на

 

рис. 4, б.

 

 

 

 

 

В этом случае

эпю­

 

ра изгибающих

момен­

 

тов примет

вид,

пред­

 

ставленный на рис. 4, г.

зп

Ординаты этой

эпюры

 

неизвестны,

однако ее

Рис. 4

26


можно представить как сумму двух эпюр

(рис.

4, д), возни­

кающих в основной системе (рис.,4, в ) :

 

 

 

М = М / + М р°.

 

 

(а)

Будем полагать, что равенства (а) записаны для ряда, на­

меченных сечений. Тогда можно говорить,

что

в

этом ра­

венстве:

 

 

 

уИу1— матрица изгибающих моментов

в основной системе

от неизвестных метода сил (изгибающие

моменты

вычисляются в недеформированном состоянии);

Мр° — матрица изгибающих моментов,

вычисляемая в ос­

новной системе в деформированном состоянии.

В свою очередь, как было показано в расчете рам методом сил [4],‘матрицу М х° можно представить в виде

M x° ~ L xX .

Матрицу М р можно записать как

Мр = А$ А.

Тогда (а) запишется в виде

Л 1 -т Л.А'+л0.\

или

M = \\LxAq\\

[ .

(б)

После того, как определена эта матрица-столбец, можно перейти к определению перемещений в данной конструкции.

Во-первых, перемещения по направлению неизвестных X t должны равняться нулю

bx = L x'B M = 0 .

(в)

Подставляя в зависимость (в) матрицу М и вынося за матрицы произведение характерных размеров и нагрузок (Р, I и £7), получим

^ . ^ ^ \ \ L / B H \ L X A 0\\

(г)

E l

Во-вторых, выбранные перемещения А определяются как

A = L/BM

или

А

РЕ \\l i b \v \\l xa q\\

X

(д)

 

EI

д

24


где матрица 7, формируется также для основной системы, по от-сил Я, =1, приложенных по направлениям выбранных

Запишем (г) и (д) совместно

0

Я / 2

I

W

B \ \ L X A 0\\

X

д

EI

 

Ц

( с )

1

 

 

 

 

 

 

Введем обозначения:

 

 

 

я /2

1

 

 

\\LxAq\\— А

Е/ ~

х

 

 

 

 

и

 

 

 

L'BA = C.

(ж)

Тогда формулу (е) можно записать как:

Л

с

X

 

Г !

 

±

 

 

Л !

 

 

 

 

 

 

(С—/. Е°)

X

= 0

(14)

А

 

 

 

 

 

где

 

0

01—вырожденная (не имеющая

 

 

0

Е

 

 

обратной) единичная матрица.

 

 

 

 

 

неоп­

Таким образом, задача'об устойчивости статически

ределимой конструкции свелась к решению

уравнения (14).

Заметим, что в нем не только матрица Е,

но и матрица С мо­

жет быть вырожденной, а поэтому

при его решении

нельзя

ограничиться отысканием .наибольшего

собственного

числа

1-шах

для матрицы С. Программа

решения

таких уравнений

приведена в Приложении.

 

 

 

 

 

 

 

Итак, для решения задачи об устойчивости статически не­

определимой конструкции необходимо сформировать

матри­

цы

L x, Lj, В, А, пользуясь выражением

(ж), получить мат­

рицу С и решить матричное уравнение (14).

 

 

 

Если в конструкции имеются

упругоподатливые опоры, го

тайную конструкцию можно рассматривать

как комбиниро­

ванную. В ней элементами, работающими только

на

осевую

нагрузку, будут эти опоры, причем

их податливость

можно.

Охарактеризовать коэффициентом

податливости

К. Тогда в

матрицы L,

В и А необходимо добавить

 

элементы с индек­

сом N [7]

 

 

 

 

7м„ х A m, i

L М, х L А ’, х

в =

Вм

=

L ' m, i

,

, Д

7д -, лг A x t«

h\i

 

Вх

 

25


Здесь у матриц, составленных для элементов, работающих на изгиб, введен индекс М.

Указанная методика хорошо реализуется на ЭВМ «Наири» для рам средней сложности. Для более сложных рам можно использовать алгоритм, указанный в работе [3].

П Р И М Е Р Ы

Пример 1. Найдем значение критической силы для рамы, показанной на рис. 5,а. Форма потери устойчивости заданной рамы при достижении силой Р критического значения, а так­

же примерный вид эпюры

в заданной

статически неопре­

делимой раме приведен на рис. 5, б,

в. Основная система после

приведения рамы к безразмерному виду

(за характерные ве­

личины приняты Р, I и EI)

показана на рис. 5, г. Форма поте­

ри устойчивости и предполагаемая

эпюра Мрп для основной

системы показана на рис. 5, д.

 

 

Охарактеризуем эту форму двумя перемещениями _М и Д2, которых достаточно для определения значений ординат эпю­ ры Мр. Для того, чтобы определить ординаты эпюры М разо­ бьем раму на 3 участка и выберем 5 сечений, которые полно­ стью определят эту эпюру (см. рис. 5, в). Сформируем исход­ ные матрицы.

Матрица податливости

 

В =

Я,

 

 

В и

 

 

 

В т

 

В ,

go

2

1

4

В и в п[

2

 

g 2

1

 

 

 

Так как эпюра М в сечениях 2 и 3 непрерывна, то

1

1 2

Для формирования матрицы 1.х построим эпюры М\,х и М2, х (рис. 5, е). Тогда

26