Femap: Расчетная длина колонны переменного сечения

02/2022

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

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

Сегодня подробно расскажем о потере устойчивости составной колонны.

Что пытаемся решить?

Вполне разумно начинать с самого начала! Прежде всего нужно признать, что придется затронуть некоторые “общие” вещи, чтобы разобрать данный пример. Но этот пример, безусловно, является базовым элементом нашего сегодняшнего обзора.

Итак, вопрос - какова расчетная длина колонны на Рисунке 1:

Рисунок. Расчетная длина колонны переменного сечения
Рисунок 1

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

Вопрос - какова расчетная длина подобной колонны, не является “случайным”. Если спроектировать подобную колонну, то КЭ программное обеспечение (с некоторыми принятыми правилами для программного кода) попросит указать расчетную длину, каждого проектируемого компонента и это наиболее важный этап.

Вопросы

Если колонна имеет переменное поперечное сечение, то у нее не будет расчетной длины. По крайней мере, не в классическом смысле этого слова. Это, конечно, не означает, что такую колонну нельзя спроектировать. Это вполне выполнимая задача, просто следует знать, как это сделать. Лучше всего – начать с простого. Следует предельно четко представлять, что и почему требуется для разработки изделий с учетом прочности. И вот с чего можно начать.

Вам не понадобится расчетная длина

Все, что нужно, это знать гибкость компонента конструкции. В настоящее время наиболее распространенный способ проектирования изделий делится на два этапа. Сначала рассчитывается допустимая пластическая нагрузка для компонента конструкции (то есть, когда изделие разрушается… буквально). Затем для этой допустимой нагрузки применяется понижающий коэффициент. С учетом прочности (потери устойчивости и т.д.), этот понижающий коэффициент вычисляется в зависимости от гибкости компонента конструкции.

Для чего нужна расчетная длина?

Использование расчетной длины – один из самых простых способов вычисления критической силы по уравнению Эйлера. Когда она известна, то и легко рассчитать гибкость. Вот почему программное обеспечение (САПР) просит указать расчетную длину!

Другими словами, чтобы спроектировать колонну, нужно знать ее гибкость. И используемое программное обеспечение потребует от нас “ввести” расчетную длину. Давайте остановимся на ней подробнее, чтобы понять принцип работы.

Как устроен процесс проектирования

Для ответа на вопрос о расчетной длине необходим более широкий контекст. Вот почему стоит сделать акцент на том, как устроен процесс проектирования. Пользователи программой, даже те, кто пользуется длительное время, могут не всегда понимать некоторые уравнения, поскольку программное обеспечение делает это сама “в фоновом режиме”. Но важно понимать, что там происходит.

Вот почему стоит рассказать о стандартном подходе к проектированию сжатой части конструкции.

Стандартный подход (основанный на EN 1993-1-1) к проектированию “простой” сжатой колонны выглядит примерно так:

  • Рассчитайте пластическую нагрузку: вы знаете, какую нагрузку можете приложить, прежде чем объект “выйдет из строя”, ввиду того, что вы просто раздавили его. Как можно видеть ниже, сделать это чрезвычайно просто. Такой “тип” предельной нагрузки обычно называют предельной нагрузкой “поперечного сечения”. Проще говоря, это величина нагрузки, которую потребуется приложить к “вырезанному фрагменту” колонны, чтобы раздавить его, поэтому очевидно, что устойчивость не является частью проблемы (вырезанный фрагмент поперечного сечения под действием усилия гидравлического пресса не потеряет устойчивость). Если бы наше поперечное сечение относилось к классу 4 (стандарта EN 1993-1-1), то даже в таком случае была бы некоторая потеря устойчивости. Предположим, что мы используем “нормальное” горячекатаное поперечное сечение, поэтому класс 4 доставит проблем.
Рисунок. Площадь поперечного сечения. Предел текучести
  • Вышеизложенного недостаточно. Было бы очень хорошо, если бы мы могли просто закончить на этом, но устойчивость нужно учитывать. Проще говоря, мы даже не будем рассчитывать предельную нагрузку потери устойчивости, мы просто уменьшим вышеуказанную предельную нагрузку поперечного сечения на “коэффициент снижения допускаемых напряжений”, как это можно увидеть ниже. Стоит обратить внимание, что коэффициент запаса прочности изменился с “типа” M0 на M1, что означает лишь то, что мы имеем дело с задачей устойчивости, оба значения коэффициента фактически равны 1.0
Рисунок. Коэффициент снижения допускаемых напряжений
  • Давайте разберемся с коэффициентом снижения допускаемых напряжений. Этот коэффициент учитывает потерю устойчивости и уменьшает предельную нагрузку поперечного сечения, как можно было видеть выше. Вот как его можно найти:
