Оценка напряженно-деформированного состояния оснований и грунтовых сооружений при статических и сейсмических воздействиях тема диссертации и автореферата по ВАК РФ 05.23.02, кандидат технических наук Мишин, Дмитрий Вячеславович

• Мишин, Дмитрий Вячеславович
• кандидат технических наук
• 2004

Оглавление диссертации кандидат технических наук Мишин, Дмитрий Вячеславович

ВВЕДЕНИЕ.

ГЛАВА 1. РАЗРАБОТКА МЕТОДОВ И АЛГОРИТМОВ РЕШЕНИЯ ОТДЕЛЬНЫХ ПРИКЛАДНЫХ ЗАДАЧ, ВХОДЯЩИХ В ПРОГРАММНЫЙ КОМПЛЕКС «ДИСК-ГЕОМЕХАНИКА».

1.1. Расчет НДС грунтовых массивов при статических нагрузках.

1.1.1. Математическая постановка задачи.

1.1.2. Модель упругопластического деформирования.

1.1.3. Метод численного решения.

1.2. Расчет НДС грунтовых массивов при сейсмических нагрузках.

1.2.1. Математическая постановка задачи.

1.2.2. Метод численного решения.

ГЛАВА 2. МОДЕЛИРОВАНИЕ СЕЙСМИЧЕСКИХ ВОЗДЕЙСТВИЙ.

2.1. Определение ускорения условной нижней границы основания.

2.2. Корректирование акселерограмм при моделировании сейсмических воздействий.

ГЛАВА 3. ПРОГРАММНЫЙ КОМПЛЕКС <4ЩСК-ГЕОМЕХАНИКА».

3.1. Общая программная архитектура и интерактивная среда комплекса «ДИСК-Геомеханика».

3.2. Краткая характеристика основных прикладных программ комплекса.

3.2.1. Пакет программ "Статика".

3.2.2. Пакет программ "Динамика".

3.3. Примеры тестовых расчетов.

3.3.1. Плоская деформация.

3.3.2. Упругопластическая плоская деформация.

3.3.3. Консолидация слоя.

3.3.4. Распространение волны ускорения.

3.3.5. Расчет ускорений нижней границы.

3.3.6. Корректирование акселерограмм.

ГЛАВА 4. ПРАКТИЧЕСКОЕ ПРИМЕНЕНИЕ ПРОГРАММНОГО КОМПЛЕКСА <<ЦИСК-ГЕОМЕХАНИКА».

4.1. Анализ сейсмической устойчивости шельфовой нефтяной платформы гравитационного типа.

4.1.1. Общая характеристика объекта исследования.

4.1.2. Статическое напряженное состояние.

4.1.3. Критерии возможности разжижения и нарушения недренированной прочности

4.1.4. Расчеты динамической реакции.

4.2. Оценка состояния системы «основные сооружения - основание» АЭС.

4.2.1. Общая характеристика объекта исследования.

4.2.2. Моделирование статического напряженного состояния.

4.2.3. Анализ сейсмостойкости основных сооружений.

4.3. Оценка сейсмостойкости плотины Ирганайского гидроузла.

4.3.1. Общая характеристика объекта исследования.

4.3.1. Анализ сейсмостойкости плотины.

Введение диссертации (часть автореферата) на тему «Оценка напряженно-деформированного состояния оснований и грунтовых сооружений при статических и сейсмических воздействиях»

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

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

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

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

Исследование напряженно-деформированного состояния энергетических сооружений и их оснований в настоящее время главным образом основано на компьютерном моделировании. Это особенно относится к сооружениям гидротехнического типа, которые строятся на водонасы-щенных основаниях сложной геологической структуры. Большую роль при расчете подобных сооружений играют различного рода статические и динамические воздействия, а также эффекты, связанные с технологией и последовательностью их возведения.

Начиная с 70-х годов, для решения рассматриваемых задач широко применяется метод конечных элементов (МКЭ). Среди преимуществ этого метода, сыгравших важную роль в его распространении, достаточно указать на его универсальность по отношению к различным классам задач и на его высокую вычислительную эффективность при использо-« вании современных вычислительных технологий.

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

