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

ANSYS Mechanical

.pdf
Скачиваний:
773
Добавлен:
31.08.2019
Размер:
61.81 Mб
Скачать

vk.com/club152685050 | vk.com/id446425943

Научно-исследовательский центр СтаДиО

Инв.№ 2009-РААСН-01/1

 

 

“Утверждаю”

 

Генеральный директор ЗАО НИЦ СтаДиО

 

 

 

А.М.Белостоцкий

 

 

“___”

___________2009г.

Верификационный отчет

по программному комплексу ANSYS Mechanical

Том 1

Руководитель работы докт.техн.наук, профессор А.М.Белостоцкий

Ответственные исполнители

С.И. Дубинский А.А. Аул

Исполнители

А.И. Нагибович И.Н. Афанасьева О.А. Козырев А.С. Павлов

Москва, 2009

vk.com/club152685050 | vk.com/id446425943

 

ANSYS Mechanical. Верификационный отчет. Том 1

 

 

 

Содержание

 

1

Список сокращений.................................................................................................................

5

2

Введение.....................................................................................................................................

7

3

Назначение и область применения ПС................................................................................

9

4

Описание ПС...........................................................................................................................

11

 

4.1

Типы решаемых задач (виды расчетов) .....................................................................

11

 

4.2

Верифицируемые возможности ПК ANSYS Mechanical ...........................................

11

4.3Реализация верифицируемых типов решаемых задач (виды расчетов) в ПК

ANSYS Mechanical......................................................................................................................

13

4.3.1

Линейная статика ....................................................................................................

13

4.3.1.1

Прямой разреженный метод..............................................................................

14

4.3.1.2 Итерационный метод сопряженных градиентов с предобуславливанием....

15

4.3.2

Стационарные задачи теплопроводности (фильтрации).................................

16

4.3.3 Собственные частоты и формы колебаний........................................................

16

4.3.3.1 Метод итераций в подпространстве..................................................................

17

4.3.3.2 Прямой блочный метод Ланцоша .....................................................................

19

4.3.3.3

Метод Ланцоша-PCG..........................................................................................

19

4.3.4

Линейная устойчивость..........................................................................................

19

4.3.5

Гармонический анализ (установившиеся вибрации) .......................................

20

4.3.5.1 Метод суперпозиций форм колебаний .............................................................

21

4.3.6

Спектральный динамический (сейсмический) анализ ....................................

22

4.3.6.1 Спектральный отклик с возбуждением по “платформенной” схеме.............

23

4.3.6.2

Спектральная плотность мощности, СПМ.......................................................

24

4.3.7

Прямой динамический расчет...............................................................................

25

4.3.7.1

Метод Ньюмарка.................................................................................................

27

4.3.7.2

Метод Хильбера–Хьюза.....................................................................................

29

4.3.8

Геометрическая нелинейность (включая устойчивость).................................

29

4.3.8.1

Большие деформации.........................................................................................

30

4.3.8.2

Большие повороты..............................................................................................

31

4.3.8.3

Ужесточение при нагружении...........................................................................

31

4.3.9Физическая нелинейность (металл, железобетон, грунт, резина, дерево,…)33

4.3.9.1

Пластичность материалов, не зависящая от времени .....................................

34

4.3.9.2

Ползучесть...........................................................................................................

41

4.3.9.3

Вязкопластичность .............................................................................................

42

4.3.9.4

Нелинейная упругость........................................................................................

43

4.3.9.5

Гиперупругость...................................................................................................

44

4.3.9.6

Моделирование бетона.......................................................................................

47

4.3.10 Нелинейные статические и динамические задачи............................................

48

4.3.10.1

Метод Ньютона-Рафсона ...................................................................................

49

4.3.10.2

Метод окаймляющих дуг...................................................................................

50

4.3.11 Нестационарный нелинейный анализ (теория поля) .......................................

52

4.3.11.1

Метод трапеций (Хьюза)....................................................................................

52

4.3.12

Суперэлементные алгоритмы................................................................................

53

4.3.12.1

Статические задачи.............................................................................................

53

4.3.12.2

Динамические задачи.........................................................................................

54

4.3.13

Контактные задачи..................................................................................................

56

4.3.13.1

Тип контакта “узел–поверхность”.....................................................................

56

4.3.13.2

