Расчет гофрированной мембраны на собственные частоты

[qeswer 4]

Всем доброго времени. Определяю собственные частоты для гофрированной мембраны. Моделирую с помощью Shell элементов и пяту с объемными элементами, кроме того пробовал и полностью с объемными. Результат получается, но моды получаются не физичными, можете посоветовать как можно устранить данную проблему. Также когда рассчитываю мембрану с пластиной, происходит внедрение тел в друг друга. хоть контакты и заданы. Как можно решить это?

[статист 610]

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

вадратная оболочка со стороной квдарата 10 мм. толщина 1 мм.

материал сталь, 2е11, 0.3

три угла закреплены по трансляционным степеням свободы, так называемое шарнирное опирание.

в четвертом углу приложен момент в плоскости обоболочки (вокруг нормали к оболочке).

величина момента 10 н*мм.

угол поворота в узле прикладывания момента: 0,85 град (надеюсь, тут нет косяка с размерностями. чет мне припоминается, какие-то проблемы мне когда-то попадались. они были решаемые, но решать их сейчас не хочется, я ведь не помню, о чем речь была)

 

силы реакций в трех опорах на картинке в правом нижнем углу.

 

только, пожалуйста, не просите энергию деформации)

Сделал тот же тест в Code_Aster.

 

И рассмотрел три варианта:

1) когда эта шестая степень свободы не задействована

COEF_RIGI_DRZ = 1e-3

2) когда фиктивная жесткость задействована, но малая

COEF_RIGI_DRZ = - 1e-3

3) когда она задействована, но большая

COEF_RIGI_DRZ = -1.0

 

 

Получил следующие результаты:

1) Реакции в закрепленным узлам по нулям. Здесь привожу реакции в одном узле

 DX           DY           DZ          
0.00000E+00  0.00000E+00  0.00000E+00

перемещения также по нулям:

https://imageup.ru/img209/4452724/disp_rigi_1.png.html

 

2) Появились реакции

 DX           DY           DZ          
6.48418E-01  3.84107E+00  0.00000E+00

максимум общих перемещений 8,5e-5

https://imageup.ru/img143/4452726/disp_rigi_2.jpg.html

 

3) Реакции изменились

DX           DY           DZ          
3.90169E-01  1.52623E+00  0.00000E+00

максимум перемещений теперь стал равен 2,9e-5

https://imageup.ru/img140/4452731/disp_rigi_3.jpg.html

 

 

В общем явно всё зависит от этой непонятной жесткости. Но как явно я проследить из этого эксперимента не могу.

 

@soklakov, может поиграетесь с параметром SECCONTROL - DRLSTIF в shell281.

https://www.mm.bme.hu/~gyebro/files/ans_help_v182/ans_cmd/Hlp_C_SECCONTROL.html

Может быть его можно и выключить в ансисе.

Изменено 30 июля 2023 пользователем статист

[ДОБРЯК 615]

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

когда уже?

Я уже показал что результаты отличаются в 56 раз. 

[soklakov 1 485]

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

может поиграетесь с параметром SECCONTROL - DRLSTIF в shell281.

зачем?

[ДОБРЯК 615]

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

В общем явно всё зависит от этой непонятной жесткости.

Естественно будет зависеть от фиктивной жесткости которую вы вводите. 

[soklakov 1 485]

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

В общем явно всё зависит от этой непонятной жесткости.

она бесконечно большая. все просто.

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

Я уже показал что результаты отличаются в 56 раз. 

как я могу из этого сделать вывод, что момент не момент?

особенно с учетом того, что я вашим результатам не доверяю?

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

Может быть его можно и выключить в ансисе.

кого выключить? связь этой степени свооды с остальными? а зачем?

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

этой непонятной жесткости

а почему она непонятная? вы же ссылки скидывали про нее

[ДОБРЯК 615]

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

она бесконечно большая. все просто.

Бесконечно большая жесткость в КЭ. Как вы потом СЛАУ будете решать? :=)

[soklakov 1 485]

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

Бесконечно большая жесткость в КЭ. Как вы потом СЛАУ будете решать? :=)

правильн понимаете) она будет какая-то конечная

[статист 610]

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

Естественно будет зависеть от фиктивной жесткости которую вы вводите. 

Ещё ее можно и выключить. Не знаю только, можно выключить ли в ансисе.

 

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

зачем?

Чтобы убедиться, что перемещения зависят от него.

DRLSTIF = Drill stiffness scale factor.

[soklakov 1 485]

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

Чтобы убедиться, что перемещения зависят от него.

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

[статист 610]

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

а почему она непонятная? вы же ссылки скидывали про нее

Потому что я ее задаю, она не высчитывается из материала и геометрии элемента.

 

 

[soklakov 1 485]

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

Потому что я ее задаю, она не высчитывается из материала и геометрии элемента.

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

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

высчитывается из материала и геометрии элемента.

думаю, дефолтное значение на них опирается. на материал так точно

[статист 610]

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

бесконечное число туда не запихнуть.

Не запихнуть =) Решение разваливается.

Я пока не очень понимаю цели этого теста. Физичность этой жесткости я прощупать не могу =)

[soklakov 1 485]

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

Не запихнуть =) Решение разваливается.

невязки страдают на скачках. можно тренд построить в пределах того, где решается. 

@ДОБРЯК а у вас элементов сколько в модели было?

[статист 610]

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

можно тренд построить в пределах того, где решается. 

Да пока не очень хочется. Не вижу конечной цели исследования. Мне этот эксперимент не понятен. Выжал из него, что есть. Лучше подожду, что @ДОБРЯК скажет.

[soklakov 1 485]

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

Я уже показал что результаты отличаются в 56 раз. 

я вот элементы на второй порядок переставил, получил поворот 0,048 рад.

результаты снова отличаются. ну и что?

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

Не вижу конечной цели исследования.

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

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

Выжал из него, что есть.

вы задали два маленьких значения там, где нужны большие.

[статист 610]

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

вы задали два маленьких значения там, где нужны большие.

В документации по Code_Aster рекомендуемое значение этого параметра 1e-5, так что не думаю, что я задал маленькую величину

[soklakov 1 485]

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

не думаю, что я задал маленькую величину

потом расскажете)

[ДОБРЯК 615]

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

 

@ДОБРЯК а у вас элементов сколько в модели было?

2502

[soklakov 1 485]

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

Лучше подожду, что @ДОБРЯК скажет.

он сказал номер элемента)

тьфу... количество

[статист 610]

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

потом расскажете)

Лучше повторите это в ансисе и посмотрите

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

он сказал номер элемента)

тьфу... количество

Ну пока он ничего не сказал о том, зачем ему это было нужно