Настройки и математика расчёта на устойчивость

[/prep7 4]

Здравствуйте, форумчане!
Всё ещё занимаюсь потерей устойчивости в APDL. Перед решением использую команду BUCOPT,SUBSP,1,20,,RANGE, что, по моему пониманию, должно означать поиск множителя нагрузки первой формы потери устойчивости в диапазоне (20;inf). 
В силу некоторых причин приходится прогонять огромное количество расчётов на разных вариантах геометрии. Для подавляющего числа вариантов прокатывает замечательно. Но вот иногда встречаются варианты, для которых ансис не может найти форму при данных настройках. Например, есть геометрия, множитель нагрузки для которой равен 106. При задании BUCOPT,SUBSP,1,20,,RANGE решение не находится (нет результатов в файле rst). Если же задать, например, BUCOPT,SUBSP,1,50,,RANGE или BUCOPT,SUBSP,1,70,,RANGE - всё шикарно. Если задавать по умолчанию поиск вокруг нуля - тоже ничего не находит. 

А теперь вопрос: мне непонятно, как может сходимость решения зависеть от заданного диапаона? Или я неправильно понимаю логику поиска критической нагрузки (перебор значений в диапазоне, фактически)?
Заранее спасибо за ответы :)
 

Изменено 17 мая 2019 пользователем /prep7

[AlexKaz 1 821]

Помогло вот это?

Цитата

It is recommended that you request an additional few buckling modes beyond what is needed in order to enhance the accuracy of the final solution. It is also recommended that you input a non zero SHIFT value and a reasonable LDMULTE value (i.e., a smaller LDMULTE that is closer to the last buckling mode of interest) when numerical problems are encountered.

When using the RANGE option, defining a range that spans zero is not recommended. If you are seeking both negative and positive eigenvalues, it is recommended that you use the CENTER option.

 

Изменено 17 мая 2019 пользователем AlexKaz

[/prep7 4]

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

[AlexKaz 1 821]

В LS-DYNA коэффициент может зависеть от начального значения нагрузки и точности решателя. Если взять singlе float и слишком маленькую нагрузку - коэффициент не будет найден. В double float проблем (вроде как) нет.

@/prep7, вкурил, подумал что это начальное значение нагрузки...

 

Изменено 17 мая 2019 пользователем AlexKaz

[/prep7 4]

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

Почему 106 а не 0.106?

потому что прорешал с другим диапазоном и получил 106.. Форма похожа на первую

Нащупал, что проблема связана с выделением памяти для решателя

Нее, не LS Dyna! Я тупал, что не указал, что это банальный линейный анализ устойчивости в чёрном

Выяснил, что проблема появляется, когда не хватает оперативки для in-core memory mode. Буду пытаться объяснить себе, как это работает)

Изменено 17 мая 2019 пользователем /prep7

[Fedor 1 597]

       !______________________Решение____________________________________
       allsel
       n=0 ! число мод для проверки устойчивости, если 0 статический анализ **************
       /SOLU   
           ANTYPE,STATIC             ! STATIC ANALYSIS
           *IF,n,GT,0,THEN  $ PSTRES,ON   $   *ENDIF      ! CALCULATE PRESTRESS EFFECTS
           SOLVE
       FINISH

       *IF,n,GT,0,THEN           ! проверка на устойчивость
         ! для строительных конструкций не ниже 1.36 коэффициент запаса по устойчивости (следует из СНИП)    
           /SOLU
              ANTYPE,BUCKLE            
              BUCOPT,LANB,n             ! USE BLOCK LANCZOS EIGENVALUE EXTRACTION METHOD
              MXPAND,n                 
              SOLVE
           FINISH
       *ENDIF

 

Я так обычно делаю и проблем не припомню...

[/prep7 4]

Плюс-минус километр выяснил. Проблема заключается в выделении памяти для решателя. В статике спаср решатель сам при нехватке оперативы при настройках по умолчанию автоматом отправляется в out-of-ore memory mode. При решении устойчивости сабспейсом такого перехода почему-то не происходит. Помогает насильное управление памятью. В моём случае, оперативы вполне хватает, и я могу себе позволить зарядить DSPOPTION,,INCORE. Теперь всё пашет, оперативы всегда выделяется чуть больше, чем надо для in core решения. От так уот

