авторефераты диссертаций БЕСПЛАТНАЯ РОССИЙСКАЯ БИБЛИОТЕКА - WWW.DISLIB.RU

АВТОРЕФЕРАТЫ, ДИССЕРТАЦИИ, МОНОГРАФИИ, НАУЧНЫЕ СТАТЬИ, КНИГИ

 
<< ГЛАВНАЯ
АГРОИНЖЕНЕРИЯ
АСТРОНОМИЯ
БЕЗОПАСНОСТЬ
БИОЛОГИЯ
ЗЕМЛЯ
ИНФОРМАТИКА
ИСКУССТВОВЕДЕНИЕ
ИСТОРИЯ
КУЛЬТУРОЛОГИЯ
МАШИНОСТРОЕНИЕ
МЕДИЦИНА
МЕТАЛЛУРГИЯ
МЕХАНИКА
ПЕДАГОГИКА
ПОЛИТИКА
ПРИБОРОСТРОЕНИЕ
ПРОДОВОЛЬСТВИЕ
ПСИХОЛОГИЯ
РАДИОТЕХНИКА
СЕЛЬСКОЕ ХОЗЯЙСТВО
СОЦИОЛОГИЯ
СТРОИТЕЛЬСТВО
ТЕХНИЧЕСКИЕ НАУКИ
ТРАНСПОРТ
ФАРМАЦЕВТИКА
ФИЗИКА
ФИЗИОЛОГИЯ
ФИЛОЛОГИЯ
ФИЛОСОФИЯ
ХИМИЯ
ЭКОНОМИКА
ЭЛЕКТРОТЕХНИКА
ЭНЕРГЕТИКА
ЮРИСПРУДЕНЦИЯ
ЯЗЫКОЗНАНИЕ
РАЗНОЕ
КОНТАКТЫ


Pages:     | 1 | 2 || 4 |

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

-- [ Страница 3 ] --

Уравнение перемещений в (14) - наиболее общая матричная форма записи
интеграла Дюамеля для ДДС. Подчеркнута отличительная особенность этого уравнения от известных интегралов Дюамеля, состоящая в том, что оно не требует построения ИПФ -. наиболее трудоемкой части анализа. Выражение подынтегральной матричной функции H(t-) = 2Rе {Ф(t-)U-1} записано в простой математической форме и содержит величины U, Ф(t-), вычисление которых основано на решении МКУ, не прибегая к его спектральному разложению. Известный аналог матрицы H(t-) (матрица функций Грина) в общем случае не имеет аналитического представления. Это является сильным препятствием при определении динамической реакции на основе ИПФ

В простейших случаях, важных для приложений строительной механики, из (14) при (t0 = 0) Y0 = 0 = 0 получены вычислительные формулы для интеграла Дюамеля (рис. 4, 5).

Внезапно приложенные силы постоянной величины, исчезающие при t = tj: Р(t) = Р0, где Р0 = [p0j] (j = 1,…,n) (рис. 4, а). После интегрирования в (14) уравнение реакции ДДС на активном участке нагружения принимает вид (tt1):

где E - единичная матрица. При tt1 система совершает свободные колебания под действием начальных условий:

Вибрационная нагрузка Р(t) = sin (t+ )Р0. Здесь = diag (1,…, n), = diag (1,…, n) - диагональные матрицы угловых частот и начальных фаз вибрационных сил, P0 = [р0j] (j = 1,..., n) - вектор амплитуд возмущающих сил. На рис. 4, б показаны параметры нагрузки, действующей в j-м узле.

Полная реакция системы определена векторами перемещений и скоростей

Вычисление (15) связано с решением непрерывного уравнения Сльвестра

Исследованы случаи разрешимости уравнения (16). При полной диссипации ДДС (det C > 0) решение уравнения (16) всегда единственно. При колебаниях консервативной системы (С =0) существуют условия для неоднозначного решения. Эго происходит при совпадении частот собственных и вынужденных колебаний (резонанс). Показано, что (16) эквивалентно n векторным уравнениям:

где Ik (t), Fk(t) - k-е столбцы матриц I(t), F(t).