На основе МКЭ был создан целый ряд программных комплексов, охватывающих широкий круг задач математической физики, теории деформируемых систем, гидромеханики, механики грунтов и т.п. К ним, например, относятся COSMOS/M, ANSYS, ABAQUS, MSC Nastran, PLAXIS, FLAC, FEMWATER и др. Среди отечественных разработок можно отметить системы SCAD, ЛИРА, СтаДиО, FEModels и др.

Характерной особенностью этих комплексов является их универсальность по отношению к рассматриваемому широкому кругу задач. Однако такая универсальность нередко оборачивается недостатком при расчете сложных задач обладающих определенной спецификой. Методы и алгоритмы, заложенные в той или иной программе, могут неэффективно работать на узком, но практически важном классе задач. Кроме того, стремление охватить широкий круг задач часто оборачивается неоправданным усложнением интерфейса с пользователем, что затрудняет освоение и практическое применение таких программ.

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

В данной работе во многом используются методические, алгоритмические и программные разработки, выполненные в течение ряда лет во Всероссийском научно-исследовательском институте им. Б.Е. Веденеева (ВНИИГ). Остановимся вкратце на этих работах.

Первой работой ВНИИГ в области конечно-элементного моделирования, воплотившейся в создание программы для ЭВМ, стала работа Ефимова Ю.Н., Сапожникова Л.Б. и Троицкого А.П. (математическая лаборатория и лаборатория сейсмостойкости) [22,24]. Этими авторами в 1969-1972гг. была разработана программа МКЭСТД для ЭВМ М-220, позволявшая решать задачи статической теории упругости и определять собственные частоты и формы для плоских тел. Программа написана в машинных кодах. Использовались треугольные элементы с линейной интерполяцией перемещений. Решения алгебраических систем выполнялось методом квадратного корня для симметричной ленточной матрицы [21]. Промежуточные данные хранились на внешних накопителях (магнитная лента и магнитный барабан). Максимальный порядок алгебраической системы в задачах статики составлял 2560 при ширине полуполосы не более 59. При решении проблемы свободных колебаний максимальный порядок системы - 512.

В 1969-1971 гг. Б.А. Шойхетом в Лаборатории оптического метода исследования напряжений разработана программа РЕШСИСТ для решения больших систем линейных алгебраических уравнений с симметричной положительно определенной матрицей [50], что является основной операцией при расчетах по МКЭ и во многом определяет их эффективность и даже выполнимость. При этом использован алгоритм матричной прогонки для блочно-трехдиагональной матрицы. Для хранения диагональных блоков матрицы используется ленточная схема, и такая же схема, но для смещенной ленты, используется для наддиагональных блоков. Программа была реализована на языке АЛГОЛ-60 для ЭВМ М-220 и позволяла решать задачи с числом неизвестных до 3600.

Алгоритм представляет интерес для решения больших систем в тех случаях, когда невозможно обойтись без использования дисковой памяти, т.к. его алгебраическая структура подсказывает естественный способ обмена с внешней памятью, доставляя удобные единицы обмена - блочные уравнения. Следует, однако, отметить, что в этом алгоритме ленточная схема хранения блоков, в особенности наддиагональных, не является достаточно эффективной т.к. в общем случае эта схема не отвечает структуре матриц, получаемых в МКЭ и их заполнению в процессе исключения неизвестных. Более адекватной была бы профильная схема хранения диагональных блоков и схема неполных столбцов для наддиагональных.

В настоящее время также используются алгоритмы скалярной факторизации, позволяющие решать системы с положительно определенными и знаконеопределенными симметричными матрицами, а также несимметричными матрицами. Для хранения используется схема разреженного профиля. Обмен с диском при этом не используется, что для ♦ двумерных задач на сегодняшний день полностью оправдано и, кроме того, позволяет получить дополнительное повышение эффективности V расчетов за счет применения специальных сортировок узлов, уменьшающих заполнение матрицы при факторизации. При реализации трехмерных задач метод матричной прогонки следует рассматривать как одну из перспективных альтернатив.

