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

Нелинейная потеря устойчивости оболочки


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

Доброго Здоровья, глубокоуважаемые форумчане!

Уже не одну неделю бьюсь с нелинейным расчётом потери устойчивости в WB, но всё никак и мозги, чувствую, скоро «поедут» :)

Прошу помощи/совета/пинка ))

Ситуация в следующем.

Необходимо смоделировать потерю устойчивости металлической дымовой трубы и определить критическую силу и критические напряжения с дефектом в виде вогнутости или просто идеальной оболочки. (к примеру Н = 60 м, Д = 2.1 м, т = 8 мм, с промежуточными рёбрами каждые 5м – по необходимости выложу файл проекта). Материал: сталь Ст. 3 Ry = 245 МПа.

Для расчёта, в место ветровой нагрузки, примем смещение верха д.т. – понимаю, что это не корректно для замены ветрового воздействия.

Что делаю:

1. Вариант. Создаю геометрию (идеальную). Прикладываю единичное смещение и делаю расчёт Linear Bukling – далее по фото1 (экспортирую геометрию через UPGEOM 1:1).

post-35507-1367783619_thumb.png

В «Р1» в настройках материала выбираю кривую упрочнения «Bilinear Isotropic» Ry = 245 MPa, Tangent Mod = 76 MPa. (пробовал и другие кривые –рез. тот же.)

Прикладываю смещение 200 мм и настраиваю решатель (фото2)

post-35507-1367783637_thumb.png

на каком-то шаге решение обрывается.

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

И смотрю реакцию в заделке д.т. (фото3),

post-35507-1367783641_thumb.png

надеюсь, что по ней смогу определить критическую силу – место, где сила в опорной реакции начнёт падать (например, см. фото4

post-35507-1367783645_thumb.png

– не знаю как у меня получилось, но это для вымышленной модели, в реальном расчёте использую те же настройки, но результат другой).

2. Вариант. Моделирую вогнутость вручную и выполняю физ./геом. нелинейный расчёт. Критическую силу так же не могу получить…..

Благодарю тех, кто дочитал до конца :) Буду очень признателен за помощь.

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


Помощь:

Второй вариант лучше.

Совет:

Сформулируйте задачу. Какая еще критическая сила? Какие критические напряжения? С силой еще можно понять, но по факту надо выходить на скорость ветра и сравнение с имеющимися ветрами в ветровом регионе. А напряжения тут при чем?

Пинок:

это ж где советуют схему варианта 1 для нелинейной потери устойчивости?

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

Благодарю soklakov :)

Под критической силой понимаю ветровое давление действующее на сооружения для каждого ветрового района, например. И если нормативное 600 Па, то потеря устойчивости будет, например при 1300 Па. (за напряжения только недавно подумал :) ведь если посмотреть формулы СНиП для расчёта местной устойчивости оболочек, то они и варьируют значениями СИГМАкр)

По схеме 1. Пример есть в инете :) Начинал свои исследования со схемы 2, но нужно результата не получил (например, как на фото 4).... поэтому другие варианты реализации начал искать...

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

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

Напомните, пожалуйста, у Вас какая цель? Глобальная/конечная. Дать разрешение на проект, выбрать оптимальный вариант конструкции, спокойно спать по ночам, написать диплом? От того, куда идти, очень сильно зависит, какую дорогу лучше выбрать.

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

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

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

Благодарю. Понимаю, что придётся крутиться))) и никуда от этого не деться...

и ещё. Какую кривую упрочнения Вы посоветуете?

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

Какую кривую упрочнения Вы посоветуете?

С максимально доступным разрешением кривой на первых полпроцента деформации. Если есть - то и до процента-двух хорошо бы. Дальше уже бессмысленно.

Еще есть вариант нагрузить до критического значения, а потом выполнить Explicit расчет. Неплохой вариант, кстати, если красивые картинки интересуют.

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

С трубами не шутят, так что обеспечьте в линейной постановке запас не меньше двойки - тройки и не мучайтесь. :unsure:

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

С трубами не шутят, так что обеспечьте в линейной постановке запас не меньше двойки - тройки и не мучайтесь. :unsure:

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

Вы правы, soklakov, это "не реальное" проектирование :) в двояком смысле этого слова как слитно, так и раздельно ))

провозолся сегодня снова с моделькой) пришёл к тому, что нужно упростить и попробовать "домучать" одну секцию.

