Знакомство с системой ИСПА

[Борман 2 396]

ключ Разреженный метод

Нет такого ключа.

[Борман 2 396]

А теперь проведите расчет разреженным прямым методом.

Время работы 0 час 0 мин 17 сек

После завершения запустите только прямой решатель.

Время работы 0 час 0 мин 22 сек

И куда смотреть ?

[Борман 2 396]

Чтобы понять что происходит приведите полный код.

Скрытый текст: Сначала прямой

В В О Д Д А Н Н Ы Х

И Н И Ц И А Л И З А Ц И Я

Читаем таблицу ПАРАМЕТРЫ

Читаем таблицу УСКОРЕНИЕ

Читаем таблицу ДИАГРАММА_ЕТ

Читаем таблицу УЗЛЫ

Читаем таблицу ХАРАКТЕРИСТИКИ_СЕЧ

Читаем таблицу ГЕОМЕТРИЯ_СЕЧ

Читаем таблицу МАТЕРИАЛЫ

Читаем таблицу ЭЛЕМЕНТЫ

Читаем таблицу ЗАКРЕПЛЕНИЕ

Читаем таблицу УРАВНЕНИЯ

Читаем таблицу УРАВНЕНИЯ_СТЕРЖНИ

Читаем таблицу ТОЧЕЧ_ВЕС

Читаем таблицу СИЛЫ

Читаем таблицу СИЛЫ

Читаем таблицу ДАВЛЕНИЕ

Читаем таблицу ТЕМПЕРАТУРЫ

Читаем таблицу ТЕМПЕРАТУРЫ

Читаем таблицу ПОТОК_ОБМЕН

Читаем таблицу ПОТОК_ОБМЕН

Время работы 0 час 0 мин 0 сек

Б Е З О Ш И Б О К

Д И А Г Н О С Т И К А М О Д Е Л И

И Н И Ц И А Л И З А Ц И Я

Р а с ч е т г е о м е т р и ч е с к и х

х а р а к т е р и с т и к с е ч е н и й

Д И А Г Н О С Т И К А М О Д Е Л И

Диагностика сечений

Диагностика материала

Диагностика ортотропного материала

Диагностика элементов

Диагностика сеток контакта

Диагностика закрепления

Диагностика уравнений

Диагностика нагружения

Диагностика нагружения края давлением

Д И А Г Н О С Т И К А Б Е З О Ш И Б О К

С о с т а в л е н и е т а б л и ц ы

Всего обработано элементов 500

О п т и м и з а ц и я ш и р и н ы л е н т ы

обработано узлов 1000

обработано узлов 2000

обработано узлов 3000

обработано узлов 3035

Начальный узел оптимизации 4329

С о с т а в л е н и е т а б л и ц ы

с т е п е н е й с в о б о д ы

Всего обработано элементов 500

Количество внешних степеней свободы 615

Количество степеней свободы 8490

Время работы 0 час 0 мин 3 сек

Б Е З О Ш И Б О К

О Б Р А Б О Т К А Н А Г Р У З К И

И Н И Ц И А Л И З А Ц И Я

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

м а с с

О б р а б о т к а с о б с т в е н н о г о

в е с а

Вес модели 1.681452E+003

Вес модели с учетом

точечного веса 1.681452E+003

К о о р д и н а т ы ц е н т р а

т я ж е с т и

X = 12.205887 Y = 24.499994 Z = 100.000023

Время работы 0 час 0 мин 0 сек

Б Е З О Ш И Б О К

С О С Т А В Л Е Н И Е М А Т Р И Ц Ы Ж Е С Т К О С Т И

И Н И Ц И А Л И З А Ц И Я

С о с т а в л е н и е м а т р и ц

ж е с т к о с т и э л е м е н т о в

Всего обработано 500 элементов

С о с т а в л е н и е м а т р и ц

з а д а н н ы х п о л ь з о в а т е л е м

С о с т а в л е н и е г л о б а л ь н о й

м а т р и ц ы ж е с т к о с т и

О б р а б о т а н о 8490 с т р о к

С о с т а в л е н и е о п о р н о й

м а т р и ц ы

О б р а б о т а н о 615 с т р о к

