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

Геометрия, конструирование и исследование напряженно-деформированного состояния оболочек в форме каналовых поверхностей Иоахимсталя

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

Целью диссертационной работы является исследование геометрии каналовых поверхностей Иоахимсталя, конструирование оболочек различных очертании на основе этих поверхностей, разработка и реализация на ЭВМ метода расчета оболочек сложной геометрии и, в частности, для каналовых поверхностей Иоахимсталя. Воздействие окружающей среды играет особую роль при выборе оптимальной формы оболочки в авиа, судо… Читать ещё >

Содержание

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

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

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

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

1. Архитектурная выразительность — при покрытии общественных зданий с учетом национальной специфики.

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

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

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

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

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

Все вышесказанное подтверждает актуальность темы диссертации.

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

Научная новизна работы заключается в следующем:

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

Предложено три способа образования каналовых поверхностей Иоахимсталя, на основании которых выполнены три макета.

Разработаны новые конструктивные и архитектурные формы на основе каналовых поверхностей Иоахимсталя.

Разработан и реализован на ЭВМ алгоритм применения вариационно-разностного метода для расчета пластин и оболочек сложной геометрии. Алгоритм и программа апробированы на тестовых задачах.

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

Научная и практическая ценность работы:

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

Внедрение работы.

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

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ.

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

2. Получено условие, при выполнении которого рассматриваемая поверхность является каналовой поверхностью Иоахимсталя.

3. Предложено три способа образования каналовых поверхностей Иоахимсталя.

4. Разработаны новые конструктивные и архитектурные формы на основе каналовых поверхностей Иоахимсталя. Выполнены макеты, реализующие три способа образования таких поверхностей.

5. Разработан алгоритм расчета пластин и оболочек сложной геометрии вариационно-разностным методом.

6. По описанному алгоритму ВРМ реализована программа расчета тонкостенных конструкций на ЭВМ. Алгоритм и программа апробированы на тестовых задачах.

7. Проведены расчеты оболочек в форме каналовых поверхностей Иоахимсталя вариационно-разностным методом на различные виды нагрузок.

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

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

10. Результаты исследований внедрены в более чем 10 дипломных работах студентов.

ЗАКЛЮЧЕНИЕ

.

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

