- •2. Определение амплитуды собственных колебаний и графическое изображение собственных форм
- •3. Проверка ортогональности собственных форм колебаний
- •5. Определение амплитудных значений инерционных сил
- •6. Определение эпюры изгибающих моментов от действия собственного веса вибраторов и амплитудных значений изгибающих моментов при вынужденном стационарном режиме колебания рамы
- •7. Построить эпюру моментов при одновременном действии статических и динамических сил и определить положение опасного сечения конструкции
- •8. Определение максимального напряжения и проверка условий прочности в наиболее опасном сечении
- •5.7. Сейсмические колебания системы с конечным числом степеней свободы
8. Определение максимального напряжения и проверка условий прочности в наиболее опасном сечении
кН/м2 = 53.2МПа <R = 190 МПа.
Следовательно, условие прочности рассматриваемой рамы обеспечено.
5.7. Сейсмические колебания системы с конечным числом степеней свободы
В теории расчета сооружений на сейсмические воздействия (теория сейсмичности), как и в других областях динамики различных механических систем, обычно применяются расчетные схемы с распределенными и дискретными параметрами (массами). Система с дискретными параметрами хотя и носит приближенный характер, но более универсальна и можно получить решение для системы любой сложности, вследствие чего наиболее часто применяются в инженерных расчетах.
Для получения динамичных расчетных схем в виде системы с конечным числом степеней свободы, фактическая распределенная масса система концентрируется в определенных местах в виде материальных точек. В итоге получается, невесомая система, несущая определенное количество сосредоточенных масс. Число степеней свободы система равно числу независимых геометрических параметров, однозначно определяющих положение сосредоточенных масс в произвольном моменте времени.
Рис. 5.13
В качестве примера рассмотрим методику расчета здания, имеющего n этажей на сейсмическом воздействии. Сконцентрировав массу здания на уровнях перекрытия и фундаментной плиты получим, систему в виде консольного стержня жестко заделанной в фундаментной плите, лежащей в условиях полного прилипания на поверхности упругого инерционного основания (рис. 5.13).
Будем рассматривать поперечные колебания стержня в плоскости (zy). Начало системы координат поместим в центре тяжести подошвы фундамента сооружения. Жесткость стержня по высоте изменяется по произвольному закону. На характер деформаций стержня не накладывается никаких ограничений, кроме требования линейной деформируемости.
Положение системы в произвольный момент времени t > 0 оп
ределяется линейными горизонтальными смещениями (yi + y0) (i = 1,2,..., n + 1) (рис. 5.13).
Так как y0(t) есть перемещение грунтов основания при землетрясении на свободной поверхности земли, в предположении отсутствия сооружения, то оно здесь принимается заранее заданной величиной. Следовательно, если нам удастся определить величины yi (t) (i = 1,2,..., n + 1), мы через значения этих величин в произвольный момент времени можем определить положение заданной системы.
Отсюда следует, что рассматриваемая система, располагая (n + + 1) количеством сосредоточенных масс, имеет (n + 1) степеней свободы.
Колебания линейной системы при заданном внешнем кинематическом воздействии y0(t) полностью определяется ее инерционными и деформативными свойствами и параметрами рассеивания энергии. Инерционные свойства рассматриваемой системы характеризуются сосредоточенными массами mi (i = 1,2,..., n + 1) и характером их распределения по высоте. Деформативные свойства системы могут быть охарактеризованы при помощи единичных перемещений dik (i, k =1,2,..., n + 1) представляют собой горизонтальное перемещение точек i от действия единичной горизонтальной силы, приложенной в точке k. Перемещение dik в рамках принятой расчетной схемы определяется:
, (5.48)
где - горизонтальные перемещения точки i от действия единичной горизонтальной силы, приложенной в точке k, обусловленные соответственно: деформациями конструктивных элементов здания; относительным сдвигом между подошвой фундаментной плиты и основанием; поворотом подошвы фундаментной плиты относительно основания.
Выражение dik можно записать в следующем виде:
, (i, k = 1,2,..., n + 1), (5.49)
Так как, фундаментная плита считается абсолютно жестким, поэтому при i = n + 1, или k = n + 1 следует принимать .
Здесь определяется по формуле Мора; с1 и сj - являются коэффициентами квазистатической жесткости основания при равномерном сдвиге и неравномерном сжатии или растяжении и их значения можно определить по следующим соотношениям а автора
[7]:
, (5.50)
где приняты следующие обозначения: - скорость распространения поперечных волн в грунтах основания; r - плотность грунтов основания; F - площадь подошвы фундаментной плиты; - момент инерции площади подошвы фундаментной плиты относительно оси х.
Для учета рассеивания энергии при колебаниях системы воспользуемся теорией Фойгта, согласно которой диссипативные силы прикладываются к сосредоточенным массам в состоянии движения системы, величина которых пропорциональна скорости движения сосредоточенных масс. Коэффициенты пропорциональности для рассматриваемой системы определяются по формуле автора [7]:
, i, k = 1,2,..., n + 1. (5.51)
Величина , где ; d - логарифмический декремент колебания, характеризует рассеивания энергии по корректированной гипотезе Фойгта за счет внутреннего неупругого сопротивления материалов конструкций при их деформации; h1 - характеризует излучение энергии в основании за счет сдвиговых деформаций, происходящих на контактной поверхности между фундаментной плитой и основанием; hj - коэффициент рассеивания энергии за счет неравномерных линейных деформаций, происходящих на контактной поверхности между фундаментной плитой и основанием.
Акустическое сопротивление основания при равномерном сдвиге h1 и неравномерного сжатия и растяжения hj определяются по известным соотношениям и автора [7]:
, (5.52)
где - скорость распространения продольных волн в грунтовом основании.
Воспользуемся методом сил и запишем величину перемещения yi (t) произвольной массы с номером i = 1,2,..., n + 1, от действия сил инерции и сил учитывающих рассеивание энергии в рассматриваемой системе:
,
(i = 1,2,..., n + 1) (5.53)
Здесь Ik(t) (k-1,2,...,n+1) сила инерции, действующая на k-ую массу и определяется по принципу Даламбера:
. (5.54)
Сила сопротивления , возникающая в k-ой массе, согласно гипотезе Фойхта прямо пропорциональна величине скорости его движения:
. (5.55)
Подставляя выражения (5.54) и (5.55) в (5.53) и после некоторых преобразований получим: дифференциальное уравнение движения заданной системы в следующем виде:
, (i = 1,2,..., n + 1). (5.56)
Для расчета сооружений на сейсмические воздействия справедливы нулевые начальные условия, т.е. предполагается что до начала землетрясения сооружение находится в состоянии покоя. При землетрясении сооружение, переходя в движение, ее состояния характеризуется системой уравнений (5.56).
Для расчета системы дифференциальных уравнений (5.56) применяется метод преобразования Лапласа, т.е. искомые функции находятся по формуле:
, (5.57)
где yi(s) является изображением функции yi(t) по Лапласу и определяется по формуле:
. (5.58)
Подставляя (5.57) в (5.56) с учетом нулевых начальных условий задачи, получим:
, (i = 1,2,..., n + 1). (5.59)
Последнее представляет систему алгебраических уравнений относительно перемещений yi(s) в изображениях Лапласа.
Решение (5.59) в изображениях представлено в виде:
, (5.60)
где Di(s) - представляет собой определитель системы неоднородных алгебраических уравнений (5.59); D(s) - определитель той же системы при неизвестных yi(s).
Применяя к выражению (5.60) операции обратного преобразования Лапласа с применением теоремы свертки, получим решение задачи в следующем виде:
, (5.61)
где ak корни определяются а - дифференциал определителя по ak при s = ak, т.е.
В традиционных методах расчета сооружения на сейсмостойкость, как правило, применяется следующее упрощающее допущение, что основание сооружения является абсолютно твердым телом, т.е. с1®¥ и сj®¥. Если исходить из условия существования полного прилипания между фундаментной плитой и основанием на их контактной поверхности, очевидно, что масса с номером n + 1, т.е. фундаментная плита, полностью повторяет закон движения основания. С другой стороны, так как закон движения основания в данном случае считается исходной известной функцией, следовательно закон движения фундаментной плиты тоже следует считать известной величиной. Поэтому число степеней свободы рассматриваемой системы (см. рис. 5.13) на одну единицу уменьшается и принимает значения равной n.
Искомыми величинами в данном случае являются перемещения сосредоточенных масс с номерами i = 1,2,...,n.
С учетом данного обстоятельства уравнение движения сооружения из (5.56) упрощается и принимает вид:
, (i = 1,2,...,n) (5.62)
Для решения системы дифференциальных уравнений (5.62) с постоянными коэффициентами применяется метод разложения колебаний по формам, основанный на методе разделения переменных, т.е.
. (5.63)
Сначала, для определения собственной частоты и собственного вектора Xiv, рассматриваются собственные колебания системы без учета сил сопротивления. В данном случае из (5.62) получим уравнения движения системы без учета сил сопротивления в свободном режиме колебаний:
= 0, (i = 1,2,...,n). (5.64)
Подставляя решение (5.63) в (5.64), с учетом условий ортогональности собственных форм колебаний, т.е.
= 0, (i ¹ k, i , k = 1,2,...,n), (5.65)
и после ряда преобразований получим:
. (5.66)
Выполнение этих равенств для произвольного значения t возможно лишь в том случае, если каждая из них в отдельности равна одной и той же постоянной при любом значении v. Обозначив эту постоянную через , получим:
= 0, (i = 1,2,...,n). (5.67)
Последние уравнения представляют собой систему n линейных однородных алгебраических уравнений относительно неизвестных Xin для каждой n = 1,2...,n формы колебаний.
Для определения ненулевых решений системы (5.67) необходимо обеспечить равенство нулю ее детерминанта:
. (5.68)
В развернутом виде это выражение представляет собой алгебраическое уравнение n-ой степени относительно . Корни этого уравнения вещественны, положительны и в большинстве случаев отличны друг от друга. Таким образом из решения (5.68) определяется n положительных значений w0, которые в возрастающем порядке: w01 > w02 >...> wn являются собственными частотами системы без учета ее диссипативных свойств.
После определения собственных частот из решения системы (5.67) определяются значения собственных векторов Хiv (i = 1,2,..., n) для каждой v-ой формы колебаний.
Для определения решения уравнений движения системы в вынужденном режиме колебаний подставим выражение (5.63) в уравнения (5.62) и с учетом условий ортогональности (5.65), получим:
, (5.69)
где приняты следующие обозначения:
. (5.70)
При нулевых начальных условиях, решая дифференциальное уравнение (5.69) и подставляя в (5.63), окончательно получим:
, (5.71)
где - называется частотой собственных колебаний с учетом диссипативных свойств системы.
Как показывает уравнение (5.69), применение принципа разложения колебаний по собственным формам, позволяет рассматривать колебания линейно-деформируемых систем по отдельным формам независимо от колебаний по другим формам, вследствии чего системы со многими степенями свободы рассчитываются как системы с одной степенью свободы для каждой отдельной формы.
Указанный подход позволяет при рассмотрении системы со многими степенями свободы оценить динамический эффект внешнего воздействия через значения коэффициента динамичности для системы с одной степенью свободы.
Сейсмические колебания системы с одной степенью свободы из (5.54) принимают вид:
, (5.72)
откуда
. (5.73)
Последнее уравнение можно записать в виде:
, (5.74)
где .
Решение последнего уравнения имеет вид:
. (5.75)
Дважды дифференцируя последнее выражение получим относительное ускорение , после суммирования с получим последнее ускорение системы с одной степенью свободы в следующем виде:
, (5.76)
где
. (5.77)
Коэффициент динамичности в данном случае определяется по формуле:
. (5.78)
Примерный вид графика ускорения колебания грунтов при сейсмических воздействиях представлен на рис. 5.14.