ПОРТЫ, ПОРТОВОЕ хозяйство И ТРАНСПОРТНАЯ ЛОГИСТИКА
DOI: 10.24143/2073-1574-2018-2-113-132 УДК 629.08.
Н. Н. Панасенко, А. В. Синельщиков, П. В. Яковлев
РАСЧЁТНЫЙ АНАЛИЗ СЕЙСМИЧЕСКОЙ БЕЗОПАСНОСТИ ГРУЗОПОДЪЁМНЫХ КРАНОВ
Выбор методов теории сейсмостойкости применительно к оценке сейсмостойкой безопасности грузоподъёмных кранов, эксплуатируемых при строительстве зданий и сооружений, зависит от наличия и вида модели сейсмического воздействия. При наличии сейсмического входа в виде сейсмических коэффициентов динамичности либо их модификации в виде сейсмических спектров ответа применяют линейно-спектральный метод, его «плоскую» версию по СП 14.13330.2014 «Строительство в сейсмических районах», а также в НП-031-01 «Нормы проектирования сейсмостойких атомных станций». Метод динамического анализа теории сейсмостойкости сооружений является развитием линейно-спектрального метода. С его помощью процедура конечно-элементного моделирования действительных конструкций упрощается, но метод динамического анализа требует использования методов численного компьютерного интегрирования больших систем дифференциальных сейсмических уравнений 2-го порядка.
Ключевые слова: сейсмическое воздействие, трёхкомпонентная акселерограмма, синтезированная (вероятностно-статистическая) акселерограмма, сейсмический спектр ответа, линейно-спектральный метод, метод динамического анализа, математическая модель сейсмического воздействия.
Введение
Землетрясения представляют собой одну из самых разрушительных природных катастроф с последствиями в виде человеческих жертв и огромных материальных убытков.
По статистическим данным, в мире ежегодно происходит более 300 тыс. землетрясений. В СНГ сейсмически активной считается 28 % территории, ежегодно регистрируется до 75 землетрясений, причём в среднем раз в три года происходит одно разрушительное землетрясение.
Задача сокращения материальных убытков и уменьшения человеческих жертв при возможных землетрясениях в ходе освоения сейсмически активных территорий обусловлена требованиями по сейсмической безопасности грузоподъёмных кранов, транспортирующих, в том числе, опасные радиоактивные грузы в портах и терминалах. Основой сейсмической безопасности являются методы расчёта подъёмных сооружений на сейсмические воздействия по оценке напряженно-деформированного состояния их несущих конструкций, связанных с обеспечением их прочности для восприятия сейсмических нагрузок. В статье показано, что требования по разработке современных проектов подъёмных сооружений удовлетворяются, если учитывается пространственный характер воздействия, а сооружения рассматриваются как единые пространственные системы, кроме того, учитываются различного характера деформации и геометрические и кинематические нелинейности.
Установлено, что расчёт и проектирование сейсмостойких подъемных сооружений стали возможными за счёт развития компьютерной техники и компьютерного математического обеспечения, основанного на развитии аппарата общей и строительной механики и численных методов решения больших систем дифференциальных уравнений сейсмических колебаний второго порядка.
Математическая модель сейсмического воздействия
Уравнение вынужденных сейсмических колебаний недемпфированной системы расчётно-динамической модели (РДМ) подъёмного сооружения с п степенями свободы имеет вид
[M ]{V(t)}+ [K ]{V (t )}= P„ + p - [M ]{äCt)} C1)
где [М] и [Ä] - матрицы масс и жёсткости РДМ конструкции; {A(t)} - вероятностно-
статистическая акселерограмма (ВСА) землетрясения; {Рст}, {Рд} - векторы статических и динамических эксплуатационных и технологических нагрузок на грузоподъёмную машину соответственно; {V(t)} и V(t)} - вектор перемещений РДМ и его вторая производная.
В соответствии с п. 5.1.1. Руководства по безопасности 006-98 «Определение исходных сейсмических колебаний грунта для проектных основ» [1] для синтезирования расчётных акселерограмм используется полученный в результате обработки набор аналоговых акселерограмм, подходящих для сходных сейсмических и грунтовых условий, представленных в работах [2-6], сводный ансамбль которых приведён на рис. 1.
¡эд «
о
R
я я
а
о и
о ^
0,20 0,15 0,10 0,05 0,00 -0.05 -0,10 -0,15 -0,20
0,0
5,0
1-,0
15,0
20,0
Т, с
о
«
и
Щ
о и
о
0,20 0,15 0,10 0,05 0,00 -0.05 -0,10 -0,15 -0,20
0,0 5,0 10,0 15,0 20,C
Т, с
0,25 tg 0,20 0,15 0,10 0,05 0,00 -0,05 -0,10 > -0,15 -0,20
ч о
Sí
я я
Л
о и
0,0
5,0
fc
... ,мш i.l Ii
рт'Г'Т
likhlil : 'ТТГилШлШ
r гг
'II r I
1
10,0
15,0
20,0
Т, с
Рис. 1. Сводный ансамбль реальных 7-балльных акселерограмм максимального расчётного землетрясения MSK-64 [2-6] на дневной поверхности
Принято считать, что расчётная акселерограмма РБ 006-98 [1]
m
aß1 = A( t Bi sin (rot + ф)
(2)
1=1
является консервативной. Здесь А(0 - огибающая, а фазовые углы фг- представляют собой равномерно распределённые в интервале от 0 до 2п случайные величины. Шаг по частоте Дю определяется из условия гладкости спектра реакции как (- ю;-1)/ю;-1 = 0,06347 .
Расчёты производятся в интервале периодов Т от 0,05 до 3 с. В качестве первого приближения для амплитуды В^ в (2) используются значения, непосредственно взятые с заданной кривой Р(Т) [7] для соответствующих значений частот юг-.
Уход от консервативной акселерограммы (2) авторам настоящей работы представляется в разработке вероятностно-статистическая акселерограммы, получаемой на основе обработки ансамбля исходных, предпочтительно реальных, акселерограмм землетрясений баллов MSK-64 [2-6].
Для построения ВСА нами приняты следующие допущения [8, 9]:
- на времени действия тэ эффективной фазы землетрясения 4 < тэ < 10 с сейсмическое воздействие (СВ) является случайным стационарным процессом с нулевым матожиданием, дисперсией о2, корреляционной функцией (КФ) - К(т) и функцией спектральной плотности (ФСП) - Ga(ю);
- сейсмическое воздействие, заданное акселерограммой а(0, имеет нормальный закон распределения;
- сейсмическое воздействие, в отличие от СНиП 11-7-81* [7], задаётся для трёх направлений пространства: двух горизонтальных X, У и вертикального Z. Причем вероятность проектного землетрясения равна Тс / 103, а максимального расчетного землетрясения (МРЗ) -Тс / 104, где Тс - нормативный срок службы проектируемого сооружения, например, согласно РД 22-01-97 [10]. Сейсмическое воздействие как случайную функцию а (/) представим ансамблем [2-6] выборочных функций {(/)}, каждая из которых не описывает всех свойств СВ.
Проведя статистическую обработку ансамбля исходных акселерограмм, получим воздействие, учитывающее все свойства ансамбля. Для этого каждая акселерограмма исходного ансамбля переоцифровывается с одним и тем же шагом Дt (0,01 < Дt < 0,03 с), и для каждой из них устанавливается длительность, соответствующая длительности тэ эффективной фазы ансамбля. В результате получается ансамбль переоцифрованных с одинаковым по времени шагом акселерограмм равной длительности, рассматриваемых как реализация случайного процесса Л(/).
Для каждого момента времени tk (с шагом оцифровки) проводится осреднение мгновенных значений процесса землетрясения: определяется (а|/к) и а2( а|<к) - математическое ожидание
и среднеквадратическое значение (дисперсия) соответственно:
1=X
< а\ь > = £( а^к)/X; (3)
1=X 2 ,
о 2% >=Х( а,\Ч -< а\Ч >)/( х -1), (4)
1=1
где X - число реализаций сейсмических процессов. Выборка а объемом X в (3), (4) имеет нормальный закон распределения. С изменением объёма выборки будет меняться дисперсия (4). В этом случае для момента времени 4 значение (а|/к) процесса является случайной величиной со всеми параметрами дисперсии
о< ,. > = о ,. л/Х (5)
< а\*к > а\*к /
и матожидания выборочного среднего, которое с вероятностью Р = F(ир) = 2Ф(ир) будет находиться в интервале
К а К) - иРо ъ/Я ) < Ц < (( а^) + иро ъ/Я ), (6)
где верхняя граница (6) является среднестатистической акселерограммой (ССА):
(а)к =(а\ь)+ир°а0 • (7)
В (6) и (7) F(UP), Ф(иР) - табулированные функции нормального распределения (нормированного и Лапласа) соответственно; ир - квантиль нормального распределения, соответствующий принятой проектировщиком доверительной вероятности Р, значения которого приведены в табл. 1.
1=1
Таблица 1
Квантили UP нормального распределения
^^^^^^ Показатель Критерий ^^^^^^ Доверительная вероятность Р
0,90 0,95 0,975 0,990 0,995 0,999
Односторонний 1,282 1,645 1,96 2,362 2,576 3,0
Двухстор онний 1,645 1,96 2,362 2,576 3,0 3,3
Заметим, что выбор доверительной вероятности Р в (6) не является математической задачей, а определяется конкретно решаемой проблемой. Величина
а = (1 - Р) (8)
называется уровнем значимости или риском пользователя (риском 1-го рода) и определяет критическую область, в которую попадает а процентов неучтённых возможных средних значений выборки X.
При малом числе реализаций процесса вместо квантиля ир в (6) и (7) следует принимать квантиль tp распределения Стьюдента (табл. 2.).
Таблица 2
Квантили ^ распределения Стьюдента
^^^^^ Показатель Объём выборки^^^^^ Доверительная вероятность Р
0,95 0,99 0,999
5 3,04 5,04 9,43
15 2,22 3,08 4,28
20 2,15 2,93 3,98
30 2,08 2,80 3,72
40 2,03 2,71 3,53
50 1,99 2,64 3,41
100 1,96 2,58 2,29
Анализ реальных акселерограмм [2-6] показывает, что у консервативных акселерограмм (2), построенных для различных ансамблей, пиковые ускорения в момент времени tk могут существенно отличаться, поэтому для повышения точности представления СВ и повышения качества сейсмических расчётов применяют ВСА. Учитывая, что среднеквадратичное значение о ^ в (6)
тоже является случайной величиной, имеющей доверительный интервал
(1 - ар ) о < о I <(1 + ар ) о I , (9)
V ^ р / а\1к а\1к \ 1 р /
где др - квантиль распределения (табл. 3), поэтому для построения вероятностно-статистической акселерограммы целесообразно принимать верхнюю границу Р-процентного интервала среднеквадратичного значения оиз (9).
Таблица 3
Значения квантилей qp
Показатель дискретных^^^ значений ансамблй--^^ акселерограмм Доверительная вероятность Р
0,95 0,99 0,999
5 1,37 2,67 5,64
10 0,65 1,08 1,80
15 0,46 0,73 1,15
20 0,37 0,58 0,88
30 0,28 0,43 0,63
40 0,24 0,35 0,50
50 0,21 0,30 0,43
100 0,14 0,20 0,27
После этого вероятностно-статистическая акселерограмма с учётом (9) представляется в виде
ak =
[(ak\h) + UP G(a\h)
+ U„
G
(1+qP)
(10)
она учитывает Р процентов свойств всей исходной информации ансамбля {аОД} [2-6] и может быть использована как модель сейсмического воздействия на грунте при расчёте зданий и сооружений на сейсмостойкость линейно-спектральным методом и методом динамического анализа. Очевидно, что ВСА (10) интенсивностью 7 баллов по шкале MSK-64 рекомендуется использовать для практических расчётов зданий и сооружений на СВ по МДА (рис. 2), причём выбор значений доверительной вероятности Р в (10) производится в зависимости от вероятности (8).
ч о
R
я я
Л
о и
о >
0,15
0,10
0,05
0,00
-0.05
-0,10
20,0
Т, с
о
«
и
5
Щ
щ
о и
о
Т, с
ч о
Sí
я я
Л
о и
о >
20,0
Т, с
k
а
б
в
Рис. 2. Расчётная ВСА максимального расчётного землетрясения 7 баллов MSK-64 на отметке +0 на дневной поверхности: а, б - горизонтальные компоненты х, y; в - вертикальная компонента z
В целях практических расчетов сооружений на СВ получение значений оцифровки любой акселерограммы a(t) с любым шагом в пределах длины полупериодов осуществляется с использованием формулы
a (t + т) = am sin (лг/L ), (11)
в которой L - длительность импульса; t - время начала импульса; т - время в пределах импульса, 0 < т < Ь ; ат - амплитуда импульса. Здесь следует учесть, что ВСА (10) обеспечены оцифровками как процесс на временной оси с шагом At (0,01 < ^ < 0,08 с). Пользователю трёхкомпонентными ВСА (10) при сейсмическом проектировании зданий и сооружений МДА следует помнить, что их качество при полноте информации большей, чем [2-6], может быть улучшено. Кроме того, их качество должно оцениваться через призму сведений [7], согласно которым максимальные ускорения пиков акселерограмм 7, 8 и 9 баллов отличаются друг от друга в 2 раза и составляют 1, 2 и 4 м/с соответственно.
В заключение следует указать, что ВСА при наличии РДМ рассчитываемого объекта служит основой для построения методом динамического анализа поэтажной ВСА, на уровне установки встроенного технологического оборудования, крановых конструкций и др. с целью их детального расчёта на сейсмостойкость.
С целью оценочного анализа ВСА 7 баллов (см. рис. 2) МРЗ MSK-64 на рис. 3 и 4 представлены её сейсмический спектр ответа (ССО) и сейсмические коэффициенты динамичности (СКД) соответственно на дневной поверхности на временной и частотной оси, в том числе, для сравнения, СКД СНиП И-781*[7], для грунтов II категории.
о
R
я я
Л
о и
о >
0,6 0,5
0,4 0,3 0,2 0,1 0
\ Г 0,02 0,04 ---0,08 ■
«Л
J ?учА Л
/л V \ , \
V Ч/\,
— - S У
10
15
20
25
30
Частота со, Гц
ч
о «
я
я я
Л
о и
о >
0,5 0,4
0,3 0,2
0,1 0
10
15
20
Л 0,02 0,04 0,08
г И1' tyj; * \ - \ у/\ ^
& v' 14 -* *■. V * J -
■ч - - ------ -----
25
30
Частота со, Гц
0,35
tg
0,3
^
(О 0,25
ч.
0,2
0,15
Я
я 0,1
а о 0,05
0
0,02 0,04 0,08
д Л л
J^jj V s _ \ \ лЧ
\ — — ' — " ^ \
10
15
20
25
30
Частота со, Гц
Рис. 3. Сейсмические спектры ответа ВСА максимального расчётного землетрясения 7 баллов М8К-64 на дневной поверхности (затухание 0,02; 0,04; 0,08): а, б - горизонтальные компоненты х, у; в - вертикальная компонента г
а
0
5
б
0
5
в
0
5
X
Р(ш)
6 5 4 3 2
1
0
-СКД СП 14.13330.2014
-0,02
---0,04
--0,08
Частота ю, Гц
У
Р(ю) 5
4
3
2 1
0
Л Л А 1 | СКД СП 14.13330.2014 0,02 0,04 0,08 -
Дг уУЛл Г\ -----
/ у \д /
/ г
//у ^ \ • о N _ Ч__
— "" "" "ч. ■ч - V __
10
15
20
25
30
Частота ю, Гц
г
Р(ю) 6
5
4
3
2
1
0
Л | 1 СКД СП 14.13330.2014 0,02 0,04 0,08
А
<:\
/ . чЛУ \ V ь' __ / . \ . ~ У V
—— л" Ч ____* ^ " N _
10
15
20
25
30
Частота со, Гц
Рис. 4. Сейсмические безразмерные коэффициенты динамичности Р(ю) ВСА максимального расчётного землетрясения 7 баллов MSK-64 на дневной поверхности (затухание 0,02; 0,04; 0,08): а, б - горизонтальные компоненты х, у; в - вертикальная компонента г
Для оценки амплитудных характеристик ВСА 7 баллов (см. рис. 2) на рис. 5 представлена её амплитудно-частотная характеристика (АЧХ).
а
б
0
5
в
0
5
g
0,004 0,0035 0,003 0,0025 0,002 0,0015 0,001 0,0005 0
10
15
20 25 30
Частота га, Гц
С
0,004
0,0035 0,003 0,0025 0,002 0,0015 0,001 0,0005 0
Частота га, Гц
N
0,0025 0,002 0,0015
0,001 0,0005 0
1 1
Will lu, f'| Uli L1 IL »
ff ' flWv (VW fVMVvJ
10
15
20
25
30
Частота га, Гц
Рис. 5. Амплитудно-частотная характеристика ВСА максимального расчётного землетрясения 7 баллов М8К-64: а, б - горизонтальные компоненты х, у; в - вертикальная компонента г
Очевидно, что пользователю ВСА 7 баллов при необходимости следует учитывать грунтовые условия, как согласно [7, 12], так и результатам обеспеченности и сейсмомикрорайониро-вания, например, согласно [3, 13]. Кроме того, оцифровку ВСА 7 баллов МРЗ MSK-64 можно найти по рекомендациям работы [11].
Необходимо указать, что затуханию = 0,02; 0,04; 0,05 и 0,08 (см. рис. 3, 4) соответствует логарифмический декремент затухания колебаний 5з, определяемый из формулы
^=5з (2п)2 +5;
=5/2л,
(12)
величины которого принимают значения 0,12; 0,25; 0,314 и 0,5 соответственно, что отвечает сведениям [14].
Математическая модель линейно-спектрального метода
В качестве расчётного базиса линейно-спектрального метода (ЛСМ) примем собственные векторы, полученные, в частности, для стержневой динамической системы сооружения из уравнений свободных колебаний «-го порядка
а
0
б
в
-0,5
[М )}+[К Г (г )}={0},
(13)
полученное из (1). Для системы со многими степенями свободы исключим время t в (13) с помощью вектора перемещений
(V (г)} = [Ф]эт(ю? + 8)
(14)
и его второй производной по времени
д \
(—} = - [Л] [Ф] 8Ь(ю г + 8) = - [Л]М,
(15)
что приводит к системе уравнений для собственных значений относительно форм колебаний
[М]-1[К][Ф] - [Ф][Л] = (0},
в которой [Ф] - фундаментальная матрица произведения матриц [М] К] порядка п х п, составленная из собственных векторов, расположенных по столбцам
[Ф] = [{Ф}-{Ф- }-{ф- }] =
чФп1 у
Ч. >
чФ пт у
'Ф1п ^
Чфпп У
(16)
характеризующая форму колебаний полной системы, у которой изменяются только амплитуды и 8 - фазовый угол; [ Л ] - диагональная матрица собственных значений, элементами которой являются квадраты собственных частот
[Л]=кю2...ют ...ю^ х п.
Подставив выражения (14) и (15) в (13), получим
-[Л] [М][Ф] 8ш(ю г + 8) + [К] [Ф] 8ш(ю г + 8) = (0} , откуда приходим к уравнению для собственных значений (однородной системе)
(-[Л] [М ] + [К ])[Ф] = (0}.
(17)
(18)
Условие нетривиальности решения однородной системы (18) приводит к частотному уравнению метода конечных элементов (МКЭ) в прямой форме
[М]-1 [К] - [Е][Л] = 0,
(19)
легко разрешаемому с помощью ЭВМ, где [Е] - единичная матрица.
Здесь всегда принимается во внимание, что решение частотного уравнения, согласно (18), удовлетворяет условию [М ]-1[К ] = [Ф][Л][Ф]-1 .
Поскольку уравнение движения (1) является линейным, связь между обобщёнными координатами {V} и нормальными координатами {у} определяется соотношением
{V}=[ФМ
(20)
Подставив в (1) выражение (20) и его вторую производную и умножив полученное слева уравнение на транспонированный вектор т-й формы колебаний (фт}т, получим распадающееся уравнение п-го порядка (векторы {Рст}, {Рд} опущены):
{фт }Т [ М ][Ф]М+К }Т [ М ][Ф]М=-{ФИ }Т [ М ]{А (г)},
(21)
которое, согласно условиям ортогональности, доказанным Рэлеем, с учётом (18) распадается на несвязанные уравнения по формам колебаний:
+ К„¥ = -(ф„)т [М]{Л(0}, (22)
в которыхМт и Кт - модальные масса и жёсткость соответственно:
Ма ={фт }Т [М]{фт }; (23)
Кт ={фт }Т [К]{фт }. (24)
Вектор в правой части (22) с учётом (23), (24) представляет собой долю сейсмической нагрузки на сооружение, вызывающую колебания его РДМ по т-й форме. Если СВ в (22) представим в виде
t )}={0S}A (t), (25)
где {cos} - трёхкомпонентный вектор направляющих косинусов трёхкомпонентной ВСА (см. рис. 2) в общей системе координат (ОСК), то уравнение (22) с учётом затухания (12) примет вид
V + 2^ т и т V + и! V = - DmAa (t), (26)
где ют - собственная т-ая частота колебаний осциллятора, определяемая из выражения (17): ®т = Кт/Мт ; Dm - коэффициент мы колебаний РДМ сооружения:
т
®т = Кт/Мт ; Dm - коэффициент влияния с учётом (23, 24) - величина постоянная для т-й фор-
От ={фт }Т [МК/Мт . (27)
Теперь покажем, что матрица жёсткости полой стержневой системы сооружения порядка п х п, состоящая из X конечных элементов определяется методом суперпозиции:
[к L„ = X
jk=i
№ [к] № [к]
14x14 J
ОСК
nxn
(28)
согласно которому матричная форма матрицы отдельного стержневого КЭ [К ^^ должна быть
ОСК
представлена в виде блока, в котором [К]'6 = ([К]6Х6 , где Т - индекс транспонирования. Известно, что матрица жёсткости масс из местной системы координат (МСК) в ОСК переводится с применением диагональной матрицы преобразования [г]'*
jk . 14x14 •
№14 = ([Т Кх14 )Т х([ККхм ]х([Т ]*14 ). (29)
ОСК ч МСК У
в которой
[т ]14х14 = Г[Л] 3х3 ^ \|3х3 ^ \1Эх3 ^ Вх3 I, (30)
а матрица [Л]3х3 в (30) - матрица направляющих косинусов углов Эйлера [8, 9]; в практических целях матрицу [К]4х14 стержневого КЭ обычно представляют в блочном виде типа (28).
4 x14
МСК
Из (22) и (26) стало очевидно, что сейсмические нагрузки по своей природе являются расчётными инерционными нагрузками. Для их определения необходимо знать абсолютные ускорения колебаний земли {Л(г)} из (25). Для многомассовых систем грузоподъёмных кранов и крановых рельсовых путей расчётные сейсмические нагрузки определяются вектором
£ (* )}=[М )}, (31)
в котором {Я^)} - вектор абсолютных ускорений земли (см. рис. 2) в обобщённых координатах
(20). На основании (22) и (31) вклад от т-й формы колебаний конструкции сооружения в искомую теории ЛСМ сейсмостойкости расчётную сейсмическую нагрузку с учётом (27) имеет вид
С )}т = М ]{Фт } Р У„ ,т (' ). (32)
Здесь уа,т (/) - абсолютное ускорение т-го осциллятора, испытывающего СВ {¿(г)}.
Принципиально важным моментом в ЛСМ пространственных конструкций зданий и сооружений с крановыми нагрузками является переход от динамической задачи (32) к квазистатической. Это осуществляется заменой временной функции уат(/) в (32) на постоянную величину уат,
которая представляет собой максимальную реакцию осциллятора с частотой ют на воздействие, заданное акселерограммой А(?) (25). При этом величину у ат (/) следует получать из сейсмического
спектра ответа (рис. 3) для заданной акселерограммы [15, 16] (см. рис. 2), поэтому расчёт т-й сейсмической силы (32), в соответствии с ЛСМ, осуществляется с использованием ССО с учётом (27) по формуле
Мт =[М ]{фт Рт^ К ,4т ) , (33)
где Щшт, 4т) - отсчёт ССО с затуханием Е,т на частоте ют, построенный для акселерограммы ) из (25) [16]. При этом у а т = V(гат,4т). Если исходная акселерограмма (в данном случае ВСА) задана трёхкомпонентным воздействием вида (25), что чаще всего имеет место на практике, сейсмическую нагрузку (33) по т-й форме колебаний на конструкцию сооружения вычисляют по формуле
М. = М ]{фт ^ ^ ^ )}, (34)
{фт } [М]{фт }
в которой {Щют, 4т)} - вектор порядка п*, построенный на основе трёхкомпонентного сейсмического спектра ответа (см. рис. 3) с использованием ВСА (9). При этом п* < п, где п - число степеней свободы РДМ сооружения, «* = 30-50 Гц, где:
{V (гаи ,4т )} =
{{V (гат,4т)},{| (гат)},{| (гат)},0,0,0}^
{{V ( гат ,4 т )},{| ( гат ,4т )},{| ( ®т ,4т )},0Д0£
(35)
Фактически в рамках ЛСМ вводится понятие квазистационарных расчётных сейсмических сил {^}т (34), каждый вектор которых по т-й форме колебаний РДМ сооружения характеризуется тем, что при приложении эквивалентного ей вектора статических сил конструкция сооружения получит сейсмические смещения { в ОСК, определяемые из уравнения [17]:
К Г }т = {*}. . (36)
Используя полученный вектор смещений {V}т из (36), согласно теории метода конечных элементов, в МСК определяется известный вектор внутренних усилий {2}т в каждом КЭд РДМ конструкции крана по т-й форме колебаний:
(Ж =[К]0Хуг ([ТКХ14 {V}т ), (37)
где [гКх14 - матрица (30) преобразования ОСК ^ МСК; [к]14х14 - матрица жёсткости КЭ;к;
МСК МСК
{V }„ - вектор перемещений узлов / и к КЭ;к 12-го порядка, выделенный из вектора {V }т п-го порядка, который получен в результате решения (36):
{v }m=
(8Х ,5, ,5 г, Фх, ф, ,ег ,0Z) (5Х ,5 , ,5 z, Фх, ф, ,0z ,0Z)
где в строках первые три параметра - линейные перемещения узла /(к) по осям Х12 ОСК, далее три угловые и производная от угла закручивания.
Результирующие внутренние усилия от действующих расчётных сейсмических сил (34) должны быть получены суммированием векторов (37) для всех учитываемых форм колебаний (как правило, до 30 Гц). Однако, поскольку их значения для разных собственных форм достигаются в различные моменты времени, они не могут быть определены в рамках ЛСМ. Поэтому суммарные расчётные сейсмические силы определяют с помощью эмпирических формул, установленных путём сопоставления расчёта по ЛСМ и методом динамического анализа (МДА) на действие акселерограммы, в частности, среднеквадратичным суммированием:
Q}s =JX(QL )2, (38)
где п* << п - число учитываемых форм колебаний. Для близко расположенных частот, таких как Ю < 1,1 Юм внутренние усилия (37), найденные методом среднеквадратичного суммирования (38), оказываются заниженными по сравнению с внутренними усилиями, найденными с помощью МДА. В этом случае авторы [18] рекомендуют алгебраическое суммирование, а для полученных сумм - среднеквадратичное суммирование:
1
ч q
X XIQL
у
(39)
где q - число групп; гч - число частот в q-й группе. Полученные среднеквадратичным суммированием (39) внутренние усилия не имеют знака и, следовательно, отсутствует информация о направлении их действия. Поэтому в рамках ЛСМ принимаются наихудшие условия нагружения, которые получаются при совпадении направления действия статических и динамических эксплуатационных и расчётных сейсмических нагрузок. Тогда для каждого КЭ/к суммарные внутренние усилия вычисляются по формуле
&}=&} + ^п({2}о )-{2} . (40)
Дополнительно отметим, что предложенная структура ЛСМ не исключает плоский расчёт согласно [7]. Для этой цели на рис. 4 приведены графики сейсмических коэффициентов динамичности (СКД) как ВСА 7 баллов МSK-64 (см. рис. 2) при затухании (11) 0,02; 0,04 для стальных конструкций и при затухании 0,08 для железобетонных конструкций зданий с крановыми нагрузками, так и строительные СКД свода правил 14.13330 [7]. Для проектировщиков сейсмостойких кранов интерес может представить сейсмический расчёт на 8-балльное землетрясение либо его долевые составляющие от 7,1 балла с шагом 0,1 до 8 баллов, поскольку, согласно картам общего сейсмического районирования ОСР-97, на территории стран СНГ и России, семибалльные районы составляют 58 % общей площади, восьмибалльные - 28 %, девятибалльные - 14 %. Тогда, как предписывает приложение 3 РБ 006-98 [1], используют единичную нормализованную акселерограмму, которую получают делением амплитуд ВСА 7 баллов на её максимальную амплитуду. С её помощью получают модифицированные (синтезированные) ВСА в границах от 7 до 8 баллов с шагом А/ = 0,1, пересчёт которых производится по формулам, принятым по [19]:
Л^ (0 = — Лпу (0-100,3/-0Д; (41)
т 1 000 п
m = 1
m=1V l=1
(() = УЖ * ^)• 10°3/" (42)
где Л^'г - нормализованная (единичная) акселерограмма; I - балльность модифицированной акселерограммы; *ту (/) - модифицированная (синтезированная) расчётная акселерограмма, в которой переходные коэффициенты, значения которых приведены в табл. 4, имеют вид
Кп = 10°'36 - °'1. (43)
Таблица 4
Переходные коэффициенты от нормализованной к модифицированным акселерограммам МРЗ MSK-64
Баллы модифицированной (синтезированной) акселерограммы Множитель Кп (б/р) перехода от нормализованной акселерограммы
7,0 0Д0000
7,1 0,10715
7,2 0,11482
7,3 0,12303
7,4 0,13183
7,5 0,14125
7,6 0,15136
7,7 0,16218
7,8 0,17378
7,9 0,18621
8,0 0,19953
* Компонента Z в (42) дополнительно умножается на коэффициент 0,5.
Для оценки спектрального состава ВСА 7 и 8 баллов MSK-64 традиционно проводят анализ Фурье [11, 22], в результате которого, прежде всего, получают их АЧХ (см. рис. 5), из которых следует необходимость качественного определения спектра составных частот (17), особенно в области 2-10 Гц (рис. 3-5).
В заключение настоящего раздела статьи укажем на достоинство предложенного развития теории ЛСМ, которое позволяет вести расчётный анализ сейсмической безопасности пространственных металлоконструкций кранов, в отличие от плоского решения СП 14.13330 [7], смоделированных на основе любых базовых КЭ. Укажем также на недостатки предложенной методики, главный из которых состоит в том, что она завышает внутренние усилия (40) в среднем на 20-25 % и, как следствие, завышает металлоёмкость кранов при сейсмическом проектировании в соответствии с требованиями ФНП [20]. Это ставит перед проектировщиками задачу дальнейшего методологического перехода от ЛСМ к МДА [7, 16, 18, 22].
Основы метода динамического анализа
При расчётном анализе сейсмической безопасности зданий с крановыми нагрузками (рис. 6) можно отказаться от иерархических построений РДМ ЛСМ и подвергнуть сейсмическому анализу комплексную РДМ, в которой представлены все расчётные элементы действительной конструкции. Расчётно-динамическая модель здания содержит степеней свободы п = 149 862, в т. ч. крана 130/32 т - 33,5 м п = 10 950, узлов и = 1 825, КЭ 5 = 1 883, а также крана 16/3,2 т -32 м п = 6 156, и = 1 026, 5 = 987.
Рис. 6. Расчётно-динамическая конечно-элементная модель промышленного здания с крановыми нагрузками, оборудованного мостовыми кранами 130/32 т (3 шт.) и 16/3,2 т (2 шт.): число степеней свободы п = 149 862, узлов и = 24 977, КЭ 5 = 28 480, в т. ч. оболочных 17 960, стержневых 10 520
Для проведения динамического расчётного анализа пространственных конструкций, аппроксимированных дискретной конечно-элементной сеткой стержней, пластин, оболочек и трёхмерных КЭ, как это следует из теории совместных конечных элементов в рамках алгоритма МКЭ, на необходимое для расчётов расчётное сочетание нагрузок можно, с учётом работы [23], сформировать математическую конечно-элементную модель кранового сооружения с п степенями свободы:
([М1 + [МI ^ )}+ Уз (М ]К Г V (')}+ [К № }= К }+ Рдин ]- (М1 + МIШ)}, (44)
где уз, согласно (12), - коэффициент потерь.
Векторы в правой части уравнения движения (44) обозначают, как в (1), внешние статические, динамические и кинематические воздействия соответственно, в которых Лту (/) может обозначать ускорение переезда краном стыков и просадок рельсовых путей, акселерограммы землетрясений [16], а выражение перед вектором скоростей {V(')} - матрица затухания Цейтлина-Мартемьянова [24]. От матричного линейного уравнения (44) можно перейти к нелинейному расчётному анализу. Для этого в левую часть (44) необходимо ввести вектор {^ (V, V)}, характеризующий силы трения заторможенных ходовых колёс, включающиеся связи реборд ходовых колёс и односторонние связи, в частности стальные канаты грузовых полиспастов, демпферы и пр. После этого становится очевидным, что комплексная РДМ промздания должна быть загружена сочетанием нагрузок согласно методу предельных состояний ГОСТ 28609-90, в котором - акселерограмма расчётного землетрясения интенсивностью согласно картам ОСР-97 [19], ОСР-97А, ОСР-97В, ОСР-97С, ОСР-97D, характеризующим разный уровень вероятности превышения сейсмической интенсивности, указанной в баллах шкалы MSK-64 на каждой из карт. Так, карта А соответствует 10 %, карта В - 5 %, карта С - 1 %, карта D - 0,5 % вероятности возможного превышения или 90, 95, 99 и 99,5 % непревышения соответственно в течение 50 лет расчётной сейсмической интенсивности, что соответствует повторяемости сейсмического эффекта на земной поверхности в среднем один раз в 500, 1 000, 5 000 и 10 000 лет.
В качестве способа интегрирования системы сейсмических нелинейных уравнений движения (44) использован жёсткоустойчивый метод Гира, как метод динамического анализа, в виде формул дифференцирования назад [23], позволяющий контролировать знак производной дf /ду на каждом шаге интегрирования и эффективно строить алгоритм интегрирования системы дифференциальных уравнений (44).
Согласно методу Гира в общем виде нелинейное матричное уравнение движения типа (44) преобразуется в систему уравнений порядка 2п вида:
$[ £ Г }=г},
'[М ]{г } = -[С ]{Г }-[ К ]{У }- ^(У V)} + {Р}, где [Е] - единичная матрица. Далее система (45) с учётом вектора фазовых координат
7={{У },{У })т ={{У }}т
(45)
(46)
записывается в виде
или в развёрнутом виде
[ Я]{У }=[ вт+т
(ж, = У, I > 1, 2, ..., п) :
(47)
[ Щ [0]
_[0] [М]
{у } {у }
[0] [Щ] -[К] -[С].
^Ц 0 А Цу }(+1{р}-{Л (у )}]
(48)
В отличие от одношаговых методов, в (48) нет необходимости вычислять матрицу, обратную матрице масс [М]Л, что, как известно, существенно повышает точность расчёта. Решение уравнения (47) ищется в виде итерационного сходящегося процесса
7(0) _ р7 . Zn _PZn-1'
znv+1) = znv) - ^ (znv))'
(49)
(50)
где 2п - вектор Нордсика [25], представленный в транспонированном (Т) виде:
^ =
¡7}', „{V}п, 2,г {У}", 6,г {Т}", .., ^ {У<»}
(51)
где к - порядок метода Гира (максимальное число первых членов ряда Тейлора полученного решения, совпадающих с точным решением дифференциального уравнения движения); Р - треугольная матрица Паскаля порядка 2п:
Р =
11111 12 3 4
1 3 6
1
к (к -1)
к (к -1) к к 1
(52)
где I - вектор вида {I}= {10, /1, ..., 1к}т с постоянными коэффициентами, зависящими от порядка метода Гира, приведёнными в табл. 5.
т
1
0
1
Таблица 5
Порядок метода Гира
Показатель Коэффициенты^ Порядок метода Гира k
2 3 4 5 6
¿0 2/3 6/11 12/25 60/137 20/49
ii 1 1 1 1 1
¿2 1/3 6/11 7/10 225/274 58/63
¿3 - 1/11 1/5 85/274 5/12
¿4 - - 1/50 15/274 25/252
¿5 - - - 1/274 1/84
¿6 - - - - 1/1 764
В формуле (50) F (2п у)) - функция невязки:
F= НЕ(Хп,({7}п)(- Н({У}. (53)
Здесь ({7}")(у) и ({7} ) - решения, полученные на у-ом шаге итерации. Кроме того, в (50) Ж - итерационная матрица вида
( дЕ (^ ) V1
W =
dZ
■I
(54)
Очевидно, что если правая часть уравнения (47) линейна вследствие того, что в уравнении типа (44) вектор {я (у, У)} = 0 , то итерационный процесс (49), (50) сходится за одну итерацию.
Кроме того, метод Гира, как и МДА, имеет средства контроля погрешностей вычислений на каждом шаге, что позволяет строить адаптивные вычислительные процессы с автоматическим выбором шага интегрирования и порядка метода интегрирования. Это позволяет решать как линейные, так и нелинейные задачи динамического анализа.
Заключение
Основным этапом расчёта сложных пространственных систем крановых сооружений произвольного вида, составленных из тонкостенных стержней (пластин, оболочек, трёхмерных КЭ), методом конечных элементов в перемещениях является формирование матриц жёсткости, масс и демпфирования отдельного конечного элемента, находящегося в условиях пространственного деформирования в МСК. Причём динамический анализ полной системы с распределёнными массами обычно требует более сложных вычислений, чем системы с сосредоточенными массами, поскольку матрица сосредоточенных масс - диагональная, а матрица распределённых масс имеет много недиагональных элементов. Это приводит к так называемой инерционной взаимосвязи масс действительных конструкций крановых сооружений, что принято считать важным для проектировочных задач. При этом важно обосновать приоритет теории МДА сейсмостойкости перед ЛСМ, что обуславливает особые требования к назначению расчётных сейсмических воздействий, как правило, представленных трёхкомпонентными сейсмическими спектрами ответа (либо СКД), акселерограммами и сейсмическими волновыми процессами.
СПИСОК ЛИТЕРА ТУРЫ
1. РБ-006-98. Определение исходных сейсмических колебаний грунта для проектных основ. М.: Госатомнадзор России, 2000. 74 с.
2. Материалы для экспериментального проектирования и технических исследований сейсмостойких сооружений. ГИПРОТИС Госстроя СССР. Таблицы ускорений грунта прошедших землетрясений интенсивностью 7 и 8 баллов. М., 1961. Вып. 1. 200 с.
3. Полевые инженерно-геологические изыскания с целью уточнения несущих свойств грунтов площадки строительства ХОТ-2. Договор № 12731-Г от 01.06.2009. Этап 2: «Геофизические измерения с целью сейсмомикрорайонирования площадки». Иркутск: ЗАО «Востсибтизис», 2009. 190 с.
4. Center for Engineering Strong Motion Data. United States Geological Survey. URL: http:// strongmo-tioncenter.org/ (дата обращения: 15.02.2017).
5. PEER-NGA. URL: https://ngawest2.berkeley.edu/ (дата обращения: 16.02.2017).
6. Сейсмическое микрорайонирование и разработка синтезированных акселерограмм с учётом характеристик грунтов в основании здания № 1 завода РТ-2. Отчёт по договору с ГХК № 10-10/3-2005 от 15.12.05 г. М.: ИГЭ РАН, 2006. В 4-х т. Т. 4. Арх. № 5794 от 18.08.2006. 103 с.
7. СП 14.1330, 2014. Строительство в сейсмических районах. Актуализированная редакция СНиП II-7-81*. М.: Минстрой РФ, 2014. 125 с.
8. Панасенко Н. Н. Динамика и сейсмостойкость подъёмно-транспортного оборудования атомных станций: дис. ... д-ра техн. наук. В 2-х ч. Новочеркасск: НГТУ, 1992. Ч. 1. 85 с.
9. Синельщиков А. В. Динамика и сейсмостойкость мостовых кранов: дис. ... канд. техн. наук. Астрахань: АГТУ, 2000. 276 с.
10. РД10-112-1-2004. Рекомендации по экспертному обследованию грузоподъёмных машин. Общие положения. М.: Ростехнадзор, 2004. 20 с.
11. Панасенко Н. Н. Вероятностно-статистическая модель расчётного сейсмического воздействия // Механики XXI веку. 2016. № 15. С. 263-277.
12. Радчикова Л. И. К вопросу о пересчёте акселерограмм сильных движений на различные грунтовые условия // Исследования по сейсмической опасности. Вопросы инженерной сейсмологии. М.: Наука, 1988. Вып. 29. С. 81-93.
13. Мацеля В. И., Сеелев И. Н., Леконцев А. В., Хафизов Р. Р., Панасенко Н. Н., Синельщиков А. В., Яковлев П. В. Сравнительный анализ параметров конечно-элементных моделей грунтов, полученных численными методами // Вестн. Астрахан. гос. техн. ун-та. 2017. № 1 (63). С. 23-31.
14. Панасенко Н. Н., Рабей В. В., Синельщикова Л. С. Конечно-элементная модель демпфирования колебаний несущих металлоконструкций грузоподъёмных кранов // Вестн. Астрахан. гос. техн. ун-та. 2013. № 2 (56). С. 41-49.
15. Панасенко Н. Н., Синельщиков А. В. Метод конечных элементов в теории сооружений из тонкостенных стержней // Строительная механика тонкостенных стержней. Волгоград: ВолГУ, 2013. С. 228-288.
16. Panasenko N. N., Sinelshicov A. V., Rabey V. V. Calculated Justification of Seismic Stability of Load-Lifting Cranes// WSEAS Transaction on Applied and Theoretical Mechanics. 2014. Vol. 9. P. 104-123.
17. Панасенко Н. Н., Синельщиков А. В. Расчётное обоснование сейсмостойких грузоподъёмных кранов. В 3-х ч. Подъёмные сооружения и специальная техника (Украина). Одесса, 2010. Ч. 1. № 10. С. 23-26.
18. Бирбраер А. Н. Прочность и надёжность конструкций АЭС при особых динамических воздействиях. М.: Энергоатомиздат, 1989. 304 с.
19. Уломов В. И. Вероятностно-детерминированная оценка сейсмических воздействий на основе карт ОСР-97 и сценарных землетрясений // Сейсмостойкое строительство. Безопасность сооружений. 2005. № 4. С. 60-68.
20. Правила безопасности ОПО, на которых используются подъёмные сооружения. СПб: Изд-во ДЕАН, 2016. 160 с.
21. НП-043-11. Правила устройства и безопасной эксплуатации грузоподъёмных кранов для ОИАЭ М.: Ростехнадзор, 2013. 15 с.
22. Синельщиков А. В., Панасенко Н. Н., Синельщикова Л. С. Спектральный анализ Фурье вероятностно-статистических акселерограмм // Вестн. Астрахан. гос. техн. ун-та. 2012. № 2 (54). С. 22-31.
23. Синельщиков А. В., Хальфин М. Н. Дискретные методы динамического анализа грузоподъёмных кранов // Изв. вузов. Северо-Кавказ. регион. Сер.: Технические науки. 2007. № 3. С. 34-38.
24. Мартемьянов А. И. Проектирование и строительство зданий и сооружений в сейсмических районах. М.: Стройиздат, 1985. 253 с.
25. Car C. W. Numerical initialvalye problems in ordinary differential equations. New Jersey: Prentice-Hall, Englewood Cliffs, 1971. 350 p.
Статья поступила в редакцию 01.11.2017
ИНФОРМАЦИЯ ОБ АВТОРАХ
Панасенко Николай Никитович - Россия, 414056, Астрахань; Астраханский государственный технический университет; д-р техн. наук, профессор; профессор кафедры техники и технологии наземного транспорта; [email protected].
Синельщиков Алексей Владимирович - Россия, 414056, Астрахань; Астраханский государственный архитектурно-строительный университет; канд. техн. наук, доцент; зав. кафедрой прикладной механики и графики; [email protected].
Яковлев Павел Викторович - Россия, 414056, Астрахань; Астраханский государственный технический университет; д-р техн. наук, профессор; профессор кафедры безопасности жизнедеятельности и гидромеханики; [email protected].
N. N. Panasenko, A. V. Sinelshikov, P. V. Yakovlev
CALCULATED ANALYSIS OF SEISMIC SECURITY OF HOISTING CRANES
Abstract. The choice of methods of the theory of seismic resistance in relation to the assessment of seismic safety of hoisting cranes operated in buildings and structures depends on the availability and the type of a seismic model. In the presence of a seismic input in the form of seismic coefficients of dynamism or their modified variants in the form of seismic response spectra, there is used a linear-spectral method, to be exact, its "flat" version, according to BD 14.13330.2014 "Building in seismic zones" and DS-031-01 "Design standard of earthquake resistant atomic power stations". Dynamic analysis of the theory of seismic resistance of structures is an extension of a linear-spectral method. It simplifies finite-element modelling of real structures, but uses methods of numerical integration of big systems of second-order seismic differential equations.
Key words: seismic impact, 3-component accelerogram, synthesized (probability-statistical) accelerogram, seismic response spectrum, linear-spectral method, dynamic analysis method, mathematical model of seismic impact.
REFERENCES
1. RB-006-98. Opredelenie iskhodnykh seismicheskikh kolebanii grunta dlia proektnykh osnov [SG-006-98. Determining starting seismic vibrations of the ground for design data]. Moscow, Gosatomnadzor Rossii Publ., 2000. 74 p.
2. Materialy dlia eksperimental'nogo proektirovaniia i tekhnicheskikh issledovanii seismostoikikh sooru-zhenii. GIPROTIS Gosstroia SSSR. Tablitsy uskorenii grunta proshedshikh zemletriasenii intensivnost'iu 7 i 8 ballov [Documents for experimental design and engineering of seismic resistant structures. State Institute of standard and experimental design of State Building Committee of the USSR. Tables of ground accelerations of earthquakes with magnitude 7.0-8.0] Moscow, 1961, vol. 1. 200 p.
3. Polevye inzhenerno-geologicheskie izyskaniia s tsel'iu utochneniia nesushchikh svoistv gruntov ploshchadki stroitel'stva KhOT-2. Dogovor № 12731-G ot 01.06.2009. Etap 2: «Geofizicheskie izmereniia s tsel'iu seismomikroraionirovaniia ploshchadki» [Field geotechnical surveys for specifying ground bearing capacity on construction site KhOT-2. Contract № 12731-G of 01.06.2009. II stage: "Geophysical measuring for determining seismic stability of the site"]. Irkutsk, ZAO «Vostsibtizis» Publ., 2009. 190 p.
4. Center for Engineering Strong Motion Data. United States Geological Survey. Available at: http://strongmotioncenter.org/ (accessed: 15.02.2017).
5. PEER-NGA. Available at: https://ngawest2.berkeley.edu/ (accessed: 16.02.2017).
6. Seismicheskoe mikroraionirovanie i razrabotka sintezirovannykh akselerogramm s uchetom kharakter-istik gruntov v osnovanii zdaniia № 1 zavoda RT-2 [Seismic microzoning and developing synthesized accelerograms subject to ground parameters in the basement of the building № 1 of the plant RT-2]. Otchet po dogovoru s GKhK № 10-10/3-2005 ot 15.12.05 g. V 4-kh tomakh. Moscow, IGE RAN Publ., 2006. Vol. 4, no. 5794, 18.08.2006. 103 p.
7. SP 14.1330, 2014. Stroitel'stvo v seismicheskikh raionakh. Aktualizirovannaia redaktsiia SNiP II-7-81* [DS 14.1330, 2014. Building in seismic zones. Updated edition of Building Standards and Rules II-7-81*]. Moscow, Minstroi RF Publ., 2014. 125 p.
8. Panasenko N. N. Dinamika i seismostoikost' pod"emno-transportnogo oborudovaniia atomnykh stantsii. Dissertatsiia ... doktora tekhnicheskikh nauk [Dynamics and seismic stability of carrying and lifting equipment at atomic power stations. Dissertation. of Ph. Doctor].Novocherkassk, NGTU, 1992. Part 1. 285 p.
9. Sinel'shchikov A. V. Dinamika i seismostoikost' mostovykh kranov. Dissertatsiia ... kandidata tekhnicheskikh nauk [Dynamics and seismic stability of bridge cranes. Dissertation ... of Candidate of technical sciences]. Astrakhan, AGTU, 2000. 276 p.
10. RD 10-112-1-2004. Rekomendatsii po ekspertnomu obsledovaniiu gruzopod"emnykh mashin. Obshchie polozheniia [RD 10-112-1-2004. Recommendations on expert examination of lifting mechanisms. General statements]. Moscow, Rostekhnadzor Publ., 2004. 20 p.
11. Panasenko N. N. Veroiatnostno-statisticheskaia model' raschetnogo seismicheskogo vozdeistviia [Probabilistic statistical model of calculated seismic impact]. MekhanikiXXIveku, 2016, no. 15, pp. 263-277.
12. Radchikova L. I. K voprosu o pereschete akselerogramm sil'nykh dvizhenii na razlichnye gruntovye usloviia [To the problem of recalculating accelerograms of intensive movements onto different conditions of the soil]. Issledovaniia po seismicheskoi opasnosti. Voprosy inzhenernoi seismologii. Moscow, Nauka Publ., 1988. Vol. 29. Pp. 81-93.
13. Matselia V. I., Seelev I. N., Lekontsev A. V., Khafizov R. R., Panasenko N. N., Sinel'shchikov A. V., Iakovlev P. V. Sravnitel'nyi analiz parametrov konechno-elementnykh modelei gruntov, poluchennykh chislennymi metodami [Comparative analysis of elements of finite-element models of grounds, obtained using numerical methods]. VestnikAstrakhanskogo gosudarstvennogo tekhnicheskogo universiteta, 2017, no. 1 (63), pp. 23-31.
14. Panasenko N. N., Rabei V. V., Sinel'shchikova L. S. Konechno-elementnaia model' dempfirovaniia kolebanii nesushchikh metallokonstruktsii gruzopod"emnykh kranov [Finite-element model of vibration damping of load-bearing metalworks of cranes]. Vestnik Astrakhanskogo gosudarstvennogo tekhnicheskogo universiteta, 2013, no. 2 (56), pp. 41-49.
15. Panasenko N. N., Sinel'shchikov A. V. Metod konechnykh elementov v teorii sooruzhenii iz tonkosten-nykh sterzhnei [Finite-element method in theory of thin-walled bar structures]. Stroitel'naia mekhanika tonkostennykh sterzhnei. Volgograd, VolGU Publ., 2013. Pp. 228-288.
16. Panasenko N. N., Sinelshicov A. V., Rabey V. V. Calculated Justification of Seismic Stability of Load-Lifting Cranes. WSEAS Transaction on Applied and ThejreticalMechanics, 2014, vol. 9, pp. 104-123.
17. Panasenko N. N., Sinel'shchikov A. V. Raschetnoe obosnovanie seismostoikikh gruzopod"emnykh kranov. V 3-kh chastiakh Pod"emnye sooruzheniia i spetsial'naia tekhnika (Ukraina) [Design feasibility analysis of seismic resistant cranes. In 3 volumes. Lifting constructions and special equipment]. Odessa, 2010, part 1, no. 10, pp. 23-26.
18. Birbraer A. N. Prochnost' i nadezhnost' konstruktsii AES pri osobykh dinamicheskikh vozdeistviiakh [Strength and stability of structures of atomic power stations in special dynamic conditions]. Moscow, Ener-goatomizdat Publ, 1989. 304 p.
19. Ulomov V. I. Veroiatnostno-determinirovannaia otsenka seismicheskikh vozdeistvii na osnove kart OSR-97 i stsenarnykh zemletriasenii [Probabilistically determined assessment of seismic impact using maps of OSR-97 and scenario earthquakes]. Seismostoikoe stroitel'stvo. Bezopasnost'sooruzhenii, 2005. no. 4, pp. 60-68.
20. Pravila bezopasnosti OPO, na kotorykh ispol'zuiutsia pod"emnye sooruzheniia [Safety rules for dangerous production facilities with lifting constructions]. Saint Petersburg, DEAN Publ., 2016. 160 p.
21. NP-043-11. Pravila ustroistva i bezopasnoi ekspluatatsii gruzopod"emnykh kranov dlia OIAE. [SR-043-11. Rules of design and operating of hoisting cranes at nuclear facilities]. Moscow, Rostekhnadzor Publ., 2013. 15 p.
22. Sinel'shchikov A. V., Panasenko N. N., Sinel'shchikova L. S. Spektral'nyi analiz Fur'e veroiatnostno-statisticheskikh akselerogramm [Fourier spectral analysis of probabilistic and statistical accelerogramsy. Vestnik Astrakhanskogo gosudarstvennogo tekhnicheskogo universiteta, 2012, no. 2 (54), pp. 22-31.
23. Sinel'shchikov A. V., Khal'fin M. N. Diskretnye metody dinamicheskogo analiza gruzopod"emnykh kranov [Discrete methods of dynamic analysis of hoisting cranes]. Izvestiia vuzov. Severo-Kavkazskii region. Seriia: Tekhnicheskie nauki, 2007, no. 3, pp. 34-38.
24. Martem'ianov A. I. Proektirovanie i stroitel'stvo zdanii i sooruzhenii v seismicheskikh raionakh [Design andconstruction of residential buildings and structures in seismic zones]. Moscow, Stroiizdat Publ., 1985. 253 p.
25. Car C. W. Numerical initialvalye problems in ordinary differential equations. New Jersey: Prentice-Hall, Englewood Cliffs, 1971. 350 p.
The article submitted to the editors 01.11.2017
INFORMATION ABOUT THE AUTHORS
Panasenko Nickolay Nikitovich - Russia, 414056, Astrakhan; Astrakhan State Technical University; Doctor of Technical Sciences, Professor; Professor of the Department of Technique and Technology of Land Transport; [email protected].
Sinelshchikov Alexey Vladimirovich - Russia, 414056, Astrakhan; Astrakhan State University of Architecture and Civil Engineering; Candidate of Technical Sciences, Assistant Professor; Head of the Department of Applied Mechanics and Graphics; [email protected].
Yakovlev Pavel Viktorovich - Russia, 414056, Astrakhan; Astrakhan State Technical University; Doctor of Technical Sciences, Professor; Professor of the Department Safety of Vital Functions and Hydromechanics; [email protected].