В дальнейшем на базе программы РЕШСИСТ разрабатывается ряд конечно-элементных программ. В той же лаборатории оптического метода исследования напряжений в 1974г. разработана программа решения плоской задачи упругости ПЗУ-М (Гордон Л.А., Заливако С.Б.) [11]. Программа позволяла определять напряжения в конструкциях, находящихся в плоскодеформированном или плосконапряженном состоянии под воздействием объемных и поверхностных нагрузок и температур. При этом материал конструкции - линейно-деформируемый, кусочно-t однородный, изотропный. В том же году в лаборатории строительной механики разработаны программа ОЗУ с аналогичными возможностями, но для задач в осесимметричных областях (Заливако С.Б., Храпков A.A., Готлиф A.A., Айзенберг В.И.) [12] и программа ПЛИТА для расчета рейснеровских пластин переменной толщины на винклеровском основании (Гордон Л.А., Фридман Е.Ш.) [13]. Языком программирования также служил АЛГОЛ-60 для ЭВМ М-220.

В 1977-1980гг. в лаборатории строительной механики разработан комплекс программ ДУПР (Вовкушевский A.B., Шойхет Б.А., Готлиф A.A., Люцко Е.А., Ремезова Е.В.). Программы комплекса реализованы на языке ФОРТРАН для ЕС ЭВМ. Комплекс позволял рассчитывать плоские и осесимметричные задачи упругости в т.ч. физически нелинейные с учетом эффектов растрескивания, последовательности приложения нагрузок, анизотропии и неоднородности свойств, линейные задачи теплопроводности и фильтрации с передачей результатов последних в задачу упругости для определения температурных и фильтрационных нагрузок.

В 1980-1982гг. в лабораториях строительной механики и оснований разработана программа ТУ ОС предназначенная для решения плоской задачи теории упругости с идеальными односторонними связями и двумерного уравнения Пуассона (Вовкушевский A.B., Зейлигер В.А.) [7,30]. В основу программы лег подход, разработанный первым из авторов [8]. Программа позволяла решать задачи теории упругости с учетом раскрытия швов и образования трещин, задачи теплопроводности и фильтрации, используя результаты последних, как воздействия для задачи теории упругости. Программа написана на языке АЛГОЛ для ЭВМ БЭСМ-6.

В 1982г. в лаборатории строительной механики разработана программа КРАК-82 для статического расчета плоских нелинейных стержневых железобетонных систем на упругом основании с учетом трещи-нообразования в бетоне, возможного нарушения контакта с основанием, последовательностей возведения и приложения воздействий (Готлиф A.A., Иванов И.В., Гордон Л.А.). Программа написана на языке ФОРТРАН для ЕС ЭВМ. Для решения систем уравнений также используется алгоритм [50], но переписанный на ФОРТРАН.

В том же 1982г. в лаборатории оснований разработана программа СУ Iii (Прокопович B.C.) [43] предназначенная для расчета дилатирую-щих упрутопластических грунтовых массивов с изменяющимися характеристиками в условиях плоской постановки. Математической основой программы послужил подход, предложенный Николаевским В.Н. [39]. Программа написана на языке АЛГОЛ для ЭВМ БЭСМ-6. Для решения нелинейных задач использован метод упругих решений. Возникающие при этом задачи упругости решались с помощью программы ТУОС.

Программа позволяла определять развитие напряжений, деформаций, смещений и пластических зон в процессе нагружения при приложении статических и кинематических нагрузок. Решение задач с 250-300 узлами занимало около 2-х часов машинного времени. Разработанный Про-коповичем B.C. алгоритм применяется и в настоящее время без существенных изменений.

В 1978-1985гг. Ефимовым Ю.Н. и Сапожниковым Л.Б. [25] в математическом отделе ВНИИГ разработан программный комплекс MFE для расчета сооружений и оснований методом конечных элементов. Языком разработки служил Фортран-IV для ЕС ЭВМ. Прототипом была программа МКЭСТД, разработанная теми же авторами в 1972г. для ЭВМ 11-го поколения. Программный комплекс позволял проводить расчеты статического и динамического напряженных состояний, стационарных температурных и фильтрационных полей. При этом предусматривалась возможность учета температурных деформаций, фильтрационных сил и процессов замерзания-оттаивания.

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

