Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Строительная механика и металлоконструкции строительных и дорожных машин

..pdf
Скачиваний:
8
Добавлен:
19.11.2023
Размер:
27.58 Mб
Скачать

 

У'

q=2>кн/м

5

Р=5,7кН

с ' l U U I U I

см

^

J

 

 

к3м

VS,2%------------------

J

 

 

 

а)

 

 

М1 [ кНм]

Г | М °Р[ *НМ] Г ж-ч

X

 

 

 

f/C ^X

 

 

 

*)

 

 

4£/

 

 

 

 

 

 

•)

 

 

 

 

 

 

 

 

1,5EJ

 

 

 

 

 

? J

^ 3EJ

J

r }

 

 

 

 

 

21

 

 

 

 

0,75 EJ

 

 

 

 

 

 

3EJ

^

*Г22

 

 

 

 

1

 

к

Р=5,7

— 5 , 5

 

 

Т

?

 

- 2 , 4

?

 

 

 

 

 

 

ж )

Рис. 1.46. Схемы для расчета плоской рамы методом перемещений

реакция направлена в сторону, противоположную соответствующему перемещению.

5. Решаем систему уравнений (1.98) и получаем Z x = 1 ,2/(£ /); Z2 = 2/(£У) м.

6. Строим в соответствии с выражением ( 1 .101) суммарную эпюру моментов Мр (кН • м) в заданной системе (рис. 1.46, з) .

7 . Производим универсальный контроль решения —проверяем рав­ новесие узлов и отсеченных частей (рис. 1.46, и). Эпюры Qp и Np в заданной сйстеме строятся по правилам, изложенным ранее (см. ме­ тод сил).

Пространственные рамы рассчитывают в той же последовательности. Особенность состоит в том, что в свободные жесткие узлы вводят прост­ ранственные плавающие заделки (изображенные в виде параллелепипе­ дов), которым задают угловые перемещения в трех плоскостях. От поворота узла одни стержни (лежащие в плоскости поворота) изгибают­ ся, а другие (лежащие в перпендикулярной плоскости) закручиваются (если имеются связи, препятствующие повороту сечения). Так как

угол закручивания в этом

случае определяется по формуле

=

= Мкр//(С /кр),то п ри^= 1

Мкр = GJKp/l.

 

Коэффициенты r.k и свободные члены Rip определяют на основе эпюр изгибающих и крутящих моментов.

На рис. 1.47, а изображена плоскопространственная система, имею­ щая упругоподатливую линейную связь в точке С, а на рис. 1.47, б — соответствующая основная система. Число отличных от нуля основных неизвестных равно трем (угол поворота узла в плоскости рамы заведо­ мо равен нулю).

При построении эпюр следует учесть, что от линейного перемещения Zi = 1 , а также от силы Р в элементах рамы возникнут только изгибаю­ щие моменты; от угла поворота Z 2 = 1 в элементах ЕС и CD возникнут изгибающие моменты, а в элементах АВ и CD — крутящие моменты; от поворота Z 3 = 1 в элементах АВ и CD возникнут изгибающие момен­ ты, а в элементе CD — крутящие моменты (кручения в элементе СЕ не будет, так как линейная опорная связь Е не препятствует поворотам сечения).

Вырезав узел С в различных состояниях системы из уравнений равновесия 2 Z = О, ЪМХ = 0 и ХМу = 0 определим коэффициенты rik и Rjp\ далее при составлении системы канонических уравнений следует учесть, что согласно выражению (1.103) Гц =r®x + 1/ф.

Особенности расчета симметричных рам. Основная система любой симметричной системы должна быть симметричной. При этом так же, как и в методе сил, нагрузку следует разложить на симметричную и кососимметричную составляющие, а неизвестные — сгруппировать. Особенно удобен метод перемещений при расчете симметричных систем на действие симметричной нагрузки. В этом случае в силу симметрии будут заведомо равны нулю основные неизвестные, определяющие линейные перемещения частей рамы (ригелей), расположенных поперек оси симметрии, а также углы поворота узлов рамы, лежащих на ее оси

Рис. 1.47. Схемы для расчета плоскопространственной рамы методом перемещений

я

Рис. 1.48. Схемы для расчета симметричной рамы методом перемещений

симметрии. Естественно, будут также равны нулю все кососимметрич­ ные неизвестные.

Изображенная на рис.

1.48, а рама 7 раз кинематически неопреде­

лима. Здесь пу = 5 и пл ^

2 (см. рис. 1.48, б) . Однако в силу сказанного

выше Z 3 = Z 4 = Z 5 = Z 6 = Z 7 = 0 (рис. 1.48, в) и, таким образом,