Тип контакта “узел–узел” ..................................................................................

57

4.3.13.3

Контактные алгоритмы ......................................................................................

58

4.3.14

Оптимизация.............................................................................................................

60

4.3.14.1

Метод аппроксимации подзадачи.....................................................................

60

4.3.14.2

Метод первого порядка......................................................................................

62

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), МГСУ (niccm@mgsu.ru), 2009

3

vk.com/club152685050 | vk.com/id446425943

 

ANSYS Mechanical. Верификационный отчет. Том 1

 

 

4.4 Библиотека конечных элементов.................................................................................

64

 

4.4.1

Элементы связей ......................................................................................................

64

 

4.4.2 Балочные (стержневые) элементы .......................................................................

64

 

4.4.3 Плоские (двумерные) элементы............................................................................

65

 

4.4.4

Оболочечные элементы ..........................................................................................

65

 

4.4.5

Объемные элементы................................................................................................

66

 

4.4.6

Суперэлементы.........................................................................................................

66

 

4.4.7

Контактные элементы.............................................................................................

67

 

4.4.8

Специальные элементы..........................................................................................

67

 

4.4.9

Комбинированные элементы.................................................................................

68

 

4.5 Создание геометрических и расчетных моделей.......................................................

68

 

4.5.1 Обзор методов создания моделей..........................................................................

68

 

4.5.2

Системы координат.................................................................................................

69

 

4.6 Приложение нагрузок, проведение вычислений, операции с результатами.......

70

 

4.6.1

Приложение нагрузок..............................................................................................

70

 

4.6.2 Шаги нагрузки, промежуточные шаги и итерации решения..........................

71

 

4.6.3 Сущность «времени» при приложении нагрузок ..............................................

71

 

4.6.4

Проведение расчета .................................................................................................

72

 

4.6.5

Обзор постпроцессоров...........................................................................................

72

5

Матрицы верификации ........................................................................................................

74

 

5.1 Матрицы верификации (Том 2. Верификационные примеры из ANSYS

 

 

Verification Manual)....................................................................................................................

74

 

5.2 Матрицы верификации (Том 3. Задачи, предложенные экспертами, и

 

 

разработанные авторами отчета)............................................................................................

81

6

Описание верификационных примеров............................................................................

86

7

Результаты верификации ПС..............................................................................................

88

8 Проект Свидетельства РААСН о верификации ПС (с приложением) ........................

89

9 Список используемой литературы (Перечень источников)..........................................

95

Приложение 1 (том 2). Верификационные примеры из ANSYS Verification Manual

Приложение 2 (том 3). Задачи, предложенные экспертами и разработанные авторами отчета

Приложение 3 (том 4). Опыт применения ANSYS в проектной и экспертной практике

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), МГСУ (niccm@mgsu.ru), 2009

4

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 1

1 Список сокращений

 

 

Сокращение

Значение

 

 

 

Матрицы, векторы

 

 

[B]

Матрица связи деформаций и перемещений

 

 

[C]

Матрица демпфирования

 

 

[D]

Матрица упругих констант

 

 

[I]

Единичная матрица

 

 

[M]

Матрица масс (инерции)

 

 

[K]

Матрица жесткости

 

 

[Kt]

Матрица теплопроводности

 

 

[L]

Нижняя матрица треугольного вида

 

 

[S]

Матрица геометрической жесткости (stress stiffness matrix)

 

 

[R]

Матрица поворота

 

 

[TR]

Матрица преобразования из ЛСК (локальной системы

 

координат) в ГСК (глобальную систему координат)

 

 

[U]

Верхняя матрица треугольного вида; матрица формы

 

 

{F}

Вектор нагрузок

 

 

P, {P}

Давление, вектор давления

 

 

{Q}

Вектор теплового потока

 

 

{T}

Вектор температур

 

 

u, v, w, {u}

Перемещения, вектор перемещений (заданных и искомых)

 

 

{φ}, {ψ}, {e}

Собственные векторы

 

 

{ω}

Круговые частоты собственных колебаний

 

 

{f}, {λ}

Частоты собственных колебаний, собственные значения

 

 

{α}

Вектор переноса поверхности текучести

 

 

 

Переменные

 

 

DV

Переменные проекта

 

 

m, n, M, N

Количество

 

 

SV

Переменные состояния

 

 

