Гибкий метод

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

Гибкость членства

Гибкость обратна жесткости . Например, рассмотрим пружину, у которой Q и q являются соответственно ее силой и деформацией:

• Соотношение жесткости пружины равно Q = kq, где k — жесткость пружины.
• Его соотношение гибкости: q = f Q , где f — гибкость пружины.
• Следовательно, f = 1/ k .

Типичное отношение гибкости члена имеет следующую общую форму:

$\mathbf {q} ^{m}=\mathbf {f} ^{m}\mathbf {Q} ^{m}+\mathbf {q} ^{om}$ ( 1 )

где

m = номер участника m .

$\mathbf {q} ^{m}$ = вектор характеристических деформаций элемента.

$\mathbf {f} ^{m}$ = матрица гибкости элемента, которая характеризует склонность элемента к деформации под действием сил.

$\mathbf {Q} ^{m}$ = вектор независимых характеристических сил стержня, которые являются неизвестными внутренними силами. Эти независимые силы порождают все силы на концах членов за счет равновесия членов.

$\mathbf {q} ^{om}$ = вектор характерных деформаций элемента, вызванных внешними воздействиями (такими как известные силы и изменения температуры), приложенными к изолированному, отсоединенному элементу (т. е. с $\mathbf {Q} ^{m}=0$).

Для системы, состоящей из множества членов, соединенных между собой в точках, называемых узлами, отношения гибкости членов можно объединить в одно матричное уравнение, опустив верхний индекс m:

$\mathbf {q} _{M\times 1}=\mathbf {f} _{M\times M}\mathbf {Q} _{M\times 1}+\mathbf {q} _{M\times 1}^{o}$ ( 2 )

где M — общее количество характерных деформаций или сил элементов в системе.

В отличие от матричного метода жесткости , где отношения жесткости элементов можно легко интегрировать с помощью условий узлового равновесия и совместимости, нынешняя гибкая форма уравнения ( 2 ) представляет серьезные трудности. С членами сил $\mathbf {Q} _{M\times 1}$ поскольку основные неизвестные, количество узловых уравнений равновесия, как правило, недостаточно для решения - если только система не является статически определенной .

Уравнения узлового равновесия

Чтобы решить эту трудность, сначала мы воспользуемся уравнениями узлового равновесия, чтобы уменьшить количество независимых неизвестных сил-членов. Узловое уравнение равновесия системы имеет вид:

$\mathbf {R} _{N\times 1}=\mathbf {b} _{N\times M}\mathbf {Q} _{M\times 1}+\mathbf {W} _{N\times 1}$ ( 3 )

где

$\mathbf {R} _{N\times 1}$: Вектор узловых сил на всех N степенях свободы системы.

$\mathbf {b} _{N\times M}$: Результирующая матрица узлового равновесия

$\mathbf {W} _{N\times 1}$: Вектор сил, возникающих от нагрузки на элементы.

В случае детерминированных систем матрица b квадратная и решение для Q можно найти сразу из ( 3 ) при условии, что система устойчива.

Основная система

Для статически неопределимых систем M > N и, следовательно, мы можем дополнить ( 3 ) уравнениями I = M − N вида:

$X_{i}=\alpha Q_{j}+\beta Q_{k}+\cdots \qquad i=1,2,\ldots ,I$ ( 4 )

Вектор X — это так называемый вектор избыточных сил, а I — степень статической неопределенности системы. Обычно мы выбираем j , k ,…, $\alpha$, и $\beta$ такой, что $X_{i}$ — это реакция поддержки или внутренняя сила, действующая на члена. При соответствующем выборе избыточных сил теперь можно решить систему уравнений ( 3 ), дополненную ( 4 ), и получить:

$\mathbf {Q} _{M\times 1}=\mathbf {B} _{R}\mathbf {R} _{N\times 1}+\mathbf {B} _{X}\mathbf {X} _{I\times 1}+\mathbf {Q} _{v\cdot M\times 1}$ ( 5 )

Замена в ( 2 ) дает:

$\mathbf {q} _{M\times 1}=\mathbf {f} _{M\times M}{\Big (}\mathbf {B} _{R}\mathbf {R} _{N\times 1}+\mathbf {B} _{X}\mathbf {X} _{I\times 1}+\mathbf {Q} _{v\cdot M\times 1}{\Big )}+\mathbf {q} _{M\times 1}^{o}$ ( 6 )

