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

книги / Метод конечных элементов в динамике сооружений

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
6.89 Mб
Скачать

iio ним описать огибающую, которая й позволит судить о величине напряжений и о положении опасного сечения. Для пояснения изложенного рассмотрим простейшую за­ дачу одного измерения.

Стержень постоянного поперечного сечения (рис. 14) заделан одним концом; свободный конец стержня на­ гружен продольным сосредоточенным мгновенным им­ пульсом S. В процессе дальнейших рассуждений будем

считать

этот

импульс

постоян­

ным. Пусть после этого длина

стержня

1\

оказалась разрушен­

ной, а оставшийся участок стер­

жня /3 стал

совершать

упругие

колебания.

Рассматривая

стер­

жень как

систему

с одной

сте­

пенью свободы, принимая в рас­

чет только

первую

частоту

ос­

новного

тона

и

соответствую­

щую ей

собственную функцию,

получим формулу для

вычисления

перемещений

 

точек

стержня

 

 

 

у = asin

 

sin ©t t.

 

 

 

 

 

 

2/3

 

Для вычисления частоты получим выражение

Из этой формулы видно, что частота изменяется об­ ратно пропорционально длине стержня. Если через ©0 обозначить частоту, соответствующую полной длине стержня I, то частота стержня уменьшенной длины 1$

Эта формула показывает, что с увеличением длины разрушенной части стержня частота свободных колеба­ ний оставшейся части возрастает, и если бы импульс, приложенный к упругой части, не изменился, то с увели­ чением воронки возросли бы напряжения. Но в действи­ тельности с увеличением воронки все большая часть энергии импульса затрачивается на местные разрушения и на долю общих деформаций остается только часть об­ щего эффекта. Чтобы учесть это, подсчитаем кинетиче-

Скую энергию системы с учетом той энергии, которая со­ общается разрушенной части стержня:

2« = j ( - f -)2М* + Й = а2в? ц (-^2- + / , ) ,

О

Приравняем это выражение удвоенной работе внеш­ него импульса; тогда получим выражение для динамиче­ ского коэффициента

(h!2 + k)

Обозначим через а0 динамический коэффициент, со­ ответствующий стержню, не имеющему разрушенной части:

2S

со0ц/

Путем сравнения полученных формул приходим следующей зависимости между а и а0

а = а0 = ДоР;

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

Несмотря на уменьшение длины стержня, динамиче­ ский коэффициент становится меньше. В этом сущест­ венная разница в выводах, которые получаются при учете и без учета местных разрушений. Значение коэф­

фициента р в

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

отношения IJI указано в

табл. 1.

 

 

 

 

 

Таблица 1

 

 

 

 

 

 

и п

Р

UH

Р

и п

Р

и п

Р

0,1

0,82

0,6

0,25

0,3

0,54

0,8

0,11

0,2

0,67

0,7

0,17

0,4

0,43

0,9

0,05

0,5

0,34

1,0

0

Данными табл. 1 воспользуемся для построения ха­ рактеристических линий напряжений. Величину нор­ мальных напряжений в поперечных сечениях получим

делением нормальной силы на площадь Поперечного сечения по формуле

оN_= E J l . = Еа---- cos —

F дх 2/3 2/3

Наибольшее

напряжение в крайнем волокне будет

в сечении при х= 0,

