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

ОСНОВНЫЕ ОБОЗНАЧЕНИЯ.

ВВЕДЕНИЕ.

ГЛАВА 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. Проведены расчеты оболочек в форме каналовых поверхностей Иоахимсталя вариационно-разностным методом на различные виды нагрузок.

136

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. Мельников Н.П. Металлические конструкции за рубежом. -М: Изд-во литературы по строительству, 1971. 400 с.

75. Милинский В.И. Дифференциальная геометрия. Л., 1934.-332 с.

76. Михайленко В.Е., Ковалев С.Н. Конструирование форм современных архитектурных сооружений. Киев: Бущвельник, 1978. - 112 с.

77. Михайленко В.Е., Обухова B.C., Подгорный А.Л. Формообразование оболочек в архитектуре. Киев: Буд1вельник, 1972. - 206 с.

78. Михлин С.Г. Вариационные методы математической физики. -М.: Наука, 1970.-360 с.

79. Молчанов H.H. Численные методы решения некоторых задач теории упругости. Киев: Наукова думка, 1979. - 316 с.

80. Наср Юнее Ахмед Аббуши. Волнообразные купола // Строительная механика инженерных конструкций и сооружений: Межвузовский сборник научных трудов. М.: АСВ, 2002. - вып.11. - с.49-58.

81. Новожилов В.В. Краткий очерк развития теории оболочек в СССР // Исследования по теории пластин и оболочек. Казань.: Изд-во Казан. ун-та, 1970. - вып.6-7. - с.5-22.

82. Новожилов В.В. Новый метод расчета тонких оболочек // Изв. АН СССР. ОТН. 1946. -№.1. -с.35-48.

83. Новожилов В.В. Теория тонких оболочек. Л.: Судостроение, 1962.-431 с.

84. Новожилов В.В., Финкельштейн P.M. О погрешности гипотезы Кирхгофа в теории оболочек //ПММ. 1943. -т.7. - вып.5. - с.331-340.

85. Новожилов В.В., Черных К.Ф., Михайловский Е.И. Линейная теория тонких оболочек. Л.: Политехника, 1991. - 656 с. ил.

86. Норден А.П. Теория поверхностей. М.: ГИТТЛ, 1956. - 260 с.

87. Павелайнен В.Я. Безмоментное напряженное состояние торойдального покрытия // Исследования по теории упругости и пластичности. -Л.: ЛГУ, 1961.-№1,-с. 110-118.

88. Паймушин В.Н. Некоторые задачи статики незамкнутых оболочек сложной формы и об одном методе их численного решения // Вопросы расчета прочности конструкций летательных аппаратов. Казань: Изд-во Казанского авиационного института, 1979. - с.67-76.

89. Паймушин В.Н., Андреев С. В. К численному исследованию напряженно-деформированного состояния однослойных и трехслойных пластин и оболочек сложной геометрии // Прикл. мех. Киев, 1983, - т.7. - с.24-30.

90. Плотников Ф. А. Расчет железобетонных спиральных камер водяных турбин круглого сечения. // Применение железобетона в машиностроении: Сб. статей. -М.: Машиностроение, 1964. с. 440-457.

91. Рекач В.Г. Развитие некоторых разделов теории расчета тонких оболочек // Труды УДН «Строительство»: Строительная механика. М.: УДН, 1967. - t.XXVIII. - вып.З. - с.3-5.

92. Рекач В.Г., Кривошапко С.Н. Расчет оболочек сложной геометрии. М.: Изд-во УДН, 1988. - 178 с.

93. Рекач В.Г., Рыжов H.H. Некоторые возможности расширения круга задач по конструированию и расчету оболочек // Труды УДН «Строительство»: Строительная механика. М.: УДН, 1970. - t.XLVIII. - вып.6.- с.3-8.

94. Розин Л.А. Задачи теории упругости и численные методы их решения. -СПб.: Изд-во СПбГУ, 1998. 532 с.

95. Рузлева Н.П. Кинематика образования циклической поверхностей // Труды УДН «Математика»: Прикладная геометрия. М.: УДН, 1967.- t.XXVI. вып.З. - с. 100-104.

96. Рыжов H.H. Каркасная теория задания и конструирования поверхностей // Труды УДН «Математика»: Прикладная геометрия. М.: УДН, 1967. - т. XXVI. - вып.З. - с. 3-16.

97. Рыжов H.H., Гершман И.П., Осипов В.А. Прикладная геометрия поверхностей // Тр. Моск. Научн.-метод, семинара по начертательной геометрии и инженерной графике. -М.: МАИ, 1972. ч.З. - вып.242. - с. 57-91.

98. Сальман Аль-Духейсат. К вопросу расчета без моментной оболочки в виде отсека циклиды Дюпена // Исследования по расчету элементов пространственных систем. М.: УДН, 1987. - с. 106-110.

99. Самуль В.И. Основы теории упругости и пластичности. М.: Высш. гик., 1982.-264 с.

100. Санчкс-Аркас М. Оболочки. М.: Изд-во литературы по строительству, 1964. -172 с.

101. Секулович М. Метод конечных элементов: Пер. с серб. М.: Стройиздат, 1993. - 664 с. ил.

102. Скидан И.А. Аналитическая теория формообразования оболочек // Труды международной научной конференции «0болочки-2001». М: Изд-во РУДН, 2001. - с.366-371.

103. Скидан И.А. Обобщенные цилиндрические координаты и их применение в прикладной геометрии // Прикладная геометрия и инженерная графика, Киев, 1971. - вып. 3. - с. 112-118.

