Файл: Фоломеев, А. А. Снижение материалоемкости железобетонных конструкций-1.pdf

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

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

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

Добавлен: 29.10.2024

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

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

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

Тогда

оо

w (х, у, t) = 2 wi (■*. у) • Ti (0 =

оо Г

оо

 

 

оо

 

= ^ w . ( x , y)Zj(t).

( 1.20)

 

y-1

 

Очевидно, что система функций

 

оо

 

 

Wj (х, у) = 2

Ф*/’ W,(X, у) (у = 1, 2, 3,...)

( 1-21)

будет фундаментальной, a ощ w2, ... •— аппроксимирующими функ­ циями. В излагаемом методе решение систем уравнений движения отдельных составляющих элементов коробки при соответствующих граничных условиях заменяется более простым с точки зрения ма­ шинного счета решением — определением собственных чисел и соб­ ственных векторов матрицы Л-1Б.

Суть приближенного метода, ранее примененного В. В. Болоти­ ным [13] и В. К. Кабуловым [36] в исследованиях динамической устойчивости и колебаний плоских рамных систем, заключается з том, что в сумме [I. 20] сохраняется конечное число членов.

§ 3. ВЫБОР АППРОКСИМИРУЮЩИХ ФУНКЦИЙ

Объем вычислений, как и точность приближенного решения за­ дачи с помощью указанного выше метода, зависит от выбора ап­ проксимирующих функций Ш|, w% w3, ..., w п. Строгих и четких пра­ вил, которыми следовало бы руководствоваться при выборе аппроксимирующих функций, нет. Однако в качестве общих рекомендаций могут быть приняты следующие замечания.

Аппроксимирующие функции надо подбирать так, чтобы они составляли полную систему функций/ удовлетворяли геометриче­ ским граничным условиям (I. 10), а также соответствовали физи­ ческому смыслу решаемой задачи. При построении этих функций нужно исходить из возможных форм колебаний системы.

Удачным выбором является такой, при котором облегчается вычисление интегралов (I. 13) и при котором матрицы А и В име­ ют большие диагональные члены в сравнении с членами, не лежа­ щими на диагонали. Иначе говоря, мы стремимся сделать аппрок­

симирующие wi (x, у) и фундаментальные w. (х, у) функции на

13


столько близкими, насколько это возможно, стремясь при этом не увеличивать объем вычислительной работы.

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

§ 4. ОПРЕДЕЛЕНИЕ

ЧАСТОТ И ФОРМ

КОЛЕБАНИЙ

КОРОБКИ С УЧЕТОМ

ЖЕСТКОСТИ ИЗГИБАЕМОЙ ПАНЕЛИ

 

 

 

 

 

 

 

 

 

Задачу будем решать при следующих исходных данных (рис. 2):

изгибаемые панели

(1,

3) из керамзитобетона

с

удельным весом

 

 

 

 

1,2 т/м»,

£ = 75 000

кг/см2,

о=0,25,

 

 

 

 

высотой

3 м,

шириной —4,

толщи*

 

 

 

 

ной —0,28 м;

панели,

работающие

 

 

 

 

на сдвиг (2, 4), из тяжелого бетона

 

 

 

 

с

удельным

 

весом

2,5

т/м3,

£ =

 

 

 

 

= 200 000

кг/см2, G'=0,4£, шириной

 

 

 

 

5 м, высотой—3, толщиной—0,12 м;

 

 

 

 

перекрытия

из тяжелого

бетона с

 

 

 

 

удельным весом 2,5 т/м3, длиной 4 м,

 

 

 

 

шириной-— 5,

 

толщиной — 0,1

м.

Рис. 2. Коробка крупнопанельно­

Временная нагрузка

на

перекры­

тия — 0,2 т/м2.

 

 

 

 

 

го здания.

 

 

 

В

качестве

аппроксимирующих

 

 

 

 

функций выбираем

 

 

 

 

 

 

то,

 

y>i +

v,

i =

1, 5,

 

 

 

(1.22)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

=

.

(2/ —1) тс

у, =

 

.

(2i —1) л

 

 

 

 

sin -------х\

Sin -------- г™——

 

 

 

 

 

 

 

 

 

 

 

 

 

•у;

 

 

 

 

.

,

,

0,862

,

3,426

 

 

 

 

°i = sm \

У\ \ = —ft— >h = — b~

 

 

 

 

, _

6,437

, ,

9,529

_i

_

12,645

 

 

 

 

 

Кл ~

 

ь ’ '4

b ’ Л5 —

b '

 

 

 

 

Функции x t, yi

подобраны исходя

из

соображения,

что из­

гибаемые панели будут совершать кососимметрические колебания. Функции vt — форма сдвига панели в плоскости действия сил

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

нагрузки. Для рассматриваемого нами случая ^ - = 1, а значени я

\ заимствованы из работы [74а] {ml — распределенная масса панели, работающая на сдвиг; М — присоединенная масса перекрытия с учетом временной нагрузки).

14


Аппроксимирующую функцию выбрали

так, чтобы она удов­

летворяла

геометрическим

условиям сопряжения. На самом деле

при у = 0,

w . = 0 — перемещение

в основании отсутствует;

при

у = Ь, щ. = v tj y^ b =

const — перекрытие не

деформируется;

при

х — 0 ( х = а),

wt =

панели

(/,

3)

принимают

на кромках

х = 0, х — а форму

сдвига панелей

(2,

4).

 

 

 

Можно

считать,

что

выбранная

функция wt удовлетворяет

требованиям, предъявляемым вариационным методом.

по­

Подставляя

(1.22)

в (1.13) и учитывая

исходные

данные,

лучаем следующие значения для элементов матрицы А, 5-1СГ6 ,

А-1 7М 0“6 :

 

 

А

 

 

1,28310

0,43391

-0,10000

0,06762

—0,04533

0,43391

0,05632

0,3684

0,06332

0,01550

-0,10000

0,03684

1,07623

0,01232

-0,00874

0,06762

0,06332

0,01232

1,07493

0,00700

-0,04533

0,01550

-0,00874

0,00700

1,07481

 

 

В -10-6

 

 

0,282494776

0,342306036

0,650811640

2,14694145

1,46216990

0,003646636

4,427522292

— 0,271800304

- 0,408197612

0,558166664

0,005814304

0,065462784

24,7308406

0,426121884

0,332712528

0,010113192

0,394649384

— 0,131621452

80,8334335

0,138635812

0,005668620

0,002170856

0,133682332

0,19612064

205,162398

 

 

А - 1 -В -10-6

 

 

 

 

0,259413074

-1,37311215

3,16401577

-

0,802190639

10,8758267

-0,103337240

4,75493863

— 2,350872232

— 4,56529326

7,81768995

0,033416972

—0,232406536

23,3632599

— 0,387952522

-

3,15376503

-0,020138065

0,176979222

— 0,4-53486584

-

75,5264522

 

1,63570024

0,018099509

-0,127459494

0,483704443

0,281006214

 

191,491565

Собственные числа (квадраты

частот)

данной матрицы равны

\ =

ш] =

