chm_3

Плоско-деформированное состояние

В плоско-деформированном состоянии находятся тела, форма и размеры поперечного сечения, условия нагружения которых не зависят от одного из направлений. Размер тела в этом направлении велик, и продольной деформацией можно пренебречь. В качестве примера рассматривается длинный брус (рис. 4 .1, а), находящийся на твердой горизонтальной площадке под действием вертикальной нагрузки, не изменяющейся вдоль оси z. Форма и размеры его поперечного сечения вдоль этой же оси не изменяются.

Г Г_p



z 0 x

а б

Рис. 4.1. Схема плоско-деформированного состояния (а) и форма

поперечного сечения тела (б)

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

, ,

, , , .

С другой стороны, из того же закона Гука следует

И, в силу ,

то есть компонента тензора напряжений не является независимой величиной. В этом случае матричное соотношение ( 4 .15) представляется в виде

, (4.19)

Учитывая, что , и решение задачи не зависит от переменной z, в системе разрешающих соотношений ( 4 .11) остаются лишь два уравнения

В матричной форме эти уравнения записываются в виде, аналогичном выражению ( 4 .13). Вид матрицы {_m} определен выражением ( 4 .19), а также используются обозначения

Связь компонент тензора деформаций и вектора перемещений определяется в виде, аналогичном выражению ( 4 .17), где

С учетом введенных обозначений окончательно система алгебраических уравнений относительно искомых коэффициентов разложения имеет вид, аналогичный выражению ( 4 .18). Необходимо отметить, что при интегрировании по области  и границе Г,

причем во всех слагаемых системы уравнений ( 4 .18) появляется общий множитель L, равный размеру рассматриваемого тела в направлении оси z, который можно сократить. Иными словами, интегрирование в дальнейшем производится только по области _p поперечного сечения и по контуру _p границы (рис. 4 .1, б).

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

где i, j, k – номера вершин произвольного конечного элемента, пробные функции _i, _j, _k определены выражениями (2.5). В матричной форме эти выражения представим в виде

Соответственно, матрицы [B_r], r = i, j, k имеют вид

, , ,

система уравнений ( 4 .18) преобразуется к форме

. (4.20)

Вычисляются значения матриц, входящих в эту систему уравнений:

;

При выводе этих выражений принято, что в пределах конечного элемента модуль упругости Е, температура T, коэффициенты Пуассона  и температурного расширения материала  постоянны; S_p – площадь p-го треугольного конечного элемента. Вводятся матричные обозначения

где {u_p} – вектор всех узловых перемещений p-го треугольного конечного элемента, [K_p] – матрица жесткости для этого же конечного элемента. Теперь система уравнений ( 4 .20) для одного конечного элемента представляется в форме (суммирование по индексу p не производится)

Пример 4.1. Рассматривается осадка длинной стальной полосы с квадратным поперечным сечением размером 22 м², зажатой между двумя гладкими горизонтальными плитами (рис. 4 .2). Каждая из плит осаживается на величину . Требуется определить, на какую величину сместятся боковые стороны этой полосы в результате деформирования.

Так как форма поперечного сечения является симметричной, можно рассматривать лишь четверть исследуемой области (рис. 4 .2, б). Кинематические граничные условия указаны на том же рисунке. Поскольку плиты, деформирующие полосу, абсолютно гладкие, трение между ними и стальной полосой отсутствует, касательные усилия на контактной поверхности равны нулю. На свободной боковой поверхности полосы отсутствуют как нормальные, так и касательные нагрузки.

Для упрощения анализа алгоритма определения напряженно-деформированного состояния объекта рассматриваемая часть поперечного сечения полосы аппроксимируется только двумя конечными элементами, как это показано на рис. 4 .2, б.

y y 



0 x

а б

Рис. 4.2. Геометрическая схема осадки длинной полосы (а) и схема кинематических граничных условий (б)

Пробные функции для первого треугольного элемента, показанного на рис. 4 .3, определены в виде

Содержание