Уравнения ( 5 ) и ( 6 ) представляют собой решение для первичной системы , которая является исходной системой, которая была статически детерминирована путем разрезов, обнажающих избыточные силы. $\mathbf {X}$. Уравнение ( 5 ) эффективно сводит набор неизвестных сил к $\mathbf {X}$.

Уравнение совместимости и решение

Далее нам нужно настроить $I$ уравнения совместимости, чтобы найти $\mathbf {X}$. Уравнения совместимости восстанавливают требуемую сплошность на сечениях путем задания относительных перемещений $\mathbf {r} _{X}$ при избыточности X к нулю. То есть, используя метод единичной фиктивной силы :

$\mathbf {r} _{X}=\mathbf {B} _{X}^{T}\mathbf {q} =\mathbf {B} _{X}^{T}{\Big [}\mathbf {f} {\Big (}\mathbf {B} _{R}\mathbf {R} +\mathbf {B} _{X}\mathbf {X} +\mathbf {Q} _{v}{\Big )}+\mathbf {q} ^{o}{\Big ]}=0$ ( 7а )

или $\mathbf {r} _{X}=\mathbf {F} _{XX}\mathbf {X} +\mathbf {r} _{X}^{o}=0$ ( 7б )

где

$\mathbf {F} _{XX}=\mathbf {B} _{X}^{T}\mathbf {f} \mathbf {B} _{X}$

$\mathbf {r} _{X}^{o}=\mathbf {B} _{X}^{T}{\Big [}\mathbf {f} {\Big (}\mathbf {B} _{R}\mathbf {R} +\mathbf {Q} _{v}{\Big )}+\mathbf {q} ^{o}{\Big ]}$

Уравнение ( 7b ) может быть решено относительно X , а затем силы на стержнях находятся из ( 5 ), тогда как узловые перемещения могут быть найдены по формуле

$\mathbf {r} _{R}=\mathbf {B} _{R}^{T}\mathbf {q} =\mathbf {F} _{RR}\mathbf {R} +\mathbf {r} _{R}^{o}$

где

$\mathbf {F} _{RR}=\mathbf {B} _{R}^{T}\mathbf {f} \mathbf {B} _{R}$ – матрица гибкости системы .

$\mathbf {r} _{R}^{o}=\mathbf {B} _{R}^{T}{\Big [}\mathbf {f} {\Big (}\mathbf {B} _{X}\mathbf {X} +\mathbf {Q} _{v}{\Big )}+\mathbf {q} ^{o}{\Big ]}$

Перемещения опор, происходящие в резервных местах, можно включить в правую часть уравнения ( 7 ), а перемещения опор в остальных местах необходимо включить в правую часть уравнения (7). $\mathbf {r} _{X}^{o}$ и $\mathbf {r} _{R}^{o}$ также.

Преимущества и недостатки

Хотя выбор избыточных сил в ( 4 ) кажется произвольным и затруднительным для автоматического расчета, это возражение можно преодолеть, перейдя непосредственно от ( 3 ) к ( 5 ), используя модифицированный процесс исключения Гаусса–Жордана . Это надежная процедура, которая автоматически выбирает хороший набор резервных сил для обеспечения численной стабильности.

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

Верхние аргументы были действительны до конца 1990-х годов. Однако недавние достижения в области численных вычислений показали возвращение метода сил, особенно в случае нелинейных систем. Были разработаны новые основы, которые позволяют «точные» формулировки независимо от типа или природы нелинейностей системы. Основным преимуществом метода гибкости является то, что ошибка результата не зависит от дискретности модели и что это действительно очень быстрый метод. Например, упруго-пластическое решение непрерывной балки с использованием силового метода требует всего 4 элемента балки, тогда как коммерческий код FEM, основанный на жесткости, требует 500 элементов, чтобы дать результаты с той же точностью. В заключение можно сказать, что в случае, когда решение задачи требует рекурсивных оценок силового поля, как в случае структурной оптимизации или идентификации системы , эффективность метода гибкости неоспорима.

См. также

• Метод конечных элементов в строительной механике
• Структурный анализ
• Метод жесткости

Внешние ссылки

• Последовательные деформации - силовой метод. Архивировано 10 августа 2008 г. в Wayback Machine.