Софт ЛОГОС - что за зверь? В наших краях не водится

[AlexKaz 1 851]

Поскольку геморрой с созданием объёмной сетки на сложных телах продолжается и конца ему не видно, есть лайфхак, который работает чуть чаще.

Скрытый текст

 

Способ 0 - это не способ, а убогое признание, что генератор (как конкретный инструмент импортозамещения, как следствие) толком не работает:

- создать сетку в платных забугорных WB, Patran, Nastran, Abaqus etc., экспортировать в формат, который переварит препост. По факту, это не решение проблемы легальным способом, поскольку лицензий на эти софтины фирмам сегодня не купить. А кто купил бессрочные лицензии раньше, тому замещение вообще зачем?

 

Способ 1:

- создать объёмную сетку в бесплатных SALOME-Meca, GMSH, FreeCAD, LS-Prepost, OpenFOAM, импортировать в препост, в том числе через ряд конверсий. Кто понимает минимальный технический английский со словарём, тот уже разобрался. Кто не хочет использовать гуглтранслейт и моск, тот ноет что ничегонеполучаицца.

 

Способ 2, самый хитрожопый на текущий момент, поэтому нигде не описан:

- вариант 1) не самый хороший. Получить корректный замкнутый STL в любом CAD, затолкать в препост-прочность, затем "Инструменты -> Сетка -> Генерация объёмной сетки -> Генерация тетраэдральной сетки от поверхностной сетки". Проверить, заполнилось ли внутреннее пространство сетки элементами, использовать инструмент по Ctrl+F8. Затем удалять поверхностную сетку после создания объёмной. 

- вариант 2) лучше. Построить качественную замкнутую поверхностную сетку в SALOME-Meca, GMSH, FreeCAD, LS-Prepost, экспортировать в STL. Повторить вариант 1).

- вариант 3) почти идеальный, но работает не всегда. Импортировать в Аэрогидро-Препост геометрию, построить сетку поверхностную/объёмную, сохранить проект. В препост-прочность импортировать efr-сетку из подпапок проекта аэрогидро. При необходимости, повторить пункт 1). Способ плох тем, что Аэрогидро - отдельный платный модуль.

- вариант 4) разные комбинации вариантов 1-3 и способов 0 и 1.

 

Способ 3) самый спорный. Продолжать использовать только мэшер препоста. 

 

Наконец, кардинальный способ 4) посмотреть в сторону опенсорц. См. ремарку про гуглтранслейт и моск.

 

 

[AlexKaz 1 851]

Подробнее про вариант 3) создания объёмной сетки в  прочность на очень сложной для препоста геометрии. Способ крайне сложный и многоэтапный, что только подтверждает, что в прочности надо конкретно улучшать ситуацию с мэшером.

Скрытый текст

 

1 сначала в Препост Аэрогидро импортируется геометрия; сложная геометрия как правило подгружается с дефектами, поэтому для неё применяется п.2 ниже
- отчасти дефекты, особенно связанные со скруглениями, могут быть устранены варьированием точности и размерности в настройках при импорте геометрии, но это не панацея в общем случае, к тому же может очень значительно увеличивать время импорта
- геометрию предпочтительнее подгружать по одному телу, и сетку соответственно строить только для одного редактируемого выбранного тела, но в целом на результат не сильно влияет, если учитывать точность при импорте

 

2 открытая в текущий момент в рабочем поле геометрия лечится:
- для начала, геометрия осматривается визуально на дефекты
- затем подгружается окно "Диагностика геометрии" (на панели зелёный кубик с жёлтым треугольником с восклицательным знаком)
- в окне наиболее значимые пункты по степени важности - это открытые контуры, затем перекрытия
- т.е. сначала пытаемся автоматикой (ПКМ по элементу дерева -> Автоисправление)  убрать все проблемные открытые контуры, затем все перекрытия - часто этого достаточно
- если недостаточно:
- удаляются битые поверхности
- при необходимости кривые (в т.ч. на поверхностях) разбиваются инструментом "разбить ребро" - очень удобная фича (в LS-PREPOST тоже имеется)
- создаются прямые/кривые/дуги/скругления по точкам
- по кривым и прямым создаются новые поверхности
- поверхности в ряде случаев требуется пересоздать: "Режим редактирования поверхностей -> Пересоздать грань". особенно это касается скруглений, которые требуется пересоздавать заново по соседним граням чтобы получить коррекную форму скругления