при расчете рамы надо решать только два уравнения с двумя неизвест­ ными. На рис. 1.48, г - е построены единичные и грузовая эпюры мо­ ментов. Следует отметить особенность характера эпюры моментов на верхнем ригеле в первом единичном состоянии — она имеет прямо­ угольное очертание, так как построена от одновременного воздействия

парных углов поворота Z\ = 1.

Определив Т ц /2 = 6EJ/l\ г12/2 = 2EJ/J\ г2212 = \ \EJ/l\ R \p!2 =

=

- q l2/3 и R 2p = 0 из решения системы уравнений, получим Z x =

=

5,5ql3/(93EJ) и Z 2 = -q l3/(93EJ). На рис. 1.48, е - з показаны

эпюры внутренних сил в заданной системе.

1.6. МАТРИЧНЫЕ МЕТОДЫ РАСЧЕТА СТЕРЖНЕВЫХ СИСТЕМ

1.6.1. ОСНОВЫ МАТРИЧНЫХ МЕТОДОВ РАСЧЕТА

Общие понятая. В связи с широким использованием ЭВМ в инже­ нерной практике многие задачи расчета конструкций решаются в мат­ ричной форме. Благодаря использованию матричного языка изложение методов расчета получило более компактную форму, существенно упростилось программирование на ЭВМ и тд. Следует отметить, что при расчетах в матричной форме континуальная задача заменяется дискрет­ ной — значения искомой функции (усилий или перемещений) опреде­ ляются не для всех точек базиса, а лишь для узловых точек. На интерва­ ле между узловыми точками значение функции интерполируется.

Предполагая, что читатель знаком с матрицами и действиями над ними, напомним смысл матричных операций на простейших примерах.

Запишем принцип независимости действия сил в обычной форме:

Si = 5цР\ + si2P2 +

+ sjkPk +

+ sinPn9

где Sf - усилие в i-й связи системы;

- сила,приложенная в к-й точке; каждый

из коэффициентов s.^ численно равен усилию 5/ от соответствующей безразмерной силы = 1.

В матричной форме эта формула записывается так:

 

S12

s/n] х

Pi

