Киевский Национальный Университет Строительства и Архитектуры.
Кафедра систем автоматизированного проектирования и управления.
КУРСОВАЯ РАБОТА.
По предмету: «Методология разработки программных продуктов и больших систем».


На тему: «Проектирование напряжённо-деформированного состояния тонкостенных (замкнутых и разомкнутых) оболочечных железобетонных конструкций переменной жёсткости».
                                                        Выполнили:

                                                           студенты группы КСП-42

                                                           Демьяненко Е.И.

                                                           Шепель В.В.

                                                        Проверил:

                                                           Яловец А.Л.
1999г.




1. АНАЛИЗ ПРОБЛЕМЫ. ПОСТАНОВКА ЦЕЛЕЙ.
     В данном курсовом проекте имеет место следующая актуальность темы. Замкнутые и разомкнутые в окружном направлении конические оболочки переменной жёсткости широко используются как конструктивные оболочечные элементы в разных отраслях машиностроения, авиастроения, судостроения, а также строительной индустрии.

1.    Улучшение технико-экономических характеристик и качества проектирования конических оболочек.

1.1           Уменьшение массы конических оболочек.

1.2           Достижение высокой жёсткости и прочности.

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

1.4           Учёт реальных факторов при изготовлении конических оболочек.

2.    Улучшение эксплуатационных характеристик конических оболочек.

2.1           Улучшение поведения конструкции при сложных условиях работы и требования предъявляемые к ним.

2.2           Повышение точности определения факторов при напряжённо-деформированном состоянии конструкции.

2.3           Исследование поведения замкнутых конических оболочек.

2.4           Исследование поведения разомкнутых конических оболочек.

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

3.1           Исследование решения двумерных краевых задач при различных граничных условиях.

3.2           Исследование различных вариационно-разностных и проекционных методов.

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

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

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

     Таким образом, проектирование и моделирование железобетонных тонкостенных замкнутых или разомкнутых оболочечных конструкций на основе сплайн функций является актуальным.
2. ВНЕШНИЕ И ВНУТРЕННИЕ ЗАДАЧИ.
1.    Построение точного решения (погрешность не более 5%) об изгибе ортотропных конических оболочках асимметричного строения под действием нормальной поверхностной нагрузке и температурного поля. Разработка методов численного решения двухмерных краевых задач для замкнутых и разомкнутых конических оболочек поворота шаровидной структуры с изотропными и ортотропными слоями, изменяемыми в двух координатных направлениях жёсткости, которые находятся под действием асимметричных силовых и температурных нагрузок, на основе сплайн аппроксимации.

1.1           Разбиение заданного отрезка исследования на N равных частей с помощью сетки точек.

1.2           Выполнение выборки N+1 точек коллокации для расчёта В-сплайнов.

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

1.4           Подстановка решения данных уравнений в заданные граничные условия на криволинейных контурах.

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

2.1           Выполнение анализа как влияет угол конусности на напряжённо-деформированное состояние замкнутой или разомкнутой конической оболочки вращения переменной толщины.

2.2           Выполнение анализа влияния ортотропии на напряжённо-деформированное состояние замкнутой или разомкнутой конической оболочки.

2.3           Расчёт деформации конических оболочек при различных способах закрепления контуров.




3. ДЕКОМПОЗИЦИЯ СИСТЕМЫ. ПОСТРОЕНИЕ ИЕРАРХИЧЕСКОЙ СТРУКТУРЫ СИСТЕМЫ. РАЗЛОЖЕНИЕ ПРОГРАММНОГО КОМПЛЕКСА НА УРОВНИ.
     Способ управления сложными системами был известен ещё в древности. При проектировании сложной программной системы необходимо разделять её на всё меньшие и меньшие подсистемы, каждую из которых можно совершенствовать независимо. В самом деле, декомпозиция вызвана сложностью программирования системы, поскольку именно эта сложность вынуждает делить пространство состояний системы. Декомпозицию системы можно разделить на два основных вида:

Ø алгоритмическая декомпозиция;

Ø объектно-ориентированная декомпозиция.

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

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

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

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

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

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

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

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

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

На рисунке 1 показана декомпозиция объекта проектирования.




                                           6                    8                   11                14




                                           7                    9                   12                15




                                                                10                  13