3 геометрия сшивается (жёлтая кнопка с молнией "Построение топологии", затем опционально на панели "Выбрать область выделения-> Прямоугольник" и выделить всю деталь) сначала с точностью по умолчанию, если не получается - вручную поварьировать параметр в сторону увеличения абсолютной величины
- делать это нужно даже если окажется, что проблем в "Диагностике геометрии" нет, а все поверхности кажутся сшитыми и автоматика не подсвечивает красных контуров (см. ниже), при этом сшивать желательно с параметром точности, выставленным вручную, значение выбрать побольше
- если геометрия не сшита - визуально наблюдаются красные контуры на геометрии; если они остаются даже при варьировании параметра при сшивке до большой величины и даже после перезапуска и обновления окна инструмента "Диагностика геометрии" - требуется вернуться на п.2. при этом, желательно долечить геометрию в окне "диагностика геометрии" по всем пунктам дерева. скорее всего, потребуется заново удалить и создать битые или некачественные поверхности, а также пересоздать их соответствующим автоматическим инструментом 
- если и сейчас не удаётся добиться сшивки - используется последний довод - поверхности нарезаются инструментами разрезки граней (их несколько) и всё возвращается на п.2 

4 создаётся регион:
- в этом пункте можно поступить двояко, либо выбрать в дереве геометрию, ПКМ -> создать регион
- либо выбрать в дереве геометрию -> извлечь тесселяцию с флагом автосшивки, что создаст в дереве новую деталь (фактически, голый stl) - выбрать её и создать для неё регион
- для тесселяции включаем видимость сетки stl - на панели "Стиль отображения геометрии -> Закрашенный с рёбрами"
- при создании региона выбрать "Разбиение региона на границы -> НЕ РАЗБИВАТЬ", что часто избавляет от генерации элементов с нулевой площадью, обнаружить которые можно только постфактум на этапе создание объемной сетки, при этом возникает ошибка "Обнаружена граница 'id=...,' с нулевым количеством треугольников - объемная сетка не может быть построена
- для тесселяции надо убедиться, что все поверхностные элементы stl сшиты, т.е. все грани стыкуются корректно по узлам; если этого где-то не произошло - возвращаемся на п.2. или используем корректировку stl инструментами работы с поверхностной сеткой ("режим редактирования фасеточной геометрии") - это крайне нежелательно делать, т.к. не позволит не вылечит саму геометрию, соответственно, и сохранить её можно будет только битой

5 строится поверхностная сетка (довольно плотная):
- в ветке генераторов отмечается "поверхностная сетка" и "ТЕТРАЭДРАЛЬНАЯ СЕТКА" (это важно)
- для выбора размера сетки необходимо определиться с размером геометрии - включить "Вид -> Масштабная линейка", затем задать ориентировочный размер элемента в "Блок генераторов -> Базовый размер"
- на панели выбирается кнопка "плоская сетка с молнией" - запуск генератора поверхностной сетки
- если поверхностная сетка не построилась, что бывает, возвращаемся на п.2.
- если поверхностная сетка построилась - ЗАПУСКАЕТСЯ ГЕНЕРАТОР ОБЪЁМНОЙ СЕТКИ
- если объёмная сетка построилась - вся мегаопупея с гарантией ~ 99.99% завершится хорошо

6 если поверхностная сетка построилась, но не построилась объёмная:
- либо тесселированная поверхностная сетка экспортируется в stl (ПКМ по названию тесселированной детали -> Экспортировать геометрию -> stl)
- либо проект Аэрогидро сохраняется, открывается Препост Прочность, в нём импортируется поверхностная сетка из подпапок проекта аэрогидро (efr-сетка с названием *surface*), затем сохраняется в stl
- сохранённый stl-файл загружается КАК ГЕОМЕТРИЯ в НОВОМ проекте Аэрогидро, затем сетка перестраивается, т.е. повторяется п.4 и 5

