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

[AlexKaz 1 821]

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

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

 

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

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

 

[ДОБРЯК 602]

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

Если сначала сместить

Если сначала сместить, то балка будет деформироваться, будет энергия деформации.

Энергия деформации = работе сил. Если правая часть = 0, то на выходе будут нулевые перемещения.)

Изменено 13 мая 2022 пользователем ДОБРЯК

[Fedor 1 597]

Понятно зачем :) 

[LazyBitch 275]

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

Зачем заранее-то?

ой, а я уже в 14-м мае живу.

[piden 2 874]

Just now, LazyBitch said:

ой, а я уже в 14-м мае живу.

а в году каком?)

[LazyBitch 275]

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

Понятно зачем :) 

фальстарт, простите

[Fedor 1 597]

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

Если сначала сместить, то балка будет деформироваться, будет энергия деформации.

Энергия деформации = работе сил. Если правая часть = 0, то на выходе будут нулевые перемещения.)

Да я одновременно модифицировал матрицу жесткости чтобы учесть заданное перемещение и прикладывал силу. Потом решал .  При нелинейных деформациях без силы происходило как и положено вращение при неизменных размерах. А вот с добавлением сил начинались проблемы. Бало начало девяностых, гайдарномика бушевала. Были проблемы и поважнее. Надо было выжить, экспортом леса занялся чисто конкретно. Так вопрос и завис  :) 

[piden 2 874]

34 minutes ago, Fedor said:

А если работу считать по общему перемещению то она получается очень большой.  Соответственно и внутренняя энергия ... 

А если силу противоположно перемещению направить? То будет ли работа отрицательной? Соответственно и внутренняя энергия ... 

 

@ДОБРЯК, признавайся, ты давеча @Fedora не кусал? Чисто на всякий случай интересуюсь...

[ДОБРЯК 602]

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

Предлагаю проверить на такой задачке на которой споткнулся лет тридцать назад. 

Стальной брусок. С одной стороны шарнир с другой задано перемещение. Синим. Просто такая задача решается и при линейных деформациях и при нелинейных деформациях. Проблемы начались когда прикладывал силы.  Все решал конечно в рамках теории упругости трехмерными  элементами . Интересно как с такими задачками справляется нынешнее ПО :) 

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

Недожали вы решение этой задачи.)

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

А если силу противоположно перемещению направить? То будет ли работа отрицательной?

Если матрица жесткости будет отрицательно определенная, то силу прикладываете в одну сторону, а перемещение в другую. Будет.

[Fedor 1 597]

Цитата

А если силу противоположно перемещению направить? То будет ли работа отрицательной? Соответственно и внутренняя энергия ...

хороший вопрос. Скалярное произведение может быть отрицательным... :) 

 

 

Цитата

Недожали вы решение этой задачи

Дожал бы не спрашивал... :) 

Но квадратичный функционал в мкэ положительно определенный как  обычно пишут. То есть (U, CU)=(U,F) >0  

Ну или не только слегка отрицательный как показывает практика, поэтому обычно используется L D Lt  разложение вместо квадратного корня.  То есть в случае мкэ произведение вектора перемещений на вектор нагрузок неотрицательно должно быть ... 

[SergeyGertsik 64]

5 часов назад, LazyBitch сказал:

Господа, предлагаю посчитать балку переменного сечения, которую предложил Сергей в Логос, ИСПА, Фидесис, АПМ, Лира Софт.

В АПМ будет считать @Jesse,  в Фидесис @Diraimer

я - в Лире. 

Добряк не согласится, он только увиливать горазд (да простит меня админ)

[piden 2 874]

8 minutes ago, SergeyGertsik said:

Добряк не согласится, он только увиливать горазд

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

 

On 5/10/2022 at 5:12 PM, SergeyGertsik said:

Более того, я каждый день буду спрашивать: "Дайте конкретную задачу, где ИСПА лучше Логос"

Пример задачи, где ИСПА справляется гораздо лучше Логос: поставить на уши раздел форума "Динамика и Прочность", генерировать лулзы и срачи более 20 лет...

Хотя в этой области оно вообще вне конкуренции... 

[SergeyGertsik 64]

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

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

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

[Борман 2 375]

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

он создал