т. е. в заделке ама1{с=Еа(п1213).

Если разрушенная

часть стержня равна нулю, то

° £ а к с = £ а ( Я / 2 / ) -

Т о Г Д а a M a . < c / a J a K C = [ l / ( l + / i / / ) ] .

Эта формула показывает, что при учете энергии, ко­ торая расходуется на местные разрушения, нормальные напряжения в опасном сечении уменьшаются, несмотря на то что упругие колебания совершает стержень уко­ роченной длины. Для наглядности закона изменения напряжений в заделке построим график (рис. 14). На­ чало координат поместим на вертикали, проведенной через свободный конец стержня, и ось абсцисс напра­ вим справа налево, параллельно оси стержня.. По оси абсцисс будем откладывать h, т. е. длину разрушенной части стержня, а по вертикали — сгМакс> возникающую в заделке при данной U. Каждая ордината построенного таким образом графика равна тому напряжению, кото­ рое возникает в заделке, если длина разрушенной части стержня равна абсциссе этой точки.

Этот график (см. рис. 14) характеризует изменение нормальных напряжений в данном сечении стержня в зависимости от величины разрушенной части, поэтому его будем называть «характеристической линией» для напряжений. Для каждого поперечного сечения стерж­ ня можно построить такую линию, поэтому напряжен­ ное состояние стержня будет описано семейством ха­ рактеристических линий; огибающая этого семейства даст представление о тех изменениях, которые возни­ кают в напряженном состоянии стержня при данном размере воронки. В данном случае огибающая этого се­ мейства сливается с характеристикой, построенной для сечения в заделке.

4.5.Определение частоты численным способом

Для вычисления основной частоты колебаний балки, имеющей воронку и изображенной на рис. 11, использу­ ем общую формулу (2.28)

Для функций X принимаем синусоиду

Х л = sin-"JE. и x ; =

- ^ - s i

ппх

г

I

Is

Жесткость EJ(x) и масса ц, (д:) балки являются те­ перь функцией х и изменяются по длине пролета.

Подставляя эти значения в формулу для частоты, ис­ пользуем конечные элементы и, заменяя интегралы сум­ мами, найдем

 

 

2

J M .

1

EJо п* л 4

(-т У - Л - \ ,

to

I*

S I ™” )1 *

J

 

 

где / о и Цо — момент инерции и

пая масса балки без учета во -

ронки.

 

 

 

 

 

Если обозначить численный коэффициент Кп, то по­

лучим

 

 

 

 

 

 

 

СОл = Кп <*>л0,

 

 

где сопо — частота балки с

постоянной жесткостью

EJ0 и погонной

массой ц.

 

 

 

 

 

Коэффициент

Кп удобно

вычислять

в

табличной']

форме

В табл. 2 дан пример вычисления коэффициента Ki для балки, имеющей симметричное ослабление в сере­

дине

пролета

глубиной 0,3

h на

длине

0,6

/.

 

 

 

 

 

 

 

Т а б л и ц а

2

 

sin —

sin* **

 

 

J ,

. ЯХ

JLsin»

**

хП

J/Jo

 

---- sin*-----

1

l

 

J,

l

Pa

l

 

 

 

0,1

0,309

0,095

,

1

0,095

0,095

0,2

0,588

0,345

1

1

0,345

0,345

0,3

0,809

0,655

0,343

0,700

0,225

0,458

0,4

0,950

0,904

0,343

0,700

0,312

0,633

0,5

1,000

1,000

0,343

0,700

0,343-0,5* 0,700-0,5-

 

 

 

 

2

1,148

1,881

* Коэффициент 0,5 введен в результате учета симметрии.

Из табл. 2 получаем /С= 1,148/1,881=0,61;

=

=0,61 wj.

 

Значение К для разных размеров воронки указано в табл.3

 

 

 

 

Т а б л и ц а 3

Глубина ослабле-

 

Ширина ослабления

 

 

 

 

 

ния

0,2/

0.4/

0.6/

0.8/

 

0,1- Н

0,94

0,88

0,84

0,81

0,3 -Л

0,87

0,72

0,61

0,51

0,5-Л

0,86

0,65

0,43

0,30

4.6.Колебания, вызванные начальной скоростью при ступенчатом изменении частоты

Рассмотрим случай, когда системе с одной степенью свободы в начальный момент сообщили ско­ рость y 'Q, затем эта система, предоставленная самой се­

бе, совершает колебания и изменяет в процессе движе-

to

V

2

 

 

-

 

1

 

0 ОМQ61 1,2 1.6 * <*>п

 

----- Внезапноеизменениечастоты

ния свою частоту согласно графику на рис. 15. Приме­ ним для решения пошаговый метод.

В течение первого интервала времени от 0 до т про­ гиб и скорость вычисляют по формулам:

у1= — Sin COj / и у\ = уо cos (й± t.

®1

Обозначим уо= yj/coj, тогда при/=T i

Уъ = У0sin со, т, и у^ = у0 ®кcos o)t т^.

В этот момент свойства системы изменяются, и но­ вая система будет иметь начальный прогиб y Xt и на­

чальную скорость y'Xi. В пределах второго интервала времени тг прогиб и скорость определяются по формулам:

Ух, =

yXisin (о21+

cos CD21,

y'Xt =

+ yXl©2C0S02t ~~ Ух, sin Ш2*•

Подставляя вместо t значение тг, найдем прогиб и ско­ рость в конце второго интервала. Так же можно перей­ ти к третьему интервалу и т. д. В конце п-го интервала времени получим:

Уп = Уп- 1 sin v +

©Л cos <»„т„,

Уп = *л-1 “ л C0S ®п Тл -

Уп- 1sin *>« V

Если наибольший прогиб будет достигнут в преде­ лах п-го интервала времени, то он будет вычисляться по формуле

Уш“с = VУп-1 + (“ ^ “ ) *

После подстановки в формулу для n-го интервала значе­ ний п-1-ro и т. д. прогиб и скорость п-го интервала бу­ дут выражены через начальную скорость

у„ = у0{2>)- »;= »o Mi R -

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

На рис. 16 дан график изменения коэффициента т], на который следует умножить наибольший прогиб, вы­ численный для системы с постоянной массой и жест­ костью для того, чтобы получить наибольший прогиб для системы с изменившейся частотой. На графике по оси абсцисс отложено отношение двух частот. В чис­ лителе стоит частота системы, вычисленная для началь­ ного момента времени, в знаменателе— частота систе­ мы с учетом изменений. Переход же от первой частоты к последней совершается непрерывно в виде линейной функции от времени. Кривая, показанная сплошной ди-

йией, относится к случаю, когда процесс изменения свойств системы (массы и жесткости) продолжается в течение одной десятой доли периода свободных колеба­ ний. За этот отрезок времени частота системы изменяет­ ся от G>i ДО ©п-

Для сравнения на том же графике пунктиром пока­ зано изменение коэффициента TJ для случая, когда про­ исходит мгновенное изменение свойств системы, т. е. система с частотой ©i в момент времени п переходит сразу в систему, имеющую частоту ©п. Как и следовало ожидать, при отношении частот, равном единице, обе кривые пересекаются. Весьма интересно отметить, что при <Bi/d)n<l кривая, которая относится к системе, по­ степенно изменяющей свои свойства, идет выше пунк­ тирной кривой. Это указывает на то, что наибольший прогиб системы, постепенно изменяющей свои свойства, будет больше такого же прогиба системы с быстрым из­ менением этих свойств. Если же ©i/©n> l , то получа­ ется обратная картина.