К 1989г. на основе ДУПР в лаборатории строительной механики ♦ ВНИИГ разрабатывается базовая версия нового комплекса, получившего название ДИСК (Заболотная В.А., Иванов И.В., Колтон Л.Г., Турчина

O.A., Шойхет Б.А.). Комплекс был предназначен для ЕС ЭВМ и в основном запрограммирован на языке PL/1. Хорошо продуманная программная архитектура, открытость по отношению к новым расширением в сочетании с широкими возможностями, заложенными авторами, предопределили его успешное внедрение и применение. Базовый вариант комплекса позволял решать плоские и осесимметричные задачи теплопроводности и упругости, плоские стационарные задачи фильтрации, статические и динамические задачи изгиба пластин и оболочек, плоские динамические задачи упругости, плоские квазистационарные задачи теплопроводности. Во всех случаях учитывались неоднородность и анизотропия. Было предусмотрено решение последовательности задач разной v физической природы на одной сетке КЭ с передачей результатов одной в другую в качестве исходных данных.

К 1993г. в лаборатории статики оснований и грунтовых сооружений комплекс ДИСК был пополнен пакетами программ, ориентированными на решение геотехнических задач и эта его модификация получила наименование ДИСК-Геомеханика (Готлиф A.A. при участии Прокоповича B.C.). В частности, в этот комплекс были включены возможности, предоставлявшиеся программой СУПТ. Новые программные блоки позволяли решать: плоские упругопластические задачи с изменяющимися характеристиками и учетом последовательности возведения, плоские задачи консолидации, плоские задачи теплопроводности с фазовым переходом, плоские и осесимметричные задачи фильтрации. »

В 1997-2000гг. в проблемной лаборатории комплекс ДИСК частично модернизируется для использования на персональных компьютерах (Вовкушевский A.B., Корсакова Л.В., Готлиф A.A., Прокопович B.C.,

Петров В.А., Озерова В.Д., Мишин Д.В.). Однако вследствие стреми-♦ тельной в то время эволюции ЭВМ и операционных систем на территории бывшего СССР жизнь этих программ оказалась недолгой.

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

С сожалением приходится констатировать, что почти все из описанных выше программ в настоящее время совершенно утрачены вслед-*•■ ствие обновления парка ЭВМ. Вместе с тем, некоторые из содержавшихся в них оригинальных алгоритмов, и по сей день весьма интересных с практической точки зрения, в силу своей специфичности не находят места в доступных универсальных системах математического моделирования.

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

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

Диссертация состоит из введения, трех глав, заключения и списка цитируемой литературы.

Заключение диссертации по теме «Основания и фундаменты, подземные сооружения», Мишин, Дмитрий Вячеславович

ЗАКЛЮЧЕНИЕ

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

Выполненные в рамках данной работы исследования позволяют сформулировать следующие основные выводы и результаты:

1. Разработаны алгоритмы решения задач статики и динамики систем "сооружение-основание" с учетом нелинейных свойств основания, последовательности возведения сооружений, изменений во времени параметров нагружения, механических характеристик и порового давления. Использование МКЭ в качестве базового численного метода позволило рассмотреть задачи в областях со сложной конфигурацией, неравномерным распределением параметров и сложным нелинейным поведением материала.

2. Разработан программный комплекс "ДИСК-Гемеханика для Windows", реализующий разработанные подходы.

3. Выполнено тестирование разработанных алгоритмов и программ путем сравнения численных результатов с аналитическими решениями ряда модельных задач. Проведенное тестирование свидетельствует о достоверности результатов, получаемых с помощью разработанных алгоритмов и программ.

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

5. Решены практические задачи:

• Оценка сейсмической устойчивости платформы для добычи нефти и газа, устанавливаемой в шельфовой материковой зоне.

• Оценка несущей способности грунтового основания энергоблока АЭС при действии статических и сейсмических нагрузок.

