Диплом, курсовая, контрольная работа
Помощь в написании студенческих работ

Пространственное напряженно-деформированное состояние дорожных конструкций при динамическом нагружении

Диссертация: Пространственное напряженно-деформированное состояние дорожных конструкций при динамическом нагружении
ДиссертацияПомощь в написанииУзнать стоимостьмоей работы

Полученная система 15 интегро-функциональных уравнений имеет довольно сложную структуру. Ее описанию и исследованию свойств интегральных операторов, входящих в уравнения системы, посвящен третий пункт второй главы. На основе проведенного исследования в следующем пункте предложен аппроксимационный метод решения этой системы при любом фиксированном значении параметра преобразования Фурье… Читать ещё >

Содержание

  • Глава 1. НАПРЯЖЕННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ ДОРОЖНОЙ КОНСТРУКЦИИ
    • 1. 1. Пространственное напряженное состояние при ф динамическом нагружении (постановка задачи)
    • 1. 2. Частные случаи
      • 1. 2. 1. Плоское напряженное состояние
      • 1. 2. 2. Напряженное состояние многослойного полупространства (грунта)
    • 1. 3. Принцип суперпозиции при решении задачи о пространственном напряженном состоянии дорожного покрытия
    • 1. 4. Решение вспомогательных задач для подобластей
    • 1. 5. Решение вспомогательной задачи для полупространства грунта)
  • Глава 2. РАСЧЕТ ПРОСТРАНСТВЕННОГО НДС СИСТЕМЫ «ДОРОЖНАЯ КОНСТРУКЦИЯ — ГРУНТ» ПРИ ДИНАМИЧЕСКОМ НАГРУЖЕНИИ
    • 2. 1. Алгоритм сведения задачи к системе интегро-функциональных уравнений для изолированной полосы
    • 2. 2. Пространственное напряженное состояние дорожного покрытия, жестко сцепленного с однородным полупространством (фунтом)
    • 2. 3. Свойства операторов системы интегро-функциональных уравнений
    • 2. 4. Решение системы интегро-функциональных уравнений (случай одной полосы)
    • 2. 5. Алгоритм решения пространственной краевой задачи в общем случае
    • 2. 6. Особенности практической реализации методики на
  • ПЭВМ
  • Глава 3. АНАЛИЗ РЕЗУЛЬТАТОВ ЧИСЛЕННОГО ЭКСПЕРИМЕНТА
    • 3. 1. Основные особенности пространственного НДС дорожной щ конструкции при динамическом нагружении
    • 3. 2. Упрощенные модели и их возможности
    • 3. 3. Особенности генерации волновых полей в грунте движущимися поверхностными нагрузками
      • 3. 3. 1. Плоское напряженно-деформированное состояние
      • 3. 3. 2. Пространственное напряженно-деформированное состояние
      • 3. 3. 3. Некоторые особенности распределения полей смещения от движущейся нагрузки
    • 3. 4. Сопоставление теории с экспериментом
    • 3. 5. Выводы

Пространственное напряженно-деформированное состояние дорожных конструкций при динамическом нагружении (реферат, курсовая, диплом, контрольная)

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

Существует целый ряд причин, в том числе и технологического характера, но одна из важнейших — несовершенство расчетных схем, принятых при расчете и проектировании дорожных конструкций, связанных с введением тех или иных упрощающих предположений. В силу этого проблема дальнейшего развития методов исследования системы «дорожная конструкция — грунт» остается актуальной. Так, например, наиболее распространенная в данное время расчетная схема определения необходимой толщины слоев дорожной одежды основана на модели слоистого упругого полупространства при статическом нагружении. Различные усовершенствования этой расчетной схемы основаны на абсолютизации свойств отдельных элементов и отбрасывании особенностей строения системы, которые на первый взгляд мало влияют на ее напряженно — деформированное состояние (НДС). Довольно часто при этом используются плоские модели (неизменность НДС в любом поперечном разрезе системы) и осесимметричные (осевая симметрия конструкции и прилагаемой к ней нагрузки).

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

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

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

Это связано с тем, что реальная область, охваченная возмущениями, достаточно велика.

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

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

РАН [Воровича И.И.| и Бабешко В.А.

Большой вклад в развитие теории контактных взаимодействий внесли Александров В. М., Арутюнян Н. Х., Бабешко В. А., Белоконь.

A.В., Бородачев Н. М., Ватульян А. О., Ворович И. И., Галин J1.A., Глушков Е. В., Горшков А. Г., Гринченко В. Т., Гузь А. Н., Калинчук.

B.В., Коваленко Е. В., Купрадзе В. Д., Лурье А. И., Малый. В.И., Мартыненко М. Д., Моссаковский В. И., Мхитарян С. М., Нуллер Б. М., Партон В. З., Перлин П. И., Попов Г. Я., Проценко B.C., Пряхина О. Д., Рвачев В. Л., Саркисян B.C., Сеймов В. М., Селезнев М. Г., Соболь Б. В., Тарлаковский Д. В., Угодчиков А. Г., Улитко А. Ф., Устинов Ю. А., Шерман Д.И.

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

Учет деформативных свойств штампа в основном связан с решением контактных задач для пластины, слоя или полупространства, подкрепленного стрингером (балкой), изложенных в работах Агабекяна П. В., Агаяна K. J1., Александрова В. М., Арутюняна Н. Х., Бабешко В. А., Воровича И. И., Григоряна Э. Х., Мхитаряна С. М., Пряхиной О. Д., Эрдогана Ф., Гупта Г. и мн. др. [2, 3, 9, 13, 47, 48, 61, 62, 121, 128, 166, 172].

Методы, используемые при решении задач механики контактных взаимодействий, условно можно разделить на аналитические (асимптотические, факторизации функций и матриц-функций, фиктивного поглощения) [4, 6, 7, 8, 14, 15, 17, 21, 46, 49, 51, 97, 107, 129, 130, 132], аналитико-численные (метод граничных интегральных уравнений (МГИУ), методы разложения решения по системам ортогональных функций и др.) [39, 59, 119, 120, 136, 143, 144, 145] и прямые численные (метод конечного элемента (МКЭ), метод конечных разностей (МКР) и др.) [26, 29, 30, 35, 57, 65, 115].

Использование прямых численных методов в настоящее время в основном связано с применением программных комплексов широкого или специального назначения, к которым можно отнести «ANSYS», «COSMOS», «POLUS», «ЛИРА», «СПРИНТ», «МИРАЖ», «Micro Fe» и Др.

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

Бабешко В.А., Белоконя А. В., Бреховских Л. М., Воровича И. И., Гетмана И. П., Глушкова Е. В., Гринченко В. Т., Космодамианского А. С., Ляпина А. А., Мелешко В. В., Молоткова Л. А., Петрашеня Г. И., Попова Г. Я Поручикова В. Б., Пряхиной О. Д., Сеймова В. М., Селезнева М. Г., Слепяна Л. И., Трофимчука А. Н., Улитко А. Ф., Устинова Ю. А. и др. [16, 18, 20, 22, 23, 24, 27, 32, 36, 37, 41, 42, 43, 44 45, 50, 52, 53, 55, 63, 64, 67, 69, 73, 74, 75, 76, 77, 86, 89, 98, 109, 117, 118, 122, 123, 139, 146, 147, 148, 149, 150, 154, 156, 162].

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

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

Остановимся теперь на основном содержании работы.

В первой главе приводится постановка пространственной динамической краевой задачи, моделирующей систему «дорожная конструкция-грунт». Дорожная конструкция описывается пакетом полос конечной толщины и ширины с плоскопараллельными границами, жестко сцепленных между собой и со слоистым полупространством, моделирующим грунт. На поверхности полосы в некоторой области Q действует система осциллирующих усилий прямолинейно и равномерно движущихся со скоростью V0.

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

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

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

Во второй главе изложены вопросы построения решения краевой задачи в пространственной постановке. Для этой цели используются принцип суперпозиции, метод интегральных преобразований и соотношения, полученные в первой главе. В результате решение задачи о контакте полосы с полупространством сводится к системе 15 интегро-функциональных уравнений с 15 неизвестными функциями напряжения Xj, Yj (j = 1,2.6), Zj (i = 1,2,3) три неизвестные функции Zj (i = 1,2,3) определяют напряжения вдоль границы контакта полосы и полупространства). Через указанные функции определяется НДС элементов исследуемой системы.

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

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