Благодаря ребятам - <noindex>http://forum.dwg.ru/showthread.php?t=93015</noindex> построил диаграмму расстяжения стали (Ст3) - файлик во вложении.

post-35507-1367868298_thumb.png

Ст3.rar

Приложил осевую силу 1МН на оболочку с вогнутостью и посчитал Linear Buckling, получил коэф. запаса 33.177.

Умножил первоначальную осевую силу на .5 от коэф. запаса. устанавливаю настройки (с дампингом получается такой же результат)

post-35507-1367868322_thumb.png

Результат

post-35507-1367868380_thumb.png

post-35507-1367868387_thumb.png

Почему эквивалентные напряжения доходят только до 277 МПа, если у меня кривая задана до 480 МПа?

С максимально доступным разрешением кривой на первых полпроцента деформации. Если есть - то и до процента-двух хорошо бы. Дальше уже бессмысленно.

Еще есть вариант нагрузить до критического значения, а потом выполнить Explicit расчет. Неплохой вариант, кстати, если красивые картинки интересуют.

- тоже осень интересно и картинки красивые нужны :) только вот пока у меня критического значения нет)

Очень признателен за советы, т.к. в универе некого спросить :)

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

Приложил осевую силу 1МН на оболочку с вогнутостью и посчитал Linear Buckling, получил коэф. запаса 33.177.

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

с дампингом получается такой же результат

Какой такой же? На что смотреть? Картинок много.

Почему эквивалентные напряжения доходят только до 277 МПа, если у меня кривая задана до 480 МПа?

С чего Вы взяли, что 277? Вот если бы Вы для разных нагрузок, разных геометрий, разных несовершенств получали бы эту цифру, то можно было бы сделать такой вывод. А так... ну потек материал, ну развалилась задача, все нормально. Спросите "что делать"? А что Вам надо? Вы ведь так и не сказали.

- тоже осень интересно и картинки красивые нужны :) только вот пока у меня критического значения нет)

Как же нет, задача ведь не просто так разваливается. Для этой КЭМ максимально выдерживаемое усилие фактически и есть критическое)
Ссылка на сообщение
Поделиться на других сайтах

В нелинейном расчете кроме пошагового статического нагружения, задания пластичности материала (желательно. консервативно - с чем-то типа площадки текучести) надо обязательно включать учет больших деформаций (NLGEOM, ON или Solution / Solution informatiom / Large Deflection/ON). Это было сделано?

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

надо обязательно включать учет больших деформаций (NLGEOM, ON или Solution / Solution informatiom / Large Deflection/ON). Это было сделано?

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

Как раз по этой теме на ютубе туториал есть

<noindex>http://www.youtube.com/watch?v=AGkE-Be3tEs</noindex>

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

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

То что вы называете линейным расчетом на устойчивость дает сущственно меньшую ошибку.

У Бормана нелинейный расчет давал цифру 122, а у меня линейный 48.

И мы долго выясняли почему. :unsure:

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

У Бормана нелинейный расчет давал цифру 122

Не вдаваясь в метод расчета скажу, что для пользователя расчет называется Linear Buckling.
Ссылка на сообщение
Поделиться на других сайтах

Не вдаваясь в метод расчета скажу, что для пользователя расчет называется Linear Buckling.

Называть можно как угодно. :unsure:

Жесткость конструкции нелинейнео зависит от приложенных сил. Это же ведь должен знать пользователь. Если берется объяснять. :unsure:

Главное какой расчет на самом деле зашит в черном ящике, а не как его называть. :unsure:

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

Как раз по этой теме на ютубе туториал есть]]>http://www.youtube.com/watch?v=AGkE-Be3tEs ]]>

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

Чтобы не было путаницы не нужно геометрически нелинейный расчет называть линейным.

Значения точек бифуркации нелинейно зависят от приложенной силы.

Для общего случая естественно. :unsure:

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

Чтобы не было путаницы не нужно геометрически нелинейный расчет называть линейным.

