
Sb98342
.pdfMP, EX, 1, 2E11 – модуль Юнга для материала со справочным номером 1 равен 2Е11;
MP, DENS, 1, 7800 – плотность для этого же материала равна 7800 ед.; MP, KXX, 3,43 – теплопроводность для материала со справочным номе-
ром 3 – 43 ед.
Кроме свойств в Y- и Z-направлениях (по умолчанию совпадают со значением в Х-направлении) многие другие свойства уже заведены, и по умолчанию предполагается, что, например, коэффициент Пуассона (NUXY) – 0.3, модуль сдвига (GXY) – EX/(2*(1+NYXY)) и излучательная способность – 1.
MP |
Для определения температурно-зависимых свойств мате- |
||||
MPTEMP |
риала можно использовать либо команду MP (в комбинации с |
||||
MPTGEN |
командами MPTEMP или MPTGEN), либо команды MPTEMP |
||||
MPDATA |
и MPDATA. Команда МР позволяет определять свойства как |
||||
|
функцию в форме полинома в зависимости от температуры. |
||||
|
Полином может быть линейным, квадратичным, кубичным |
||||
|
или 4-й степени: |
|
|
|
|
|
свойство C +C T+C T 2+C T 3+C T 4 , |
||||
|
0 |
1 |
2 |
3 |
4 |
где Cn – коэффициенты; Т – температура.
Коэффициенты вводятся с использованием полей C0 , C1, C2 , C3 , C4 командой МР. Если определен только коэффициент C0 , свойство материала постоянно; если определены C0 и C1, то свойство меняется линейно, и т. д. Когда задается температурно-зависимое свойство, программа сама вычислит полином в дискретных точках температуры с линейной интерполяцией вне точек экстремума. Для задания свойств 2-го и высшего порядка при определении соответствующих температурных шагов используются команды MPTEMP или MPTGEN, которые должны быть заданы перед командой MP.
Второй метод определения температурно-зависимых свойств материалов заключается в использовании комбинации команд MPTEMP и MPDATA. MPTEMP (MPTGEN) задает серию температур, а MPDATA определяет соответствующие значения свойств материала. Например, следующие команды определяют температурно-зависимую энтальпию для материала 4:
MPTEMP,1, 1600, 1800, 2000, 2325, 2326, 2335, –
заданы 6 температур из положения 1;
MPTEMP, 7, 2345, 2355, 2365, 2374, 2375, 3000, –
21

6 следующих температур из положения 7;
MPDATA, ENTH, 4, 1, 53.81, 61.23, 68.83, 81.51, 81.55, 82.31 –
соответствующие значения ENTH;
MPDATA, ENTH, 4, 7, 84.48, 89.53, 99.05, 112.12, 113.00, 137.40
Если число точек данных свойств и температуры не совпадают, для таблицы свойств функции используются только те положения, которые имеют обе заданные точки. Чтобы определить различный набор температур для следующего свойства материала, следует сначала исправить текущую таблицу температур, используя команду MPTEMP (без аргументов), и затем задать новые температуры (используя дополнительные команды MPTEMP
или MPTGEN).
MPPLOT Команда MPPLOT рисует график зависимости свойств матеMPLIST риала от температуры. Рис. 4.1 показывает нарисованную командой MPPLOT кривую энтальпии как функцию от темпера-
туры. Команда MPLIST печатает свойства материалов.
ENTH
ENTH
100
TEMP
50
1500 |
2000 |
2500 |
3000 TEMP |
Рис. 4.1
О температурно-зависимых свойствах материалов можно сказать следующее:
1. Для изменения точки данных свойств на существующей кривой переопределите, используя команду MPDATA, точку с желаемыми данными на точку с соответствующим номером положения. Например, чтобы изменить значение ENTH в позиции 6 на кривой зависимости «энтальпия – температура» с 82.31 на 83.09, команда должна быть такой:
MPDATA, ENTH, 4, 6, 83.09
22
MPDRES |
Чтобы изменить точку с данной температурой на сущест- |
|
вующей кривой, вам понадобятся две команды: MPTEMP с |
|
соответствующим номером положения для задания новой |
|
температуры и, чтобы связать таблицу с новой температурой |
со свойством материала.
Например: чтобы изменить температуру в позиции 7 на кривой зависимости «энтальпия – температура» с 2345 на 2340, команда должна быть:
MPTEMP, 7, 2340 – меняет позицию 7;
MPDRES, ENTH, 4 – связывает ENTH для материала 4 с новой температурой.
Всякий раз, определяя температурно-зависимое свойство, команда MPDRES немедленно записывает пару данных (температура – свойство) в базу данных. Изменение точки данных температуры может изменить только будущее свойство материала, не влияя на уже существующее. Команда MPDRES изменяет уже записанное в базе данных свойство материала. Два дополнительных поля команды MPDRES позволяют изменить записанное свойство и записать его под новой меткой или новым справочным номером материала.
MPTRES Команда MPTRES позволяет восстанавливать таблицу температур существующих свойств материала из базы данных. В дальнейшем эти же точки данных температуры можно использовать для различных свойств.
2. Температурно-зависимые свойства материалов извлекаются из базы данных для счета по мере формирования матрицы элементов. Программа сперва вычисляет температуру в центре каждого элемента (для тепловых элементов – в точках интегрирования), определяет соответствующее значение свойства материала при помощи линейной интеграции таблицы «свойство – температура» и затем использует это значение для вычисления матрицы элементов. Если температура элементов выпадает из пределов табличных данных, допускается постоянное значение для свойства материала при предельном значении интервала.
MPWRITE Линейные свойства материала (либо температурноMPREAD зависимые, либо постоянные) могут быть записаны из базы CDWRITE данных в кодированный файл (MPWRITE) или читаться из кодированного файла в файл базы данных (MPREAD). Команда CDWRITE, MAT может быть использована для записи
как линейных, так и нелинейных свойств материала в файл.
23