184303,074; \

= «в* = 4760914, 16;

Х3 =

ш23 =

23384750,2; \

= со2 =

75514503,7;

1 5=

о)2 =

19151157.

 

 

Матрица собственных векторов следующая:

0,99975

—0,27969

0,13991

-0,00941

0,05736

0,02205

0,96001

—0,12567

—0,06440

-0,04166

- 0

00121

0,01234

0,98211

-0,00700

0,01880

0,00021

-0,00239

0,00893

0,99785

—0,01421

—0

00008

0,00065

—0,00292

0,00238

0,99720

15


Для вычисления элементов матриц А~1, Л_1-Б-10_6, собственных чисел и соответствующих им собственных векторов в настоящей и последующих задачах мы использовали стандартные программы.

*77777777

Рис. 3. Формы свободных колебаний.

а для вычисления элементов матриц А и В — составили стандарт­ ные программы, блок-схема которых приведена в приложении 1.

Периоды свободных колебаний определяются по формуле

Т. = -Щ = .

VsЛ ■

Период основного тона настоящей задачи равен 0,0146 сек. Формы свободных колебаний, вычисляются по формуле

5

Wj = 2

*»,(■*, у),

i=i

 

а для случая х = 0 (х = а) показаны на рис. 3.

§ 5. ОПРЕДЕЛЕНИЕ ЧАСТОТ И ФОРМ КОЛЕБАНИЙ КОРОБКИ БЕЗ УЧЕТА ЖЕСТКОСТИ ИЗГИБАЕМОЙ ПАНЕЛИ

В этом случае массу изгибаемых панелей перераспределим на перекрытия и панели, работающие на сдвиг. Жесткостью изгибае­ мой панели будем пренебрегать. Тогда в выражениях (I. 13) отпа­ дут члены, содержащие поверхностный интеграл и жесткость D. Аппроксимирующими функциями будут

 

 

=

sin >-гу; г =

1,5,

 

 

(1.23)

где

0,784 л

3,370

л

6,405

>

9,507

_

12,628

~~Ъ

b I ' 3 ~

Ъ ’ 4 —

Ь , 5 ~

Ь

 

16

Рис. 4. Формы свободных колебаний без учета жесткости изгибаемых панелей.

Для рассматриваемого случая -^- = 0,78125 и \ вычислены

на основании работы [74 а]. Подставляя (1.23) в (1.13) и учитывая указанные выше замечания, получаем:

 

 

А

 

 

0,80750

0,06235

-0,03336

0,02262

-0,01702

0,06235

0,63274

0,01048

—0,00706

0,00527

-0,03336

0,01048

0,61862

0,00389

—0,00296

0,02262

—0,00706

0,00389

0,61547

0,00202

—0,01702

0,00527

—0,00296

0,00202

6,61438

1

 

 

в - к г 6

 

 

 

 

0,19319400

—0,00012672

0,00012672

— 0,00009216

 

0,00059904

0,0001152

2,319390072

-0,00031104

— 0,00011520

 

0,00021888

0,00004608

-0,00021888

8,02485504

 

0,00006912

— 0.00160128

0,00003456

—0,00001152

0,00020736

 

17,5027737

 

0,0023616

0,00008064

-0,00003456

-0,00043776

 

0,00119808

 

30,7692288

 

 

Л - 1 • В • 10“ 6

 

 

 

 

0,242101187

-0,291211396

0,563326081

- 0,841166514

1,11700233

-0,024426047

3,69657365

■0,273306637

0,405264826

- 0,534602980

0,013635630

-0,079212296

13,0082468

 

— 0,231746338

0,307837118

-0.С0Р3429С8

0,053720271

■0,105982131

28,4754822

 

— 0,209787742

0,006381541

-0,040390156

0,080257935

 

— 0,119568936

50,1194587

2-207


Собственные числа ^(квадраты частот) данной матрицы равны

V =

(0? =239113,571;

\ =

а>2 =3694911,47; \

= ш* =130088^7,1;

Х4 =

= 28477054,0; \

= ш52 = 50121915,4.

 

Матрица собственных векторов следующая:

 

 

0,99997

-0,08186

0,04407

-0,03001

0,02271

 

0,00693

0,99640

-0,02899

0,01641

-0,01166

 

-0,00102

0,00851

0,99858

—0,01497

0,00838

 

0,00031

-0,00215

0,00694

0,99928

—0,00976

 

—0,00012

0,00086

—0,00217

0.С0561

0,99959

Период основного тона 0,0128 сек. Формы свободных коле­ баний, вычисленные по формуле

5

vj = Z tfvi (У)'

приведены на рис. 4.

§ 6. ОПРЕДЕЛЕНИЕ ЧАСТОТ И ФОРМ КОЛЕБАНИЙ С УЧЕТОМ СДВИГА ПЕРЕКРЫТИЙ В СВОЕЙ ПЛОСКОСТИ

Как уже было отмечено, проведенные в Советском Союзе и за рубежом натурные исследования колебаний зданий показывают, что деформации вертикальных и горизонтальных элементов (пере­ крытий) по длине здания неодинаковы. В ряде случаев существен­ ны деформации в своей плоскости не только сборных, но и моно* литных перекрытий.

Сучетом деформации сдвига перекрытий элементы матрицы А

иВ вычисляются по формулам

Аи =

J j 4 , w k wfixdy + j'mc v k v.dy + j mn uk иf i x ,

Bu =

^ Dl^ wiw bdxdy -

~d w vk dy -

(1.24)

В качестве аппроксимирующих берем функции

18