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

[AlexKaz 1 850]

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

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

 

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

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

 

[SergeyGertsik 64]

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

@SergeyGertsik на сайте сказано, что

и

В связи с этим возникает вопрос. Какое геометрическое ядро используется в вашей программе?

Геом. ядро Аскона, как в Компасе

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

Если вы не разработчик. Зачем придумываете свои термины?)

 

Какие изгибные заклинивания) в оболочке Тимошенко, с которыми Голованов успешно борется?

А в стержне Тимошенко с такими же изгибными заклиниваниями боретесь?)

 

Какие термины я придумал? Заклинивание? Это в отчетах разработчиков термин в других пакетах это называется locking.

 

Это не костыли, абсолютно те же алгоритмы применяются и в настране и в ансисе и т.д.

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

elibrary не выдаёт статьи без регистрации.

Если статьи были выпущены, скорее всего, они есть в свободном доступе на сайте СЯЦ.

http://book.sarov.ru/wp-content/uploads/2021/05/Works-RFNC-VNIIEF-v23-1-11.pdf

Решатель явный со всеми вытекающими типа hourglass instability.

Без регистрации статьи журнала ВАНТ можно посмотреть на сайте http://book.sarov.ru/

Оболочка Голованова применяется для статического и вибрационного анализов. Для явной динамики реализованы оболочечные элементы BLT и Хьюса-Лю

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

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

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

Как раз полноинтегрируемые элементы в явной динамике подвержены объемному заклиниванию.

[ДОБРЯК 611]

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

Какие термины я придумал? Заклинивание? Это в отчетах разработчиков термин в других пакетах это называется locking.

В теории МКЭ нет такого термина. И в теории Тимошенко нет такого термина. 

И в документации Ансис а разделе где описывается теория построения КЭ нет такого термина.

Коммерческие пакеты могут вводить любые термины. Вы понимаете что означает данный термин для 3-х узловой оболочки или 4-х узловой оболочки.

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

Как раз полноинтегрируемые элементы в явной динамике подвержены объемному заклиниванию.

А в линейной статике их не заклинивает?)

На формулах покажите что такое объемное заклинивание 3-х узловой оболочки или 4-х узловой. И расскажите как разработчики вашей программы борются с объемным заклиниванием.)

Поймите, что явная или неявная динамика не имеет никого отношения к методу конечных элементов.

Вы уходите от теории построения 3- и 4-х узловых тонких оболочек в вашей программе. Хоть функции формы (базисные функции) можете написать для этих элементов и объяснить что означает термин полное интегрирование КЭ для этих функций формы.?

[SergeyGertsik 64]

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

В теории МКЭ нет такого термина. И в теории Тимошенко нет такого термина. 

И в документации Ансис а разделе где описывается теория построения КЭ нет такого термина.

Коммерческие пакеты могут вводить любые термины. Вы понимаете что означает данный термин для 3-х узловой оболочки или 4-х узловой оболочки.

Не хочу затевать споры и тем более искать для кого-то информацию. Но если вы посмотрите в книгу Bathe K.J. Finite Element Procedures. New Jersey: Prentice-Hall, Inc. 1996 на термины в конце книги со ссылка ми на страницы, то увидите что locking (заклинивание или запирание в наших источниках) имеет место быть в теории МКЭ. И руководстве ансиса тоже есть про заклинивания и методы борьбы с ними как для оболочек, так идля объемных элементов В теории Тимошенко его естественно нет, потому что он относится к МКЭ. Для оболочек это означает излишнюю жесткость на изгиб или мембранные деформации.

[ДОБРЯК 611]

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

И руководстве ансиса тоже есть про заклинивания и методы борьбы с ними как для оболочек, так идля объемных элементов

Дайте эту ссылку, где сказано, что в Ансис заклинивает) все КЭ оболочки и все объемные КЭ.

Объясните своими словами, что означает термин заклинивает. Я сам термин не понимаю. Я сделал КЭ модель из тонких оболочек. Составил матрицу жесткости [K]. Решил СЛАУ. {F} = [K]{V}. Получил перемещения {V}. В моей программе ничего не заклинивает, почему с ваших слов в Ансис что-то заклинивает, если КЭ в моей программе и в Ансис один в один?

[SergeyGertsik 64]

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

Дайте эту ссылку, где сказано, что в Ансис заклинивает) все КЭ оболочки и все объемные КЭ.

Объясните своими словами, что означает термин заклинивает. Я сам термин не понимаю. Я сделал КЭ модель из тонких оболочек. Составил матрицу жесткости [K]. Решил СЛАУ. {F} = [K]{V}. Получил перемещения {V}. В моей программе ничего не заклинивает, почему с ваших слов в Ансис что-то заклинивает, если КЭ в моей программе и в Ансис один в один?

