Р асчет зданий и сооружений на нестационарные воздействия

ЛИТЕРАТУРА:

1. Потапов а.Н. Динамический анализ дискретных диссипативных систем при нестационарных воздействиях. – Челябинск, 2003. – 167 с.

2. Чернов ю.Т. Вибрации строительных конструкций. – м.: асв, 2011. – 384 с.

3. Тексты лекций

Основные определения

Дискретная система. Это система с конечным числом степеней свободы

Консервативная система. Система, колебания которой происходят без учета сил сопротивления

Диссипативная система. Это система, при колебаниях которой учитываются силы сопротивления

Диссипация – рассеяние части энергии во внешнюю среду

Внешние динамические параметры системы: массы, характеристики жесткости, характеристики демпфирования (затухания), внешние силы

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

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

Нестационарные воздействия (процессы, нагрузки). Это такие воздействия, при которых происходят изменения внутренних динамических параметров конструкции

Лекция 1: УПРУГИЕ КОЛЕБАНИЯ ДИСКРЕТНЫХ СИСТЕМ

1. Дифференциальное уравнение движения дискретной диссипативной системы (ддс)

Решение задачи колебаний упругой ДДС требует рассмотрения матричного неоднородного дифференциального уравнения движения

M (t) + C (t) + KY(t) = P(t), (1)

и начальных условий: Y(t₀) = Y₀, (t₀) = ₀ (2)

где M = diag (m₁, ... , m_n), C = С^Т = [c_ij], K = K^T = [r_ij]  М_n (R) (i, j = 1, ... , n) – вещественные симметрические матрицы масс, демпфирования и жесткости;

Y(t) = [y_i (t)], P(t) = [p_i (t)]  M_n,1 (R) (i = 1, ... , n) – векторы перемещений и заданных нагрузок; n – число степеней свободы ДДС.

Уравнение (1) описывает колебания произвольной упругой ДДС в рамках теории вязкого трения (рис. 1).

Математически: ур. (1) – система неоднородных ОДУ 2-го порядка с постоянными коэффициентами (матрицы M, C, K – постоянные величины).

Основная система для построения матрицы жесткости K (рис. 2).

ОБЩЕЕ РЕШЕНИЕ ОДНОРОДНОГО ОДУ

Однородное дифференциальное уравнение с постоянными матрицами-коэффициентами, соответствующее уравнению (1)

M (t) + C (t) + KY(t) = 0, (3)

Общий интеграл однородного ОДУ

Y(t) = Re [(t) A], (4)

где (t)  фундаментальная матрица однородного ОДУ

(t) = e^St ; (5)

A – вектор постоянных интегрирования, зависящий от начальных условий.

Дифференцируем матричную функцию (5) и подставляем в левую часть уравнения (3) вместо соответствующих векторов

(t) = S(t), (t) = S²(t)

M S²(t) + C S(t) + K (t) = 0

[M S² + C S + K] (t) = 0

M S² + C S + K = 0

Лекция 2

Х АРАКТЕРИСТИЧЕСКОЕ

М АТРИЧНОЕ КВАДРАТНОЕ УРАВНЕНИЕ

И ЕГО АНАЛИЗ

Структура решения МКУ

Решение динамической задачи с учётом упруго-вязкого сопротивления связано с анализом матричного квадратного уравнения.

Пусть S  M_n(С) – некоторая квадратная матрица. Матрица (5):

(t) = e^St

является фундаментальной матрицей уравнения (3) тогда и только тогда, когда матрица S удовлетворяют характеристическому нелинейному уравнению, имеющему вид МКУ

MS² + CS + K = 0 (6)

где: M, C, K – матрицы масс, демпфирования и жёсткости:

M = M^T, C = C^T, K = K^T  M_n(R); (7)

S  М_n(C) – матрица внутренних динамических параметров ДДС.

Решение МКУ (6) представляется в аналитической форме записи, имеющей вид корневой пары, структура которой определяется формулой:

S_1,2 = M^1(–C + V  U)/2. (8)

где V, U – соответственно кососимметрическая и симметрическая матрицы.

Следует отметить, что в отличие от обычного (скалярного) квадратного уравнения, которое на комплексной плоскости всегда имеет только два решения, МКУ допускает множество решений, то есть множество корневых пар вида (8):

