Расчет вертикальной жесткости фундамента с учетом взаимного влияния свай
- Авторы: Колесников А.О.1, Костюк Т.Н.1, Попов В.Н.2
-
Учреждения:
- Новосибирский государственный архитектурно-строительный университет (Сибстрин)
- Институт теоретической и прикладной механики имени С.А. Христиановича СО РАН
- Выпуск: Том 15, № 3 (2019)
- Страницы: 229-236
- Раздел: Динамика конструкций и сооружений
- URL: https://journals.rudn.ru/structural-mechanics/article/view/21415
- DOI: https://doi.org/10.22363/1815-5235-2019-15-3-229-236
Цитировать
Полный текст
Аннотация
Цели. В статье описаны исследования по определению влияния расстояния между сваями в кусте при вертикальных колебаниях фундаментов на величину динамической жесткости и сравниваются результаты, полученные при натурных испытаниях, в рамках волновой модели и согласно СП 26.13330.2012. Методы. Рассчитаны собственные частоты колебаний фундаментов из девяти свай диаметром d с высоким ростверком в зависимости от расстояний между сваями - 2 d , 3 d и 5 d . Полученные результаты сравнены с данными, полученными после серии импульсных динамических нагружений, выполненных в полунатурных условиях на моделях фундаментов. Значения скоростей поперечных волн определены непосредственно на опытном полигоне по результатам измерений. Результаты. Установлено, что увеличение расстояния между сваями в кусте приводит к повышению частот собственных колебаний свайных фундаментов. При изменении расстояния между сваями с 2 d до 5 d при вертикальных колебаниях фундаментов частота возрастает в 1,32 раза. Показано преимущество результатов расчетов в рамках волновой модели по сравнению с методом из СП 26.13330.2012 и их высокое совпадение с величинами, полученными в ходе экспериментов, что позволяет достоверно определять амплитудно-частотные характеристики фундаментов. Показатели, определяемые согласно СП 26.13330.2012, имеют значительно более низкие значения частот собственных колебаний с максимальным расхождением с экспериментальными данными в 2,7 раза и не в полной мере отражают изменение расстояния между сваями.
Полный текст
Введение 1 В настоящее время при оценке характеристик колебаний свайных фундаментов под машины с динамическими нагрузками существуют требования к точности получаемых результатов [1]. Ди- намический отклик системы «машина - фундамент - грунт» определяется ее собственной частотой и амплитудой колебаний при нормальных условиях работы машины, являющимися двумя наиболее важными параметрами, которые необходимо определить при проектировании фундамента. Если для увеличения собственной частоты и уменьшения амплитуды колебаний фундамента требуется использование свай, то для такого расчета необходимо определение жесткости и демпфирования свайного основания с использованием теории колебаний. Сваи почти всегда используются в виде кустов или групп, и эффект их взаимодействия существенно влияет на динамические характеристики фундаментов. Взаимодействие свайного фундамента с окружающим грунтом представляет собой одну из наименее изученных задач динамики фундаментов [2]. Дальнейшее усложнение этой задачи происходит от взаимодействия между сваями в группе. Если расстояние между сваями очень велико, жесткость группы может быть оценена просто путем суммирования жесткостей от одиночных свай. Сваи, находящиеся на близком расстоянии друг от друга, будут оказывать значительное взаимное влияние при динамическом нагружении за счет волны, распространяющейся по грунту от каждой сваи. В результате групповой эффект при динамическом возбуждении может сильно меняться, а отсутствие оценки взаимодействия свая - грунт - свая может сделать результат расчета свайного фундамента абсолютно неточным [3]. Концепция взаимодействия между сваями была введена в [4]. В [5] в решении для группы свай указано, что при близком расстоянии свай смещение одной сваи увеличивается за счет смещения всех других свай, а жесткость и демпфирование группы уменьшаются. В отсутствии аналитического решения для динамического взаимодействия свай, было предположено, что факторы взаимодействия могут быть определены из статического решения [4]. Однако, в стандарте Американского института бетона [6] указано, что динамический групповой эффект значительно отличается от статического группового эффекта, а групповая динамическая жесткость и демпфирование могут быть уменьшены или увеличены взаимодействием системы свая - грунт - свая. Точное теоретическое решение динамического взаимодействия сваи с грунтом затруднено, поэтому используются различные приближенные методы [5; 7-13]. Однако существует несколько достоверных подтверждений предлагаемых методов результатами динамических испытаний как на маломасштабных группах свай [5; 14-17], так и на группах полномасштабных свай [18; 19], поэтому дальнейшие исследования являются актуальными и полезными. В [13; 20] была рассмотрена математическая модель и ее решения, описывающие динамическое взаимодействие сваи с грунтом при вертикальном сдвиге в зависимости от ее расположения в составе свайной группы. Из представленных результатов следует, что динамическая жесткость при взаимодействии боковой поверхности сваи с грунтом может изменяться более чем на 20 % в зависимости от ее расположения в группе и расстояния до соседних свай. Однако какого-либо экспериментального подтверждения полученных результатов не приводится. Несмотря на расширяющееся распространение компьютерных технологий с привлечением больших пакетов программ при оценке характеристик колебаний свайных фундаментов, на практике при инженерных расчетах желательно получать надежные результаты с использованием простых выражений [21]. В связи с этим целью настоящего исследования является изучение динамического поведения группы свай при вертикальном возбуждении и проверка достоверности модели, рассмотренной в [13; 20], для описания взаимодействия висячих свай в кусте с учетом их взаимного расположения. Для решения поставленной задачи проведено комплексное исследование (с применением теоретического анализа), включающее определение собственных частот при вертикальных колебаниях свайных фундаментов в ходе экспериментов на опытном полигоне. Испытания проводились на модельных группах свай 2×2, представляющих собой стальные трубы, что снижало величину сопротивления грунта под их нижним концом и позволяло более точно оценивать жесткость боковой поверхности вертикальному сдвигу. Экспериментально определены собственные частоты фундаментов для различных расстояний между сваями. Полученные результаты сравнивались с теоретическими решениями. 1. Методика исследований С целью оценки взаимовлияния свай в ростверке на собственные частоты при вертикальных колебаниях фундаментов были проведены экспериментальные исследования, схема которых представлена на рис. 1. Испытывались три специально изготовленных свайных фундамента с высоким ростверком с различными вариантами расстояний между девятью сваями в кусте: фундамент Ф-1 - расстояние между сваями 2 диаметра, фундамент Ф-2 - расстояние между сваями 3 диаметра, фундамент Ф-3 - расстояние между сваями 5 диаметров. Все опытные фундаменты Ф-1, Ф-2, Ф-3 были выполнены в виде монолитного железобетонного ростверка с габаритными размерами 1,0×1,0 м и толщиной 20 см, опирающегося на девять жестко закрепленных свай диаметром d = 76 мм с рабочей длиной 1,4 м, изготовленных из металлических труб с толщиной стенки 3,5 мм. Опалубкой при бетонировании опытных фундаментов Ф-1, Ф-2, Ф-3 служили выровненные стенки котлована. Грунт опытного полигона до глубины 9,3 м состоял из лессовидной маловлажной твердой супеси с плотностью r = 1,70 т/м3 и модулем деформации Е = 14 МПа, подстилаемой полутвердыми суглинками. Грунтовые воды на площадке отсутствовали. Масса каждого из фундаментов равнялась 690 кг. Рис. 2. Диапазон измерений применяемой аппаратуры: - область измеряемых параметров; - виброметр АВМ-1 [Figure 2. Measuring range of the vibration equipment: - area of measuring parameters; - vibrometer АVM-1] Рис. 1. Схема экспериментальных исследований опытных фундаментов Ф-1, Ф-2, Ф-3: m - количество диаметров свай d Собственные частоты колебаний фундаментов lz при вертикальных колебаниях фундаментов определяются согласно [Figure 1. Scheme of testing foundations F-1, F-2, F-3 experimental researches: m - quantity of pile diameters d] λz = Kz , m По одинаковой методике на каждом фундаменгде Kz - вертикальная жесткость фундамента; m - те Ф-1, Ф-2, Ф-3 выполнялись серии импульсных вертикальных динамических нагружений (по 20 раз). Импульсное нагружение производилось стальным грузом, имеющим форму параллелепипеда массой 6 кг, который свободно падал на поверхность каждого испытуемого фундамента с высоты 0,5 м. масса всей установки. Вертикальная жесткость незаглубленных опытных фундаментов определяется суммированием жесткостей каждой из девяти свай в основании с учетом их волнового взаимодействия в зависимости от расстояния между их осями: I Для определения собственных частот колебаний использовалась аппаратура АВМ-1, включаz K = α-1 å i=1 Kzi , (1) ющая пять трехкоординатных датчиков виброакселерометров марки Analog Devices, аналого-цифгде α - коэффициент взаимного влияния свай; Kz - ровой преобразователь (АЦП) и персональный компьютер, позволяющая регистрировать колебания и одновременно обрабатывать получаемую информацию. вертикальная жесткость одиночной сваи, I = 9. Коэффициент взаимного влияния свай α определяется по формуле S Вибродатчик представлял собой металлический куб с установленными внутри интегральными двухканальными акселерометрами, имеющий возмож- α = w10 , Sw1 j (2) ность закрепления на верхней грани фундамента. где Sw10 - коэффициент жесткости одиночной сваи Основные характеристики виброакселерометра следующие: измерительный диапазон ±5g, спектральбез учета взаимного влияния; Sw1 j - среднее знаная плотность шума 250 μg/√Hz rms. Датчики подключались к АЦП с помощью соединительных кабелей. В качестве АЦП использовался внешний модуль L-CARD Е14-140. Диапазон измерений чение коэффициента жесткости свай в кусте. Коэффициент жесткости одиночной сваи без учета взаимного влияния определяется согласно [12] по формуле применяемой аппаратуры представлен на рис. 2. Данные измерений для каждого из фундаментов Ф-1, Ф-2, Ф-3 сравнивались с результатами Sw10 (a0 ) = 2πa0 J0 (a0 )J1 (a0 ) + Y0 (a0 )Y1 (a0 ) . 0 0 0 0 J 2 (a ) + Y 2 (a ) расчетов, полученных при использовании аналитических выражений, предложенных в [1; 13; 20]. Среднее значение коэффициента жесткости свай в кусте, согласно [13; 20], связано с безразмерной частотой колебаний a0 и расположением свай в S = 4Sw14 (a0 ) + 4Sw15 (a0 ) + Sw13 (a0 ) , ростверке. Например, для рассматриваемого в данной рабоw1 j 9 те фундамента из девяти свай Sw1 j определяется исгде Sw13 , Sw14 , Sw15 - коэффициенты жесткости, приходя из схемы расположения свай (рис. 3) по формуле веденные в [20], и равные S (a ) = 2πa J0 (a0 )J1 (a0 ) + Y0 (a0 )Y1 (a0 ) - 3 πa J0 (a0 )J1 (ma0 ) + Y0 (a0 )Y1 (ma0 ) C, w13 0 0 J 2 (a ) + Y 2 (a ) 2 0 J 2 (a ) + Y 2 (a ) 0 0 0 0 0 0 0 0 S (a ) = 2πa J0 (a0 )J1 (a0 ) + Y0 (a0 )Y1 (a0 ) - 5 πa J0 (a0 )J1 (ma0 ) + Y0 (a0 )Y1 (ma0 ) C, w14 0 0 J 2 (a ) + Y 2 (a ) 4 0 J 2 (a ) + Y 2 (a ) 0 0 0 0 0 0 0 0 S (a ) = 2πa J0 (a0 )J1 (a0 ) + Y0 (a0 )Y1 (a0 ) - 7 πa J0 (a0 )J1 (ma0 ) + Y0 (a0 )Y1 (ma0 ) C, w15 0 0 J 2 (a ) + Y 2 (a ) 8 0 J 2 (a ) + Y 2 (a ) 0 0 0 0 0 0 0 0 C = N N å J2n-1 (a0 )[Y2n-2 (a0 ) - Y2n (a0 )]- Y2n-1 (a0 )[J2n-2 (a0 ) - J2n (a0 )] n=1 , å J2n-1 (a0 )[Y2n-2 (ma0 ) - Y2n (ma0 )]- Y2n-1 (a0 )[J2n-2 (ma0 ) - J2n (ma0 )] n=1 где a0 = ωr0 / Vs ; ω - угловая частота колебаний; ния металлической сваи-трубы; rp - внешний ра- Vs - скорость поперечных волн в грунте; m - кодиус сваи; f18,1 § коэффициент, определяемый по личество диаметров свай между их осями; Jn ,Yn - табл. 1, согласно [23], зависящий от соотношений 2 функции Бесселя первого и второго рода, N = 2. L / rp и Ep / G , где L - длина сваи, G = Vs ρ - модуль сдвига грунта, ρ - плотность грунта, который прорезает свая. Коэффициенты вертикальной жесткости Таблица 1 f для висячих свай 18,1 [Table 1. Vertical stiffness coefficients 18,1 f for friction piles] Ep / G L / rp 10 000 2500 1000 500 250 10,8696 0,0021 0,0052 0,0104 0,0187 0,0332 21,7391 0,0031 0,0083 0,0166 0,0301 0,0509 32,6087 0,0042 0,0104 0,0218 0,0364 0,0571 43,4783 0,0042 0,0125 0,0260 0,0405 0,0582 46,7391 0,0052 0,0135 0,0270 0,0416 0,0582 54,3478 0,0052 0,0145 0,0281 0,0416 0,0582 65,2174 0,0062 0,0166 0,0291 0,0416 0,0582 76,0870 0,0062 0,0177 0,0301 0,0416 0,0582 86,9565 0,0073 0,0187 0,0301 0,0416 0,0582 100,0000 0,0083 0,0197 0,0301 0,0416 0,0582 Рис. 3. Схема коэффициентов жесткости опытного свайного фундамента [Figure 3. Scheme of stiffness ratios of the trial pile foundation] Вертикальная жесткость одиночной сваи определяется согласно [22] по формуле Kzi Kzi E A = p p rp f18,1 , 2. Результаты исследований в которой p E = 2,1×105 МПа · модуль упругости Результаты измерений и расчетов частот собматериала сваи; Ap - площадь поперечного сечественных колебаний фундаментов Ф-1, Ф-2, Ф-3 приведены в табл. 2 и 3 соответственно. Для каждого испытуемого фундамента производилась регистрация частот собственных колебаний и одновременный расчет по формулам (1), (2) и согласно [1] для различных вариантов размещения свай в ростверке m = 2, 3, 5. Необходимо отметить, что при расчетах использовались значения скоростей поперечных волн Vs , определенных непосредственно на опытном полигоне согласно методу, описанному в [22]. Были использованы две погруженные в грунт на расстоянии 2 м друг от друга металлические сваи «источник» и «приемник» диаметром 76 мм и длиной 1500 мм, имеющие на верхних торцах металлические пластины с закрепленными трехкоординатными вибродатчиками. При вертикальном ударе по свае-источнику регистрировались импульсы на ней и на свае-приемнике. По времени прохождения импульса между сваями высчитывалась скорость поперечных волн в грунте. Значение скорости поперечных волн для грунта опытного полигона составило Vs = 146 м/с. ет уменьшению частот вертикальных колебаний в 1,32 раза, при этом коэффициент взаимного влияния увеличивается от своего самого минимального значения, равного 1,00, до 1,32. Можно отметить, что увеличение расстояния между сваями с 2 до 3 диаметров более значительно влияет на увеличение частот собственных колебаний и коэффициента взаимного влияния и составляет около 19 % в обоих случаях. Дальнейшее увеличение расстояния между сваями приводит к менее существенному увеличению частоты собственных вертикальных колебаний и коэффициента взаимного влияния свай и равняется 11 %. 100 lz 80 60 40 20 Таблица 2 Частоты собственных колебаний фундаментов Ф-1, Ф-2, Ф-3 [Table 2. Natural frequencies of the foundations F-1, F-2, F-3] 1,4 2 3 m 4 5 а a Фундамент [Foundation] Измеренная частота, Гц [Measured frequency, Hz] Расчетная частота, Гц [Calculated frequency, Hz] СП Формулы 26.13330.2012 (1), (2) [SP [Formulas 26.13330.2012] (1), (2)] Ф-1 [F-1] 76,90 29,39 89,22 Ф-2 [F-2] 91,02 33,73 96,80 Ф-3 [F-3] 101,69 41,15 100,13 1,2 Коэффициенты взаимного влияния свай фундаментов Ф-1, Ф-2, Ф-3 Таблица 3 1,0 2 3 m 4 5 б [Table 3. Coefficients of the piles mutual influence of foundations F-1, F-2, F-3] Фундамент [Foundation] Измеренный коэффициент, a [Measured coefficient, a] Расчетный коэффициент, a [Calculated coefficient, a] СП 26.13330.2012 [SP 26.13330.2012] Формула (2) [Formula (2)] Ф-1 [F-1] 1,32 1,41 1,30 Ф-2 [F-2] 1,11 1,22 1,10 Ф-3 [F-3] 1,00 1,00 1,03 Из результатов экспериментов, проведенных на опытном полигоне, следует, что уменьшение расстояния между сваями в ростверке приводит к понижению частот собственных вертикальных колебаний свайных фундаментов и коэффициента взаимовлияния (рис. 4). Уменьшение расстояния между сваями в кусте с 5 диаметров до 2 способству- Рис. 4. Зависимость резонансной частоты lz (а) и коэффициента взаимного влияния a (б) от расстояния между сваями m: штриховые линии - данные измерений авторов (D); сплошные линии - результаты расчетов по формулам (1) и (2) ( ) и согласно СП 26.13330.2012 ( ) [Figure 4. Correlation of the resonant frequency lz (а) and the piles mutual ratio a (б) from the distance between piles m: dashed lines - the results of measurements of the authors (D); solid lines - the results of calculations according to formulas (1) and (2) ( ) and SP 26.13330.2012 ( )] Приведенные в табл. 2 и 3 результаты, полученные согласно [1], демонстрируют значительно более низкие (в 2 раза) значения рассчитанных частот собственных колебаний в сравнении с данными опытных замеров. Однако полученные в результате расчетов значения коэффициентов взаимного влияния практически полностью совпадают с опытными данными. Из чего следует, что эффект взаимного влияния свай в кусте согласно выражениям [1] для определения жесткости свайных фундаментов учитывается достаточно полно, но существуют некоторые проблемы при определении жесткости. В табл. 2, 3 и на рис. 4 отображены результаты проведения инженерных расчетов с использованием формулы (2) для определения жесткостей свайных фундаментов с различным расстоянием между сваями - случаи расположения 2d, 3d, 5d. При вычислении частот собственных колебаний свайных фундаментов с различным размещением 1.Rashidifar M.A., Rashidifar A.A., Abertavi A. Nonlinear characteristics of the pile soil system under vertical vibration // Universal Journal of Engineering Science. 2016. Vol. 4. No. 4. Pp. 59-65. 2.Novak M., El Sharnouby B. Evaluation of dynamic experiments on pile group // Jl. of Geotech. Eng. 1984. Vol. 110. Issue 6. Pp. 738-756. 3.Poulos H.G. Analysis of the settlement of pile groups // Geotechnique. 1968. Vol. 18. Issue 4. Pp. 449-471. 4.Novak M., Grigg R.F. Dynamic experiments with small pile foundations // Canadian Geotech. J. 1976. Vol. 13. Issue 4. Pp. 372-385. 5.ACI Committee 351. Foundations for dynamic equipment. American Concrete Institute, 2004. 63 p. 6.Matlock H., Foo S.H.C., Bryant L.M. Simulation свай в ростверке (рис. 3) коэффициенты Sw1 j для of lateral pile behaviour under earthquake motion // Proc. разных случаев размещения свай в группе ( j =3, 4, 5) определялись согласно [20]. Представленные значения частот подтверждают, что расчеты в рамках используемых приближений позволяют получить хорошее совпадение с экспериментальными данными. Рассматривая результаты, полученные с использованием формулы (1), различие с опытными данными оценивается от 2 до 16 %. Большая величина расхождения зафиксирована при расстоянии между сваями m = 2 для фундамента Ф-1 как следствие неточности изготовления опытного ростверка в полевых условиях. Из оценки результатов, полученных с использованием формулы (2), следует, что различий с опытными данными практически не существует, разброс для случаев при m = 2, 3, 5 составляет от 1 до 3 %. Выводы Исследование показало, что соотношения (1) и (2), учитывающие взаимное влияние свай в кусте, полученные в рамках волновой модели и применяемые для расчета динамических жесткостей при рассмотрении вертикальных колебаний свайных фундаментов, дают более точные результаты при определении резонансных частот по сравнению с методом, предложенным в [1], что подтверждается сравнением результатов с данными экспериментальных исследований, выполненных в полунатурных условиях. Частоты колебаний, найденные с использованием рассмотренных формул, показывают лучшее приближение к частотам, определенным опытным путем, по сравнению с существующими методами действующих норм.
Об авторах
Алексей Олегович Колесников
Новосибирский государственный архитектурно-строительный университет (Сибстрин)
Автор, ответственный за переписку.
Email: ao_kolesnikov@mail.ru
SPIN-код: 1225-7017
кандидат технических наук, доцент, кафедра инженерной геологии, оснований и фундаментов
Российская Федерация, 630008, Новосибирск, ул. Ленинградская, 113Татьяна Николаевна Костюк
Новосибирский государственный архитектурно-строительный университет (Сибстрин)
Email: ao_kolesnikov@mail.ru
SPIN-код: 5302-2974
магистрант, кафедра инженерной геологии, оснований и фундаментов
Российская Федерация, 630008, Новосибирск, ул. Ленинградская, 113Владимир Николаевич Попов
Институт теоретической и прикладной механики имени С.А. Христиановича СО РАН
Email: ao_kolesnikov@mail.ru
SPIN-код: 9483-9689
доктор физико-математических наук, главный научный сотрудник
Российская Федерация, 630090, Новосибирск, ул. Институтская, 4/1Список литературы
- СП 26.13330.2012. Фундаменты машин с динамическими нагрузками. Актуализированная редакция СНиП 2.02.05-87. Введ. 2013-01-01. ASCE Geotech. Eng. Div. Spec. Conf. on Earthq. Eng. and Soil Dyn. 1978. Vol. 1. Pp. 600-619.
- Sheta M., Novak M. Vertical vibration of pile groups // J. Geotech. Engng. 1982. Vol. 108. Pp. 570-590.
- Veletsos A.S., Dotson K.W. Vertical and torsional vibration of foundations in inhomogeneous media // Jl. of Geotech. Eng. 1988. Vol. 114. Issue 9. Pp. 1002-1021.
- El Naggar M.H., Novak M. Non-linear model for dynamic axial pile response // Jl. of Geotech. Eng. 1994. Vol. 120. Issue 2. Pp. 308-329.
- Gazetas G., Makris N. Dynamic pile-soil-pile interaction. Part I. Analysis of axial vibration // Earthquake Engng. ASCE. 1991. Vol. 20. Pp. 115-132.
- Баранов В.А. О расчете вынужденных колебаний заглубленного фундамента // Вопросы динамики и прочности: тр. Рижского ПИ. 1967. № 14. С. 195-209.
- Колесников А.О., Попов В.Н. Оценка динамических реакций на контурах нескольких круглых вырезов при колебаниях пластины // Строительная механика инженерных конструкций и сооружений. 2014. № 3. С. 37-43.
- Hassini S., Woods R.D. Dynamic experiments with model pile foundations // Proc. 12th ICSMFE. 1989. Vol. 2. Pp. 1135-1138.
- Han Y., Vaziri H. Dynamic response of pile groups under lateral loading // Soil Dyn. and Earthq. Eng. 1992. Vol. 11. Pp. 87-99.
- Manna B., Baidya D.K. Dynamic nonlinear response of pile foundations under vertical vibration. Theory versus experiment // Soil Dyn. and Earthq. Eng. 2010. Vol. 30. Pp. 456-469.
- Burr J.P., Pender M.J., Larkin T.J. Dynamic response of laterally excited pile groups // Jl. of Geotech. & Geoenviron. Eng. 1997. Vol. 123. No. 1. Pp. 1-8.
- Blaney G.W., Muster G.L., O'Neill M.W. Vertical vibration test of a full-scale pile group // Geotechnical Special Publication. 1987. No. 11. Pp. 149-165.
- Crouse C.B., Cheang L. Dynamic testing and analysis of pile-group foundations // Geotechnical Special Publication. 1987. No. 11. Pp. 79-98.
- Колесников А.О., Попов В.Н. Динамические реакции на контурах круговых вырезов с учетом их взаимного расположения при колебаниях пластины // Известия ВНИИГ имени Б.Е. Веденеева. 2017. Т. 283. C. 3-10.
- Prakash S., Sharma H.D. Pile Foundations in Engineering Practice. John Wiley&Sons, Inc., 1990. 784 p.
- Chowdhury I., Dasgupta S.P. Dynamics of Structure and Foundation - A Unified Approach: 1. Fundamentals. CRC Press, 2008. 882 p.
- Novak M., El Sharnouby B. Stiffness and damping constants for single piles // Jl. of Geotech. Eng. 1983. Vol. 109. Pp. 961-974.