О б р а б о т к а с и л д л я п р о б ы

р а в н о в е с и я

Время работы 0 час 0 мин 4 сек

Реальное время работы 0 час 0 мин 8 сек

Б Е З О Ш И Б О К

Р Е Ш Е Н И Е У Р А В Н Е Н И Й

Матрица жесткости не будет записана на диск

Допуск диагонального члена 1.000000E-012

И Н И Ц И А Л И З А Ц И Я

П е р е п и с ы в а е м м а т р и ц у

ж е с т к о с т и

Прогнозируемое время решения

для Intel I7 2900 МГц 0 час 0 мин 1 сек

переписано 5000 степеней

Всего переписано 8490 степеней

Максимальная ширина ленты 897 в строке 4411

Средняя ширина ленты 601

максимальный диагональный член 2.803872E+008

минимальный диагональный член 6.088548E+006

средний диагональный член 6.242097E+007

Р а з л о ж е н и е м а т р и ц ы

обработано 613 строк

обработано 1236 строк

обработано 2214 строк

обработано 2999 строк

обработано 3504 строк

обработано 4009 строк

обработано 4572 строк

обработано 5113 строк

обработано 5734 строк

обработано 6295 строк

обработано 6858 строк

обработано 7536 строк

обработано 8490 строк

О б р а б о т к а п р а в о й ч а с т и

О б р а т н ы й х о д

О п р е д е л е н и е о п о р н ы х

р е а к ц и й

Время работы 0 час 0 мин 24 сек

Реальное время работы 0 час 0 мин 28 сек

Б Е З О Ш И Б О К

В Ы В О Д Р Е З У Л Ь Т А Т О В

Будут записаны данные для постпроцессора

Будут рассчитаны напряжения в плоских элементах

Будут рассчитаны напряжения в объемных элементах

Будут записаны напряжения в осесимметричных элементах

Будут рассчитаны напряжения в стержневых элементах

И Н И Ц И А Л И З А Ц И Я

Н о м е р в а р и а н т а 1

З а п и с ь п е р е м е щ е н и й

З а п и с ь о п о р н ы х р е а к ц и й

Проба равновесия нагружающих и опорных сил

0.000E+000 0.000E+000 0.000E+000 0.000E+000 0.000E+000 -4.300E+003

-3.802E-010 9.418E-010 1.935E-010 -4.399E-008 2.408E-008 4.300E+003

------------------------------

-3.802E-010 9.418E-010 1.935E-010 -4.399E-008 2.408E-008 -1.295E-007

Проба равновесия положительных и отрицательных нагрузок

4.689E+002 3.140E+002 7.933E+002 9.242E+004 9.300E+004 1.953E+004

-4.689E+002 -3.140E+002 -7.933E+002 -9.242E+004 -9.300E+004 -1.953E+004

------------------------------

-3.802E-010 9.418E-010 1.935E-010 -4.399E-008 2.408E-008 -1.295E-007

Р а с ч е т н а п р я ж е н и й

Всего обработано 500 элементов

О б р а б о т к а э л е м е н т о в

п о л ь з о в а т е л я

О С Р Е Д Н Е Н И Е Н А П Р Я Ж Е Н И Й

Время работы 0 час 0 мин 1 сек

Б Е З О Ш И Б О К

Реальное время проведенного расчета 0 час 0 мин 55 сек

Скрытый текст: Оптимизация разреженности

В В О Д Д А Н Н Ы Х

И Н И Ц И А Л И З А Ц И Я

Читаем таблицу ПАРАМЕТРЫ

Читаем таблицу УСКОРЕНИЕ

Читаем таблицу ДИАГРАММА_ЕТ

Читаем таблицу УЗЛЫ

Читаем таблицу ХАРАКТЕРИСТИКИ_СЕЧ

Читаем таблицу ГЕОМЕТРИЯ_СЕЧ

Читаем таблицу МАТЕРИАЛЫ

Читаем таблицу ЭЛЕМЕНТЫ

Читаем таблицу ЗАКРЕПЛЕНИЕ

Читаем таблицу УРАВНЕНИЯ

Читаем таблицу УРАВНЕНИЯ_СТЕРЖНИ