104. Слицкоухов Ю.В., Буданов В.Д., Гапоев М.М. и др. Конструкции из дерева и пластмасс. -М.: Стройиздат, 1986. 543 с. ил.

105. Стеблянко В.Т. Об одном методе задания частного вида эпитро-хоидальных поверхностей // Прикладная геометрия и инженерная графика. -Киев, 1975. -вып.20. -с.89-91.

106. Сыроежкин М.Н. Железобетонные камеры кольцевого сечения // Энергомашиностроение, М., 1959. - №10. - с.26-32.

107. Тимошенко С.П., Войновский-Кригер С. Пластины и оболочки. М.: Физматгиз, 1963. - 635 с.

108. Фортран 90. Международный стандарт: Пер. с англ. -М.: Финансы и статистика, 1998. 378 с.

109. Хемминг Р.В. Численные методы для научных работников и инженеров: Пер. с англ. М.: Наука, 1972. - 400 с.

110. Церпа В. Исследование работы в области теории оболочек в ФРГ // Второй Международный конгресс по тонкостенным оболочкам-покрытиям. -Осло, 1957.-М., 1960.

111. Чаттерджи Б.К. Некоторые оболочки, возведенные в Индии // Большепролетные оболочки: Международный конгресс в Ленинграде. М.: Стройиздат, 1969. - т. 1. - с. 207-220.

112. Черных К.Ф. Линейная теория оболочек. Л.: Изд-во ЛГУ, 1962. -Т.1.-374 е.; 1964.-Т.2.-395 с.

113. Шуликовский В.И. Классическая дифференциальная геометрия.- М.: ГИФМЛ, 1963. 540 с.

114. Эльсгольц Л.Э. Дифференциальные уравнения и вариационное исчисление. -М: Наука, 1965. 279 с.

115. Якубовский A.M. Исследование аналитического метода задания циклид Дюпена при выделении их из конгруэнции окружностей // Прикладная геометрия. М.: УДН, 1971. - вып.4. - с.26-40.

116. Якубовский A.M. Циклические каркасы из линий кривизны // Труды УДН «Математика»: Прикладная геометрия. М.: УДН, 1967.- т.XXVI. вып.З. - с.66-78.

117. Якупов Н.М. Об одном методе расчета оболочек сложной геометрии // Исслед. по теории оболочек: Тр. семинара. Казань: Казанский физико-технический институт, 1984. - вып. 17. - ч.2. - с.4-17.

118. Якупов Н.М. О некоторых работах по расчету оболочек сложной геометрии // Тр. семинара. Казань: Казанский физико-технический институт, 1990. - №25. - с.86-94.

119. Якупов Н.М. Фрагменты оболочек сложной геометрии в тороидальной системе координат // Исследования по теории оболочек. Труды семинара. Казань: Казанский физико-технический институт, 1988. - вып.21, -ч.1.-с. 130- 137.

120. Якупов Н.М., Галимов Ш.К., Хисматуллин Н.И. От каменных глыб к тонкостенным конструкциям. Казань: Изд-во «SOS», 2001. - 96 с.

121. Бранков Г.Й. Вълнообразни черупкови конструкции. София: Изд-во на Българската академия на науките, 1961. - 80 с.

122. Бранков Г.Й. Някои въпроси от моментната и безмоментната теория на черупките. Трудове на Научноизследователския строителен институт София, 1959. - №2.

123. Argyris J.H. Matrix methods of structural analysis // Proc. 14-th meeting of AGARD. AGARDograph. 1962. - 72.

124. 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.

125. 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.

126. Caley A. On The Cyclide // Quarterly Journal of pure and applied mathematics. 1873. -№12. - p. 148.

127. Darboux G. System ortogonoux et les coordonees curilinges. 2nd adn. Paris, Guatvier-Villars. -1910.

128. Dean W.R. The distortion of curved tube due to internal pressure. Phil. Mag VII ser. Vol. 28. - 1939.

129. Dixon R. Asymmetric shells a new approach // Bulletin of international Association for Shell and Spatial Structures. -1991. - Vol. 32. -№3. - p. 133-137.

130. Doehlemann K. Geometrische trasformation, II, Leipzig. 1908.

131. Dupin Ch. Applications de geometry et de mechanique, Paris, Bachelier. 1822.

132. Dupin Ch. Develepment de geometrie. Paris. 1813.

133. Forsyth A.R. Lectures on the differential geometry of curves and surfaces. Cambridge.-1920

134. 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.

135. Maxwell J.C. On the Cyclide // Quarterly Journal of pure and applied mathematics, -1868. №9.

136. Nelson R.L. Stresses in shell structures // J. Sound and Vibr. 1981. -79. -N3. -p.397-414.

137. 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.

138. Jiirgen Joedicke. Shell Architecture. London: Alec Tiranti Ltd., 1962.-304 c.152

139. КАНАЛОВЫЕ ПОВЕРХНОСТИ ИОАХИМСТАЛЯ С РАЗЛИЧНЫМИ ОЧЕРТАНИЯМИ В ПЛАНЕ1. г2{а) = 1 + ¡исо$>{а) эпитрохоидальная поверхность• С2= 0, //-0.5.155• С2 =0, а = \, ju= 0.8, к = 2.25.156

140. ОБОЛОЧКИ НА ОСНОВЕ КАНАЛОВЫХ ПОВЕРХНОСТЕЙ ИОАХИМСТАЛЯ, ПРОЕКТИРУЕМЫЕ В ДИПЛОМНЫХ РАБОТАХ СТУДЕНТОВ СТРОИТЕЛЕЙ РУДН

