Борман

Сравнение автоматических TETRA-генераторов.

409 сообщений в этой теме

В теме будет произведен сравнительный анализ различных автоматических тетра-негераторов (HyperMesh, ANSYS, Nastran, ABAQUS, ИСПА, Cosmos, ... ). Кроме того приглашаю коллег CFD-шников присоединиться к обсуждению.

Господа, ну вы все поняли.

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

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

Обращаюсь тем, кто работает с ANSYS наиболее плотно. Какие показатели качества можно снять с тетра-элемента? Наспех нашел только это...

*GET Preprocessing Items, Entity = ELEM

Item1=SHPAR

У кого набита рука, напишите пожалуйста код для расчета мат. ожидания и дисперсии этих показателей (и тех, что еще найдете). Кто общается с APDL ежедневно, тому делов на 5-10 минут.

В общем жду инициативы.

PS. Рассматривам любые пакеты.. Жалко, конечно, что ИСПА не делает экспорт никуда, кроме как в формат ИСПА.

Экспорт только в формате ИСПА.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

на apdl написать быстрее чем 5 минут. и готов это сделать, НО

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

Форма не все. одинаковые элем в разных программах дадут разное.

а разные могут дать примерно одинаковое.

в ансис элементы очень вытянутые дают приемлемые результаты.

за счет улучшенных функций формы.

*get,ke,elem,,count

ne=0

*dim,shp,,ke,1

*do,j,1,ke

ne=elnext(ne)

! здесь подставить запрос нужного

*get,shp(j),elem,ne,....

*enddo

*vscfun,mea,mean,shp(1)

*vscfun,std,stdv,shp(1)

*msg,ui,mea,std

mean value=%13.5e%/standdeviation=%13.5e

друзья=вам явно скучно...

могу предложить намного более полезные задачи...

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Форма не все. одинаковые элем в разных программах дадут разное.

а разные могут дать примерно одинаковое.

в ансис элементы очень вытянутые дают приемлемые результаты.

за счет улучшенных функций формы.

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

В общем неконструктивно. Ждем заинтересованных...

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

цель - не красивая сетка, а верные результаты. и быстро.

удачи.Зенкевич в помощь.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

2 Испа

с вами даже неинтересно дискутировать.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

как можно крупную тетра сетку

ИСПА, это не формализуемый критерий. У генераторов есть вполне четкие параметры сетки. Так же огласите правила игры... иначе у вас преимущество.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Если у генератора есть четкие параметры, то кто-то их формализует

Их формализует пользовательский интерфейс.

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

Я не могу строить сетки с таким условием. Галочки (максимизировать размер элемента) такой нет, видите ли..

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

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Что же будем тренироваться на более простых областях. Давайте пример. Язык показывать не буду, не имею такой привычки.

Да можно и на этой.. Готов строить сетки на этой области с заданным размеров элемента (ребра).

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Пустая болтовня пошла...

1) Фиксируем геом. модель.

2) Фиксируем параметры сетки

3) Договаривается что и как сравниваем.

4) Нажимаем нужную кнопку.

5) Сравниваем.

6) Убеждаемся что...

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Написать переходник из формата ИСПА в любую другую систему например АНСИС это дело одного дня.

По-упражнялся. :sport_box:

Кстати, блин, ваш STEP отвратительно импортировался в ГМ. У него всегда был проблемы с импортом таких (тонкостенных) вещей. Чувствую, перейду все таки на более молодую версию ГМ.. там все это полечили.

post-5875-1278536940_thumb.jpg

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Какие мои действия по гекса сеткам? Руководите!!!

Вроде тетра обуждаем.

Предстваьте, что вы ничего не знаете в ИСПА, и вам нужно построить тетра-сеть. Перед вами три объекта - поля для ввода MIN EDGE, MAX EDGE, и кнопка "ПУСК". Делайте на кубе 10х10х10 сеть тетров с ребром min=max=1. И кидайте сеть в своем формате. Я сделаю то же самое в АНСИС и ГМ.

Вечером я чего-нибудь насчитаю. Все равно начинать надо с простого..

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Буду тренироваться на кубиках (на кошках)

До Федора вам еще далеко. Федор может из интеграла Лебега пельменей налепить и сварить в кастрюле в неевклидовой метрикой. А вы пока тренеруйтесь на кубиках.. :drinks_drunk:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Ладно, буду играть в открытую.

Взял куб 10х10х10 и покрывал тетра-сеткой размером 1.0. Настроек вообще никаких.. Специально для ИСПА код ANSYS:

ET,1,45 - элемент

BLOCK,0,10,0,10,0,10, - объем

ESIZE,1,0, - размер элемента =1.0, 0 - неиспользуемый параметр.

MSHAPE,1,3D - 3D, triangle-shaped сетка.

MSHKEY,0 - Use free meshing

VMESH,1

В HM тоже единственная настройка - размер = 1.0 без всяких галочек.

В слову, для "неопытных" пользователей АНСИС выдает такое сообщение:

<noindex>Изображение</noindex>