Синусоидальный импульс P(t) = sin (t)P0, где = E/t1; P0 = [p0j] (j = 1,…,n) (рис. 4, в). Реакция ДДС на активном участке (tt1) вычисляется в соответствии с (15), где матричная функция I(t) определена при = 0; I(t) = [(S)2 + 2]-1 F(t). Реакция системы при tt1 выражена уравнениями:

где Z(t) = [Ф(t) + Ф(t-t1)][U(S2 + 2)]-1.

Периодические импульсы. На рис. 5 показаны импульсы сил, действующие в j-м узле конструкция. Рассмотрено действие периодических импульсов прямоугольной и синусоидальной форм постоянной длины (t0= t’i – ti-1), периодичности (T = ti-ti-1) и амплитуды p0j.

Вычисление параметров реакции ДДС от действия последующих импульсов обеспечивается на основе информации (в виде начальных условий: Y0, 0) о кинематических характеристиках узлов системы, вызываемых предыдущим импульсным воздействием. Получена система уравнений, определяющих реакцию при вынужденных колебаниях от действия i-й группы импульсов (ti-1tt’i):

Для случая прямоугольной и синусоидальной форм импульсов матричная функция (t) представлена соответствующими выражениями:

После исчезновения i-й группы импульсов ДДС совершает свободные колебания на интервале времени (t’itti):

под действием начальных условий, назначаемых на основе (17) в конце предыдущего интервала времени: Y0=Y(t’i), 0=(t’i)

Приведен анализ реакции каркасного здания, изображенного на рис. 2, на действие периодических синусоидальных импульсов и вибрационной нагрузки. В соответствии с РДМ (рис. 2, б), вектор Y(t) (17) имеет следующую структуру:

Y(t) = [xi(t), yi(t), x2(t), y2(t), x3(t), y3(t), 1(t), 2(t), 3(t)]T где xi(t), yi(t) - поступательные перемещения центра тяжести перекрытия i-го этажа вдоль осей х и y соответственно; i(t)- угол поворота перекрытия вокруг центра жесткости упругих связей i-го этажа.

Векторы сил действуют в уровнях перекрытий под углом к оси х (рис. 2. а). Амплитуды импульсного воздействия на каркас вычислялись исходя из нормативного значения ветрового давления на поверхность здания, равного q = 3.5 105 кН/см2. Для вектора амплитуд Р0 = [F0, M0]T при = 45° имеем F0 = [33, 33, 22, 22, 10, 10] кH; M0 = [13608, 0, 0] кHсм. При длине импульсов t = 0.15 с рассмотрено действие на каркас одиночных ударов и периодических импульсов с периодичностью, равной половине периода: T = 0.5 T1 = 0.3332 с и периоду основного тона колебаний T= T1 = 0.6663 с, где T1 = 2 /1, 1 = 9.4298 с-1.

На рис. 6-9 приведены осциллограммы параметров динамической реакции каркаса при периодичности импульсов T1/2. Перемещения (рис. 6) и скорости (рис. 7) центров тяжести перекрытий на осциллограммах представлены обеими линейными составляющими в направлении осей х (a), у (б) и угловой составляющей (в); восстанавливающие (рис. 8) и диссипативные (рис. 9) силы, действующие в перекрытиях этажей, - линейной составляющей по оси х (а) и угловой составляющей (б). Сравнительный анализ реакции здания оценивался с помощью модели А.И. Цейтлина C = M (пунктир).

Проведен анализ наиболее загруженных колонн каркаса при варьировании ряда параметров периодического импульса. При циклическом изменении параметров ta и строились поверхности максимальных нормальных напряжений в зависимости от периодичности импульсов. Число повторений импульсов во всех случаях принималось, равным 5.

На рис. 10, а приведена поверхность нормальных напряжений, построенная на сетке из 45х33 узлов при периодичности импульсов T = 2ta. Интервалы варьируемых величин: ta [0.02, 0.9] с (при шаге ta= 0.02), [0, /2] рад (при шаге = /64 рад). Шаг временного анализа на активном участке составлял 0.01 с, на участке свободных колебаний - 0.016 с. Рассматриваемый режим нагружения (при T = 2ta) характерен для случая действия ветровой нагрузки, пульсационная составляющая которой может быть моделирована в виде периодического импульса. Поверхность нормальных напряжений, построенная при периодичности импульсов T = 0.3332 с (рис. 10, 6), имеет следующие характеристики. Сетка содержит 21х17 узлов, параметры нагружения: ta [0.01, 0.31] с; [0, /2] рад. Шаг временного анализа - 0.01 c.