141. ВОЗМОЖНОСТЬ КОНСТРУИРОВАНИЯ КОРПУСА ПРИТОЧНОЙ НАСАДКИ В ФОРМЕ КАНАЛОВЫХ ПОВЕРХНОСТЕЙ1. ИОАХИМСТАЛЯ

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

143. Проточные насадки или воздухораспределители. I) статический (улиточный) центробежный конструкции ВНИИГС. II) конструирование корпуса в форме каналовых поверхностей Иоахимсталя. а приточный, б - приточно-вытяжной.

144. ВЛИЯНИЕ ПАРАМЕТРА С НА ВИД КАНАЛОВОЙ ПОВЕРХНОСТИ

145. ИОАХИМСТАЛЯ ПРИ г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

146. РАСЧЕТ ТОНКОСТЕННЫХ КОНСТРУКЦИЙ ВАРИАЦИОННО-РАЗНОСТНЫМ МЕТОДОМ Авторы: проф. ИВАНОВ В.Н.

147. НАСР ЮНЕС (Палестина) ПРОГРАММА ПОДГОТОВЛЕНА НА КАФЕДРЕ СОПРОТИВЛЕНИЯ МАТЕРИАЛОВ РОССИЙСКОГО УНИВЕРСИТЕТА ДРУЖБЫ НАРОДОВ МОСКВА, 2000 г.

148. М 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

149. PROGRAM SHELLVRM USE GEOMETRY1.PLICIT REAL*8 (A-H,0-Z)

150. REAL,PARAMETER :: PI=3.1415926535897932384626433832795 DIMENSION AH(30000000),MN(6),MH(300000),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

151. PRINT *, ' Последнее обновление программы 00:43 19.09.2001' PRINT*

152. PRINT*,'ВВОДИТЕ ИМЕНА ФАЙЛОВ ИСХОДНЫХ ДАННЫХ И РЕЗУЛЬТАТОВ' READ(*,*) NMR,NMW !ввод имен файлов

153. OPEN(UNIТ=5,FILE=NMR,STATUS='OLD 1) !c исходными данными

154. OPEN(UNIT=6,FILE=NMW,STATUS='NEW') !и результатами.

155. WRITE(*,*) ' „АТА' READ(*,*) СDAT WRITE(6,1000) СDAT

156. WRIТЕ(б, * ) ' !!!!!!!!!!!!!!!!!!!!!!!!!!!! 1 !!!!!! 1 1 !!!!!!!! 1 1 1 1 !!! 1

157. M I I | I I | I I | | | M ¡I | I i

158. WRITE(б,*)' !! РАСЧЕТ ТОНКОСТЕННЫХ КОНСТРУКЦИЙ ВАРИАЦИОННО-РАЗН *ОСТНЫМ МЕТОДОМ !!'

159. WRITE(б,*)' !! Авторы: проф. ИВАНОВ В.Н., НАСР ЮНЕС (Палестина)

160. 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

161. X | I I I I I I I И II I I I I I

162. WRITE(6,1010) NMR С NTP-номер типа задачи: 1-ИЗГИБ ПЛАСТИН,

163. С 2-ПЛОСКАЯ ЗАДАЧА ТЕОРИИ УПРУГОСТИ,1. С 3-ОТСЕКИ ОБОЛОЧКИ,1. С 4-ЗАМКНУТЫЕ ОБОЛОЧКИ.

164. С NGEO-номер геометрии задачи: 1-ПРЯМОУГОЛЬНАЯ СИСТЕМА КООРДИНАТ,

165. С 2-ПОЛЯРНАЯ СИСТЕМА КООРДИНАТ,

166. С 3-ЦИЛИНДРИЧЕСКИЕ ИЛИ ПОЛОГИЕ ОБОЛОЧКИ,1. С 4-ОБОЛОЧКИ ВРАЩЕНИЯ,1. С 5-СФЕРИЧЕСКАЯ ОБОЛОЧКА,

167. С 6-КАНАЛОВЫЕ ПОВЕРХНОСТИ ИОАХИМСТАЛЯ

168. 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

169. WRITE(6 *)' ***************************************************** ******************** '1. ПЛОСЕАЯ ДЕФОРМАЦИЯ

170. РАСЧЕТ ОТСЕКОВ (СЕКТОРОВ) ОБОЛОЧКИ РАСЧЕТ ЗАМКНУТОЙ ОБОЛОЧКИ1. SELECT CASE(NTP) CASE(1)

171. WRITE(б,1018) ' РАСЧЕТ ПЛАСТИНКИ НА ИЗГИБ1. CASE(2)

172. WRITE(6,1018) ' РАСЧЕТ ПЛОСКОЙ ЗАДАЧИ ТЕОРИИ УПРУГОСТИ1. SELECT CASE(KL) CASE(0)

173. WRITE(б,1017) ' ПЛОСКОЕ НАПРЯЖЕННОЕ СОСТОЯНИЕ

174. CASE(1) WRITE(б,1017) END SELECT CASE(3) WRITE(6,1018) CASE(4) WRITE(6,1018) CASE DEFAULT

175. STOP 'HE ПРАВИЛЬНЫЙ ВВОД НОМЕРА ТИПА ЗАДАЧИ' END SELECT SELECT CASE (NGEO) CASE(1)

176. WRITE(6,1019) ' ПРЯМОУГОЛЬНАЯ СИСТЕМА КООРДИНАТ

177. 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)