Заклинивание это излишняя жесткость. В ансисе не заклинивает и  в Логосе не заклинивает, именно потому что и там и там реализованы специальные методы подавления заклиниваний.

[ДОБРЯК 611]

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

Заклинивание это излишняя жесткость. В ансисе не заклинивает и  в Логосе не заклинивает, именно потому что и там и там реализованы специальные методы подавления заклиниваний.

Теперь понятно, что заклинивание - это излишняя жесткость в ваших терминах. 

Какие специальные методы подавления заклиниваний реализованы в Логосе?

Только не рассказывайте про то что реализовано в Ансис. Пока вы не можете пояснить, что сделано в вашей программе но как-бы знаете, что сделано в Ансис.)

Расскажите про то, что сделано в Логосе, например для тетраэдра, или для 3-х узловой оболочки.  Как реализованы специальные методы подавления заклиниваний, например для этих элементов?

 

[SergeyGertsik 64]

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

Теперь понятно, что заклинивание - это излишняя жесткость в ваших терминах. 

Какие специальные методы подавления заклиниваний реализованы в Логосе?

Только не рассказывайте про то что реализовано в Ансис. Пока вы не можете пояснить, что сделано в вашей программе но как-бы знаете, что сделано в Ансис.)

Расскажите про то, что сделано в Логосе, например для тетраэдра, или для 3-х узловой оболочки.  Как реализованы специальные методы подавления заклиниваний, например для этих элементов?

 

Реализованы ANS, EAS, b-bar, MITC. Как бы то, что сделано в ансис, настран и т.д. тоже самое и реализуется в Логос. В этом и есть импортозамещение, реализовывать требуемый предприятиям функционал популярных расчетных пакетов по их запросам.

Для тетраэдров сразу не могу сказать что реализовано без обращения к документации, для 8 узловых гексаэдров применяются EAS (в т.ч. упрощенный) и b-bar. Для 3-х узловой оболочки применяется  MITC.

[AlexKaz 1 850]

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

Какие специальные методы подавления заклиниваний

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

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

[ДОБРЯК 611]

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

Реализованы ANS, EAS, b-bar, MITC. Как бы то, что сделано в ансис

Опять какие-то буковки. Что означают эти буковки? 

И опять вы прикрываетесь Ансис.)

Вот элемент который есть в Ансис. 

3-х узловая оболочка. Прогиб w по каждому из ребер меняется по кубическому закону. И в этом элементе нет этих буковок ANS, EAS, b-bar, MITC. 

И в этом элементе нет никакого заклинивания.)) Это совместная 3-х узловая оболочка.

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

В этом и есть импортозамещение

Именно такие КЭ надо импортозамещать. А вы придумываете какую-то химию в своих КЭ в Логосе и говорите, что это импортозамещание.)

Вы вместо Мерседеса предлагаете Запорожец. Это аналогия понятная всем.)

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

[SergeyGertsik 64]

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

Опять какие-то буковки. Что означают эти буковки? 

И опять вы прикрываетесь Ансис.)

Вот элемент который есть в Ансис. 

3-х узловая оболочка. Прогиб w по каждому из ребер меняется по кубическому закону. И в этом элементе нет этих буковок ANS, EAS, b-bar, MITC. 

И в этом элементе нет никакого заклинивания.)) Это совместная 3-х узловая оболочка.

Именно такие КЭ надо импортозамещать. А вы придумываете какую-то химию в своих КЭ в Логосе и говорите, что это импортозамещание.)

Вы вместо Мерседеса предлагаете Запорожец. Это аналогия понятная всем.)

Я не могу понять что вы от меня хотите. Расписывать всю математику Логоса? Я этого не буду делать.

Ну и что вы мне показали? Функции формы в L координатах и что? В Логосе тоже самое, и в Логосе и в Ансисе борются с заклиниваниями, читайте документацию.

Чтобы Понять что это за буковки вот вам ссылка на пару статей:

Andelfinger, U. & Ramm, E. (1993). EAS-elements for two-dimensional, three-dimensional, plate and shell structures and their equivalence to HR-elements. International Journal for Numerical Methods in Engineering. 36, 1311-1337. ;

Pramin Norachan, Songsak Suthasupradit, Ki-Du Kim. A co-rotational 8-node degenerated thin-walled element with assumed natural strain and enhanced assumed strain // Finite Elements in Analysis. No. 50. 2012. pp. 70-85.

Плюс книга которую я привел.

При чем тут прикрываться ансисом я вообще не понял.

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

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

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

Во первых, речь идет про статику, а не явную динамику, во вторых эти алгоритмы применяются во ВСЕХ пакетах численного моделирования, это не Логосу костыли придумали, это принятые алгоритмы ВСЕХ пакетов. Советую почитать литературу по МКЭ потяжелее Зинкевича (при всем уважении к этому автору).

