На рис. 3 показаны эпюры изгибающих моментов и прогибов для сталежелезобетонной балки при действии распределенной нагрузки q и сосредоточенной силы F , приложен-
ной в середине второго пролета (
x
94,5
м).
Рис. 3. Эпюры изгибающих моментов и прогибов
Изгибающие моменты в опорных сечениях отрицательные, следовательно, в верхних волокнах поперечных сечений возникают растягивающие напряжения.
На рис. 4 представлены графики изменения изгибающих моментов в 1-м и 2-м опорных сечениях в зависимости от положения сосредоточенной силы.
Рис. 4. Графики изменения изгибающего момента в 1-м и 2-м опорных сечениях
Отрицательный изгибающий момент вызывает растягивающие напряжения в верхних волокнах опорных поперечных сечений балки, которые могут достигать предела прочности бетона на растяжение. В качестве предельных значений изгибающих моментов примем следующие: пред к м; пред к м. Этим значениям изгибающих моментов отвечают значения растягивающих нормальных напряжений:
Значения предельных изгибающих моментов отмечены на рис. 4
65
горизонтальными линиями. Красным цветом на них выделены участки балки, при расположении на которых силы F , в опорных поперечных сечениях могут возникать предельные растягивающие напряжения.
В последующих расчетах мы будем использовать понятие коэффициента снижения
изгибной жесткости в окрестности точки поперечного сечения (0 1) . Коэффициент |
|
|||
учитывает частичное вовлечение бетона в работу на растяжение [1]. Суть коэффициента |
в |
|||
следующем. Если при определенном положении силы F на балке в верхних волокнах попе- |
||||
речного сечения с координатой |
x |
i |
напряжение превышает предельное значение, |
то |
|
|
|
|
|
ослабляем соответствующее сечение, меняя его изгибную жесткость в окрестности точки
EJ (xi ) EJ (xi ) (скорректированная на участке изгибная жесткость при новом поло- |
||||||||
жении силы |
остается неизменной). В расчетах |
значение коэффициента |
принималось рав- |
|||||
ным |
0,7 |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
2. Расчет неразрезной балки переменного сечения методом конечных разностей |
|||||||
|
Дифференциальные уравнения равновесия балки переменного сечения имеют вид |
|||||||
|
|
|
|
|
, |
|
|
(1) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2) |
|
|
|
|
|
, |
|
|
|
где V (x) , |
M (x) - искомые функции прогиба и |
изгибающего момента, |
EJ (x) - изгибная |
|||||
жесткость, |
q |
- интенсивность распределенной нагрузки, Dirac(x xi ) |
- |
дельта функция |
||||
Дирака, фиксирующая положение сосредоточенной силы F в точке |
xi . |
|
|
|||||
|
Для шарнирно опертой по концам балки прогибы и изгибающие моменты в концевых |
|||||||
точках равны нулю: |
|
|
|
|
|
|||
|
|
|
V (0) 0, M (0) 0,V (L) 0, M (L) 0, |
|
|
(3) |
||
где L |
- длина балки. |
|
|
|
|
|
||
|
Решение дифференциальной задачи (1)-(3) выполняется методом конечных разностей. |
|||||||
Отрезок интегрирования |
разбивается на |
участков длиной h L / n . Координаты уз- |
||||||
ловых точек полученной сетки связаны формулой |
|
|
|
|||||
|
|
|
|
) |
. |
|
|
(4) |
|
Значение вторых производных в точке заменяем приближенно с помощью формул |
|||||||
численного дифференцирования: |
|
|
|
|
||||
|
|
|
|
|
|
|
|
(5) |
|
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
|
|
(6) |
|
|
|
|
|
, (i 1,...,n 1) . |
|||
Соотношения (5), (6) справедливы с точность до величин порядка |
. При увеличении числа |
|||||||
дискретных элементов (участков разбиения) погрешность стремится к нулю |
). |
|||||||
|
Краевые условия (3) принимают вид |
|
|
|
|
|||
|
|
|
V (x0 ) 0, M (x0 ) 0,V (xn ) 0, M (xn ) 0 . |
|
|
(7) |
||
66
После подстановки выражений (6) в (5), а затем в (1) приходим к системе линейных алгебраических уравнений относительно неизвестных прогибов V (xi ) в узлах сетки:
(8)
.
На рис. 5 представлена структура матрицы жесткости балки переменного сечения при разбиении балки на n=9 элементов. Матрица симметричная, пятидиагональная, положительно определенная. На рис. 5 показан также фрагмент программного кода ее создания в системе символьной математики Maplе.
Выполняется генерация трех списков: d1 - элементы главной диагонали; d2 - элементы 2-й диагонали (расположенной под главной); d3 - элементы 3-й диагонали. Сформированная таким образом матрица зависит от списка изгибных жесткостей балки EJ[k] в узлах сетки. При изменении значений изгибных жесткостей EJ[k] (на предмет появления трещин в
бетоне) автоматически корректируется матрица A0.
Рис. 5. Структура матрицы жесткости
Вектор узловых сил корректируется при каждом новом положении силы F. Сила F приводится в узлы дискретного элемента, на котором она находится по правилу рычага. Соб-
ственный вес моделируем распределенной нагрузкой |
q . |
Матрица жесткости неразрезной балки получается из матрицы жесткости однопролетной балки ( L 189 м). Сетка вдоль балки наводится так, чтобы координаты промежуточных опор совпадали с узлами сетки. Далее обнуляются строки и столбы матрицы A0 с номерами узлов, отвечающих промежуточным опорам. Элементам главной диагонали матрицы с этими номерами присваиваем значения, равные единице, а соответствующие компоненты вектора правой части делаем равными нулю.
В последующих расчетах шаг дискретизации назначался по результатам сравнения
конечно-разностного |
и точного (построенного в Maple) решений. При шаге |
h L / n 189 / 189 1 м |
погрешность M в определении изгибающих моментов не превы- |
шала
0, 01
%
.
Систему линейных алгебраических уравнений решаем методом Холецкого. Определяем прогибы, после чего по значениям функции прогибов в узлах сетки вычисляем значения изгибающих моментов по формулам (6).
67
3. Алгоритм расчета неразрезной сталежелезобетонной балки на подвижную нагрузку с учетом раскрытия трещин
Расчет неразрезной сталежелезобетонной балки на подвижную нагрузку с учетом раскрытия трещин выполнялся по следующему алгоритму:
1. |
Задаем интенсивность постоянной нагрузки q=5.772 к /м и величину перемещаемой силы |
|||
|
F=400 |
к ; предельные изгибающие моменты |
|
; |
|
|
и коэффициент снижения изгибной жесткости |
0,7 |
; шаг дискретизации |
|
|
|
||
|
балки |
h 1 м ; |
|
|
2. |
Перемещаем нагрузку F по балке с шагом h и вычисляем перемещения и изгибающие |
|||
|
моменты в узлах сетки; |
|
|
|
3.При каждом положении нагрузки сравниваем вычисленные изгибающие моменты с предельным значением изгибающего момента. Если изгибающий момент превышает предельное значение, то ослабляем соответствующее сечение, меняя его изгибную жесткость
(
EJ |
k |
|
|
|
EJk
) (скорректированная на участке изгибная жесткость при новом положе-
нии силы остается неизменной);
(
пред
4.Определяем координаты положений сил, при которых происходит выключение бетона из работы, и координаты ослабленных сечений, отвечающих данному положению сил, строим графики.
На рис. 6 представлена схема раскрытия трещин при перемещении силы F по балке
.1 3500 к м ). Первый отрезок на рисунке отмечен сноской 12 – 64. Первое число
указывает на координату ослабленного сечения.
x 12 м
приложения силы
F
, а второе на координату
x 64 м
Рис. 6. Схема раскрытия трещин по длине балки при пред |
кН м |
Ниже приведен список координат приложения силы и соответствующие координаты ослабленных сечений:
[[12, 64], [15, 65], [17, 63], [18, 66], [22, 67], [25, 2],[27, 68], 34, 69], 81,127], [86, 128], [88, 126], [92,129], [137, 125], [140, 124], [143, 123], [150, 122]].
68
F
На рис. 7 показан график изменения изгибной жесткости при перемещении нагрузки по балке.
Рис. 7. График изменения изгибной жесткости при пред к м
Раскрытие трещин в окрестности первой промежуточной опоры происходило на участках балки с координатами [64, 65. 63, 66, 67, 62, 68, 69] (координата опорного сечения х 63 м ). Длина зоны раскрытия трещин составила 7 м. Последовательность координат сече-
ний появления трещин не является возрастающей. Раскрытие трещин происходит при расположении силы в первом пролете слева от опорного сечения.
Раскрытие трещин в окрестности второй промежуточной опоры (координата опорного сечения х 126 м ) происходило на двух участках балки с координатами: [127, 128, 126,
129] (сила во втором пролете); [125, 124, 123, 122] (сила в третьем пролете). Длина зоны раскрытия трещин составила 7 м. Последовательность координат сечений появления трещин не является возрастающей. Раскрытие трещин происходит при расположении силы слева и
справа от опорного сечения. |
|
Аналогичный расчет выполняем для предельного момента |
. |
Схема раскрытия трещин представлена на рис. 8. |
|
Рис. 8. Схема раскрытия трещин по длине балки при пред |
кН м |
69