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

[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

 

[Jesse 979]

Был на инженерной конфе в Кронштадте 5 окт, представитель АО "ЦИФРА" гордо заявил, что они закупили ЛОГОС, обучают работе в нём молодых сотрудников и даже проводят соревнования в нём. На вопрос из зала в чём считали конструкции которые есть в презентации, Дмитрий Михалюк гордо заявил, что в Ansys, ибо по уровню автоматизации и удобству ему пока нет равных

[Radon 17]

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

Был на инженерной конфе в Кронштадте 5 окт, представитель АО "ЦИФРА" гордо заявил, что они закупили ЛОГОС, обучают работе в нём молодых сотрудников и даже проводят соревнования в нём. На вопрос из зала в чём считали конструкции которые есть в презентации, Дмитрий Михалюк гордо заявил, что в Ansys, ибо по уровню автоматизации и удобству ему пока нет равных

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

[athrips 1]

Скажите пожалуйста, в Логос-прочность есть возможность активировать\деактивировать элементы по ходу расчета?

[AlexKaz 1 851]

@athrips , сейчас нет. Сколько-то времени назад запрос на реализацию был отправлен в техподдержку. Когда именно доработают, разработчики пока не сообщали.

Если совсем заморачиваться, в Логос такое можно при большом желании реализовать через USERMAT, по аналогии с ANSYS - делая Юнг околонулевым или принудительно зануляя тензор напряжений по ходу расчёта опционально в нужных конечных элементах.

4 часа назад, Radon сказал:

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

Есть операция "Морфинг узлов", переносит узлы с одной поверхности на другую и присваивает узлы геометрии чтобы дальше поверхностный и объёмный мэшер смог построить согласованную сетку. Работает на простой геометрии, в т.ч. циклической, но на сложной отказывается. Аналогично, ждём доработки.

4 часа назад, Radon сказал:

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

Да, при этом критерий "искажённости" не всегда понятен. Это могут быть как элементы с острыми углами, так и с обратным порядком узлов. В общем случае, "лечится" выставлением флага "Расчёт с искажёнными ячейками" в дереве проекта. Вроде счёт с этой фичей корректный. Но хотелось бы, конечно, без сюрпризов с искажениями.

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

[AlexKaz 1 851]

3 часа назад, athrips сказал:

активировать\деактивировать элементы по ходу расчета?

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

[AlexKaz 1 851]

12 часов назад, AlexKaz сказал:

В явных расчётах в материалах есть критерии разрушения.

О, оказывается, в неявных тоже есть. Ни разу не юзал в неявных.

[Jesse 979]

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

О, оказывается, в неявных тоже есть. Ни разу не юзал в неявных.

Явный решатель тоже есть в Логосе?)

Ну,раз такая пляска пошла, то чего там по набору КЭ и моделей материалов? Гексы есть? Гиперупругие, модели пластичности навороченные тоже имеются?)

[athrips 1]

Есть труба сложной формы, в которую засыпают однотипные детали и нужно проверить, застрянут они там или нет. В первую очередь мне в этой задаче непонятно, что делать с деталями, которые сквозь трубу прошли - выглядит целесообразным их либо телепортировать обратно на вход, либо как-нибудь устранить из расчета. Удалять или перемещать подобласти макросом можно, но как отбирать их по координате, да ещё проделывать это в процессе расчета - не понимаю. Или я вообще не с той стороны подхожу к проблеме?

[AlexKaz 1 851]

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

Или я вообще не с той стороны подхожу к проблеме?

Разве что в явном решателе перейти в SPH-постановку, где можно указать габаритную коробку. За пределами коробки SPH-частицы в расчёт не берутся. Но я не проверял будет ли сохраняться корректно геометрия деталей в SPH в процессе.

1 час назад, Jesse сказал:

по набору КЭ и моделей материалов? Гексы есть? Гиперупругие, модели пластичности навороченные тоже имеются?

Да. Да. Да. Универу не хватает разве что связки "гиперупругий с вязкоупругостью" чтобы посмотреть релаксацию резин. Но? при желании? универ в принципе должен сам осилить написать подходящий usermat.

[AlexKaz 1 851]

1 час назад, athrips сказал:

Удалять или перемещать подобласти макросом можно, но как отбирать их по координате, да ещё проделывать это в процессе расчета - не понимаю.

Если бы речь шла про LS-DYNA, в ней есть рестарты, в которых можно гасить и перемещать part'ы произвольно.

Если про Логос - без понятия. М.б. разбить процесс на стадии, результат расчёта каждой i-й стадии сохранять, деформированную сетку сэйвить в k-файл (или efr, что проще, т.к. экспорт в k-файл должен работать только в свежих релизах 2024 или в 2025 г.), а его уже модить как угодно, затем закидывать сетку в стадию i+1. Точно не проверял, можно ли при таком подходит сохранять в нужных part'ах скорости и НДС чтобы перекидывать их от стадии к стадии.

Изменено 17 октября пользователем AlexKaz

[athrips 1]

@AlexKaz , спасибо! Насчет SPH-постановки: в свойствах подобластей есть параметры "множитель длины сглаживания", "коэффициент квадратичной вязкости" и "множитель минимальной плотности" - вы не подскажете, за что они отвечают? Предполагаю, что это что-то про размер области влияния сглаживающей функции и про искуственную вязкость, но в документации при описании теории эти параметры вообще не упоминаются.

[AlexKaz 1 851]

@athrips , ориентируйтесь на Хэлп к LS-DYNA, в нём к SPH-картам есть соответствующие описания.

И ещё по удалению тел в Логос: в справке есть упоминание, что (произвольные?) тела можно удалять указанием наборов в command.xml, плюс потребуется в дереве расчёта указать необходимость сохранения файла рестарта. Порядок запуска на счёт в рестарте в справке был описан.

[AlexKaz 1 851]

Видео «Определение НДС лопатки авиационного ГТД» в Прочность от Адванс Инжиниринг.
https://youtu.be/9Aq55JcNSnE

Источник https://advengineering.ru/ru/news/vebinar-logos-2024-10-31/ . 

Часть нюансов осталась за кадром, но в целом нормальный уровень.

[ДОБРЯК 615]

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

Часть нюансов осталась за кадром, но в целом нормальный уровень.

Это статическая задача. Какой контроль невязок?

Автор ролика поставил в очередь решение на супер-пупер компьютер у которрого 3 000 ядер.  И сразу остановил решение. Эти ядра сдаются в аренду. Каждый год нужно доплачивать...:=)

Это же смешно...:=)

[AlexKaz 1 851]

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

Какой контроль невязок?

Настроить можно по силе, перемещению, пластической деформации, числу итераций и ещё немного в глубоком погружении. Тут не очень принципиальны настройки сходимости, всё равно сойдётся, и сетка для обучения, а не для расчётов для КД.

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

И сразу остановил решение.

Задача с пластикой и в ANSYS за секунду не решится. Цель очевидна, показать типовые часто используемые инструменты в Препост, а не медитировать перед выдачей в лог.

Нормально же получилось =)

[ДОБРЯК 615]

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

Цель очевидна, показать типовые часто используемые инструменты в Препост

Естественно типовая задача. Все решают лопатки в этой постановке. И именно в пластике.:=)

Лопатки работают в упругой зоне. Такие ответственные конструкции работают в упругой зоне.

Это уже много раз обсуждали на форуме...

 

[AlexKaz 1 851]

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

Все решают лопатки в этой постановке. И именно в пластике.:=)

Конструкторам и прочнистам ГТД виднее, когда в такой постановке, а когда в другой. За всех сказать нельзя. Главное, сейчас инструмент позволяет переходить к более сложным случаям, чем на видео.

[ДОБРЯК 615]

13.11.2024 в 06:07, AlexKaz сказал:

Конструкторам и прочнистам ГТД виднее

Если конструкция или деталь работает в линейной зоне, то зачем реальным расчетчикам в реальной жизни, считать пластику? Или глубокую пластику... Вы себе ответьте на этот вопрос. :=)

Речь не идет о рекламе, или ливерной или докторской диссертации.

Речь идет о реальных расчетах. Интересно посмотреть как отличаются результаты при 20 град и при 200  град.

На обычном персональном компьютере... Тогда будет понятно, что это за зверь Логос...

 

 

[Fedor 1 630]

При высоких температурах проявляется ползучесть. Это и является ограничением для кпд по Карно...