Шестой пункт посвящен описанию практической реализации на ПЭВМ предложенного алгоритма. Приведена блок-схема алгоритма расчета амплитудно-частотной характеристики (АЧХ) точки конструкции и дано ее подробное описание.

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

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

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

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

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

В четвертом пункте третьей главы сопоставляются расчеты с экспериментом. Приведены примеры АЧХ и амплитудно-временных характеристик (АВХ) и показано, что основные качественные характеристики расчетный АЧХ во всех случаях хорошо коррелируют с данными эксперимента.

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

Научная новизна исследования заключается в следующем: — в работе поставлена и решена новая динамическая задача строительной механики в пространственной постановке, моделирующая систему «дорожная конструкция-грунт» — разработан аналитико-численный метод решения динамической задачи строительной механики для полосы конечной ширины (дорожной конструкции), жестко сцепленной с полупространством (грунтом), в пространственной постановке;

— разработан алгоритм расчета напряжений и деформаций системы на ПЭВМ;

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

Достоверность полученных в работе результатов обеспечена:

1. Обоснованным и корректным применением математических методов исследования модельных задач при получении аналитических и численных решений.

2. Сопоставлением результатов расчета для частных и предельных случаев с полученными другими методами.

3. Сопоставлением ряда расчетных характеристик с данными натурного эксперимента.

Основные результаты, полученные в настоящей диссертационной работе, сводятся к следующему:

1. Поставлена и решена новая динамическая задача строительной механики в трехмерной постановке, которая описывает НДС элементов системы «дорожная конструкция — грунт» при воздействии движущегося автотранспорта.

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

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

4. Проведен подробный численный эксперимент для изучения основных динамических характеристик НДС элементов системы под воздействием движущейся осцилирующей нагрузки.

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

6. На основе асимптотического решения построены диаграммы направленности продольных и поперечных волн от движущейся нагрузки в полупространстве (грунте) и изучена их зависимость от скорости движения транспорта.

7. Результаты численного эксперимента были использованы при планировании широкомасштабных натурных экспериментальных исследований в ДорТрансНИИ РГСУ, а также при проведении госбюджетных и хоздоговорных НИР ДорТрансНИИ РГСУ.

8. Данные натурных исследований иллюстрируют хорошее совпадение с материалами численного эксперимента и подтверждают все основные качественные и количественные выводы, полученные в результате теоретических исследований, проведенных в работе.

9. Полученные результаты дают теоретическую основу для разработки уточненных инженерных методик расчета и проектирования автодорог с учетом основных особенностей динамического воздействия транспортного потока.

ЗАКЛЮЧЕНИЕ

.

Показать весь текст