Читаем таблицу ТОЧЕЧ_ВЕС

Читаем таблицу СИЛЫ

Читаем таблицу СИЛЫ

Читаем таблицу ДАВЛЕНИЕ

Читаем таблицу ТЕМПЕРАТУРЫ

Читаем таблицу ТЕМПЕРАТУРЫ

Читаем таблицу ПОТОК_ОБМЕН

Читаем таблицу ПОТОК_ОБМЕН

Время работы 0 час 0 мин 0 сек

Б Е З О Ш И Б О К

Д И А Г Н О С Т И К А М О Д Е Л И

Оптимизация разреженности

И Н И Ц И А Л И З А Ц И Я

Р а с ч е т г е о м е т р и ч е с к и х

х а р а к т е р и с т и к с е ч е н и й

Д И А Г Н О С Т И К А М О Д Е Л И

Диагностика сечений

Диагностика материала

Диагностика ортотропного материала

Диагностика элементов

Диагностика сеток контакта

Диагностика закрепления

Диагностика уравнений

Диагностика нагружения

Диагностика нагружения края давлением

Д И А Г Н О С Т И К А Б Е З О Ш И Б О К

С о с т а в л е н и е т а б л и ц ы

Всего обработано элементов 500

О П Т И М И З А Ц И Я

С о с т а в л е н и е т а б л и ц ы

с т е п е н е й с в о б о д ы

Всего обработано элементов 500

Количество внешних степеней свободы 615

Количество степеней свободы 8490

Время работы 0 час 0 мин 1 сек

Б Е З О Ш И Б О К

О Б Р А Б О Т К А Н А Г Р У З К И

И Н И Ц И А Л И З А Ц И Я

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

м а с с

О б р а б о т к а с о б с т в е н н о г о

в е с а

Вес модели 1.681452E+003

Вес модели с учетом

точечного веса 1.681452E+003

К о о р д и н а т ы ц е н т р а

т я ж е с т и

X = 12.205887 Y = 24.499994 Z = 100.000023

Время работы 0 час 0 мин 0 сек

Б Е З О Ш И Б О К

С О С Т А В Л Е Н И Е М А Т Р И Ц Ы Ж Е С Т К О С Т И

И Н И Ц И А Л И З А Ц И Я

С о с т а в л е н и е м а т р и ц

ж е с т к о с т и э л е м е н т о в

Всего обработано 500 элементов

С о с т а в л е н и е м а т р и ц

з а д а н н ы х п о л ь з о в а т е л е м

С о с т а в л е н и е г л о б а л ь н о й

м а т р и ц ы ж е с т к о с т и

О б р а б о т а н о 8490 с т р о к

С о с т а в л е н и е о п о р н о й

м а т р и ц ы

О б р а б о т а н о 615 с т р о к

О б р а б о т к а с и л д л я п р о б ы

р а в н о в е с и я

Время работы 0 час 0 мин 4 сек

Реальное время работы 0 час 0 мин 8 сек

Б Е З О Ш И Б О К

Р Е Ш Е Н И Е У Р А В Н Е Н И Й

Матрица жесткости не будет записана на диск

Допуск диагонального члена 1.000000E-012

И Н И Ц И А Л И З А Ц И Я

Решение проводится по разреженному алгоритму

максимальный диагональный член 2.803872E+008

минимальный диагональный член 6.088548E+006

средний диагональный член 6.242097E+007

П р я м о й х о д

Всего обработано 8490 строк

О б р а б о т к а п р а в о й ч а с т и

О б р а т н ы й х о д

О п р е д е л е н и е о п о р н ы х

р е а к ц и й

Время работы 0 час 0 мин 17 сек

Реальное время работы 0 час 0 мин 19 сек

Б Е З О Ш И Б О К

В Ы В О Д Р Е З У Л Ь Т А Т О В

Будут записаны данные для постпроцессора

Будут рассчитаны напряжения в плоских элементах

Будут рассчитаны напряжения в объемных элементах

Будут записаны напряжения в осесимметричных элементах

Будут рассчитаны напряжения в стержневых элементах

И Н И Ц И А Л И З А Ц И Я

Н о м е р в а р и а н т а 1

