SOLVING THE PROBLEM OF LONGITUDINAL BENDING OF A ROD BY METHODS OF DISCRETE CALCULUS OF VARIATIONS

Abstract


The question of the longitudinal bending of a rod consisting of rigid links connected by hinges is being studied. It is shown that, as in the classical version of a solid rod, the problem can be posed as a variational minimum energy problem, but given by a functional defined on a class of functions with a discrete domain of definition. The main provisions of the classical calculus of variations are transferred to functionals of this type: the formula of variation is found, a generalization of the main lemma of the calculus of variations is proved, an analogue of the Euler equation, which is a difference equation, is obtained. Applying the results obtained and the known properties of classical difference equations, we succeeded in solving an analogue of the Euler problem for two types of a hinge rod: for a rod consisting of links of the same length, and for an arbitrary choice of lengths of links. In both cases we find the critical Euler force, as well as the equation and the form of the deflection curve.

Full Text

Введение Вариационные принципы при решении задач механики используются давно: хорошо известен, например, принцип наименьшего действия, который утверждает, что любой оптимальный физический процесс (т.е. реализуемый природой) происходит с наименьшими затратами энергии. Соответствующая математическая постановка означает, что функционал, выражающий суммарную энергию, для оптимального процесса должен достигать минимума. Развитие вариационного исчисления как самостоятельной математической дисциплины позволило решить ряд задач, которые сейчас считаются классическими: в частности, задача Эйлера о продольном изгибе стержня может быть поставлена и решена как вариационная задача. Большинство вариационных задач ставится на множестве непрерывных функций и, как следствие, большинство исследуемых функционалов имеют вид интегралов. Заметим, однако, что такой выбор пространства функций и вида функционалов не является единственно возможным. Есть задачи, в которых исходный аргумент естественнее считать дискретным, что автоматически приводит к замене непрерывной функции последовательностью, а интеграла - суммой. В данной работе рассматривается задача, которую можно считать дискретным аналогом задачи Эйлера о продольном изгибе стержня. Для решения этой задачи оказалось удобным изначально считать аргумент дискретным, а вместо интеграла энергии записывать сумму. Исследование этого функционала на экстремум потребовало переноса на дискретный случай известных утверждений классического вариационного исчисления, что составило первую часть настоящей работы. Во второй части работы полученные результаты применяются к исследованию задачи о продольном изгибе шарнирного стержня. 1. Описание объекта и постановка задачи Напомним формулировку и решение классической задачи Эйлера о продольном изгибе стержня. Задача 1. Стержень длиною опирается своими концами и подвержен давлению P (рис. 1). При определённом значении P (критическая сила Эйлера) происходит продольный изгиб стержня. Требуется определить наименьшую величину силы P, дающую продольный изгиб. Рис. 1. Непрерывный стержень Как известно из [1], решение задачи 1 сводится к определению минимума интеграла энергии , где E - модуль упругости, I - наименьший момент инерции поперечных сечений стержня, - угол касательной кривой прогиба с горизонтальной осью. Функционал дополнен граничными условиями . Решая данную задачу как вариационную, находим наименьшее значение критической силы Эйлера . Уравнение кривой прогиба выглядит следующим образом: . Теперь рассмотрим задачу, которую можно считать дискретным аналогом задачи 1. Задача 2. Стержень длиною разделён шарнирами на единичные участки. Стержень опирается своими концами и подвержен давлению P (рис. 2). Шарниры, находящиеся на одинаковом расстоянии друг от друга, считаем идеальными, т.е. они дают стержню свободно изгибаться. При определённом значении P (критическая сила Эйлера) происходит продольный прогиб стержня. Требуется определить наименьшую величину силы P, дающую продольный изгиб. Рис. 2. Дискретный стержень Как будет показано ниже, задачу 2 тоже можно рассматривать как вариационную, но возникающий там функционал энергии не является интегральным, так как область определения искомых функций в такой постановке удобно считать целочисленной. Для решения этой (дискретной) задачи следует разработать соответствующий математический аппарат - дискретное вариационное исчисление. 2. Дискретное вариационное исчисление 2.1. Вариация функционала. Необходимое условие экстремума Напомним основные теоремы и определения из классического вариационного исчисления. Пусть - произвольное линейное нормированное пространство. Если каждому элементу ставится в соответствие по определенному закону единственное вещественное число , то говорят, что на определен функционал . Кратко этот факт отображают соотношением . Пространство называют областью определения функционала. Так как значения функционала - числа, то на функционалы легко переносятся понятия наибольшего и наименьшего значения, которые в случае их достижимости называют максимумом и минимумом. Для краткости речи эти понятия обозначают одним термином - экстремумы. Определение 1. Точка называется точкой локального максимума (соответственно - минимума) функционала , если для некоторого при всех , для которых , выполнено неравенство (соответственно - ). Определение 2 [2]. Вариацией функционала в точке называется число , определяемое по формуле (1) Здесь , принимает значения в некотором интервале, содержащем точку Теорема 1 [2]. Пусть в точке существует вариация . Тогда, если - точка экстремума функционала , то Теорема 1 является аналогом необходимого условия экстремума для функций. При переходе от функций к функционалам роль производной начинает играть вариация функционала. Эта теорема лежит в основе классического вариационного исчисления. Заметим, что в теореме 1 не накладывается никаких условий ни на пространство , ни на вид функционала , но в большинстве учебников и монографий по вариационному исчислению выбирается в качестве пространство дифференцируемых до некоторого заданного порядка функций, на которых рассматриваются функционалы, имеющие вид определенных интегралов (однократных, двойных, тройных и т. д.). Простейшей из вариационных задач считается так называемая задача с закрепленными границами: найти в указанном классе функций такую, на которой функционал принимал бы наибольшее или наименьшее из всех возможных значений, т.е. (2) Теорема 2 [2]. Если - точка экстремума вариационной задачи (2), то она является решением уравнения Эйлера: . 2.2. Дискретная вариационная задача и аналог уравнения Эйлера Очевидно, что выбирать в качестве пространство непрерывных на отрезке функций, а в качестве функционала определенный интеграл - это только одна из многих возможностей. Предположим, что множество, на котором определяются допустимые функции - не отрезок , а набор целочисленных точек , на которых определяются вещественные числа , т.е. задаются функции целочисленного аргумента. Для таких функций в задаче (2) интеграл естественно заменить суммой, и мы получаем дискретный аналог вариационной задачи (2): (3) где - конечная разность. Покажем, что на задачу (3) можно перенести ряд основных теорем из классического вариационного исчисления. Из теоремы 1 следует, что для определения точек экстремума функционала следует найти его вариацию по формуле (1) и приравнять ее к нулю. Обозначим , , и найдём вариацию функционала (3), используя свойства конечных разностей и формулы дифференцирования: Для дальнейших преобразований выведем формулу «суммирования по частям», для этого используем формулу произведения конечных разностей [3]: . (4) Рассмотрим подробнее слагаемое Преобразуем второе слагаемое в правой части (4) Здесь через S обозначен «оператор сдвига»: Из задачи (3) и определения получаем: следовательно, вариация принимает вид (6) Далее нам понадобится следующий дискретный аналог основной леммы вариационного исчисления. Лемма 1. Пусть - заданный набор чисел. Если при любом числовом наборе таком, что справедливо равенство , то . Доказательство проведем методом от противного. Предположим, что утверждение леммы неверно, т.е. найдется такое число , что . Построим числовой набор по следующему правилу Тогда , что несовместимо с условиями леммы. Мы пришли к противоречию, следовательно, наше допущение неверно и при всех Теперь, применяя лемму 1, приведём задачу (3) к разностному уравнению с заданными граничными условиями. Теорема 3. Если - точка экстремума вариационной задачи (3), то она является решением разностного уравнения Эйлера: = 0. (7) Доказательство. Из теоремы 1 следует, что , следовательно, сумма в формуле (6) равна 0 при любом числовом наборе таком, что . Следовательно, мы находимся в условиях применимости леммы 1, откуда следует выполнение равенства (7) при всех Напомним, что при и заданы граничные условия (см. задачу (3)), используя которые находим частное решение уравнения (7). 3. Дискретная задача Эйлера 3.1. Шарнирный стержень с равными звеньями Вернёмся к задаче 2 и, применяя полученные выше формулы, найдём её решение. Пусть E - модуль упругости, - наименьший момент инерции поперечных сечений стержня, - угол между кусками балки и в n-й точке c шарниром. Требуется определить наименьшую величину силы P, дающую продольный прогиб. По аналогии с задачей 1 решение задачи будет сводиться к поиску минимума потенциальной энергии после деформации. Потенциальная энергия изгиба определяется формулой: , где - момент прогиба в . При сжатии конца стержня на величину потенциальная энергия стержня уменьшается на . Если потенциальная энергия до деформации была равна нулю, то после деформации она выразится формулой . Так как , а для малых значений справедливо , то . В случае равновесия потенциальная энергия принимает минимальное значение. Поэтому решение задачи сводится к определению минимума функционала с граничными условиями , . Применяя к функционалу теорему 3, получаем: Следовательно, уравнение Эйлера имеет вид: (8) или , (9) где . Заметим, что . Решение (9), как известно по данным [3], зависит от дискриминанта характеристического уравнения , знак которого определяют три случая, приводящие к существенно различным решениям. Случай 1. Если , то, с учетом сделанного выше замечания, . Характеристическое уравнение имеет два вещественных различных корня и . По теореме Виета , , следовательно, . Будем считать, что , тогда , а решение уравнения (9) имеет вид . Найдем : . Подставляя первое граничное условие, имеем: . Так как , то , следовательно, . Из второго граничного условия имеем равенство , которое выполняется только при , что соответствует случаю тривиального положения равновесия. Итак, случай или не дает устойчивого прогиба. Случай 2. Пусть . Так как , то . Характеристическое уравнение принимает вид откуда следует , а решение (9) имеет вид . Так как , то из граничных условий получаем откуда следует, что . Значит, случай или также не дает устойчивого прогиба. Случай 3. Остается случай . В этом случае характеристическое уравнение имеет корни . Эти корни будут комплексно-сопряжёнными, т.е. . По теореме Виета , следовательно, , значит, , а решение уравнения (9) имеет вид где . Применяя граничное условие получим: . (10) Тогда из условия, что найдём минимальную силу: (11) Так как нас интересует ненулевой прогиб, то . Функция обращается в нуль при , так что (11) не будет выполнено, если , и наименьшее значение , при котором оно выполнено, - это . Из определения получаем , следовательно, Подставляя и , получаем, что наименьшее значение критической силы Эйлера равно , (12) так как при прогиб стержня равен нулю. Заметим, что найденная формула реализуется только при , т.е. при и при . Для полноты картины найдем вид кривой прогиба. Так как при малых значениях имеем , то будет задаваться кусочно-линейной функцией, угловой коэффициент каждого звена которой равен . Найдём первую прямую. Так как при имеем , то . Из граничного условия получаем, что , тогда . Рассуждая аналогичным образом, найдем уравнения всех остальных звеньев: (13) Так как звенья непрерывно стыкуются, из (13) получаем для следующие формулы: …, Найдём , используя равенство (10) и формулу для суммы косинусов [4]: , где . Для того, чтобы построить кривую прогиба, достаточно теперь отметить на плоскости точки с координатами и соединить их ломаными. 3.2. Шарнирный стержень со звеньями разной длины Рассмотрим обобщение задачи 2: пусть длины звеньев стержня не одинаковы (но рационально соизмеримы). Задача 3. Стержень длиною разделён шарнирами на участки произвольной целочисленной длины. Стержень опирается своими концами и подвержен давлению P (рис. 3). Шарниры, находящиеся на расстоянии от начала стержня, считаем идеальными, т.е. они дают стержню свободно изгибаться. При определённом значении P (критическая сила Эйлера) происходит продольный прогиб стержня. Требуется определить наименьшую величину силы P, дающую продольный изгиб. Рис. 3. Дискретный стержень c шарнирами и звеньями разной длины Решение. Пусть E - модуль упругости, - наименьший момент инерции поперечных сечений стержня, - угол между кусками балки и в n-й точке c шарниром. Требуется определить наименьшую величину силы P, дающую продольный прогиб. Чтобы использовать результаты раздела 3.1, можно применить следующий прием: разделим стержень на одинаковые звенья так, чтобы существующие шарниры попали в деления, а пустые деления заменим на виртуальные шарниры (на рис. 4 они обозначены цветом). Рис. 4. Стержень с виртуальными шарнирами Таким образом, решение задачи снова сводится к поиску минимума функционала, выражающего потенциальную энергию после деформации: Для этого функционала были найдены экстремали (формула (10)) и (формула (12)). Отличие результата только в том, что равенство будет интересовать нас не при всех n, a лишь при Далее, так как виртуальные шарниры не дают прогибаться стержню и по сути являются точками отрезка, составляющего кусок стержня, то значения угла в этих точках будут равны между собой и равны значению угла реального шарнира, находящегося в начале отрезка. Найдем вид кривой прогиба. Рассуждения, аналогичные тем, что были проведены при решении задачи 2, приводят к следующим уравнениям: Из этих равенств следует, что , а так как формула для уже найдена выше, то можно легко построить кривую прогиба, которая будет кусочно-линейной.

About the authors

I. A Aksenenko

Perm National Research Polytechnic University

References

  1. Ванько В.И. Очерки об устойчивости элементов конструкций. - М.: Изд. МГТУ им. Н.Э. Баумана, 2015.
  2. Эльсгольц Л.Э. Дифференциальные уравнения и вариационное исчисление. - М.: Эдиториал УРСС, 2000.
  3. Романко В.К. Разностные уравнения: учебное пособие. - М.: БИНОМ. Лаборатория знаний, 2006.
  4. Цыпкин А.Г., Цыпкин Г.Г. Математические формулы. Алгебра. Геометрия. Математический анализ: справочник. - М.: Наука, 1985.

Statistics

Views

Abstract - 97

PDF (Russian) - 73

Refbacks

  • There are currently no refbacks.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies