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

книги из ГПНТБ / Морозов, В. А. Регулярные методы решения некорректно поставленных задач-1

.pdf
Скачиваний:
23
Добавлен:
22.10.2023
Размер:
14.58 Mб
Скачать

а)

I г I J MJ M 3 \г \

'

___m

.2,7

 

2,7

 

44,15

81 SO/2 = 40 ,8 0

6)

В )

/ .

3 o,os

 

L->-------L _

 

 

 

*к ..I___

 

 

 

Ч

? -

 

\S'75\ | V f | .

32,70

j

1,35

1,35

 

1

Рис. VIII.6. Изменение поперечной жесткости по длине пролет­ ного строения:

а — участки постоянного сечення; б — фактическая эпюра моментов инерции; в — расчетная эпюра

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

№ участков . . . . .

1

2

3

4

5

Моменты инерции, м4

964

1058

1328

1428

4750

Изменение моментов инерции по длине пролетного строения изо­ бражено на рис. VIII.6, б. Как видим, моменты инерции участков 1 и 2, а также участков 3 и 4 незначительно отличаются друг от дру­ га. Поэтому указанные участки объединены попарно и для них при­ няты средние значения моментов инерции, что дает погрешность в пределах 5%. Кроме этого, для упрощения дальнейших расчетов граница между участками 2 и 3 переносится на 1,35 м в сторону средней опоры с тем, чтобы она совпала с точкой приложения со­ средоточенного груза. Таким образом, в расчете пролетное строе­ ние рассмотрено состоящим из трех участков постоянного сечения (на половине длины) с моментами инерции Iо, Ii, /ц. Принятая в расчете эпюра изменения моментов инерции по длине пролетного строения приведена на рис. VIII.6, в.

221

Расчетные моменты инерции в лг4 и долях Iо имеют следующие значения:

№ участков .................

 

О

I

II

Моменты

инерции,

лИ

1010

1378

4750

Моменты

инерции

в

/0

1.37 /о

4,70/0

долях /о ......................

 

Жесткости участков равны

соответственно

Е10, Е1%, Е1ц, где

Е = 3 500 000 Т/м2— модуль упругости 'бетона марки 400.

Жесткости сечений опор определены по их проектным размерам (см. рис. VIII.4). Столбы опор в поперечном направлении не связа­ ны друг с другом и работают независимо. Моменты инерции /п, /в верхнего и нижнего сечений двух столбов относительно центральной оси, параллельной оси моста, соответственно равны 5, 9 и 10,4 м4.

Эти величины приводим к бетону марки 400 (путем умножения на отношение модулей упругости бетонов марок 300 и 400) и выра­ жаем в долях момента инерции /о пролетного строения. Оконча­ тельно будем иметь:

 

3 150 000

5,9

0,0053/0;

в—

3 500 000

' 1010

 

 

3 150 000

1 М =

0 0 0 93/

3 500 000

1010

и

Жесткости верхнего и нижнего сечений опор в поперечном на­ правлении соответственно равны E Iв и Е1п.

У п р о щ е н и е д и н а м и ч е с ко й расчетной схем ы с учетом усл о ви й с и м м е тр и и . О пределение л и ш н е й неизвестной . При вычислении сил инерции поперечного направления ввиду симметрии системы доста-

а} -*|А л - А

Рис. VI11.7. Упрощенная динамическая расчетная схема:

а — фасад; б — план

222

точно

рассмотреть

только

 

план

 

 

 

формы

симметричных ко­

 

 

 

С

А^А

лебаний, для которых ка­

 

 

 

 

 

 

7 В 3 10 II

R

сательная

к

упругой

ли­

А>0 1 2

3 4 S

6

нии прогибов

в середине

 

t

 

■12

Р=1

 

-13

пролетного

строения

па­

 

 

 

РЧ ■W

раллельна оси моста.

Для

 

 

\R

 

77777

■определения этих форм и

 

 

 

 

 

 

 

 

 

соответствующих

собст­

 

 

 

 

 