Коэффициент снижения допускаемых напряжений
  • Допустим, что греческие буквы в уравнениях могут вызывать беспокойство, но на самом деле это все довольно просто. Как ни странно, сказанное выше (и ниже) может быть получено непосредственно из дифференциальных уравнений деформации сжатой колонны. Эту формулу, при необходимости можно вывести самостоятельною В любом случае, для приведенного выше выражения требуются 2 вещи: гибкость, которую мы рассмотрим далее, и “загадочный” параметр, который просто сокращает представленное выше уравнение:
Коэффициент неоднородности
  • Обычно в МКЭ это подразумевает большой объем работы, но в приведенном выше описании все предельно просто. Сечения разделены на группы, и вы всего лишь берете значение этого параметра из таблицы (в зависимости от группы, которой принадлежит сечение). Конечно, в нашем случае это не так просто, поскольку мы имеем колонну с двумя поперечными сечениями - так что вы можете взять “худшее” из двух значений коэффициента, или просто “худшее” из имеющихся. Теперь можем вычислить гибкость.
Предельная пластическая нагрузка
  • Как видно, гибкость зависит только от предельной нагрузки поперечного сечения и критической силы. Итак, по сути мы знаем “все”, что нам нужно для расчета колонны… кроме, конечно, критической силы.

Критическая сила по Эйлеру

Если мы опустим детали того, как все это работает, то критическая сила Эйлера может быть сведена к такому предложению:

Критическая нагрузка на часть конструкции – это такая нагрузка, при действии которой эта часть выйдет из строя из-за упругой потери устойчивости.

Точка бифуркация

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

Femap. Модуль Юнга

Теперь появилось понимание, зачем нужна расчетная длина. Которая нужна, чтобы вычислить критическую нагрузку (которая требуется для расчета гибкости). Но из того же уравнения (и с учетом иллюстраций выше) получается, что нельзя получить расчетную длину для нашей колонны. И это потому, что:

Допущения для критической силы Эйлера:

  • Стержень идеально прямой
  • Поперечное сечение постоянно (то есть постоянный момент инерции по длине)
  • Стержень сжимается “только” вдоль оси
  • Материал колонны изотропный и однородный.

В нашем случае колонна, как единое целое не имеет критической силы Эйлера, поскольку она не “вписывается” в допущения.

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

Если не по Эйлеру, тогда что?

Тот факт, что для нашей колонны нет критической силы Эйлера, не означает, что она не потеряет устойчивость, это просто означает, что она не “вписывается” в предложенное Эйлером решение.

По-прежнему существует сила, которая при ее приложении вызовет упругую потерю устойчивости колонны, но в нашем случае мы не можем использовать уравнение Эйлера для вычисления этой силы. Такая сила будет являться критической, так как она вызовет упругую потерю устойчивости, но не будет критической силой по Эйлеру воспринимаемой как сила, вызывающая упругую потерю устойчивости, которую можно рассчитать по формуле Эйлера.

Примечательно, что в стандарте EN 1993-1-1 не говорится, что вы должны использовать “Критическую силу Эйлера” для расчета гибкости, в стандарте написано, что нужно использовать “Критическую силу”. Однако приведены допущения, указывающие, что вся методика применяется для прямых стержней с постоянным поперечным сечением и т.д. Но в действительности вам нужна критическая сила для расчета гибкости.

Если мы не можем использовать уравнение Эйлера, то нужно посмотреть что можно сделать? Конечно, мы можем применить МКЭ, а точнее LBA (Linear Bifurcation Analysis - Линейный анализ потери устойчивости). Благодаря этому инструменту мы сможем получить Критическую силу для любого случая.

Посмотрим принцип работы

В нашем случае - это колонна из профиля HEB 360 длиной 8 м. Почему так? Конечно потому, что такая задача идеально подходит для уравнения Эйлера. Таким образом, мы можем использовать как уравнение Эйлера, так и LBA для вычисления Критической силы. Конечно, предполагается, что мы получим одинаковое решение в обоих случаях!

Задача довольно проста, если вы хотите произвести вычисления вручную. Хотя колонне проще потерять устойчивость относительно ее слабой оси, я просто приму, что она закреплена таким образом, чтобы не деформироваться в этом направлении. Это ничего не меняет в самой методике вычислений – просто мне кажется, что потеря устойчивости относительно “сильной оси” выглядит лучше.

HEB 360

То же самое можно сделать с помощью LBA.

Поскольку это “балочная конструкция”, я буду использовать RFEM, программу, которую я по умолчанию использую для такого рода задач.