По итогам работы можно сформулировать следующие основные результаты и выводы:

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

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

  1. Н.П., Андреев Н. П., Дерюга А. П. Вариационные принципы теории упругости и теории оболочек. М.: Наука, 1978. — 288 с.
  2. A.A., Дубинский Ю. А., Копченова Н. В. Вычислительные методы для инженеров: Учеб. пособие. -М.: Высш. шк., 1994. 544 с. ил.
  3. В.Н., Хампе Э., Рауэ Э. Проектирование железобетонных тонкостенных пространственных конструкций. М.: Стройиздат, 1990. -232 с.
  4. Н.С. Численные методы. М.: Наука, 1975. — 631 с.
  5. И.К. Геометрия циклид Дюпена и их применение в строительных объектах // Расчет оболочек строительных конструкций. М.: УДН, 1982. -С.161- 129.
  6. И.Н., Семендяев К. А. Справочник по математике для инженеров и учащихся ВТУЗОВ. М.: Наука, 1986. — 544 с. 7. Вайнберг Д. В., Вайнберг Е. Д. Расчет пластин. К.: Буд1вельник, 1970.-436 с.
  7. П.М., Варвак Л. П. Метод сеток в расчетах строительных конструкций. М.: Строийиздат, 1977. — 160 с.
  8. К. Вариационные методы в теории упругости и пластичности: Пер. с англ. М.: Мир, 1987. — 544 с.
  9. В.З. Общая теория оболочек и ее приложение в технике. М.: Гостехиздат, 1949. — 784 с.
  10. В.З. Строительная механика оболочек. М.: Стройиздат, 1936.- 263 с.
  11. В.З. Тонкостенные упругие системы. М.: Госстройиз-дат, 1958.-502 с.
  12. В.З., Леонтьев H.H. Балки, плиты и оболочки на упругом основании. — М., i960.
  13. М. Я. Дифференциальная геометрия. М., Л.: ГИТТЛ, 1949.-512 с.
  14. К.З., Паймушин В. Н. Теория оболочек сложной геометрии (геометрические вопросы теории оболочек). Казань: Изд-во Казанского университета, 1985. — 164 с.
  15. Д., Кон.-Фоссен С. Наглядная геометрия. М.: Наука, 1981. -344 с.
  16. Дж., Ван Лоун Ч. Матричные вычисления: Пер. с англ.- М.: Мир, 1999.-548 е., ил.
  17. А.Л. Качественное исследование напряженного состояния тонкой оболочки // ПММ. 1945. — т.9. — вып. 6. — с. 14−20.
  18. А.Л. Теория упругих тонких оболочек. М.: Гос-техиз-дат, 1953. — 512 с.
  19. А.Н. Гидроэнергетические оборудование гидроэлектростанции и его монтаж. М.: Энергия, 1972. — 320 е., ил.
  20. С.А. и др. Конструкции гидротурбин и расчеты их деталей. -М.: Машгиз, 1956. 178 с.
  21. Я.М. Решение задач теории оболочек методами численного анализа // Прикладная механика. 1984. — т.20. — № 10. — с.32−40.
  22. Я.М., Гуляев В. И., Гоцуляк Е. А., Ашури К. Напряженно- деформированное состояние трубчатых оболочек под действием равномерно распределенного давления // Прикладная математика. 1983. -т.Х1Х. -№ 3. — с. 11−18.
  23. Я.М., Мукоед А. П. Решение задач теории оболочек на ЭВМ. Киев: Высш. шк., 1979. — 280с.
  24. М.Ф. Отсасывающие трубы гидроэлектростанций. -М.: Энергия, 1970.-270 с.
  25. В.И., Баженов В. А., Гоцуляк Е. А., Гайдайчук В. В. Расчет оболочек сложной формы. Киев: Буд1вельник, Библиотека проектировщика, 1990.- 192 с.
  26. А.П. Двойственность вариационно-разностных схем расчета оболочек // Пространственные конструкции в красноярском крае.- Красноярск, 1982. с. 19−32.
  27. Ю.А., Жуковский Э. З., Ермолов В. В. и др. Современные пространственные конструкции (железобетон, металл, дерево, пластмасс): Справочник. М: Высш. шк., 1991. — 543 с. ил.
  28. Жиль-Улбе Матье. Алгоритм расчета эпитрохоидальной оболочки по без моментной теории в перемешенях // Строительная механика инженерных конструкций и сооружений. Волгоград: Перемена, 1996. — Вып.7.- с.18−21.
  29. Жиль-Улбе Матье. Разрешающее дифференциальное уравнение эпитрохоидальной оболочки в перемещениях // Строительная механика инженерных конструкций и сооружений. М.: Изд-во РУДН, 1996. — Вып.6. -с.56−59.
  30. Жиль-Улбе Матье. Расчет эпитрохоидальной оболочки в усилиях и в перемешенях: Дис.. канд. техн. Наук. М.: РУДН, 1997. — 134 с.
  31. Э. З. Шилобреев Ю.А. Шевченко О. В. Большепролетные пространсвеные конструкции в мировой практике // Строительная механика и расчет конструкций. 1990. — № 5. — с.87−91.
  32. О. Метод конечных элементов в технике: Пер. с англ. -М.: Мир, 1975.-542 с.
  33. В.Н. Каналовые поверхности Иоахимсталя с плоской линией центров // Исследования пространственных систем: Материалы семинара кафедры сопротивления материалов РУДН. М.: Изд-во РУДН, 1996.- с.32−36.
  34. В.Н. Вариационные принципы и методы решения задач теории упругости: Учеб. пособие. М.: Изд-во РУДН, 2001. — 176 с. ил.
  35. В.Н. Геометрия циклических поверхностей // Сборник научных работ аспирантов инженерного факультета. М.: УДН, 1971. — вып.8. — с.137−142.
  36. В.Н. Конструирование оболочек на основе каналовых поверхностей Иоахимсталя // Вестник РУДН: специальный выпуск серии «Инженерные исследования». -М.: Изд-во РУДН, 2000. № 1. — с.57−61.
  37. В.Н. Некоторые аспекты геометрии циклид Дюпена // Вестник РУДН: Серия «Инженерные исследования», спец. выпуск «геометрия и расчет тонкостенных пространственных конструкций». М.: Изд-во РУДН, 2002. -№ 1.- с. 12−21.
  38. В.Н. Некоторые вопросы теории поверхностей с семейством плоских координатных линий // Расчет оболочек строительных конструкций. М.: УДН, 1977. — вып. 10. — с. 37−48.
  39. В.Н. Об одном классе каналовых поверхностей // Вопросы прочности пространственных систем: Материалы XXVIII научной конференции инженерного ф-та. М.: РУДН, 1992. — с. 54−62.
  40. В.Н. Об уравнениях безмоментной теории оболочек в форме каналовых поверхностей Иоахимсталя // Строительная механика инженерных конструкций и сооружений: Межвузовский сборник научных трудов. М.: МБК «Биоконтроль», 1994. — вып.4. — с. 83−85.
  41. В.Н. Расчет оболочек в форме циклических поверхностей: Дис.. канд. техн. наук. -М.: УДН, 1970. 117 с.
  42. В.Н. Теория расчета оболочек в форме циклических поверхностей // Доклады научно-технической конференции инженерного факультета. М.: УДН, 1971. — с. 27−29.
  43. В.Н. Условия образования каналовых поверхностей //
  44. Строительная механика инженерных конструкций и сооружений: Межвузовский сборник научных трудов. М.: МБК «Биоконтроль», 1995. — вып.5.- с.7−16.
  45. В.Н. Циклические поверхности: геометрия, классификация, конструирование оболочек // Труды международной научной конференции «Оболочки- 2001». М: Изд-во РУДН, 2001. — с. 127−134.
  46. В.Н., Махмуд Хуссейн Аль-Хадж. Алгоритм расчета эпи-трохоидальной оболочки по безмоментной теории // Вопросы прочности пространственных систем: Материалы XXVIII научной конференции инженерного ф-та. М.: РУДН, 1992. — с.58−63.
  47. В.Н., Махмуд Хуссейн Аль-Хадж. Координатная сеть линий кривизны эпитрохоидальной поверхности // Исследования по строительной механике пространственных систем. М.: УДН, 1990. — с.38−44.
  48. В.Н., Махмуд Хуссейн Аль-Хадж. Расчет эпитрохоидальной оболочки по безмоментной теории // Расчет и проектирование гражданских и гидротехнических сооружений: Межвузовский сборник научных трудов. М.: МБК «Биоконтроль», 1994. — вып.4. — с. 21−26.
  49. И.Е. Справочник по гидравлическим расчетам. М.: Машиностроение, 1975. — 559 е., ил.
  50. В.А., Позняк Э. Г. Аналитическая геометрия. М.: Наука, 1968. -232 с.
  51. В.Ф. Основы теории поверхностей. М.: ОГИЗ, Гостехиз-дат, 1947.-Т.1.-512 е.- 1948. — т.2. — 408 с.
  52. В.М. Каркас поверхностей изогнутой отсасывающей трубы гидротурбин // Труды УДН, Математика, Прикладная геометрия, М., 1971. — t.LIII. -вып.4.-с. 98−102.
  53. H.H. Численные методы. М.: Наука, 1978. — 318 с.
  54. М.Д. Численные методы строительной механики. -М.: Наука, 1976.-520 с.
  55. JI.B. Приближенный расчет железобетонных спиральных камер высоконапорных ГЭС // Применение железобетона в машиностроении. Сб. статей. -М.: Машиностроение, 1964. с.458−472.
  56. П.А. Примеры расчета упругих оболочек. М., 1966., -132 с.
  57. C.B. Параметрические уравнения некоторых спирои-дальных поверхностей // Кибернетика, графика и прикладная геометрия поверхностей. Труды МАИ. -М., 1974. вып.296. — с. 112−124.
  58. В.А. Плоская задача теории упругости. М.: Высш. шк, 1976, — 152 с.
  59. Н.В. Основы расчета упругих оболочек: Учеб. пособие. М.: Высш. шк., 1972. — 296 с.
  60. М.С., Паймушин В. Н., Андреев C.B. К вариационным методам исследований устойчивости тонких оболочек сложной геометрии // Тр. XII Всесоюзной конф. по теории оболочек и пластин. Ереван: Изд-во Ереванского университета, 1980. — т.1. — с.67−72.
  61. М.С., Паймушин В. Н., Снигирев В. Ф. Вычислительная геометрия в задачах механики оболочек. М.: Наука, 1989. — 208 с.
  62. М.С., Паймушин В. Н., Фирсов В. А. К решению двухмерных задач механики деформирования оболочек сложной геометрии // Вопросы вычислительной и прикладной математики. Ташкент: Изд-во АН УзССР, 1980. — вып.60. — с.70−76.
  63. М.С., Файзуллина М. А. Обзор работ по расчету на изгиб и устойчивость пластин и оболочек сложного очертания. Казан, физ.-техн. ин-т. Казань., 1986. — 36 с. — Рук. деп. в ВИНИТИ N8071- В86 от 1.12.86г.
  64. И.И. Об одном методе исследования циклических поверхностей // Труды ВЗЭИ, -М., 1958. вып. 13. — с. 51−60.
  65. Р. Строительство оболочек в Австрии // Большепролетные оболочки: Межд. конгресс в Ленинграде. М.: Стройиздат, 1969. -т.2.-с.319−320.
  66. С.Н. Основные сведения из дифференциальной геометрии поверхностей. М.: РУДН, 1992. — 32 с.
  67. Кришна Редди Г. В. Безмоментная теория оболочек в форме цик-лид Дюпена // Исследования по теории сооружений. М., 1967. — вып. 15.- с.90−99.
  68. Кришна Редди Г. В. Расчет оболочек в форме циклид Дюпена: Дис.. канд. техн. наук. -М.: УДН, 1966. 157 с.
  69. В.А. Тонкостенные зонтичные оболочки. -Л.: Госстрой-издат, 1958.- 172 с.
  70. М.Е. Определение усилий в кольцевой оболочке, образованной вращением вокруг вертикальной оси части окружности // Пространственная работа промышленных зданий и инженерных сооружений. -Л., 1982. с.60−73.
  71. А. И. Общая теория упругих тонких оболочек // ПММ.- 1940. т.4. — вып.2. — с.7−34.
  72. Ляв А. Математическая теория упругости. М., Л.: ОНТИ, 1935.- 674 с.
  73. Махмуд Хуссейн Аль-Хадж. Решение задач расчета тонких упругих оболочек в форме эпитрохоидальных поверхностей: Дис.. канд. техн. наук. М.: РУДН, 1992. — 157 с.
  74. В.И. Дифференциальная геометрия. Л., 1934.-332 с.
  75. В.Е., Ковалев С. Н. Конструирование форм современных архитектурных сооружений. Киев: Бущвельник, 1978. — 112 с.
  76. В.Е., Обухова B.C., Подгорный А. Л. Формообразование оболочек в архитектуре. Киев: Буд1вельник, 1972. — 206 с.
  77. С.Г. Вариационные методы математической физики. -М.: Наука, 1970.-360 с.
  78. H.H. Численные методы решения некоторых задач теории упругости. Киев: Наукова думка, 1979. — 316 с.
  79. Наср Юнее Ахмед Аббуши. Волнообразные купола // Строительная механика инженерных конструкций и сооружений: Межвузовский сборник научных трудов. М.: АСВ, 2002. — вып.11. — с.49−58.
  80. В.В. Краткий очерк развития теории оболочек в СССР // Исследования по теории пластин и оболочек. Казань.: Изд-во Казан. ун-та, 1970. — вып.6−7. — с.5−22.
  81. В.В. Новый метод расчета тонких оболочек // Изв. АН СССР. ОТН. 1946. -№.1. -с.35−48.
  82. В.В. Теория тонких оболочек. Л.: Судостроение, 1962.-431 с.
  83. В.В., Финкельштейн P.M. О погрешности гипотезы Кирхгофа в теории оболочек //ПММ. 1943. -т.7. — вып.5. — с.331−340.
  84. В.В., Черных К. Ф., Михайловский Е. И. Линейная теория тонких оболочек. Л.: Политехника, 1991. — 656 с. ил.
  85. А.П. Теория поверхностей. М.: ГИТТЛ, 1956. — 260 с.
  86. В.Я. Безмоментное напряженное состояние торойдального покрытия // Исследования по теории упругости и пластичности. -Л.: ЛГУ, 1961.-№ 1,-с. 110−118.
  87. В.Н. Некоторые задачи статики незамкнутых оболочек сложной формы и об одном методе их численного решения // Вопросы расчета прочности конструкций летательных аппаратов. Казань: Изд-во Казанского авиационного института, 1979. — с.67−76.
  88. В.Н., Андреев С. В. К численному исследованию напряженно-деформированного состояния однослойных и трехслойных пластин и оболочек сложной геометрии // Прикл. мех. Киев, 1983, — т.7. — с.24−30.
  89. Ф. А. Расчет железобетонных спиральных камер водяных турбин круглого сечения. // Применение железобетона в машиностроении: Сб. статей. -М.: Машиностроение, 1964. с. 440−457.
  90. В.Г. Развитие некоторых разделов теории расчета тонких оболочек // Труды УДН «Строительство»: Строительная механика. М.: УДН, 1967. — t.XXVIII. — вып.З. — с.3−5.
  91. В.Г., Кривошапко С. Н. Расчет оболочек сложной геометрии. М.: Изд-во УДН, 1988. — 178 с.
  92. В.Г., Рыжов H.H. Некоторые возможности расширения круга задач по конструированию и расчету оболочек // Труды УДН «Строительство»: Строительная механика. М.: УДН, 1970. — t.XLVIII. — вып.6.- с.3−8.
  93. Л.А. Задачи теории упругости и численные методы их решения. -СПб.: Изд-во СПбГУ, 1998. 532 с.
  94. Н.П. Кинематика образования циклической поверхностей // Труды УДН «Математика»: Прикладная геометрия. М.: УДН, 1967.- t.XXVI. вып.З. — с. 100−104.
  95. H.H. Каркасная теория задания и конструирования поверхностей // Труды УДН «Математика»: Прикладная геометрия. М.: УДН, 1967. — т. XXVI. — вып.З. — с. 3−16.
  96. H.H., Гершман И. П., Осипов В. А. Прикладная геометрия поверхностей // Тр. Моск. Научн.-метод, семинара по начертательной геометрии и инженерной графике. -М.: МАИ, 1972. ч.З. — вып.242. — с. 57−91.
  97. Сальман Аль-Духейсат. К вопросу расчета без моментной оболочки в виде отсека циклиды Дюпена // Исследования по расчету элементов пространственных систем. М.: УДН, 1987. — с. 106−110.
  98. В.И. Основы теории упругости и пластичности. М.: Высш. гик., 1982.-264 с.
  99. М. Метод конечных элементов: Пер. с серб. М.: Стройиздат, 1993. — 664 с. ил.
  100. И.А. Аналитическая теория формообразования оболочек // Труды международной научной конференции «0болочки-2001». М: Изд-во РУДН, 2001. — с.366−371.
  101. И.А. Обобщенные цилиндрические координаты и их применение в прикладной геометрии // Прикладная геометрия и инженерная графика, Киев, 1971. — вып. 3. — с. 112−118.
  102. Ю.В., Буданов В. Д., Гапоев М. М. и др. Конструкции из дерева и пластмасс. -М.: Стройиздат, 1986. 543 с. ил.
  103. В.Т. Об одном методе задания частного вида эпитро-хоидальных поверхностей // Прикладная геометрия и инженерная графика. -Киев, 1975. -вып.20. -с.89−91.
  104. М.Н. Железобетонные камеры кольцевого сечения // Энергомашиностроение, М., 1959. — № 10. — с.26−32.
  105. С.П., Войновский-Кригер С. Пластины и оболочки. М.: Физматгиз, 1963. — 635 с.
  106. Фортран 90. Международный стандарт: Пер. с англ. -М.: Финансы и статистика, 1998. 378 с.
  107. Р.В. Численные методы для научных работников и инженеров: Пер. с англ. М.: Наука, 1972. — 400 с.
  108. В. Исследование работы в области теории оболочек в ФРГ // Второй Международный конгресс по тонкостенным оболочкам-покрытиям. -Осло, 1957.-М., 1960.
  109. .К. Некоторые оболочки, возведенные в Индии // Большепролетные оболочки: Международный конгресс в Ленинграде. М.: Стройиздат, 1969. — т. 1. — с. 207−220.
  110. К.Ф. Линейная теория оболочек. Л.: Изд-во ЛГУ, 1962. -Т.1.-374 е.- 1964.-Т.2.-395 с.
  111. В.И. Классическая дифференциальная геометрия.- М.: ГИФМЛ, 1963. 540 с.
  112. Л.Э. Дифференциальные уравнения и вариационное исчисление. -М: Наука, 1965. 279 с.
  113. A.M. Исследование аналитического метода задания циклид Дюпена при выделении их из конгруэнции окружностей // Прикладная геометрия. М.: УДН, 1971. — вып.4. — с.26−40.
  114. A.M. Циклические каркасы из линий кривизны // Труды УДН «Математика»: Прикладная геометрия. М.: УДН, 1967.- т.XXVI. вып.З. — с.66−78.
  115. Н.М. Об одном методе расчета оболочек сложной геометрии // Исслед. по теории оболочек: Тр. семинара. Казань: Казанский физико-технический институт, 1984. — вып. 17. — ч.2. — с.4−17.
  116. Н.М. О некоторых работах по расчету оболочек сложной геометрии // Тр. семинара. Казань: Казанский физико-технический институт, 1990. — № 25. — с.86−94.
  117. Н.М. Фрагменты оболочек сложной геометрии в тороидальной системе координат // Исследования по теории оболочек. Труды семинара. Казань: Казанский физико-технический институт, 1988. — вып.21, -ч.1.-с. 130- 137.
  118. Н.М., Галимов Ш. К., Хисматуллин Н. И. От каменных глыб к тонкостенным конструкциям. Казань: Изд-во «SOS», 2001. — 96 с.
  119. Г. Й. Вълнообразни черупкови конструкции. София: Изд-во на Българската академия на науките, 1961. — 80 с.
  120. Г. Й. Някои въпроси от моментната и безмоментната теория на черупките. Трудове на Научноизследователския строителен институт София, 1959. — № 2.
  121. Argyris J.H. Matrix methods of structural analysis // Proc. 14-th meeting of AGARD. AGARDograph. 1962. — 72.
  122. Barony S.Y., Totlenham H. The analysis of rotational shells using a curved ring element and the mixed variational formulation // Int. J. Numer Meth. Eng. 1976.- 10.-N4.-p.861−872.
  123. Brebbia C.A., Hadid H.A. Analysis of plates and shells using finite elements //Pev. roum. sci techn. ser. mec. appl.- 1973. 18. -N15. — p.939−962.
  124. Caley A. On The Cyclide // Quarterly Journal of pure and applied mathematics. 1873. -№ 12. — p. 148.
  125. Darboux G. System ortogonoux et les coordonees curilinges. 2nd adn. Paris, Guatvier-Villars. -1910.
  126. Dean W.R. The distortion of curved tube due to internal pressure. Phil. Mag VII ser. Vol. 28. — 1939.
  127. Dixon R. Asymmetric shells a new approach // Bulletin of international Association for Shell and Spatial Structures. -1991. — Vol. 32. -№ 3. — p. 133−137.
  128. Doehlemann K. Geometrische trasformation, II, Leipzig. 1908.
  129. Dupin Ch. Applications de geometry et de mechanique, Paris, Bachelier. 1822.
  130. Dupin Ch. Develepment de geometrie. Paris. 1813.
  131. Forsyth A.R. Lectures on the differential geometry of curves and surfaces. Cambridge.-1920
  132. Kikuchi F., Ando Y. A new variational functional for the finite element method and its application to plate and shell problems // Nucl. Eng. Design. -1972. -N25. -p.95−113.
  133. Maxwell J.C. On the Cyclide // Quarterly Journal of pure and applied mathematics, -1868. № 9.
  134. Nelson R.L. Stresses in shell structures // J. Sound and Vibr. 1981. -79. -N3. -p.397−414.
  135. Skopinsky V.N. Stress analysis of shell intersections with torus transition under internal pressure leading. // Trans. ASME. J. Pressure Vessel Technol. 1997. — 119, № 3. — p.288−292.
  136. Jiirgen Joedicke. Shell Architecture. London: Alec Tiranti Ltd., 1962.-304 c.152
  137. КАНАЛОВЫЕ ПОВЕРХНОСТИ ИОАХИМСТАЛЯ С РАЗЛИЧНЫМИ ОЧЕРТАНИЯМИ В ПЛАНЕ1. г2{а) = 1 + ?исо$>{а) эпитрохоидальная поверхность• С2= 0, //-0.5.155• С2 =0, а =, ju= 0.8, к = 2.25.156
  138. ОБОЛОЧКИ НА ОСНОВЕ КАНАЛОВЫХ ПОВЕРХНОСТЕЙ ИОАХИМСТАЛЯ, ПРОЕКТИРУЕМЫЕ В ДИПЛОМНЫХ РАБОТАХ СТУДЕНТОВ СТРОИТЕЛЕЙ РУДН
  139. ВОЗМОЖНОСТЬ КОНСТРУИРОВАНИЯ КОРПУСА ПРИТОЧНОЙ НАСАДКИ В ФОРМЕ КАНАЛОВЫХ ПОВЕРХНОСТЕЙ1. ИОАХИМСТАЛЯ
  140. Одним из примеров, когда геометрия конструкции задается технологическим требованием, является предложенный нами вариант замены существующей конструкции корпуса приточной насадки, имеющей прямоугольную форму сечения (рис. 2.21).
  141. Проточные насадки или воздухораспределители. I) статический (улиточный) центробежный конструкции ВНИИГС. II) конструирование корпуса в форме каналовых поверхностей Иоахимсталя. а приточный, б — приточно-вытяжной.
  142. ВЛИЯНИЕ ПАРАМЕТРА С НА ВИД КАНАЛОВОЙ ПОВЕРХНОСТИ
  143. ИОАХИМСТАЛЯ ПРИ г2 (а) = д//Г соз2 ка +р1 соз2 ка + а4 М- 0.5 < С2 <0.5 а = 1, // = 0.8,? = 2.5с2 =с2 = с2 = 00:43 19.09.2001 г. ! м м м м | м I м I м — I | м | м м I I м м м м ! ! м ! м м м | | I м I м м I I м I I
  144. РАСЧЕТ ТОНКОСТЕННЫХ КОНСТРУКЦИЙ ВАРИАЦИОННО-РАЗНОСТНЫМ МЕТОДОМ Авторы: проф. ИВАНОВ В.Н.
  145. НАСР ЮНЕС (Палестина) ПРОГРАММА ПОДГОТОВЛЕНА НА КАФЕДРЕ СОПРОТИВЛЕНИЯ МАТЕРИАЛОВ РОССИЙСКОГО УНИВЕРСИТЕТА ДРУЖБЫ НАРОДОВ МОСКВА, 2000 г.
  146. М II II II II II I II I I I I II II I II I II II I I I II II II I I I II II I II II II I II II II I I I
  147. PROGRAM SHELLVRM USE GEOMETRY1. PLICIT REAL*8 (A-H, 0-Z)
  148. REAL, PARAMETER: PI=3.1 415 926 535 897 932 412 023 218 372 608 DIMENSION AH (30 000 000), MN (6), MH (300 000), MZ (3) CHARACTER*12 NMR, NMW, СDAT, NP1, NP2,NP3,NP4,NP5 COMMON NTP, NGEO, N1, N2,KU, KW, KUS, KN, LS, NT1, NT2,NTR, PA (12) PRINT 800 PRINT 900
  149. PRINT *, ' Последнее обновление программы 00:43 19.09.2001' PRINT*
  150. PRINT*,'ВВОДИТЕ ИМЕНА ФАЙЛОВ ИСХОДНЫХ ДАННЫХ И РЕЗУЛЬТАТОВ' READ (*,*) NMR, NMW! ввод имен файлов
  151. OPEN (UNIТ=5,FILE=NMR, STATUS='OLD 1) !c исходными данными
  152. OPEN (UNIT=6,FILE=NMW, STATUS='NEW') !и результатами.
  153. WRITE (*,*) ' &bdquo-АТА' READ (*,*) СDAT WRITE (6,1000) СDAT
  154. WRIТЕ (б, *) ' !!! 1 !!! 1 1 !!! 1 1 1 1 !!! 1
  155. M I I | I I | I I | | | M? I | I i
  156. WRITE (б,*)' !! РАСЧЕТ ТОНКОСТЕННЫХ КОНСТРУКЦИЙ ВАРИАЦИОННО-РАЗН *ОСТНЫМ МЕТОДОМ !!'
  157. WRITE (б,*)' !! Авторы: проф. ИВАНОВ В.Н., НАСР ЮНЕС (Палестина)
  158. WRITE (б, *) ' I mi I м I м I I I I м м I м I I I I м м I м I I I I I I м I I I I I I м I м I I
  159. X | I I I I I I I И II I I I I I
  160. WRITE (6,1010) NMR С NTP-номер типа задачи: 1-ИЗГИБ ПЛАСТИН,
  161. С 2-ПЛОСКАЯ ЗАДАЧА ТЕОРИИ УПРУГОСТИ, 1. С 3-ОТСЕКИ ОБОЛОЧКИ, 1. С 4-ЗАМКНУТЫЕ ОБОЛОЧКИ.
  162. С NGEO-номер геометрии задачи: 1-ПРЯМОУГОЛЬНАЯ СИСТЕМА КООРДИНАТ,
  163. С 2-ПОЛЯРНАЯ СИСТЕМА КООРДИНАТ,
  164. С 3-ЦИЛИНДРИЧЕСКИЕ ИЛИ ПОЛОГИЕ ОБОЛОЧКИ, 1. С 4-ОБОЛОЧКИ ВРАЩЕНИЯ, 1. С 5-СФЕРИЧЕСКАЯ ОБОЛОЧКА,
  165. С 6-КАНАЛОВЫЕ ПОВЕРХНОСТИ ИОАХИМСТАЛЯ
  166. READ (5,*) NTP, NGEO, N1, N2,KL С WRITE (б,*) NTP, NGEO, N1, N2,KL1. NTR=NTP1.(NTP, EQ.4)NTR=3 NT1 = 1 NT2 = 21.(NTP.EQ.2)NT2=1 IF (NTP.EQ.1)NT1=2
  167. WRITE (6 *)' ***************************************************** ******************** '1. ПЛОСЕАЯ ДЕФОРМАЦИЯ
  168. РАСЧЕТ ОТСЕКОВ (СЕКТОРОВ) ОБОЛОЧКИ РАСЧЕТ ЗАМКНУТОЙ ОБОЛОЧКИ1. SELECT CASE (NTP) CASE (1)
  169. WRITE (б, 1018) ' РАСЧЕТ ПЛАСТИНКИ НА ИЗГИБ1. CASE (2)
  170. WRITE (6,1018) ' РАСЧЕТ ПЛОСКОЙ ЗАДАЧИ ТЕОРИИ УПРУГОСТИ1. SELECT CASE (KL) CASE (0)
  171. WRITE (б, 1017) ' ПЛОСКОЕ НАПРЯЖЕННОЕ СОСТОЯНИЕ
  172. CASE (1) WRITE (б, 1017) END SELECT CASE (3) WRITE (6,1018) CASE (4) WRITE (6,1018) CASE DEFAULT
  173. STOP 'HE ПРАВИЛЬНЫЙ ВВОД НОМЕРА ТИПА ЗАДАЧИ' END SELECT SELECT CASE (NGEO) CASE (1)
  174. WRITE (6,1019) ' ПРЯМОУГОЛЬНАЯ СИСТЕМА КООРДИНАТ
  175. CASE (2) WRITE (6,1019) CASE (3) WRITE (6,1019) CASE (4) WRITE (6,1019) CASE (5) WRITE (6,1019) CASE (6) WRITE (6,1019)
  176. ПОЛЯРНАЯ СИСТЕМА КООРДИНАТ ЦИЛИНДРИЧЕСКАЯ (ИЛИ ПОЛОГАЯ) ОБОЛОЧКА ОБОЛОЧКА ВРАЩЕНИЯ СФЕРИЧЕСКАЯ ОБОЛОЧКА КАНАЛОВАЯ ПОВЕРХНОСТЬ ИОАХИМСТАЛЯ1. CASE DEFAULT
  177. STOP 'НЕ ПРАВИЛЬНЫЙ ВВОД НОМЕРА ГЕОМЕТРИИ ЗАДАЧИ' END SELECT
  178. WRITE (6 ' * ^ * * * * ^ * * ^ * * * * * * * * * ^ * * ^ ^ ^ к -к к к к к к к к к к к к к к -к к к к к * к к к кк к *к-ккккккккккккккк к к~кк I
  179. NLl=NHl+Nl+2 NH2=NL1+N1+1 NL2=NH2+N2+2 NAS=NL2+N2+1 NMNH=11. AH (NL1)=PA (3)
  180. KUS= (N1 + 3)* (N2 + 3)+NTR1. D=', PA (7)
  181. Адрес (@) каждого массива в едином массиве АН ЫН1−0 Н1(м.шагов сетки по X) NLl-@ НЬ1(м.длин (коорд.)по X) NH2-@ Н2(м.шагов сетки по У) NL2-@ HL2(м.длин (коорд.)по Y) количество общих узловыхг
  182. WRITE (6,*) 'KUS=', KUS WRITE (6,1040) DLSUX-0
  183. READ (5,*) I, J, A! c номера сетки I no J
  184. WRITE (6,1070) 60 READ (5,*) MN, MZ
  185. WRITE (6,*) 'KN=', KN С IK=11. С DO 120 IW=-1,N1+1
  186. С WRITE (6,82) (MH (К), K=IK, IK+(N2+3)*NTR-1)
  187. С 120 IK=IK+(N2+3)*NTR CALL LSI (MH (NMNH)) WRITE (6,*) 'LS=', LS NAR=NAS+KN*(LS+1) NAP=NAR+KN NMX=NAP+LS+1 WRITE (*,*) 'NMX=', NMX WRITE (*,*) 'KN=', KN
  188. CALL SHLLESF (AH (NH1), AH (NH2), AH (NL1), AH (NL2), MH (1), AH (NAS)) WRITE (*,*) 'esf'
  189. CALL SHLLOAD (AH (NH1), AH (NL1), AH (NH2), AH (NL2), MH (1), AH (NAR)) WRITE (*,*) 'load' NP=0 NT = 1
  190. CALL DSLT (AH (NAS), AH (NAR), AH (NAP), KN, LS, NT) WRITE (*,*) 'dslt' С WRITE (6,*) (AH (I), I=NAR, NAR+KN-1)
  191. CALL SHLLSORT (AH (NAR), AH (NAP), MH (1)) WRITE (*,*) 'sort'
  192. CALL SHLLRES (AH (NAP), AH (NH1), AH (NH2), AH (NL1), AH (NL2)) WRITE (*,*) 'result'
  193. SUBROUTINE SHLLESF (HI, Н2, HL1,HL2, MNH, AS) USE GEOMETRY1. PLICIT REAL*8 (A-H, 0-Z)
  194. DIMENSION HI (0:N1+1), H2(0:N2+1), HL1(0:N1), HL2(0:N2)
  195. DIMENSION AS (1), DF (6,9,4), HN (3,3), DS (4), DO (6,9,4)
  196. DIMENSION CF (12), EG (3, 6), AP (200)
  197. DIMENSION MNH (KU:KW,-1:N2+1,-1:N1+1)
  198. DIMENSION RM (9,9), RM1(9,9), KI (4), KPSL (7), PAS (3)
  199. DIMENSION EGN (3,6), EGT (6,6), EGDF (6,9)
  200. COMMON NTP, NGEO, N1, N2,KU, KW, KUS, KN, LS, NT1, NT2,NTR, PA (12) WRITE (6,*) 'PA=' WRITE (6,*) PA
  201. HU=PA (2) !Nu-коэф. Пуассона1. HN (1, 1)=11. HN (1, 2)=HU1. HN (1, 3)=01. HN (2,1)=HU1. HN (2,2)=11. HN (2,3)=01. HN{3,1)=01. HN (3,2)=01. HN (3,3)=(1.-HU)/21. DO 4 K=l, 41. DO 4 1=1,б1. DO 4 J=l, 91. DO (I, J, K) =01. DO DO DO DO DO DO DO DO1. CONTINUE
  202. DO 200 1=0,N1 DO 200 J=0, N2 X = HL1 (I) Y=HL2(J)
  203. CALL COEF (X, Y, CF, PA, NGEO) AB=CF (1)*CF (2)/4 C WRITE (*,*) 'X=', X, ' Y= ', Y
  204. DS (1)=H1(1+1)*H2(J+l)*AB DS (2)=H1(I)*H2(J+l)*AB DS (3) =H1(I)*H2(J)*AB DS (4)=H1(1+1)*H2(J)*AB IF (I.EQ.O.OR.J.EQ.O) DS (3)=0 IF (I.EQ.O.OR.J.EQ.N2) DS (2)=0 IF (I.EQ.Nl.OR.J.EQ.O) DS (4)=0 IF (I.EQ.Nl.OR.J.EQ.N2) DS (1)=0
  205. C IF (I.EQ.1.AND.J.EQ.1) THEN1. C WRITE (6,*)'DS='1. C WRITE (6,*) DS1. C ENDIF1. C WRITE (*,*) PA
  206. C WRITE (*,*) AO, Al, A2, A3
  207. C WRITE (*,*) BO, Bl, B2, B3
  208. C WRITE (*,*) A1B1, A1B2, A2B1, A2B21. C END IF1. DO 8 K=1,4
  209. DF (4 2, K) =DF (4,2,K)*A3/(1+A0)
  210. DF (4 5, K) =DF (4,5,K)*A3/A0
  211. DF (4 8, K) =DF (4,8,K)*A3/(AO*(1+A0))
  212. DF (6 4, K) =DF (6,4,K)*B3/(1+B0)
  213. DF (6 5, K) =DF (6,5,K)*B3/B0
  214. C 301 WRITE (6,302) (DF (L1,L3,L), L3=1,9)302 FORMAT (9E9.2) C END IF
  215. DO 16 K=KU, KW DO 15 L=KU, KW
  216. DO 11 11=1,9 DO 11 Jl=l, 9 11 RM (II, J1)=01. DO 12 NT=NT1,NT2
  217. CALL EGFMATRIX (CF, EG, NT, L)
  218. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1) THEN
  219. C WRITE (6,*) 'L=', L, ' NT=', NT, ' EG', (NT-1)* 3 + L1. C DO 11=1,3
  220. C WRITE (6,302)(EG (II, Jl), Jl=l, 6)1. C ENDDO1. C ENDIF
  221. CALL MULTMATRIX (HN (1,1), EG (1,1), AP (1), 3,3,3,6)
  222. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1) THEN
  223. C WRITE (6,*) 'L=', L, ' NT=', NT, ' NEG', (NT-1) *3 + L1. C DO IL=1,31. C KJ=IL1. C DO JL=1,61. C EGN (IL, JL)=AP (KJ)1. C KJ=KJ+31. C ENDDO1. C ENDDO1. C DO IP=1,3
  224. C WRITE (6,302) (EGN (IP, JP), JP=1, 6)1. C ENDDO1. C ENDIF
  225. CALL EGFMATRIX (CF, EG, NT, K)
  226. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1) THEN
  227. C WRITE (6,*)'K=', K,' NT=', NT,' EG',(NT-1)*3+K1. C DO 11=1,3
  228. C WRITE (6,302)(EG (II, Jl), Jl=l, 6)1. C ENDDO1. C ENDIF
  229. CALL MULTMATRIX (EG (1,1), AP (1), AP (19), 3,6,3,6)
  230. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1) THEN
  231. C WRITE (6,*) 'L=', L, ' NT= ', NT, 1 EGT*NEG', (NT-1) *3 + L1. C DO IL=1,61. C KJ=IL1. C DO JL=1,61. C EGT (IL, JL)=AP (KJ+18)1. C KJ=KJ+61. C ENDDO! JL1. C ENDDO ! IL1. C DO IP=1,6
  232. C WRITE (6,302)(EGT (IP, JP), JP=1,6)1. C ENDDO IIP C ENDIF
  233. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1) THEN
  234. C WRITE (6,*)'L=', L,' NT=', NT, 1 EGT*NEG',(NT-1)*3+L1. C KI = 1 C DO IW=1,6
  235. C WRITE (6,302) (AP (JI), JI=KI, KI+5)1. C KI=KI+61. C ENDDO1. C ENDIF1. DO 12 IDF=1,4
  236. CALL MULTMATRIX (AP (19), DF (1,1,IDF), AP (55), 6,6,6,9)
  237. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1.AND.L.EQ.2.AND.NT.EQ.1) THEN
  238. C WRITE (6,*)'L=', L,' NT=', NT,' EGT*NEG*DF', IDF1. C DO IL=1,61. C KJ=IL1. C DO JL=1,91. C EGDF (IL, JL)=AP (KJ+54)1. C KJ=KJ+61. C ENDDO1. C ENDDO1. C DO IP=1,6
  239. C WRITE (6,302)(EGDF (IP, JP), JP=1,9)1. C ENDDO1. C ENDIF
  240. CALL MULTMATRIX (DF (1,1,IDF), AP (55), RM1(1,1), 6,9,6,9)
  241. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1.AND.L.EQ.2.AND.NT.EQ.1) THEN
  242. C WRITE (6,*)'L=', L,' NT=', NT,' RM1', IDF1. C DO IP=1,9
  243. C WRITE (6,302)(RM1(IP, JP), JP=1,9)1. C ENDDO C ENDIF1. DO 11=1,9 DO Jl=l, 9
  244. RM (II, Jl)=RM (II, Jl)+RM1(II, Jl)*DS (IDF)*PA (NT+5)1. ENDDO1. ENDDO
  245. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1.AND.L.EQ.2.AND.NT.EQ.1) THEN
  246. C WRITE (6,*)'L=', L,' NT=', NT,' RM',' IDF=', IDF1. C DO IP=1,9
  247. C WRITE (6,302)(RM (IP, JP), JP=1,9)1. C ENDDO1. C ENDIF
  248. CONTINUE C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1) THENc c c c c c
  249. WRITE (6,*) 'K=', K, ' L= ', L WRITE (6, *)'RM' DO IP=1,9
  250. WRITE (6,302)(RM (IP, JP), JP=1,9) ENDDO ENDIF1.(K.NE.3.AND.L. NE.3) GOTO 130 LK0 = 111 = 0
  251. IF (I.NE.0.AND.I.NE.N1) GOTO 90 KJ0=11.(I.EQ.Nl) KJ0=-1 LK2=31. DO 80 J1=J-1,J+11.(MNH (3,Jl, I-KJ0).NE.0) GOTO 80 LK2=LK2+1 LK1 = 1
  252. KI (3-LKO)=LK2+KJ0*3 KI (5-LKO)=LK2-KJ0*3 DO 70 Kl=l, 9 KI (LK0)=K1 KI (LKO+2)=K17 0 RM (KI (1), KI (2))=RM (KI (1), KI (2))+RM (KI (3), KI (4)) 80 CONTINUE
  253. IF (J.NE.0.AND.J.NE.N2) GOTO 120 KJ0=11.(J.EQ.N2) KJ0=-1 LK2=-11. DO 110 11=1−1,1+11.(MNH (3,J-KJ0,II).NE.0) GOTO 1101.2=LK2+311 = 1
  254. KI (3-LKO)=LK2+KJ0 KI (5-LKO)=LK2-KJ0 DO 100 Kl=l, 9 KI (LK0)=K1 KI (LK0+2)=K1
  255. RM (KI (1), KI (2))=RM (KI (1), KI (2))+RM (KI (3), KI (4)) 110 CONTINUE 120 LK0=LK0+11.(LK0.EQ.2.AND.LK1.EQ.1) GOTO 60 130 L1=0
  256. DO 14 11=1−1,1+1 DO 14 J1=J-1,J+1 L1=L1+11. NU=MNH (K, Jl, II)
  257. C IF (I.EQ.1.AND.J.EQ.2.AND.L.EQ.1) THEN
  258. C WRITE (6,*)'K=', K,' 11=', II,' J1=', J1,' NU=', NU1. C ENDIF1.(NU.LT.l) GOTO 14 L2 = 0
  259. DO 13 12=1−1,1+1 DO 13 J2=J-1,J+l L2=L2+11. NH=MNH (L, J2,12)1.(NH.LT.1.OR.NH.LT.NU) GOTO 13 NK=(NU-1)*(LS+1)+NH-NU+1 AS (NK)=AS (NK)+RM (LI, L2)13 CONTINUE14 CONTINUE
  260. С IF (I.EQ.l.AND.J.EQ.l) THEN1. С WRITE (6,*) 'RM='1. С DO 319 13=1,9
  261. С 319 WRITE (6,302) (RM (13,14), 14=1,9) С ENDIF401 CONTINUE15 CONTINUE16 CONTINUE 200 CONTINUE
  262. ПРОЦЕДУРА УЧЕТА УПРУГОГО ОСНОВАНИЯ ПОД КОНСТРУКЦИЕЙ
  263. СПЛОШНОЕ УПРУГОЕ ОСНОВАНИЕ' ' УПРУГИЕ ОПОРЫ ВДОЛЬ ЛИНИИ' ' УПРУГИЕ УЗЛОВЫЕ ОПОРЫ'1. DO I=KPSL (2), KPSL (3)
  264. DO K=l, 3 NH=MNH (K, J, I) NK=(NH-1)*(LS+1)+11.(NH.GT.0.AND.PAS (K).NE.0) AS (NK)=AS (NK)+PAS (K)*DSA ENDDO ENDDO ENDDO GOTO 500 400 KP=KPSL (2)
  265. DO I=KPSL (3), KPSL (4), KPSL (5)
  266. DO J=KPSL (6), KPSL (7) 11 = 1 J1=J1.(KP.EQ.2) THEN I1 = J J1 = I ENDIF X=HL1(II)1. Y=HL2(Jl)
  267. DO I=KPSL (2), KPSL (3), KPSL (4) DO J=KPSL (5), KPSL (6), KPSL (7) DO K=l, 31. NH=MNH (K, J, I)
  268. NK=(NH-1)*(LS+1)+1 IF (NH.GT.0.AND.PAS (K).NE.0) AS (NK)=AS (NK)+PAS (K) ENDDO ENDDO ENDDO GOTO 500 510 CONTINUE
  269. FORMAT (15,5X, 615,5X, 3E11.3) C NH=11. C WRITE (6,*) 'AS='1. C DO NU=1,KN
  270. C WRITE (6,*) (AS (N1), NI=NH, NH+LS)1. C NH=NH+LS+11. C END DO1. C WRITE'(6, *) ' EG= '1. C DO 310 L0=1,21. C DO 310 Ll=l, 31. C WRITE (6,*) L0, Ll1. C DO 310 L2=l, 3
  271. C 310 WRITE (6,302) (EG (L0,L1,L2,L3), L3=1,6)1. RETURN END1. SUBROUTINE MHHl (MNH)
  272. DIMENSION MNH (KU:KW,-1:N2+1,-1:N1+1)
  273. COMMON NTP, NGEO, Nl, N2, KU, KW, KUS, KN, LS, NT1, NT2,NTR, PA (12) DO 3 K=KU, 2 DO 1 I=-l, N2+1 MNH (K, I, -1) =0 1 MNH (K, I, N1+1)=0 DO 2 J=-1,N1+1 MNH (K,-1,J)=02 MNH (K, N2 + 1, J)=0
  274. CONTINUE IF (KW.NE.3)GOTO 4 MNH (KW,-1,-1)=-1 MNH (KW,-1,N1+1)=-1 MNH (KW, N2+1,-1)=-1 MNH (KW, N2+1,N1+1)=-14 RETURN END1. SUBROUTINE LSI (MNH)
  275. DIMENSION MNH (KU:KW,-1:N2+1,-1:N1+1)
  276. COMMON NTP, NGEO, N1, N2,KU, KW, KUS, KN, LS, NT1, NT2,NTR, PA (12) LS=1
  277. DO 104 KS=KW, KU,-1 K2=MNH (KS, J+l, I) IF (K2.LT.1) GOTO 104 K3=K2-I<11.(K3.GT.LS) LS=K3104 CONTINUE105 CONTINUE106 CONTINUE RETURN END
  278. SUBROUTINE DSLT (AS, AR, AP, NS, LP, NT) IMPLICIT REAL*8(A-H, O-Z) DIMENSION AS (1), AR (1), AP (1) N=LP+11.(NT.NE.l) GOTO 5 DO 4 1=1,NS-1 M=(1−1)*N+11. DO 1 J=2,N
  279. AP (J) =AS (M+J-1) /AS (M) DO 3 K=I+1,I+LP IF (K.GT.NS) GOTO 4 M=K-I
  280. K1=(I-1)*N+M+1 IF (AS (KI).EQ.0) GOTO 3 LI-(K-l)*N DO 2 J=1,N-M
  281. AS (Ll+J)=AS (Ll+J)-AS (KI)*AP (J+M)3 CONTINUE4 CONTINUE5 DO 7 1=1,NS-11.(AR (I).EQ.0) GOTO 71. M=(1−1)*N + 11. A1=AR (I)/AS (M)1. DO 6 K=I+1,I+LP1.(K.GT.NS) GOTO 71. M=K-I+(1−1)*N1. AR (K)=AR (K)-AS (M+l)*A16 CONTINUE7 CONTINUE I=NS+18 1=1−11. Ml=(1−1)*N+1 M=M1
  282. DO 9 K=I+1,I+LP IF (K.GT.NS) GOTO 10 M=M+1
  283. AR (I)=AR (I)-AR (K)*AS (M) 9 CONTINUE 10 AR (I)=AR (I)/AS (Ml) IF (I.GT.l) GOTO 8 RETURN END
  284. SUBROUTINE MULTMATRIX (A, B, C, M, N, N1, K) IMPLICIT REAL*8(A-H, O-Z)
  285. DIMENSION A (1), B (1), C (1) M1=M ID=1 JD=M1.(N.EQ.Nl) GOTO 11. M1=N1.=M1. JD=11 K1=01. DO 2 J=1,K1. DO 2 1=1,Ml1. K1=K1+11. C (K1)=01. DO 2 11=1,N11.=(1−1)*ID+(I1−1)*JD+11.=N1*(J-l)+11
  286. C (KI)=C (KI)+A (IA)*B (IB) RETURN1. END
  287. SUBROUTINE EGFMATRIX (CF, EG, IP, K)1.PLICIT REAL*8(A-H, 0-Z)
  288. DIMENSION EG (3,6), CF (12) DO 17 Jl=l, 3 DO 17 Kl=l, 6 17 EG (J1, K1) =01. A12=CF (1)*CF (1)1. A22=CF (2)*CF (2)1. A1A2=CF (1)*CF (2)
  289. GOTO (1,2,3,4,5,6) (IP-1)*3+K
  290. EG (1,2)=1./CF (1) EG (2,1)=CF (7)/A1A2 EG (3,1)=-CF (6)/A1A2 EG (3,3)=1./CF (2)1. RETURN
  291. EG (1,1)=CF (6)/A1A2 EG (2,3)=1./CF (2)
  292. EG (3,1)=-CF (7)/A1A2 EG (3,2)=1./CF (1) RETURN
  293. EG (1,1)=CF (3) EG (2,1)=CF (4) RETURN
  294. EG (1,1)=CF (9)/CF (1) EG (1,2)=CF (3)/CF (1)
  295. EG (2,1)=CF (7)*CF (3)/A1A2 EG (3,1)=CF (10)*CF (1)/CF (2) EG (3,3)=CF (3)/CF (2) RETURN
  296. EG (1,1)=CF (6)*CF (4)/A1A2 EG (2,1)=CF (11)/CF (2)
  297. EG (2,3)=CF (4)/CF (2) EG (3,1)=CF (12)* CF (2)/CF (1) EG (3,2)=CF (4)/CF (1) RETURN
  298. EG1=CF (6)/A1A2 EG2=CF (7)/A1A21. EG (1,2)=CF (5)/(A12*CF (1))1. EG (1,3)=-EGl/CF (2)1. EG (1,4)=-1./A121. EG (2,2)=-EG2/CF (1)
  299. EG (2, 3)=CF (8)/(A22*CF (2))1. EG (2,6)=-1./A221. EG (3,2)=2*EG1/CF (1)1. EG (3,3)=2*EG2/CF (2)1. EG (3,5)=-2./A1A21. RETURN1. END
  300. SUBROUTINE SHLLOAD (HI, HL1, H2,HL2,MNH, AR) USE GEOMETRY1. PLICIT REAL*8 (A-H, O-Z)
  301. DO 2 0 I=MH (2), MH (3) DO 2 0 J=MH (4), MH (5) X=HL1(I) Y=HL2(J)
  302. HL1(MH (2)))+QL (3)*(Y-HL2(MH (4))))*DS30 KP=MH (2)1. DO I=MH (3), MH (4), MH (5)
  303. DO J=MH (6), MH (7) 11 = 1 J1=J1.(KP.EQ.2) THEN I1 = J J1 = I ENDIF X=HL1(II) Y=HL2(J1)
  304. CALL COEF (X, Y, CF, PA, NGEO) NH=MNH (QL (1), Jl, II) IF (NH.LT.l) CYCLE DS = 01. SELECT CASE (KP) CASE (1)1. AB=CF (2)/21.(Jl.GT.O) DS=DS+AB*H2(Jl) IF (Jl.LT.N2) DS=DS+AB*H2(Jl+1)
  305. AR (NH)=AR (NH)+(QL (2)+QL (3)*(Y-HL2(MH (6))))*DS CASE (2)1. AB=CF (1)/21.(I1.GT.0) DS=DS+AB*H1(II) IF (I1.LT.N1) DS=DS+AB*H1(11+1)
  306. AR (NH)=AR (NH)+(QL (2)+QL (3)*(X-HL1(MH (6))))*DS END SELECT ENDDO ENDDO1. GOTO 104 0 KP=MH (2)1. DO I=MH (3), MH (4), MH (5)
  307. DO J=MH (6), MH (7) 11 = 1 J1=J1.(KP.EQ.2) THEN I1=J J1 = I ENDIF X=HL1(I1) Y=HL2(Jl)
  308. CALL COEF (X, Y, CF, PA, NGEO) DS = 01. SELECT CASE (KP) CASE (1)
  309. A=H2(Jl+1)/(H2(Jl)) AK (1)=-A/((1+A)*H2(Jl)) AK (2)=(A-l)/(A*H2(Jl)) AK (3)=1./(A*(1+A)*H2(Jl+1)) AB=CF (1)/21.(Il.GT.O) DS=DS+AB*H1(II) IF (Il.LT.Nl) DS=DS+AB*H1(I1+1) DO K=l, 3 DO 12=1,3 NH=MNH (K, Jl, II—2+12) IF (NH.LT.l) CYCLE
  310. AR (NH)=AR (NH)+(QL (2)+QL (3)*(X-HL1(MH (6))))*DS*AK (I2) ENDDO ENDDO END SELECT ENDDO1. ENDDO • •1. GOTO 105 0 DO I=MH (2), MH (3), MH (4)
  311. DO J=MH (5), MH (6), MH (7) DO K=l, 31. NH=MNH (K, J, I)1.(NH.GT.0.AND.QL (K).NE.0) AR (NH)=AR (NH)+QL (K) ENDDO ENDDO ENDDO GOTO 1060 DO I=MH (2), MH (3), MH (4)
  312. DO J=MH (5), MH (6), MH (7) A=Hl (1+1)/(HI (I)) AK (1)=-A/((1+A)*H1(I)) AK (2)=(A-1)/(A*H1(I)) AK (3)=1./(A*(1+A)*H1(1+1)) B=H2(J+l)/(H2 (J)) BK (1)=-B/((1+B)*H2(J)) BK (2)=(B-l)/(B*H2(J)) BK (3)=1./(B*(1+B)*H2(J+l)) DO K=l, 3 DO 12=1,3
  313. NH=MNH (K, J, I-2+I2) IF (NH.LT.l) CYCLE AMX=QL (2)*AK (I2) AMY=QL (3)*BK (I2) AR (NH)=AR (NH)+AMX+AMY ENDDO ENDDO ENDDO ENDDO GOTO 107 0 G0=QL (1) DO K=l, 3 GXYZ (K)=0 GRIV (K)=01. GRIV (K)=G0 *MH (K+1) ENDDO
  314. DO 1=0,N1 DO J=0,N2 X=HL1 (I) Y=HL2(J)
  315. FORMAT (15,5X, 615,5X, 3F10.2) RETURN END
  316. SUBROUTINE SHLLSORT (AR, ARS, MNH) IMPLICIT REAL*8(A-H, O-Z)
  317. DIMENSION AR (1), ARS (1), MNH (KU:KW,-1:N2 + 1, -1:N1 + 1) COMMON NTP, NGEO, N1, N2,KU, KW, KUS, KN, LS, NT1, NT2,NTR, PA (12) DO 10 1=1,KUS 10 ARS (I)=0 N=0
  318. C WRITE (6,*) (ARS (N), N=K, K+NTR-1)1. C 21 K=K+NTR RETURN END
  319. SUBROUTINE SHLLRES (AR, HI, H2, HLl, HL2)1. USE GEOMETRY1. PLICIT REAL*8(A-H, O-Z)
  320. WRITE (6,*) 'X Y W Mx My H'
  321. WRITE (6,*) '=====================================================11. CASE (2)
  322. WRITE (6,*) ' X Y U V Nx Ny1. S'
  323. WRITE (6,*) 1 ================================================== ===1. CASE DEFAULT1. WRITE (6,*) X Y U V W Nx1. Ny S '
  324. WRITE (6,*) '=====================================================1. END SELECT1. CASE (2)1. SELECT CASE (NTP) CASE (1)
  325. WRITE (6,*) 'X Y W SigMx SigMy TawXY
  326. WRITE (6,*) '=====================================================1. CASE (2)
  327. WRITE (6,*) X Y U V SigNx SigNy1. TawXY'
  328. WRITE (6,*) '======================================================1. CASE DEFAULT
  329. WRITE (6,*) X Y U V W SigNx1. SigNy TawXY'
  330. WRITE (6,*) '======================================================
  331. END SELECT END SELECT JH=1 GOTO 11 C 10 IF (NTP.GE.3) THEN10 WRITE (6, '(/)')
  332. SELECT CASE (NPRZ) CASE (1)1. WRITE (6,*) X Y U V W Mx1. My H'
  333. WRITE (6,*) '=====================================================1. CASE (2)
  334. WRITE (6,*) X Y U V W SigMx1. SigMy TawXY'
  335. WRITE (6,*) '=====================================================11. END SELECT1. C ENDIF
  336. C 11 N=NXY (2) 11 N=NXY (2) KM=31.(N.EQ.2) KM=6
  337. DO 4 0 I1=NXY (KM), NXY (KM+1), NXY (KM+2) DO 35 I2=NXY (9-KM), NXY (10-KM), NXY (11-KM)1 = 11 J=I21.(N.EQ.l) GOTO 21 = 121. J=I1
  338. RES (1)=HL1(I) RES (2)=HL2(J) DO L=1,3 EPS (L)=0 EKS (L)=0 ENDDO1. C JH=1
  339. C GOTO 11 C 10 WRITE (6, '(/)')1. C WRITE (6,*) X Y U V W Mx1. C * My H'
  340. C WRITE (6,*) '====================================================.1. C 11 DO L=1,31. C EPS (L)=01. C EKS (L)=01. C ENDDO1. DO KF=KU, KW1. RES (KF+2)=AR (KF, J, I)1. ENDDO
  341. CALL COEF (RES (1), RES (2), CF, PA, NGEO)1. DO L=1,121. CF1(L)=CF (L)1. CF2(L)=CF (L)1. ENDDO
  342. AB=CF (1)*CF (2) SELECT CASE (NTP) CASE (1)
  343. EKS (1)=RES (3)*CF (9)/CF (1)+RES (4)*CF (4)*CF (6)/AB EKS (2)=RES (3)*CF (7)*CF (3)/AB+RES (4)*CF (11)/CF (2) EKS (3)=RES (3)*CF (10)*CF (1)/CF (2)+RES (4)*CF (12)*CF (2)/CF (1) CASE (2)
  344. EPS (1)=CF (6)*RES (4)/AB EPS (2)=CF (7)*RES (3)/AB
  345. EPS (3)=-CF (6)*RES (3)/AB-CF (7)*RES (4)/AB CASE (3,4)
  346. EPS (1)=CF (6)*RES (4)/AB+RES (5)*CF (3)
  347. W (1)=AR (K, J, I-KI) W (2)=AR (K, J, I-KI+1) W (3)=AR (K, J, I-KI+2) CALL DIF (HX1,HX2,W, X, DFX) KJ=11.(K.NE.3.AND.J.EQ.0) KJ=0 HY1=H2(J-KJ+1) HY2=H2(J-KJ+2) Y=HY1*KJ1.(K.NE.3.AND.J.EQ.N2) THEN KJ=21. Y=HY1+HY2 ENDIF1. W (1)=AR (K, J-KJ, I)
  348. W (2)=AR (K, J-KJ+1,I) W (3)=AR (K, J-KJ+2,I) CALL DIF (HY1,HY2,W, Y, DFY) IF (K.NE.3)THEN1.(I.EQ.0.OR.I.EQ.N1) CALL COEF (HL1(I-KI+1), RES (2), CF1, PA, NGEO) IF (J.EQ.0.OR.J.EQ.N2) CALL COEF (RES (1), HL2{J-KJ+1), CF2, PA, NGEO) ENDIF
  349. CALL COEF (HL1(I + KI), HL2(J+KJ), CF, PA, NGEO) ENDIF
  350. CALL DIF (HY1,HY2,W2,Y, DFX)
  351. EKS (3)=EKS (3)-2*DFX (1)/(CF (1)*CF (2))1. END SELECT1. ENDDO! K=KU, KW1. (NTP.NE.1) THEN
  352. RES (6)=PA (6)*(EPS (1)+PA (2)*EPS (2)) RES (7)=PA (6)*(PA (2)*EPS (1)+EPS (2)) RES (8)=PA (6)*(1.-PA (2))*EPS (3) /2 RES (12) =RES (6) /PA (5) RES (13)=RES (7)/PA (5) RES (14)=RES (8)/PA (5) ENDIF1.(NTP.NE.2) THEN
  353. RES (9)=-PA (7)*(EKS (1)+PA (2)*EKS (2)) RES (10)=-PA (7)*(PA (2)*EKS (1)+EKS (2))
  354. RES (11)=PA (7)*(1.-PA (2))*EKS (3)/21. BH2=PA (5)* PA (5)/61. RES (15)=RES (9)/BH21. RES (16)=RES (10)/BH21. RES (17)=RES (11)/BH21. ENDIF1. SELECT CASE (JH)
  355. CASE (1) SELECT CASE (NPRZ) CASE (1)1. SELECT CASE (NTP) CASE (1)
  356. WRITE (6,100) RES (1), RES (2), RES (5), RES (9), RES (10), RES (11) CASE (2)
  357. WRITE (6,100) RES (1), RES (2), RES (3), RES (4), RES (6), RES (7), RES (8) CASE (3,4)
  358. WRITE (6,100) RES (1), RES (2), RES (3), RES (4), RES (5), RES (6), *RES (7), RES (8) END SELECT CASE (2) SELECT CASE (NTP) CASE (1)
  359. WRITE (6,100) RES (1), RES (2), RES (5), RES (15), RES (16), RES (17) CASE (2)
  360. WRITE (6,100) RES (1), RES (2), RES (3), RES (4), RES (12), RES (13), RES (14) CASE (3,4)
  361. WRITE (6,100) RES (1), RES (2), RES (3), RES (4), RES (5), RES (12) ,
  362. RES (13), RES (14) END SELECT1. END SELECT CASE (2)
  363. SELECT CASE (NPRZ) CASE (1)
  364. WRITE (6,100) RES (1), RES (2), RES (3), RES (4), RES (5), RES (9), RES (10), *RES (11) CASE (2)
  365. WRITE (6,100) RES (1), RES (2), RES (3), RES (4), RES (5), RES (15), RES (16), *RES (17) END SELECT END SELECT JH=JH+11.(JH.EQ.2.AND.NTP.GE.3) GOTO 101. CONTINUE IF (NTP.GE.3)THEN
  366. WRITE (6,*) '---------------------------------------------------1. ELSE IF (NTP.EQ.2) THEN
  367. WRITE (6,*) '---------------------------------------------------1. ELSE
  368. WRITE (6, *) '---------------------------------------------------i1. ENDIF1. CONTINUE JH=JH+11.(JH.EQ.2.AND.NTP.GE.3) GOTO 101. MODULE GEOMETRY
  369. CF (2)=PA (8) CF (4)=1./PA (8) RETURN5 CF (1)=PA (8)1. CF (2)=PA (8)* DCOS (X)1. C F (3) = 1. / P A (8)1. CF (4)=CF (3)1. CF (7)=-CF (1)*DSIN (X)1. RETURN
  370. ARN=PA (8) c AL1=PA (9) c CE=PA (10) c AL2=PA (11)
  371. AMU=PA (9) CE=PA (10) AKN=PA (11) c DCS2=DCOS (X)*DCOS (X)c DSN2=DSIN (X)*DSIN (X)
  372. ПРЯМОУГОЛЬНАЯ система! координат
  373. ПОЛЯРНАЯ система координат
  374. ЦИЛИНДРИЧЕСКИЕ оболочки !{ Цилиндр с радиусом R=PA (8) }
  375. Оболочки вращения в ! ЦИЛИНДРИЧЕСКОЙ системе координат ! { Цилиндр с радиусом R=PA (8) }
  376. Оболочки вращения в ! СФЕРИЧЕСКОЙ системе координат !{ Сфера с радиусом В=РА (8) }
  377. ОБОЛОЧКИ СЛОЖНОЙ ГЕОМЕТРИИ
  378. Выше приводятся основные параметры рассматриваемой оболочки Например: Начальный радиус, количество волн, амплитуда волнистых оболочек и т. д.
  379. А9=-1*1*(3*1*12*1*11 + 4*:Е*12*К+1*2*1*)+ 3*1*Б*1*1*11*Б*Б*11-Е**Б*111*31С2 АЮ=02*0<1*(К12-Е*2*:Е*2*Е2-Е**:Е*11*0)-02.1*0*(2*ЗЮ2-Е**0*Е*1*11) А121= ((2*316(2*1*1*1−1*2*31 601) /ЭЮг) (01*01) А122=((01*01−2*0*011) / (01*01*01)) *Е*2/ЗЮ<2
  380. Матрица геметрических коэффициентов, записана в общем случае для ! каналовых поверхностей Иоахимсталя!
  381. СБ (1) =АЕШ*02 * ЭЮО/О СБ (2)=АЕ*Ы*О1*Е0"/0
  382. CF (3)=1.*(2.*SIG2-R*RRII* D)/(ARN * G 2 * SIG 2 * SIGQ)
  383. CF (4)=2. *R2/(ARN*G1*SIGQ)
  384. CF (5)=ARN*G2 *(D*SIGQI-SIGQ*DI)/D2
  385. С F (6)=ARN * SIGQ *(D*G2J-G2*DJ)/D2
  386. CF (7)=ARN* FJ*(D*G1I-G1*DI)/D2
  387. CF (8)=ARN*G1*(D* FJJ-FJ*DJ)/D2
  388. CF (9)=A9*D/(ARN*G2 * SIG2 *SIG2 * SIGQ)
  389. CF (10)=2.*A10/(ARN*ARN*SIG2*SIG2*G2*G2*G2)1. CF (11)=0
  390. CF (12)=2.*(A121+A122)/(ARN*ARN*FJ) RETURN END SUBROUTINE COEF
  391. ПРОЦЕДУА УЧЕТА СОБСТВЕННОГО ВЕСА КОНСТРУКЦИИ
  392. SUBROUTINE GRAVITY (X, Y, GRIV, GXYZ) IMPLICIT REAL*8 (A-H, 0-Z) DIMENSION GRIV (3), GXYZ (3)
  393. COMMON NTP, NGEO, N1, N2,KU, KW, KUS, KN, LS, NT1, NT2,NTR, PA (12)1. SELECT CASE (NGEO)
  394. NGEO=l ПРЯМОУГОЛЬНАЯ СИСТЕМА КООРДИНАТ1. CASE (1) DO K=l, 31. GXYZ (K)=GRIV (K) ENDDO
  395. NGEO=2 ПОЛЯРНАЯ СИСТЕМА КООРДИНАТ1. CASE (2)
  396. GXYZ (1) =DCOS (Y)*GRIV (1)+DSIN (Y)*GRIV (2) GXYZ (2)=-DSIN (Y)*GRIV (1)+DCOS (Y)*GRIV (2) GXYZ (3)=GRIV (3)
  397. NGEO=3 ЦИЛИНДРИЧЕСКАЯ ОБОЛОЧКА1. CASE (3)1. GXYZ (1)=GRIV (1)
  398. GXYZ (2)=-DSIN (Y/PA (8))*GRIV (2)+DCOS (Y/PA (8))*GRIV (3) GXYZ (3)=DCOS (Y/PA (8))*GRIV (2)+DSIN (Y/PA (8))*GRIV (3)
  399. NGEO=4 ЦИЛИНДРИЧЕСКАЯ СИСТЕМА КООРДИНАТ1. CASE (4)1. GXYZ (1)=GRIV (3)
  400. GXYZ (2)=-DSIN (Y)*GRIV (1)+DCOS (Y)*GRIV (2) GXYZ (3)=DCOS (Y)*GRIV (1)+DSIN (Y)*GRIV (2)
  401. NGEO=5 СФЕРИЧЕСКАЯ СИСТЕМА КООРДИНАТ1. CASE (5)
  402. NGEO=6 ОБОЛОЧКИ СЛОЖНОЙ ГЕОМЕТРИИ (КАНАЛОВЫЕ ПОВЕРХНОСТИ ИОАХИМСТАЛЯ)1. CASE (6)с ARN=PA (8)с AL1=PA (9)с СЕ=РА (10)с AL2=PA (11)
  403. С RI1=(DCS2 + 2 * DSN2)/(DCS2*DCOS (X))
  404. С RIII=DSIN (X)*(5*DCS2+6*DSN2)/(DCS2*DCS2)
  405. R=1+AMU*(1+DCOS (AKN*X)) RI=-AMU*AKN*DSIN (AKN*X) RII=-AMU*AKN*AKN*DCOS (AKN*X) RIII=AMU*AKN*AKN*AKN* DSIN (AKN*X)1. С R=2*(l+AMU*DCOS (X))1. С RI=-2*AMU*DSIN (X)1. С RII=-2*AMU*DCOS (X)1. С RI11 = 2 *AMU*DSIN (X)
  406. С F=0.75*DTAN (Y) С FJ=0.75/(DCOS (Y)*DCOS (Y))
  407. С FJJ=1.5* DTAN (Y)/(DCOS (Y)*DCOS (Y))1. F=DTAN (Y)
  408. FJ=1./(DCOS (Y)*DCOS (Y)) FJJ=2.* DTAN (Y)/(DCOS (Y)*DCOS (Y))
  409. МЕТОДИКА ВВОДА ИСХОДНЫХ ДАННЫХ1. И ВЫВОДА РЕЗУЛЬТАТОВ
  410. Ввод необходимых исходных данных в файл «INput.dat» осуществляется построчно, строго в следующем порядке:1. ЫТР ЫСЕО N1 N2 Пплоск.1. Е V XI о Х20 Н
  411. РА08=' 0 'РАО9=' 0 'РА10=' О 'РА11=' О 'РА12=' О1. Л1 3−2 А13 211 12 А1 11 32 А1 и v VI -2 0 0 0 0 0 0 0 1 11 12 11 12 0 0 Кх1 Кх2 Кш
  412. ПНапр 11 1 г А1 11 12 Кх1 Кх2 Кш3 11 12 А1 11 1 2 А1 Кх1 Кх2 Кш-2 0 0 0 0 0 0 0 0 1 11 12 11 12 Ндейст 0 <7о дх1 4x2
  413. ПНапр 11 1 г А1 11 12 Ндейст Ро Рпр
  414. ПНапр 11 1 г А1 11 12 0 м0 Мпр4 11 12 А1 11 12 А1 Рх1 Рх2 Рт5 11 12 А1 11 12 А1 0 мх1 Мх2б Нх ну Нг 0 0 0 до 0 0−2 0 0 0 0 0 0 0 0 0
  415. Пус ПНапр 11 1 г А1 11 12 А1−2 0 0 0 0 0 0 0 ввоаитыЕ паралпетры
  416. ОБЩИЕ ДАННЫЕ ПО ГЕОМЕТРИИ И МЕХАНИКЕ МА ТЕРИАЛАКОНСТРУКЦИИ
  417. ШТ номер типа задачи: ЫТР=1 — изгиб пластин, ШТ=2 — плоская задача теории упругости, ЫТР=3 — отсеки (секторы) оболочки.
  418. N1, N2 — число шагов сетки в направлении осей X/ и х2 соответственно.
  419. Пплоск показатель вида плоской задачи теории упругости при ЫТР=2: ПГщоск= 0 — задача о плоском напряжении, Пплоск- 1 — задача о плоской деформации. Значение Пплоск всегда принимает нулевое значение, когда ЫТР^ 2.
  420. Е модуль упругости материала конструкции.
  421. V коэффициент Пуассона материала.
  422. XIо и Х20 смещение сетки по координатным осям х}, х2 соответственно (координаты начала сетки, которая смещена по осям х-, х2 на величины XI0 и Х20 соответственно).1. Н толщина конструкции.
  423. РАО8- О 'РА09- О 'РАЮ- О 'РА11='0 'РА 12='0 эта строка используется при вводе дополнительных коэффициентов, используемых для конкретной конструкции, например: радиус сферы, радиус цилиндра и т. д.1. РАЗБИВКА СЕТКИ
  424. А/ разбивка сетки в направлении оси X/ от номера сетки до номера сетки /2 с шагом А/.71 72 4/ ~~ разбивка сетки в направлении оси х2 от номера сетки 7−1 до номера сетки у2 с шагом А/.
  425. Данные построчно вводятся необходимое число раз, до тех пор пока не будут указаны все закрепления узлов. Номера закрепленных узлов могут повторяться в разных строках для удобства организации ввода. Шаги А/, А/ -не могут равняться нулю.
  426. При жестком защемлении краев, помимо контурных узлов закрепляются законтурные узлы с координатами по х, (-1, N1+1), и по (-1, N2+1).-2 0 0000 0 0 признак конца ввода данных граничных условий.
  427. Строка вводится после окончания ввода всех данных о закреплении узлов.
  428. УЧЕТ УПРУГОВО ОСНОВАНИЯ И УПРУГИХ ОПОР
  429. При упругом опирании конструкции возможны следующие виды:
  430. Упругое опирание вдоль линии.
  431. Пнапр ц ь Л/ у'. 72 Кх, Кх2 К т
  432. Точечное (Узловое) упругое опирание.3 /. /2 А/ 71 72 А/ Кх, Кх2 Кт
  433. Если упругие опоры работают только вдоль одной из осей х, х2, т, то остальные коэффициенты упругости равняются нулю.-2 0 0000 0 0 показатель конца ввода данных об упругих опираниях.1. НАГРУЗКИ
  434. Рассмотрим записи этих нагрузок в файле исходных данных: 1. Линейно распределенная нагрузка по поверхности.
  435. Линейно распределенная нагрузка вдоль линии.
  436. Пнапр ?1 /2 А/. /2 Ндейст р0 рпр
  437. Узловые сосредоточенные нагрузки.4 /1 12 А/ 71 72 А/ Рх! Рх2 Рт
  438. Узловые сосредоточенные изгибающие моменты.5 ц к А/ 71 72 А/ 0 Мх1 Мх2
  439. Если собственный вес направлен противоположно одной из осей глобальной системы, то go принимается со знаком (-).-2 0 0000 0000 показатель конца ввода информации о нагрузках.
  440. ВЫВОД РЕЗУЛЬТАТОВ РАСЧЕТА НА ПЕЧАТЬ
  441. Пус показатель вида выводимых результатов: Пус=1 — в усилиях, Пус=2- в напряжениях- Пнапр показывает направление вывода результатов:
  442. Этой строкой завершается файл с исходными данными, и является завершающей информацией работы программы. запусн програтты
  443. По окончании работы программы результаты расчета могут быть просмотрены в файле результатов расчета и при необходимости выданы на печать.202
Заполнить форму текущей работой

На отлично!