178. ПОЛЯРНАЯ СИСТЕМА КООРДИНАТ ЦИЛИНДРИЧЕСКАЯ (ИЛИ ПОЛОГАЯ) ОБОЛОЧКА ОБОЛОЧКА ВРАЩЕНИЯ СФЕРИЧЕСКАЯ ОБОЛОЧКА КАНАЛОВАЯ ПОВЕРХНОСТЬ ИОАХИМСТАЛЯ1. CASE DEFAULT

179. STOP 'НЕ ПРАВИЛЬНЫЙ ВВОД НОМЕРА ГЕОМЕТРИИ ЗАДАЧИ' END SELECT

180. WRITE ( 6 ' * ^ * * * * ^ * * ^ * * * * * * * * * ^ * * ^ ^ ^ к -к к к к к к к к к к к к к к -к к к к к * к к к кк к *к-ккккккккккккккк к к~кк I

181. NLl=NHl+Nl+2 NH2=NL1+N1+1 NL2=NH2+N2+2 NAS=NL2+N2+1 NMNH=11. AH(NL1)=PA(3)

182. KUS= (N1 + 3)* (N2 + 3)+NTR1. D=',PA(7)

183. Адрес (@) каждого массива в едином массиве АН ЫН1-0 Н1(м.шагов сетки по X) NLl-@ НЬ1(м.длин(коорд.)по X) NH2-@ Н2(м.шагов сетки по У) NL2-@ HL2(м.длин(коорд.)по Y)количество общих узловыхг

184. WRITE(6,*) 'KUS=',KUS WRITE(6,1040) DLSUX-0

185. READ(5,*) I, J,A !c номера сетки I no J

186. WRITE(6,1070) 60 READ(5,*) MN, MZ

187. WRITE(6,*) 'KN=',KN С IK=11. С DO 120 IW=-1,N1+1

188. С WRITE(6,82) (MH(К),K=IK,IK+(N2+3)*NTR-1)

189. С 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

190. CALL SHLLESF(AH(NH1) ,AH(NH2),AH(NL1) , AH(NL2) ,MH(1) ,AH(NAS) ) WRITE(*,*) 'esf'

191. CALL SHLLOAD(AH(NH1),AH(NL1) , AH(NH2) , AH(NL2) , MH(1) ,AH(NAR) ) WRITE(*,*) 'load' NP=0 NT = 1

192. CALL DSLT(AH(NAS),AH(NAR),AH(NAP),KN,LS,NT) WRITE(*,*) 'dslt' С WRITE(6,*) (AH(I),I=NAR,NAR+KN-1)

193. CALL SHLLSORT(AH(NAR),AH(NAP),MH(1)) WRITE(*,*) 'sort'

194. CALL SHLLRES(AH(NAP),AH(NH1),AH(NH2),AH(NL1),AH(NL2)) WRITE(*,*) 'result'

195. SUBROUTINE SHLLESF(HI,Н2,HL1,HL2, MNH,AS) USE GEOMETRY1.PLICIT REAL*8 (A-H,0-Z)

196. DIMENSION HI(0:N1+1),H2(0:N2+1),HL1(0:N1),HL2(0:N2)

197. DIMENSION AS(1),DF(6,9,4),HN(3,3),DS(4),DO(6,9,4)

198. DIMENSION CF(12 ) , EG(3, 6) , AP(200)

199. DIMENSION MNH(KU:KW,-1:N2+1,-1:N1+1)

200. DIMENSION RM(9,9),RM1(9,9),KI(4) , KPSL(7),PAS(3)

201. DIMENSION EGN(3,6),EGT(6,6),EGDF(6,9)

202. COMMON NTP, NGEO,N1,N2,KU,KW,KUS,KN,LS,NT1,NT2,NTR,PA (12) WRITE(6,*) 'PA=' WRITE(6,*) PA

203. 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

204. DO 200 1=0,N1 DO 200 J=0, N2 X = HL1 (I) Y=HL2(J)

205. CALL COEF(X,Y,CF,PA,NGEO) AB=CF(1)*CF(2)/4 C WRITE(*,*) 'X=',X, ' Y= ' , Y

206. 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

207. C IF (I.EQ.1.AND.J.EQ.1) THEN1. C WRITE(6,*)'DS='1. C WRITE(6,*) DS1. C ENDIF1. C WRITE(*,*) PA

208. C WRITE(*,*) AO, Al, A2, A3

209. C WRITE(*,*) BO, Bl, B2, B3

210. C WRITE(*,*) A1B1, A1B2, A2B1, A2B21. C END IF1. DO 8 K=1,4

211. DF (4 2, K) =DF(4,2,K)*A3/(1+A0)

212. DF (4 5, K) =DF(4,5,K)*A3/A0

213. DF ( 4 8, K) =DF(4,8,K)*A3/(AO*(1+A0))

214. DF ( 6 4, K) =DF(6,4,K)*B3/(1+B0)

215. DF ( 6 5, K) =DF(6,5,K)*B3/B0

216. C 301 WRITE(6,302) (DF(L1,L3,L),L3=1,9)302 FORMAT(9E9.2) C END IF

217. DO 16 K=KU,KW DO 15 L=KU,KW

218. DO 11 11=1,9 DO 11 Jl=l,9 11 RM(II,J1)=01. DO 12 NT=NT1,NT2

219. CALL EGFMATRIX(CF,EG,NT,L)

220. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1) THEN

221. C WRITE(6,*) 'L=',L, ' NT=',NT, ' EG', (NT-1)* 3 + L1. C DO 11=1,3

222. C WRITE(6,302)(EG(II,Jl),Jl=l,6)1. C ENDDO1. C ENDIF