[sm9h 3]

Решил не плодить новых тем. Имеется вопрос для спецов.

Рассчитываю устойчивость тонкостенного цилиндра под внешним давлением. Элементы shell281, 7 точек интегрирования по толщине; (secoffset, mid).

Упругий расчет. Тип анализа buckle с предварительным расчетом при внешнем давлении 1 Па. С одного конца цилиндр заделан по ux, uy, uz, с другого - по ux, uy. Ось z  совпадает с осью вращения цилиндра.  Проблема в том, что с увеличением длины цилиндра критическое давление постоянно уменьшается. Увеличивал длину от 4 до 32 метров при диаметре оболочки 1,6м. Размер элементов постоянный, меняется только количество элементов.

Может дело в элементах 281?

Кто сталкивался с таким анализом? Отпишитесь, буду признателен.

 

Изменено 10 февраля 2020 пользователем sm9h

[karachun 2 038]

@sm9h Вроде так и должно быть - закрепления на торцах жесткие, значит чем цылиндр длиннее тем он податливее. Сделайте так чтобы торцы могли свободно скользить, тогда критическое давление должно быть постоянным.

[sm9h 3]

Есть значение критической длины, например, L>0,815D*корень(D/(D-c)) в одной из книг. В моем случае правая часть равна 13м, если взять диаметр 1,6м, а толщину 16 мм. Считается, что если длина больше критической, то заделки не должны влиять на критическое давление. Почему-то в расчете этого не наблюдается.

Изменено 10 февраля 2020 пользователем sm9h

[karachun 2 038]

Материал - сталь?

[sm9h 3]

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

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

С одного торца закрепления в осевом направлении не было задано, то есть цилиндр имеет возможность свободно деформироваться вдоль оси.

[karachun 2 038]

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

деформироваться вдоль оси

Но не поперек.

Еще раз.

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

Материал - сталь?

И приведите результаты, которые Вы у себя насчитали для разных длин.

[sm9h 3]

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

Материал - сталь?

Сталь. Модуль упругости 200000Мпа, кэф-т Пуассона 0,3.

 

[karachun 2 038]

@sm9h

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

И приведите результаты, которые Вы у себя насчитали для разных длин.

 

[Jesse 958]

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

Есть значение критической длины, например, L>0,815D*корень(D/(D-c)) в одной из книг

а можно конкретно узнать откуда эта формула?)
быть может в ней загвоздка, что-то оговорено очень ситуативное..

быть может это формула для случая длинной оболочки, когда она теряет устойчивость как стержень?

Изменено 10 февраля 2020 пользователем Jesse

[sm9h 3]

Результаты выложу чуть позже - бегаю по работе)).

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

а можно конкретно узнать откуда эта формула?)
быть может в ней загвоздка, что-то оговорено очень ситуативное..

Формула из книги по конструированию хим. оборудования автора Поникарова И.И. В других источниках формула схожая, может отличаться значением коэффициента. Да, в формуле ошибся - в знаменателе должна быть толщина цилиндра s: L>0,815D*корень(D/(s-c)). 

Изменено 10 февраля 2020 пользователем sm9h
Доп. Информация.

[Jesse 958]

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

Да, в формуле ошибся - в знаменателе должна быть толщина цилиндра s: L>0,815D*корень(D/(s-c))

хз честно сказать.. такая простая задачка, можно сказать даже элементарная. А посчитать в Ансис Меканикал целое дело. В СВ просто 2 кнопки жать и всё.) пару минут делов

[karachun 2 038]

@sm9h

1. Посмотрите на форму потери устойчивости для жесткой заделки (L=32 м).

2. Посмотрите на форму потери устойчивости для скользящей заделки (L=32 м).

3. Делайте выводы сами.

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

в знаменателе должна быть толщина цилиндра

Что это такое?

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

А посчитать в Ансис Меканикал целое дело

Зачем меканикал, есть воркбенч. А еще лучше в Femap.

[sm9h 3]

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

хз честно сказать.. такая простая задачка, можно сказать даже элементарная. А посчитать в Ансис Меканикал целое дело. В СВ просто 2 кнопки жать и всё.) пару минут делов

А результат в вашем случае совпадает с теорией?