На рис. 17 показано изменение прогиба с течением времени для случая, когда в результате изменения свойств системы происходит увеличение частоты. В точках, где меняется частота, кривая, показывающая изменение прогиба, опускается на графике ниже. Наи­ больший прогиб получается меньше того, который был бы у системы с постоянной частотой. Если бы в процес­ се движения происходило постепенное уменьшение час­ тоты системы, то максимальный прогиб был бы больше.

4.7. Внезапное нагружение

Пусть внешняя сила Р была внезапно приложена к

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

/ = П, затем

она прекратила свое действие, и в момент

времени

то изменились свойства системы, т. е. ее мас­

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

на рис. 18. Для момента времени

t, удовлетворяющего

неравенству

то, можно написать

 

 

у = ух%cos ©21 -f

sin ©21

и для всех последующих интервалов

уп = Ух cos at + Утд 1 sin <ап t.

П 1

(Вд

В этой формуле г/То

и у^ — прогиб и скорость систе­

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

обычные формулы, но сохранив в них только то слагае­ мое, которое зависит от внешних сил, и отбросив сла­ гаемые, относящиеся к начальным условиям, так как в данном случае при /= 0 , Уо=0 и у'0 = 0

Для скорости получим формулу

У \~ -

sin

COS ©!^т0----J- j

Если вместо силы ввести в формулу импульс S| = =PiTi и обозначить

Уо = OOi/П! CDJTI /2

то формулы примут вид:

Ух. = Уоsin“1(•'о--1-)*?;. = У„«1 cos (Т0 - -S-) -

Чтобы оценить влияние продолжительности интервалов времени т0 и на результат, ограничимся несколькими ступенями изменения частоты и подставим полученные значения в формулу (с. 58)

Уп = Уоj 2 si и—1(V -i — cos u j-h

+ ^ t cos^ - ‘{rn- 1~ T -)sin““']-

Наибольший прогиб t/макс= г)Уо. Обозначим т)— чис­ ленный коэффициент. Тогда т| будет характеризовать изменение наибольшего прогиба системы при разных •соотношениях частот и интервалов времени. Сначала рассмотрим предельный случай, когда a>i = 002 H T 0= T I , ■тогда коэффициент rj= 1 и наибольший прогиб равен лрогибу системы с постоянной массой и жесткостью, как этого и следовало ожидать. Исследуем изменение вели­ чины г]; для этого выразим частоту через период коле­ баний. Тогда получим

Пользуясь этой формулой, составим табл. 4 измене­ ния -л в функции от соотношений: ©1/(02и TI/TY

При помощи табл. 4 составлен график (рис. 19). По оси абсцисс отложено отношение частот, по оси орди-

Т а б л н ц а 4

то/т,=2

Ш,/Оа

*1/Тг

0

0

0,445

0,707

1,0

2,0

0,10

0,456

0,603

0,779

1,0

1,8

0,25

0,925

0,936

0,957

1,0

1,2

1,00

1,00

1,00

1,00

1,0

1,0

нат — значение т]. Каждая кривая относится к опреде­ ленному отношению Т1/Г1

Соседние файлы в папке книги