• Оценка сейсмической устойчивости каменно-набросной плотины Ирганайского гидроузла.

6. Разработанные методы и программы могут быть использованы при определении статического и динамического НДС оснований энергетических сооружений и грунтовых сооружений.

Список литературы диссертационного исследования кандидат технических наук Мишин, Дмитрий Вячеславович, 2004 год

1. Белкова И.Н., Глаговский В.Б., Готлиф A.A., Мишин Д.В. Оценка динамической реакции основания энергоблока при сейсмических воздействиях. // Известия ВНИИГ им. Б.Е. Веденеева. 2001. Т.239. С. 135-143.

2. Беллендир E.H., Глаговский В.Б., Готлиф A.A., Прокопович B.C. Математическое моделирование грунтовых сооружений и оснований. Изв. ВНИИГ им. Б.Е. Веденеева, Т. 231,1996, С.272 286.

3. Беллендир E.H., Глаговский В.Б., Мишин Д.В., Финагенов О.М. Оценка надежности основания энергоблока при сейсмических воздействиях. // Труды Международной конференции: Геотехника. Оценка состояния оснований и сооружений. Т.1., 2001.С. 273-280.

4. Бирбраер А.Н. Расчет конструкций на сейсмостойкость. СПб.: Наука, 1998.-255с.

5. Вержбицкий В.М. Основы численных методов. М.: Высш. шк., 2002, 840с.

6. Виноградов В.В., Мишин В.А. О коррекции акселерограмм сотрясения грунта при воздействиях природной и техногенной сейсми-ки. Сейсмостойкое строительство, 1998, №6, с. 12-15.

7. Вовкушевский A.B., Зейлигер В.А. Программа решения задачи упругости с односторонними связями и уравнения Пуассона методом конечных элементов для ЭВМ БЭСМ-6 (шифр ТУ ОС). JL, Изд. ВНИИГ, 1980,128с.

8. Вовкушевский A.B., Шойхет Б.А. Расчет массивных гидротехнических сооружений с учетом раскрытия швов.-М.:Энергоиздат, 1981,136 с.

9. Голицын Б.Б. Избранные труды, т.2. М.:Изд. АН СССР, 1960. -492с.

10. Гордон JI.A., Заливако С.Б., Шойхет Б.А. Решение плоской задачи теории упругости методом конечных элементов для ЭВМ типа "М-220". Л., Изд. ВНИИГ, 1974, 44с.

11. Гордон Л.А., Заливако С.Б., Храпков A.A. Решение осесимметрич-ной задачи теории упругости методом конечных элементов для ЭВМ типа "М-220". Л., Изд. ВНИИГ, 1974, 80с.

12. Гордон Л.А., Фридман Е.Ш. Программа расчета пластин на упругом основании винклеровского типа методом конечных элементов для ЭВМ типа "М-220". Л., Изд. ВНИИГ, 1974, 76с.

13. Горелик Л.В. Расчеты консолидации оснований и плотин из грунтовых материалов. Л.:Энергия 1975.

14. Гинзбург С.М., Добрецова И.В., Корсакова Л.В. Расчет бетонных сооружений с учетом изменения реологических свойств бетона в процессе строительства и эксплуатации. — Известия ВНИИГ им. Б.Е. Веденеева, 2000г., т. 237, с. 83 88.

15. Л.В.Горелик, Б.М.Нуллер. О расчете деформаций упругих тел с изменяющимися характеристиками.- Известия ВНИИГ им. Б.Е.Веденеева. Сборник научных трудов, т. 136, 1980, с.20-22.

16. Л.В.Горелик, Б.М.Нуллер, Б.А.Шойхет. Расчетные модели грунтов, подверженных замораживанию и оттаиванию.- Известия ВНИИГим. Б.Е.Веденеева. Сборник научных трудов, т. 151, 1981, с.66-70.126

17. Грайзер В.М. «Истинное» движение почвы в эпицентральной зоне. М.: ИФЗ АН СССР, 1984. - 198с.

18. А.Джоррдж, Д.Лю. Численное решение больших систем уравнений. М.: Мир, 1984, 333с.