ну т.е. "WARNING !!! Покупай меня не глядя !!!"

Простейшие параметры получившихся сеток:

Пока смотрел только объем (мат. ожидание, дисперсия, узлов, элементов)

ANSYS 10

MEAN =

STDEV =

Number of Defined Nodes =

Number of Defined Elements =

HM

MEAN =

STDEV =

Number of Defined Nodes =

Number of Defined Elements =

К слову.. Объем правильного тетраедра со стороной 1.0 равен 0.1178511.

Очень интересно получилось, однако. Средний объем отличается ровно в 2 раза. Вроде правильно все задавал... MEAN*NE - примерно одинаковые..

Перепроверил - все верно. Что скажете, ИСПА ?

sergeyd

Спасибо, скрипт работает.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

кстати, для 45 это избыточная диагностика=перестраховка.

ничего особенного не замечал.

хотя предпочитаю 186 и 187.

еще попробуйте 285 тип= тетраэдр 4-узловой.

(свои специальные функции формы, а не вырожденные из гексаэдра)

тестировать особо времени не было= но напряжения на вид "глаже".

и в нелине он "робастнее"

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Number of Selected Nodes = 1635

Maximum Element Number = 7660

См. сюда (Ansys, все по умолчанию):

/prep7

ET,1,SOLID45

!*

BLOCK,10,,10,,10,,

ESIZE,1,0,

MSHAPE,1

MSHKEY,0

VMESH,all

!

! Number

! Nodes . . . . . . . . . . . 1597

! Elements. . . . . . . . . . 7431

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

! Nodes . . . . . . . . . . . 1597

! Elements. . . . . . . . . . 7431

Копирую ваш код в только что открытый АНСИС - и получаю свои числа. У вас какая версия? У меня v10.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Копирую ваш код в только что открытый АНСИС - и получаю свои числа. У вас какая версия? У меня v10.

Ansys v12.0.1

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Кубик как у Бормана.

Препрцессор Medina:

ничего не менял, всё автоматом генерировалось.

Number of TET4 elements = 6094

Number of tetra nodes = 1336

Shortest tetra edge = 0.3104

Longest tetra edge = 2.1036

Сетка в Настране

Препрцессор Ansa:

ничего не менял, всё автоматом генерировалось.

Number of TET4 elements = 6129

Number of tetra nodes = 1353

Shortest tetra edge = 0.466

Longest tetra edge = 1.857

Сетка в Настране

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Чудны дела твои, о Господи...

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Спасибо Vova

MEDINA

MEAN =

STDEV =

Number of Defined Nodes =

Number of Defined Elements =

ANSA

MEAN =

STDEV =

Number of Defined Nodes =

Number of Defined Elements =

PS. Сетка ANSA cодержит плоские элементы по поверхности.

Чудны дела твои, о Господи...

Сделайте экспорт сети.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

[/color]PS. Сетка ANSA cодержит плоские элементы по поверхности.

Поправил..

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Vova, подскажи пожалуйста как здесь файлы красиво прикреплять. Что-то не нашел я на форуме файловый архив.

***

Кинул зипованый файл на <noindex>http://ifolder.ru/18472089</noindex>

В архиве находится cdb в нейтральном формате Ансиса (и геометрия и сетка). Закачать в Ансис можно как через Реад Инпут, так и через Реад Архив.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Islander

Он как VIP наверное... Мне тоже понравилось..

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Vova, подскажи пожалуйста как здесь файлы красиво прикреплять. Что-то не нашел я на форуме файловый архив.

Нажимаешь ответить.

Внизу стоит

Выберите файл - нажимаеш серую кнопку у меня там стоит что типа "искат" выбираешь файл и жмешь кнопку загрузить.

в маленьком меню можно ещё файл в текст вставлять.

всё

Islander

Он как VIP наверное... Мне тоже понравилось..

Может быть, я то думал зачем мне этот статус нужен. Полезная вещь однако.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Спасибо, Vova! Вроде получилось - архивы цепляет (rar и zip):)

as.zip

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Islander тоже подключился...

ANSYS 12.1

MEAN =

STDEV =

Number of Defined Nodes =

Number of Defined Elements =

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Вот теперь ИСПА дождемся, и узнаем, у кого он списал свой генератор. :lamo:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

тогда надо ещё инфу по Патрану, Femapу и ICEM. Тогда может и определим. :smile:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Воркбенч.

По умолчанию бьет такую модель на кирпичи и на КЭ со срединным узлом. Выставил метод - "на тетраэдры", размер сторон элемента - 1 мм и опцию "без срединных узлов". В итоге получил сетку:

узлов - 1846, КЭ - 8907.

В прикрепленном файле - командный файл Ансиса для построения сетки (в классике).

post-1320-1278666603_thumb.jpg

aw12_0.zip

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

По умолчанию он бьет на кубы, добавляет диагонали под 45 гр.

А то что видишь ты это сглаживание идет вледом за построением сетки.

Если сглаживание отключить, то увидишь что он делает.

Для чистоты эксперимента возьми ICEM и там поупражняйся.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