7 строится объёмная сетка в Препост Прочность:
- открывается Препост Прочность, в нём импортируется поверхностная сетка из подпапок проекта аэрогидро (efr-сетка с названием *surface*) либо импортируется поверхностный stl
- "Инструменты -> Сетка -> Генерация объёмной сетки -> Генерация тетраэдральной сетки от поверхностной сетки". Проверить, заполнилось ли внутреннее пространство сетки элементами, использовать инструмент по Ctrl+F8. Затем удалять поверхностную сетку после создания объёмной.

 

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

 

Дополнительно в процессе как вспомогательный инструмент можно использовать LS-PREPOST, например для проверки корректности stl-сетки, но в общём случае это не обязательно. В LS-Prepost удобно визуально наблюдать, есть ли в сетки незакрытые контуры, но корректировать, закрывая контуры, бесполезно - препосты логоса объёмную сетку строить всё равно не захотят, ссылаясь на свои непонятные проблемы.
Также в LS-PREPOST удобно строить объёмную сетку из тетраэдров по качественной поверхностной сетки из п.6, при этом удобно быстро понижать количество элементов (метод 2 при создании тетрасетки) при довольно хорошем сохранении границ сетки, но не идеально всё же; затем сохранить в k-файл и открыть объёмную сетку в Препост Прочность.
Также в LS-PREPOST в ряде случаев удаётся долечить геометрию и сшить поверхности до солида, если это всё же не получается в Аэрогидро-Препост.
Ещё один лайфхак - битая stl-сетка импортируется в SALOME-Meca, близкие узлы сшиваются, линейная сетка заменяется на квадратичную, сохраняется в stl и загружается в ls-prepost - далее можно получить объёмную сетку из тетраэдров.

 

 

Ранее в сериале:

Скрытый текст

07.11.2023 в 08:34, AlexKaz сказал:

Поскольку геморрой с созданием объёмной сетки на сложных телах продолжается и конца ему не видно, есть лайфхак, который работает чуть чаще.

 

Способ 0 - это не способ, а убогое признание, что генератор (как конкретный инструмент импортозамещения, как следствие) толком не работает:

- создать сетку в платных забугорных WB, Patran, Nastran, Abaqus etc., экспортировать в формат, который переварит препост. По факту, это не решение проблемы легальным способом, поскольку лицензий на эти софтины фирмам сегодня не купить. А кто купил бессрочные лицензии раньше, тому замещение вообще зачем?

 

Способ 1:

- создать объёмную сетку в бесплатных SALOME-Meca, GMSH, FreeCAD, LS-Prepost, OpenFOAM, импортировать в препост, в том числе через ряд конверсий. Кто понимает минимальный технический английский со словарём, тот уже разобрался. Кто не хочет использовать гуглтранслейт и моск, тот ноет что ничегонеполучаицца.

 

Способ 2, самый хитрожопый на текущий момент, поэтому нигде не описан:

- вариант 1) не самый хороший. Получить корректный замкнутый STL в любом CAD, затолкать в препост-прочность, затем "Инструменты -> Сетка -> Генерация объёмной сетки -> Генерация тетраэдральной сетки от поверхностной сетки". Проверить, заполнилось ли внутреннее пространство сетки элементами, использовать инструмент по Ctrl+F8. Затем удалять поверхностную сетку после создания объёмной. 

- вариант 2) лучше. Построить качественную замкнутую поверхностную сетку в SALOME-Meca, GMSH, FreeCAD, LS-Prepost, экспортировать в STL. Повторить вариант 1).

- вариант 3) почти идеальный, но работает не всегда. Импортировать в Аэрогидро-Препост геометрию, построить сетку поверхностную/объёмную, сохранить проект. В препост-прочность импортировать efr-сетку из подпапок проекта аэрогидро. При необходимости, повторить пункт 1). Способ плох тем, что Аэрогидро - отдельный платный модуль.

- вариант 4) разные комбинации вариантов 1-3 и способов 0 и 1.

 

Способ 3) самый спорный. Продолжать использовать только мэшер препоста. 

 

Наконец, кардинальный способ 4) посмотреть в сторону опенсорц. См. ремарку про гуглтранслейт и моск.

Изменено 7 ноября пользователем AlexKa

 

[ДОБРЯК 615]

6 часов назад, Jesse сказал:

НЕ ВЕРЮ!