венных

частот

 

(см.

 

 

 

 

 

§ VII.3) можно использо­

Рис.

VIII.8.

Основная система

 

вать упрощенную

расчет­

 

 

 

 

 

ную схему

(рис. VIII.7). Она включает в себя только половину мо­

ста. По оси симметрии установлена

жесткая

скользящая

опора.

Масса, укрепленная по оси симметрии, делится пополам, т. е. при­ нимаем т ц = 0,545т1 (см. табл. V III.1). Число степеней свободы рассматриваемой схемы /г = 14. В отношении симметричных форм собственных колебаний эта расчетная схема динамически эквива­ лентна исходной, приведенной на рис. VIII.2.

Принятая расчетная схема в статическом отношении является однажды статически неопределимой системой. В качестве лишней неизвестной удобно принять горизонтальную силу R взаимодейст­ вия между пролетным строением и промежуточной опорой. Основ­ ная система, полученная в результате устранения связи между указанными элементами, приведена на рис. VIII.8; там же показа­ ны положительные направления перемещений и сил.

Для определения единичных перемещений расчетной схемы предварительно строим линию влияния лишней неизвестной; в каче­ стве внешней нагрузки принимаем горизонтальную единичную силу Р = 1, прикладываемую последовательно в точках прикрепления со­ средоточенных грузов (рис. VIII.8). Ординаты линии влияния лиш­

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

 

R k= - J ^ ,

 

(VIII.1)

°RR

 

 

где 8д>р — обобщенное перемещение основной

системы

по направ­

лению силы R от единичной внешней силы

в точке

/г*; блн —

то же, перемещение от силы R = 1.

 

 

Указанные перемещения определяем с учетом только изгибных деформаций. Схемы для вычисления изгибающих моментов от еди­ ничных сил Р, R приведены на рис. VIII.9, VIII. 10.

При расположении единичной силы в пределах пролетного стро­ ения (£ = 1 — 11) его реакции равны А = Р 1, В = 0 и для изгибаю-

* Физически эти перемещения представляют собой расхождение точки 6 про­ летного строения и верхней точки опоры.

223

Рис. VIII.9. Схема построения единич­ ных эпюр моментов по пролетному строению

л»;

ничных эпюр моментов по опоре

 

 

р=/

 

 

о 1 г 3 4 5 В 7 8 9 ЮЦ1

 

 

с |d

ж

 

 

' R-1

 

 

hi.

 

 

 

 

3,65

 

 

 

11,75

 

 

 

19,85

 

 

 

27,35

 

 

 

35,05

 

 

 

45,15

 

 

 

52,25

 

 

 

60,35

 

 

 

55,45

 

 

 

75,55

 

 

 

«54,45

 

 

 

. 43,15

Рис.

V III.11.

Эпюры М р ,

M r в о с н о в ­

н о й

системе

при расположении еди­

ничной силы

в пределах

пролетного

 

 

строения

 

щих моментов в сечениях пролетного строения имеем формулы

(см. рис. VIII.9):

__

при x ^ a k

М р= х ;

при л :> а *

M p ~ a k,

где а* — расстояние силы Р = 1 от левой опоры. Координату х так­ же отсчитывают от левой опоры; M r совпадает с Мр при k = 6.

Единичные эпюры МР, построенные вышеуказанным способом при расположении силы Я = 1 в точках 1И , приведены нэ

224

рис. VIII. 11; опора при этом не деформируется и в

ее пределах

А /р= 0 . Эпюра M r приведена на том же рисунке.

 

Единичные перемещения определяем по формуле

Максвелла —

Мора. Интеграл записываем только в пределах длины пролетного

строения;

сообразно

с

участками

постоянной жесткости (см.

рис. VIII.6,

в) он разбит на пять интегралов и окончательно можно

записать:

 

 

 

 

 

 

 

•*(*)_ 1

X5 --

--

 

xd ------

Г

MpMpdx

, Г MpMpdx . г

MpMpdx