t

Время, толщина

 

 

ε

Деформации

 

 

σ

Напряжения

 

 

x, y, z

Координаты конечных элементов

 

 

X, Y, Z

Узловые координаты (обычно в декартовой ГСК)

 

 

V

Объем

 

 

δU

Вариация работы внутренних сил

 

 

δV

Вариация работы внешних сил

 

 

κ

Пластичная работа

 

 

Λ

Растяжение

 

 

θ

Параметр интегрирования по времени

 

 

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), МГСУ (niccm@mgsu.ru), 2009

5

vk.com/club152685050 | vk.com/id446425943

 

ANSYS Mechanical. Верификационный отчет. Том 1

 

 

 

 

 

 

Сокращение

 

 

Значение

 

 

 

 

 

 

 

 

 

Константы материалов

 

 

 

 

 

 

E

 

 

Модуль упругости (Юнга)

 

 

 

 

 

 

ν

 

 

Коэффициент Пуассона

 

 

 

 

 

 

α

 

 

Коэффициент температурного расширения

 

 

 

 

 

 

β

 

 

Демпфирование, связанное с жесткостью

 

 

 

 

 

 

ξ

 

 

Коэффициент затухания

 

 

 

 

 

Надстрочные и подстрочные индексы

 

 

 

 

 

 

ac

 

 

Узловые эффекты, обусловленные ускорением

 

 

 

 

 

 

a

 

 

Внешние (приложенные) эффекты (например, внешняя

 

 

 

 

нагрузка)

 

 

 

 

 

 

cr

 

 

Ползучесть (creep)

 

 

 

 

 

 

d

 

 

Динамический

 

 

 

 

 

 

e

 

 

Элементная основа в глобальных координатах

 

 

 

 

 

 

el

 

 

Упругий

 

 

 

 

 

 

f

 

 

Не закрепленные (свободные) (степени свободы)

 

 

 

 

 

 

g

 

 

Внутреннее теплообразование(internal heat generation)

 

 

 

 

 

 

i

 

 

Номер итерации для достижения равновесия

 

 

 

 

 

 

 

 

 

Элементная основа в элементных координатах

 

 

 

 

 

 

m

 

 

Основные (главные)(Master)

 

 

 

 

 

 

max

 

 

Максимальное значение

 

 

 

 

 

 

n

 

 

Номер подшага (шага по времени)

 

 

 

 

 

 

nd

 

 

Эффекты, примененные непосредственно к узлу

 

 

 

 

 

 

nr

 

 

Восстанавливающий эффект

 

 

 

 

 

 

pl

 

 

Пластичность

 

 

 

 

 

 

pr

 

 

Давление

 

 

 

 

 

 

r

 

 

Реактивные эффекты (например, реактивные усилия);

 

 

 

 

закрепленные (степени свободы)

 

 

 

 

 

 

s

 

 

Зависимые (slave), псевдостатические

 

 

 

 

 

 

sw

 

 

Разбухание (swelling)

 

 

 

 

 

 

t, th

 

 

Термический, температурный (thermal)

 

 

 

 

 

 

tot

 

 

Полный

 

 

 

 

 

 

^

 

 

Обозначение редуцированной матрицы или вектора

 

 

 

 

 

 

.

 

 

Дифференцирование по времени

 

 

 

 

 

 

 

 

 

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), МГСУ (niccm@mgsu.ru), 2009

6

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 1

2 Введение

Численное (в обиходе – “компьютерное”) моделирование статического, температурного и динамического напряженно-деформированного состояния зданий, сооружений и конструкций без каких-либо существенных ограничений основано на различных формулировках метода конечных элементов (МКЭ), реализованного в известных программных комплексах. В отличие от машиностроения, где широко распространены т.н. “тяжелые” универсальные зарубежные (ABAQUS, ANSYS, NASTRAN, COSMOS/M, …) и

отечественные (СТАДИО, ЗЕНИТ,…) пакеты, в промышленном, гражданском и транспортном строительстве расчеты, как правило, проводятся специализированными строительно-ориентированными пакетами (Лира, SCAD, MicroFE, Stark ES, Robot Millennium и др.).

