Перейти к публикации

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


Рекомендованные сообщения

Здравствуйте, форумчане!
Всё ещё занимаюсь потерей устойчивости в 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 - всё шикарно. Если задавать по умолчанию поиск вокруг нуля - тоже ничего не находит. 

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

Изменено пользователем /prep7
Ссылка на сообщение
Поделиться на других сайтах


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

Цитата

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.

 

Изменено пользователем AlexKaz
Ссылка на сообщение
Поделиться на других сайтах

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

Ссылка на сообщение
Поделиться на других сайтах

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

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

 

Изменено пользователем AlexKaz
Ссылка на сообщение
Поделиться на других сайтах
6 минут назад, AlexKaz сказал:

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

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

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

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

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

Изменено пользователем /prep7
Ссылка на сообщение
Поделиться на других сайтах

       !______________________Решение____________________________________
       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

 

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

Ссылка на сообщение
Поделиться на других сайтах

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

Ссылка на сообщение
Поделиться на других сайтах
  • 8 месяцев спустя...

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

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

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

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

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

 

Изменено пользователем sm9h
Ссылка на сообщение
Поделиться на других сайтах

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

Ссылка на сообщение
Поделиться на других сайтах

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

Изменено пользователем sm9h
Ссылка на сообщение
Поделиться на других сайтах
41 минуту назад, karachun сказал:

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

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

Ссылка на сообщение
Поделиться на других сайтах
2 минуты назад, sm9h сказал:

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

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

Еще раз.

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

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

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

Ссылка на сообщение
Поделиться на других сайтах

@sm9h

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

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

 

Ссылка на сообщение
Поделиться на других сайтах
15 минут назад, sm9h сказал:

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

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

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

Изменено пользователем Jesse
Ссылка на сообщение
Поделиться на других сайтах

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

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

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

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

Изменено пользователем sm9h
Доп. Информация.
Ссылка на сообщение
Поделиться на других сайтах
10 минут назад, sm9h сказал:

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

:smile:

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

Ссылка на сообщение
Поделиться на других сайтах

@sm9h

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

Fixed.thumb.png.935d5f97df6e4d7b20ff401e64f76d93.png

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

Sliding.thumb.png.b27365b7463b6a4601aea2086c56efc9.png

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

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

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

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

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

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

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

Ссылка на сообщение
Поделиться на других сайтах
2 минуты назад, Jesse сказал:

:smile:

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

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

Ссылка на сообщение
Поделиться на других сайтах

Присоединяйтесь к обсуждению