Трудно поверить, что отечественный генератор работает четко, но это факт. Два элемента по всей длине. 

[ДОБРЯК 615]

@Jesse теперь вы решите тестовую задачу. Верхнюю и нижнюю полку нужно также разбить на 2 элемента.)

Читаете STEP формат. Ничего лечить не надо. Только генерация.

https://disk.yandex.ru/d/KH8VeJbQgXLSEw

@AlexKaz и вы проверьте Логос. Если у вас получится сгенерировать этот тест, то проверим скорость решения задач Динамики и прочности.

[AlexKaz 1 851]

Да можно то можно, было бы желание, и, главное, время. Ниже 3 элемента в полке, можно и с двумя, но потребуется заново перестраивать мелкую сетку на поверхностях отверстий. Даже если выделять их прямоугольником через Ctrl + ЛКМ, то выбираются, то не выбираются. Видимо, требуется скрывать поверхность спинки. В WB, LS-Prepost и SALOME я б так и делал, а тут...

 

Ага, Ctrl+F8. И тогда получается.

Изменено 16 февраля 2023 пользователем AlexKaz

[Jesse 979]

1 час назад, ДОБРЯК сказал:

@Jesse теперь вы решите тестовую задачу. Верхнюю и нижнюю полку нужно также разбить на 2 элемента.)

 

Читаете STEP формат. Ничего лечить не надо. Только генерация.

https://disk.yandex.ru/d/KH8VeJbQgXLSEw

 

а нормальные тестовые примеры придумать слабо? как я выше 2Д змейку.
Вот делать мне больше нечего на 10 лям элементов генерить тестовый пример

[AlexKaz 1 851]

Ну и вот, например,

[ДОБРЯК 615]

22 минуты назад, AlexKaz сказал:

но потребуется заново перестраивать мелкую сетку на поверхностях отверстий. Даже если выделять их прямоугольником через Ctrl + ЛКМ, то выбираются, то не выбираются.

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

8 минут назад, Jesse сказал:

Вот делать мне больше нечего на 10 лям элементов генерить тестовый пример

А меньше 10 000 000 никак не получается? ;=)

У меня получилась модель 

Количество узлов=215519

Количество элементов = 857434
Нормальный тест для современных генераторов и компьютеров.

Изменено 16 февраля 2023 пользователем ДОБРЯК

[AlexKaz 1 851]

4 минуты назад, ДОБРЯК сказал:

Поясните пожалуйста...

Я писал, что он требует доработки, нужны ручные манипуляции 100 студентов за еду, потеря времени. По крайней мере, в версии от 2022. Тем не менее, пашет.

[ДОБРЯК 615]

9 минут назад, AlexKaz сказал:

Тем не менее, пашет.

И какая сетка получилась? Сколько узлов, элементов? И время генерации этой сетки...

[AlexKaz 1 851]

В Логос с моими настройками ~155 сек, 1.68 млн элементов, 520 тыс. узлов.

 

Вообще говоря, ИСПАвский генератор мне напоминает Gmsh. Возможно, это даже одно и тоже? =)

Вот сетка, созданная в 1 клик во FreeCAD, куда прикручен Gmsh.

Изменено 16 февраля 2023 пользователем AlexKaz

[ДОБРЯК 615]

4 минуты назад, AlexKaz сказал:

Вообще говоря, ИСПАвский генератор мне напоминает Gmsh. Возможно, это даже одно и тоже?

А вы сделайте сетку в FEMAP или в NX и сравните. :=) Я сравнивал сетку в NX и в ИСПА практически один в один.

 А генератор Gmsh в исходных кодах где-то можно скачать? И на каком языке он написан?

[AlexKaz 1 851]

Прочнистов устроит любой адекватный вариант сетки. К слову, возможен импорт таких сеток: STL *.stl, VTK *.vtk и *.vtu, LS-DYNA *.k, ABAQUS *.inp, CGNS *.cga и *.cgns, NASTRAN *.bdf и *.nas, ANSYS *.cdb, COSMOS *.geo, EFR *.efr, NGEOM *.ngeom, MESH *.msh, FLUENT MSH *.msh, TEC *.tec.

Может быть, ещё что-нибудь прикрутили в 2023-м.

 