Нелинейные свойства материала
Preprocess Нелинейные свойства материалов – это обычно таблич-
ные данные: пластичность (кривые – при различных зако-
ElemenType |
нах упрочнения), магнитные поля (В–Н кривые), ползучесть, |
RealCons |
расширение, гиперпластичность и т. д. |
Material |
Первый шаг для определения нелинейного свойства ма- |
териала – активизация таблицы данных с помощью команды
|
ТВ. Например: |
|
Linear |
||
|
||
- |
ТВ, ВН, 2 |
|
|
||
TB_Table |
активизируется таблица В-Н для материала 2. |
|
|
||
|
Сейчас, чтобы ввести табличные данные, можно исполь- |
|
|
||
|
-зовать панель для ввода данных или набор команд. Рис. 4.2
TBEdit> |
показывает оба метода определения примера В–Н кривой. |
|
Правый блок показывает, как задавать кривую с помощью |
||
TBTEMP |
||
команд TBPT |
||
|
|
|
Start |
– |
Preproc |
|
Finish |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Main Cmds |
Elemtype |
RealConst |
Mater |
SolMod |
DirecGen |
DOF |
Link |
Other |
Flotran |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Util Cmds |
Linar |
TBTable |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Reference |
T = |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Array Edit |
|
|
|
|
|
R |
|
B |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
150.0 |
0.2100 |
|
TBPT, DEFI, 150, .21 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
300.0 |
0.5500 |
|
TBPT, DEFI, 300, .55 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
460.0 |
0.8000 |
|
TBPT, DEFI, 460, .80 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Insert Row* |
|
|
4 |
640.0 |
0.9500 |
|
TBPT, DEFI, 640, .95 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Insert Col* |
|
|
5 |
720.0 |
1.000 |
|
TBPT, DEFI, 720, 1.0 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Delete Row* |
|
|
6 |
890.0 |
1.100 |
|
TBPT, DEFI, 890, 1.1 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Delete Col* |
|
|
7 |
1020.0 |
1.150 |
|
TBPT, DEFI, 1020, 1.15 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Up/Down* |
|
|
8 |
1280.0 |
1.250 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Left/Right* |
|
|
9 |
1900.0 |
1.400 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(Rate=Page) |
|
|
10 |
0 |
0 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Move To* |
|
|
11 |
0 |
0 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Help |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Save |
Reset |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Quit |
Apply |
Pick inside the BUTTON-MENU |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 4.2
24