Учитывая контекст, имеете ли Вы в виду, что расчет, который за бугром принято называть Linear Buckling, является расчетом, учитывающим геометрическую нелинейность, то есть НЕ предполагающим малость перемещений?
Ссылка на сообщение
Поделиться на других сайтах

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • Viktor2004
      Уважаемые коллеги Чего-то я туплю страшно Пытаюсь в P-CODE писать программу редактирования барабана инструментов Макропеременная @MDIKEYI #8549 /* MDI Key image должна мне давать номер нажатой кнопки Кнопка F1 под экраном имеет код 249 Кнопка <INPUT> имеет код 152 При нажатии F1 у меня должна включаться D5001.0 а при нажатии <INPUT> выключаться Но в макропеременной #8549 код нажатой клавиши лишь кратковременно мигает. И надо много раз нажимать нужную клавишу что бы она сработала Макропеременная #8549 у меня выводится в диагностический D5000 Кто знает, может #8549 как-то настраивается, что бы программа успела на него отреагировать IF [#8549 EQ 249] THEN #548 = 1 ENDIF IF [#8549 EQ 152] THEN #548 = 0 ENDIF https://www.youtube.com/watch?v=jvsv0YoTy-8   Я как-то пример показывал с нажатием кнопки, но там программа была короткая и она успевала схватывать нажатие кнопки. А тут у меня сначала круги рисует и не успевает
    • malvi.dp
      SW2020 работает. Пожелание - уменьшить размеры окна, а то перекрывает весь чертеж, исправить ошибку в слове детели и, если получится, сделать выбор толщины и типа линии.  
    • sippovich
      Добрый вечер Уважаемые Форумчане,нет ли ни у кого инструкции по восстановлению нулевых точек револьверной головы на станки серии TF и TD серии Fanuc? Заранее благодарен.
    • Heroend
    • ДОБРЯК
      А выкинуть строчку и столбец  это не закрепление? :=)
    • nicomed
      @Kelny  То что,Вы, описали это варианты, а как я понял, @Евдоким Жаркий требуется именно комплекты. Как я понимаю различие: Балка может быть и катанного профиля или гнутого, но она может быть только одна. - это варианты , и, они, не исключают комплекты сменных частей. А, к примеру у заказчика есть три домкрата, и, в процессе эксплуатации зацепа, в разных ситуациях, может использоваться разный тип домкрата. У каждого домкрата свое посадочное место, и для него уже заказчик, будет менять пятки зацепа. Т.е. у него их (пяток) столько, сколько типов домкратов. И это у же комплекты сменных частей. Резюмируя - варианты это еще на уровне изготовителя, а комплекты на уровне заказчика (эксплуатации).       
    • maxx2000
      пусть так это тебе и продолжает казаться, для меня например всё ясно как божий день. Есть вещи которые требуют понимания. Я не ставлю своей задачей выдать набор кнопок, а лишь указать возможный вектор решения задачи. Артас обрабатывает по нормали к элементу(поверхности), а надо обрабатывать вдоль оси Z РП элемента которая должна совпадать с осью Z станка. Он намного ближе к понимаю 5 осевой обработки чем может казаться со стороны.   а зачем его разворачивать? разворачивать надо с конкретной целью. Боком режет, и пусть себе режет.
    • Heroend
      Добрый день! Сел за изучение FS и решил смоделировать настоящую АВО типа БАВО-К. В процессе были некоторые проблемы, но все решалось постепенно. Однако, так и не смог исправить одну ошибку, которая портит весь расчет. Ахиллесовой пятой моего моделирования стала 3-х ходовая теплообменная секция АВО, при входных данных P=3,2 МПа, T газа = 160 С (взятых с работающего оборудования) расчет заканчивается досрочно, с ошибкой "Твердое тело расплавилось". Во время самого расчета система набирает безумное давление (смотрел по дополнительному окну), причем как только я не игрался с расходом, динамика ошибок остается та же.  Может кто-то подсказать из-за чего может быть такая проблема?  
    • Jesse
      скорей всего ДК - это какое-то "общее" требование для всех роторных машин, где есть вал..) вот и копипастят, сами не пойми чего. Ну я уже на 90% уверен, что ДК - это про гироскопику, когда достаточно большой диск кренится вдоль оси и создаёт "гироскопику" - момент, которые гнёт вал и искажает частоты. Наблюдается две ветви - т.н. "вилка". Если эти гироскопич-е эффекты пренебрежимо малы, то это уже не ДК, а просто зав-ть СЧ от оборотов..) Пока так.   да там и от скорости вращения ротора должна зависеть added mass; когда замедляется/ускоряется, везде разная масса получится.  Интересно, а можно из эксперимента определить added mass? ну типа тот же ротор постепенно разгоняешь, и смотришь как меняется основная акустическая мода по времени..
    • Orchestra2603
      @ДОБРЯК: Решил вашу задачку с одним элементом без факторизации...   Частоты:     Формы:  
×
×
  • Создать...