Хитрые задачки МКЭ и МДТТ. Вопросы/обсуждения..

[Jesse 960]

07.05.2021 в 13:04, ДОБРЯК сказал:

Как такими моментами инерции можно пренебречь? 

В моём примере я пренебрегал собственными моментами инерции маленького кубика со стороной 0,16мм. Вы мне доказываете, что у кубика 80х80х80 большие моменты инерции и ими нельзя пренебречь. Я вам про одно, Вы мне про ддругое

[ДОБРЯК 605]

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

В моём примере я пренебрегал собственными моментами инерции маленького кубика со стороной 0,16мм. Вы мне доказываете, что у кубика 80х80х80 большие моменты инерции и ими нельзя пренебречь. Я вам про одно, Вы мне про ддругое

Масса кубика 80х80х80 4 кг. Вы пытаетесь заменить большой кубик на маленький, но при этом масса маленького) кубика 4 кг. 

Почему вы сохраняете  массу большого кубика, а моментами инерции большого кубика  пренебрегаете? Можете объяснить?

[Jesse 960]

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

Можете объяснить?

уже объяснил. лень повторяться

[ДОБРЯК 605]

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

Вы мне доказываете, что у кубика 80х80х80 большие моменты инерции и ими нельзя пренебречь.

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

В этом ошибка вашего моделирования.

[ДОБРЯК 605]

@Jesse вот ваша тестовая задача. Размер кубика 80х80х80.

Вы большой кубик моделируете маленьким. И не можете задать моменты инерции маленького кубика как у большого.

А 3 крутильная собственная форма зависит не от массы а от момента инерции большого кубика вокруг оси стержня. 

А вы этот момент инерции обнулили.

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

[ДОБРЯК 605]

Цитата

Моме́нт ине́рции — скалярная физическая величина, мера инертности во вращательном движении вокруг оси, подобно тому, как масса тела является мерой его инертности в поступательном движении. Характеризуется распределением масс в теле: момент инерции равен сумме произведений элементарных масс на квадрат их расстояний до базового множества (точки, прямой или плоскости).

Единица измерения в Международной системе единиц (СИ): кг·м².

Обозначение: I или J.

@Jesse вы обнулили меру инертности во вращательном движении вокруг оси для большого кубика. Поэтому и собственные  частоты неправильные.

 

[ДОБРЯК 605]

Обнулил массу, а моменты инерции оставил. Собственная частота крутильной формы не изменилась. А продольные собственные частоты исчезли. 

[Jesse 960]

кто как такие рёбра жёсткости моделирует в больших листовых сборках? 

Даже для конечных элементов оболочек жадноватенько выходит, учитывая протяжённость этих рёбер жёсткости..
Пока что мысли заменить эти штуки балочными КЭ без массы, то бишь лист сделать тупо плоским, а в местах где рёбра жёсткости прилепить балочку.. ну и потестить отдельно на основную воспринимаемую нагрузку (изгиб, кручение) так, чтобы суммарные перемещения совпали для случая когда лист целиком смоделирован оболочечными КЭ. 

У кого какие мысли? Кто как делает? 
 

[Jesse 960]

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

[karachun 2 038]

@Jesse Если бы были доступны прямоугольные оболочки. то все можно было бы сделать ими.

А так, как вариант - увеличить толщину.

[Jesse 960]

никак руки не доходили о сравнении метода Ланцоша и метода итераций подпространств при определении собственных частот применительно к конкретным задачам. Из очевидного пока что ясно, прямой метод Ланцоша, как правило, более шустрый, но "капризный": более требовательный к сетке, не поддерживает связи балочных и оболочечных КЭ (SW Simulation), реже даёт результат (опция "отсечение частоты" иногда помогает).
Итерационный метод подпространств более надёжный, но ооочень долгий, особенно когда запрашиваешь большое кол-во частот в большой сборке, где частоты получаются сильно сгруппированными друг отн-но друга.

 

То есть вывод пока что простой: исп-ть в лоб метод Ланцоша. Не работает - переходи на итерационный и терпи.
Но быть может в каких то случаях тот или иной способ погрешность сильную даёт? Или итерационный выгоднее, к примеру, в большой сборке на грубой сетке с отн-но искажёнными КЭ, когда запрошено большое кол-во СЧ?
Мб в сети методичка какая найдётся об исп-и данных методов именно в прикладном смысле, на разных задачках/конструкциях, где что выгоднее и т.п. и т.д.? беглый гуглинг плодов не дал
 

[Graf Kim 217]

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

Или итерационный выгоднее, к примеру, в большой сборке на грубой сетке с отн-но искажёнными КЭ, когда запрошено большое кол-во СЧ?
Мб в сети методичка какая найдётся об исп-и данных методов именно в прикладном смысле, на разных задачках/конструкциях, где что выгоднее и т.п. и т.д.? беглый гуглинг плодов не дал

 

Обратись для начала к хелпу больших пакетов.

[Jesse 960]

4 часа назад, Graf Kim сказал:

Обратись для начала к хелпу больших пакетов.

понял, благодарствую.
з.ы.: в Ансисе ж вроде Арнольди был ещё?!)

[ДОБРЯК 605]

16 часов назад, Jesse сказал:

прямой метод Ланцоша

Это большая тема. в одном сообщении всей информации не написать.

Прямой метод Ланцоша - итерационный метод.

Есть 

Прямой метод Ланцоша

Обратный метод Ланцоша

Прямой метод итераций подпространства

Обратный метод итераций подпространства

Метод итераций подпространства с уточнением (например по Якоби)

 

Это итерационные методы.

Но скорость решения зависит от реализации метода. 

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

Если в вашей программе действительно прямой метод Ланцоша, то это медленный алгоритм для определения низших собственных частот.

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

[Jesse 960]

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

Прямой метод Ланцоша - итерационный метод.

Цитата

The Lanczos algorithm is a direct algorithm

https://en.wikipedia.org/wiki/Lanczos_algorithm

[ДОБРЯК 605]

@Jesse почитайте здесь.

Алгоритм Ланцоша

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

[Jesse 960]

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

@Jesse почитайте здесь.

Алгоритм Ланцоша

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

в моём output файле не нашёл итераций. Прямой метод в СВ

[soklakov 1 475]

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

в моём output файле не нашёл итераций.

а попробуйте по той же ссылке, что привели, найти на странице слово iteration

по этой:

https://en.wikipedia.org/wiki/Lanczos_algorithm

[ДОБРЯК 605]

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

в моём output файле не нашёл итераций. Прямой метод в СВ

А для метода итераций подпространства в этом файле есть итерации?

[Jesse 960]

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

А для метода итераций подпространства в этом файле есть итерации?

вот output файл, полученный после решением итераций подпространств.

При решении Методом Ланцоша такого не наблюдается в out файле

@ДОБРЯК @soklakov, хорошо... тогда можно сделать вывод, что по крайней мере для линейной статики/модальника/баклинга итерационные методы в целом всегда быстрее по сравнению с прямыми методами?