Sj = [ i l l

Рг

Рп

Если необходимо вычислить не одно, а т усилий 5; от действия той же группы сил, то аналогично выражению (1.104) записывают:

51

s11

*12

* 1 Я

Px

52

S21

*22

Sin

Рг

 

 

 

 

(1.104)

 

*m i

*m 2

Smn

Pn

или сокращенно:

 

 

 

-►

-*

 

 

(1.105)

S =

LSP ,

 

 

где S - матрица-столбец (вектор) внутренних сил; L - матрица коэффициентов линейного преобразования; Р - матрица-столбец (вектор) нагрузки.

В общем случае т Ф п. Например, в ферме (рис. 1.1, б), состоящей

из

И

стержней и нагруженной силами в пяти узлах, т = 11 и п = 5.

 

Аналогичный вид

имеет формула

перемещений А к от

действия

нагрузки:

 

 

 

 

 

А , '

’ « и

б 1 2

5

Pi

 

 

Д 2

« 2 1

622

&2П

Рг

 

 

 

 

=

 

X

1______

(1.106)

 

А ш

5 m 1

f>m 2

I

 

 

&тп

 

 

или сокращенно:

 

 

 

 

 

Д

=

D Р ,

 

 

 

(1.107)

где

D

=

Ьд - матрица

коэффициентов линейного преобразования,

называемая

иногда матрицей влияния.

 

 

 

 

 

Чаще всего в расчетах необходимо вычислить перемещения Д^ всех

узловых

точек, где приложены обобщенные силы Рк . В этом случае

т = п и матрица коэффициентов линейного преобразования будет квад­ ратной размером п х п.

В строительной механике решают и обратные задачи: по известным перемещениям Д^ или внутренним силам Sk определяют внешние силы

Рк. В матричной форме это записывается так:

 

Pi

'bi 1

b it

Ь\т

Ai "

Рг

Ьг 1

b22

^ 2 т

д2

 

 

 

X

 

 

 

 

Рп

Ьщ

Ьп2

Ьпт

Ат

-

 

 

 

- -

или сокращенно:

 

 

 

Р =

В А ,

 

 

(1.109)

где матрица линейного преобразования В в общем случае так же, как и матрицы L или D, - прямоугольная (m X л ) .

Если т = п, то матрицы В и D будут взаимно обратными

(DB =

= BD = Е) и

 

В = D ”1

(1.110)

Операция вычисления матрицы В по имеющейся матрице D (или на­ оборот) называется обращением матрицы.

Выражение действительной работы внешних сил в матричной форме:

1

Т -►

( L i n )

А =

Р А,

2

 

Т

 

 

где матрица-столбец Р записана в виде матрицы-строки: Р

(

Л

)•

Рис. 1.49. Схемы, поясняющие податливость и жесткость элемента

С учетом выражения (1.107) эту же формулу можно записать так:

“►т

А = 1 Р D P.

( 1.112)

Операция замены столбцов строками (и наоборот) называется транспонированием матрицы.

Матрицы податливости и жесткости системы. Квадратная матрица коэффициентов линейного преобразования, входящая в формулу (1.106), называется матрицей податливости:

'

« 11

5i 2

&1П

 

D =

52 1

522

П

(1.113)

 

1

&П2

Snn

 

Элемент бу матрицы D представляет собой перемещение точки / по направлению силы P j = 1 , вызванное действием единичной силы Pj. Так как 5 .. = 5.-, матрица податливости является симметричной. Отметим, что матрица (1.113) является невырожденной — ее опреде­ литель не равен нулю: Det D Ф 0.

Понятию податливость противостоит понятие жесткость. Например, для пружины (рис. 1.49, а) податливость 6 равна перемещению, выз­ ванному силой Р= 1, а жесткость г равна силе Р (рис. 1 .49, б) , вызываю­

щей единичное удлинение А =

1 пружины. Эти коэффициенты входят

во взаимно обратные зависимости:

 

А = ЬР

и

Р =

гД,

 

(1.114)

где 6 = г ' 1 иг = 6 _1.

 

 

(1.115)

Наряду с понятием матрицы податливости системы аналогично

зависимостям

(1.115) вводится

понятие матрицы жесткости или мат­

рицы реакций:

 

 

 

 

 

г 11

т1 2

Г\П

 

R =

?2 1

Г22

Г2П

( 1 6)

 

Гп\

Гп2

Гпп

 

где Гу - реакция во введенной /-й связи от перемещения А. = 1

(см. метод пере­

мещений) .

 

С помощью матрицы R вектор сил Р можно выразить через вектор

перемещений А аналогично зависимости (1.114) :

 

Р = R А .

(1.117)

Отметим, что матрицы жесткости и податливости взаимно обратны:

R = D 1 и D = R -1.

(1.118)

Существование матрицы R следует из условия Det D Ф 0. Из симмет­ рии матрицы D следует, что при ее обращении и матрица R будет сим­ метричной, т.е. г- = г..\ это соответствует выражению (1.46).

1.6.2.РАСЧЕТЫ В МАТРИЧНОЙ ФОРМЕ

Определение перемещений. Запишем в матричной форме интеграл Мора для /-го участка при определении перемещений, зависящих только от деформаций изгиба в одной плоскости:

 

Mk Mpd*

^

(1.119)

% " I

--------M/ * A/ V

 

m e ы/к транспонированная матрица-столбец единичных моментов на/-м участке:

»/к |м% Кк м/к 1;

- диагональная матрица податливости на /-м участке:

4

0

0

4а?

 

0

о

6EJ0

 

1

0

О

 

( 1.120)

здесь а.. = J0/J>у - относительная податливость ( / 0 - любой фиксированный мо­ мент инерции, принятый за основной, J.. - момент инерции.в /-м сечении элемента l.)\ Mjp - грузовая матрица-столбец на/-м участке:

JP

м/р -

В этих формулах верхние индексы л, с и п указывают, что величины соответствуют левому, среднему и правому сечениям элемента /..

Легко убедиться в том, что перемножение матриц в формуле (1.119) полностью совпадает с вычислением интеграла Мора с использованием формулы Симпсона (1.59).

Суммируя результаты вычислений по всем п участкам системы, по­

лучаем:

 

г < а, м

(1.121)

**Р

 

/= 1

Это выражение можно записать в блочной форме:

 

- Т

м

Т

. -Т

 

 

 

Акр =

£м и

2к

 

 

х

 

A i

0

 

 

о"

 

’ м

 

0

а 2

 

 

0

X

М2р

( 1.122)

0

0

 

 

Ап

 

Мпр

 

Для

вектора перемещений Ар, с компонентами А 1р,

Д2р,

Атру вместо одной строки

надо написать т строк:

 

Ар

= МТАМр.

 

 

(1.123)

Если на /-м участке EJ. = const, то диагональная матрица податли­

вости ( 1 .120) имеет вид:

 

 

 

 

1

0

0

 

 

'/

4

0

(1.124)

А / =

0

6 E J .

0

1

 

 

0

 

Если на /-м участке обе эпюры моментов М^ и Мр линейны, то их средние ординаты можно выразить через крайние. В этом случае фор­ мулу (1.119) можно записать так:

М М dz

ij.

E J

6 E J 0

(« ;+ ф

<7С

 

 

+

 

1

[M?k M *]

х

м *

 

 

IP

 

(1.125)

X

1______

Я

 

При EJ. = const матрица податливости в выражении (1.125) имеет

вид:

 

'/

(1.126)

А/ =

6 E J,J L

 

 

 

Она может быть записана еще проще, если на одном конце участка имеется шарнир и нет сосредоточенного момента:

 

'/

(1.127)

А/ = -

[ 2 ].

6 E J ;

 

Естественно, что при этом моменты М и Мр в узловой точке для шарнира не записываются, так как они заведомо равны нулю.

Отметим, что уменьшение порядка матрицы податливости играет существенную роль при вычислениях на ЭВМ.

Если при расчете системы эпюры внутренних сил и податливостей непрерывны, то нет смысла дважды вписывать каждую их ординату на границах участков в соответствующих матрицах. Можно каждую орди­ нату в матрицах моментов вводить один раз, но при этом в блочной матрице податливостей наложить углами соседние блоки, т.е. сложить значения ’’пограничных” элементов. Например, при EJ = const и четырех участках одинаковой длины I. = d блочная матрица податливости соглас­ но выражению (1.126) будет иметь вид:

1

А =

I 1

~ i~ !

6 EJ

1___

 

Рассмотрим пример. Определим перемещение точки 4 балки переменного сечения (J0/Jz = 1 + z / 0 , нагруженной в узловых точках силами Р (рис. 1.50,а).

Вычислив узловые и промежуточные значения моментов Мр и М а также относи­ тельных податливостей о (рис. 1.50, б —г ) , запишем их в блочной форме и прове­ дем в соответствии с выражениями (1.119) и (1.120) вычисления. Учитывая при этом сложение пограничных элементов, получим:

P d 3

 

3

2,5

2

1,5

1

0,5

0

0

0 ] X

 

6 EJ.

 

 

 

 

 

 

 

-1

г

 

 

 

 

0

0

0

 

0

 

1

0

0

0

0

10

 

0

1,125

0

0

0

0

0

0

0

8

 

0

0

2,5

0

0

0

0

0

0

6

 

0

0

0

1,375

0

0

0

0

0

4,5

 

0

0

0

0

3

0

0

0

0

X 3

 

0

0

0

0

0

1,625

0

0

0

2

 

0

0

0

0

0

0

3,5

0

0

1

 

0

0

0

0

0

0

0

1,875

0

0,5

 

0

0

0

0

0

0

0

0

2

0

 

PI3

 

 

 

 

 

 

 

 

 

=

17,07

 

 

 

 

 

 

 

 

 

 

~еГ0

 

 

 

 

 

 

 

 

 

где J0

- момент инерции сечения балки в точке 1

(при 2 = 0).

 

 

 

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

I Расчет статически неопределимых рам методом сил. Система кано­ нических уравнений метода сил в матричной форме описывается выра­ жением (1.87).

Рис. 1.50. Определение перемещения в мат­ ричной форме

При рассмотрении ряда комбина-

ций нагружений векторы X и D р заме­ няют на матрицы X и Dp :

Х ц

X 12

X \k

X = х 21

Х 22

Хгк

Хщ

Хп2

Хпк

А \рк

&2рк

&прк

где к - число рассматриваемых вариантов загружения.

При расчете системы на температурные или дислокационные воз­ действия надо в матрице Dp заменить индексы р на t или d.

Изложим порядок расчета статически неопределимых систем в мат­ ричной форме~на'примере плоской рамы, изображенной на рис. 1.51, а, 1. Выбираем основную систему (рис. 1.51, б) и составляем исходные матрицы М, М° и А, разбивая длины стержней по участкам, как это по­

казано на рис. 1.51, а, В соответствии с рис. 1.51 ,б - д получаем

 

-1

0 "

’о "

 

 

-1

0

0

1

М

-1

0,5/ ;

м ; - 0

 

-1

0,5/

М

 

-1

/

м

 

 

F

0

0

0

0

 

0

2

1

0

0

А

1 / 2

1

2

0

0

0

 

6EJ

0

0

2

1

 

0

 

0

0

0

1

2

Вматрице А первый блок, состоящий из одного элемента, относится

кучастку 0 - 1 (длиной /, см. рис. 1.51, а), второй —к участку 2 - 3 и третий —к участку 4 - 5 ,