З а п и с ь п е р е м е щ е н и й

З а п и с ь о п о р н ы х р е а к ц и й

Проба равновесия нагружающих и опорных сил

0.000E+000 0.000E+000 0.000E+000 0.000E+000 0.000E+000 -4.300E+003

1.108E-009 4.311E-010 1.215E-010 -2.037E-010 2.013E-007 4.300E+003

------------------------------

1.108E-009 4.311E-010 1.215E-010 -2.037E-010 2.013E-007 2.890E-007

Проба равновесия положительных и отрицательных нагрузок

4.689E+002 3.140E+002 7.933E+002 9.242E+004 9.300E+004 1.953E+004

-4.689E+002 -3.140E+002 -7.933E+002 -9.242E+004 -9.300E+004 -1.953E+004

------------------------------

1.108E-009 4.311E-010 1.215E-010 -2.037E-010 2.013E-007 2.890E-007

Р а с ч е т н а п р я ж е н и й

Всего обработано 500 элементов

О б р а б о т к а э л е м е н т о в

п о л ь з о в а т е л я

О С Р Е Д Н Е Н И Е Н А П Р Я Ж Е Н И Й

Время работы 0 час 0 мин 1 сек

Б Е З О Ш И Б О К

Реальное время проведенного расчета 0 час 0 мин 40 сек

Скрытый текст: Опять прямой

В В О Д Д А Н Н Ы Х

И Н И Ц И А Л И З А Ц И Я

Читаем таблицу ПАРАМЕТРЫ

Читаем таблицу УСКОРЕНИЕ

Читаем таблицу ДИАГРАММА_ЕТ

Читаем таблицу УЗЛЫ

Читаем таблицу ХАРАКТЕРИСТИКИ_СЕЧ

Читаем таблицу ГЕОМЕТРИЯ_СЕЧ

Читаем таблицу МАТЕРИАЛЫ

Читаем таблицу ЭЛЕМЕНТЫ

Читаем таблицу ЗАКРЕПЛЕНИЕ

Читаем таблицу УРАВНЕНИЯ

Читаем таблицу УРАВНЕНИЯ_СТЕРЖНИ

Читаем таблицу ТОЧЕЧ_ВЕС

Читаем таблицу СИЛЫ

Читаем таблицу СИЛЫ

Читаем таблицу ДАВЛЕНИЕ

Читаем таблицу ТЕМПЕРАТУРЫ

Читаем таблицу ТЕМПЕРАТУРЫ

Читаем таблицу ПОТОК_ОБМЕН

Читаем таблицу ПОТОК_ОБМЕН

Время работы 0 час 0 мин 0 сек

Б Е З О Ш И Б О К

Д И А Г Н О С Т И К А М О Д Е Л И

И Н И Ц И А Л И З А Ц И Я

Р а с ч е т г е о м е т р и ч е с к и х

х а р а к т е р и с т и к с е ч е н и й

Д И А Г Н О С Т И К А М О Д Е Л И

Диагностика сечений

Диагностика материала

Диагностика ортотропного материала

Диагностика элементов

Диагностика сеток контакта

Диагностика закрепления

Диагностика уравнений

Диагностика нагружения

Диагностика нагружения края давлением

Д И А Г Н О С Т И К А Б Е З О Ш И Б О К

С о с т а в л е н и е т а б л и ц ы

Всего обработано элементов 500

О п т и м и з а ц и я ш и р и н ы л е н т ы

обработано узлов 1000

обработано узлов 2000

обработано узлов 3000

обработано узлов 3035

Начальный узел оптимизации 4329

С о с т а в л е н и е т а б л и ц ы

с т е п е н е й с в о б о д ы

Всего обработано элементов 500

Количество внешних степеней свободы 615

Количество степеней свободы 8490

Время работы 0 час 0 мин 3 сек

Б Е З О Ш И Б О К

О Б Р А Б О Т К А Н А Г Р У З К И

И Н И Ц И А Л И З А Ц И Я

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

м а с с

О б р а б о т к а с о б с т в е н н о г о

в е с а

Вес модели 1.681452E+003

Вес модели с учетом

точечного веса 1.681452E+003

