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

[Борман 2 398]

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

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

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

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

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

*GET Preprocessing Items, Entity = ELEM

Item1=SHPAR

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

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

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

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

[sergeyd 3]

на 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

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

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

[Борман 2 398]

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

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

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

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

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

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

[sergeyd 3]

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

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

[sergeyd 3]

2 Испа

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

[Борман 2 398]

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

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

[Борман 2 398]

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

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

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

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

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

[Борман 2 398]

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

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

[Борман 2 398]

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

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

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

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

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

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

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

[Борман 2 398]

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

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

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

[Борман 2 398]

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

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

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

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

[Борман 2 398]

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

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

[Борман 2 398]

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

Взял куб 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

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

[sergeyd 3]

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

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

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

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

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

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

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

[Islander 2]

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

[Борман 2 398]

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

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

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

[Islander 2]

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

Ansys v12.0.1

[Vova 419]

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

Препрцессор 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

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

[Islander 2]

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

[Борман 2 398]

Спасибо 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одержит плоские элементы по поверхности.

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

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