А С++-сорцы gmsh тут https://gmsh.info/src/

[Jesse 979]

2 часа назад, ДОБРЯК сказал:

Нормальный тест для современных генераторов и компьютеров.

мой пример на 1700 оболочечных КЭ. у вас почти на 3 порядка больше. Разницу ощущаете?)

[ДОБРЯК 615]

2 часа назад, AlexKaz сказал:

К слову, возможен импорт таких сеток

А экспорт сетки в каких форматах?

[AlexKaz 1 851]

17 часов назад, ДОБРЯК сказал:

А экспорт сетки в каких форматах?

STL *.stl, VTK *.vtk, LS-DYNA *.k, ABAQUS *.inp, NASTRAN *.bdf, ANSYS *.cdb, EFR *.efr, NGEOM *.ngeom, MESH *.msh, FLUENT MSH *.msh.

[ДОБРЯК 615]

40 минут назад, AlexKaz сказал:

NASTRAN *.bdf

Запишите в этот формат, то что вы сгенерировали в Логосе. Если вам не сложно...

[AlexKaz 1 851]

27 минут назад, ДОБРЯК сказал:

Запишите в этот формат

Скрытый текст

 

Простая окружность

BEGIN BULK
$
GRID    1               1.0000000.0000001.000000
GRID    2               0.8660250.5000001.000000
GRID    3               0.5000000.8660251.000000
GRID    4               0.0000001.0000001.000000
GRID*   5                               -0.500000       0.866025        
*       1.000000        
GRID*   6                               -0.866025       0.500000        
*       1.000000        
GRID*   7                               -1.000000       0.000000        
*       1.000000        
GRID*   8                               -0.866025       -0.500000       
*       1.000000        
GRID*   9                               -0.500000       -0.866025       
*       1.000000        
GRID*   10                              -0.000000       -1.000000       
*       1.000000        
GRID*   11                              0.500000        -0.866025       
*       1.000000        
GRID*   12                              0.866025        -0.500000       
*       1.000000        
GRID    13              0.1058930.2056381.000000
GRID*   14                              -0.088728       -0.576253       
*       1.000000        
GRID*   15                              0.389229        -0.039023       
*       1.000000        
GRID*   16                              -0.514938       -0.007600       
*       1.000000        
GRID*   17                              -0.253269       0.370476        
*       1.000000        
GRID*   18                              -0.345289       -0.325252       
*       1.000000        
CQUAD4         1       1       2       3      13      15
CQUAD4         2       1       6       7      16      17
CQUAD4         3       1      14       9      10      11
CQUAD4         4       1      11      12      15      14
CQUAD4         5       1      12       1       2      15
CQUAD4         6       1       8       9      14      18
CQUAD4         7       1      14      15      13      18
CQUAD4         8       1      13      17      16      18
CQUAD4         9       1      16       7       8      18
CQUAD4        10       1       6      17       4       5
CQUAD4        11       1       3       4      17      13
ENDDATA

Простой куб