Если @Fedor не пропил какие нибудь свои органы, то он сможет вспомнить свою точку зрения 15-ти летней давности на процесс "создания" ИСПА.

Ссылался на труды какой-то алкоконференции, но ссылку так и не дал.

[AlexKaz 1 821]

Речь про ЗИТУ, видимо (м.б. даже про бесплатную сегодня z88 https://en.wikipedia.org/wiki/Z88_FEM_software ). Это за 2012-й. Были и свежее диалоги.

https://cccp3d.ru/topic/50351-неоднородное-поле-напряжений/?do=findComment&comment=453955

https://cccp3d.ru/topic/42282-суперкомпьютер-своими-руками-для-расчетов-в-ansys-и-3d-приложениях/page/46/?tab=comments#comment-514702

https://cccp3d.ru/topic/58966-сейсмика-резервуара-с-частичным-наполнением/page/8/

 

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

 

Где-то я ловил Д. на полном непонимании самой примитивной математики МКЭ, можно поискать тут же на форуме. Чувак потом самостоятельно познал, что есть функции формы и куда они втыкаются.

Был и копипаст Д. из справки на немецком, переведённом гуглом. Тут даже фанатам И. стало интересно, "зачем?" =) В прочем, коммерсам виднее как продавать и что с кодом делать.

 

 

Изменено 13 мая 2022 пользователем AlexKaz

[LazyBitch 275]

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

Речь про ЗИТУ, видимо (м.б. даже про бесплатную сегодня z88 https://en.wikipedia.org/wiki/Z88_FEM_software ). Это за 2012-й. Были и свежее диалоги.

Я в восхищении. Вы всегда так легко вспоминаете события 10-ти летней давности?

[статист 608]

13 часов назад, Fedor сказал:

Интересно не на балках... 

Решается без проблем. С заданием одновременно силы и перемещений.
 

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

 

Линейное.

 

Нелинейное.

 

 

Федор, с Днем рождения!

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

[ДОБРЯК 602]

Борман, спасибо за статьи.

В Логосе я не нашел решателей для решения задач прочности. 

Аэрогидромеханика. Ничего нет для решения плохообусловленных линейной статики.  Что тут сравнивать? Ответы на форуме это детский сад.

Решают совершенно другие задачи. Миллиарды уравнений. 

На 4000 ядрах. 

Изменено 14 мая 2022 пользователем ДОБРЯК

[SergeyGertsik 64]

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

Борман, спасибо за статьи.

В Логосе я не нашел решателей для решения задач прочности. 

Аэрогидромеханика. Ничего нет для решения плохообусловленных линейной статики.  Что тут сравнивать? Ответы на форуме это детский сад.

Решают совершенно другие задачи. Миллиарды уравнений. 

На 4000 ядрах. 

Ну так искать надо лучше. Если задачи решаются, то решатели таки есть, я полагаю. А то что его можно распараллелить на 4000 ядер,  таки плюс. Это вы в ИСПА не можете ни одного решения показать.

[ДОБРЯК 602]

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

А то что его можно распараллелить на 4000 ядер

У вас нет своих решателей. Вы бы хоть сами прочитали статьи.

В статьях решаются задачи аэрогидромеханики. 

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

Это вы в ИСПА не можете ни одного решения показать.

А вы такой же специалист по ИСПА как piden. Вы нашли себе собеседника и на пару с ним доказывайте, что не может решить ИСПА.

Мне этот детский сад не интересен.

12.05.2022 в 18:56, SergeyGertsik сказал:

12.05.2022 в 07:34, ДОБРЯК сказал:

1. Есть ли автоматическое закрепление при решении СЛАУ.

 

Расширить  

1. Нет, да и зачем?

У вас нет своих решателей для решения задач прочности. Поэтому вы так и ответили. И прямых и итерационных решателей у вас своих нет. 

Поэтому у меня и нет интереса, что то сравнивать.

В Ансис серьезные решатели стоят, для локальной версии. Поэтому и сравнивали скорость и точность решения.

 

[ДОБРЯК 602]

 

Присоединяйтесь к pidenу и доказывайте, что ИСПА ничего не может решить.

Вы нашли друг друга на этом форуме.

Ставьте лайки, смейтесь, кривляйтесь как обезьянки. Двоим веселее будет.