Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Sb98342

.pdf
Скачиваний:
4
Добавлен:
13.02.2021
Размер:
1.09 Mб
Скачать

MP, 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

Full>
Reduced>

В 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

 

 

 

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]