TBPLOT Команды TBPLOT и TBLIST позволяют проверить табTBLIST лицу данных с помощью экрана или листинга. Пример графиTBCOPY ка, построенного с помощью команды TBPLOT, показан на
рис. 4.3.
11,5. |
|
|
|
|
T1 |
|
|
|
|
|
|
|
|
|
|
|
9 |
|
|
|
|
8 |
|
11,0. |
|
6 |
7 |
|
|
|
|
|
|
||
|
4 |
5 |
|
|
|
|
3 |
|
|
|
|
00,5. |
2 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
00,0. |
|
|
|
|
|
0 |
500 |
1000 |
1500 |
10002000 |
Рис. 4.3
Некоторые типы элементов могут иметь материал с анизотропными упругими свойствами, которые обычно вводятся в форме матрицы. Заметим, что эти свойства отличаются от свойства анизотропной пластичности, которая требует задавать различные кривые – в различных направлениях. Типы элементов, для которых может быть применимо свойство анизотропной упругости – SOLID64 (трехразмерные анизотропные тела), PLANE13 (двухразмерные тела со связанными полями), SOLID5 и SOLID98 (трехразмерные тела со связанными полями).
ТВ Процедура для задания анизотропных упругих свойств TBDATA материала подобна процедуре задания нелинейных свойств. TBLIST Сначала следует активизировать таблицу данных командой ТВ (с меткой Lab=ANEL) и затем определить члены матрицы коэффициентов упругости, либо при помощи панели, либо командой TBDATA. После заполнения таблицы командой TBDATA, можно проверить правильность введенных данных
командой TBLIST.
4.1.5. Создание геометрии модели
Цель этого этапа – создание конечноэлементной модели (узлов и элементов), которая адекватно описывает геометрию модели.
25

Рис. 4.4
Существует два метода создания конечноэлементной модели: моделирование сплошных тел и прямая генерация. При моделировании сплошных тел описываются геометрические границы модели и затем отдается приказ ANSYS для автоматического построения сетки (геометрия с узлами и элементами). Размер и форму элементов, которые создает программа, можно контролировать. При прямой генерации положение каждого узла и связанных с узлами элементов определяется вручную. При этом доступны несколько удобных операций – копирование модели существующих узлов и элементов, симметричное отражение и т. п.
4.2. Задание источников и граничных условий и получение решений
/SOLU |
На этом этапе команда SOLUTION (/SOLU) используется |
|
для определения типа анализа и его опций, прикладываемых |
|
нагрузок, опций задаваемой пошаговой нагрузки и затем про- |
|
грамма запускается на счет. Далее объясняется, как выполня- |
|
ется каждая из этих задач. |
Выбор типа анализа обосновывается условиями нагружения и ответом. Например, если требуется вычислить собственные частоты и модальные формы, следует выбрать модальный анализ.
26

В ANSYS доступны следующие типы анализа: статика (установившийся режим), анализ переходных процессов, гармонический, модальный, спектральный, устойчивость и подструктуры. Не все типы анализа доступны для любой дисциплины. Например, модальный анализ недоступен для тепловой модели.
Опции анализа позволяют осуществлять настройку типа анализа. Существуют следующие типичные (общие) опции анализа: метод решения, включение (выключение) матрицы жесткости, опции Ньютона–Рафсона.
Solution Для определения типа анализа и его опций можно ис-
пользовать либо контрольную панель либо команду ANTYPE
AnalyTyp |
и связанные с ней опции (TRNOPT, HROPT, MODOPT, |
Masters |
SSTIF, NROPT и т. д.). |
GapCond |
Необходимое меню для выбора подходящего анализа по- |
казано слева. Заметим, что при команде By Panel появляется
|
возможность выбора нового анализа (NewAnaly), рестарта |
|
ByPanel |
||
|
||
ByCommand |
(Restart>) или наращивания расширения (ExtPass>). |
|
|
В большинстве случаев представляется возможным но- |
|
|
||
вый анализ. Рестарт доступен только в статическом (устойчи- |
||
|
||
NewAnalys |
||
вый режим) и в динамических анализах. Наращивание расши- |
||
Restart> |
||
рения – опция, действительная только для некоторых типов |
||
ExpPass> |
||
анализа. |
||
|
||
|
||
Пример контрольной панели с вводом соответствующих |
||
|
||
Static> |
||
команд для структурного кратковременного анализа показан |
||
Modal> |
||
на рис. 4.5 (а – контрольная панель, б – команды). Нужно |
||
Spectrm> |
||
помнить, что дисциплины (структурный, тепловой, магнит- |
||
Harmonic |
||
ный и т. д.) задаются типом элемента, используемого в модели. |
||
Transient |
||
Сначала определяются тип анализа и его опции, следую- |
||
|
щий шаг – приложение нагрузок. Некоторые типы структурно-
го анализа требуют перед этим установки других составляющих – таких, как мастер степеней свободы, условия зазора.
Затем выполняется задание источников. Слово «источник», используемое в этом руководстве, предполагает как внешние, так и внутренние нагрузки, а также граничные условия (ограничения, опоры, описание определения граничных полей).
Нагрузки в ANSYS разделены на 6 категорий: DOF Concstration, «Forces», Surface Loads, Body Loads, Inertia Loads, Coupled-field Loads.
27

