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

Сравнение автоматических 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 пользователей

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


  • Сообщения

    • Bot
      Информационное моделирование: кто «платит и заказывает музыку»?

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

      Об авторах: Мария Золотова — BIM-специалист строительной компании, по первому образованию — театральный режиссер, Софья Минина — BIM-менеджер, по второму образованию — психолог, Владимир Талапов — BIM-эксперт, всего одно образование — математик. См. также BIM и Театр: искусство возможного. Часть 1, BIM и Театр: искусство возможного. Часть 2, BIM и Театр: искусство возможного. Часть 3. Информационное моделирование зданий и сооружений, как творческий процесс, имеет много общего с не менее творческим процессом постановки театральных спектаклей. Поэтому изучение театрального опыта, его наработанных взаимосвязей участников и методик их взаимодействия может существенно помочь становлению и использованию BIM в строительной отрасли. В этой части речь пойдет о тех, кто вроде бы непосредственно не участвует в постановке спектакля или информационном моделировании, но от кого эти процессы существенно зависят. «Кто платит, тот и заказывает музыку» Итак, кто же [...] View the full article
    • Guhl
      g71

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

      Понятия не имею
    • Александр 36
      g71

      Автор: Александр 36 · Опубликовано:

      А где можно скачать русифицированную версию
    • BASH_HD
      Подключение 5 оси(индексной) к 31i-A

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

      В Москве) но вообще по мнению японцев станок в Германии) но это не точно...   Аукцон и BAFA знали куда продали.
    • Alexandr97
      Проблема с обновлением созданной библиотеке в tflex docs.

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

      Спасибо!
    • Andrey_kzn
      Sinumerik 840 DSL Ошибка 21612

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

      Логика работы у них та же самая. И интерфейсные сигналы те же. Вы упускаете из виду сигнал готовность привода. Это DB3x.DBX93.5. Выстраивайте логическую цепь событий. Формирование сигнала DDX2.1 может быть завязано на готовность привода.  Не будет готовности привода, можете сколько угодно выставлять DBX2.1, ось никуда не поедет. А готовность может пропадать из-за дефекта энкодера, кабеля, разъёмов и т.п.  Надо сделать трассировку всех этих приводных сигналов,  заодно посмотрите также DBX93.7 и DBX21.7.  И понять, что какой сигнал пропадает первым.
    • Guhl
      Fanuc CNC GUIDA 25.0 x86

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

      Опять я забыл про это меню. Оно в корне лежит, а я туплю Спасибо  
    • gudstartup
      Fanuc CNC GUIDA 25.0 x86

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

      @Guhl все это настраивается здесь
    • Guhl
      Fanuc CNC GUIDA 25.0 x86

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

      У меня 25я версия, но может ревизия не та. Скачаю с рутрэкера, посмотрю 
    • mrVladimir
      Fanuc CNC GUIDA 25.0 x86

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

      Странно. Сейчас попробовал выборочно разные ЧПУ - везде есть такие пункты. А у Вас на других также нет этих панелей? Может что-то при установке не так пошло...  
×
×
  • Создать...