223. CALL MULTMATRIX(HN(1,1),EG(1,1),AP(1),3,3,3,6)

224. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1) THEN

225. 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

226. C WRITE(6,302) (EGN(IP,JP),JP=1, 6)1. C ENDDO1. C ENDIF

227. CALL EGFMATRIX(CF, EG, NT, K)

228. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1) THEN

229. C WRITE(6,*)'K=',K,' NT=',NT,' EG',(NT-1)*3+K1. C DO 11=1,3

230. C WRITE(6,302)(EG(II,Jl),Jl=l,6)1. C ENDDO1. C ENDIF

231. CALL MULTMATRIX(EG(1,1),AP(1),AP(19),3,6,3,6)

232. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1) THEN

233. 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

234. C WRITE(6,302)(EGT(IP,JP),JP=1,6)1. C ENDDO IIP C ENDIF

235. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1) THEN

236. C WRITE(6,*)'L=',L,' NT=',NT,1 EGT*NEG',(NT-1)*3+L1. C KI = 1 C DO IW=1,6

237. C WRITE ( 6,302) (AP(JI) ,JI=KI,KI+5)1. C KI=KI+61. C ENDDO1. C ENDIF1. DO 12 IDF=1,4

238. CALL MULTMATRIX(AP(19),DF(1,1,IDF),AP(55),6,6,6,9)

239. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1.AND.L.EQ.2.AND.NT.EQ.1) THEN

240. 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

241. C WRITE(6,302)(EGDF(IP,JP),JP=1,9)1. C ENDDO1. C ENDIF

242. CALL MULTMATRIX(DF(1,1,IDF),AP(55),RM1(1,1),6,9,6,9)

243. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1.AND.L.EQ.2.AND.NT.EQ.1) THEN

244. C WRITE(6,*)'L=',L,' NT=',NT,' RM1', IDF1. C DO IP=1,9

245. C WRITE(6,302)(RM1(IP,JP),JP=1,9)1. C ENDDO C ENDIF1. DO 11=1,9 DO Jl=l,9

246. RM(II,Jl)=RM(II,Jl)+RM1(II,Jl)*DS(IDF)*PA(NT+5)1. ENDDO1. ENDDO

247. C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1.AND.L.EQ.2.AND.NT.EQ.1) THEN

248. C WRITE(6,*)'L=',L,' NT=',NT,' RM',' IDF=',IDF1. C DO IP=1,9

249. C WRITE(6,302)(RM(IP,JP),JP=1,9)1. C ENDDO1. C ENDIF

250. CONTINUE C IF (I.EQ.1.AND.J.EQ.1.AND.K.EQ.1) THENc c c c c c

251. WRITE(6,*) 'K=',K, ' L= ' , L WRITE(6, *)'RM' DO IP=1,9

252. WRITE(6,302)(RM(IP,JP),JP=1,9) ENDDO ENDIF1.(K.NE.3.AND.L . NE.3) GOTO 130 LK0 = 111 = 0

253. 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

254. 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

255. 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

256. KI (3-LKO)=LK2+KJ0 KI(5-LKO)=LK2-KJ0 DO 100 Kl=l,9 KI(LK0)=K1 KI(LK0+2)=K1

257. 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 60130 L1=0

258. DO 14 11=1-1,1+1 DO 14 J1=J-1,J+1 L1=L1+11. NU=MNH(K,Jl,II)

259. C IF (I.EQ.1.AND.J.EQ.2.AND.L.EQ.1) THEN

260. C WRITE(6,*)'K=',K,' 11=',II,' J1=',J1,' NU=',NU1. C ENDIF1.(NU.LT.l) GOTO 14 L2 = 0

261. 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

262. С IF (I.EQ.l.AND.J.EQ.l) THEN1. С WRITE(6,*) 'RM='1. С DO 319 13=1,9

263. С 319 WRITE(6,302) (RM(13,14),14=1,9) С ENDIF401 CONTINUE15 CONTINUE16 CONTINUE 200 CONTINUE

264. ПРОЦЕДУРА УЧЕТА УПРУГОГО ОСНОВАНИЯ ПОД КОНСТРУКЦИЕЙ

265. СПЛОШНОЕ УПРУГОЕ ОСНОВАНИЕ' ' УПРУГИЕ ОПОРЫ ВДОЛЬ ЛИНИИ' ' УПРУГИЕ УЗЛОВЫЕ ОПОРЫ'1. DO I=KPSL(2),KPSL(3)

266. 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)

267. DO I=KPSL(3),KPSL(4),KPSL(5)

268. 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)

269. 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)

270. 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

271. FORMAT(15,5X,615,5X,3E11.3) C NH=11. C WRITE(6,*) 'AS='1. C DO NU=1,KN

272. 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

273. C 310 WRITE(6,302) (EG(L0,L1,L2,L3),L3=1,6)1. RETURN END1. SUBROUTINE MHHl(MNH)

274. DIMENSION MNH(KU:KW,-1:N2+1,-1:N1+1)

275. 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

276. 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)

277. DIMENSION MNH(KU:KW,-1:N2+1,-1:N1+1)

278. COMMON NTP,NGEO,N1,N2,KU,KW,KUS,KN,LS,NT1,NT2,NTR,PA(12) LS=1

279. 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

280. 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

281. AP (J) =AS (M+J-1) /AS (M) DO 3 K=I+1,I+LP IF(K.GT.NS) GOTO 4 M=K-I

282. K1=(I-1)*N+M+1 IF(AS(KI).EQ.0) GOTO 3 LI-(K-l)*N DO 2 J=1,N-M

