Министерство образования и науки Республики Казахстан
Карагандинский Государственный Технический Университет
Кафедра ____САПР______
ПОЯСНИТЕЛЬНАЯ ЗАПИСКА
к курсовой работе
По дисциплине: "Моделирование систем и комплексов"
Тема: "Автоматизированный анализ проектирования на микроуровне"
Руководитель
(подпись) (дата)
Студент
(подпись) (дата)
2009
Содержание
Введение
1. Задание
2. Разработка концептуальной модели и расчетной схемы объекта анализа
3. Выбор метода автоматизированного анализа объекта проектирования
3.1 Метод конечных разностей
3.2 Метод конечных элементов
4. Выбор и краткое описание программных и технических средств автоматизированного анализа
4.1 COSMOS
4.2 ANSYS
5. Построение дискретной модели объекта анализа и программной модели
6. Планирование машинных экспериментов и их реализация
7. Анализ и интерпретация результатов моделирования
Заключение
Список используемой литературы
Приложения
Введение
В данном курсовом проекте требуется спроектировать напряженно-деформированное состояние исследуемого объекта при заданных граничных условиях и установить зависимости для него от геометрических параметров на основе машинного эксперимента.
В ходе выполнения данного курсового проекта нам необходимо будет разработать концептуальную модель объекта, построить расчетную схему объекта моделирования, разработать модель для автоматизированного анализа, создать командный файл, поучить регрессионные зависимости исследуемых данных. Кроме того, по указанию преподавателя необходимо провести качественную и количественную оценки исследуемых характеристик моделируемого объекта для сравнения с результатами имитационного эксперимента.
После ряда прогонов полученной модели на ПЭВМ и получения результатов машинного эксперимента с требуемой точностью необходимо будет провести их интерпретацию и анализ в терминах объекта моделирования, а затем оформить пояснительную записку и графическую часть курсового проекта.
Концептуальная модель необходимо будет проверить на непротиворечивость и полноту описания, для чего надо убедиться, что задание содержит всю необходимую качественную и количественную информацию об объекте. Далее модель проверяют на отсутствие неформализованных и нереализуемых описаний, ситуаций, соотношений. Затем необходимо провести (согласно заданию) качественную и количественную оценки характеристик процесса функционирования объекта.
Итогом курсового проекта является техническая документация в виде разработанных схем и программного обеспечения моделирования объекта, результаты машинного эксперимента с моделью объекта, включая выводы и рекомендации по их использованию при исследовании и разработке реального объекта, пояснительная записка, содержащая документацию по выполнению автоматизированного анализа объекта моделирования.
1. Задание
Требуется спроектировать напряженно-деформированное состояние исследуемого объекта при заданных граничных условиях и установить зависимости для него от геометрических параметров на основе машинного эксперимента.
Материал - сталь;
Модуль упругости E=200000 МПа;
Коэффициент Пуассона v=0,3.
2. Разработка концептуальной модели и расчетной схемы объекта анализа
Целью данного этапа является построение концептуальной (содержательной) модели процесса функционирования объекта с проведением его формализации.
Внешние воздействия на объект проектирования выражены в виде силы нагружения 100 МПа. Силы распределены равномерно.
Таким образом, нам необходимо построить фигуру, изображенную на рисунке 2.1
Рисунок 2.1 - Объект анализа
3. Выбор метода автоматизированного анализа объекта проектирования
Моделирование многочисленных физических, биологических и химических явлений часто приводит к решению линейных и нелинейных уравнений или систем уравнений в частных производных. Существуют традиционные математические средства, позволяющие получить решение в определенных случаях, но для решения конкретных проблем, возникающих в науке и технике, невозможно обойтись без использования численных методов. С ростом производительности ЭВМ численное моделирование приобретает особое значение, т.к позволяет дополнить или даже заменить прямой эксперимент. Последний часто дорог, его постановка бывает трудоемкой или вообще невозможной. В настоящее время существует ряд различных методов автоматизированного анализа. Среди них наиболее популярными являются методы конечных элементов и конечных разностей.
3.1 Метод конечных разностей
Метод конечных разностей относится к методам сеток и содержит три этапа:
1 этап - Дискретизация - на данном этапе область решения разбивается на сетку, как правило, регулярную, в дальнейшем в узлах сетки будут определены узловые значения искомой функции.
2 этап - Аппроксимация - здесь осуществляется переход от дифференциального оператора к разностному аналогу.
3 этап - Алгебраизация - подстановка разностных аналогов дифференциальных операторов в исходное уравнение и получение системы алгебраических уравнений.
3.2 Метод конечных элементов
Этот метод в настоящее время достиг такого уровня, что многие часто сомневаются - может ли появиться лучший метод. Диапазон применимости МКЭ, их эффективность и сравнительная легкость реализации, делают их серьезными соперниками для любого метода. Достоинствами МКЭ являются гибкость и разнообразие сеток, стандартные приемы построения дискретных задач для произвольных областей, простота учета естественных краевых условий. Кроме того, математический анализ МКЭ является более простым, его методы применимы к более широкому классу задач, а оценки погрешностей приближенных решений получаются при менее жестких ограничениях, чем в методе конечных разностей. Слабой стороной является то, что он по идее представляет собой схему дискретизации всего тела, а это неизбежно ведет к большому количеству КЭ, особенно в трехмерных задачах с удаленными границами, в пределах каждой из которых не все неизвестные переменные изменяются непрерывно.
МКЭ в общем случае состоит из трех этапов:
Дискретизация - физическая область решения задачи разбивается на некоторые подобласти или конечные элементы.
Аппроксимация - искомая функция аппроксимируется функцией специального вида на каждом конечном элементе. Коэффициенты аппроксимации становятся основными неизвестными в задачах.
Алгебраизация - подстановка аппроксимирующих коэффициентов в определяющее уравнение позволяет получить систему алгебраических уравнений относительно узловых значений искомой функции.
Сравнивая эти два метода можно сказать, что каждый из методов имеет свои преимущества. Но у метода конечных разностей больше недостатков, например, сложность решения задач в случае неоднородной области решения, наличие нелинейностей.
Отсутствие решения при недостаточных условиях дают право использовать для реализации поставленной задачи МКЭ. МКЭ, в целях упрощения расчета задачи, может быть реализован по средствам программы ANSYS.
4. Выбор и краткое описание программных и технических средств автоматизированного анализа
В настоящее время существует широкий ряд систем инженерных расчетов. В данном курсовом проекте для сравнения были выбраны две программы - ANSYS и COSMOS.
4.1
COSMOS
COSMOS сочетает надежность, точность и быстроту расчетов с простотой использования, а наличие широких расчетных возможностей - с доступностью (цена системы невысока для продукта такого класса).
COSMOS сориентирован на решение задачи сквозной автоматизации проектирования, позволяя интегрировать расчетный анализ в единую программную систему.
COSMOS имеет ряд преимуществ:
Интерфейс пользователя - позволяет быстро осваивать технологию расчетов.
Система команд естественно отражает стадии расчетного анализа.
Возможность углубленного анализа по мере освоения системы.
Простейшие расчеты доступны с первого же сеанса работы.
Оперативные средства графического отображения исходных данных и полученных результатов.
Геометрическое моделирование возможно как внутренними средствами системы, так и на основе данных, импортируемых из популярных пакетов моделирования.
Модульное строение способствует формированию системы в зависимости от потребностей и финансовых возможностей пользователя.
Средства автоматизации подготовки расчетов позволяют проводить оценочные расчеты проектировщикам, не имеющим глубокой подготовки: автоматическое построение пространственных сеток конечных элементов с измельчением ее в местах концентрации напряжений, наличие модуля оптимизации конструкций.
FFE-методы скоростного расчета сокращают время расчета в десятки раз, при этом использование системных ресурсов (оперативная память, дисковое пространство) уменьшается на порядок. Повышение производительности расчетов при этом качественно изменяет подход к использованию COSMOS: в результате отбрасывается необходимость упрощения сложных пространственных моделей, а это уже - значительный скачок в производительности и оперативности процесса конструкторской доводки.
Функциональные возможности:
Определение параметров напряженно-деформированного состояния
Тепловые расчеты
Динамика
Физическая и геометрическая нелинейность
Ламинарные и турбулентные потоки
Электромагнитные расчеты
Оптимизация конструкций
Библиотека материалов
Подготовка данных для других программ МКЭ
4.2
ANSYS
Многоцелевой конечно-элементный пакет для проведения анализа в широком круге инженерных задач. Особенностью программы является файловая совместимость всех членов семейства ANSYS для всех используемых платформ. Многоцелевая направленность программы позволяет использовать одну и ту же модель для решения таких связанных задач, как прочность при тепловом нагружении, влияние магнитных полей на прочность конструкции, тепломассоперенос в электромагнитном поле.
Внутренний язык программирования APDL позволяет программировать любые процедуры, параметризовать построение модели, расчет и вывод результатов. ANSYS позволяет подключать пользовательские модели, написанные на FORTRAN, C++.
Программа ANSYS обладает многими возможностями конечно-элементного анализа - от простого линейного статического до сложного нелинейного динамического (нестационарного). Процедура типового расчета может быть разделена на три основных этапа:
Построение модели.
Приложение нагрузок (включая и граничные условия) и получение решения.
Просмотр и анализ результатов.
Построение модели (Preprocessor)
Этот этап включает определение типов конечных элементов, их констант, свойств материала и геометрии модели.
Задание типов элементов (ElementType)
Библиотека конечных элементов программы ANSYS содержит более 100 типов, каждый из которых определяет, среди прочего, применимость элемента к той области расчетов (прочностной, тепловой, магнитный и электрический анализ, движение жидкости или связанные задачи), характерную форму элемента (линейную, плоскую, в виде бруска и т.д.), а также двумерность и трехмерность элемента как геометрического тела.
Задание констант элементов (RealConstants)
Константы элемента - это свойства специфические для данного типа элемента, такие как параметры поперечного сечения балочного элемента.
Задание свойств материалов (MaterialProps)
Свойства материала требуются для большинства типов элементов. В зависимости от области приложения свойства могут быть линейными, нелинейными иили анизотропными.
Создание геометрической модели (-Modeling-)
Основной целью на этапе разработки геометрической модели является создание адекватной конечно-элементной модели, состоящей из узлов и элементов. При создании конечно-элементной модели используются три метода: твердотельное моделирование, каркасное моделирование и прямая генерация сетки.
Задание нагрузок и получение решения (Solution)
На этом этапе выбирается тип анализа и установление его опций, прикладываются нагрузки, определяются опции для выбора шага по нагрузке и инициируется решение.
Приложение нагрузок (-Loads-)
Под нагрузками понимаются как внешние и внутренние усилия, так и граничные условия в виде ограничений на перемещения. В программе ANSYS нагрузки разделены на 6 категорий:
ограничения степени свободы;
силы;
поверхностные нагрузки;
объемные силы;
инерционные нагрузки;
нагрузки для связных задач.
Большинство этих нагрузок может быть приложено или к твердотельной модели, или к конечно-элементной модели.
Запуск на счет
По команде SOLVE программа обращается за информацией о модели и нагрузках к базе данных и выполняет вычисления. Результаты записываются в специальный файл и базу данных. При этом в базе данных может храниться только один набор результатов, тогда как в файл могут быть записаны результаты для всех шагов решения.
Просмотр результатов
Постпроцессор предназначен для визуализации результатов задачи. При этом как правило используют:
Наложение деформированной сетки на недеформированную.
Построение эпюр (перемещений, напряжений и деформаций) для заданной области решений.
В программе ANSYS доступны следующие виды прочностного анализа:
Статический анализ - вычисление перемещений, напряжений и т.д. в условиях статического нагружения;
Модальный анализ - определение собственных частот и форм колебаний;
Гармонический анализ - определение отклика конструкции на гармонические составляющие возмущающей нагрузки;
Динамический анализ - определение отклика конструкции на действие произвольной нагрузки как функции времени;
Спектральный анализ - расширение модального анализа для вычисления напряжений и деформаций при действии спектра частот или случайной вибрации;
Анализ устойчивости - расчет критических нагрузок и определение форм потери устойчивости.
Таким образом, сравнивая две программы, учитывая преимущества и недостатки каждой для данного курсового проекта был выбран программный пакет ANSYS.
5. Построение дискретной модели объекта анализа и программной модели
Прежде всего, необходимо проанализировать свою деталь. Если, выданная в качестве задания, деталь симметрична, то следует строить только одну часть детали. В данной работе мы построим только четверть детали (Рисунок 5.1), так как деталь можно разбить на 4 одинаковые части.
Рисунок 5.1 - Вид четверти детали
В данной задаче параметрически заданными (то есть их можно будет менять в программе) будут следующие параметры: r, R, h (рисунок 5.2).
r- внутренний радиус сквозного отверстия;
R- радиус внутренней выемки;
h- высота внутренней выемки;
|
|
|
Рисунок 5.2 - Базовые размеры
Нагрузку приложим на область S (см. рисунок 5.3).
|
Рисунок 5.3 - Вид области с нагрузкой
После решения выше изложенных вопросов можно переходить к построению.
Запуск программы ANSYS
1. Первичный запуск (запускаем первый раз за сеанс работы)
Пуск>Программы>ANSYSED55>Interactive.
В поле WorkingDirectory ввести свою рабочую директорию (имя директории и все пути к ней должны быть на английском языке).
Например: С: BASKETPLATE
Далее нажать RUN
2. Последующие запуски.
Пуск>Программы>ANSYSED55>RunInteractiveNow
Примечание: в скобках будут указаны команды в блокноте, соответствующие командам в программе ANSYS.
1. Определение имени файла
Верхнее меню File>ChangeJobname …
Ввести имя файла (на английском языке)
OK
(/FILNAM, <имя файла>)
2. Определение заголовка задачи
Верхнее меню File>Change Title …
Ввести заголовок задачи (на англ)
OK
(/TITLE, <заголовок>)
3. Вход в препроцессор
Выбрать пункт Preprocessor главного меню
(/PREP7)
4. Задание типа элемента
Main Menu >Preprocessor> Element Type> Add/Edit/Delete
Выбрать Add…
Solid10node 92
Close
(ET,1,92)
5. Задание свойств материала
Main Menu>Preprocessor>Material Props>Polynomial
EX - модуль упругости - ввести 200000
(MP,EX,1, 200000)
NUXY - коэффициент Пуассона - ввести 0.3
(MP,NUXY,1,0.3)
6. Задание базовых размеров
(r=6)
(R=13)
(c=6)
7. Построение контрольных точек (рисунок 5.4):
Main menu>Preprocessor>Modeling>Create>Keypoints>In Active Cs
(k,1,0,0,0! задание точек для основания
k,2,60,0,0
k,3,60,-50,0
k,4,15,-50,0
k,5,15,-40,0
k,6,0,-40,0)
Рисунок 5.4 - Задание контрольных точек
8. Построение линий через контрольные линии (рисунок 5.5):
Main menu>Preprocessor>Modeling>Create>Lines>Lines>In Active Coord
(L,1,2! задание линий для основания
L,2,3
L,3,4
L,4,5
L,5,6
L,6,1)
Рисунок 5.5 - Построение линий
9. Объединение линий в область (рисунок 5.6):
Main Menu>Preprocessor>Modeling>Create>Areas>Arbitrary>By Lines
выбрать мышкой все линии, затем нажать ОК.
(al,all) ! объединение линий в область
Рисунок 5.6 - Объединение линий в область
10. Полученную область необходимо вытянуть (рисунок 5.7):
Modeling-Operate - Extrude - Areas - By XYZ offset
(voffst,1,25) ! вытягивание области на высоту
Рисунок 5.7 - Вытягивание области на высоту
11. Замена локальной системы координат:
WorkPlane-Local Coordinate System - Create Local CS - At WP Origin…
(local,100,0,0,0,-25,0,0,0,0,,
wpcsys,1,100)
12. Построение контрольных точек и соединение их линиями в новой системе координат, затем объединяем их в область и вытягиваем на высоту. Результат этого этапа приведен на рисунке 5.8.
(k,13,0,0,0
k,14,25,0,0
k,15,0,-25,0
L,13,14
L,14,15
L,15,13
al, 19, 20,21
voffst,9,35)
Рисунок 5.8 - Построение детали
13. Затем опять меняем систему координат и строим контрольные точки, соединяем их в линии.
(local,100,0,0,0,-60,0,0,0,0,,
wpcsys,1,100
k, 19,0,-13,0
k, 20,13,0,0
L,16, 19
L,16, 20)
14. Строим дугу, затем две линии и дугу объединяем в область и вытягиваем ее на высоту.
(larc, 19, 20,16,13! построение дуги
al,28,29,30
voffst,14,6)
15 Вычитание из одного объёма другого (Рисунок 5.9):
Main menu>Preprocessor>Modeling>Operate>Subtract>Volumes
(vsbv,2,3! вычитание объёмов)
Рисунок 5.9 - Вычитание объема
16. Затем опять меняем систему координат и строим контрольные точки, соединяем их в линии, строим дугу, все это объединяем в область и вытягиваем на высоту и таким образом получаем цилиндр.
(local,100,0,0,0,-54,0,0,0,0,,
wpcsys,1,100
k,24,0,-6,0! построение контрольных точек
k,25,6,0,0
L,21,24! построение линий
L,21,25
larc,24,25,21,6! построение дуги
al,22,24,25! объединение в область
voffst,10,54) ! вытягивание на высоту
17. Объединение объёмов:
Main Menu>Preprocessor>Modeling>Operate>Add>Volumes
Выбрать мышкой объёмы и нажать ОК.
(vadd,1,4) ! объединение объёмов
18 Затем вычитаем из одного объема полученный ранее цилиндр (рисунок 5.10).
(vsbv,3,2) ! вычитание объемов
Рисунок 5.10 - Вычитание объемов
19. Задание размера конечных элементов:
Main Menu>Preprocessor>Meshing>Size Cntrls>ManualSize>Global>Size
В поле Size ввести размер конечного элемента = 5.
(smrtsize,all,5) ! задание размера конечного элемента
20. Разбиение объема на конечные элементы (рисунок 5.11):
Main Menu>Preprocessor>Meshing>Mesh>Volumes>FreeMain
(vmesh,all)
Рисунок 5.11 - Разбиение на конечные элементы
21. Выход из препроцессора:
(FINISH)
22. Вход в процессор:
(/SOLU)
23. Закрепление областей по всем степеням свободы:
Main Menu>Solution>Define Loads>Apply>Structural>Displacement>On AreasMain
Выделить мышкой нужные области и нажать ОК.
(da,26,all)
24. Закрепление симметрии областей:
Menu>Solution>Define Loads>Apply>Structural>Displacement>Symmetry B. C. >On Areas
Выделить мышкой нужные области и нажать ОК.
(da,8,uy
da,21,uy
da,3,ux
da, 20,ux)
25. Приложение нагрузки к областям (рисунок 5.12):
Main Menu>Solution>Define Loads>Apply>Structural>Pressure>On Areas
Выбрать мышкой нужные области, в поле Valueнабрать 100, и нажать ОК.
(sfa,9,,pres,100)
Рисунок 5.12 - Приложение нагрузки
26. Выполнение расчета (SOLVE)
27. Выход из процессора (FINISH)
28. Вход в постпроцессор общего пользования
(/POST1)
29. Просмотр результатов (рисунок 5.13).
Main Menu>General Postproc>Plot Results
Выбрать режим просмотра.
(plnsol,u,sum,0,1)
Рисунок 5.13 - Просмотр результатов
6. Планирование машинных экспериментов и их реализация
Данный машинный эксперимент будет включать в себя изменение параметров модели с отклонением ±5 мм.
Диапазон параметров изменяется следующим образом:
r от 4 до 9;
R от 11 до 15;
H от 1 до 10;
Изменяя параметры модели, мы посредством программы ANSYS, снимали значения σx
, σy
, σz
. Так как изменяемых параметров было 3 и нужно было провести опыты во всех комбинациях, мы получили 27 опытов, которые представлены в таблице 6.1.
Таблица 6.1 - Результаты практических опытов
| № опыта | r | R | H | σx
|
σy
|
σz
|
| 1 | 6 | 13 | 6 | 61.84 | 62.746 | 107.573 |
| 2 | 6 | 13 | 1 | 49.494 | 47.841 | 96.654 |
| 3 | 6 | 13 | 10 | 54.293 | 64.308 | 97.772 |
| 4 | 6 | 11 | 6 | 73.727 | 68.827 | 108.093 |
| 5 | 6 | 15 | 6 | 58.683 | 54.223 | 106.188 |
| 6 | 6 | 11 | 1 | 52.754 | 53.011 | 107.432 |
| 7 | 6 | 11 | 10 | 44.734 | 56.755 | 105.237 |
| 8 | 6 | 15 | 1 | 46.811 | 48.589 | 108.531 |
| 9 | 6 | 15 | 10 | 47.373 | 52.555 | 97.67 |
| 10 | 4 | 13 | 6 | 46.682 | 48.097 | 81.004 |
| 11 | 4 | 11 | 6 | 48.93 | 41.278 | 82.495 |
| 12 | 4 | 15 | 6 | 53.186 | 47.617 | 80.986 |
| 13 | 4 | 11 | 1 | 38.02 | 44.192 | 95.672 |
| 14 | 4 | 11 | 10 | 49.812 | 41.77 | 91.161 |
| 15 | 4 | 15 | 1 | 37.227 | 37.348 | 95.51 |
| 16 | 4 | 15 | 10 | 41.7 | 44.667 | 91.307 |
| 17 | 9 | 11 | 6 | 101.439 | 114.478 | 131.362 |
| 18 | 9 | 11 | 1 | 168.002 | 159.947 | 125.308 |
| 19 | 9 | 11 | 10 | 92.06 | 90.598 | 126.622 |
| 20 | 9 | 13 | 6 | 111.459 | 110.351 | 120.552 |
| 21 | 9 | 13 | 1 | 148.791 | 163.141 | 118.298 |
| 22 | 9 | 13 | 10 | 102.164 | 115.427 | 120.214 |
| 23 | 9 | 15 | 6 | 104.629 | 146.047 | 127.884 |
| 24 | 9 | 15 | 1 | 158.741 | 148.498 | 129.646 |
| 25 | 9 | 15 | 10 | 107.96 | 115.379 | 129.244 |
| 26 | 4 | 13 | 1 | 41.357 | 38.541 | 93.753 |
| 27 | 4 | 13 | 10 | 38.244 | 35.417 | 94.541 |
Анализ полученных зависимостей проводился с помощью программы ANETR.
Созданный и реализованный нетрадиционный метод программы ANETR базируется на разработках профессора М.М. Протодьяконова. Он позволяет работать в данными управляемого и неуправляемого экспериментов.
Метод универсален как в области применения, так и в вариантах построения модели: модели могут иметь вид суммы, произведения частных зависимостей, их комбинации, с последовательной нейтрализацией влияния приоритетных аргументов или без таковой. Важным преимуществом метода является возможность с должным уровнем надежности оценить влияние тех аргументов, которые в традиционных методах будут исключены из анализа как не влияющие на целевую функцию существенно (например, по критерию Фишера). Это обеспечивается процедурой последовательной нейтрализации влияния на функцию сильно действующих аргументов, при которой и начинает проявляться влияние более "слабых" аргументов. Метод детерминирован, т.е. позволяет ограничивать выбор модели известными данными по теории процесса (тип уравнений, пределы и т.д.).
На выходе программа, кроме аналитического характера частных связей, дает их графики, оценку надежности. Ранжирует факторы по силе воздействия на результат, позволяет выделить существенно влияющие факторы. Графики частных связей позволяют проводить оптимизацию процесса без применения специального математического аппарата.
Разработка моделей сводится к последовательному выполнению следующих этапов:
Планирование эксперимента;
Проведение эксперимента;
Составление заказа на модели;
Анализ модели.
Для начала мы формируем файл для работы с программой ANETR:
Для σx
:
B=6; E=2; G=1; L=11; P=1;
4=F (1; 2;3);
1 6 13 6 61.84 62.746 107.573
2 6 13 1 49.494 47.841 96.654
3 6 13 10 54.293 64.308 97.772
4 6 11 6 73.727 68.827 108.093
5 6 15 6 58.683 54.223 106.188
6 6 11 1 52.754 53.011 107.432
7 6 11 10 44.734 56.755 105.237
8 6 15 1 46.811 48.589 108.531
9 6 15 10 47.373 52.555 97.67
10 4 13 6 46.682 48.097 81.004
11 4 11 6 48.93 41.278 82.495
12 4 15 6 53.186 47.617 80.986
13 4 11 1 38.02 44.192 95.672
14 4 11 10 49.812 41.77 91.161
15 4 15 1 37.227 37.348 95.51
16 4 15 10 41.7 44.667 91.307
17 9 11 6 101.439 114.478 131.362
18 9 11 1 168.002 159.947 125.308
19 9 11 10 92.06 90.598 126.622
20 9 13 6 111.459 110.351 120.552
21 9 13 1 148.791 163.141 118.298
22 9 13 10 102.164 115.427 120.214
23 9 15 6 104.629 146.047 127.884
24 9 15 1 158.741 148.498 129.646
25 9 15 10 107.96 115.379 129.244
26 4 13 1 41.357 38.541 93.753
27 4 13 10 38.244 35.417 94.541
/*
Для σy
:
Меняем только верхнюю часть файла. Куда вводим данные:
B=6; E=2; G=1; L=11; P=2;
5=F (1; 2;
3);
Для σz
:
Меняем только верхнюю часть файла. Куда вводим данные:
B=6; E=2; G=1; L=11; P=1;
6=F (1; 2;
3);
После обработки мы получаем следующие результаты
Для σx
:
Аппроксимирующие уравнения:
По 1-му фактору:
Y1= A * X**2 + B * X + C
Где A= 3.42944E+00 B= - 2.90415E+01 C= 1.05229E+02
По 2-му фактору:
Y2= A * X**2 + B * X + C
Где A= 2.27839E-01 B= - 6.26915E+00 C= 1.15724E+02
По 3-му фактору:
Y3= A + B * X
Где A= 8.46799E+01 B= - 2.00161E+00
Вид обобщенного уравнения
Y=Y (1) *Y (2) *Y (3) / 5.378386E+03
Подставив в уравнение Y уравнения по всем трем факторам с соответствующими коэффициентами получим:
Данное уравнение является зависимостью σx
от базовых размеров.
СКО%= 32.11
Для σy
:
Аппроксимирующие уравнения:
По 1-му фактору:
Y1= 1/ (A + B * X)
Где A= 3.67374E-02 B= - 3.21271E-03
По 2-му фактору:
Y2= A + B / X
Где A= 8.46504E+01 B= - 1.10833E+02
По 3-му фактору:
Y3 = A + B * X
Где A= 8.45668E+01 B= - 1.51408E+00
Вид обобщенного уравнения
Y= Y (1) *Y (2) *Y (3) / 5.774019E+03
Подставив в уравнение Y уравнения по всем трем факторам с соответствующими коэффициентами получим:
Данное уравнение является зависимостью σy
от базовых размеров.
СКО%
Для σz
:
Аппроксимирующие уравнения:
По 1-му фактору:
Y1= A + B * X
Где A= 6.08987E+01 B= 7.17218E+00
По 2-му фактору:
Y2= A * X**2 + B * X + C
Где A= 1.10477E+00 B= - 2.88999E+01 C= 1.86047E+02
По 3-му фактору:
Y3 = A * X**2 + B * X + C
Где A= 8.42139E-02 B= - 1.13545E+00 C= 2.58845E+00
Вид обобщенного уравнения
Y= Y (1) +Y (2) +Y (3)
Подставив в уравнение Y уравнения по всем трем факторам с соответствующими коэффициентами получим:
Данное уравнение является зависимостью σz
от базовых размеров.
СКО%= 28.69
7. Анализ и интерпретация результатов моделирования
В результате пробного расчета при параметрах r = 5 R=12 H=3, с помощью программы ANSYS было получено:
σx
= 48,517 σy
= 59,969 σz
= 89, 192
При расчете с помощью уравнений полученных через ANETR мы получили:
σx
= 49.06472 σy
= 50.55656 σz
= 95.03471
Сравнивая значения пробного расчета, полученные с помощью ANSYS можно сказать, что расхождения совсем незначительные.
Данные значения СКО показывают, что уравнения для модели были найдены правильно и их возможно использовать в практических задачах в заданном диапазоне.
Заключение
В ходе выполнения данного курсового проекта были решены все поставленные задачи, а именно:
1. Создана расчетная схема анализа на основании сравнительного анализа численных методов, а также программных и технических средств их осуществления;
2. Создан выбор метода автоматизированного анализа объекта проектирования;
3. Спланирован и проведен эксперимент, анализируя результаты которого, приходим к выводу, что данная модель может использоваться с параметрами: r = 5 R = 12 H=3
Список используемой литературы
1. Нургужин М.Р., Степанов П.Б. Прикладная теория систем. Микро - и макромоделирование: учебное пособие. - Караганда: КарПТИ, 1994. -132с.
2. Нургужин М.Р., Даненова Г.Т., Кацага Т.Я. Математическое моделирование: курсовое проектирование. - Караганда: КарГТУ, 2001. -102с.
3. Нургужин М.Р., Даненова Г.Т. Моделирование систем и комплексов: Лабораторный практикум: Учебное пособие. - Караганда: КарГТУ, 2006. -97с.
Приложения
Приложение А
Результаты ANETR для σx
B=6; E=2; G=1; L=11; P=2;4=F (1; 2;3);
1 6.18400E+01 6.00000E+00 1.30000E+01 6.00000E+00
2 4.94940E+01 6.00000E+00 1.30000E+01 1.00000E+00
3 5.42930E+01 6.00000E+00 1.30000E+01 1.00000E+01
4 7.37270E+01 6.00000E+00 1.10000E+01 6.00000E+00
5 5.86830E+01 6.00000E+00 1.50000E+01 6.00000E+00
6 5.27540E+01 6.00000E+00 1.10000E+01 1.00000E+00
7 4.47340E+01 6.00000E+00 1.10000E+01 1.00000E+01
8 4.68110E+01 6.00000E+00 1.50000E+01 1.00000E+00
9 4.73730E+01 6.00000E+00 1.50000E+01 1.00000E+01
10 4.66820E+01 4.00000E+00 1.30000E+01 6.00000E+00
11 4.89300E+01 4.00000E+00 1.10000E+01 6.00000E+00
12 5.31860E+01 4.00000E+00 1.50000E+01 6.00000E+00
13 3.80200E+01 4.00000E+00 1.10000E+01 1.00000E+00
14 4.98120E+01 4.00000E+00 1.10000E+01 1.00000E+01
15 3.72270E+01 4.00000E+00 1.50000E+01 1.00000E+00
16 4.17000E+01 4.00000E+00 1.50000E+01 1.00000E+01
17 1.01439E+02 9.00000E+00 1.10000E+01 6.00000E+00
18 1.68002E+02 9.00000E+00 1.10000E+01 1.00000E+00
19 9.20600E+01 9.00000E+00 1.10000E+01 1.00000E+01
20 1.11459E+02 9.00000E+00 1.30000E+01 6.00000E+00
21 1.48791E+02 9.00000E+00 1.30000E+01 1.00000E+00
22 1.02164E+02 9.00000E+00 1.30000E+01 1.00000E+01
23 1.04629E+02 9.00000E+00 1.50000E+01 6.00000E+00
24 1.58741E+02 9.00000E+00 1.50000E+01 1.00000E+00
25 1.07960E+02 9.00000E+00 1.50000E+01 1.00000E+01
26 4.13570E+01 4.00000E+00 1.30000E+01 1.00000E+00
27 3.82440E+01 4.00000E+00 1.30000E+01 1.00000E+01
Min 3.72270E+01 4.00000E+00 1.10000E+01 1.00000E+00
Max 1.68002E+02 9.00000E+00 1.50000E+01 1.00000E+01
Ћ‘ = 7.334E+01 ‘ЉЋ = 3.899E+01 ‘Љ’ = 27
---------------------------------------------------------------
3_
3.72E+01 6.34E+01 8.95E+01 1.16E+02 1.42E+02 1.68E+02
X Y YP KT 4.00E+00+---------+---------+---------+---------+---------+-
4.00000E+00 4.39064E+01 4.39337E+01 9! 3 2@211
! *
! *
! *
4.50E+00+ *
! *
! *
! *
! *
5.00E+00+ *
! *
! *
! *
! *
5.50E+00+ *
! *
! *
! *
! *
6.00E+00+ *
6.00000E+00 5.44121E+01 5.44394E+01 9! 1211@11 1
! *
! *
! *
6.50E+00+ *
! *
! *
! *
! *
7.00E+00+ *
! *
! *
! *
! *
7.50E+00+ *
! *
! *
! *
! *
8.00E+00+ *
! *
! *
! *
! *
8.50E+00+ *
! *
! *
! *
! *
9.00E+00+ *
9.00000E+00 1.21694E+02 1.21639E+02 9! 1 2111 @ 1 1 1
3$
‚€„ 11 “ђЂ‚Ќ…Ќ€џ ЏЋ 1 ”ЂЉ’Ћђ“ Y = A * X**2 + B * X + C
A= 3.42944E+00 B= - 2.90415E+01 C= 1.05229E+02
CKO= 1.69351E+01 CKO%= 43.4 CKO CP. = 3.92639E-02
3_
3.72E+01 6.34E+01 8.95E+01 1.16E+02 1.42E+02 1.68E+02
X Y YP KT 1.10E+01+---------+---------+---------+---------+---------+-
1.10000E+01 7.43864E+01 7.43320E+01 9! 1 1111 @ 1 1 1
! *
! *
! *
1.14E+01+ *
! *
! *
! *
! *
1.18E+01+ *
! *
! *
! *
! *
1.22E+01+ *
! *
! *
! *
! *
1.26E+01+ *
! *
! *
! *
! *
1.30E+01+ *
1.30000E+01 7.27027E+01 7.27299E+01 9! 1 1 11 1 1 @ 1 1 1
! *
! *
! *
1.34E+01+ *
! *
! *
! *
! *
1.38E+01+ *
! *
! *
! *
! *
1.42E+01+ *
! *
! *
! *
! *
1.46E+01+ *
! *
! *
! *
! *
1.50E+01+ *
1.50000E+01 7.29233E+01 7.29506E+01 9! 1 1 2 1 1 @ 11 1
3$
‚€„ 11 “ђЂ‚Ќ…Ќ€џ ЏЋ 2 ”ЂЉ’Ћђ“ Y = A * X**2 + B * X + C
A= 2.27839E-01 B= - 6.26915E+00 C= 1.15724E+02
CKO= 3.89802E+01 CKO%= 100.0 CKO CP. = 3.92639E-02
3_
3.72E+01 6.34E+01 8.95E+01 1.16E+02 1.42E+02 1.68E+02
X Y YP KT 1.00E+00+---------+---------+---------+---------+---------+-
1.00000E+00 8.23552E+01 8.26783E+01 9! 2 1 111 @ 1 1 1
! *
! *
! *
1.90E+00+ *
! *
! *
! *
! *
2.80E+00+ *
! *
! *
! *
! *
3.70E+00+ *
! *
! *
! *
! *
4.60E+00+ *
! *
! *
! *
! *
5.50E+00+ *
! *
! *
6.00000E+00 7.33972E+01 7.26703E+01 9! 2 1 11 @ 11 1
! *
6.40E+00+ *
! *
! *
! *
! *
7.30E+00+ *
! *
! *
! *
! *
8.20E+00+ *
! *
! *
! *
! *
9.10E+00+ *
! *
! *
! *
! *
1.00E+01+ *
1.00000E+01 6.42600E+01 6.46639E+01 9! 1 1111 1 @ 1 1 1
3$
‚€„ 1 “ђЂ‚Ќ…Ќ€џ ЏЋ 3 ”ЂЉ’Ћђ“ Y = A + B * X
A= 8.46799E+01 B= - 2.00161E+00
CKO= 3.82575E+01 CKO%= 98.1 CKO CP. = 5.24895E-01***********************************************************************************************************************
‘ЉЋ = 1.25204418E+01 R = 9.39903264E-01
***********************************************************************************************************************
Ћ ’ Љ ‹ Ћ Ќ … Ќ € џ
27 - 1.72726E-01 9 - 3.74577E-01 25 1.27272E+00 16 3.16672E+00 10 3.50867E+00 7 - 3.91772E+00
22 - 4.20056E+00 11 4.80567E+00 5 5.02351E+00 3 6.68985E+00 26 - 7.76209E+00 20 - 8.07520E+00
1 8.34282E+00 6 - 9.45141E+00 12 9.88168E+00 14 1.05491E+01 2 - 1.13707E+01 15 - 1.20411E+01
13 - 1.21811E+01 21 1.27948E+01 8 - 1.42384E+01 23 - 1.52679E+01 19 - 1.66475E+01 4 1.90514E+01
17 - 2.07282E+01 24 2.23321E+01 18 2.90101E+01
***********************************************************************************************************************
‚€„ ЋЃЋЃ™…ЌЌЋѓЋ “PЂ‚Ќ…Ќ€џ
Y (4) =Y (1) *Y (2) *Y (3) / 5.378386E+03
***********************************************************************************************************************
ѓђЂ”€Љ €‡Њ…Ќ…Ќ€џ ‘ Љ Ћ
3_
0.00E+00 7.80E+00 1.56E+01 2.34E+01 3.12E+01 3.90E+01
OT ‘ђ…„Ќ…ѓЋ ‡ЌЂ-…Ќ€џ ”“ЌЉ-€€ ‘ Љ Ћ = 3.89876382E+01+---------+---------+---------+---------+---------*
! 20 40 60 80 100%
!
ЏЋ‘‹… *** ”ЂЉ’ЋђЂ ‘ Љ Ћ = 1.25204418E+01+ *
3$
***********************************************************************************************************************
”ЋPЊ“‹џP ЊЋ„…‹€
„Ђ’Ђ: 28/4/2008 Ћ’PЂЃЋ’ЂЌЋ 27 ЋЏ›’Ћ‚ ‡ЂЉЂ‡-€Љ:
B=6; E=2; G=1; L=11; P=2;
4=F (1; 2;
3);
------------------------------------------------------------------------------------------------------------------
! N! N! KЋќ””€-€…Ќ’› “PЂ‚Ќ…Ќ€‰! CKO%!! ‘ЉЋ. ‘ђ%!
! ЂPѓ! “P-џ! - ----------------------------------------------------------! - -------------! R! - -------------!
!!! A! B! C! D! 100!! 100!
! - ---! - ---! - -------------! - -------------! - -------------! - -------------! - -------------! - -------------! - -------------!
!!!!!!!!!!
! 1! 11! 3.4294362E+00! - 2.9041529E+01! 1.0522883E+02!! 43.4!! 0.1!
!!!!!!!!!!
! 2! 11! 2.2783876E-01! - 6.2691519E+00! 1.1572414E+02!! 100.0!! 8.2!
!!!!!!!!!!
! 3! 1! 8.4679913E+01! - 2.0016056E+00!!! 98.1!! 109.6!
!!!!!!!!!!
! ***! ***!!!!! 32.1! 9.3990326E-01!!
!!!!!!!!!!
------------------------------------------------------------------------------------------------------------------
DMAX = 2.9010124E+01 ( 18)
CKO% = 32.11
ЉЋќ””€-€…Ќ’ ”€?…PЂ F = 9.69649; F (001) = 2.52846; F (005) = 1.92428
‚›ЏЋ‹Ќ€‹
***********************************************************************************************************************
Приложение Б
Результаты ANETR для σy
B=6; E=2; G=1; L=11; P=2;
5=F (1; 2;3);
1 6.27460E+01 6.00000E+00 1.30000E+01 6.00000E+00
2 4.78410E+01 6.00000E+00 1.30000E+01 1.00000E+00
3 6.43080E+01 6.00000E+00 1.30000E+01 1.00000E+01
4 6.88270E+01 6.00000E+00 1.10000E+01 6.00000E+00
5 5.42230E+01 6.00000E+00 1.50000E+01 6.00000E+00
6 5.30110E+01 6.00000E+00 1.10000E+01 1.00000E+00
7 5.67550E+01 6.00000E+00 1.10000E+01 1.00000E+01
8 4.85890E+01 6.00000E+00 1.50000E+01 1.00000E+00
9 5.25550E+01 6.00000E+00 1.50000E+01 1.00000E+01
10 4.80970E+01 4.00000E+00 1.30000E+01 6.00000E+00
11 4.12780E+01 4.00000E+00 1.10000E+01 6.00000E+00
12 4.76170E+01 4.00000E+00 1.50000E+01 6.00000E+00
13 4.41920E+01 4.00000E+00 1.10000E+01 1.00000E+00
14 4.17700E+01 4.00000E+00 1.10000E+01 1.00000E+01
15 3.73480E+01 4.00000E+00 1.50000E+01 1.00000E+00
16 4.46670E+01 4.00000E+00 1.50000E+01 1.00000E+01
17 1.14478E+02 9.00000E+00 1.10000E+01 6.00000E+00
18 1.59947E+02 9.00000E+00 1.10000E+01 1.00000E+00
19 9.05980E+01 9.00000E+00 1.10000E+01 1.00000E+01
20 1.10351E+02 9.00000E+00 1.30000E+01 6.00000E+00
21 1.63141E+02 9.00000E+00 1.30000E+01 1.00000E+00
22 1.15427E+02 9.00000E+00 1.30000E+01 1.00000E+01
23 1.46047E+02 9.00000E+00 1.50000E+01 6.00000E+00
24 1.48498E+02 9.00000E+00 1.50000E+01 1.00000E+00
25 1.15379E+02 9.00000E+00 1.50000E+01 1.00000E+01
26 3.85410E+01 4.00000E+00 1.30000E+01 1.00000E+00
27 3.54170E+01 4.00000E+00 1.30000E+01 1.00000E+01
Min 3.54170E+01 4.00000E+00 1.10000E+01 1.00000E+00
Max 1.63141E+02 9.00000E+00 1.50000E+01 1.00000E+01
Ћ‘ = 7.599E+01 ‘ЉЋ = 4.165E+01 ‘Љ’ = 27
---------------------------------------------------------------
3_
3.54E+01 6.10E+01 8.65E+01 1.12E+02 1.38E+02 1.63E+02
X Y YP KT 4.00E+00+---------+---------+---------+---------+---------+-
4.00000E+00 4.21030E+01 4.18644E+01 9! 122@12
! *
! *
! *
4.50E+00+ *
! *
! *
! *
! *
5.00E+00+ *
! *
! *
! *
! *
5.50E+00+ *
! *
! *
! *
! *
6.00E+00+ *
6.00000E+00 5.65394E+01 5.72698E+01 9! 2 3@* 2 1
! *
! *
! *
6.50E+00+ *
! *
! *
! *
! *
7.00E+00+ *
! *
! *
! *
! *
7.50E+00+ *
! *
! *
! *
! *
8.00E+00+ *
! *
! *
! *
! *
8.50E+00+ *
! *
! *
! *
! *
9.00E+00+ *
9.00000E+00 1.29318E+02 1.27827E+02 9! 1 1 3 @ 11 11
3$
‚€„ 7 “ђЂ‚Ќ…Ќ€џ ЏЋ 1 ”ЂЉ’Ћђ“ Y = 1/ (A + B * X)
A= 3.67374E-02 B= - 3.21271E-03
CKO= 1.49386E+01 CKO%= 35.9 CKO CP. = 9.87110E-01
3_
3.54E+01 6.10E+01 8.65E+01 1.12E+02 1.38E+02 1.63E+02
X Y YP KT 1.10E+01+---------+---------+---------+---------+---------+-
1.10000E+01 7.45396E+01 7.45746E+01 9! 21 11 1 @ 1 1 1
! *
! *
! *
1.14E+01+ *
! *
! *
! *
! *
1.18E+01+ *
! *
! *
! *
! *
1.22E+01+ *
! *
! *
! *
! *
1.26E+01+ *
! *
! *
! *
! *
1.30E+01+ *
1.30000E+01 7.62077E+01 7.61248E+01 9! 11 2 2 @ 1 1 1
! *
! *
! *
1.34E+01+ *
! *
! *
! *
! *
1.38E+01+ *
! *
! *
! *
! *
1.42E+01+ *
! *
! *
! *
! *
1.46E+01+ *
! *
! *
! *
! *
1.50E+01+ *
1.50000E+01 7.72137E+01 7.72615E+01 9! 1 12 2 @ 1 11
3$
‚€„ 3 “ђЂ‚Ќ…Ќ€џ ЏЋ 2 ”ЂЉ’Ћђ“ Y = A + B / X
A= 8.46504E+01 B= - 1.10833E+02
CKO= 4.16392E+01 CKO%= 100.0 CKO CP. = 5.99823E-02
3_
3.54E+01 6.10E+01 8.65E+01 1.12E+02 1.38E+02 1.63E+02
X Y YP KT 1.00E+00+---------+---------+---------+---------+---------+-
1.00000E+00 8.23453E+01 8.30527E+01 9! 2 1 2 1 @* 1 11
! *
! *
! *
1.90E+00+ *
! *
! *
! *
! *
2.80E+00+ *
! *
! *
! *
! *
3.70E+00+ *
! *
! *
! *
! *
4.60E+00+ *
! *
! *
! *
! *
5.50E+00+ *
! *
! *
6.00000E+00 7.70738E+01 7.54823E+01 9! 1 2 1 1 1 @ 1 1 1
! *
6.40E+00+ *
! *
! *
! *
! *
7.30E+00+ *
! *
! *
! *
! *
8.20E+00+ *
! *
! *
! *
! *
9.10E+00+ *
! *
! *
! *
! *
1.00E+01+ *
1.00000E+01 6.85418E+01 6.94259E+01 9! 1 1 1 11 1 @ 1 2
3$
‚€„ 1 “ђЂ‚Ќ…Ќ€џ ЏЋ 3 ”ЂЉ’Ћђ“ Y = A + B * X
A= 8.45668E+01 B= - 1.51408E+00
CKO= 4.12651E+01 CKO%= 99.1 CKO CP. = 1.14916E+00***********************************************************************************************************************
‘ЉЋ = 1.16635160E+01 R = 9.54656513E-01
***********************************************************************************************************************
Ћ ’ Љ ‹ Ћ Ќ … Ќ € џ
11 4.64553E-01 9 - 6.47581E-01 13 - 7.14781E-01 22 - 1.57465E+00 27 - 2.90206E+00 25 - 3.36981E+00
5 - 3.62067E+00 14 4.23122E+00 12 5.33308E+00 7 5.40261E+00 1 5.75339E+00 16 5.77573E+00
10 6.43520E+00 24 6.44152E+00 26 - 7.29922E+00 6 - 8.42069E+00 15 - 9.17674E+00 17 - 1.01399E+01
3 1.18882E+01 4 1.29949E+01 2 - 1.48676E+01 8 - 1.50560E+01 20 - 1.68572E+01 23 1.69392E+01
18 2.28307E+01 21 2.31746E+01 19 - 2.40212E+01
***********************************************************************************************************************
‚€„ ЋЃЋЃ™…ЌЌЋѓЋ “PЂ‚Ќ…Ќ€џ
Y (5) =Y (1) *Y (2) *Y (3) / 5.774019E+03
***********************************************************************************************************************
ѓђЂ”€Љ €‡Њ…Ќ…Ќ€џ ‘ Љ Ћ
3_
0.00E+00 8.33E+00 1.67E+01 2.50E+01 3.33E+01 4.17E+01
OT ‘ђ…„Ќ…ѓЋ ‡ЌЂ-…Ќ€џ ”“ЌЉ-€€ ‘ Љ Ћ = 4.16542733E+01+---------+---------+---------+---------+---------*
! 20 40 60 80 100%
!
ЏЋ‘‹… *** ”ЂЉ’ЋђЂ ‘ Љ Ћ = 1.16635160E+01+ *
3$
***********************************************************************************************************************
”ЋPЊ“‹џP ЊЋ„…‹€
„Ђ’Ђ: 4/5/2008 Ћ’PЂЃЋ’ЂЌЋ 27 ЋЏ›’Ћ‚ ‡ЂЉЂ‡-€Љ:
B=6; E=2; G=1; L=11; P=2;
5=F (1; 2;
3);
------------------------------------------------------------------------------------------------------------------
! N! N! KЋќ””€-€…Ќ’› “PЂ‚Ќ…Ќ€‰! CKO%!! ‘ЉЋ. ‘ђ%!
! ЂPѓ! “P-џ! - ----------------------------------------------------------! - -------------! R! - -------------!
!!! A! B! C! D! 100!! 100!
! - ---! - ---! - -------------! - -------------! - -------------! - -------------! - -------------! - -------------! - -------------!
!!!!!!!!!!
! 1! 7! 3.6737450E-02! - 3.2127082E-03!!! 35.9!! 2.5!
!!!!!!!!!!
! 2! 3! 8.4650388E+01! - 1.1083328E+02!!! 100.0!! 11.7!
!!!!!!!!!!
! 3! 1! 8.4566750E+01! - 1.5140801E+00!!! 99.1!! 224.5!
!!!!!!!!!!
! ***! ***!!!!! 28.0! 9.5465651E-01!!
!!!!!!!!!!
------------------------------------------------------------------------------------------------------------------
DMAX = - 2.4021168E+01 ( 19)
CKO% = 28.00
ЉЋќ””€-€…Ќ’ ”€?…PЂ F = 12.75440; F (001) = 2.52846; F (005) = 1.92428
‚›ЏЋ‹Ќ€‹
***********************************************************************************************************************
Приложение В
Результаты ANETR для σz
B=6; E=2; G=1; L=11; P=1;
6=F (1; 2;3);
1 1.07573E+02 6.00000E+00 1.30000E+01 6.00000E+00
2 9.66540E+01 6.00000E+00 1.30000E+01 1.00000E+00
3 9.77720E+01 6.00000E+00 1.30000E+01 1.00000E+01
4 1.08093E+02 6.00000E+00 1.10000E+01 6.00000E+00
5 1.06188E+02 6.00000E+00 1.50000E+01 6.00000E+00
6 1.07432E+02 6.00000E+00 1.10000E+01 1.00000E+00
7 1.05237E+02 6.00000E+00 1.10000E+01 1.00000E+01
8 1.08531E+02 6.00000E+00 1.50000E+01 1.00000E+00
9 9.76700E+01 6.00000E+00 1.50000E+01 1.00000E+01
10 8.10040E+01 4.00000E+00 1.30000E+01 6.00000E+00
11 8.24950E+01 4.00000E+00 1.10000E+01 6.00000E+00
12 8.09860E+01 4.00000E+00 1.50000E+01 6.00000E+00
13 9.56720E+01 4.00000E+00 1.10000E+01 1.00000E+00
14 9.11610E+01 4.00000E+00 1.10000E+01 1.00000E+01
15 9.55100E+01 4.00000E+00 1.50000E+01 1.00000E+00
16 9.13070E+01 4.00000E+00 1.50000E+01 1.00000E+01
17 1.31362E+02 9.00000E+00 1.10000E+01 6.00000E+00
18 1.25308E+02 9.00000E+00 1.10000E+01 1.00000E+00
19 1.26622E+02 9.00000E+00 1.10000E+01 1.00000E+01
20 1.20552E+02 9.00000E+00 1.30000E+01 6.00000E+00
21 1.18298E+02 9.00000E+00 1.30000E+01 1.00000E+00
22 1.20214E+02 9.00000E+00 1.30000E+01 1.00000E+01
23 1.27884E+02 9.00000E+00 1.50000E+01 6.00000E+00
24 1.29646E+02 9.00000E+00 1.50000E+01 1.00000E+00
25 1.29244E+02 9.00000E+00 1.50000E+01 1.00000E+01
26 9.37530E+01 4.00000E+00 1.30000E+01 1.00000E+00
27 9.45410E+01 4.00000E+00 1.30000E+01 1.00000E+01
Min 8.09860E+01 4.00000E+00 1.10000E+01 1.00000E+00
Max 1.31362E+02 9.00000E+00 1.50000E+01 1.00000E+01
Ћ‘ = 1.063E+02 ‘ЉЋ = 1.588E+01 ‘Љ’ = 27
---------------------------------------------------------------
3_
8.10E+01 9.11E+01 1.01E+02 1.11E+02 1.21E+02 1.31E+02
X Y YP KT 4.00E+00+---------+---------+---------+---------+---------+-
4.00000E+00 8.96032E+01 8.95875E+01 9! 21 @2 211
! *
! *
! *
4.50E+00+ *
! *
! *
! *
! *
5.00E+00+ *
! *
! *
! *
! *
5.50E+00+ *
! *
! *
! *
! *
6.00E+00+ *
6.00000E+00 1.03906E+02 1.03932E+02 9! 12 @1122
! *
! *
! *
6.50E+00+ *
! *
! *
! *
! *
7.00E+00+ *
! *
! *
! *
! *
7.50E+00+ *
! *
! *
! *
! *
8.00E+00+ *
! *
! *
! *
! *
8.50E+00+ *
! *
! *
! *
! *
9.00E+00+ *
9.00000E+00 1.25459E+02 1.25448E+02 9! 1 2 @1 12 1
3$
‚€„ 1 “ђЂ‚Ќ…Ќ€џ ЏЋ 1 ”ЂЉ’Ћђ“ Y = A + B * X
A= 6.08987E+01 B= 7.17218E+00
CKO= 5.17057E+00 CKO%= 32.6 CKO CP. = 1.90567E-02 R= 9.432861E-01
***********************************************************************************************************************
3_
8.10E+01 9.11E+01 1.01E+02 1.11E+02 1.21E+02 1.31E+02
X Y YP KT 1.10E+01+---------+---------+---------+---------+---------+-
1.10000E+01 1.08154E+02 1.08147E+02 9! 1 12@ 2 2
! *
! *
! *
1.14E+01+ *
! *
! *
! *
! *
1.18E+01+ *
! *
! *
! *
! *
1.22E+01+ *
! *
! *
! *
! *
1.26E+01+ *
! *
! *
! *
! *
1.30E+01+ *
1.30000E+01 1.03373E+02 1.03377E+02 9! 1212 @ 21
! *
! *
! *
1.34E+01+ *
! *
! *
! *
! *
1.38E+01+ *
! *
! *
! *
! *
1.42E+01+ *
! *
! *
! *
! *
1.46E+01+ *
! *
! *
! *
! *
1.50E+01+ *
1.50000E+01 1.07441E+02 1.07444E+02 9! 1 1 @21211
3$
‚€„ 11 “ђЂ‚Ќ…Ќ€џ ЏЋ 2 ”ЂЉ’Ћђ“ Y = A * X**2 + B * X + C
A= 1.10477E+00 B= - 2.88999E+01 C= 1.86047E+02
CKO= 4.70435E+00 CKO%= 29.6 CKO CP. = 4.40933E-03 R= 9.512971E-01
***********************************************************************************************************************
3_
8.10E+01 9.11E+01 1.01E+02 1.11E+02 1.21E+02 1.31E+02
X Y YP KT 1.00E+00+---------+---------+---------+---------+---------+-
1.00000E+00 1.07867E+02 1.07860E+02 9! 2 1 @121 1
! *
! *
! *
1.90E+00+ *
! *
! *
! *
! *
2.80E+00+ *
! *
! *
! *
! *
3.70E+00+ *
! *
! *
! *
! *
4.60E+00+ *
! *
! *
! *
! *
5.50E+00+ *
! *
! *
6.00000E+00 1.05126E+02 1.05130E+02 9! 11 1 1@ 21 1 1
! *
6.40E+00+ *
! *
! *
! *
! *
7.30E+00+ *
! *
! *
! *
! *
8.20E+00+ *
! *
! *
! *
! *
9.10E+00+ *
! *
! *
! *
! *
1.00E+01+ *
1.00000E+01 1.05974E+02 1.05978E+02 9! 1 11 @1 1 1
3$
‚€„ 11 “ђЂ‚Ќ…Ќ€џ ЏЋ 3 ”ЂЉ’Ћђ“ Y = A * X**2 + B * X + C
A= 8.42139E-02 B= - 1.13545E+00 C= 2.58845E+00
CKO= 4.55717E+00 CKO%= 28.7 CKO CP. = 5.29090E-03 R= 9.523351E-01
***********************************************************************************************************************
Ћ ’ Љ ‹ Ћ Ќ … Ќ € џ
14 9.33221E-02 6 1.38079E-01 7 - 1.75046E-01 19 - 3.06599E-01 20 - 7.57778E-01 16 9.43032E-01
24 1.53924E+00 8 1.94079E+00 22 - 1.94367E+00 5 2.32755E+00 23 2.50700E+00 13 2.72245E+00
3 - 2.86911E+00 25 3.01911E+00 15 3.26416E+00 18 - 3.50247E+00 4 3.52884E+00 10 - 4.44486E+00
17 5.28129E+00 26 5.57438E+00 21 - 5.74154E+00 2 - 5.86899E+00 9 - 7.03834E+00 11 - 7.72479E+00
1 7.77977E+00 27 8.24425E+00 12 - 8.53008E+00
***********************************************************************************************************************
‚€„ ЋЃЋЃ™…ЌЌЋѓЋ “PЂ‚Ќ…Ќ€џ
Y (6) =Y (1) +Y (2) +Y (3)
***********************************************************************************************************************
ѓђЂ”€Љ €‡Њ…Ќ…Ќ€џ ‘ Љ Ћ
3_
0.00E+00 3.18E+00 6.35E+00 9.53E+00 1.27E+01 1.59E+01
OT ‘ђ…„Ќ…ѓЋ ‡ЌЂ-…Ќ€џ ”“ЌЉ-€€ ‘ Љ Ћ = 1.58833364E+01+---------+---------+---------+---------+---------*
! 20 40 60 80 100%
!
ЏЋ‘‹… 1 ”ЂЉ’ЋђЂ ‘ Љ Ћ = 5.17056730E+00+ *
!
!
ЏЋ‘‹… 2 ”ЂЉ’ЋђЂ ‘ Љ Ћ = 4.70435362E+00+ *
!
!
ЏЋ‘‹… 3 ”ЂЉ’ЋђЂ ‘ Љ Ћ = 4.55717450E+00+ *
3$
***********************************************************************************************************************
”ЋPЊ“‹џP ЊЋ„…‹€
„Ђ’Ђ: 4/5/2008 Ћ’PЂЃЋ’ЂЌЋ 27 ЋЏ›’Ћ‚ ‡ЂЉЂ‡-€Љ:
B=6; E=2; G=1; L=11; P=1;
6=F (1; 2;
3);
------------------------------------------------------------------------------------------------------------------
! N! N! KЋќ””€-€…Ќ’› “PЂ‚Ќ…Ќ€‰! CKO%!! ‘ЉЋ. ‘ђ%!
! ЂPѓ! “P-џ! - ----------------------------------------------------------! - -------------! R! - -------------!
!!! A! B! C! D! 100!! 100!
! - ---! - ---! - -------------! - -------------! - -------------! - -------------! - -------------! - -------------! - -------------!
!!!!!!!!!!
! 1! 1! 6.0898722E+01! 7.1721842E+00!!! 32.6! 9.4328610E-01! 0.1!
!!!!!!!!!!
! 2! 11! 1.1047693E+00! - 2.8899930E+01! 1.8604702E+02!! 29.6! 9.5129706E-01! 3.1!
!!!!!!!!!!
! 3! 11! 8.4213889E-02! - 1.1354500E+00! 2.5884491E+00!! 28.7! 9.5233506E-01! 3.7!
!!!!!!!!!!
------------------------------------------------------------------------------------------------------------------
DMAX = - 8.5300788E+00 ( 12)
CKO% = 28.69
ЉЋќ””€-€…Ќ’ ”€?…PЂ F = 12.14765; F (001) = 2.52846; F (005) = 1.92428
‚›ЏЋ‹Ќ€‹
***********************************************************************************************************************
Приложение Г
Командный файл
/filnam, massiv! Определение имени файла
/title,lab_3! Определение заголовка задачи
/prep7 ! Вход в препроцессор
a=6! h
b=13! R
c=6! r
et,1,92! Задание типа элемента
mp,ex,1, 200000Задание свойств материала
mp,ey,1, 200000
mp,ez,1, 200000
mp,nuxy,1,0.3
mp,nuxz,1,0.3
mp,nuyz,1,0.3
k,1,0,0,0! Построение точек
k,2,60,0,0
k,3,60,-50,0
k,4,15,-50,0
k,5,15,-40,0
k,6,0,-40,0
L,1,2! Построение линий
L,2,3
L,3,4
L,4,5
L,5,6
L,6,1
al,all! Объединение в область
voffst,1,25! Вытягивание на высоту
local,100,0,0,0,-25,0,0,0,0,,! Смещение системы координат
wpcsys,1,100
k,13,0,0,0
k,14,25,0,0
k,15,0,-25,0
L,13,14
L,14,15
L,15,13
al, 19, 20,21
voffst,9,35
local,100,0,0,0,-60,0,0,0,0,,
wpcsys,1,100
k, 19,0,-13,0
k, 20,13,0,0
L,16, 19
L,16, 20
larc, 19, 20,16,13! Построение дуги
al,28,29,30
voffst,14,6
vsbv,2,3! Вычитание объемов
local,100,0,0,0,-54,0,0,0,0,,
wpcsys,1,100
k,24,0,-6,0
k,25,6,0,0
L,21,24
L,21,25
larc,24,25,21,6
al,22,24,25
voffst,10,54
vadd,1,4! Объединение объемов
vsbv,3,2
smrtsize,all,5! Задание размера конечного элемента
vmesh,all! Разбиение на конечные элементы
/solu! Вход в процессор
da,26,all! Закрепление областей
da,8,uy
da,21,uy
da,3,ux
da, 20,ux
sfa,9,,pres,100! Приложение нагрузки
solve! Запуск задачи на решение
/post1! Вход в постпроцессор
plnsol,u,sum,0,1! просмотр напряжений