Рис.1. Декомпозиция объекта проектирования.
     0.- Тонкостенная железобетонная оболочечная конструкция.

     1.- Подсистема исследования состояния спокойствия.

     2.- Подсистема проверки на наличие дефекта.

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

     4.- Подсистема проверки оболочки на прочность в упругом состоянии.

     5.- Подсистема расчёта околоарматурных напряжений.

     6.- Процедура исследования бездефектного околоарматурного состояния.

     7.- Процедура исследования конструкции с околоарматурным состоянием.

     8.- Процедура исследования напряженно-деформированного состояния от температурной нагрузки.

     9.- Процедура исследования напряженно-деформированного состояния от поверхностной нагрузки.

     10.- Процедура исследования напряженно-деформированного состояния от комбинированной нагрузки.

     11.- Процедура расчёта деформированных сред и выражения углов поворота нормали.

     12.- Процедура соотношения упругости при поверхностных и температурных нагрузках.

     13.- Процедура получения уравнения равновесия.

     14.- Процедура расчёта физико-механических характеристик.

     15.- Процедура расчёта В-сплайнов для получения точного решения.

На рисунке 2 показано разложение на уровни программного комплекса.





Рис.2. Разложение программного комплекса на уровни.




4. ПОСТРОЕНИЕ ЛОГИЧЕСКИХ СХЕМ.
     При построении логических схем пользуются таким термином как логическая ячейка, которую можно представить в следующем виде:




                                                    S_n.

	Mn. Tn. Kn.

                             A_n.




                                                                         R_n.

                             C_n.
     Где её компонентами являются следующие значения:

S_n.- Поставленная задача;

A_n.- Исходные данные;

C_n.- Ограничения;

M_n.- Модуль;

T_n.- Решаемая процедура;

K_n.- Оценка качества;

R_n.- Проектное решение.

     Рассмотрев, что собой представляет логическая ячейка перейдём к рассмотрению задач в данном курсовом проекте.


Логическая ячейка 1.




                                                    S_1.

	M1. T1. K1.

                                 A_1.

                                                                      R_1.

                                 C_1.





S_1.   -Вычисление всякого В-сплайна N–ой степени.

A_1.  -B_nⁱ(x) -наиболее употребительный базис В-сплайн N-ой степени.

       - b_i -некоторые постоянные коэффициенты.

C_1.  - i = -n, … N-1.

T_1.  - S(x) -всякий сплайн N-ой степени.

                            

R_1. - S(x) –всякий сплайн N-ой степени.
Логическая ячейка 2.




                                                      S_2.

	M2. T2. K2.

                                  A_2.

                                                                        R_2.

                                  C_2.





S_2.  - Вычисление В-сплайна нулевой степени.

A_2. - Х - координата функции по оси абсцисс.

C_2. - ограничения имеют вид:

                     

T_2. - В-сплайн нулевой степени.

               

R_2. - В-сплайн нулевой степени.
Логическая ячейка 3.




                                                      S_3.

	M3. T3. K3.

                                  A_3.

                                                                        R_3.

                                  C_3.






S_3. - Получение нового разбиения  и для него рекуррентное соотношение.

A_3. - Входные данные x, x_i, x_i+n, x_i+1, x_i+n+1, Bⁱ_n-1(x).

C_3. - Новое разбиение ,

        (x_-n<x_-n+1<…<x_-1<x₀<x₁<…<x_N<x_N+1<…<x_N+n),

                           n=1,2,…;  ;

T_3. - Рекуррентное соотношение.

           

R_3. - Bⁱ_n(x) –рекуррентное соотношение на выходе ячейки.
Логическая ячейка 4.




                                                      S_4.

	M4. T4. K4.

                                  A_4.

                                                                        R_4.

                                  C_4.





S_4.  - Хранение предыдущего значения сплайна Bⁱ_n-1(x).

A_4. - Предыдущее значение сплайна Bⁱ_n-1(x).

M_4. - Модуль хранения значения со входа ячейки.

T_4.  - Передача хранимого значения на выход.

R_4. - Предыдущее значение сплайна Bⁱ_n-1(x).
Логическая ячейка 5.




                                                      S_5.

	M5. T5. K5.

                                  A_5.

                                                                        R_5.

                                  C_5.