К о о р д и н а т ы ц е н т р а

т я ж е с т и

X = 12.205887 Y = 24.499994 Z = 100.000023

Время работы 0 час 0 мин 0 сек

Б Е З О Ш И Б О К

С О С Т А В Л Е Н И Е М А Т Р И Ц Ы Ж Е С Т К О С Т И

И Н И Ц И А Л И З А Ц И Я

С о с т а в л е н и е м а т р и ц

ж е с т к о с т и э л е м е н т о в

Всего обработано 500 элементов

С о с т а в л е н и е м а т р и ц

з а д а н н ы х п о л ь з о в а т е л е м

С о с т а в л е н и е г л о б а л ь н о й

м а т р и ц ы ж е с т к о с т и

О б р а б о т а н о 8490 с т р о к

С о с т а в л е н и е о п о р н о й

м а т р и ц ы

О б р а б о т а н о 615 с т р о к

О б р а б о т к а с и л д л я п р о б ы

р а в н о в е с и я

Время работы 0 час 0 мин 4 сек

Реальное время работы 0 час 0 мин 8 сек

Б Е З О Ш И Б О К

Р Е Ш Е Н И Е У Р А В Н Е Н И Й

Матрица жесткости не будет записана на диск

Допуск диагонального члена 1.000000E-012

И Н И Ц И А Л И З А Ц И Я

П е р е п и с ы в а е м м а т р и ц у

ж е с т к о с т и

Прогнозируемое время решения

для Intel I7 2900 МГц 0 час 0 мин 1 сек

переписано 5000 степеней

Всего переписано 8490 степеней

Максимальная ширина ленты 897 в строке 4411

Средняя ширина ленты 601

максимальный диагональный член 2.803872E+008

минимальный диагональный член 6.088548E+006

средний диагональный член 6.242097E+007

Р а з л о ж е н и е м а т р и ц ы

обработано 613 строк

обработано 1236 строк

обработано 2214 строк

обработано 2999 строк

обработано 3504 строк

обработано 4009 строк

обработано 4572 строк

обработано 5113 строк

обработано 5734 строк

обработано 6295 строк

обработано 6858 строк

обработано 7536 строк

обработано 8490 строк

О б р а б о т к а п р а в о й ч а с т и

О б р а т н ы й х о д

О п р е д е л е н и е о п о р н ы х

р е а к ц и й

Время работы 0 час 0 мин 23 сек

Реальное время работы 0 час 0 мин 27 сек

Б Е З О Ш И Б О К

В Ы В О Д Р Е З У Л Ь Т А Т О В

Будут записаны данные для постпроцессора

Будут рассчитаны напряжения в плоских элементах

Будут рассчитаны напряжения в объемных элементах

Будут записаны напряжения в осесимметричных элементах

Будут рассчитаны напряжения в стержневых элементах

И Н И Ц И А Л И З А Ц И Я

Н о м е р в а р и а н т а 1

З а п и с ь п е р е м е щ е н и й

З а п и с ь о п о р н ы х р е а к ц и й

Проба равновесия нагружающих и опорных сил

0.000E+000 0.000E+000 0.000E+000 0.000E+000 0.000E+000 -4.300E+003

-3.802E-010 9.418E-010 1.935E-010 -4.399E-008 2.408E-008 4.300E+003

------------------------------

-3.802E-010 9.418E-010 1.935E-010 -4.399E-008 2.408E-008 -1.295E-007

Проба равновесия положительных и отрицательных нагрузок

4.689E+002 3.140E+002 7.933E+002 9.242E+004 9.300E+004 1.953E+004

-4.689E+002 -3.140E+002 -7.933E+002 -9.242E+004 -9.300E+004 -1.953E+004

------------------------------

-3.802E-010 9.418E-010 1.935E-010 -4.399E-008 2.408E-008 -1.295E-007

Р а с ч е т н а п р я ж е н и й

Всего обработано 500 элементов

О б р а б о т к а э л е м е н т о в

п о л ь з о в а т е л я

О С Р Е Д Н Е Н И Е Н А П Р Я Ж Е Н И Й

Время работы 0 час 0 мин 1 сек

Б Е З О Ш И Б О К