Я просто смоделировал 8-метровую колонку HEB 360 и сверху нагрузил ее свободный край сжимающей силой 1 кН. Так как я использовал коэффициент расчетной длины потери устойчивости 1.0, это отчасти приводит к “шарнирному характеру опирания”. Внизу имеется шарнирная опора (которая также запрещает вращение колонны вокруг оси, чтобы сделать ее устойчивой). Сверху запрещены только горизонтальные перемещения.

Результат анализа выглядит следующим образом (линейный статический расчет слева и LBA справа):

линейный статический расчет слева и LBA справа

Но более важно (в случае LBA), что коэффициент критической нагрузки, который получился, составил 13986. Это значение соответствует второй форме потери устойчивости (поскольку первая форма была для “слабой оси”, которую мы решили не учитывать). Поскольку приложенная нагрузка 1 кН, то это означает, что критическая сила составляет 1 кН x 13986 = 13986 кН.

Как видно, есть возможность рассчитать критическую силу с помощью LBA. Как это было наглядно показано, данный метод “согласуется” с уравнением Эйлера, поэтому методика “работает”. Но есть еще одна важная деталь. LBA работает для любой рассматриваемой модели, даже если она не соответствует допущениям по Эйлеру.

Вернемся к нашей задаче

Будем подробно моделировать соединение без учета эксцентриситет. Это не идеальное решение, поэтому, если хотите, можете поэкспериментировать и сделать более точную модель. Этого будет “достаточно” для большинства инженеров, которые решают подобные задачи.

Famap lpha_sr

Модель можно совершенствовать. Слева – случай “простой” 8-метровой балки из профиля HEB 360, затем добавили по 1 м профиля HEA 160 за один прием, получив в конце “простую” 8-метровую балку из профиля HEA 160, представленному справа. Над рисунками вы можете видеть коэффициент критической нагрузки для каждого случая (это, фактически, есть критическая сила в кН, поскольку нагрузка всегда составляла 1 кН).

Как видите, получение критической силы для нашего случая не является чем-то “невозможным”, даже несмотря на то, что мы не могли использовать формулу Эйлера.

Вот, что требовалось сделать, но мы еще не закончили.

График стоит 1000 слов

Поскольку у нас была простая модель, то показалось, что можно провести небольшое исследование. Как видно выше, вместо того, чтобы просто показать, как рассчитать критическую силу, мы фактически проанализировал различные колонны. Давайте посмотрим, как изменяется критическая сила, когда мы еще и еще будем сверху добавлять профиль HEA 160 к нашей 8-метровой колонне из профиля HEA 360:

Femap HEA 160 HEA 360

Мы добавили несколько расчетных точек (которые не показаны на изображениях выше), чтобы кривая выглядела более гладкой. Похоже, что очень короткий участок профиля HEA 160 (примерно до 500 мм) не вносит серьезных изменений. Проще говоря, он слишком короткий, чтобы оказывать влияние, а поскольку опора “шарнирно закреплена” наверху, в любом случае ничего особенного “там не происходит”.

Но все же, после достижения профилем HEA 160 высоты примерно 750 мм, “более слабый участок” начинает устраивать хаос! На приведенных выше изображениях даже можно видеть, что максимальные деформации возникают там, где имеется “более слабое” поперечное сечение, даже если оно составляет всего 25% от общей длины. Когда мы достигаем высоты колонны в 4 м (половины) с учетом замены на профиль HEA 160, это практически соответствует состоянию, при котором вся колонна состояла бы из более слабого поперечного сечения. Конечно, имеет место негативное воздействие (поскольку “половинчатая” колонна HEB 160/HEB 360 вдвое прочнее, чем “цельная” колонна HEB 160), но на самом деле мы в 10 раз уменьшили несущую способность, за счет замены поперечного сечения колонны на более слабое вплоть до половины длины.

Мы сделали представленный выше график, поскольку следует проверить, как все “меняется”, если выполнять такие преобразования. Я думаю, что периодическое проведение таких исследований может чему-то меня научить (я надеюсь, что и вас тоже!). Мне кажется, что если “верхняя”, более слабая, часть может быть “короткой”…, то хорошо бы делать ее как можно короче. Рентабельность усилий в таком случае наибольшая!

Преобразование критической силы в расчетную длину

Нам нужна расчетная длина для вычисления критической силы. Конечно, если бы мы делали расчет вручную, то могли бы просто получить критическую силу. Можно использовать RFEM вместо Femap для балочных моделей, потому что в RFEM уже реализована вся методика проектирования по Еврокоду. Поэтому можно отобразить процентное соотношение несущей способности в соответствии со стандартом EN 1993-1-1 в виде результирующего графика на балках/колоннах.

