Рис. П8.5. Поля эквивалентных напряжений по Мизесу
Для представления полей напряжений, не являющихся инвариантными по отношению к выбранной системе координат, в силу осевой симметрии НДС удобнее перейти к цилиндрической СК. Для этого необходимо изменить СК представления результатов (по умолчанию, глобальная декартова СК, в данном случае ось 0Z глобальной цилиндрической СК совпадает с осью симметрии):
Main Menu General PostProc Options for Outp выбрать в поле [RSYS] Result coord system Global cylindric.
Отобразить радиальные напряжения σr (ðèñ. Ï8.6):
Рис. П8.6. Поля радиальных напряжений
Main Menu |
General PostProc |
Plot Results |
Contour Plot |
Element |
Solu |
Stress |
X-direction SX. |
|
|
|
|
Отобразить окружные напряжения σθ (ðèñ. Ï8.7): |
|
Main Menu |
General PostProc |
Plot Results |
Contour Plot |
Element |
Solu |
Stress |
Y-direction SX. |
|
|
|
Рис. П8.7. Поля окружных напряжений
Для контроля характера распределения напряжений по толщине можно использовать известное решение задачи Ламе для составных толстостенных цилиндров.
Создание анимации результатов расчета в терминах эквивалентных
напряжений по Мизесу: |
|
|
|
Utility Menu |
PlotCtrls |
Animate |
Deformed Results |
Stress von |
Mises SEQV |
OK |
|
|
|
После просмотра анимации выйдите из медиаплеера: |
|
Media player |
File – Exit. |
|
|
|
Для удобства просмотра результатов расчета для элементов группы «контактная поверхность» необходимо выбрать только элементы соответствующего типа (номер 3 согласно порядку ввода в начале сеанса работы):
Utility Menu |
Select Elements By attributes Element type Num |
введите 3 |
OK. |
Отображение контактных напряжений (рис. П8.8):
Main MenuGeneral PostProcPlot ResultsContour PlotElement Solu Contact Total Stress STOT OK.
Рис. П8.8. Поля контактных напряжений
Перед выходом из программы сохранить свойства материала для по-
следующего анализа: |
|
Main Menu |
Preprocessor Material Props Write to File... |
Выход из ANSYS-программы: |
|
Utility Menu |
File Exit… |
|
Выход без сохранения базы данных: Quit – No Save! |
OK. |
Выход с сохранением базы данных: Save Everything |
OK. |
8. Для сравнения проведем анализ в осесимметричной постановке (ось 0Y − ось симметрии) с использованием плоских твердотельных КЭ и линейных контактных КЭ. Если Вы не вышли из программы, задайте новый расчет и определите новое имя проекта:
Utility Menu File Clear & Start New...
Utility Menu → File → Change Jobname...
9. Задание типов и опций конечных элементов. Формирование списка типов конечных элементов: Solid – для твердотельного моделирования, Contact и Target – для создания контактной пары:
Main Menu Preprocessor Element Type Add/Edit/Delete Add Structural
Solid → Quad 8 node 82 → Apply → Contact → 2D Target 169 → Apply → Contact → 3 nd surf 172 → OK.
Main Menu General PostProc.
Для более плавного отображения цветовых переходов внутри КЭ открыть окно настроек вывода:
Main Menu General PostProc Options for Outp в поле [/EFACET] Facets/element edge выбрать 4 facets/edge.
Отобразить эквивалентные напряжения по Мизесу (рис. П8.10) и сравнить с ранее полученным результатом (см. рис. П8.5).
Рис. П8.10. Поля эквивалентных напряжений по Мизесу
Рис. П8.11. Поля радиальных напряжений
Рис. П8.12. Поля окружных напряжений
Отобразить радиальные (рис. П8.11) и окружные (рис. П8.12) напряжения и сравнить с ранее полученными результатами (см. рис. П8.6,
Ï8.7). |
|
17. Выход из ANSYS-программы: |
|
Utility Menu File Exit … |
|
Выход без сохранения базы данных: Quit – No Save! |
OK. |
Выход с сохранением базы данных: Save Everything |
OK. |
Приложение 9
Постановка условий совместности деформаций для несовместного разбиения на примере задачи изгиба бруса*
Постановка задачи
Рассматривается задача исследования НДС бруса сечением прямоугольной формы при изгибе. Особенность задачи: обеспечение совместности деформаций на стыке подобластей с различной плотностью сети конечных элементов при несовпадении граничных узлов. Такая ситуация может возникнуть, например, при импорте и сборке сеточных моделей элементов конструкции.
1. Задание типа конечного элемента:
Main Menu Preprocessor Element Type Add/Edit/Delete Add Structural Solid Brick 20 node 95 OK Close.
2. Задание свойств материала (сталь) в единицах СИ. Задание значе- ний модуля упругости (ЕХ, Па) и коэффициента Пуассона (Poisson s Ratio):
Main Menu Preprocessor Material Props Constant Isotropic OK
EX |
2.1E11 |
Poisson s Ratio (minor) |
0.3 |
3. Создание геометрической модели:
Main Menu Preprocessor Modeling Create Volumes Block By Dimensions.
Введите следующие значения: |
|
X1, X2 |
0 |
1.0 |
|
Y1, Y2 |
0 |
2.0 |
|
Z1, Z2 |
0 |
1.0 |
OK |
Скопировать созданное тело со смещением вдоль оси 0Y:
Main Menu Preprocessor Modeling Copy Volumes Pick All.
Ввести в поле DY значение 2.0 OK.
* Автор – К. Жеков, технический специалист CAD-FEM GmbH.
|
|
|
|
|
|
|
|
4. Создание сеточной модели. |
|
|
|
|
Перед созданием сеточной модели отобразить линии: |
|
Utility Menu |
Plot |
Lines. |
|
|
|
|
Показать изометрию: |
|
|
|
|
|
Utility Menu |
Plot Ctrls |
Pan, Zoom, Rotate |
Obliq. |
|
Назначить число разбиений КЭ на линиях: |
|
|
|
Main Menu |
Preprocessor |
Meshing |
Mesh |
|
Lines Picked Lines; |
– выбрать два перпендикулярных ребра (линии) в основании нижне- |
го параллелепипеда |
Apply |
|
|
|
|
NDIV No. of element divisions |
|
1 |
|
|
Space Spacing Ratio |
|
|
1 |
Apply |
|
– выбрать два перпендикулярных ребра (линии) в основании верхне- |
го параллелепипеда |
Apply |
|
|
|
|
NDIV No. of element divisions |
|
3 |
|
|
Space Spacing Ratio |
|
|
1 |
Apply |
|
– выбрать любое длинное ребро верхнего параллелепипеда |
Apply |
NDIV No. of element divisions |
|
3 |
|
|
Space Spacing Ratio |
|
|
1 |
Apply |
|
– выбрать любое длинное ребро нижнего параллелепипеда |
Apply |
NDIV No. of element divisions |
|
2 |
|
|
Space Spacing Ratio |
|
|
1 |
OK |
|
Построение сеточной модели с данными параметрами и использова- |
нием брикетной формы КЭ (рис. П9.1, à): |
|
|
|
|
Main Menu |
Preprocessor |
Meshing |
Mesh |
|
Volumes |
Mapped |
4to 6 Sided Pick All.
5.Описание условий совместности деформаций на стыке двух параллелепипедов. Для этого необходимо отобрать узлы нижней грани верхнего параллелепипеда и прилегающие к ним элементы нижнего параллелепипеда.
При отборе узлов используются операции пересечения разнородных множеств из группы операций селекции объектов:
Utility MenuSelect EntitiesVolumes By Num Pick Apply.
Выбрать верхний параллелепипед OK. Отобразить поверхности:
Utility Menu Plot Areas.
Отобрать и отобразить все поверхности селектированного тела:
Utility Menu |
Select |
Entities |
Areas |
Attached to |
Volumes |
Apply, |
|
|
|
|
|
|
|
|
|
Utility Menu → Plot → Replot. |
|
|
|
|
|
|
Из этих поверхностей отобрать нижнюю: |
|
|
|
|
Utility Menu |
Select |
Areas |
By Num Pick |
|
Reselect |
Apply. |
Выбрать нижнюю поверхность (основание параллелепипеда) OK. |
Utility Menu |
Plot |
Replot. |
|
|
|
|
|
|
Отобрать все узлы на этой поверхности: |
|
|
|
|
Utility Menu |
Select |
Nodes |
Attached to |
|
Areas, All |
From |
Full Apply. |
|
|
|
|
|
|
|
|
|
Для контроля отобразить узлы: |
|
|
|
|
|
Utility Menu |
Plot |
Nodes. |
|
|
|
|
|
|
Нижний параллелепипед разбит на два КЭ. Отобрать верхний из |
íèõ: |
|
|
|
|
|
|
|
|
|
Utility Menu |
Select |
Elements |
By Num Pick |
From Full |
OK. |
Выбрать верхний элемент |
OK. |
|
|
|
|
|
Применить к отобранным объектам функцию согласования (совме- |
стности) деформаций: |
|
|
|
|
|
|
|
|
Main Menu |
Preprocessor |
Coupling/Seqn |
|
Adjacent Regions |
OK. |
Восстановите активность всех объектов: |
|
|
|
|
Utility Menu |
Select |
Everything. |
|
|
|
|
|
Отобразить поверхности: |
|
|
|
|
|
|
Utility Menu |
Plot |
Areas. |
|
|
|
|
|
|
6. Приложение граничных условий: |
|
|
|
|
|
– заделка по верхней поверхности: |
|
|
|
|
|
Main Menu |
Preprocessor (Solution) |
Loads |
Apply |
Structural |
Displacement |
On Areas |
выбрать верхнюю грань OK |
All DOF |
OK; |
|
|
|
|
|
|
|
|
|
– давление на двух смежных боковых гранях нижнего параллелепипеда (косой изгиб):
Main MenuPreprocessor (Solution)LoadsApplyPressureOn Areas выбрать любые две смежные боковые грани нижнего параллелепипеда OK в поле Value ввести значение 1.Е7 OK.
7. Запуск решателя:
Main Menu Solution Solve Current LS