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

[AlexKaz 1 822]

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

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

 

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

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

 

[ДОБРЯК 603]

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

Под прогибом вы понимаете Uz?

XYZ - глобальная система координат. UVW - локальная система координат. U V перемещения в плоскости элемента. w - перемещение из плоскости элемента - прогиб.

Будем считать, что элемент плоский. Чтобы не переходить в криволинейную систему координат.)

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

К сожалению, ГОСТа на подобные термины нет

ГОСТа не термины нет, но есть общепринятые термины на русском языке. 

[SergeyGertsik 64]

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

И вы доказывали что такие же элементы кто-то сегодня использует в Ансис, Настран и Абакус. Можете не отвечать. Посмеялись и проехали.)

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

Я уточню ответ на ваш вопрос когда доберусь до рабочего места. Чтобы проверить точность расчетов на изгиб оболочками прошу предоставить мне постановку теста и я покажу результаты расчета.

[KorovnikovAV 509]

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

Распространенней вырожденный.

И это правильный термин.

[ДОБРЯК 603]

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

Распространенней вырожденный.

Дело не только в терминах. В Ансис я не нашел ни одного вырожденного элемента.)

[SergeyGertsik 64]

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

Дело не только в терминах. В Ансис я не нашел ни одного вырожденного элемента.)

SHELL181 Element Description

SHELL181 is suitable for analyzing thin to moderately-thick shell structures. It is a four-node element with six degrees of freedom at each node: translations in the x, y, and z directions, and rotations about the x, y, and z-axes. (If the membrane option is used, the element has translational degrees of freedom only). The degenerate triangular option should only be used as filler elements in mesh generation.

[ДОБРЯК 603]

Найдите здесь дегенеративный КЭ.

https://www.mm.bme.hu/~gyebro/files/ans_help_v182/ans_thry/thy_shp5.html#shp3d3ntrishell

Если вы уверены, что в вашей программе элемент 181 напишите функции формы для этого элемента.

[SergeyGertsik 64]

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

Найдите здесь дегенеративный КЭ.

https://www.mm.bme.hu/~gyebro/files/ans_help_v182/ans_thry/thy_shp5.html#shp3d3ntrishell

Если вы уверены, что в вашей программе элемент 181 напишите функции формы для этого элемента.

Исходя из имеющейся у меня информации, эквивалент реализованного элемента это shell181, keyopt(3)=2

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

[ДОБРЯК 603]

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

Исходя из имеющейся у меня информации, эквивалент реализованного элемента это shell181

Вот функции формы элемента shell181.

Билинейные функции формы. Элемент построенный по теории Тимошенко. Это элемент жесткий на изгиб. Если брать интеграл 2х2. По ребру прогиб меняется по линейному закону. А если вы понижаете жесткость этого элемента, то расскажите как вы это делаете.

 

Но это уже будет химия и к теории МКЭ это уже не имеет отношения.)

30.04.2022 в 16:10, SergeyGertsik сказал:

и оболочечный элемент Голованова (теория Тимошенко)

Если вы говорите, что 4-х узловая оболочка в Логосе это аналог shell181. То почему это элемент Голованова???

Что нового придумал Голованов?)

[SergeyGertsik 64]

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

Вот функции формы элемента shell181.

Билинейные функции формы. Элемент построенный по теории Тимошенко. Это элемент жесткий на изгиб. Если брать интеграл 2х2. По ребру прогиб меняется по линейному закону. А если вы понижаете жесткость этого элемента, то расскажите как вы это делаете.

 

Но это уже будет химия и к теории МКЭ это уже не имеет отношения.)

Если вы говорите, что 4-х узловая оболочка в Логосе это аналог shell181. То почему это элемент Голованова???

Что нового придумал Голованов?)

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

[ДОБРЯК 603]

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

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

В элементе 181 две опции

 

[ДОБРЯК 603]

Консольно закрепленная балка с приложенной на конце силой.

размер балки 1x1x10 см, сила 100 кГ, модуль упругости  Е = 2000000 КГ/см2, коэффициент Пуассона 0.3. Количество элементов по длине 10.

14 часов назад, SergeyGertsik сказал:

эквивалент реализованного элемента это shell181, keyopt(3)=2

Решите это задачу и вы увидите что элемент жесткий если keyopt(3)=2.

[SergeyGertsik 64]

Максимальный прогиб 0,19 см

Даже 0,1994 мм

*0,1994 см

[ДОБРЯК 603]

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

Максимальный прогиб 0,19 см

 

15 часов назад, SergeyGertsik сказал:

эквивалент реализованного элемента это shell181, keyopt(3)=2

У вас получится 0.14. 

Вы кого пытаетесь обмануть?. У меня же исходные коды этого элемента.)

А если решать элементами построенными по теории Кирхгофа, будет 0.198. Без химии. Это чистые элементы по теории МКЭ.

[SergeyGertsik 64]

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

 

У вас получится 0.14. 

Вы кого пытаетесь обмануть?. У меня же исходные коды этого элемента.)

А если решать элементами построенными по теории Кирхгофа, будет 0.198. Без химии. Это чистые элементы по теории МКЭ.

В смысле обмануть? Как есть так и получилось. Вам видео записать с решением? Напишите куда высылать.

[ДОБРЯК 603]

15 часов назад, SergeyGertsik сказал:

Исходя из имеющейся у меня информации, эквивалент реализованного элемента это shell181, keyopt(3)=2

Значит у вас не этот элемент.

 

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

Не являюсь разработчиком

Поэтому можете написать все что угодно. 

Значит у вас интегрирование 1х1. Вот такая химия. Уменьшают жесткость элемента.

Я этот элемент сделал еще в 1987 году. 

[SergeyGertsik 64]

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

Значит у вас не этот элемент.

Как я говорил уточнение смогу сделать как доберусь до рабочего места

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

Значит у вас интегрирование 1х1. Вот такая химия. Уменьшают жесткость элемента.

Я этот элемент сделал еще в 1987 году. 

респект и уважуха )

[ДОБРЯК 603]

Если моделировать 3-х узловыми элементами, 20 элементов, то перемещение = 0.197

@SergeyGertsik решите 3-х узловыми элементами.

 

[SergeyGertsik 64]

[ДОБРЯК 603]

Значит у вас какая-то другая 3-х узловая оболочка. Вы же не знаете по какой теории построена оболочка.

 

Решите задачу, когда силы в плоскости.

[SergeyGertsik 64]