По умолчанию он бьет на кубы, добавляет диагонали под 45 гр.

А то что видишь ты это сглаживание идет вледом за построением сетки.

Если сглаживание отключить, то увидишь что он делает.

Для чистоты эксперимента возьми ICEM и там поупражняйся.

А вот нужно-ли сглаживание отключать, которое по умолчанию?

Про сглаживание - не новость, поэтому и привел спецом картинку с развернутой вкладкой Details of "Mesh".

Что мне в Ансисе нравиться, так это толково расставленные "по умолчанию" опции.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

В ICEM все задается как хочешь. При этом есть режимы генерации, где защита на дурака введена очень агрессивная.

Ну и также нужно помнить, что процедуры сглаживания сетки там стоят по умолчанию.

А вот нужно-ли сглаживание отключать, которое по умолчанию?

От задачи зависит. В данном случае сглаживание не нужно.

Товарищи же хотят как можно крупнее сетку.

Что мне в Ансисе нравиться, так это толково расставленные "по умолчанию" опции.

Насколько я понял, во всех FEA пакетах такие параметры заданы.

post-1864-1278668493_thumb.jpg

post-1864-1278668501_thumb.jpg

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Насколько я понял, во всех FEA пакетах такие параметры заданы.

Есть, но по разному заданы. Об этом надо помнить, когда речь идет о сравнении.

Да и при сравнении решений нужно про это помнить, так например для модального анализа в Ансис и Настране выставлены разные опции по умолчанию (Настран - LUMPM - ON (lumped mass matrix), Ansis - LUMPM - OFF).

От задачи зависит. В данном случае сглаживание не нужно.

Товарищи же хотят как можно крупнее сетку.

Это было на первых порах, но потом вроде ограничились сравнением КЭ сеток, получаемых в пакетах на автомате.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

На автомате не получается, т.к. придется все равно какие-то параметры же начальные вводить. Например ограничения по максимальному размеру.

Ну раз такое дело, то я выставил максимальный размер элемента, превышающий сторону фигуры с запасом.

Статистика:

Element types :

NODE : 16

LINE_2 : 456

TETRA_4 : 19966

TRI_3 : 7666

И по кубику

Element types :

NODE : 8

LINE_2 : 127

TETRA_4 : 15655

TRI_3 : 1396

Кстати, а если применить Флюентовскую сеточную решалку, то результат следующий:

Element types :

NODE : 8

LINE_2 : 120

TETRA_4 : 6627

TRI_3 : 1368

post-1864-1278671209_thumb.jpg

post-1864-1278671617_thumb.jpg

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Воркбенч.

WB12

MEAN =

STDEV =

Number of Defined Nodes =

Number of Defined Elements =

Средний объем элемента меньше объема правильного тетраедра со стороной 1.0. Думаю это приговор для WB. Как все сетки появятся - сделаю summary.

Господа, прошу никакого шаманства с сеткой. Настроек две: Размер=1.0, сеть - ТЕТРА. Ожидаем сетки ИСПА, Femap, Патран, далее опционально.

Прошу выкладывать остальные сетки. В выходные буду недоступен...

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Хоть убей не пойму, чего привязались к правильному тетраэдру. Почему не к тетраэдру со сторонами 1 мм (образующие прямой угол между собой. На гранях куба межузловое расстояние 1 мм.) Объемчик элемента в этом случае - 0,16667 мм3.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

ISPA

Все что "самое кравивое" - вне конкурса. Как я нонял без настроек строится сетка "Вот такой вариант кубика".

Если сегодня вечером не упаду от усталости - сделаю сводную табличку. А то что-то не привык с ночного поезда сразу на работу идти...

Если не переругаемся, перейдем поотм и к самым красивым сеткам.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

В ICEM CFD тоже 4 генератора.

:)

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

в ИСПА четыре автоматических генератора тетра сеток. И все они (четыре) имеют право принимать участие в конкурсе

Но какой-то из них стоит по-умолчанию...

Кстати не увидел сетки ГМ и ANSYS. Выложите .NAS плиз.

Шэф, усё будет сделано у лучшем виде. "Все будет хорошо, даже лучше".

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

ISPA

ну чтож вы меня подставляете? Народ мне верит, а тут вы со своими тремя сетками...

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Создайте аккаунт или войдите для комментирования

Вы должны быть пользователем, чтобы оставить комментарий

Создать аккаунт

Зарегистрируйтесь для получения аккаунта. Это просто!


Зарегистрировать аккаунт

Войти

Уже зарегистрированы? Войдите здесь.


Войти сейчас

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу

  • реклама

  • Реклама

  • Ближайшие события

    Предстоящих событий не найдено
  • Дни рождения сегодня

    1. _Anna_
      _Anna_
      (31 год)
    2. Alecksey_K
      Alecksey_K
      (32 года)
    3. Astrabawerynd
      Astrabawerynd
      (32 года)
    4. Bobikos
      Bobikos
      (42 года)
    5. enenainureTat
      enenainureTat
      (49 лет)
    Просмотреть все