Список литературы

  1. М., Стиган И. Справочник по специальным функциям. М.: Наука, 1979. — 830 с.
  2. П.В. Контактная задача для упругой полуплоскости, усиленной на своей границе конечной и полубесконечной накладками. Сб. Механика деформируемого твердого тела. Ереван: изд. АН Арм. ССР, 1989, С. 17−23.
  3. К.Л. Периодическая контактная задача для бесконечной пластины с упругими накладками. Изв. АН Арм. ССР, Механика 1975, т.28, N 3, С. 3−11.
  4. С.М., Александров В. М. Асимптотические решения контактных задач теории упругости для полупространства и полуплоскости, неоднородных по глубине.//Изв. АН Арм. ССР. Механика. 1986. Т.39. N 3. С.13−28.
  5. Аки К., Ричарде П. Количественная сейсмология. М.: Мир, Т.1, 1983.-519 с.
  6. В.М. Асимптотические методы в смешанных задачах теории упругости // Развитие теории контактных задач в СССР. М.: Наука, 1976. С.96−100.
  7. В.М., Мхитарян С. М. Контактные задачи для тел с тонкими покрытиями и прослойками.// М.: Наука, 1983, 488с.
  8. В.М., Пожарский Д. А. Неклассические пространственные задачи механики контактных взаимодействий упругих тел. М.: Факториал, 1998. 288 с.
  9. В.М., Солодовник М. Д. Эффективный метод решения задачи о взаимодействии накладки (стрингера) с упругой полуплоскостью и некоторые новые качественные результаты. -Труды
  10. X Всесоюзной конференции по теории оболочек и пластин., Т.1, Кутаиси, 1975, С. 10−17.
  11. Ю.Алексеев А. С., Бабич В. М., Гельчинский Б. Я. Лучевой метод вычисления интенсивности волновых фронтов.//Вопросы динамической теории распространения сейсмических волн. в.5, 1961. С.5−24.
  12. И.В. К теории колебания сред с непрерывно меняющимися характеристиками.//Изв. СКНЦ ВШ,/Естеств. науки/, № 4, 1976.
  13. Н.Х., Мхитарян С. М. Периодическая контактная задача для полуплоскости с упругими накладками. ПММ, 1969, т. ЗЗ, N5, С.813−843.
  14. В.А. Обобщенный метод факторизации в пространственных динамических смешанных задачах теории упругости. М.: Наука, 1984.-256 с.
  15. В.А. О вибрации систем штампов.//Изв. АН СССР, МТТ. № 6, 1990. С. 72−78.
  16. В.А., Белянкова Т. Н., Калинчук В. В. О решении одного класса смешанных задач для слоистого полупространства.// Докл. РАН, 2001, 380, N5, С. 619−622.
  17. В.А., Ворович И. И., Селезнев М. Г. Вибрация штампа на двухслойном основании.//ПММ, т.41, в. 1, 1977. С. 166−173.
  18. В.А., Ворович И. И., Селезнев М. Г. Распространение в упругом слое волн, возникающих при колебании штампа.//Алма
  19. Ата: Наука Каз. ССР, Сб. Распространение упругих и упруго-пластических волн. 1973. С. 339−342.
  20. В.А., Глушков Е. В., Глушкова Н. В. Об особенностях в угловых точках пространственных штампов в контактных задачах.// ДАН СССР. 1981. Т 257. N2. С.289−294.
  21. В.А., Глушков Е. В., Зинченко Ж. Ф. Динамика неоднородных линейно-упругих сред. М.: Наука, 1989. — 344 с.
  22. В.А., Румянцев А. Н. Колебания штампа, частью поверхности сцепленного с упругим слоем.// ПММ, т.41, в.4, 1977.
  23. В.А., Селезнев М. Г. и др. Об одном методе исследования установившихся колебаний упругого полупространства, содержащего сферическую или горизонтальную цилиндрическую полость. //ПММ. В.1, т. 47, 1983. С. 115−121.
  24. В.А., Селезнев М. Г., Шагинян А. С. Об одном методе уточненного учета реакции упругой среды при гармоническом воздействии.-М.: Недра, Прикладная геофизика, в.89, 1981.С.79−88.
  25. В.А., Селезнев М. Г., Шагинян А. С. Способ определения параметров смещения упругой среды при гармоническом воздействии.- М.: Недра, Прикладная геофизика, в. 106, 1983. С. 32−39.
  26. В.М., Кирпичникова Н. Я. Метод пограничного слоя в задачах дифракции коротких волн. Ленинград: ЛГУ, 1974. — 125 с.
  27. Н.С. Численные методы. Т.1. М.: Наука, 1975. -632 с.
  28. А.В. К теории динамических задач с подвижными возмущениями для неоднородной упругой полосы.//- Докл. АН СССР, т.261, № 5, 1981. С. 1079−1082.
  29. М.В., Казарян К. Б. К вопросу существования поверхностных сдвиговых волн в неоднородном упругом полупространстве. -Изв. Нац. АН Армении, Мех. 2000.-53, N1, С. 6−12.
  30. П., Баттерфилд Р. Методы граничных элементов в прикладных науках. -М: Мир, 1984. 494 с.
  31. И.С. и Жидков Н.П. Методы вычислений, т.т. 1,2. М.: Наука, 1966.
  32. С.И., Полякова И. Б. Об ограниченных В-резонансах в системе массивный штамп слоистое основание.//Изв. АН СССР, МТТ. № 6, 1990. С. 67−71.
  33. С.И., Румянцев А. Н., Селезнев М. Г. Решение задачи о возбуждении волн в упругом двухслойном полупространстве.//Сб. «Методы расширения частотного диапазона вибросейсмических колебаний», Новосибирск, ИГ и Г СО АН СССР, 1987.
  34. А.Н., Калинин И. И. Распространение волн в неоднородном слое.// Изв. вузов Сев.- Кавк. регион, Техн. н., 2000, N2, С. 1518,121.
  35. .М., Столяров С. Н. О принципах излучения в среде с дисперсией.//Проблемы теоретической физики (Сб. памяти В.Е. Тамма). М., 1972. С. 267−280.
  36. К., Теллес Ж., Вроубел JI. Методы граничных элементов. М.: Мир, 1987.-524 с.
  37. Бреховских J1.M. Волны в слоистых средах. М., 1957. 502 с.
  38. .Г. Принципы излучения, предельного поглощения и предельной амплитуды в общей теории уравнений с частными производными.// Успехи математических наук. Т.21, № 3, 1966. С. 115−194.
  39. Г. Р. Теория Бесселевых функций. М.: Иноиздат, 1949. -798 с.
  40. А.О., Шамшин В. М. Новый вариант граничных интегральных уравнений и их применение к динамическим пространственным задачам теории упругости / ПММ. -1998. -62, № 3. -С.462−469.
  41. А.С. Гибкие пластинки и оболочки. М.: Гостехиздат, 1956.
  42. Вопросы возбуждения волн вибрационными источниками.- Новосибирск, ИГ и Г СО АН СССР, 1976.
  43. И.И., Александров В. М., Бабешко В. А. Неклассические смешанные задачи теории упругости. М.: Наука, 1974.
  44. И.И., Бабешко В. А. Динамические смешанные задачи теории упругости для неклассических областей. М.: Наука, 1979.
  45. И.И., Бабешко В. А., Пряхина О. Д. Динамика массивных тел и резонансные явления в деформируемых средах. Научный мир, 1999,246 с.
  46. И.И., Белянкова Т. И., Калинчук В. В. К проблеме низкочастотных резонансов при взаимодействии упругого тела с полуограниченной средой.// Докл. РАН. 1998, Т.358, N5, C.624−626.
  47. Е.И., Пряхина О. Д. Аналитический метод определения В-резонансов//Изв. АН СССР. МТТ. № 3, 1987. С. 101−106.
  48. Е.И., Пряхина О. Д. Динамическая контактная задача для упругой системы балка-слой. // Изв. АН СССР, МТТ, 1989, N1, С. 144−148.
  49. Е.И., Пряхина О. Д. Об одном эффективном методе решения задачи о колебании упругой балки на упругом слое. // Изв. АН СССР, МТТ, 1989, N4, С. 96−101.
  50. Е.И., Пряхина О. Д., Селезнев М. Г., Тукодова О. М. Исследование взаимного влияния двух штампов при гармоническом нагружении. //Сб. Исследование по расчету пластин и оболочек. -Ростов-на Дону, РИСИ, 1987. С. 83−88.
  51. Е.И., Пряхина О. Д., Тукодова О. М. Динамические свойства упругой полуограниченной среды, контактирующей с упругим инерционным элементом. //Изв. АН СССР, МТТ. № 2, 1996. С. 128−133.
  52. JI.A. Контактные задачи теории упругости и вязкоупруго-сти.// М.: Наука, 1980, 304 с.
  53. А., Образцов М. Б. Об условиях затухания решений и принципе излучения для одного дифференциального уравнения с операторными коэффициентами на полуоси.//Дифференциальные уравнения. Т. 19., в.6, 1983. С. 944−954.
  54. И.П., Устинов Ю. А. Математическая теория нерегулярных твердых волноводов. Ростов-на-Дону, Изд-во РГУ, 1993. 143 с.
  55. Г. И. Расчет сооружений, заглубленных в грунт— М.: Стройиздат, 1977.-265с.
  56. Е.В., Кириллова Е. В. Динамическая смешанная задача для пакета упругих слоев.// Приклладная маттематика и механика (Москва),-1998,-62, N3, C.455−461.
  57. Н.В., Глушков Е. В., Хофф Р. Сингулярность напряжений в многогранных угловых точках упругих разномодульных соединений. // Докл. РАН. 1999. Т. 370. N2.
  58. О., Маршан Ж., Муратидис А. Метод конечных элементов применительно к процессу трещинообразования в дорожных покрытиях и расчету времени раскрытия трещин. //Бюллетень ЦЛДМ, № 125, 1983.
  59. A.M., Гринченко В. Т. Метод однородных решений в случае негладких нагрузок. //Теоретическая и прикладная механика. Харьков, вып. 19, 1988. С.111−116.
  60. А.Г., Тарлаковский Д. В. Динамические контактные задачи для абсолютно жестких тел и упругого полупространства.// Препринт, МАИ, 1989,49 с.
  61. И.С., Рыжик И. М. Таблицы интегралов, сумм, рядов и произведений. М.: Наука, 1971. — 1108 с.
  62. Э.Х. О двух динамических контактных задачах для полуплоскости с упругими накладками. Изв. АН СССР, МТТ, 1972, N5, С. 101−116.
  63. Э.Х. О динамической контактной задаче для полуплоскости, усиленной упругой накладкой конечной длины. ПММ, 1974, Т.38, N2, С. 321−330.
  64. В.Т. Равновесие и установившиеся колебания упругих тел конечных размеров. Киев: Наукова Думка, 1978. — 264 с.
  65. В.Т., Мелешко В. В. Гармонические колебания и волны в упругих телах. Киев: Наукова думка, 1981.- 283 с.
  66. .П., Марон И. А., Шувалов Э. З. Численные методы анализа.- М.: Наука, 1967.
  67. В.А., Прудников А. П. Интегральные преобразования и операционное исчисление. М.: Физматгиз, 1961.-367 с.
  68. Дьелесан Э, Руайе Д. Упругие волны в твердых телах. М.: Наука, 1982.-424 с.
  69. М.А. Асимптотические оценки и целые функции. М.: Физматгиз, 1962.
  70. О.Ю., Улитко А. Ф. Введение в механику нестационарных колебаний и волн. Киев: Вища шк. Головное изд-во, 1989. — 184 с.
  71. В.Б. Об асимптотическом решении плоских и осесим-метричных нестационарных динамических задач.// Современные проблемы мех. сплош. ср. Труды 6 Междунар. конф., Ростов-на-Дону, 12−14 июня, 2000, Т.2, Изд-во СКНЦВШ, 2001, С. 74−77.
  72. В.В. Теория приближенных методов. Киев: Наукова думка, 1968.
  73. Излучение и регистрация вибросейсмических сигналов.- Новосибирск, ИГ и Г СО АН СССР, 1986.
  74. С.К., Ляпин А. А. Особенности расчета напряженно-деформированного состояния конструкции дорожной одежды при динамическом нагружении. //Изв. высших учебных заведений. Северо Кавказский регион. Технические науки, № 4, 1997. С. 63−66.
  75. С.К., Селезнев М. Г. О разработке новых современных методов расчета и конструирования дорожных одежд. //М.: Дороги. Наука и техника в дорожной отрасли. № 1, 2000. С.7−11.
  76. С.К., Селезнев М. Г., Ляпин А. А., Углова Е. В. Распределение энергии воздействия движущегося транспорта в элементах системы «дорожная конструкция грунт». //М.: Дороги. Наука и техника в дорожной отрасли. № 4, 2001. С.8−10.
  77. С.К., Селезнев М. Г., Углова Е. В. Динамика дорожных конструкций. Ростов н/Д: РГСУ, 2002. — 258 с.
  78. С.К., Селезнев М. Г., Углова Е. В. Необходимо разработать новые критерии расчета и конструирования дорожных одежд. //М.: Дороги. Наука и техника в дорожной отрасли. № 3, 2000. С. 1315.
  79. Исследование Земли невзрывными источниками. М.: Наука, 1981.
  80. В.Д. Направления научных исследований в связи с концепцией развития дорожной отрасли. Приложения к журналу «Автомобильные дороги», № 1, 1997. С. 2−4.
  81. В.Д. Задачи совершенствования теории и практики расчета и конструирования дорожных одежд. //Автомобильные дороги, № 1, 1992. С. 8−10.
  82. В.В., Белянкова Т. И. К проблеме исследования особенностей динамического контактного взаимодействия штампа с полупространством, ослабленным наличием трещины.// Сев.-Кав. регион. Естеств.н. 2001, Спец.вып., С. 83−85, 171.
  83. В.В., Селезнев М. Г. Некоторые особенности возбуждения и распространения упругих волн в неоднородных средах. //Сб. «Разработка и исследование источников сейсмических сигналов» -М.: ВНИИОЭНГ, 1986. С. 61−66.
  84. Г., Корн Т. Справочник по высшей математике для научных работников и инженеров. М.: Наука 1970.
  85. Л.Я. О распространении упругих волн в двухкомпо-нентных средах.// ПММ, 23, № 6, 1959. С. 1115−1123.
  86. А.С., Сторожев В. И. Динамические задачи теории упругости для анизотропных сред. Киев: Наукова думка, 1985.- 176 с.
  87. П.Е. О возбуждении нормальных и присоединенных волн в бесконечной слоистой упругой полосе.// ПММ. 1979. т.43, № 5, с. 877−886.
  88. В.Д. Нестационарное взаимодействие элементов конструкций со средой. Киев: Наукова думка, 1979. — 184 с.
  89. В.Д., Гегелиа Т. Г., Башейлешвили М. О., Бурчуладзе Т. В. Трехмерные задачи математической теории упругости и термоупругости. М.: Наука, 1976. — 603 с.
  90. М.А., Шабат Б. В. Методы теории функций комплексного переменного. М.: Наука, 1965.
  91. Н.Н. Специальные функции и их приложения. М.: Физ-матгиз. 1963.-358 с.
  92. А.И. Пространственные задачи теории упругости. М.: Гос-техиздат, 1955,-491 с.
  93. А.И. Теория упругости. М.: Наука, 1970.
  94. А.А., Румянцев А. Н., Селезнев М. Г. Особенности нестационарного воздействия массивного штампа на двухслойное полупространство с заглубленной полостью.// Изв. АН СССР, МТТ, № 6, 1990.
  95. А.А., Селезнев М. Г., Собисевич Л. Е., Собисевич А. Л. Механико-математические модели в задачах активнойсейсмологии. М. ТНИЦ ПГК Минобразования России, 1999, 291 с.
  96. Т.А. Особенности расчета кратных слабосходящихся несобственных интегралов, возникающих в пространственных динамических задачах теории упругости для полуограниченных областей.// Сборник «Известия РГСУ», № 3, г. Ростов-на-Дону, 1998, С. 194.
  97. Т.А. Динамическая задача контакта упругой полосы с полупространством.// Тезисы Международной юбилейной научно-практической конференции «Строительство-99», РГСУ, г. Ростов-на-Дону, 1999, С. 50.
  98. Т.А. Некоторые особенности распределения полей смещения от движущейся нагрузки.// Материалы Международной научно-практической конференции «Строительство-2002», РГСУ, г. Ростов-на-Дону, 2002, С. 67−68.
  99. Т.А., Селезнев М. Г. Об исследовании волновых полей от движущихся нагрузок в многослойном полупространстве.// Тезисы Международной научно-практической конференции «Строительство-98», РГСУ, г. Ростов-на-Дону, 1998, С. 165−166.
  100. Т.А., Селезнев М. Г. Об анализе волновых полей от движущейся нагрузки.// Межвузовский сборник научных трудов «Технологические процессы в транспортном машиностроении», РГУПС, г. Ростов-на-Дону, 1997, С. 94−97.
  101. Т.А., Селезнев М. Г. Об особенностях решения одной системы интегральных уравнений.// Тезисы Международной научно-практической конференции РГСУ, г. Ростов-на-Дону, 1997.
  102. Т.А., Селезнев М. Г., Собисевич A.JI. Об анализе волновых полей от движущейся нагрузки в слоистом полупространстве.// М.: Доклады Академии Наук, т.367, № 1,1999, С. 52−55.
  103. Методы расширения частотного диапазона вибросейсмических колебаний. Новосибирск, ИГ и Г СО АН СССР, 1987.
  104. Механика контактных взаимодействий. (Под редакцией Егоровича И. И., Александрова В.М.) -М.: Физматлит, 2001,-672с.
  105. С.Г., Смолицкий X.JI. Приближенные методы решения дифференциальных и интегральных уравнений.- М.: Наука, 1965.
  106. JI.A. Матричный метод в теории распространения волн в слоистых упругих и жидких средах. М.: Наука. 1984.-202с.
  107. B.C. Задачи теории упругости для неоднородных сред. Сообщения по прикладной математике. М.- ВЦАН СССР, вып.4, 1976.
  108. В.Н. и др. Механика насыщенных пористых сред. М.: Наука, 1970.-336 с.
  109. В. Теория упругости. М.: Мир, 1975. — 872 с.
  110. ИЗ. Новацкий В. Динамика сооружений. М.: Госстройиздат, 1968.-376 с.
  111. П.М. Изгиб, устойчивость и колебания пластинок. Изд-во МГУ, 1958.
  112. Дж. Конечные элементы в нелинейной механике сплошных сред.- М.: Мир, 1976. 464 с.
  113. B.C. Особенности расчета многослойных дорожных одежд на изгиб//Автомобильные дороги, № 1, 1991. 11−12с.
  114. Г. И. Распространение волн в анизатропных упругих средах.-JI.: Наука, 1980, 175 с.
  115. Г. И., Молотков JI.A., Крауклис П. В. Волны в слоисто-однородных изотропных упругих средах. Ленинград: Наука, 1982.-289 с.
  116. Г. Я. О методе ортогональных многочленов в контактныхзадачах теории упругости .//ПММ, т. ЗЗ, в. З, 1969. С. 518−531.
  117. Г. Я. Контактные задачи для линейно-деформируемого основания.-Киев-Одесса: Вища школа, 1982,-168 с.
  118. Г. Я. Концентрация упругих напряжений возле штампов, разрезов, тонких включений и подкреплений.-М.: Наука, 1982,-344с.
  119. Г. Я. К решению задач механики и математической физи-Ф ки для слоистых сред. // Изв. АН СССР. Механика. -1978. -Т.31.2.
  120. В.Б. Методы динамической теории упругости. -М.: Наука, 1986.-328 с.
  121. И.И. Введение в теорию функций комплексного переменного. М.: Наука, 1967.
  122. А.К. Пространственная деформация многослойно-Ф го основания.//В сб. Устойчивость и прочность элементов конструкции, Днепропетровск, 1973.
  123. Проектирование нежестких дорожных одежд. ОДН 218.046 -01//Информавтодор, 2001.- 144с.
  124. B.C., Николаев А. Г. Решение пространственных задач теории упругости с помощью формул переразложения // ПМ. -1986.-22. -№ 7.
  125. О.Д. Нестационарные колебания упругой балки на вяз-# коупругом основании.// Изв. АН СССР, МТТ, № 1, 1992. С. 164 169.
  126. О.Д., Фрейгейт М. Р. О динамических свойствах системы: массивное тело-полуограниченная среда. // Докл. РАН. 1998, Т. 358, N 1, С.48−50.
  127. О.Д., Фрейгейт М. Р. О связи решений нестационарныхконтактных задач с резонансными свойствами исследуемых систем. // Докл. РАН. 1998, Т. 360, N 3, С.346−348.
  128. .С. Поведение дорожной конструкции как слоистой вязкоупругой среды под действием подвижной нагрузки.//Известия вузов. Строительство и архитектура, 1975. С. 78−83.
  129. Развитие контактных задач в СССР./ Под ред. Галина JI.A./ М.: Наука, 1976.-493 с.
  130. Развитие вибрационных исследований Земной коры в Сибири.// Новосибирск, ИГ и Г СО АН СССР, 1989.
  131. P.M. Задача Буссинеска для слоистого упругого по-лупространства.//Труды Ленинградского политехнического института. Ленинград, 1948. С. 3−18.
  132. Х.А., Демьянов Ю. А. Прочность при интенсивных кратковременных нагрузках. М.: Физматгиз, 1961. -400 с.
  133. В.Л., Проценко B.C. Контактные задачи теории упругости для неклассических областей. Киев: Наукова Думка. 1977. -235 с.
  134. В.Г. Руководство к решению задач теории упругости. -М.: Высшая школа. 1977. — 215 с.
  135. Дж. Теория звука. М.: Гостехиздат, т.2, 1955. — 476 с.
  136. Т.Г., Селезнев М. Г., Селезнева Т. Н. Пространственная задача об установившихся колебаниях упругого полупространства со сферической полостью.//ПММ, т.50, в.4, 1986. С. 651−656.
  137. М.Ю. О дисперсии волн в бесконечной упругой трехслойной пластине.//Изв. РАН Мех. тверд, тела-1998, N1, С. 166−172.
  138. А.О. Простейшие методы учета многослойности дорожных одежд при оценке их напряженно деформируемого состоя-ния.//Труды Союздорнии. Исследования по механике дорожныходежд.-М.- 1985. С. 4−13.
  139. П. Ортогональные многочлены. М.: Физматгиз, 1962. -500 с.
  140. В.М. Динамические контактные задачи. Киев: Наукова думка, 1970,-283 с.
  141. В.М., Ермоленко Н. П., Зайцева Е. А. Неосесимметрич-ные периодические и нестационарные колебания круглого штампа на упругом полупространстве.// Прикладная механика. (Киев),-1997−33, N5, С. 41−48.
  142. В.М., Трофимчук А. Н., Савицкий О. А. Динамические контактные задачи для слоистых сред.// Изв. вузов Сев.-Кав. регион. Естеств. н.- 2001, Спец. вып., С. 138−140.
  143. В.М., Трофимчук А. Н., Савицкий О. А. Колебания и волны в слоистых средах. Киев: Наукова думка, 1990. — 224 с.
  144. М.Г. Возбуждение волн в двухслойной среде колеблющимся штампом.//ПММ, т.39, в.2, 1975. С. 381−384.
  145. Л.И. Нестационарные упругие волны. Л.: Судостроение, 1972.-371 с.
  146. Л.И., Яковлев Ю. С. Интегральные преобразования в нестационарных задачах механики. -Л.Судостроение, 1980.-343 с.
  147. А.В., Динамика дорожных одежд автомобильных дорог. Омск, Запсибиздат, 1975.
  148. А.В., Малофеев А. Г. Э5кспериментальное исследование волн колебаний дорожных покрытий при движении автомобилей.// Прикладная механика. T. IX, в.1, 1973.
  149. Собисевич J1.E., Шумейко В. И., Селезнев М. Г., Ляпин А. А., Со-бисевич А. Л, Корабельников Г. Я. Локальные резонансы в слоистых средах.// М.: ОИФЗ РАН, Московский филиал ГНИЦ ПГК при КубГУ Министерства образования РФ, 2000. 178 с.
  150. Д.В., Федоров С. Н. Нестационарные колебания упругой анизатропной полуплоскости.// Материалы 4-го Между-нар. симп. «Динам, и технол. пробл. мех. конструкций и сплош. сред.», Ярополец, 16−20 фев. 1998.- М.-С 24.
  151. С.П. Пластинки и оболочки. М.: Гостехиздат, 1948.
  152. Е.С. Введение в теорию интегралов Фурье. М.: Гостехиздат, 1948. — 479 с.
  153. Е. Разложения по собственным функциям, связанные с дифференциальными уравнениями второго порядка. М.: Иноиз-дат, 1960.
  154. А.Н., Самарский А. А. Уравнения математической физики,— М.: Наука, 1972.
  155. К. Дж. Интегральные преобразования в математической физике. М.: Гостехиздат, 1956. — 204 с.
  156. А.Ф. Метод собственных векторных функций в пространственных задачах теории упругости .//Киев: Наукова думка, 1979.-261 с.
  157. Я.С. Интегральные преобразования в задачах теории упругости.-Ленинград: Наука, 1967.
  158. М.В. Метод перевала. М.: Наука, 1977.
  159. А.П. Колебания деформируемых систем. М.: Машиностроение, 1970. — 736 с.
  160. А.Р. Контактная задача для упругой кусочно-однородной полуплоскости с конечным стрингером. Межвузовский сб. научных трудов, Механика, Ереван: изд. ЕГУ, 1990, N8, С.97−106.
  161. Л.П., Довгалюк А. В. Колебания полуограниченых и ограниченных слоистых тел стохастической структуры.// Прикладная механика (Киев)-1997−33, N5, С.13−19.
  162. А.И. Прикладные методы решения краевых задач строительной механики. М.: Стройиздат, 1984. — 334 с.
  163. И.С. Вибрационное излучение сейсмических волн.-М.: Недра, 1984.-220 с.
  164. Г. С., Никишин B.C. Пространственные задачи теории упругости для многослойных сред. М.: ВЦАН СССР, 1970.
  165. Ф., Гупта Г. Задача о полуплоскости с упругой накладкой. ПМ. Труды Амер. общ. инж.-мех., сер. Е, 1971, Т.38, N4,1. С.323−327.
  166. Э., Эмде Ф., Леш Ф. Специальные функции (формулы, графики, таблицы) М.: Наука, 1964.
  167. Abdelkarim А.М.А.М., Vrouwenvelder A.C.W.M., Verweij M.D. Analysis of the dynamic response of layered, elastic media by means of the fast Fourier transform.// Heron.-1999−44, N2, p. l09−125.
  168. Alblas J.B., Kuypers W.I.I. On the Diffusion of Load from a Stiff-ener into an Infinite Wedge Shaped Plate.// Applied Scienntific Research, Series A, Vol.15, N6, 1965−1966, p.429.
  169. Ivanov Ts., Savova R. A note on Rayleigh waves.// J.Teor. and Appl. Mech.-1996−26, N3, p.57−61.
  170. Jones D.V., Le Houedec D., Peplow A.T., Petyt M. Around vibration in the vicinity of a moving harmonic rectangular load on a half-space.// Eur. J. Mech.A.-1998−17, N1, p. 153−166.
  171. Muravskii G. On time-harmonic problem for non-homogeneous elastic half-space with shear modulus limited at infinite depth.// Eur. J. Mech. A.-1997−16, N2, p.277−294.
  172. Nath Sisir, Sengupta P.R. Steady-state response to moving loads in an elastic solid media.// Indian J. Pure and Appl. Math.-1999−30, N3, p.317−327.
  173. Zhang Bixing, Yu M., Lan C.Q., Xiong Wei. Elastic wave and excitation mechanism of surface waves in multilayered media.// J. Acoust. Soc. Amer.-1996−100, N6, p.3527−3538.
  174. Описание программы расчета влияния скорости движения нагрузки на распространение внутренних волн. Плоский случай.
  175. Программа PLZ предназначена для построения диаграмм направленности поля внутренних волн в однородной полуограниченной среде под воздействием движущейся нагрузки в случае плоской постановки задачи.
  176. R —расстояние от источника колебаний до рассматриваемой точки в среде-
  177. Выходными данными программы PLZ являются диаграммы направленности, посчитанные для различных значений скорости движения нагрузки и представленные в п. 3.3.3 (рис. 3.41).vO := 50 со := 20 a := 2 tO := 8000 R:=100 400 i uO т
  178. J := 1 2 v (j) :=- p := 1800 pO := v (2) -p и := —— u = 4.5×103j 2-t0m := 0.58 ф (ш) := + m— k (j) :=60 60 v (j) v (j)
  179. A (j, m) := (k (j)2 l)-(k (j)2-sinU (m))2 — l) B (j, m) := 6(j)-k (j) (l — k (j)2-sin (4.(m))2)
  180. Up (m) := -Jul (m)2 + wl (m)2w2(m) := exp^-R-S (2,m) i-^fw2(m)kr V|S1 l (2,m)|1. Us (m) := >/u2(m)2 + w2(m)aUp (m) := >/lm (Up (m))2 + Re (Up (m))2aUs (m) := >/lm (Us (m))2 + Re (Us (m))26.104.5−10aUs (m) 3−101.5−102.101.5−10aUp (m) 1-Ю5.10
  181. Описание программы расчета влияния скорости движения нагрузки на распространение внутренних волн в однородном полупространстве. Пространственный случай.
  182. Программа PRZ предназначена для построения диаграммнаправленности поля внутренних волн в однородномполупространстве в дальней от области нагружения зоне подвоздействием движущейся нагрузки в пространственном случае.
  183. R — расстояние от источника колебаний до рассматриваемой точки в среде-
  184. S (j) := — A (j, m) := ?k (j))2 -T. Ck (j))2.p (m) -Цv (j)
  185. REAL*8 у, x, aa, h, c, dd, 11, ml, pi, 10, m0,p0,lm, w, wO, hw, wn, tetl, tet2, yi, xO, yO, * zO, P, GLUB, SHAG, ABSERR, RELERR, EPS, AMAXDL, DLINA, ABSERR1, RELERR1, EPS1, AMAXDL1 COMPLEX*1 б B, BX, BB, DKONTUR1, ci, XXl
  186. SHAG=1.2 ABSERR=1.E-4 ABSERR1=1.E-3 RELERR=1.E-4 RELERR1=1.E-3 EPS=1.E-3 EPS1=1.E-2 AMAXDL=1.E+3 AMAXDL1=1.E+2 DLINA=1.l*tet2+0.5
  187. XX1=DK0NTUR1(FXX1,GLUB, DLINA, SHAG, ABSERR1, RELERR1,EPS1,AMAXDL1) write (*,*) w, XXI1 continuecall gettim (il, i2, i3,i4) write (*,*) il, i2, i3,i4 write (1,*) il, i2, i3,i4 stop end
  188. Заполнение матрицы в квадрате (1−18,1−18)--------------do 2 k=0,15,3do 3 i=l, 3 do 4 j=l, 3 .4 В (k+i, k+j)=c (i, j) 3 continue 2 continue
  189. Заполнение матрицы в квадрате (19−36,19−36)----------------do 8 k=18,33,3do 9 i=l, 3 do 10 j = l, 3
  190. B (i, 19+(k-1)*3)=DKONTUR (FLl, GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL) В (i, 20+(k-1)*3)=DKONTUR (FL2,GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL) В (i, 21+(k-1)*3)=DKONTUR (FL3,GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL) 18 continuedo 20 i=10,12 yi=y (i-9)
  191. В (i, 19+(k-1)*3)=DKONTUR (FLl, GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL) В (i, 20+(k-1)*3)=DKONTUR (FL2,GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL) В (i, 21+(k-1)*3)=DKONTUR (FL3,GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL)21 continuedo 22 i=13,15 yi=y (i-12)
  192. B (i, 19+(k-1)*3)=DKONTUR (FL4,GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL)
  193. В (i, 20+(к-1)*3)=DKONTUR (FL5,GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL) В (i, 21+(k-1)*3)=DKONTUR (FL6,GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL)22 continue flag=3do 23 ±=1,9 yi=y{i-6)
  194. B (i, 19+(k-1)*3)=DKONTUR (FLl, GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL) B (i, 20+(k-1)*3)=DKONTUR (FL2,GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL) B (i, 21+(k-1)*3)=DKONTUR (FL3,GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL)23 continuedo 24 i=16,18 yi=y (i-15)
  195. BX (i)=BX (i)+Bl (i, j)*BB (j) 51 continuewrite (1,*) BX (i) 50 continue
  196. FXX1=(BX (1)/sqrt (aa**2-yO**2)+BX (2)*yO+BX (3)*y0**2) *exp (alf*z0)returnendfunction FLl (xl) complex PARA, PARB complex*16 FL1, SXT, gam, ci real*8 xl, h, aa common /SI/ h, aa, ci common /S12/ PARA, PARB external SXTgam=PARA*x1+PARB
  197. FLl=ci/gam*(1-exp (ci*gam*h))*SXT (gam)returnendfunction FL2(x2) complex PARA, PARB complex*16 FL2, SXT, gam, ci real*8×2,h, aa common /SI/ h, aa, ci common /S12/ PARA, PARB external SXTgam=PARA*x2+PARB
  198. FL2=(exp (ci*gam*h)-ci*h*gam-l)/gam**2*SXT (gam)returnendfunction FL3(x3) complex PARA, PARB complex*16 FL3, SXT, gam, ci real*8×3,h, aa common /SI/ h, aa, ci common /S12/ PARA, PARB external SXTgam=PARA*x3+PARB
  199. FL3=2*ci/gam**3*(exp (ci*gam*h)+h*h*gam**2/2+h*gam/ci-1)* *SXT (gam) return endfunction FL4(x4) complex PARA, PARB complex*16 FL4, SXT, gam, ci real*8×4,h, aa common /SI/ h, aa, ci common /S12/ PARA, PARB external SXTgam=PARA*x4+PARB
  200. FL4=ci/gam*(exp (-ci*gam*h)-1)*SXT (gam)returnendfunction FL5(x5) complex PARA, PARB complex*16 FL5, SXT, gam, ci real*8×5,h, aa common /SI/ h, aa, ci common /S12/ PARA, PARB external SXTgam=PARA*x5+PARB
  201. FL5=(1-(ci*h*gam+l)/exp (ci*gam*h))/gam**2*SXT (gam)returnendfunction FL6(x6) complex PARA, PARB complex*16 FL6, SXT, gam, ci real*8×6,h, aa common /SI/ h, aa, ci common /S12/ PARA, PARB external SXTgam=PARA*x6+PARB
  202. FMl=pi*CJN (uf, 0)*exp (ci*bet*aa)*SYT (bet, 1) returnendfunction FM2(y2) complex PARA, PARB complex*16 FM2, SYT, bet, ci real*8 y2, h, aa common /SI/ h, aa, ci common /S12/ PARA, PARB external SYTbet=PARA*y2+PARB
  203. FM2=(exp (2*ci*bet*aa)*(l/bet-aa*ci)-aa*ci-l/bet)/bet* *SYT (bet, 1) return endfunction FM3(у3> complex PARA, PARB complex*16 FM3, SYT, bet, ci real*8 y3, h, aa common /SI/ h, aa, ci common /S12/ PARA, PARB external SYTbet=PARA*y3+PARB
  204. FM4=pi*CJN (uf, 0)/exp (ci*bet*aa)*SYT (bet, 1) returnendfunction FM5(y5) complex PARA, PARB complex*16 FM5, SYT, bet, ci real*8 y5, h, aa common /SI/ h, aa, ci common /S12/ PARA, PARB external SYTbet=PARA*y5+PARB
  205. FM5=(l/bet-aa*ci-(aa*ci+l/bet)/exp (2*ci*bet*aa))/bet* *SYT (bet, 1) return endfunction FM6(y6) complex PARA, PARB complex*16 FM6, SYT, bet, cireal*8 y6, h, aa common /SI/ h, aa, ci common /S12/ PARA, PARB external SYTbe t = PARA* у 6 + PARB
  206. FXl=pi*CJN (uf, 0)*SYT (bet, 0) returnendfunction FX2(y2)complex PARA, PARBcomplex*16 FX2, SYT, bet, ci, EAreal*8 y2, h, aacommon /SI/ h, aa, cicommon /S12/ PARA, PARBexternal SYTbet=PARA*y2+PARB1. EA=exp (ci*bet*aa)
  207. С DESCRIPTION OF PARAMETERS
  208. С A INPUT MATRIX, DESTROYED IN COMPUTATION AND REPLACED BY1. С RESULTANT INVERSE.1. С N ORDER OF MATRIX A1. С D RESULTANT DETERMINANT
  209. С L WORK VECTOR OF LENGTH N
  210. С M WORK VECTOR OF LENGTH NС1. С REMARKS
  211. С MATRIX A MUST BE A GENERAL MATRIXС
  212. С SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED1. С NONEС1. С METHOD
  213. С THE STANDARD GAUSS-JORDAN METHOD IS USED. THE DETERMINANT
  214. С IS ALSO CALCULATED. A DETERMINANT OF ZERO INDICATES THAT1. С THE MATRIX IS SINGULAR.1. С .
  215. SUBROUTINE MINV (A, D, L, M, N) DIMENSION A (l), L (1), M (1) complex*16 A, D, BIGA, HOLD1. С .
  216. С IF A DOUBLE PRECISION VERSION OF THIS ROUTINE IS DESIRED, THE
  217. С С IN COLUMN 1 SHOULD BE REMOVED FROM THE DOUBLE PRECISION
  218. С STATEMENT WHICH FOLLOWS. С
  219. С DOUBLE PRECISION A, D, BIGA, HOLD, DABSС
  220. С THE С MUST ALSO BE REMOVED FROM DOUBLE PRECISION STATEMENTS
  221. С APPEARING IN OTHER ROUTINES USED IN CONJUNCTION WITH THIS1. С ROUTINE. С
  222. С THE DOUBLE PRECISION VERSION OF THIS SUBROUTINE MUST ALSO
  223. С CONTAIN DOUBLE PRECISION FORTRAN FUNCTIONS. ABS IN STATEMENT
  224. С 10 MUST BE CHANGED TO DABS.1. С .
  225. С SEARCH FOR LARGEST ELEMENTС1. D=1.0 NK=-N
  226. DO 80 K=1,N NK=NK+N L (K)=K M (K)=K KK=NK+K BIGA=A (KK) DO 20 J=K, N IZ=N*(J-l) DO 20 I=K, N IJ=IZ+I
  227. IF (cABS (BIGA) — cABS (A (IJ))) 15,20,20 15 BIGA=A (IJ) L (K)=1 M (K)=J 20 CONTINUE
  228. С---------- INTERCHANGE ROWS----------1. J=L (K)1.(J-K) 35,35,25 25 KI=K-N
  229. DO 30 1=1,N KI=KI+N HOLD=-A (KI) JI=KI-K+J A (KI)=A (JI) 30 A (JI) =HOLD
  230. С--------- INTERCHANGE COLUMNS---------35 I=M (K)1.(I-K) 45,45,38 38 JP=N*(1−1) DO 4 0 J=1,N JK=NK+J JI=JP+J HOLD=-A (JK) A {JK)=A (JI) 40 A{JI) =HOLD
  231. С----- DIVIDE COLUMN BY MINUS PIVOT (VALUE OF PIVOT ELEMENT IS CONTAINED IN BIGA)45 ABIGA=cabs (BIGA)1.(ABIGA) 48,46,48 4 6 D=0.01. RETURN 48 DO 55 1=1,N1.(I-K) 50,55,50 50 IK=NK+I
  232. A (IK)=A (IK)/(-BIGA) 55 CONTINUE
  233. C---------- REDUCE MATRIX-------------
  234. DO 65 1=1,N IK=NK+I HOLD=A (IK) IJ=I-N DO 65 J=1,N IJ=IJ+N1.(I-K) 60,65, 60 60 IF (J-K) 62,65,62 62 KJ=IJ-I+K
  235. A (IJ)=HOLD*A (KJ)+A (IJ) 65 CONTINUE
  236. C----------DIVIDE ROW BY PIVOT-----------
  237. KJ=K-N DO 75 J=1,N KJ=KJ+N1.(J-K) 70,75,70 70 A (KJ)=A (KJ)/BIGA 75 CONTINUE
  238. С---------- PRODUCT OF PIVOTS--------------------1. D=D*BIGA
  239. С--------- REPLACE PIVOT BY RECIPROCAL-----------
  240. A (KK)=1.0/BIGA 8 0 CONTINUE
  241. С---------- FINAL ROW AND COLUMN INTERCHANGE------1. K=N 100 K=(K-1)1.(K) 150,150,105 105 I=L (K)1.(I-K) 120,120,108 108 JQ=N*(K-l) JR=N*(1−1) DO 110 J=1,N JK=JQ+J HOLD=A (JK) JI=JR+J A (JK)=-A (JI) 110 A (JI) =HOLD 120 J=M (K)1.(J-K) 100,100,125 125 KI=K-N
  242. DO 130 1=1,N KI=KI+N HOLD=A (KI) JI=KI-K+J A (KI)=-A (JI) 130 A (JI) =HOLD1. GO TO 100 150 RETURN ENDс------------------------------------------------------------------
  243. FUNCTION DKONTUR (F, GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL)
  244. С ПРОГРАММА ИНТЕГРИРОВАНИЯ ПО КОРЫТУ
  245. С F- функция типа COMPLEX*16 от аргумента типа REAL*8
  246. С * в ней необходим COMMON блок /S12/PARA, PARB1. С GLUB- глубина корыта1. С DLINA- его длина
  247. С SHAG- шаг разбиения контура вне корыта
  248. С ABSERR, RELERR- абс. и отн. погрешности вычисления интегралов разб.
  249. С EPS- отн. погрешность оценки хвоста
  250. С AMAXDL- максимальная длина контура по оси
  251. COMPLEX*16 DKONTUR, F, RESULT COMPLEX CI, PARA, PARB
  252. REAL*8 A, B, ABSERR, RELERR, ERREST, FLAG, AK, AR, EPS, AMAXDL, GLUB, DLINA, SHAG
  253. COMMON /S12/PARA, PARB EXTERNAL F1. CI=(0., 1.)
  254. PARA=(0.,-1.) PARB=(0., 0.) A=0.DO B=GLUB
  255. CALL CDQUAN8(F, A, B, ABSERR, RELERR, RESULT, ERREST, NOFUN, FLAG)1. DKONTUR=-CI*RESULT1. PARA=(1., 0.)1. PARB=-CMPLX (0., GLUB)1. A=0.DO1. B=DLINA
  256. CALL CDQUAN8(F, A, B, ABSERR, RELERR, RESULT, ERREST, NOFUN, FLAG)
  257. DKONTUR=DKONTUR+RESULT PARA=(0., 1.) PARB=CMPLX (DLINA, 0.) A=-GLUB B=0.DO
  258. CALL CDQUAN8(F, A, B, ABSERR, RELERR, RESULT, ERREST, NOFUN, FLAG)1. DKONTUR=DKONTUR+RESULT*CI1. PARA=(1., 0.)1. PARB=(0., 0.)1. A=DLINA1 B=A+SHAG
  259. CALL CDQUAN8{F, A, B, ABSERR, RELERR, RESULT, ERREST, NOFUN, FLAG)1. DKONTUR=DKONTUR+RESULT1. AK=CDABS (DKONTUR)1. AR=CDABS (RESULT)1.(AR/AK.LE.DBLE (EPS)) RETURN1.(B.GE.AMAXDL) GO TO 21. A=B1. GO TO 12 PRINT 50, AMAXDL
  260. FORMAT (' Выход no MAX длине=', E15.5) RETURN END
  261. FUNCTION DKONTUR1(F, GLUB, DLINA, SHAG, ABSERR, RELERR, EPS, AMAXDL)
  262. С ПРОГРАММА ИНТЕГРИРОВАНИЯ ПО КОРЫТУ
  263. С F- функция типа COMPLEX*16 от аргумента типа REAL*8
  264. С * в ней необходим COMMON блок /S12/PARA1,PARB11. С GLUB- глубина корыта1. С DLINA- его длина
  265. С SHAG- шаг разбиения контура вне корыта
  266. С ABSERR, RELERR- абс. и отн. погрешности вычисления интегралов разб.
  267. С EPS- отн. погрешность оценки хвоста
  268. С AMAXDL- максимальная длина контура по оси
  269. COMPLEX*16 DKONTUR1, F, RESULT1 COMPLEX CI1, PARA1,PARB1
  270. REAL*8 Al, Bl, ABSERR, RELERR, ERREST, FLAG, AK, AR, EPS, AMAXDL, GLUB, SHAG
  271. COMMON /S12/PARA1,PARB1 EXTERNAL F
  272. CI1=(0., 1.) PARA1=(0.,-1.) PARB1=(0., 0.) A1=0.DO B1=GLUB
  273. CALL CDQUAN81(F, Al, Bl, ABSERR, RELERR, RESULT1, ERREST, NOFUN, FLAG)1. DKONTURl=-CI1 *RESULT11. PARA1=(1., 0.)1. PARB1=-CMPLX (0., GLUB)1. A1=0.DO1. B1=DLINA
  274. CALL CDQUAN81(F, Al, Bl, ABSERR, RELERR, RESULTl, ERREST, NOFUN, FLAG)1. DKONTURl=DKONTURl+RESULT11. PARAl=(0., 1.)1. PARB1=CMPLX (DLINA, 0.)1. Al=-GLUB1. B1=0.DO
  275. CALL CDQUAN81(F, Al, Bl, ABSERR, RELERR, RESULTl, ERREST, NOFUN, FLAG) DKONTURl=DKONTURl+RESULTl*CI1 PARA1=(1., 0.) PARB1=(0., 0.) A1=DLINA 1 B1=A1+SHAG
  276. FUNCTION CJN (U, N) С *** Функция Бесселя JN (U) одинарной точности INTEGER N, К, M, NN4 COMPLEX*8 CJN, U, Y, Т, P, Q
  277. REAL*4 R, PI4/0.7 853 982/, SQ4/2.256 758/, AY,
  278. GAN (0:2)/5.72,5.8,5.9/, Е/0.1 192 093Е-0б/ С *** Блок 1. Формула Ганкеля?1.(ABS (U).GT.GAN (N)) GO TO 30 С *** Блок 2. Сумма ряда Q=1.01.(N.EQ.O) GO TO 22 Y=0.5*U DO 21 K=l, N21 Q=Q*Y/K
  279. Y=-0.25*U*U T=Q*Y/(N+l) CJN=Q+T1. K=1
  280. K=K+1 M=K*(N+K) T=T*Y/M CJN=CJN+T1.(ABS (T).GT.E) GO TO 23 С *** Выход из CJN RETURN
  281. С *** Блок 3. Формула Ганкеля
  282. M=N+N NN4=M*M Y=0.125/U AY=ABS (Y) P=1.01. Q=Y*(NN4−1)1. T=Q1. M=11. K=11. С * * * Сумма для «P»
  283. M=M+1 K=K+2 R=K*K-NN4 R=R/M1.(AY*R.GT.1.0) GO TO 321. T=T*Y*R1. P=P+T1.(ABS (T).LE.E) GO TO 32 С *** Сумма для «Q» M=M+1 K=K+2 R=NN4-K*K R=R/M1.(AY*R.LT.-1.0) GO TO 321. T=T*Y*R1. Q=Q+T1.(ABS (T).GT.E) GO TO 31 32 T=U-(N+N+l)*PI4
  284. CJN=SQ4 *SQRT (Y)*(P*COS (T)-Q*SIN (T)) С *** Выход из CJN RETURN END
  285. FUNCTION RIN (X, N) С *** Функция IN (X) одинарной точности INTEGER N, К, M
  286. REAL*4 RIN, X, Y, T, E/0.1 192 093Е-06/ С *** Сумма ряда Y=0.25*X*X T=Y/(N+l) RIN=1.0+T K=1
  287. K=K+1 M=K*(N+K) T=T*Y/M RIN=RIN+T1.(T.GT.E) GO TO 11 IF (N.EQ.O) GO TO 13 Y=0.5*X DO 12 K=1,N
  288. RIN=RIN*Y/K C *** Выход из RIN13 RETURN END
  289. SUBROUTINE CDQUAN8(FUN, A, B, ABSERR, RELERR, RESULT, ERREST, NOFUN, FLAG)
  290. COMPLEX*16 FUN, RESULT, F0, F (16), CORll
  291. COMPLEX*1 6 QPREV, QNOW, QDIFF, QLEFT, QRIGHT (31), FSAVE (8,30)
  292. REAL*8 A, B, ABSERR, RELERR, ERREST, FLAG, ESTERR, TOLERR,
  293. REAL*8 WO, Wl, W2, W3,W4,AREA, XO, STONE, STEP, TEMP, X (16), XSAVE (8, 30)
  294. TEGER LEVMIN, LEVMAX, LEVOUT, NOMAX, NOFIN, LEV, NIM, I, J, NIMA, NOFUN
  295. VMIN=1 LEVMAX=30 LEVOUT=6 NOMAX=5000
  296. NOFIN=NOMAX-8 *(LEVMAX-LEVOUT+2**(LEVOUT+1))1. W0=3956.0/14 175.01. Wl=23 552.0/14 175.01. W2=-3712.0/14 175.01. W3=41 984.0/14 175.01. W4=-18 160.0/14 175.01. FLAG=0.0
  297. RESULT=(0.0,0.0) COR11=0.0 ERREST=0.0 AREA=0.0 NOFUN=01.(A.EQ.B) RETURN1. V=01. NIM=11. X0=A1. X (16)=B1. QPREV=0.01. F0=FUN (X0)1. TONE=(B-A)/16.01. X (8)=(X0+X (16))/2.01. X (4)=(X0+X (8))/2.01. X (12)=(X (8)+X (16))/2.01. X (2)=(X0+X (4))/2.0
  298. X (6) = (X (4)+X (8)) /2. 025
Заполнить форму текущей работой

На отлично!