Overall stability of compressed compound bars in variable cross section

Full Text

Введение. Расчет на общую устойчивость стержней сплошного или составного переменного сечения на основе технической или уточненной (балка Тимошенко) теорий представляет интерес и в настоящее время, если иметь в виду построение замкнутых (или математическое обоснование) приближенных методов ее решения [1-9]. Следует отметить практическую значимость континуальной расчетной модели центрально-сжатого составного стержня переменного сечения, применяемой при расчете на общую устойчивость гибких дискретных стержневых систем башенного типа, в частности радио и телевизионных мачт, стрел кранов, буровых вышек, опор ЛЭП и эстакад [5, 10]. Вопросы обеспечения общей устойчивости таких систем при их проектировании играют важную роль в связи с применением в качестве материала для них высокопрочных сталей, а также легких сплавов, обладающих высокими деформативными свойствами. Условно сдвигово-изгибная расчетная модель стержня может успешно использоваться и при исследовании устойчивости плоских рам [11]. Наиболее полно устойчивость сплошных стержней переменного сечения со степенным законом изменения жесткости исследовалась А. Н. Динником [1], а в статье [4] приведены численные результаты расчетов для таких стержней. Подобная задача рассматривалась также в [2, 3] и некоторых других. Вместе с тем следует отметить, что результаты этих исследований проведены без учета деформаций поперечного сдвига, играющих заметную роль именно в стержнях составного сечения. Исключением являются статьи одного из соавторов [5, 10, 12], в которых получены замкнутые решения задач устойчивости для составных стержней с квадратичным законом изменения моментов инерции при их загружении полярными силами, а также приближенные решения для балки Тимошенко сплошного сечения [6-7]. В отличие от приведенных выше исследований, настоящая работа посвящена построению замкнутого (в пределах принятых допущений) решения задачи общей устойчивости для составного стержня переменного сечения. Для этой цели предлагается гипотетическая модель стержня сплошного сечения, эквивалентная по деформациям изгиба и единичным углам сдвига исследуемой сквозной стержневой конструкции. Существенным представляется и то, что в дальнейшем используется уточненная авторами модель С. П. Тимошенко. Такой подход позволил сформулировать критерий, при котором оказывается справедливой и традиционная схема учета деформаций сдвига. Следует также подчеркнуть предлагаемую авторами настоящей работы корректную постановку рассматриваемой краевой задачи, основанную на дифференциальных уравнениях как 2-го, так и 4-го порядков, что при соответствующих граничных условиях позволило исследовать различные способы закрепления составного стержня, и не только симметричные, но и несимметричные формы потери устойчивости. В работе для этой цели используется прием повышения порядка разрешающего дифференциального уравнения. Вместе с тем, как отмечалось выше, авторы обычно пользовались лишь разрешающим уравнением 2-го порядка, что не позволило им исследовать все возможные формы потери устойчивости равновесия. 1. Постановка задачи. Рассматривается устойчивость центрально сжатого стержня переменного сечения (рис. 1). Сжимающая сила P считается приложенной строго вдоль центральной оси стержня, которая остается прямолинейной вплоть до момента потери устойчивости (критического состояния Pcr ). Таким образом, потеря устойчивости сопровождается разветвлением форм равновесия (бифуркацией). Поскольку рассматриваемая ниже механическая система является консервативной, используется статический критерий потери устойчивости, основанный на уравнениях равновесия, составленных в деформированном состоянии стержня. При этом многопанельный составной стержень заменяется эквивалентным по единичным углам сдвига и изгибной жесткости стержнем сплошного сечения. Деформации поперечного сдвига в соответствии с гипотезой Энгессера-Тимошенко определяются законом из Рис. 1. Общая схема к задаче исследования устойчивости стержня переменного сечения [Figure 1. The general scheme to the stability task of solving a rod of variable section] менения перерезывающих сил, а деформации сжатия, ввиду их малости, не учитываются. Как правило, ввиду постоянства поперечного сечения ветвей составных стержней с непрерывно изменяющейся по степенному закону изгибной жесткостью наибольший практический интерес представляет случай, когда Ax = A = const, Jx (x) = Jmax (x/b)α , (1) где Ax , Jx (x) - соответственно площадь и момент инерции поперечного сечения составного стержня. Остальные обозначения приведены на рис. 1. Принимая во внимание сформулированные допущения, рассматриваемая задача представляет задачу о бифуркации изгибных форм равновесия. Исходными соотношениями при этом являются уравнение совместности деформаций: dy = ψ + γ, (2) dx физические уравнения, связывающие усилия и деформации: dψ = -Mx , dx µ µ γ= Qx = γ Qx , γ = ¯ ¯ , GA GA EJx (3) (4) а также уравнения равновесия в деформированном состоянии стержня (рис. 1): Mx = P y + H(x - a) + M0 , Qx = P ψ + H. (5) (6) Здесь y, ψ, γ - соответственно прогибы, углы поворота и углы сдвига в том же сечении стержня; Mx , Qx - изгибающий момент и перерезывающая сила; H, M0 - реактивные усилия в сечении x = a; E, G - модули упругости и сдвига; µ, γ - соответственно коэффициент неравномерности распределения ¯ касательных напряжений по поперечному сечению и единичный угол сдвига от Q = 1. В силу малости деформаций при составлении уравнений равновесия было принято (рис. 1, b): sin ψ ∼ ψ, cos ψ ∼ 1. Дифференцируя равенство (2), = = находим dψ dγ d2 y + . (7) EJx 2 = EJx dx dx dx В результате подстановки в (7) продифференцированного соотношения (4), имея при этом в виду (6), получаем EJx d2 y dψ = EJx (1 + γ P ) ¯ . 2 dx dx (8) Выражение (3) с учетом (5) принимает следующий вид: EJx dψ = - P y + H(x - a) + M0 . dx (9) Последнее равенство с учетом зависимости (1) представляет собой дифференциальное уравнение устойчивости составного стержня переменного сечения со степенным законом изменения жесткости при наличии деформаций сдвига (случай I): x b EJmax (1 + γ P )-1 ¯ α d2 y(x) dx2 + P y(x) = -H(x - a) - M0 . (10) Следует отметить, что приведенный выше вывод дифференциального уравнения устойчивости (10) основан на независимых (5) и (6) уравнениях равновесия. Вместе с тем в обычной трактовке для этой цели вместо (6) используется известная зависимость Д. И. Журавского. С учетом (5) имеем Qx = dMx dy =P + H. dx dx (11) Повторяя рассуждения (7)-(9), при наличии (11) получаем в этом случае такое дифференциальное уравнение устойчивости (случай II): EJmax (1 - γ P ) ¯ 344 x b α d2 y(x) dx2 + P y = -H(x - a) - M0 . (12) Добавляя к (10) или (12) граничные условия, определяемые соответствующими условиями закрепления концов стержня, получаем математическую формулировку рассматриваемой краевой задачи. Легко установить связь между уравнениями (10) и (12). Действительно, из степенного разложения (1 + γ P )-1 = 1 - γ P + (¯ P )2 - (¯ P )3 + . . . , ¯ ¯ γ γ замечаем, что (1+ γ P )-1 ∼ 1- γ P при γ P ¯ ¯ ¯ 1. Таким образом, когда γ P ¯ 1, = уравнение (10) тождественно совпадает с (12). Следует обратить внимание, что уравнениям (10) и (12) соответствует одно однородное дифференциальное уравнение b2-α xα d2 y + ωy = 0, dx2 в котором используются различные критические параметры ω: ω = P b2 (1 + γ P ) (EJmax )-1 ¯ ¯ ω = P b2 (1 - γ P )-1 (EJmax ) -1 (в случае I), (в случае II). (13) (14) 2. Метод решения. Решение задачи осуществляется непосредственным интегрированием уравнения (10) или (12). Известно, что общее решение линейного неоднородного дифференциального уравнения с переменными коэффициентами может быть представлено в виде суммы: y(x) = y1 (x) + y2 (x), (15) где y1 (x) - общее решение соответствующего (10) или (12) однородного уравнения, а y2 (x) - какое-либо частное решение неоднородного уравнения (10). Общий интеграл однородного уравнения d2 y1 (x) + bα-2 ωx-α y1 (x) = 0 dx2 (16) √ √ y1 (x) = C1 xJν (¯) + C2 xYν (¯), x x (17) можно получить в виде где x= ¯ 1 1/2 x ω k b k , k= 2-α , 2 ν = |2 - α|-1 . (18) Jν (¯), Yν (¯) - цилиндрические функции (функции Бесселя) 1-го и 2-го рода x x порядка ν, а C1 , C2 - постоянные интегрирования. Частное решение y2 (x) определяется по виду правой части (10): y2 (x) = - 1 H(x - a) + M0 . P (19) Складывая (17) и (19) в соответствии с (15), окончательно имеем √ √ 1 H(x - a) + M0 , y(x) = C1 xJν (¯) + C2 xYν (¯) - x x P (20) при этом H и M0 определяются условиями закрепления верхнего конца стержня. Следует отметить особый случай α = 2, соответствующий квадратичному закону изменения жесткости (1), при котором согласно (18) ν → ∞ и, естественно, полученное выше общее решение (17) и (20) не проходит. Рассмотрим этот случай отдельно. При α = 2 уравнение (16) представляет собой уравнение Эйлера d2 y1 (x) + ωx-2 y1 (x) = 0, dx2 решение которого записывается в виде √ ¯ ¯ y1 (x) = x C1 cos x1 + C2 sin x1 , ¯ где (21) √ 4ω - 1 ln x. (22) 2 При этом критический параметр ω > 0.25, поскольку случай ω < 0.25 соответствуют нереально большим гибкостям стержня. Частное решение неоднородного уравнения (10) в случае α = 2 определяется по формуле (19), а его общее решение в соответствии с (15), (19), (21) - равенством x1 = ¯ y(¯) = x √ x C1 cos x1 + C2 sin x1 - ¯ ¯ 1 H(x - a) + M0 . P (23) В дальнейшем искомым является критический параметр ωcr , зная который, сжимающая критическая сила может быть определена из соотношений (13) или (14) (случаи I и II). Случай I. Из равенства (13) следует квадратное уравнение [13], корни которого 1 1 EJmax Pcr = - ± + ωcr , 2¯ γ 4¯ 2 γ γ b2 ¯ или в более удобной форме: Pcr = π 2 EJmax , (β1 l)2 (24) где l - длина стержня (см. рис. 1); β1 - коэффициент приведенной длины стержня: β1 = β 346 ωcr EJmax b-2 1 ωcr EJmax + 2 4¯ γ γ b2 ¯ 1/2 - 1 2¯ γ -1 . (25) Общая устойчивость сжатых составных стержней. . . Здесь β= 1 π √ , 1 - k ωcr k= a , b (26) причем β - безразмерный коэффициент приведенной длины, определяемый без учета деформаций сдвига решетки [4]. Вид ωcr = ωcr (l) будет определен ниже. Случай II. Аналогично из (14) имеем Pcr = ωcr EJmax 2 + γ EJ b ¯ max ωcr где β1 = или Pcr = π 2 b + ωcr γ EJmax ¯ l π 2 EJmax , (β1 l)2 -1/2 . (27) Здесь также ωcr = ωcr (l). 3. Решение в случае симметричной формы потери устойчивости. Поскольку кривая, определяющая форму потери устойчивости в зависимости от условий закрепления стержня, должна удовлетворять четырем граничным условиям (по два на каждом конце стержня), полученное выше решение (20) или (23) в общем случае при наличии двух произвольных постоянных C1 и C2 не может это обеспечить, поэтому оно должно быть расширено. Следует отметить, что равенства (20) или (23) можно использовать при исследовании симметричных форм потери устойчивости стержней переменного сечения. Естественно, сам стержень при этом должен быть конструктивно симметричным (рис. 2). В дальнейшем понадобится выражение для углов наклона касательной к деформированной оси стержня, определяющей форму потери устойчивости. Дифференцируя выражения (17) и (21) в случае α = 2, получим 1 n d¯ x dy1 (x) = C1 x-1/2 Jn (¯) + x1/2 Jn (¯) - Jn+1 (¯) x x x + dx 2 x ¯ dx d¯ x 1 n x x x , (28) + C2 x-1/2 Yn (¯) + x1/2 Yn (¯) - Yn+1 (¯) 2 x ¯ dx где d¯ x x (1-2n)/2n 1 = ω 1/2 . dx a a А в случае α = 2 получим dy1 (x) 1 d¯1 x = C1 x-1/2 cos x1 - x1/2 sin x1 ¯ ¯ + dx 2 dx + C2 x-1/2 sin x1 + x1/2 cos x1 ¯ ¯ d¯1 x , (29) dx причем d¯1 x 1 1 = (4ω - 1)1/2 . dx 2 x (30) Рис. 2. Расчетная схема стержня при симметричных формах потери устойчивости [Figure 2. The design scheme for a rod with symmetric cases buckling] Рассмотрим ниже некоторые примеры расчета. 1. Симметричный составной стержень при произвольном α = 2 шарнирно закреплен по концам. Соответствующими граничными условиями при этом являются условия симметрии формы потери устойчивости (рис. 2): dy1 = 0. (31) dx x=b Подстановка равенств (17), (28) в соотношения (31) формирует следующую однородную систему алгебраических уравнений для C1 и C2 : y1 (a) = 0; C1 Jn (¯) + C2 Yn (¯) = 0, a a 1 -1/2 n C1 b + b1/2 ¯ Jn (¯ - b1/2 Jn+1 (¯ + b) b) 2 b 1 -1/2 n + C2 b + b1/2 ¯ Yn (¯ - b1/2 Yn+1 (¯ = 0. b) b) 2 b 348 (32) Общая устойчивость сжатых составных стержней. . . Здесь a= ¯ 2 ω 1/2 , 2-α ¯= b 2 ω 1/2 (b/a)(2-α)/2 , 2-α n= 1 . |2 - α| (33) Разыскивая нетривиальные решения для C1 и C2 , приравниваем главный детерминант системы (32) к нулю. В результате получаем Jn (¯) a 1 -1/2 n b + b1/2 ¯ Yn (¯ - b1/2 Yn+1 (¯ = b) b) 2 b n 1 -1/2 b + b1/2 ¯ Jn (¯ - b1/2 Jn+1 (¯ . (34) = Yn (¯) a b) b) 2 b Трансцендентное уравнение устойчивости (34) позволяет определить итерационным путем наименьший корень ωmin - критический параметр, а затем по формулам (26), (24) - коэффициент приведенной длины β1 и критическую силу Pcr . 2. Аналогичное решение задачи в случае α = 2 при наличии равенств (22), (29) и (31) приводит к соотношениям C1 cos a1 + C2 sin a1 = 0 при x = a, ¯ ¯ 1 1 b b b C1 cos ¯1 - (4ω - 1)1/2 sin ¯1 + C2 sin ¯1 + (4ω - 1)1/2 cos ¯1 b 2 2 при x = b, где a1 , ¯1 определяются по формулам (33). ¯ b Трансцендентное уравнение для определения ωmin принимает вид 1 1 cos a1 sin ¯1 + (4ω - 1)1/2 cos ¯1 = sin a1 cos ¯1 - (4ω - 1)1/2 sin ¯1 ¯ b b ¯ b b 2 2 и уточняет соответствующее ошибочное уравнение, приведенное в монографии [3]. 3. Рассмотрим составной стержень, жестко защемленный по концам при α = 2. Из-за симметрии формы потери устойчивости граничные условия запишутся следующим образом: dy1 dx x=a = dy1 dx x=b = 0. (35) В результате подстановки выражения (28) в (35) получаем однородную систему уравнений: C1 C1 1 -1/2 n a + a1/2 Jn (¯) - a1/2 Jn+1 (¯) + a a 2 a ¯ 1 -1/2 n +C2 a + a1/2 Yn (¯) - a1/2 Yn+1 (¯) = 0, a a 2 a ¯ 1 -1/2 n b + b1/2 Jn (¯ - b1/2 Jn+1 (¯ + b) b) 2 a ¯ 1 -1/2 n +C2 b + b1/2 ¯ Yn (¯ - b1/2 Yn+1 (¯ = 0, b) b) 2 b 349 С е н и ц к и й Ю. Э., И ш у т и н А. С. откуда затем следует трансцендентное уравнение для определения ωmin : = 1 -1/2 1/2 n 1 -1/2 1/2 n a +a Jn (¯)-a1/2 Jn+1 (¯) a a b +b ¯ Yn (¯ b)-b1/2 Yn+1 (¯ = b) 2 a ¯ 2 b 1 -1/2 1/2 n 1 -1/2 1/2 n b +b ¯ Jn (¯ a +a Yn (¯)-a1/2 Yn+1 (¯) . a a b)-b1/2 Jn+1 (¯ b) 2 2 a ¯ b Здесь a, ¯ и n определяются по формулам (33). ¯ b 4. В случае α = 2 по той же схеме c использованием равенств (30), (35) (соответственно при x = a и x = b) получаем 1 1 C1 cos a1 - (4ω - 1)1/2 sin a1 + C2 sin a1 + (4ω - 1)1/2 cos a1 = 0, ¯ ¯ ¯ ¯ 2 2 1 1 C1 cos ¯1 - (4ω - 1)1/2 sin ¯1 + C2 sin ¯1 + (4ω - 1)1/2 cos ¯1 = 0. b b b b 2 2 Приравнивая главный детерминант этой системы к нулю, получаем уравнение для определения ωmin : 1 1 cos a1 - (4ω - 1)1/2 sin a1 sin ¯1 + (4ω - 1)1/2 cos ¯1 = ¯ ¯ b b 2 2 1 1 b ¯ ¯ = cos ¯1 - (4ω - 1)1/2 sin ¯1 sin a1 + (4ω - 1)1/2 cos a1 . (36) b 2 2 Равенство (36) уточняет соответствующее ошибочное уравнение, приведенное в монографии [3]. 4. Обобщение решения рассматриваемой краевой задачи. На основе решений (17), (21) разрешающего дифференциального уравнения втрого порядка (10), (12), были сформулированы (п. 3) трансцендентные уравнения устойчивости для критического параметра ωcr = ωmin центрально сжатого составного стержня со степенным законом изменения жёсткости. При этом получены результаты для заранее (априори) известных симметричных форм потери устойчивости. Однако в случае произвольных несимметричных форм исследовать устойчивость составного стержня, пользуясь соотношениями (17), (21), не представляется возможным из-за переопределённости системы граничных условий (по два условия на каждом конце для всех схем закрепления стержня). Для этой цели предложена методика, связанная с повышением порядка разрешающего дифференциального уравнения и дальнейшей корректной постановкой рассматриваемой краевой задачи. Дифференцируя дважды уравнения (10) и (12), получаем y (4) (x) + 2α α(α - 1) ω y (x) + + α y (x) = 0. 2 x x x (37) Здесь в соответствии с (13), (14) ω представляет собой критический параметр преобразованных уравнений (10) и (12): ω= 350 P bα EJmax (1 + γ P ) ¯ (1 - γ P )-1 ¯ В дальнейшем будем рассматривать квадратичный закон изменения жёсткости стержня (α = 2). Построим для него новое замкнутое решение задачи. При α = 2 дифференциальное уравнение (37) принимает вид y (4) (x) + 4 2+ω y (x) = 0. y (x) + x x2 (38) Введём новую функцию y (x) = z(x) (39) и запишем уравнение (38) так: z (x) + 4 2+ω z (x) + z(x) = 0. x x2 (40) Общее решение уравнения (40) имеет вид z(¯) = exp(-n¯)[C1 cos m¯ + C2 sin m¯], x x x x (41) где x = ln x, ¯ 3 n= , 2 1 m = (4ω - 1)1/2 , 2 (42) ω > 1/4. Случай ω < 1/4 соответствует нереально большим гибкостям стержня и в дальнейшем не рассматривается. Из дифференциального уравнения (39) находим y (x) = exp -(1 + n)¯ -C1 (n cos m¯ + m sin m¯)+ x x x + C2 (cos m¯ - n sin m¯) ; x x y (x) = 1 + n)2 -1 (1 - m2 x m-1 C1 sin m¯ + 1-n cos m¯ + x m 1-n sin m¯ - cos m¯ + C3 , x x m -1 -1 (1 - n)2 (2 - n)2 y(x) = 1 + 1- m-1 × m2 m2 2-n (sin m¯ + cos m¯) + x x × C1 (sin m¯ - cos m¯) + x x m 2-n +C2 (sin m¯ - cos m¯) + x x (cos m¯ + sin m¯) + C3 x + C4 . x x m + C2 (43) Соотношения (41) и (43) представляют собой решение рассматриваемой задачи, которое справедливо для произвольных условий закрепления составного стержня с квадратичным законом изменения жёсткости. Добавляя в каждом конкретном случае однородные граничные условия, получаем соответствующую краевую задачу на собственные значения. Ниже сформулируем уравнения устойчивости для различных условий закрепления составного стержня при α = 2. 1. Шарнирное закрепление обоих концов. В этом случае на обоих концах перемещения и моменты равны нулю. Принимая во внимание соотношения (3) и (8), имеем граничные условия y(a) = y(b) = 0, y x=a =y x=b = 0. (44) Подставляя выражения (41), (43) в равенства (44), получаем однородную систему алгебраических уравнений относительно постоянных C1 , C2 , C3 , C4 . Разыскивая затем её нетривиальные решения, приравниваем нулю главный детерминант. В результате получаем [B(ω, a) + D(ω, a)] [B(ω, a) + D(ω, a)] E(ω)a E(ω) ¯ ¯ ¯ ¯ cos m¯ a sin m¯ a 0 0 = 0. B(ω, ¯ + D(ω, ¯ b) b) B(ω, ¯ + D(ω, ¯ b) b) E(ω)b E(ω) cos m¯ b sin m¯ b 0 0 (45) Здесь 2-n (sin m¯ + cos m¯), x x m (1 - n)2 (2 - n)2 E(ω) = 1 + 1+ m, m2 m2 где n, m вычисляются по формулам (42), а a и ¯ находятся из равенств ¯ b B(ω, x) = sin m¯ - cos m¯, ¯ x x D(ω, x) = ¯ ¯ = ln b. b a = ln a, ¯ (46) Трансцендентное уравнение устойчивости (45) решается итерационным путём. Разыскивается минимальный корень ωmin = ωcr , а затем из равенств (25), (26) и (24) - коэффициенты приведенной длины и критическая нагрузка. 2. Жёсткое защемление обоих концов. В этом случае граничные условия формулируются следующем образом: y(a) = y(b) = 0, y (x) x=a = y (x) x=b = 0. (47) Подставляя выражения (43) в равенства (47), получаем однородную систему, приравниваем нулю главный детерминант и получаем трансцендентное уравнение устойчивости для ωmin = ωcr : [B(ω, a) + D(ω, a)] [B(ω, a) + D(ω, a)] E(ω)a E(ω) ¯ ¯ ¯ ¯ F (ω, a) ¯ G(ω, a) ¯ H(ω) 0 = 0, ¯ + D(ω, ¯ ¯ + D(ω, ¯ B(ω, b) b) B(ω, b) b) E(ω)b E(ω) F (ω, ¯ b) G(ω, ¯ b) H(ω) 0 где F (ω, x) = sin m¯ + ¯ x 1-n cos m¯, x m H(ω) = 1 + 352 G(ω, x) = ¯ 1-n sin m¯ - cos m¯, x x m (1 - n)2 m. m2 (48) Общая устойчивость сжатых составных стержней. . . 3. Консольный стержень. Граничные условия задачи для этого случая имеют вид y (a) = y (a) = 0, y(b) = y (b) = 0. (49) Используя выражения (42), (43) и (49), в соответствии с описанной выше процедурой получаем следующее трансцендентное уравнение устойчивости для критического параметра ωmin = ωcr : cos m¯ a K(ω, a) ¯ B(ω, ¯ + D(ω, ¯ b) b) F (ω, ¯ b) sin m¯ a L(ω, a) ¯ B(ω, ¯ + D(ω, ¯ b) b) G(ω, ¯ b) 0 0 0 0 = 0. E(ω)b E(ω) H(ω) 0 (50) Здесь K(ω, x) = n cos m¯ + m sin m¯, L(ω, x) = n sin m¯ - m cos m¯. ¯ x x ¯ x x 4. Случай, когда один конец защемлён, а другой шарнирно закреплён. Граничные условия формулируются следующим образом: y(a) = y (a) = 0, y(b) = y (b) = 0. Применяя описанный выше алгоритм, получаем следующее уравнение устойчивости для ωmin = ωcr : [B(ω, a) + D(ω, a)] [B(ω, a) + D(ω, a)] E(ω)a E(ω) ¯ ¯ ¯ ¯ cos m¯ a sin m¯ a 0 0 = 0. F (ω, ¯ + D(ω, ¯ b) b) G(ω, ¯ + D(ω, ¯ b) b) E(ω)b E(ω) F (ω, ¯ b) G(ω, ¯ b) H(ω) 0 (51) Таким образом, в зависимости от условий закрепления стержня критический параметр ωcr является наименьшим корнем одного из полученных трансцендентных уравнений (45), (48), (50), (51). 5. Численный анализ результатов. Исходными данными для проведения расчётов устойчивости центрально-сжатого составного стержня переменного сечения с квадратичным законом изменения жёсткости α = 2 являются модуль упругости материала E, площадь поперечного сечения ветви Ab , а также соотношения площадей поперечных сечений горизонтальных поясов и раскосов An /Ab и Ap /Ab решётчатого стержня. Известными считаются Jmax , l (рис. 1) размеры dH , db и число панелей i многопанельного стержня. С использованием формулы (1) при x = a и соотношения (46) вычисляются параметры a, b, a, ¯ ¯ b: a=l db , dH - db b = l + a, a = ln a, ¯ ¯ = ln b. b Единичные углы сдвига определяются для каждой панели стержня (рис. 3) из соотношения δ11 γi ≈ tg γi = ¯ ¯ , (l/i) ¯ где δ11 - единичное перемещение по направлению Q. За расчётное γi прини¯ мается среднее арифметическое значение. Рис. 3. Схема к определению единичных углов сдвига для панели [Figure 3. The design scheme for the definition of angular shift for a panel] В зависимости от способов закрепления обоих концов стержня итерационным методом решается соответствующее трансцендентное уравнение (45), (48), (50) или (51). Для этой цели задаёмся ωk и, варьируя его, добиваемся удовлетворения указанных уравнений. Определив наименьший корень ωkp = ωmin по формулам (25) или (27), вычисляем коэффициент приведённой длины β или β1 , а затем из выражения (24) определяем критическую силу Pcr . В качестве модельного примера рассматривается десятипанельный стальной стержень решётчатого типа пролётом l = 2.0 м, d = 0.2 м и db = 0.045 м. Каждая из четырёх ветвей стержня трубчатого сечения диаметром 0.01 м с толщиной стенки ∆ = 0.001 м. Другие исходные данные: E = 2.1·1011 Н/м2 , Ab = 0.126 · 10-3 м2 , Jmax = 0.126 · 10-5 м4 , Ap /Ab = 0.1, An /Ab = 0.2 (Ab , Ap , An - соответственно площади поперечных сечений ветвей, раскосов и поясов). Для указанного типа решётки единичный угол сдвига составил γ = (4.51/E) · 103 Н-1 . ¯ На рис. 4 приведены графики коэффициента приведенной длины β для составного стержня (α = 2) при различных условиях закрепления его концевых сечений, полученные в результате варьирования деформаций сдвига решётки γ . Кривые 1-4 соответствуют консольному, шарнирному, жёстко за¯ щемлённому и смешанному способам закрепления соответственно. Рис. 4. Значения коэффициента приведенной длины β для составного стержня (α = 2) для различных условий закрепления: 1 - консольное; 2 - шарнирное; 3 - жёсткое; 4 - смешанное [Figure 4. The values of the coeffcient β for the composite rod (α = 2) under different conditions of fixing: 1 - single-sided support; 2 - hinge support; 3 - rigid fixing; 4 - combined fixing] Следует отметить, что значения β превышают аналогичные значения коэффициентов приведенной длины, подсчитанные при γ = 0. При этом вли¯ яние деформаций сдвига решётки на β возрастает с увеличением жёсткости закрепления концов стержня и становится максимальным для жёстко защемлённого стержня. В отличие от стержней сплошного сечения, при оценке устойчивости составных стержней учёт деформаций сдвига решётки является необходимым. В заключение можно сделать следующие выводы. 1. Сформулирована краевая задача устойчивости для центрально-сжатых составных стержней переменного сечения со степенным законом изменения жесткости, основанная на уточненных авторами соотношениях теории С. П. Тимошенко. 2. Получены впервые обобщенные трансцендентные уравнения устойчивости справедливые для исследования как симметричных, так и несимметричных форм равновесия составных стержней переменного сечения. 3. Показано, что для рассматриваемых стержней деформации поперечного сдвига решетки оказывают заметное влияние на окончательные результаты и должны учитываться при проведении расчётов в прикладных задачах.

About the authors

Yuri E Senitsky

Samara State University of Architecture and Civil Engineering

Email: senitskiy@mail.ru
Dr. Techn. Sci.; senitskiy@mail.ru; Corresponding Author), Head of Dept., Dept. of Strength of Materials and Structural Mechanics 194, Molodogvardeyskaya st., Samara, 443001, Russian Federation

Alexander S Ishutin

Samara State University of Architecture and Civil Engineering

Student, Dept. of Strength of Materials and Structural Mechanics 194, Molodogvardeyskaya st., Samara, 443001, Russian Federation

References

Supplementary files