hRP~ —

 

 

 

 

 

 

Е J

 

Т0

h J

7 Г - + J

 

 

х 7

 

 

-VИ

л>

 

 

+ . f

 

 

 

 

 

 

М p M p d x ,

Г А'\рЬАpd x

(A = l,

2 ,.'.., 11),

 

 

 

 

 

 

где х5, хс, ха, Ху, хц — границы участков с постоянной жесткостью (см. рис. VIII.9). Моменты инерции Iо, /г, /ц были приведены выше.

Указанные интегралы вычисляем способом Верещагина (пере­

множением эпюр). Результаты

вычислений (значения

E I ^ p )

приведены в табл. VIII.2.

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а VII 1.2

Местоположе­

' - ' О

°RP

* *

Местоположе­

Ы 0о р р

 

точек

ние единичной

точек

ние единичной

 

 

 

F ,

S ( *>

 

 

силы Р=1

 

 

 

СИЛЫ P e l

 

 

 

 

 

 

1

На пролет-

9 229

—0,0412

8

На пролет-

120 728

—0,5390

 

НОМ

 

 

 

 

НОМ

 

 

2

строении

 

 

 

 

строении

 

 

То же

29 475

— 0,1316

9

То же

128 082

- 0 ,5 7 1 9

3

я

48 987

—0,2187

10

 

132 537

— 0,5918

4

69 089

— 0,2995

11

 

134 099

—0,5987

5

п

83 420

— 0,3725

12

На опоре

—92 818

0,4144

6

97 748

—0,4364

13

То же

—34 531

0,1542

7

 

110 480

—0,4933

14

— 4 360

0,0195

 

 

 

 

 

 

 

 

При расположении единичной силы в пределах высоты опоры (точки k = 12, 13,14) пролетное строение не деформируется и момен­ ты Мр в нем равны нулю. Эпюры моментов Мр, MR для опоры пока­ заны на рис. VIII.10. Единичные перемещения от внешней силы определяются формулой Максвелла — Мора:

1 f p= Г~MPMpdx

13) щ

J Б1Х

 

2 2 5

где Ix — переменный момент инерции сечений опоры; h — высота опоры.

Указанный интеграл можно вычислить путем замены непрерыв­ ного закона изменения момента инерции ступенчатым и интегриро­ вания по участкам. Мы используем для вычислений формулы при­ ложения I.

Предварительно определяем вспомогательные коэффициенты:

 

0,0053 =

0,5699;

0,0093

 

1 _ с = 0 ,4 3 0 1 ;

- 1 - 3 0

1,9183;

2 (1 —с)2

2с — 1

0,7557;

( 1 - с ) 2

1,1625:

(1 — с)

с-

4,0823;

( 1 - с ) 3

-----------= 3,0807. ( 1 - с ) 2

Кроме того, учитывая соотношение /н = 0,0093 /0, будем иметь

/г3А = 1 1 ^ ^ 3 3 4 600.

/н 0,0093

k

1

2

3

4

5

6

7

V

 

 

 

 

 

 

 

 

E I 05

Л О Г

3

 

 

 

 

 

 

 

1

0,6401

 

 

 

 

 

 

2

1,9219

5,9728

 

 

 

 

 

3

2,9930

9,4407

15,2831

 

 

 

 

4

3,8800

12,3172

20,2080

27,2041

 

 

 

5

4,6022

14,6587

24,2246

32,9938

40,5990

 

 

6

5,2086

16,6016

27,5600

37,8095

47,0103

55,0847

 

7

5,7284

18,3139

30,5007

42,0535

52,6601

62,2552

70,9798

8

6,1532

19,6926

32,8669

45,4701

57,2095

68,0292

78,0503

9

6,4583

20,6813

34,5637

47,9429

60,4713

72,1698

83,1220

10

6,6428

21,2803

35,5923

49,4049

62,4486

74,6784

86,1943

11

6,7075

21,4891

35,9514

49,9238

63,1398

75,5543

87,2681

12

3,8241

12,2648

20,2874

27,8015

34,5686

40,5060

45,7810

13

1,4228

4,5436

7,5484

10,3440

12,8618

15,0710

17,0335

14

0,1796

0,5737

0,9531

1,3061

1,6240

1,9027

2,1507

Тогда на основе формулы (4) приложения I для искомых еди­ ничных перемещений получим

E I J > $ = E I 0 Г м р Щ ^ х _ _ з з 4 б о д Г _ 1 > 9 1 8 3 _ j_ 0 7 5 5 7 £ » +

J

Е*х

 

I

Л

о

 

 

 

 

+ 1,1625 ( ^ ) 2~

(4,0823 +

3,0807 ^

In ^0,5699 + 0,4301 ^

 

=

12, 13,

14),

(VIII.2)

где Xk — ордината точек k от верха опоры (см. рис. УШЛО).

Результаты вычислений

по указанной

формуле

(значения

E I qR^p) 'приведены в табл. VI11.2.

При

вычислении

единичного

перемещения 6RR следует иметь в виду,

что

сила

R —1

вызывает

изгибающие моменты как в пролетном строении,

так

и

в опоре.

Соответствующие единичные

эпюры

приведены

на рис.

VIII. 10-

VIII.И.

 

 

 

 

 

 

 

Полное перемещение $RR состоит из двух частей, соответст­ вующих деформациям соответственно пролетного строения и опоры:

 

 

8я я = 8w?+<W

 

 

 

 

Очевидно,

b'RR

совпадает с перемещением

8

$

при £ = 6 (см.

табл. VIII.3):

E / 0b'RR= 97748.

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

VIII.3

 

 

 

ft

 

 

 

 

3

9

10

ii

12

 

13

*14

 

 

 

E I q5 . 1 0 - 3

 

 

 

 

 

 

 

vft

 

 

 

 

86,4439

 

 

 

 

 

 

92,5261

99,6251

 

 

 

 

 

96,2118

103,9815

108,9804

 

 

 

 

97,5012

105,5047

110,7576

112,7229

 

 

 

50,0276

53,0744

54,9200

55,5662

31,81222

 

 

18,6136

19,7472

20,4339

20,6743

12,2829

6,6073

 

2,3502

2,4934

2,58005

2,6104

1,7530

1,3487

0,40031

226

Рис. V III . 12. Линия влияния лишней неизвестной

Вторую часть перемещения Sw вычисляем по формуле (VIII.2) при Xh= 0; учитывая противоположные положительные направле­ ния сил Р и R по рис. V III.10, результат следует принять с обрат­

ным знаком. Будем иметь Д708я/? — 126 218. Тогда окончательно

E I $ r r = 97 748 + 126 218 = 223 966.

Далее непосредственно вычисляем ординаты линии влияния лишней неизвестной по формуле (VIII.1). Результаты приведены в табл. VIII.2. Линия влияния R графически изображена на рис. V III.12.

О пределение е д и н и ч н ы х перем ещ ений расче тно й схем ы . Для вычисления периодов и форм собственных колебаний предваритель­ но нужно определить единичные перемещения принятой расчетной схемы для точек k прикрепления сосредоточенных грузов. Их вычис­ ляют по формуле

Sv* = S v * —(--

 

( V III . 3)

где бvfc — единичные перемещения для

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

изобра­

женной на рис. VIII.8 (6„ь — перемещение точки v основной

систе­

мы поперек оси моста от силы Р = 1,

приложенной в точке k );

о — единичные перемещения для основной системы от силы Р = 1 ; Rh — значения лишней неизвестной при^диничной силе Р в точке /г.

Значения единичных перемещений

 

основной

системы при

расположении

внешней силы

на пролетном строении

(точки

k = l , 2, ...,

11)

вычисляют по формуле Максвелла — Мора.

Разде­

ляя, как и выше, полулролет на

пять участков с постоянной жест­

костью, для

8vft

получим выражение

 

 

 

 

 

 

 

г мM^Mpdxет рах

сп титрахM.,Mpdx

d

-Т7 - п

 

 

 

 

r>

M^Mpdx

 

 

Ъ.,ь= -

Т0

Ч

/,

 

4

 

 

 

 

 

\

 

 

 

 

г7 титрахM 4M pdx

XRг»

M.,Mpdx

k = l ,

2 , . . . ,

11).

(VIII.4)

 

 

 

 

 

(v,

Л

228

Знамения

единичных

перемещений

8vft

при

расположении

единичной

силы на опоре

(£=12,

13,

14) вычисляем по форму­ ле приложения I. Учиты­ вая соотношение между /и и /о и значение 1/1—с

= 2,325, по формуле (3)

получим

£/0ifvft= 3 3 4 600 X

X { — 1,9183 + 0,7557 X

X^ - 2 , 3 2 5 - ^ - X

кк

X

+ 1,1625 X

X ( + -

4,0823 +

+ 3’080 7 ( + + Т ) +

+ 2 , 3 2 5 ^ ' ^ 1

J

X

 

 

к

к

 

 

 

X In 0,5699 + 0,4301^-'

)■

 

 

 

 

 

 

 

 

 

(VII .5)

 

Значения

перемещений 6,R для точек v = l, 2, ...,

11 совпадают

со значениями о,*,

 

-вычисленными по формуле (VII 1.4) при k = 6.

Значения

blR

для точек v= 12, 13, 14 вычисляем

по формуле

(VII 1.5) при хн 0.

 

 

 

Значения Rk принимаем по табл. VIII.2.

Результаты вычислений по формулам (VIII.3) — (VIII.5) приве­ дены в табл. VIII.3 *.

В ней даются значения £706уй. Учитывая симметричность матри­ цы единичных перемещений (6й, = 8„ й), в таблице даны перемеще­

ния только для точек v.^'/г.

 

 

 

Полученные

единичные

перемещения

изображены

на

рис. VIII. 13. Данные на этом

рисунке графически представляют

* Эти вычисления весьма трудоемки. Все приведенные выше вычисления еди­ ничных перемещений могут быть выполнены на ЭЦВМ с помощью стандартных программ для расчета стержневых систем.

2 2 9

упругие линии прогибов пролетного строения и опоры системы, изо­ браженной на рис. VIII.7 от силы Р = 1, расположенной в точках v. Их можно рассматривать так же, как упругие линии исходной сис­ темы (рис. VIII.2) от двух сил Р = 1, расположенных в точках v симметрично относительно середины пролета.

Общее очертание упругих линий, приведенных на рис. VIII. 13, соответствует физической картине деформирования моста в попе­ речном направлении, а их плавность свидетельствует о правильно­ сти вычисления единичных перемещений.

Определение периодов и форм собственных поперечных колеба­

ний. Для определения периодов и форм собственных колебаний вычисляем собственные числа Хг и координаты собственных векто­ ров Xik матрицы вида

[ /72A 8 vft] 1 ,

где nik — сосредоточенные массы по табл.

V III.1

(как было указа­

но выше, массу т и нужно принять

с коэффициентом 0,5);

Svfc-— единичные перемещения по табл. VIII.3;

д=14 — порядок

матрицы (число степеней свободы системы).

 

Определение собственных чисел и координат собственных векто­ ров произведено с помощью ЭЦВМ «Раздан-2» по программе опре­ деления сейсмических нагрузок по матрицам единичных перемеще­ ний и масс, которая скорректирована Тбилисским зональным науч­ но-исследовательским институтом типового и экспериментального проектирования жилых и общественных зданий (ТбилЗНИИЭП) *.

Исходная информация для вычислений на машине была пред­ ставлена в виде матрицы-столбца масс в Т~1 сек2 н матрицы еди­

ничных перемещений

[о,Ц? в м Г-1.

 

Для получения

элементов

этой матрицы

числа, данные в

табл. VIII.3, умножали на величину

 

= ---------

1---------

= 0,283 •1СГ6 Г -1

м~\

Е10 3500 000-1010

 

 

Машиной были определены первые три собственные числа: А,|, Х% Хз и координаты соответствующих собственных векторов X ih, X2h, Xah (k — l, 2, ..., 14). Напомним, что они соответствуют сим­ метричным формам собственных колебаний исходной расчетной схемы (см. рис. VIII.2).

* Вышеуказанная программа требует, чтобы точки системы были занумерова­ ны по условию б11<622 < 6пп. Для этого в ходе вычислений была введена но­ вая временная нумерация точек, которую мы здесь не приводим. При других про­ граммах указанное условие не обязательно.

230

Собственные числа, опре­ деленные машиной, приве­ дены в табл. VIII.4. Там же даны круговые частоты и пе­ риоды колебаний первых трех собственных форм, вычисленные по форму­ лам:

 

 

Т а б л и ц а V I I I . 4

S

 

Круговые

Собственные

о.

Собственные

частоты

периоды Т-}

о

•&

числа Л.

о.1

сек

 

2

 

 

 

1

0,00718283

11,7992

0,5325

2

0,00015029

81,5713

0,07703

3

0,000019689

225,3675

0,02788

?/ = — L r

2, 3);

Vh

 

 

Т t = 2n У хТ (г = 1 , 2 , 3 ) .

Координаты собственных векторов, вычисленные машиной, при­ ведены в табл. VIII.5. Они определяют формы собственных колеба­

ний. Напомним,

что Xik есть ордината упругой линии г-го собствен­

ного

колебания в точке k. Все

формы

нормированы по условию

Kj,n =

l (с точностью до знака). Полученные формы колебаний изо­

бражены на рис. V III.14.

 

 

 

 

 

 

 

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

условий ортогональности

(см.

§

II.1). С этой целью сопоставляем

значения ЪпиХ^Х^ и X1mh\XihXjk\ для всех

пар

форм колебаний.

 

К

К

сумм производим

по условным

массам

Для простоты вычисление

 

 

 

 

 

 

 

 

Т а б л и ц а VIII.5

 

Ординаты форм собственных

 

 

Значения

Xj k

 

 

колебаний

 

 

 

 

 

 

 

 

Условные

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

массы

 

/= 1

/=1

 

 

i=i

1=2

/=3

 

 

i = 2

 

 

 

 

 

 

 

 

7-2

i= 3

j = з

 

 

 

 

 

 

 

1

0,06514

0,18280

0,33612

1,00

0,01191

0,02189

0,06144

2

0,20808

0,54866

0,91324

1,00

0,11416

0,19003

0,50106

3

0,34610

0,80000

1,02857

1,00

0,27688

0,35599

0,82286

4

0,47653

0,88075

0,63928

1,00

0,41970

0,30464

0,56305

5

0,59590

0,76786 — 0,01881

1,19

0,54451 —0,01334 —0,01718

6

0,70411

0,49274 —0,56681

2,04

0,70776

—0,81416 —0,56975

7

0,80339

0,13730 —0,81051

1,19

0,13127

—0,77488 —0,13243

8

0,88653

—0,26444 —0,54480

1,00

—0,23443 —0,48298 +0,14407

9

0,94824

—0,63516

0,10073

1,00

—0,60228

0,09552 —0,06398

10

0,98647

—0,90042

0,72652

1 ,00

—0,88824

0,71669 —0,65418

11

1,00000

— 1,00000

1,00000

0,545

—0,54500

0,54500 —0,54500

12

0,51838

0,38513 —0,75214

0,27

0,05390 —0,10527 —0,07821

13

0,19317

0,15411 —0,59209

0,33

0,00982 —0,03774 —0,03011

14

0,02449

0,02311 —0,16337

0,39

0,00022 —0,00156 —0,00147

 

 

 

YmkX ikX jk

0,00018 —0,00017

0,00017

 

 

 

2mk 1 X ikX jk |

4,54008

4,45969

4,18479

231

Соседние файлы в папке книги из ГПНТБ