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

[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

 

[Fedor 1 630]

Для линейного тензора деформаций правильно. Еще нарисуйте картинку деформаций. Например первого инварианта тензора деформаций.  Или напряжений  ... 

[SergeyGertsik 66]

У вас зафиксирован Х, у меня на референсной картинке нет. Просто опустил по Y, естественно, линейка не захотела растягиваться по Х

[статист 610]

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

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

Я знаю.

[SergeyGertsik 66]

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

Для линейного тензора деформаций правильно. Еще нарисуйте картинку деформаций. Например первого инварианта тензора деформаций.  Или напряжений  ... 

Зачем? они будут околонулевыми

[ДОБРЯК 615]

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

Еще нарисуйте картинку деформаций.

Я вам уже третий день объясняю, что и деформации и напряжения нулевые. Зачем рисовать?

[Fedor 1 630]

ЧТД.  Длина линейки при линейном тензоре деформаций увеличивается, а деформаций в дифференциальном смысле нет :) 

[SergeyGertsik 66]

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

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

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

Все верно, чем мельче сетка, тем меньше шаг - критерий Куранта-Фридриха-Леви.

Я считаю не в ПНС, а в трехмере

[ДОБРЯК 615]

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

Просто опустил по Y

Ну термины у вас. Как на зоне.)

Надо уравнение записать w = 200 мм.

[SergeyGertsik 66]

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

Я знаю.

А вот мне сейчас лень. У меня же статус "утомительно", да, я ленив

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

ЧТД.  Длина линейки при линейном тензоре деформаций увеличивается, а деформаций в дифференциальном смысле нет :) 

Не увеличится там ничего, я задал перемещение только по Y, по Х конец мог гулять, и естественно он шагнул так чтобы энергия деформирования была минимальной, в этом случае ~0

[ДОБРЯК 615]

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

Я считаю не в ПНС, а в трехмере

В Логосе и оболочка трехмерная?

Или в солидах считаете?

[SergeyGertsik 66]

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

Ну термины у вас. Как на зоне.)

Надо уравнение записать w = 200 мм.

Хм, для меня опустил - это сделал ниже по уровню. Если у вас ассоциации с зоной, то чтож...

Только что, ДОБРЯК сказал:

В Логосе и оболочка трехмерная?

Или в солидах считаете?

В локальных координатах элемента 2D, в рамках трехмерного расчета 3D

[статист 610]

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

Все верно, чем мельче сетка, тем меньше шаг - критерий Куранта-Фридриха-Леви.

Наверное, поэтому в МКЭ, где сеточная сходимость краеугольный камень, пользуются неявной схемой, чтобы решать задачи без оглядки на критерий шага.

Ну то есть чтобы решить простую задачу с помощью Логоса типа такой:

https://scadsoft.com/tests_scad/index/test/40

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

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

[ДОБРЯК 615]

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

В локальных координатах элемента 2D

Так вы и считаете данную задачу в плоскости. 

[SergeyGertsik 66]

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

Наверное, поэтому в МКЭ, где сеточная сходимость краеугольный камень, пользуются неявной схемой, чтобы решать задачи без оглядки на критерий шага.

Ну то есть чтобы решить простую задачу типа такой:

https://scadsoft.com/tests_scad/index/test/40

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

Явные методы применяются в своей области - падение, взрыв, пробитие. Посчитайте разрыв контейнера с ВВ при его детонации в неявной постановке

Да в статике (как у вас в примере) придется долго ждать, или применять демпфирование и динамическую релаксацию

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

Так вы и считаете данную задачу в плоскости. 

Ну и что? формулировка элемента не знает считаю ли я в плоскости или с изгибами

[ДОБРЯК 615]

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

Наверное, поэтому в МКЭ, где сеточная сходимость краеугольный камень, пользуются неявной схемой, чтобы решать задачи без оглядки на критерий шага.

В неявной схеме шаг должен быть. Зависит от метода решения.

Только в неявной схеме шаг по времени принципиально больше. Гораздо быстрее решается одна и та же задача.

[статист 610]

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

Явные методы применяются в своей области - падение, взрыв, пробитие. Посчитайте разрыв контейнера с ВВ при его детонации в неявной постановке

Да в статике (как у вас в примере) придется долго ждать, или применять демпфирование и динамическую релаксацию

Ну то есть подытоживая, для большинства задач Динамики и прочности приобретение Логоса нерентабельно.

[ДОБРЯК 615]

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

Ну и что? формулировка элемента не знает считаю ли я в плоскости или с изгибами

Оболочка - это двухмерный элемент. Стержень - одномерный. Вам как специалисту тех поддержки это нужно знать.

[SergeyGertsik 66]

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

В неявной схеме шаг должен быть. Зависит от метода решения.

Только в неявной схеме шаг по времени принципиально больше. Гораздо быстрее решается одна и та же задача.

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

[ДОБРЯК 615]

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

Ну то есть подытоживая, для большинства задач Динамики и прочности приобретение Логоса нерентабельно.

Естественно. Задачи будут решаться с 10-100 раз дольше. С той же точностью.

 

[SergeyGertsik 66]

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

Ну то есть подытоживая, для большинства задач Динамики и прочности приобретение Логоса нерентабельно.

С чего это? Статика, модальный, динамика есть. Неявный тоже есть, но только для гексов информация у меня о нем двухлетней давности, может уже и подразвили.