Реальное время проведенного расчета 0 час 0 мин 51 сек

Модель: ANM.rar

[Борман 2 396]

Это я уже понял... кнопки я нажимаю хорошо..

[Борман 2 396]

У меня какие то идеологические непонятки... О чем этот шаг ?

Постпроцессор - Вывод результатов.

Это рестарт решения уравнений в удобной для пользователя форме. Еще одна фишка.

Ну как бы не особо фишка.. Жаль что итерации опять с начала отситывает, непонятно продолжает ли вообще.

[Борман 2 396]

Ну ладно, ИСПА, устраивайте экзамен. Надо же мне выпуститься из младшей школы.

[Борман 2 396]

Запустите GPROF. Загрузите эту модель.

Установите деформированное состояние

Постпроцессор - Деформации - деформированное состояние

Затем

Постпроцессор - Анимация - Создание анимации (модификация)

И управляйте положением модели мышкой.

Ну да.. этого то мне как раз и не хватало...

Можно ли вывести разом реакции во всех опорах, а не отпрашивать узлы по одному? Если я хочу узнать сумму реакций по некоторым узлам ?

Можно ли строить эпюры по твердотельным элементам ?

[Борман 2 396]

А что вы понимаете под твердртельным элементом? Вы создали МЭ. По ним получили КЭ. Откуды появились твердые тела?

Ну вы и буквоед.. Просто эпюру сквозь партию 28-х элементов.

[Борман 2 396]

Поэтому вопрос буквоеда, что в данном случае ЭПЮРА?

Есть некая воображаемая линия, которая проходит сквозь тело (из 28-х элементов). Нужно нарисовать/вывести численно и т.п. значение того или иного поля по этой линии с заданным шагом по длине.

[Борман 2 396]

Нам бы схемку, иль чертеж. Мы б затеяли вертеж.

Да все просто... Есть у вас тело, одна грань у него закреплена, например. Само тело нагружено. Необходимо узнать усилия (три силы, три момента) относительно какой-либо точки по закрепленной грани. Для сложности считаем, что крепиться не грань, а узлы на ней..

[Борман 2 396]

А если в теле много граней закреплено?

В общем случае - по наперед заданным граням, или узлам.

[Борман 2 396]

Скажу, как это делается в Ансисе..

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

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

Видео карта считает 2 сек.

Я думаю, что Ансис тоже не будет с 8000 степенями свободы долго париться..

[Борман 2 396]

Борман вы проверьте копирование 1 МЭ.

Работает.

[Борман 2 396]

Причина в предобуславливателе.

А прямые методы не дружат с GPU ?

[a_schelyaev 370]

Добавлю, что по этой же причине и Абакус не может перейти на GPU.

В Абакусе все силы брошены на ускорение факторизации силами GPU. И ускорение уже есть.

[Fedor 1 625]

" Там реализован блочно-разреженный предобуславливатель" - наилучший из возможных предобуславливателей - обратная матрица

Я вообще-то видал статьи по алгебре. Эта уж очень не похожа, поэтому так и считаю. Коль занимаетесь математикой, то будте любезны придерживаться принципов которые заложили многие поколения классиков. Основное - цель понимание того, что происходит и как. А я этого в этой рекламе не увидел. Но это мнение постмодерниста, а не модерниста

[Борман 2 396]

ИСПА, а не сравнить ли нам "ИСПА" и АНСИС по скорости на ЦПУ моей тачки?

Вы скажете самые оптимальные настройки решателя "ИСПА", я оставлю все по-умолчанию в АНСИС...

[a_schelyaev 370]

Так я и говорю о факторизации матрицы.

Не будет ускорения на GPU. На ЦПУ быстрее будет и дешевле железо.

Я то говорю, что ускорение есть.

Вот логи повнимательнее изучу и будет отчетик по Абакусу.

[a_schelyaev 370]

ИНТЕЛ тоже не стоит на месте. Сравнивать скорость уже нужно с 12 ядрами от ИНТЕЛ.

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

А на Абакусе, пока что, GPU позиционируется строго под большие задачи.

[a_schelyaev 370]

От перемены мест слагаемых сумма не меняется.

Ускорение есть.

:)