19. Ефимов Ю.Н. Программа решения системы линейных алгебраических уравнений с симметричной матрицей ленточной структуры для ЭВМ типа "М-220" (Шифр РСПС). Л.: Изд. ВНИИГ, 1972, 36с.

20. ЗО.Зенкевич O.K. Метод конечных элементов в технике.- М.: Мир 1975, 541 с.

21. Качанов Л.М. Основы теории пластичности. М.: Наука 1969, 420с.

22. Колмогоров А.Н., Фомин C.B. Элементы теории функций и функционального анализа. М.: Наука, 1989, 623с.

23. Кольский Г. Волны напряжений в твердых телах.: Изд. иностр. лит., 1955.

24. Ломбардо В.Н., Муравьева Н.Г. Корректировка приборных записей смещений и ускорений при землетрясениях. Известия ВНИИГ, 1974, т. 105, с. 156-161.

25. Мишин Д.В. Корректирование акселерограмм при моделировании сейсмических воздействий. Сейсмостойкое строительство, 2001, №2, с. 17-22.

26. Миышн Д.В. Программная архитектура и интерактивная среда конечно-элементного расчетного комплекса ДИСК-Геомеханика. // Известия ВНИИГ им. Б.Е. Веденеева. 2002. Т.241. С. 193-196.

27. Мишин Д.В. Определение параметров теплопереноса при оценке надежности энергетических сооружений. // Известия ВНИИГ им. Б.Е. Веденеева. 2002. Т.241. С. 196-199.

28. Николаевский В.Н. Механические свойства грунтов и теория пластичности. М. ¡ВИНИТИ, 1972 - 88с.

29. Николаевский В.Н. Механика пористых и трещиноватых сред.-М.: Недра, 1984.-232 с.

30. Перепелица В.А. Определение истинного вида смещения почвы по сейсмограмме методом почленного интегрирования. Известия АН СССР. Физика Земли, 1972 №2. - с. 66-72.

31. Перепелица В.А. Определение истинного вида смещения «почвы» по сейсмограмме. М.: Наука, 1974. - 76с.

32. Писсанецки С. Технология разреженных матриц. М.: Мир, 1988, 410с.

33. Прокопович B.C. Расчеты на ЭВМ упругопластических деформаций грунтовых массивов. Известия ВНИИГ им. Б.Е. Веденеева. 1982. Т. 156. с. 3-6.

34. Саваренский Е.Ф., Косарев Г.Л., Садиков Ф.С. О восстановлении истинного движения почвы из записи землетрясения. Известия АН СССР. Физика Земли, 1977 №1. - с. 77-83.

35. Трапезников Л.П. Температурная трещиностойкость массивных бетонных сооружений.-М.:Энергоатомиздат, 1986 г, 272 с.

36. Уздин М.А., Сандович Т.А., Аль-Насер-Мохомад Самих Амин. Основы теории сейсмостойкости и сейсмостойкого строительства зданий и сооружений. СПб.: ВНИИГ, 1993. - 176с.

37. Флорин В.А. Основные уравнения динамики грунтовой массы. -Известия ВНИИГ. 1939. -Т.25. -с.190-196.

38. Форсайт Дж., Моулер К. Численное решение систем линейных алгебраических уравнений. — М., Мир, 1969.

39. Учет сейсмических воздействий при проектировании гидротехнических сооружений. (Пособие к разделу 5: Гидротехнические сооружения СниП II-7-81) СПб.: ВНИИГ, 1986. - 310с.

40. Шойхет Б.А. Программа решения системы алгебраических уравнений с симметричной блочно- трехдиагональной матрицей для ЭВМ типа "М-220". Л., Изд. ВНИИГ, 1974,28с.

41. Bathe К J., Wilson E.L. Numerical methods in finite element analysis. -New-Jersey, 1975.

42. Biot M.A. General theory of three-dimentional consolidation // J. Appl.Phys. 1941. - v.12, N 2.

43. Lysmer J., Kuhlemeyer L. Finite dynamic model for ifinite media. // J. Eng. Mech. 1969, с 859-877.