283. 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

284. DO 9 K=I+1,I+LP IF(K.GT.NS) GOTO 10 M=M+1

285. 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

286. SUBROUTINE MULTMATRIX(A,B,C,M,N,N1,K) IMPLICIT REAL*8(A-H,O-Z)

287. 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

288. C(KI)=C(KI)+A(IA)*B(IB) RETURN1. END

289. SUBROUTINE EGFMATRIX(CF,EG, IP,K)1.PLICIT REAL*8(A-H,0-Z)

290. 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)

291. GOTO (1,2,3,4,5,6) (IP-1)*3+K

292. 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

293. EG(1,1)=CF(6)/A1A2 EG(2,3)=1./CF(2)

294. EG(3,1)=-CF(7)/A1A2 EG(3,2)=1./CF(1) RETURN

295. EG(1,1)=CF(3) EG(2,1)=CF(4) RETURN

296. EG(1,1)=CF(9)/CF(1) EG(1,2)=CF(3)/CF(1)

297. 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

298. EG(1,1)=CF(6)*CF(4)/A1A2 EG(2,1)=CF(11)/CF(2)

299. EG(2,3)=CF(4)/CF(2) EG(3,1)=CF(12)* CF(2)/CF(1) EG(3,2)=CF(4)/CF(1) RETURN

300. 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)

301. 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

302. SUBROUTINE SHLLOAD(HI,HL1,H2,HL2,MNH,AR) USE GEOMETRY1.PLICIT REAL*8 (A-H,O-Z)

303. DO 2 0 I=MH(2),MH(3) DO 2 0 J=MH(4),MH(5) X=HL1(I) Y=HL2(J)

304. HL1(MH(2) ) )+QL(3)*(Y-HL2(MH (4) ) ) )*DS30 KP=MH(2)1. DO I=MH(3) ,MH(4) ,MH (5)

305. 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)

306. 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)

307. 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)

308. 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)

309. 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)

310. CALL COEF(X,Y,CF,PA,NGEO) DS = 01. SELECT CASE(KP) CASE(1)

311. 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

312. 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)

313. 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)

314. 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

315. 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

316. DO 1=0,N1 DO J=0,N2 X=HL1 (I) Y=HL2(J)

317. FORMAT(15,5X,615,5X,3F10.2) RETURN END

318. SUBROUTINE SHLLSORT(AR,ARS,MNH) IMPLICIT REAL*8(A-H,O-Z)

319. 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

320. C WRITE(6,*) (ARS(N),N=K,K+NTR-1)1. C 21 K=K+NTR RETURN END

321. SUBROUTINE SHLLRES(AR,HI,H2,HLl,HL2)1. USE GEOMETRY1.PLICIT REAL*8(A-H,O-Z)

322. WRITE(6,*) 'X Y W Mx My H'

323. WRITE(6,*) '=====================================================11. CASE(2)

324. WRITE(6,*) ' X Y U V Nx Ny1. S'

325. WRITE(6,*) 1 ================================================== ===1. CASE DEFAULT1. WRITE(6,*) X Y U V W Nx1. Ny S '

326. WRITE(6,*) '=====================================================1. END SELECT1. CASE(2)1. SELECT CASE(NTP) CASE(1)

327. WRITE(6,*) 'X Y W SigMx SigMy TawXY

328. WRITE(6,*) '=====================================================1. CASE(2)

329. WRITE(6,*) X Y U V SigNx SigNy1. TawXY'

330. WRITE(6,*) '======================================================1. CASE DEFAULT

331. WRITE(6,*) X Y U V W SigNx1. SigNy TawXY'

332. WRITE(6,*) '======================================================

333. END SELECT END SELECT JH=1 GOTO 11 C 10 IF(NTP.GE.3) THEN10 WRITE(6, '(/)')

334. SELECT CASE (NPRZ) CASE(1)1. WRITE(6,*) X Y U V W Mx1. My H'

335. WRITE(6,*) '=====================================================1. CASE(2)

336. WRITE(6,*) X Y U V W SigMx1. SigMy TawXY'

337. WRITE(6,*) '=====================================================11. END SELECT1. C ENDIF

338. C 11 N=NXY(2) 11 N=NXY(2) KM=31.(N.EQ.2) KM=6

339. 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

340. RES(1)=HL1(I) RES(2)=HL2(J) DO L=1,3 EPS(L)=0 EKS(L)=0 ENDDO1. C JH=1

341. C GOTO 11 C 10 WRITE(6, '(/)')1. C WRITE(6,*) X Y U V W Mx1. C * My H'

342. 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

343. 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

344. AB=CF(1)*CF(2) SELECT CASE(NTP) CASE(1)

345. 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)

346. EPS(1)=CF(6)*RES(4)/AB EPS(2)=CF(7)*RES(3)/AB

347. EPS(3)=-CF(6)*RES(3)/AB-CF(7)*RES(4)/AB CASE(3,4)

348. EPS(1)=CF(6)*RES(4)/AB+RES(5)*CF(3)

349. 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)

350. 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

351. CALL COEF(HL1(I + KI) ,HL2(J+KJ) ,CF,PA,NGEO) ENDIF

352. CALL DIF(HY1,HY2,W2,Y,DFX)

353. EKS(3)=EKS(3)-2*DFX(1)/(CF(1)*CF(2))1. END SELECT1. ENDDO !K=KU,KW1. (NTP.NE.1) THEN

354. 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