S_5. - Вычисление В-сплайна третей степени Bⁱ₃(x).

A_5. - t -размер шага движения по оси абсцисс.

C_5. - ограничения:

       ;            x < x_i-2;

                                      x_i-2 < x < x_i-1;

                                      x_i-1 < x < x_i;

                                      x_i < x < x_i+1;

                                      x_i+1 < x < x_i+2;

                                      x_i+2 < x.

T_5. - Выражение В-сплайна третей степени.

               

R_5. – Выражение В-сплайна третей степени Bⁱ₃(x).
Логическая ячейка 6.




                                                      S_6.

	M6. T6. K6.

                                  A_6.

                                                                        R_6.

                                  C_6.





S_6. - Вычисление В-сплайна пятой степени Bⁱ₅(x).

A_6. - t -размер шага движения по оси абсцисс.

C_6. - ограничения:

       ;           x < x_i-3;

                                      x_i-3 < x < x_i-2;

                                      x_i-2 < x < x_i-1;

                                      x_i-1 < x < x_i;

                                      x_i < x < x_i+1;

                                      x_i+1 < x < x_i+2;

                                      x_i+2 < x < x_i+3;

                                      x_i+3 < x.

T_6. - Выражение В-сплайна пятой степени.

  

R_6. - Выражение В-сплайна пятой степени Bⁱ₅(x).
Логическая ячейка 7.




                                                      S_7.

	M7. T7. K7.

                                  A_7.

                                                                        R_7.

                                  C_7.





S_7. - Определение количества шагов на заданном интервале.

A_7. - h - шаг разбиения по оси абсцисс.

        Рассматриваемый интервал: (x_i, x_i+1).

T_7. - Количество шагов на заданном интервале t.

                                       ;

R_7. - Количество шагов на заданном интервале h.

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




6. СХЕМА АНАЛИЗА ПРОЕКТА.
     Программный продукт составлен на модульном принципе. Он состоит из отдельных модулей (подпрограмм) реализующих определённые части вычислительного процесса. Модули по своему назначению делятся на модули специального и общего назначения.

Модули специального назначения (модули 1 и 6) зависят от вида граничных условий, нагрузок, структуры пакета слоёв, геометрических характеристик конической оболочки, способа выбора точек коллокации и т.д.

     Структурная схема программы с указанием связей между её модулями показана на рисунке 4.

     Выполнение программы начинается с головной (модуль 1) где задаётся вся исходная информация с клавиатуры пользователем или из базы данных. После задания исходной информации подключаем к работе специальные модули 2 и 3, вычисляющие, соответственно, сплайны 3-й и 5-й степени на заданных отрезках. После этого подключаем к работе специальный модуль 4, который, вычисляет сплайн нулевой степени сплайны 3-й и 5-й степени находит конечный В-сплайн. Специальный модуль 5 производит всю необходимую работу с матрицами. Данные, полученные после работы модулей 4 и 5, используются в программе вычисления правых частей системы обыкновенных дифференциальных уравнений (модуль 6). После задания и получения всей исходной информации управление передаётся модулю 7, в котором реализована вычислительная схема метода дискретной ортогонализации.

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








Рис.4. Структурная схема программы.




К модулям общего назначения относятся:

1)          Модуль 6. В нём производится вычисление вектора правых частей системы обыкновенных дифференциальных уравнений по заданным геометрическим и механическим параметрам оболочки.

2)          Модуль 7. Этот модуль является основной частью программы. В нём реализуется алгоритм решения линейной краевой задачи, производится управление работой остальных модулей.

3)          Модуль 8. Производит вывод на печать значений разрешающих функций в точках выдачи результатов.

4)          Модуль 9. В этом модуле осуществляется решение систем линейных алгебраических уравнений с несколькими правыми частями методом Гаусса с выбором главного элемента.

5)          Модуль 10. В теле данного модуля по заданной информации (в матричном виде) головной программы производится формирование расширенных матриц левых и правых граничных условий.

При этом возможно получение следующих результатов:

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

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

3)       На основе разрабатываемого подхода к решению двумерных краевых задач построен алгоритм численного решения задач о деформации ортотропных замкнутых и разомкнутых конических оболочек вращения переменной жёсткости.

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