Вы можете опубликовать сообщение сейчас, а зарегистрироваться позже. Если у вас есть аккаунт, войдите в него для написания от своего имени.
Примечание: вашему сообщению потребуется утверждение модератора, прежде чем оно станет доступным.

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу.




  • Сообщения

    • gudstartup
      а вы хоть станок проверяли по программе на изделии на точность прежде чем товарищей этих выгнать? если нет то грешите на самих себя! система в наших краях еще не распространенная поэтому и тем тут нет надо в поднебесную писать
    • AlexArt
      Ну допустим, ты и на другом ресурсе это опубликовал. А не коммуниздил. Но вот продвигать воровство от государства, ворующее из Вики, это верх мерзости.
    • maxx2000
      Ах, да. Фильтры выбора добавили. Теперь можно выбрать только то что видно на первом плане, а не вместе с тем что с обратной стороны детали. В общем надо обновляться. Как раз работёнка на прессформу нарисовалась 
    • maxx2000
      Причина того - Кроилово. Кроилово всегда приводит к попадалову. Месяц простоял сколько мильонов деревянных потеряли? Вопрос риторический. И ещё будет стоять. Как памятник человеческой глупости и жадности.
    • AlexKaz
      "9 июля 1968 года на мышах был проведен самый знаменитый эксперимент американского ученого-этолога Джона Кэлхуна «Вселенная-25». Суть опыта заключалась в создании идеальных условий, где мыши могли бы жить и размножаться, не ведая никаких забот, вдали от хищников и в отсутствие эпидемий и заболеваний. Для этих целей ученый построил специальный загон, куда были помещены четыре пары белых мышей (самцов и самок). В распоряжении мышей всегда была чистая вода и еда в изобилии, специальные гнезда, где можно обустроить себе жилище ― гнезд в загоне хватало для проживания нескольких тысяч мышей. Температура в загоне в среднем составляла около 20 ℃ и была комфортной для мышей. Животные не подвергались никаким влияниям извне и жили в идеальных условиях в свое удовольствие. А дальше началось самое интересное. На первом этапе эксперимента мыши хорошо размножались, вели активный образ жизни, охотно играли. На следующей фазе эксперимента мыши стали есть меньше, перестали наедаться до отвала. На третьей фазе эксперимента, когда в загоне были уже сотни мышей, произошло распределение социальных ролей, стала ярко выраженной иерархия, клановость. Появились так называемые отверженные ― молодые особи, которых другие, взрослые мыши сгоняли в центр загона, не давали им вести нормальный образ жизни, причиняли физический вред. В природе такое, наверное, было бы невозможно, ведь эти мыши-агрессоры просто не дожили бы до старости: их бы съели хищники. Но в загоне Кэлхуна хищников не было, и взрослые мыши начали попросту издеваться над молодняком. Образовались две большие группировки: самцы-одиночки и самки-одиночки. При этом самки-одиночки отказывались спариваться <с менее статусными многочисленными молодыми самцами и с оставшимися старыми статусными> и отвергали ухаживания самцов. У мышей стал проявляться тотальный индивидуализм, мыши не стремились создать семью. На последней, четвертой стадии мышиная популяция стала сокращаться. Появились самцы, которых сам Кэлхун назвал «красивыми» (англ. beautiful ones), из-за отсутствия ран и рубцов. <В оригинале: They never engaged in sexual approaches toward females, and they never engaged in fighting, and so they had no wound or scar tissue. Thus their pelage remained in excellent condition. - Дословный перевод: Они никогда не прибегали к сексуальным подходам к самкам, и они никогда не участвовали в боях, и поэтому у них не было ран или рубцовой ткани. Таким образом, их шерсть сохранилась в отличном состоянии.> Эти мыши не вступали в борьбу за самок и территорию, не проявляли активности к размножению и только питались, спали и чистили шёрстку. У мышей стали проявляться различные формы девиантного поведения, вспышки агрессии. Самки стали проявлять агрессию, защищать себя сами, стали умерщвлять своих детенышей, а затем окончательно отказались размножаться. На пике эксперимента в загоне одновременно проживало чуть более двух тыс. мышей. Еды и гнезд было достаточно для дальнейшего роста популяции, но через четыре года после начала эксперимента Кэлхун остановил свой опыт, потому что в загоне осталось чуть более сотни мышей, и все они уже вышли из репродуктивного возраста. По итогам эксперимента Кэлхун пришел к выводу, что достижение определенной плотности населения и заполнение социальных ролей в популяции приводит к распаду общества" https://physicsoflife.pl/dict/pic/calhoun/calhoun.. https://scientificrussia.ru/articles/utopiya-dlya-mys.. https://ru.wikipedia.org/wiki/Кэлхун,_Джон_(этолог)
    • gudstartup
      @Koels вот в чем дело пока ds609 это предупреждение поэтому F может и не появится если sv601 это значит ошибка. возможно при нагреве радиатора серво определяет это как предупреждение или ваш вентилятор крутиться медленнее чем оригинальный и серва думает что он встал хотяпри этом обычно на экране в строке состояния FAN.мигает больше у меня вариантов нет....  
    • ДОБРЯК
      Решите любым алгоритмом. Тогда будет конструктивный разговор. :=)
    • Fedor
      https://en.wikipedia.org/wiki/List_of_numerical_analysis_topics#Eigenvalue_algorithms     :) 
    • Юрий К.Ф.
      Добрый день. Не нашёл тут тему по стойке Китайско Китайской)) Lynuc N3ME. Видать мне так повезло с её наличием)) Приобрели 5-ти осевой Китаец. В б/у состоянии после удара по оси Z. Отремонтировали по механике, заменили батареи на драйверах, выставили лимиты. Всё Ок. Пригласили со стороны людей которые бы разобрались по операторской части. Те два выходных ковырялись, после сказали покажут расскажут, но за огромные деньги. Не сошлись. После месяц станок простоял, когда включили перестал реагировать на регулировку скорости шпинделя. То есть в режиме Jog, включаем обороты, которые стандартно 2140-2149 об/мин. При регулировке процетности не меняются (сама процентность показывает на мониторе). Так же при включении оборотов через команду M03S300 или другое значение, скорость так же показывает 2140-2149 об/мин. Грешить на тех товарищей с которыми не сошлись по деньгам для обучения, как то не хочется. Поковырялся в настройках шпинделя, вроде всё в норме. Проводку на шпинделе прозвонил, целая. В чём причина, не понятна. Кто нибудь сталкивался с подобным, или с подобной стойкой? Может подсказать варианты причины подобного?
    • ДОБРЯК
×
×
  • Создать...