Большинство из этих нагрузок может быть применимо либо к геометрической модели «solid model» (точки, линии, поверхности), либо к конечноэлементной модели (узлы и элементы).
Пользователю необходимо знать о двух понятиях, относящихся к нагружению: это шаг нагружения и шаг приращения. Шаг нагружения – это просто конфигурация нагрузок, для которых будет получено решение. В структурном анализе, например, можно на первом шаге нагружения приложить нагрузку ветра и на втором – гравитационную. Шаги нагружения также полезны для деления на несколько сегментов кривой кратковременной нагрузки.
|
|
ByPanel |
By Comand |
|
|
|
|
|
|
|
|
|
|
||
ANSYS |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Main Cmds |
|
|
|
|
|
|
Analysis Options |
|
|
|
|
|
|
||
Util Cmds |
|
|
New, Restart, or Expansion Pass: |
|
New Analysis |
||||||||||
Reference |
|
|
Discipline (based on active DOF) |
|
STRUCTURAL |
||||||||||
Full> |
|
|
|
|
|
|
Analysis type |
|
TRANSIENT |
||||||
Reduced> |
|
|
|
|
|
|
Analysis method |
|
|
|
FULL |
||||
ModeSup> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Prestress or stress stiff. effects |
|
INCLUDE STRESS |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
STIFFNESS |
||||
|
|
Large deformation effects |
|
|
|
|
INCLUDE LARGE |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
DEFORM |
|||
|
|
Newton-Raphson option |
|
|
|
|
PROGRAM CHOOSES |
||||||||
|
|
Newton-Raphson adaptive descent |
DO NOT USE ADAPT DE- |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
SCENT |
|||
|
|
Mass Matrix formulation |
|
|
|
|
DEFAULT ELEM MASS |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
MATRIX |
|||
|
|
Equation solver to be used |
|
|
|
|
FRONTAL EQUATION |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
SOLVER |
|||
|
|
JCG Solver tolerance value |
|
|
|
|
0.1000E-07 |
||||||||
|
|
Difference (in degrees) between absolute zero and |
|
|
|
|
|
|
|||||||
|
|
the temperature system being used |
|
|
|
|
|
|
|||||||
|
Help |
To make change, pick |
buttons |
|
|
RESET |
APPLY |
QUIT |
|||||||
Save |
Reset |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Quit |
Apply |
Working . . . |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
а |
|
|
|
|
|
|
ANTYPE, TRANS, TRNOPT, FULL, SSTIF, ON, NLGEOM, ON,
б
Рис. 4.5
28
Шаг приращения – часть шага нагружения. Они используются главным образом для точности и сходимости в кратковременном и нелинейном анализе. Шаг приращения также известен как шаг, взятый за определенный промежуток времени.
Опции шагов нагружения – число шагов нагружения, время окончания шага нагружения и контроль выхода – могут меняться от шага к шагу. В зависимости от задаваемого типа анализа опции шага нагружения могут быть необходимы или нет.
|
4.3. Инициализация решений и просмотр результатов |
SOLVE |
Команда для начала вычислений. При ее запуске ANSYS |
|
читает информацию для создания модели и прикладываемых |
|
нагрузках из базы данных. Результаты записываются в резуль- |
|
тирующий файл (jobname.RST, jobname.PTH, jobname.RMG) и |
|
базу данных. Различие между организацией результирующего |
|
файла и фала базы данных заключается в следующем: в базу |
|
данных может помещаться только один набор результатов, в то |
|
время как в результирующий файл могут быть записаны все |
|
наборы результатов для всех шагов приращения. |
SOLVE |
Другая действующая команда позволяет решать множе- |
|
ственные шаги нагружения соответствующим способом. |
Когда результаты будут вычислены, для просмотра результатов ANSYS |
|
использует постпроцессор. Доступны два постпроцессора: POST1 и POST26. |
|
/POST1 |
Основной постпроцессор. Используется для просмотра |
|
результатов на 1-м шаге приращения (шаге времени) всей мо- |
|
дели. Команда для ввода – /POST1 имеет силу только на на- |
|
чальном уровне. Есть возможность получить контур экрана, |
|
деформированные формы и табличный список для просмотра |
|
и интерпретации результатов анализа. В POST1 доступны |
|
многие другие возможности – вычисление ошибки, комбина- |
|
ция случаев нагружения, вычисление результирующих данных |
|
и операции доступных путей и т. п. |
/POST26 |
Основной постпроцессор. Используется для просмотра |
|
результатов определенных точек модели для всех временных |
|
шагов. Команда для ввода – /POST26, имеет силу только на |
29

