книги / Метод конечных элементов в динамике сооружений
..pdfiio ним описать огибающую, которая й позволит судить о величине напряжений и о положении опасного сечения. Для пояснения изложенного рассмотрим простейшую за дачу одного измерения.
Стержень постоянного поперечного сечения (рис. 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 т^.
5Д
В этот момент свойства системы изменяются, и но вая система будет иметь начальный прогиб 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