Несомненными достоинствами последних являются большая практика применения в массовом проектировании, наличие сертификатов Госстроя РФ, подтверждающих учет методик СНиП, умеренная цена, делающая их доступными даже небольшим проектным организациям. Применяемые алгоритмы прошли испытание временем и их возможностей достаточно для типовых практических задач.

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

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

Трудоемкость подготовки моделей и совершенствование алгоритмов автоматической генерации сеток.

Геометрическая нелинейность в поведении большепролетных сооружений

Учет разнообразного физической нелинейности поведения строительных материалов (металлы и сплавы, бетон и железобетон, каменная и кирпичная кладка, деревянные и резиноподобные, грунты и скальные массивы)

Необходимость совместного учета грунта и конструкции в нелинейной

постановке

Сложный характер ветрового нагружения и его большая динамическая составляющая

Необходимость анализа “нештатных” ситуаций в высоконелинейной динамической постановке с учетом разрушения (например, прогрессирующего разрушения)

Необходимость решения нетиповых задач теплообмена и вентиляции

Оптимизация проектных параметров нетиповых конструкций

Программные продукты ANSYS разрабатываются с конца 60-гг. прошлого столетия. Компания ANSYS, Inc. основана в 1970 г. Дж. Суонсоном (Dr. John Swanson) и до 1994 года именовалась как SAS Inc. (Swanson Analysis Systems).

Первая коммерческая версия выпущена в 1971 году. На 2009 год насчитывается около 200 тысяч коммерческих и более 200 тысяч университетских инсталляций по всему миру. В СНГ ANSYS приобрели и используют уже около 300 организаций.

Программный комплекс аттестован Госатомнадзором РФ (рег. номер ПС в ЦОЭП при РНЦ КИ №490 от 10.09.2002, рег. номер паспорта аттестации №145 от 31.10.2002, Рис. 2.1). Ранее и позже сертифицирован многими авторитетными зарубежными органами (ISO-9001, ISO 9000-3, British standard BS 5750, Lloyd’s Register’s software certification, NAFEMS QA certification и др.). Сертификат Lloyd’s Register’s эквивалентен сертификации Морским Регистром РФ.

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), МГСУ (niccm@mgsu.ru), 2009

7

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 1

Рис. 2.1 Аттестационный паспорт Госатомнадзора РФ (2002 г.)

Основная цель разработанного верификационного отчета заключается в демонстрации широкого круга конечных элементов и других возможностей ПК ANSYS Mechanical в задачах, имеющих отношение к строительной отрасли и имеющих “классическое” или легко получаемое теоретическое решение. Кроме того, тесное согласование решений ПК ANSYS Mechanical с теоретическими результатами в данном отчете позволяют пользователю быть уверенным в достоверности результатов ПК ANSYS Mechanical. Была предпринята попытка включить большинство типов элементов и основные возможности ПК ANSYS Mechanical в набор тестовых примеров. Эти примеры могут служить основанием для дополнительной проверки и уточнения возможностей ПК ANSYS самим пользователем.

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), МГСУ (niccm@mgsu.ru), 2009

8

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 1

3 Назначение и область применения ПС

ANSYS Mechanical – универсальный многоцелевой программный комплекс, позволяющий решать стационарные и нестационарные задачи теории поля (теплопроводность, фильтрация и др.), выполнять расчет статического, температурного и динамического напряженно-деформированного состояния и оптимизации произвольных пространственных комбинированных конструкций, зданий и сооружений с учетом эффектов физической, геометрической, структурной (контакты с трением) и генетической (история возведения и нагружения) нелинейностей на основе метода конечных элементов.

Широкий спектр возможностей ПК ANSYS позволяет инженерам-исследователям моделировать произвольные пространственные комбинированные (массивно-оболочечно- плитно-стержневые из различных материалов-сред) строительные конструкции, здания и сооружения, включая:

наземные и подземные, высотные и большепролетные, монолитные и панельные, вантовые;

металлические (стальные, чугунные, сплавы), бетонные, железобетонные, сталежелезобетонные, каменные и кирпичные, деревянные, резиновые и резинометаллические, грунтовые;

сложные конструктивные узлы.

Основной задачей конечноэлементного расчета конструкций является исследование ее поведения в условиях приложения разнообразных внешних нагрузок и воздействий. ПК ANSYS дает возможность приложения следующих видов нагрузок (применительно к строительной тематике):

заданные тепловые и фильтрационные потоки (теория поля);

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