Вибрационное воздействие на каркас осуществлялось с помощью двух сил Fi(t) = F0 sin(it + i) (i = 1, 3), действующих в уровнях перекрытий 1-го и 3-го этажей. Причем вектор силы Fi(t) совпадает с направлением оси у и приложен к центру тяжести перекрытия 1-го этажа, а вектор F3(t) совпадает с направлением оси х, действуя по линии i-j на расстоянии l от центра тяжести перекрытия 3-го этажа (рис. 2,а). Параметры нагрузки:

F01 = 15 кН, 1 = 96 с-1, 1 = 0; F03 = 20 кH, 3 = 120 c-1, 3 0.

Вектор амплитуд P0 = [F0, M0]T имеет следующие значения: F0 = [0, 15, 0, 0, 20, 0] кН, М0 = [500, 0, 12000] кHсм. Моментная составляющая M01 = 500 кHсм вектора М0, действующая в перекрытии 1-го этажа, получена вследствие несовпадения положения центра тяжести С1 перекрытия с центром жесткости O1 упругих связей (рис. 2, а).

Временной анализ реакции каркаса проводился при изменении фазы 3 силы F3 в интервале [0, 2} с шагом = /36 рад (5°) относительно нулевой начальной фазы 1 нагрузки F1. По результатам анализа наибольшие отклонения в максимальных значениях нормальных напряжений и относительных перемещений концов наиболее нагруженной колонны каркаса соответствуют значениям 3 = 1,484 рад (z = 0.084 МПа: 3-й этаж, колонна № 1; = 0.0224 мм: 2-й этаж. колонна № 13) и 3 = 3,142 рад (z = 0.181 Мпа, = 0,049 мм: 3-й этаж, колонна № 1) и отличаются друг от друга более, чем а 2 раза.

На рис. 11 приведены осциллограммы линейной составляющей параметров реакции, действующих вдоль оси х: перемещений (а), скоростей (б), восстанавливающих (в) и диссипативных (г) сил для значения 3 =3.142 рад.

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

Внешняя нагрузка представлена в виде вектора

где f(t) - безразмерная скалярная функция времени t; P0 = [р0j] (j= 1,…, n) - вектор постоянных усилий. На основе (18) введен вектор импульсов сил

Доказательство соотношений взаимности в ДДС базируется на двух положениях. Одно - связано с формой записи систем разрешающих уравнений (14), в которых величины Y0, 0 для удобства приняты нулевыми. Показано, что эти системы обладают симметричной структурой:

Осуществлен переход (путем обращения матриц D(t) и V(t) к матрицам динамической жесткости L(t) и импульсов H(t) также симметричного вида):