355. RES(9)=-PA(7)*(EKS(1)+PA(2)*EKS(2) ) RES(10)=-PA(7)*(PA(2)*EKS(1)+EKS(2))

356. 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)

357. CASE(1) SELECT CASE (NPRZ) CASE(1)1. SELECT CASE(NTP) CASE(1)

358. WRITE(6,100) RES(1),RES (2) ,RES(5) ,RES (9) ,RES(10) ,RES(11) CASE(2)

359. WRITE(6,100) RES(1),RES(2) ,RES(3) ,RES(4) ,RES(6),RES (7) ,RES (8) CASE(3,4)

360. 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)

361. WRITE(6,100) RES(1) ,RES(2),RES(5) ,RES(15),RES (16) ,RES(17) CASE(2)

362. WRITE(6,100) RES(1),RES(2),RES(3),RES(4),RES(12),RES(13),RES(14) CASE(3,4)

363. WRITE(6,100) RES(1) ,RES(2) ,RES(3) ,RES(4) , RES(5) , RES(12) ,

364. RES(13) ,RES(14) END SELECT1. END SELECT CASE(2)

365. SELECT CASE (NPRZ) CASE(1)

366. WRITE(6,100) RES(1),RES(2),RES(3),RES(4),RES(5),RES(9),RES(10), *RES(11) CASE(2)

367. 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

368. WRITE(6,*) '---------------------------------------------------1. ELSE IF(NTP.EQ.2) THEN

369. WRITE(6,*) '---------------------------------------------------1. ELSE

370. WRITE(6, *) '---------------------------------------------------i1. ENDIF1. CONTINUE JH=JH+11.(JH.EQ.2.AND.NTP.GE.3) GOTO 101. MODULE GEOMETRY

371. 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

372. ARN=PA(8) c AL1=PA(9) c CE=PA(10) c AL2=PA(11)

373. AMU=PA(9) CE=PA(10) AKN=PA(11) c DCS2=DCOS(X)*DCOS(X)c DSN2=DSIN(X)*DSIN(X)

374. ПРЯМОУГОЛЬНАЯ система !координат

375. ПОЛЯРНАЯ система координат

376. ЦИЛИНДРИЧЕСКИЕ оболочки !{ Цилиндр с радиусом R=PA(8) }

377. Оболочки вращения в ! ЦИЛИНДРИЧЕСКОЙ системе координат ! { Цилиндр с радиусом R=PA(8) }

378. Оболочки вращения в ! СФЕРИЧЕСКОЙ системе координат !{ Сфера с радиусом В=РА(8) }

379. ОБОЛОЧКИ СЛОЖНОЙ ГЕОМЕТРИИ

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

381. А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*31601) /ЭЮг ) (01*01) А122=( (01*01-2*0*011) / (01*01*01) ) *Е*2/ЗЮ<2

382. Матрица геметрических коэффициентов, записана в общем случае для ! каналовых поверхностей Иоахимсталя!

383. СБ (1) =АЕШ*02 * ЭЮО/О СБ(2)=АЕ*Ы*О1*Е0"/0

384. CF(3)=1.*(2.*SIG2-R*RRII* D)/(ARN * G 2 * SIG 2 * SIGQ)

385. CF(4)=2. *R2/(ARN*G1*SIGQ)

386. CF(5)=ARN*G2 *(D*SIGQI-SIGQ*DI)/D2

387. С F(6)=ARN * SIGQ *(D*G2J-G2*DJ)/D2

388. CF(7)=ARN* FJ*(D*G1I-G1*DI)/D2

389. CF(8)=ARN*G1*(D* FJJ-FJ*DJ)/D2

390. CF(9)=A9*D/(ARN*G2 * SIG2 *SIG2 * SIGQ)

391. CF(10)=2.*A10/(ARN*ARN*SIG2*SIG2*G2*G2*G2)1. CF(11)=0

392. CF(12)=2.*(A121+A122)/(ARN*ARN*FJ) RETURN END SUBROUTINE COEF

393. ПРОЦЕДУА УЧЕТА СОБСТВЕННОГО ВЕСА КОНСТРУКЦИИ

394. SUBROUTINE GRAVITY(X,Y,GRIV,GXYZ) IMPLICIT REAL*8 (A-H,0-Z) DIMENSION GRIV(3),GXYZ(3)

395. COMMON NTP,NGEO,N1,N2,KU,KW,KUS,KN,LS,NT1,NT2,NTR,PA(12)1. SELECT CASE(NGEO)

396. NGEO=l ПРЯМОУГОЛЬНАЯ СИСТЕМА КООРДИНАТ1. CASE(1) DO K=l,31. GXYZ(K)=GRIV(K) ENDDO

397. NGEO=2 ПОЛЯРНАЯ СИСТЕМА КООРДИНАТ1. CASE(2)

398. GXYZ(1) =DCOS(Y)*GRIV(1)+DSIN(Y)*GRIV (2) GXYZ(2)=-DSIN(Y)*GRIV(1)+DCOS(Y)*GRIV(2) GXYZ(3)=GRIV(3)

399. NGEO=3 ЦИЛИНДРИЧЕСКАЯ ОБОЛОЧКА1. CASE(3)1. GXYZ(1)=GRIV(1)

400. 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)

401. NGEO=4 ЦИЛИНДРИЧЕСКАЯ СИСТЕМА КООРДИНАТ1. CASE(4)1. GXYZ(1)=GRIV(3)

402. GXYZ(2)=-DSIN(Y)*GRIV(1)+DCOS(Y)*GRIV(2) GXYZ(3)=DCOS(Y)*GRIV(1)+DSIN(Y)*GRIV(2)