[ДОБРЯК 611]

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

Ну и что вы мне показали? Функции формы в L координатах и что? В Логосе тоже самое, и в Логосе и в Ансисе борются с заклиниваниями, читайте документацию.

Функции формы в L координатах для этого элемента известны только для перемещений в плоскости элемента. А для прогиба их не знает никто. Читайте еще раз.)

В этом элементе нет никакого заклинивания. От слова совсем нет.)

Читайте еще раз. 

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

Прогиб w по каждому из ребер меняется по кубическому закону.

Вы не слышите то что вам объясняют.)

[SergeyGertsik 64]

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

Функции формы в L координатах для этого элемента известны только для перемещений в плоскости элемента. А для прогиба их не знает никто. Читайте еще раз.)

В этом элементе нет никакого заклинивания. От слова совсем нет.)

Читайте еще раз. 

Вы не слышите то что вам объясняют.)

Как вы хотели увидеть в апроксимации поля перемещений заклинивание?

[ДОБРЯК 611]

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

Во первых, речь идет про статику, а не явную динамику, во вторых эти алгоритмы применяются во ВСЕХ пакетах численного моделирования, это не Логосу костыли придумали, это принятые алгоритмы ВСЕХ пакетов.

Вот вы опять про все пакеты. Это смешно уже читать. Вы не можете объяснить какая  3-х узловая оболочка в вашей программе. И как вы ее делаете мягче, но постоянно говорите про все пакеты.)

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

Как вы хотели увидеть в апроксимации поля перемещений заклинивание?

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

То что этот элемент жесткий на изгиб я написал уже давно.)

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

Если это элементы Тимошенко 1-го порядка аппроксимации, то они очень жесткие на изгиб.

Давно уже написал, что это элементы жесткие на изгиб. Как вы в своей программе с этим боретесь для 3-х узлового элемента. Функции формы я вам написал.)

[SergeyGertsik 64]

Аналог трехузловой оболочки Логос ANSYS: SHELL181, Full with Incompatible , NASTRAN: CTRIA3, ABAQUS Standard: S3 алгоритм MITC. Как я уже говорил, математику я расписывать не собираюсь.

 

Алгоритм MITC используется для заклиниваний опять же говорю. То что этот элемент жесткий также как тетраэдр первого порядка все знают, и я не говорил, что в Логосе мягче чем в ансисе. 4-х узловые элементы Голованова показывают нормальные результаты. В треугольниках первого порядка оболочки крайне не рекомендуется их вообще применяют только для построения сетки как переходные места. Вы это имели ввиду? Так сразу и надо было говорить что  "элементы Тимошенко 1-го порядка аппроксимации, то они очень жесткие на изгиб" имелось ввиду для треугольников.

[ДОБРЯК 611]

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

В треугольниках первого порядка оболочки крайне не рекомендуется их вообще применяют только для построения сетки как переходные места.

Откуда такая уверенность. Откуда такая информация?

Читайте еще раз.

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

Вот элемент который есть в Ансис. 

3-х узловая оболочка. Прогиб w по каждому из ребер меняется по кубическому закону.

Вы не читаете, то что вам объясняют?  Прогиб w по каждому из ребер меняется по кубическому закону в этом элементе. Это что элемент первого порядка???

В вашей программе в 3-х узловой оболочке прогиб w по какому закону меняется? 

Только не рассказывайте мне про Ансис, Настран и Абакус. Очень вас прошу.

 

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

Так сразу и надо было говорить что  "элементы Тимошенко 1-го порядка аппроксимации, то они очень жесткие на изгиб" имелось ввиду для треугольников.

4-х узловые элементы Тимошенко, то же жесткие на изгиб. Если я буду говорить про 4-х узловые элементы вы совсем запутаетесь.)

Перемещения u и v к теории Тимошенко не имеют никакого отношения. Только прогиб w  и углы поворота. По какому закону меняется прогиб w в ваших оболочках Тимошенко? 

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

[SergeyGertsik 64]

У нас под 3-х узловой оболочкой подразумевается дегенеративный 4-х угольник поэтому сетки должны быть преимущественно четырехугольные.

[ДОБРЯК 611]

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

У нас под 3-х узловой оболочкой подразумевается дегенеративный 4-х угольник поэтому сетки должны быть преимущественно четырехугольные.

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

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

[SergeyGertsik 64]

В ансис и т.д. тоже есть вырожденные треугольники. Под прогибом вы понимаете Uz?

[KorovnikovAV 510]

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

дегенеративный

в термине уверены?

[SergeyGertsik 64]

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

в термине уверены?

Нет. Распространенней вырожденный. Но и под вырожденным можно понимать как представление оболочки как вырожденного трехмерного тела. К сожалению, ГОСТа на подобные термины нет