= M^1(–C + V⁽ⁱ⁾  U⁽ⁱ⁾)/2 (i = 1, 2, … , q).

При этом, если уравнение (6) обладает конечной разрешимостью, то число корневых пар (8) будет равно

q = /2,

где – число сочетаний из 2n по n (n – порядок матрицы S).

Вопрос о конечной (или бесконечной) разрешимости МКУ (6) связан с анализом вспомогательной спектральной задачи, которая здесь не рассматривается. В дальнейшем будем полагать, что МКУ всегда имеет конечную разрешимость. Тогда для его решения достаточно получить какую-либо одну корневую пару при фиксированном индексе «i». Все остальные корневые пары уже не будут содержать новой информации с точки зрения построения общего решения дифференциального уравнения (3). Поэтому верхний индекс «i» у матриц S, V и U может быть опущен. В то же время, наличие нижнего индекса у матричных корней:

S_k, S_l (k, l = 1, 2; k  l)

всегда свидетельствует о принадлежности данных корней только общей корневой паре (8).

Формулу (8) иногда удобнее записывать с помощью одного нижнего индекса (вместо двух) у матрицы S:

S_k = M^1[–C + V – (–1)^kU]/2 (k = 1, 2) (8 а)

либо в более предпочтительном варианте, сохраняя общность выражения:

S_k = M^1(–C + V + U)/2,

S_l = M^1(–C + V – U)/2 (k, l = 1, 2; k  l). (8 б)

По аналогии с решением «как обычно» (для скалярного квадратного уравнения) матричную структуру (8) также можно представить через дискриминант:

M^–1U = , D = 4K – (C + V )M^–1(C – V ). (9)

Здесь D – матрица-дискриминант МКУ

Помимо формул (9) матрицы V и U связаны между собой соотношением

(C + V)M^1U = UM^1(C – V), (10)

представляющим собой условие симметрии матрицы (C + V)M^1U.

Зависимости (9), (10) образуют систему разрешающих уравнений относительно искомых матриц V, U, значения которых определяются с помощью построенной итерационной схемы.

Свойства матричных соотношений

Матричные корни S_k (k = 1, 2) в (8) обладают многими специальными свойствами. Эти свойства проявляются в виде связей между матричными коэффициентами M, C, K, корнями S_k, а также матрицами V и U, входящими в структуру решения (8). Ниже приведены некоторые из этих свойств, играющие ключевую роль в последующем анализе.

Важнейшим свойством матричных корней, принадлежащих общей корневой паре (8), являются соотношения

S_k^ТM + MS_l = –C, S_k^ТMS_l = K (k, l = 1, 2; k  l). (11)

Нетрудно видеть, что эти зависимости имеют сходный вид с формулами Виета для обычного квадратного уравнения. Действительно, первое соотношение в (11) в обобщенном смысле можно трактовать как сумму матричных корней, равную матрице C в линейном члене МКУ, взятой с обратным знаком. Второе соотношение также в обобщенном смысле представляет собой произведение матричных корней, равное свободному члену МКУ (матрице K).

Наряду с формулами (11) имеют место более простые матричные соотношения, устанавливающие свойства симметрии для характеристик решения в структуре корневой пары:

V = –V^T, U = U^T; (12)

US = S^ТU. (13)

Согласно (12), матрица V является кососимметрической, а U симметрической. Свойства (12) находятся в тесной взаимосвязи с соотношениями (11). Действительно, если в левые части формул (11) подставить выражения матричных корней (8 а), то первое соотношение в (11) приводит к тождеству, а второе соотношение дает две независимые формулы: одна из них связана с матрицей-дискриминантом (9), другая представляет равенство (10). Уравнение (13) выражает условие симметрии матрицы US. Это условие выполняется независимо от номера корня в (8), поэтому нижний индекс у S опущен.

Соотношения (11)–(13) по своей природе являются фундаментальными, так как они получены чисто алгебраическим путем из анализа МКУ не прибегая к использованию каких-либо физических свойств системы кроме тех, которые оговорены требованиями (7). Следовательно, формулы (11)–(13) справедливы для любой колебательной системы с проявлением таких физико-механических свойств материала как упругость, вязкость, пластичность, ползучесть и др. и выполняются тогда и только тогда, когда параметры расчетной динамической модели удовлетворяют условиям симметрии и вещественности матриц M, C, K в МКУ (6).