пульсационная (динамическая) составляющая ветровой нагрузки;

сейсмические, заданные трехкомпонентными спектрами ответа и акселерограммами;

силовые динамические, заданные временной реализацией;

вибрационные, заданные амплитудами и частотами воздействия;

случайные динамические, заданные спектрально (PSD).

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

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

Ниже приведены общие соображения по требуемому объему оперативной памяти для расчетов по ПК ANSYS:

1) Общим правилом для решателей ANSYS может быть 1 GB оперативной памяти на миллион степеней свободы (DOFs). Это правило подходит как для итерационных решателей (таких как PCG), так и для прямых (таких как direct sparse solver). Трехмерные модели, состоящие из элементов второго порядка, требуют больше памяти при расчете (требования к памяти возрастают в 2–3 раза, в случае применения sparse solver), а также несколько широко используемых элементов в ANSYS (к примеру, SOLID92 и SOLID95) используют наоборот

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), МГСУ (niccm@mgsu.ru), 2009

9

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 1

гораздо меньше памяти при использовании итерационного решателя PCG и опции MSAVE,ON (сейчас она используется по умолчанию везде, где это только возможно).

2)Требования итерационного решателя к системе ввода/вывода те же, что и к памяти за исключением некоторых видов расчетов, включающих множество шагов нагружения, которые могут сгенерировать очень большие файлы результатов. Запуск “Sparse solver” требует 10 GB файлов на миллион степеней свободы. Это число может увеличиваться в 2–3 раза при использовании 3D моделей с элементами высшего порядка. Блочный метод Ланцоша (Block Lanczos) является, пожалуй, наиболее требовательным к вычислительным ресурсам и операциям ввода/вывода решателем ПК ANSYS.

3)Средства ввода/вывода не могут ускоряться процессорами. Это проблема является также проблемой управления итерационными решателями и может привести к потере производительности при использовании решателей Sparse. Большинство современных операционных систем, включая Windows, используют некоторые объемы оперативной памяти для большего кэша файлов. Таким образом, аппаратные системы с большим объемом памяти, в общем случае более предпочтительны.

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

5)Уже существуют примеры решения задач с использованием ПК ANSYS, где было задействовано от 8 до 16 процессоров с более чем 100 Гб оперативной памяти. На таких аппаратных системах ПК ANSYS может решать модели размером в более чем 100 миллионов степеней свободы с лидирующей в своем классе инженерных программных комплексов скоростью. В меньших системах ПК ANSYS показывает предельно допустимую для них пиковую производительность, только если они правильно сконфигурированы, что достигается имеющимися в программном комплексе средствами оптимального управления памятью.

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), МГСУ (niccm@mgsu.ru), 2009

10

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 1

4 Описание ПС

4.1Типы решаемых задач (виды расчетов)

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

1)решение стационарных и нестационарных задач теории поля (теплопроводность, фильтрация и т. п.);

2)статические расчеты конструкций (в т.ч., “монтаж” и расчет на прогрессирующее обрушение);

3)расчет собственных частот и форм колебаний в энергетически значимом частотном диапазоне и решение задач линейной устойчивости (частичная проблема собственных значений);

4)гармонический анализ (установившиеся вынужденные колебания);

5)решение задач линейно-спектральной теории сейсмостойкости (на спектры ответа);

6)спектральный динамический анализ (с разложением по собственным формам колебаний);

7)спектральный анализ случайных колебаний;

8)переходные динамические процессы (прямое интегрирование уравнений

движения);

9)решение задач оптимизации геометрической формы и структуры.

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

4.2Верифицируемые возможности ПК ANSYS Mechanical

Верифицируемые возможности ПК ANSYS Mechanical можно структурно записать следующим образом:

9“Прочностные” задачи

Статика

Собственные частоты и формы колебаний

Гармонический анализ

Переходные процессы

Спектральный анализ и случайные колебания

Линейная устойчивость конструкций

9Нелинейные факторы

Геометрические нелинейности

Большие деформации

Большие прогибы

Напряженная жесткость

Физическая нелинейность

Пластичность

Ползучесть

Вязкопластичность

Гиперупругость…

9Контактные задачи

Типы контакта

Элементы “поверхность-поверхность”

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), МГСУ (niccm@mgsu.ru), 2009

11