403. NGEO=5 СФЕРИЧЕСКАЯ СИСТЕМА КООРДИНАТ1. CASE(5)

404. NGEO=6 ОБОЛОЧКИ СЛОЖНОЙ ГЕОМЕТРИИ (КАНАЛОВЫЕ ПОВЕРХНОСТИ ИОАХИМСТАЛЯ)1. CASE(6)с ARN=PA(8)с AL1=PA(9)с СЕ=РА(10)с AL2=PA(11)

405. С RI1=(DCS2 + 2 * DSN2)/(DCS2*DCOS(X) )

406. С RIII=DSIN(X)*(5*DCS2+6*DSN2)/(DCS2*DCS2)

407. 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)

408. С F=0.75*DTAN(Y) С FJ=0.75/(DCOS(Y)*DCOS(Y))

409. С FJJ=1.5* DTAN(Y)/(DCOS(Y)*DCOS(Y) )1. F=DTAN(Y)

410. FJ=1./(DCOS(Y)*DCOS(Y)) FJJ=2.* DTAN(Y)/(DCOS(Y)*DCOS(Y))

411. МЕТОДИКА ВВОДА ИСХОДНЫХ ДАННЫХ1. И ВЫВОДА РЕЗУЛЬТАТОВ

412. Ввод необходимых исходных данных в файл «INput.dat» осуществляется построчно, строго в следующем порядке:1. ЫТР ЫСЕО N1 N2 Пплоск.1. Е V XI о Х20 Н

413. РА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 Кш

414. ПНапр 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

415. ПНапр 11 1г А1 11 12 Ндейст Ро Рпр

416. ПНапр 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

417. Пус ПНапр 11 1г А1 11 12 А1-2 0 0 0 0 0 0 0 ввоаитыЕ паралпетры

418. ОБЩИЕ ДАННЫЕ ПО ГЕОМЕТРИИ И МЕХАНИКЕ МА ТЕРИАЛАКОНСТРУКЦИИ

419. ШТ номер типа задачи: ЫТР=1 - изгиб пластин, ШТ=2 - плоская задача теории упругости, ЫТР=3 - отсеки (секторы) оболочки.

420. N1, N2 — число шагов сетки в направлении осей X/ и х2 соответственно.

421. Пплоск показатель вида плоской задачи теории упругости при ЫТР=2: ПГщоск= 0 - задача о плоском напряжении, Пплоск- 1 - задача о плоской деформации. Значение Пплоск всегда принимает нулевое значение, когда ЫТР^ 2.

422. Е модуль упругости материала конструкции.

423. V коэффициент Пуассона материала.

424. XIо и Х20 смещение сетки по координатным осям х}, х2 соответственно (координаты начала сетки, которая смещена по осям х;, х2 на величины XI0 и Х20 соответственно).1. Н толщина конструкции.

425. РАО8- О 'РА09- О 'РАЮ- О 'РА11='0 'РА 12='0 эта строка используется при вводе дополнительных коэффициентов, используемых для конкретной конструкции, например: радиус сферы, радиус цилиндра и т.д.1. РАЗБИВКА СЕТКИ

426. А/ разбивка сетки в направлении оси X/ от номера сетки до номера сетки /2 с шагом А/.71 72 4/ ~~ разбивка сетки в направлении оси х2 от номера сетки 7-1 до номера сетки у2 с шагом А/.

427. Данные построчно вводятся необходимое число раз, до тех пор пока не будут указаны все закрепления узлов. Номера закрепленных узлов могут повторяться в разных строках для удобства организации ввода. Шаги А/, А/ -не могут равняться нулю.

428. При жестком защемлении краев, помимо контурных узлов закрепляются законтурные узлы с координатами по х, (-1, N1+1), и по (-1, N2+1).-2 0 0000 0 0 признак конца ввода данных граничных условий.

429. Строка вводится после окончания ввода всех данных о закреплении узлов.

430. УЧЕТ УПРУГОВО ОСНОВАНИЯ И УПРУГИХ ОПОР

431. При упругом опирании конструкции возможны следующие виды:

432. Упругое опирание вдоль линии.

433. Пнапр ц ь Л/ у'. 72 Кх, Кх2 К т

434. Точечное (Узловое) упругое опирание.3 /. /2 А/ 71 72 А/ Кх, Кх2 Кт

435. Если упругие опоры работают только вдоль одной из осей х\,х2, т , то остальные коэффициенты упругости равняются нулю.-2 0 0000 0 0 показатель конца ввода данных об упругих опираниях.1. НАГРУЗКИ

436. Рассмотрим записи этих нагрузок в файле исходных данных: 1. Линейно распределенная нагрузка по поверхности.

437. Линейно распределенная нагрузка вдоль линии.

438. Пнапр ¿1 /2 А/ . /2 Ндейст р0 рпр

439. Узловые сосредоточенные нагрузки.4 /1 12 А/ 71 72 А/ Рх! Рх2 Рт

440. Узловые сосредоточенные изгибающие моменты.5 ц к А/ 71 72 А/ 0 Мх1 Мх2

441. Если собственный вес направлен противоположно одной из осей глобальной системы, то go принимается со знаком (-).-2 0 0000 0000 показатель конца ввода информации о нагрузках.

442. ВЫВОД РЕЗУЛЬТАТОВ РАСЧЕТА НА ПЕЧАТЬ

443. Пус показатель вида выводимых результатов: Пус=1 - в усилиях, Пус=2- в напряжениях; Пнапр показывает направление вывода результатов:

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

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