BEGIN BULK
$
GRID    1               0.2480420.0000000.755814
GRID    2               0.0000000.0000001.000000
GRID    3               0.0000000.0000000.515423
GRID    4               0.5382110.0000000.487699
GRID    5               1.0000000.0000000.000000
GRID    6               1.0000000.0000000.535381
GRID    7               0.5000000.0000001.000000
GRID    8               1.0000000.0000001.000000
GRID    9               0.0000000.0000000.000000
GRID    10              0.4955910.0000000.000000
GRID    11              0.0000001.0000000.508034
GRID    12              0.0000001.0000000.000000
GRID    13              0.0000000.5000000.000000
GRID    14              0.0000000.4922710.507265
GRID    15              0.0000000.4646191.000000
GRID    16              0.0000001.0000001.000000
GRID    17              0.6052930.6341840.000000
GRID    18              0.2483440.2494460.000000
GRID    19              0.4870331.0000000.000000
GRID    20              1.0000001.0000000.000000
GRID    21              0.7010920.2812640.000000
GRID    22              1.0000000.5094900.000000
GRID    23              1.0000000.5403480.481343
GRID    24              1.0000000.4758061.000000
GRID    25              1.0000001.0000000.531377
GRID    26              1.0000001.0000001.000000
GRID    27              0.4611551.0000000.473785
GRID    28              0.7368781.0000000.756286
GRID    29              0.4704401.0000001.000000
GRID    30              0.5276550.4595561.000000
GRID    31              0.2463800.7283011.000000
CTETRA         1       1      27      17      19      20
CTETRA         2       1      21      18       4      17
CTETRA         3       1      23      21      17      22
CTETRA         4       1      12      13      19      27
CTETRA         5       1      23      25      17      27
CTETRA         6       1      25      27      20      17
CTETRA         7       1      27      14      16      11
CTETRA         8       1      23      22      17      20
CTETRA         9       1       4      24       8       6
CTETRA        10       1      23       4       5       6
CTETRA        11       1       4      23       5      21
CTETRA        12       1      14      27      16      31
CTETRA        13       1      23      27      28      25
CTETRA        14       1      29      31      16      27
CTETRA        15       1      30      27      29      28
CTETRA        16       1      10       4      18       9
CTETRA        17       1      18      17      14       4
CTETRA        18       1      29      26      30      28
CTETRA        19       1      27      30      29      31
CTETRA        20       1      23      28      26      25
CTETRA        21       1      12      27      11      13
CTETRA        22       1      15       2       1      30
CTETRA        23       1      10      21      18       4
CTETRA        24       1      13      17      14      18
CTETRA        25       1      10       5      21       4
CTETRA        26       1      15      14      30       1
CTETRA        27       1      15      14      16      31
CTETRA        28       1      27      31      14      30
CTETRA        29       1      23      26      30      24
CTETRA        30       1      23       4      24      30
CTETRA        31       1      26      23      30      28
CTETRA        32       1      23      14      17       4
CTETRA        33       1      23      17      14      27
CTETRA        34       1       7       1      30       4
CTETRA        35       1       4      23      24       6
CTETRA        36       1       2       7       1      30
CTETRA        37       1      23      30      14       4
CTETRA        38       1       1       4      14      30
CTETRA        39       1       9      13      14      18
CTETRA        40       1       3       2       1      15
CTETRA        41       1      24       4       8      30
CTETRA        42       1      14       9       4       3
CTETRA        43       1      27      14      11      13
CTETRA        44       1       3       1      14      15
CTETRA        45       1      15      30      14      31
CTETRA        46       1       3       4      14       1
CTETRA        47       1      27      19      17      13
CTETRA        48       1      23      14      30      27
CTETRA        49       1      27      17      14      13
CTETRA        50       1       9      14       4      18
CTETRA        51       1      23       4      17      21
CTETRA        52       1      23       5      21      22
CTETRA        53       1      23      20      17      25
CTETRA        54       1       8      30       4       7
CTETRA        55       1      23      27      30      28
ENDDATA

 

 

[ДОБРЯК 615]

22 часа назад, AlexKaz сказал:

Ну и вот, например,

Я просил записать эту модель.

[AlexKaz 1 851]

Нет, побережём интернет и не будем плодить сущности.

[soklakov 1 485]

15.02.2023 в 16:49, Jesse сказал:

ну тогда попробуйте у себя также сделать.

Вот пример "змейка" шириной 2 мм

А че не гексами?))

 

15.02.2023 в 07:40, AlexKaz сказал:

Единственное, хочется, чтобы она работала на солидах и не только мельчила сетку, но и укрупняла. Нужно это в том числе для ручной борьбы с бессмысленным мельчением сетки в прочности.

Не выглядит ли более логичным глобально бить крупно и мелчить где надо?

Ну то есть отменять это бессмысленное измельчение ещё на входе 

Догадываюсь ... Невозможно?

[Jesse 979]

19 минут назад, soklakov сказал:

А че не гексами?))

нетю)
а так не против лишний раз посмотреть  на криволинейной геометрии в Ансисе гекс-сетку строят..)
Вперёд, Соклаков. Отряхни уже свой АНсис от нафталина..))

 

20 минут назад, soklakov сказал:

Не выглядит ли более логичным глобально бить крупно и мелчить где надо?

а представь у тебя много разных мелких фичеров. Хочется глобально на них сделать мелкую сетку (чтоб не выбирать), а выборочно большие фичеры сделать крупную сетку..