Нам нужно получить свой “коэффициент запаса”, то необходимо выполнить расчет по нормативной методике. И мы должен “сказать” программе, какова критическая сила для каждой части конструкции. В действительности дело обстоит иначе. Программное обеспечение, которое мы будем использовать, исходит из того, что мы не заинтересованы в самостоятельном вычислении критической силы. Оно просто запросит величину расчетной длины, а затем вычислит “остальное”. Но поскольку у нас есть критическая сила, и мы хотим ввести данные в программу, то нужно сначала “преобразовать их” в расчетную длину.

Но у нашей колонны нет расчетной длины

Верно, что наша колонна имеет критическую силу (но не критическую силу по Эйлеру, поскольку она не укладывается в рассмотренные допущения). Но программу не “волнуют” такие нюансы. Лучше всего выяснить это в руководстве к программе, но наше программное обеспечение хочет вычислить расчетную длину каждой балки, используя уравнение Эйлера. В конце концов, именно поэтому оно в первым делом спрашивает нас о расчетной длине. То что нам нужно сделать, это “видоизменить” уравнение Эйлера, чтобы получить:

Femap. Модуль Юнга, расчетная длина

Прекрасно, мы знаем критическую силу благодаря LBA. Теперь мы можем ввести все числа и вычислить то, что нам нужно. Однако, расчетная длина не будет иметь никакого “физического” смысла, поскольку наша колонна не соответствует допущениям о колонне Эйлера. Но это небольшая проблема! Программа возьмет нашу “бестолковую” расчетную длину и введем ее в уравнение Эйлера, которое в любом случае не подходит для нашего случая. Но результат, который получится по этому уравнению (критическая сила), в нашем случае будет правильным.

Конечно, прежде чем мы начнем использовать данный способ, нужно сделать один учебный пример и проверить, работает ли наша программа подобным образом.

В уравнении для расчетной длины необходимо использовать момент инерции. Он зависит от поперечного сечения, так что следует ли вам использовать маленькое сечение из верхнего профиля HEB 160 или большое сечение из профиля HEB 360. Конечно, нам нужно использовать оба сечения.

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

Каждая часть колонны имеет одну и ту же критическую силу, но проверяется независимо друг от друга

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

Как с этим справиться? На самом деле все очень просто

Давайте на секунду поразмыслим о верхней части.

Femap

Если бы мы сделали 8-метровую колонну только из профиля HEB 160, то она бы легче теряла устойчивость, чем рассмотренная выше. Иначе говоря, профиль HEB 360 (нижняя часть) более жесткий, и он “помогает” верхней части оставаться устойчивой. Это означает, что если мы сделаем колонну только из профиля HEB 160, имеющую такую же несущую способность по устойчивости (критическую силу), то она будет короче, чем указанная рассмотренная выше колонна. Она была бы высотой 5.15 м, вот вычисление:

Момент инерции профиля HEB 160

Конечно, мы можем сделать, то же самое и для нижней части. Но на этот раз это будет “хуже” 8 м. Иначе говоря, если бы вся колонна была из профиля HEB 360, то расчетная длина составила бы 8 м. Но поскольку верхняя часть “слабее”, она фактически делает всю конструкцию менее устойчивой. Таким образом, колонна, изготовленная исключительно из профиля HEB 360, будет иметь расчетную длину больше, чем те 8 м, чтобы обеспечить идентичную несущую способность по устойчивости.

А если точнее, расчетная длина “фиктивной” колонны из профиля HEB 360 должна быть: 21.43 м. Расчеты просты:

Момент инерции профиля HEB 360

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

Подводя итог, вот что мы сделали в этом примере шаг за шагом:

  • Смоделировали колонну и выполнили LBA, таким образом, мы можем получить множитель критической нагрузки, который даст нам критическую силу для нашей колонны.
  • Общая критическая сила. Все части колонны потеряют устойчивость “одновременно”. Вне зависимости от поперечных сечений все компоненты колонны всегда будут иметь общую критическую силу. Это потому, что колонна разрушается как одно целое, поэтому решающую роль играет самое слабое звено. Проще говоря, мы не можем выгнуть “часть” колонны. Колонна всегда выходит из строя из-за потери устойчивости, как единое целое от опоры до опоры (или от узла крепления до узла крепления).
  • Вычислили расчетную длину для каждой части колонны. Когда мы понимаем, что все компоненты колонны выгибаются “одновременно”, все становится намного проще. Для обеих частей мы можем обратным пересчетом получить длину “гипотетической” колонны постоянного поперечного сечения, которая должна иметь такую же несущую способность по устойчивости, как и наша исследуемая колонна. Благодаря этому можно спроектировать нашу колонну как “совокупность” “гипотетических” колонн, используя программный алгоритм, реализованный в КЭ программном обеспечении.