Эти результаты можно считать расширением известных результатов (теорема Максвелла о взаимности перемещений: D(t) = D(t)T и теоремы Рэлея о взаимности реакций: L(t) = L(t)T и импульсов: H(t) = H(t)T для случая нестационарного процесса, протекающего в ДДС. Соотношение V(t) = V(t)T, по своей сути, есть теорема о взаимности скоростей масс от действия единичных импульсов, прикладываемых к узлам ДДС. В литературных источниках данный закон не выявлен, хотя не исключено, что для частных случаев задачи он известен.

Второе положение относится к алгебраической трактовке принципа взаимности, впервые данной в 1927 г. проф. П.Л. Пастернаком. Согласно этому положению. свойство взаимности присуще любой системе n линейных уравнений с n неизвестными, обладающей симметричной структурой коэффициентов.

На основании изложенного получены соотношения взаимности:

выражающие теорему взаимности Бетти в форме произведения перемещений и сил (первое соотношение) и в форме произведения скоростей и импульсов сил (второе соотношение) в произвольной упругой ДДС. Векторы P(t)’, Y(t)’, (t)’ и

Z(t)’ представляют новые системы соответственно сил, перемещений, скоростей и импульсов в исходной ДДС

Соотношения (19) выполняются для любой системы сил, определяемой вектором нагрузки (18), и являются обобщением результатов Рэлея, доказавшим первый закон взаимности в (19) со всеми его следствиями в ДДС для частного случая системы сил, гармонически изменяющейся во времени (то есть при условии f(t) = sin(t + )). Второй закон в (19) и его частный случай H(t) = H(t)T, были доказаны Рэлеем для консервативной системы, находящейся под действием мгновенных импульсов.

Результаты обобщены в виде теоремы о предпосылках закона взаимности в произвольной упругой ДДС: Пусть характер воздействия динамической нагрузки в узлах ДДС определяется вектором (18). Тогда, если матрицы М, С, К дифференциального оператора уравнения движения (7) обладают симметрией, то: (а) полная система уравнений динамической реакции (14) также имеет симметричную структуру; (б) к данной упругой системе применимы законы взаимности как в форме общих (19), так и частных теорем.

Показано, что динамическая реакция Y(t) = D(t)P(t) выражается через её статическую составляющую Yст = K-1 P0 посредством матричной функции

характеризующей учет динамического эффекта в произвольной конечномерной системе при нестационарных воздействиях, вследствие чего Y(t) = (t)Yст

Для матриц (t), D(t), L(t), V(t), H(t) приведены конечные формулы для случая действия внезапно приложенной нагрузки.

Последующие главы диссертации, с шестой по восьмую, посвящены упругопластическому анализу ДДС при действии кратковременной нагрузки.

В шестой главе «Теоремы состояний и анализ внутренних динамических параметров системы» предложены математические модели упругопластического расчета и доказаны теоремы, характеризующие качественные уровни состояний конструкции в процессе ее неупругого деформирования.

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

t [ti, ti+1] (I = 0, 1,…), на которых динамические параметры системы неизменны. Время ti соответствует открытию или закрытию шарниров пластичности. Это позволило обобщить временной анализ ДДС на случай движения конструкции с неупругой восстанавливающей силой, используя для этой цели метод временного анализа, разработанный для упругой системы.

Условия динамического равновесия ДДС с идеально упругопластическими восстанавливающими силами (вектор R(t), см. рис. 12) представлены в виде

Математические модели расчета включают в себя физические соотношения

и комплекс условий: упругости, текучести в j-м элементе конструкции при t = tm и разгрузки в том же элементе при t = ti (ti > tm):

Здесь K(ti)Y(t) - квазиупругая составляющая вектора (21); K(ti), K(j) - матрицы жесткости квазиупругой системы и j-го элемента конструкции. Составляющие вектора (21): R0(ti) - вектор предельных значений и R*(ti) - вектор остаточных усилий определяются в упругопластических пружинах при текучести и разгрузке соответственно; Y*(ti) - вектор остаточных перемещений ДДС, значения которого вызваны пластическими деформациями в j-м элементе конструкции.

Упругие колебания происходят при условии, когда вдоль всех степеней свободы ДДС значения относительных перемещений не превышают их предельно упругих значений y0j (j = 1,…,n)

Сформулирована задача неупругого анализа ДДС, в узлах которой действует нагрузка P(t) = [pj(t)] (j = 1,…,n) (рис. 13). Для вычисления динамической реакции системы на любом i-м интервале времени t [ti, ti+1] необходимо удовлетворить уравнению движения (20) физическими соотношениями (21) так, чтобы выполнялись условия упругости, текучести и разгрузки (22).

Далее проведен анализ собственных колебаний квазиупругой системы на интервале t [ti, ti+1], требующий рассмотрения характеристического МКУ (8) и соотношений Виета в (13) при K = K(ti), где Sk и Sl, из обшей корневой пары (9).

Рост текучести (tt1) вызывает снижение коэффициентов жесткости ДДС и, вследствие этого, изменение внутренних динамических параметров. Этот факт отражен пятью теоремами состояний, устанавливающими критерии соответствия между определителями матриц в равенстве (13): SkT MSl = K(ti)(kl).

Теорема 6.1 (об условии невырожденного состояния квазиупругой системы). Матрица жесткости К(ti,) квазиупругой системы (i 0) невырожденна тогда и только тогда, когда невырожденны обе матрицы внутренних динамических параметров в корневой паре (9) характеристического МКУ, т. е

Теорема 6.2 (об условии вырожденного состояния квазиупругой системы). Пусть ДДС обладает полной диссипацией. Матрица жесткости К(ti) квазиупругой системы (i> 0) вырожденна на интервале t [ti, tm],тогда и только тогда, когда одна из матриц Sk в (9) вырожденна, а другая - нет, т. е.

Если пластические шарниры образуются во всех опасных сечениях конструкции (при t = tq [ti, tm]), то K(tq) = 0. Такое деформированное состояние квазиупругой системы названо предельно вырожденным состоянием (ПВС).

Показано, что в процессе пластического деформирования частотно-демпфированный спектр системы становится подвижным. Характер кривых собственных частот неупругой конструкции показан на графиках, отвечающих условиям теоремы 6.1 (на всем интервале реагирования (рис. 14)), теорем 6.2, 6.3 (на интервале t [tq, tq+1] (рис. 15)) и теорем 6.4, 6.5 (при t [tq, tq+1] (рис. 16))

Для спектральных норм матриц коэффициентов демпфирования и собственных частот в отдельных состояниях квазиупругой системы построены двухсторонние априорные оценки:

нижние (1, 1, 01 ) и верхние (2, 2, 02) грани спектральных норм ||G||, ||W|| и ||W0|| (W0 - матрица собственных частот соответствующей консервативной системы), равных максимальным значениям внутренних динамических параметров ДДС: ||G|| = m = max (1,…, n), ||W|| = n, ||W0|| = 0n

Оценки (23) соответствуют упругим колебаниям ДДС при пропорциональном демпфировании. Оценки (24) - движению упругопластической системы, состояние которой удовлетворяет условию теоремы 6.4 при t [tq, tq+1] [ti, tm]

На рис. 17, 18 дана графическая интерпретация двухсторонних оценок в зависимости от частотно-демпфированных уровней упругопластической ДДС. Показаны нижние (1(t), 1(t), 01(t)) и верхние (2(t), 2(t), 02(t)) грани норм ||G(t)||, ||W(t)|| и ||W0(t)|| на всем участке упругопластического нагружения, когда выполняются условия теоремы 6.1 (t [ti, tm]), теорем 6.2 и 6.3 (t [ti, tm]) (рис. 17), а также теорем 6.4 и 6.5 (ПВС при t [tq, tq+1]) (рис. 18).

Седьмая глава «Неупругий временной анализ: обобщение интеграла Дюамеля» посвящена разработке метода динамического расчета диссипативных конструкций за пределом упругости при кратковременном нестационарном воздействии, построению и реализации разрешающих уравнений неупругих колебаний ДДС при различных состояниях квазиупругой системы.

На основе предложенных математических моделей построен шаговый процесс, в котором упругопластический расчет сведен к последовательности квазиупругих решений на интервалах времени t [ti, ti+1] (i = 0, 1,…). В результате интегрирования уравнения движения (20), с учетом (21), получены уравнения полной реакции квазиупругой системы

Приведенный результат есть обобщение интеграла Дюамеля для физически нелинейной системы с идеально упругопластическим поведением материала. Уравнения реакции ДДС (25) обеспечивают получение замкнутого решения в рамках принятой модели деформирования. Первые два члена в уравнениях, стоящие под знаком суммы, выражают реакцию системы при свободных колебаниях. Последний член - при вынужденных колебаниях. При этом уравнения реакции при свободных колебаниях включают реакцию от упругопластических смещений узлов ДДС при текучести и разгрузке.

Решены вопросы практической реализации уравнений неупругих колебаний ДДС. В зависимости от условий состояния квазиупругой системы получены расчетные формулы для вычисления второго интеграла в (26). При выполнении условий теоремы состояния 6.1 его вычисление проводится по формуле

где YR(ti) = K(ti)-1 [-R0(ti) + R*(ti)] - вектор упругопластических смещений узлов ДДС, накопленных к моменту времени t = ti; 0 = Y0 – YR(ti) - новый вектор начальных условий

При условии вырожденного состояния системы (теорема состояния 6.2) для интеграла.Jk(t) в (26) получено выражение в виде функционального ряда

При выполнении условия теоремы состояния 6.3 (С = 0) показано, что параметры динамической реакции в уравнениях (25) становятся неопределенными, ввиду det U=0 и неограниченного возрастания величины U-1 .

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



Pages:     | 1 | 2 || 4 |
 





 
© 2013 www.dislib.ru - «Авторефераты диссертаций - бесплатно»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.