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

[AlexKaz 1 828]

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

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

 

Способ 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 828]

Подробнее про вариант 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

 

[ДОБРЯК 606]

19.05.2022 в 11:52, Fedor сказал:

Да мне математику плевать на какое то жалкое приложение математики вроде физики. Мы же не физические эксперименты обсуждаем, а решение чисто математических задач

@SergeyGertsik вот почитайте. Гениальный математик придумал математическую задачу и хочет чтобы ее решили. Сам он не может решить.

Правильного решения с граничными условиями которые придумал математик нет. А то что у задачи нет физического смысла гениальному математику плевать

[статист 609]

7 минут назад, SergeyGertsik сказал:

Да вы знаток как все устроено во ВНИИЭФ, знаете что сделают, что не сделают и в какие сроки

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

Я на вскидку даю, что в лучшем случае эта опция появиться года через два, если вы зададитесь целью. А если не зададитесь, то она так и будет болтаться в обсуждениях, а тех поддержка, как вы, будет доказывать, что да оно и не надо =)

 

15 минут назад, SergeyGertsik сказал:

Эту, только без W. Могу посчитать в статике с W. Только вопрос: F следящая нагрузка, или мертвая?

Посчитайте, пожалуйста.

[Fedor 1 604]

Давайте считать что на веревочке подвешен груз создающий силу... 

[SergeyGertsik 64]

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

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

Я на вскидку даю, что в лучшем случае эта опция появиться года через два, если вы зададитесь целью. А если не зададитесь, то она так и будет болтаться в обсуждениях, а тех поддержка, как вы, будет доказывать, что да оно и не надо =)

А у кого препост допилен настолько, что перестали его дорабатывать?

Решение задачи с геом нелинейностью

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

@SergeyGertsik вот почитайте. Гениальный математик придумал математическую задачу и хочет чтобы ее решили. Сам он не может решить.

Правильного решения с граничными условиями которые придумал математик нет. А то что у задачи нет физического смысла гениальному математику плевать

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

Даже не подрастянули, по Х не фиксировал.

Это на этапе приложения кинематики только

Вот так, если зафиксировать, как в постановке задачи

[статист 609]

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

А у кого препост допилен настолько, что перестали его дорабатывать?

Решение задачи с геом нелинейностью

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

Даже не подрастянули, по Х не фиксировал.

Это на этапе приложения кинематики только

 

 

Вы с помощью пластин решали? Если не сложно, посчитайте с солидами.

[SergeyGertsik 64]

С пластинами. Как дойдут руки , посчитаю.

Результат без геом. нелинейности

[Fedor 1 604]

Цитата

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

Это не спортивно. Интересно когда сразу все в рамках одной задачи. Когда-то это считал тестом для больших деформаций... 

[статист 609]

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

С пластинами. Как дойдут руки , посчитаю.

Года через два? :)

[ДОБРЯК 606]

16 минут назад, SergeyGertsik сказал:

Не, ну у этой задачи все понятно, мы наклонили и растянули линейку

Вы же не читали предыдущие десять страниц.)

Если вы ее растянули почему нулевые деформации?

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

Вы с помощью пластин решали? Если не сложно, посчитайте с солидами.

Разницы не будет. Это ПНС задача. Две степени свободы в узле. И если стержнями моделировать 2 степени свободы в узле.

[SergeyGertsik 64]

11 минут назад, Fedor сказал:

Это не спортивно. Интересно когда сразу все в рамках одной задачи. Когда-то это считал тестом для больших деформаций... 

Пф, одновременно

[статист 609]

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

Разницы не будет. Это ПНС задача. Две степени свободы в узле. И если стержнями моделировать 2 степени свободы в узле.

Разница во времени. Если считать ПНС в явной постановке, то это несколько минут. Если Солидами в явной постановке - до часа.

[SergeyGertsik 64]

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

Года через два? :)

Возможно. Я волен сам выбирать себе сроки )

[статист 609]

1 минуту назад, SergeyGertsik сказал:

Возможно. Я волен сам выбирать себе сроки )

У Вас же 3000 ядер под рукой. Это совсем не сложно.

[SergeyGertsik 64]

Только что, статист сказал:

Разница во времени. Если считать ПНС в явной постановке, то это несколько минут. Если Солидами в явной постановке - до часа.

Не согласен. Шаг по времени зависит от характерного линейного размера, если он одинаков, то не шибко будет отличаться. В оболочечном 4 узла с 6 степенями свободы - 24 DOF, в солиде 8 узлов и 3 степени свободы - 24 DOF

Только что, статист сказал:

У Вас же 3000 ядер под рукой. Это совсем не сложно.

Вы полагаете, что количество вычислительных ресурсов влияет на время создания расчетной модели? Мол если специалист делает модель на пентиум 2 с одним ядром на час, то специалист с 3000 ядер под рукой должен сделать модель за 1,2 секунды?

[Fedor 1 604]

 

Видите нулевые напряжения, соответственно и нулевые деформации. 

 

"Пф, одновременно" - симпатично, спасибо :) 

Изменено 20 мая 2022 пользователем Fedor

[статист 609]

5 минут назад, SergeyGertsik сказал:

Не согласен. Шаг по времени зависит от характерного линейного размера, если он одинаков, то не шибко будет отличаться. В оболочечном 4 узла с 6 степенями свободы - 24 DOF, в солиде 8 узлов и 3 степени свободы - 24 DOF

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

В ПНС две степени свободы, Добряк об этом выше писал.

 

7 минут назад, SergeyGertsik сказал:

Вы полагаете, что количество вычислительных ресурсов влияет на время создания расчетной модели? Мол если специалист делает модель на пентиум 2 с одним ядром на час, то специалист с 3000 ядер под рукой должен сделать модель за 1,2 секунды?

Вы эту модель сделаете за час?

 

[SergeyGertsik 64]

2 минуты назад, Fedor сказал:

 

Видите нулевые напряжения, соответственно и нулевые деформации. 

Так откуда им там взяться, если в той задаче я X не зафиксировал? Я на картинке линейные размеры привел, они не поменялись. Просто поворот.

[ДОБРЯК 606]

11 минут назад, SergeyGertsik сказал:

Просто поворот.

Как же просто поворот. Смотрите внимательно картинку. Коэф. деформаций = 1.

@SergeyGertsik вы опять на скорость решаете задачу?)

[SergeyGertsik 64]

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

Как же просто поворот. Смотрите внимательно картинку. Коэф. деформаций = 1.

Что такое коэффициент деформаций? 

[ДОБРЯК 606]

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

Вы эту модель сделаете за час?

Если есть модель из оболочек модель делается за 1 мин.