Jump to content

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

Борман

By Борман,
in Построение расчетных сеток

 Share Followers 0

Recommended Posts

Борман

Борман 2,286

Posted
Posted

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

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

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

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

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

*GET Preprocessing Items, Entity = ELEM

Item1=SHPAR

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

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

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

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

Link to post
Share on other sites


sergeyd

sergeyd 3

Posted
Posted

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

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

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

Link to post
Share on other sites
Борман

Борман 2,286

Posted
  • Author
Posted

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

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

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

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

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

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

Link to post
Share on other sites
sergeyd

sergeyd 3

Posted
Posted

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

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

Link to post
Share on other sites
Борман

Борман 2,286

Posted
  • Author
Posted

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

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

Борман 2,286

Posted
  • Author
Posted

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

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

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

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

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

Link to post
Share on other sites
Борман

Борман 2,286

Posted
  • Author
Posted

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

Да можно и на этой.. Готов строить сетки на этой области с заданным размеров элемента (ребра).
Link to post
Share on other sites
Борман

Борман 2,286

Posted
  • Author
Posted

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

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

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

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

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

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

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

Link to post
Share on other sites
Борман

Борман 2,286

Posted
  • Author
Posted

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

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

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

post-5875-1278536940_thumb.jpg

Link to post
Share on other sites
Борман

Борман 2,286

Posted
  • Author
Posted

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

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

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

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

Link to post
Share on other sites
Борман

Борман 2,286

Posted
  • Author
Posted

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

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

Борман 2,286

Posted
  • Author
Posted

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

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

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

Link to post
Share on other sites
sergeyd

sergeyd 3

Posted
Posted

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

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

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

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

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

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

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

Link to post
Share on other sites
Islander

Islander 2

Posted
Posted

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

Link to post
Share on other sites
Борман

Борман 2,286

Posted
  • Author
Posted

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

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

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

Islander 2

Posted
Posted

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

Ansys v12.0.1

Link to post
Share on other sites
Vova

Vova 411

Posted
Posted

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

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

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

Link to post
Share on other sites
Борман

Борман 2,286

Posted
  • Author
Posted

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

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

Сделайте экспорт сети.
Link to post
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.
Note: Your post will require moderator approval before it will be visible.

Guest
Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

×
 Share
Followers 0
Go to topic listing

  • Recently Browsing   0 members

    No registered users viewing this page.


  • Сообщения

    • Kelny
      Как на листе детали поставить ссылку на номер позиции, под которой она записана в сборке?

      By Kelny · Posted

      Возможно вам при таких запросах надо выпускать чертежи иначе - в виде каталогов и в этом возможно поможет SOLIDWORKS Composer.   Ещё можно глянуть: ЕСКД Каталог изделия. ЕСКД Структура изделия.   А что бы всё было интерактивно выпускать не бумажную КД, а интерактивный документ, что бы клик по позиции в таблице просто приводил к чертежу выбранной позиции без всяких доп.пометок в документе чертежа детали.
    • Россиянин1975
      UnigraphicsNX help

      By Россиянин1975 · Posted

      надо обмозговать...если деталь разработана конструктором и она пришла мне в обработку как технологу, на стадии конструирования я уже сделать ничего не смогу. А вот если я делаю то вы конечно правы. Ну например я затянул пов-ть паза, а потом мне же ведь придется эту поверхность подавить для того чтоб я сделал паз (это касается не только паза)? или как-то я недопонимаю. спасибо , немного нехватило времени допереть поэкспериментировать самому. Вот так вот реально сложный вопрос будет чтоб я без таких как вы гуру делал...Тут намекают что ютуб скоро в россии закроют, без этого форума совсем трудно будет осваивать...
    • Суслик
      Какой параметр отвечает за энергосбережение и где его искать?

      By Суслик · Posted

      С трудом разобрался как копировать. Надеюсь то, что надо... https://drive.google.com/file/d/1TLpVBbRKm-L3hGqLGIYoPvdW3xUdfY8m/view?usp=drivesdk К сожалению был на выходных, сейчас в ночную работаю. Теоретически по инструкции нашел параметр, отвечающий за авто отключение. Но он это или нет, без понятия. Кроме того, в самом окне управления реле нашел и таймер на авто отключение гидросистемы (по времени примерно похоже), и само реле. Но почему-то не дает отредактировать, ссылаясь на "незавершенность записи". Но это только теоретически, так как возможности заняться этим в плотную нет. Постоянные наладки на 1000+ деталей...
    • Krusnik
      Траектория по координате с градусом, нужна помощь!

      By Krusnik · Posted

      Он пишет, что R*sin(α) - это координата Z R*cos(α) - это координата Y Где R - радиус детали координата X написана в таблице   Алгебру в школе не учили?   Берёте эксель, забиваете в него цифры из таблицы и формулы выше. Полученные значения вставляете в кривую по точкам.
    • 01aaddss
      DMU 70 TNC 426/430M

      By 01aaddss · Posted

      merhaba, dmu50 makine ve post dosyasının bir simülasyonuna ihtiyacım var, bana yardımcı olabilir misiniz?
    • soklakov
      Определение нагрузок на фундамент

      By soklakov · Posted

      10^18 частиц. Да, мощнее имеющихся.
    • Golem
      Hockert cw 1250

      By Golem · Posted

      Спасибо большое, я с таким не сталкивался.
    • VOleg
      Как на листе детали поставить ссылку на номер позиции, под которой она записана в сборке?

      By VOleg · Posted

      Я бы сказал "Странный вопрос". За 40+ лет ни разу на чертеж детали не ставил номера позиций т.к. ни разу не нужно было. Но вот последние несколько лет люди начали спрашивать... Видимо на производство пришли "новые люди"... Вот, помогаю, как могу. Еще надо на листах металла на складе писать на какие детали куплен этот металл... И на трубы, что ли? А на сверлах писать какие отверстия нужно будет просверлить...
    • NGM
      Проверка чертежа технологом

      By NGM · Posted

      Смотрю стандарт, и понимаю - сколько человек его прочтёт, столько трактовок и будет. Сплошные абстрактные формулировки; целый раздел про аттестацию, начинающийся с того, что она "не носит обязательный характер"... Зато, если буквально читать требования к специалисту к нормоконтролю - какой-то терминатор получается: и национальную систему стандартизации знает (всю!), и опыт по разработке КД/ТД имеет, и английским языком владеет. Вот для чего нужны такие ГОСТ-ы?! Заготовка для локального СТО, не более.  
    • NGM
      Проверка чертежа технологом

      By NGM · Posted

      В идеальном мире розовых единорогов блюющих радугой подпись "Т. контр" в чертеже должен ставить технолог, который будет разрабатывать техпроцесс на проектируемое изделие. На практике я такого не встречал никогда. Предприятиям, где такое практикуется - моё заочное уважение. @Фарита, расслабьтесь. Лучшее, что Вы можете сделать у себя - организовать собственную группу конструкторов, которые будут: а) перечерчивать полученную КД под требования вашего производства, б) выносить мозги разработчику бесконечными ПИ. Поверьте, так работает абсолютное большинство предприятий промышленности - почти везде есть подразделения "конструкторского сопровождения производства", по факту являющиеся живыми примерами того, что ни ЕСКД, ни ЕСТД не работают и никогда не работали (во всяком случае, первое слово из обеих абревиатур).
×
×
  • Create New...