Перейти к публикации

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

Борман

Автор: Борман,
в Построение расчетных сеток

 Рассказать Подписчики 0

Рекомендованные сообщения

Борман

Борман 2 375

Опубликовано:
Опубликовано:

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

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

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

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

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

*GET Preprocessing Items, Entity = ELEM

Item1=SHPAR

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

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

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

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

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


sergeyd

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 375

Опубликовано:
  • Автор
Опубликовано:

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

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

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

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

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

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

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

Борман 2 375

Опубликовано:
  • Автор
Опубликовано:

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

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

Борман 2 375

Опубликовано:
  • Автор
Опубликовано:

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

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

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

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

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

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

Борман 2 375

Опубликовано:
  • Автор
Опубликовано:

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

Да можно и на этой.. Готов строить сетки на этой области с заданным размеров элемента (ребра).
Ссылка на сообщение
Поделиться на других сайтах
Борман

Борман 2 375

Опубликовано:
  • Автор
Опубликовано:

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

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

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

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

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

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

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

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

Борман 2 375

Опубликовано:
  • Автор
Опубликовано:

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

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

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

post-5875-1278536940_thumb.jpg

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

Борман 2 375

Опубликовано:
  • Автор
Опубликовано:

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

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

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

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

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

Борман 2 375

Опубликовано:
  • Автор
Опубликовано:

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

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

Борман 2 375

Опубликовано:
  • Автор
Опубликовано:

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

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

sergeyd 3

Опубликовано:
Опубликовано:

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

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

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

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

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

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

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

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

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 375

Опубликовано:
  • Автор
Опубликовано:

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

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

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

Islander 2

Опубликовано:
Опубликовано:

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

Ansys v12.0.1

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

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

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

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

Борман 2 375

Опубликовано:
  • Автор
Опубликовано:

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

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

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

Присоединяйтесь к обсуждению

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

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

×
 Рассказать
Подписчики 0
Перейти к списку тем

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

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


  • Сообщения

    • gudstartup
      Вероятно слетел драйвер ключа DMU-50

      Автор: gudstartup · Опубликовано:

      считывание происходит при помощи вх\вых сигналов контроллера plc 
    • gudstartup
      Подпрограмма пересчёта положения системы координат через sin, cos для fanuc. Кинематика 3+2, поворотные оси B,C.

      Автор: gudstartup · Опубликовано:

      @Maks Horhe так все таки скиньте бэкап эмулируем ваше чпу в cncguide и посмотрим куда поедет?  можете снять видео с фиксацией координатных позиций после каждого кадра. Выложу вашу программу пусть программисты посмотрят все ли в ней ок. %O0002 G40 G17 G94 G90 G49 G80 N1 G91 G28 Z0.0 N2 G91 G28 X0.0 Y0.0 N3 G91 G28 B0.0 C0.0 N4 M03 S200 N5 G90 G0 G53 B0.0 C0.0 N6 G54 N7 X0.0 Y0.0 N8 G90 G43 H01 N9 G90 G0 X0.0 Y0.0 N10 G90 G0 Z200.0 N11 G01 Z10.0 F1500. N12 M00 N13 G00 Z200.0 N14 G40 G49 G69 N15 G00 G53 Z0.0 N16 G00 G54 B0.0 C0.0 N17 G68.2 X0.0 Y0.0 Z0.0 1135. J39.2044 K-129.2315 N18 G53.1 N19 G01 X0.0 Y0.0 F1500 N20 G90 G43 H01 N21 G90 G01 X0.0 Y0.0 F1500 N22 G90 G01 Z200. F1500 N24 M00 N25 G00 Z200. N26 G40 G49 G69 N27 G91 G28 Z0.0 N28 G28 X0.0 Y0.0 N29 G91 G28 B0.0 C0.0 N30 M5 N31 M30
    • gudstartup
      Вероятно слетел драйвер ключа DMU-50

      Автор: gudstartup · Опубликовано:

      @karlf 530 считывает ключ по специальному протоколу при помощи plc и получает его серийный номер а из него определяет возможные режимы доступа. там нет драйвера а есть plc модуль или несколько эти модули написаны на питоне  надпись smartkey исчезает с экрана при запуске чпу??
    • ДОБРЯК
      Расчет на сейсмику

      Автор: ДОБРЯК · Опубликовано:

      Для того, чтобы получить правильные высшие) формы при виртуальном эксперименте, нужно сделать грамотную КЭ модель. От разговора на эту тему вы постоянно уклоняетесь.  То нет компьютера под рукой, сделать простейший тест, то теряете интерес. :=) Сходимости энергии деформации при расчетах статики, недостаточно для точного определения высших собственных форм и частот.  Для того, чтобы грамотно использовать метод конечных элементов, нужно сделать много-много тестов в статике, динамике и ... Одной кнопки и двух конечных элементов в 3Д программе недостаточно для определения высших собственных форм...  У вас в качестве инструмента всего два конечных элемента, шести узловая несовместная оболочка Тимошенко и десяти узловой тетраэдр. И еще контакты при решении задачи на собственные числа. Вам ли говорить про правильность определения высших собственных форм для сложных изделий... :=)    
    • vad0000
      Delem 65. Нужна консультация, не работает ось X.

      Автор: vad0000 · Опубликовано:

      Покажите схему с разрешением на движение
    • vad0000
      Delem 65. Нужна консультация, не работает ось X.

      Автор: vad0000 · Опубликовано:

      Вход, а не выход Вытащить Аналоговый вход и все, как будто туда ничего не подключено И если мы подключим сигнал к энкодеру оси Х, то он стнтет одинаковый с аналоговым входом, который не подключен?
    • Snake 60
      Как зафиксировать деталь в "воздухе" во время расчета в simulation?

      Автор: Snake 60 · Опубликовано:

      @waze4534  Посмотрите вверх и прочитайте текст на красной полоске...
    • kkk
      Помощь новичку Solidworks

      Автор: kkk · Опубликовано:

      Я так понимаю, что предупреждение про "касательные" не просто так выскакивает. Если скруглить прямую стыковку отрезков эскиза (минимальным радиусом) то все работает даже без объединенной кривой, достаточно эскиза.
    • karlf
      Вероятно слетел драйвер ключа DMU-50

      Автор: karlf · Опубликовано:

      Подскажите пожалуйста, может кто сталкивался. Станок DMU-50 на стойке TNC 530, перестал определяться ключ доступа. Сам ключ вроде работает, если переключать на нём режимы, то в шкафу на соответствующих блоках лампочки тоже переключаются. Но изначально был уровень доступа 4, а теперь уровни доступа не активны. Ключ только один, запасных нет. Есть какой-то старый бэкап, пробовал его накатить, но какой-то он непонятный - станок грузится, но почти в конце загрузки выдаёт какую-то ошибку по параметрам. Может кто знает в каком из разделов и в какой папке искать установленные драйверы ключа?
    • YuriySt54
      Ограничение лишнего при построении массива

      Автор: YuriySt54 · Опубликовано:

      В данной конструкции можно сразу одним телом делать и обрезать. Но меня интересовало, можно ли сразу при построении массива ограничить лишнее.
×
×
  • Создать...