|
|
начальном уровне. Есть возможность получить графики ре- |
|||||||
|
|
зультатов данных в зависимости от времени (или частоты) и |
|||||||
|
|
табличные листинги. Другие возможности POST26 включают |
|||||||
|
|
арифметические вычисления и комплексную алгебру. |
|
|
|||||
5. АНАЛИЗ МАГНИТНОГО ПОЛЯ И ПОТЕРЬ В ДВУХ ДЛИННЫХ |
|||||||||
|
|
ПАРАЛЛЕЛЬНЫХ ПРОВОДНИКАХ |
|
|
|||||
В системе двух длинных прямолинейных шин протекают переменные |
|||||||||
токи, сдвинутые по фазе на 180º друг относительно друга. Вследствие боль- |
|||||||||
шой протяженности проводников рассматривается двумерная задача в декар- |
|||||||||
товой системе координат (см. рисунок). Характеристики системы представ- |
|||||||||
лены в таблице. |
|
|
|
|
|
|
|
||
y |
|
|
|
|
|
|
|
|
|
y3 |
Воздух |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Материал № 1 |
|
|
|
|
Az = 0 + j0 |
|
||
|
|
|
|
|
|
|
|
||
|
Тип элемента |
|
|
|
|
|
|
|
|
|
Plane53 |
|
|
|
|
|
|
|
|
|
Степени свободы |
|
|
|
|
|
|
|
|
|
A |
|
|
|
|
|
|
|
|
y2 |
|
|
|
|
|
|
Шины 1 и 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I2 |
Материалы № 2 и № 3 |
|
||
|
|
a |
|
|
|
|
Тип элемента |
|
|
|
|
|
|
|
|
Plane53 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Степени свободы |
|
|
y1 |
|
|
|
|
|
|
A, VOLT |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b |
|
c |
b |
|
|
|
|
|
|
|
|
|
|
f = 50 Гц |
|
|
|
|
Az = 0 + j0 |
|
|
|
|
|
|
|
y0 |
|
|
|
|
|
|
|
|
|
|
x0 |
|
x1 |
x2 |
x3 |
x4 |
|
x5 |
x |
Номер |
Номер типа |
Описание |
Удельное |
Относительная |
Степени |
||||
материала |
элемента |
материала |
электрическое |
магнитная |
свободы |
||||
|
|
|
|
|
сопротивление, |
проницаемость |
|
|
|
|
|
|
|
|
Ом·м |
|
|
|
|
1 |
|
1 |
Воздух |
|
1E+25 |
1 |
A |
|
|
|
(PLANE53) |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
2 |
Проводник |
2E-8 |
|
1 |
A, V |
|
|
|
(PLANE53) |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
3 |
|
2 |
Проводник |
2E-8 |
|
1 |
A, V |
|
|
|
(PLANE53) |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
30 |
|
|
|
|