Геометрия изгиба: Вторая производная и теория кривизны плоских кривых

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

Геометрический смысл второй производной и скорость изменения наклона касательной

Геометрический смысл второй производной и скорость изменения наклона касательной

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

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

От координаты к уклону: препарирование первой производной

Для начала зафиксируем базовый механизм работы дифференциального исчисления. Пусть задана непрерывная и гладкая функция y=f(x)y = f(x), где yy — зависимая переменная (ордината), а f(x)f(x) — функция от независимой переменной xx (абсциссы). Значение f(x)f(x) (где f(x)f(x) — значение функции в точке xx) определяет высоту графика.

Первая производная f(x)f'(x) (где f(x)f'(x) — первая производная функции) отвечает на вопрос: «С какой скоростью меняется высота при смещении вправо?». Геометрически значение f(x0)f'(x_0) (где f(x0)f'(x_0) — значение первой производной в точке x0x_0) равно тангенсу угла наклона касательной к графику в этой точке:

f(x0)=tan(α)f'(x_0) = \tan(\alpha)

В этой формуле f(x0)f'(x_0) обозначает значение первой производной в точке x0x_0, tan\tan — функцию тангенса, а α\alpha — угол между касательной и положительным направлением оси абсцисс.

Если первая производная больше нуля (f(x)>0f'(x) > 0, где f(x)f'(x) — первая производная), график идет на подъем, если меньше нуля (f(x)<0f'(x) < 0, где f(x)f'(x) — первая производная) — на спуск. Однако знание одного лишь уклона не дает информации о форме самой горы. Подъем может становиться все более крутым, может выполаживаться, а может оставаться неизменным (как на идеальном пандусе). Чтобы описать характер этого подъема, требуется инструмент, измеряющий скорость изменения самой скорости.

Вторая производная как «руль» функции

Вторая производная f(x)f''(x) (где f(x)f''(x) — вторая производная функции) вычисляется как производная от первой производной:

f(x)=ddx[f(x)]f''(x) = \frac{d}{dx}[f'(x)]

Здесь f(x)f''(x) — вторая производная функции, ddx\frac{d}{dx} — оператор дифференцирования по переменной xx, а f(x)f'(x) — первая производная.

Если первая производная — это уклон траектории, то вторая производная показывает, как быстро этот уклон меняется при продвижении вдоль оси xx (где xx — независимая переменная). В терминах автомобильного движения: функция f(x)f(x) (где f(x)f(x) — функция траектории) — это дорога, f(x)f'(x) (где f(x)f'(x) — первая производная) — положение компаса, указывающего направление капота, а f(x)f''(x) (где f(x)f''(x) — вторая производная) — это скорость вращения рулевого колеса.

Рассмотрим механику этого вращения детально. Угол наклона касательной α\alpha (где α\alpha — угол наклона касательной) может меняться по-разному. Если вторая производная больше нуля (f(x)>0f''(x) > 0, где f(x)f''(x) — вторая производная) на некотором интервале, это означает, что первая производная f(x)f'(x) (где f(x)f'(x) — первая производная) возрастает. Тангенс угла наклона становится все больше.

  • Если мы двигались вниз (отрицательный уклон), спуск начинает выполаживаться (уклон растет от 5-5 до 1-1, стремясь к нулю).
  • Если мы двигались вверх (положительный уклон), подъем становится все круче (уклон растет от 11 до 55).

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

Если же вторая производная меньше нуля (f(x)<0f''(x) < 0, где f(x)f''(x) — вторая производная), первая производная убывает. Касательная поворачивается по часовой стрелке. Подъем сменяется выполаживанием, а спуск становится еще более обрывистым. Кривая загибается вниз, образуя купол. Это выпуклость вверх.

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

  1. Если на интервале вторая производная больше нуля (f(x)>0f''(x) > 0, где f(x)f''(x) — вторая производная), то график функции лежит строго выше любой своей касательной, проведенной на этом интервале (не считая самой точки касания).
  2. Если на интервале вторая производная меньше нуля (f(x)<0f''(x) < 0, где f(x)f''(x) — вторая производная), то график функции лежит строго ниже любой своей касательной.

Анализ базовых моделей изгиба

Разберем, как вторая производная формирует геометрию на классических примерах.

Пример 1: Квадратичная парабола

Рассмотрим функцию y=x2y = x^2, где yy — значение функции, а xx — независимая переменная. Ее первая производная y=2xy' = 2x (где yy' — первая производная функции, а xx — независимая переменная) показывает, что уклон линейно зависит от координаты. В точке x=2x = -2 (где xx — абсцисса) уклон равен 4-4 (крутой спуск), в нуле он равен 00 (дно), в точке x=2x = 2 (где xx — абсцисса) уклон равен 44 (крутой подъем).

Вторая производная:

y=2y'' = 2

В этой формуле yy'' — вторая производная функции y=x2y = x^2.

Она положительна и постоянна на всей числовой оси. Это означает, что скорость приращения уклона неизменна. На каждый шаг вдоль оси xx (где xx — независимая переменная) касательная увеличивает свой тангенс ровно на 2 единицы. График везде лежит выше своих касательных, образуя идеальную чашу. Постоянство второй производной делает параболу кривой с равномерно распределенным «усилием изгиба» относительно оси абсцисс.

Пример 2: Кубическая парабола

Функция y=x3y = x^3, где yy — значение функции, а xx — независимая переменная. Первая производная y=3x2y' = 3x^2 (где yy' — первая производная, а xx — независимая переменная). Заметим, что первая производная неотрицательна (y0y' \geq 0, где yy' — первая производная) при любых xx (где xx — независимая переменная). График никогда не идет вниз, он только растет (или на мгновение замирает в нуле).

Вторая производная:

y=6xy'' = 6x

Здесь yy'' — вторая производная функции y=x3y = x^3, а xx — независимая переменная.

Знак второй производной зависит от координаты xx (где xx — независимая переменная):

  • При x<0x < 0 (где xx — независимая переменная) имеем y<0y'' < 0 (где yy'' — вторая производная). График выпуклый вверх (купол), он лежит ниже своих касательных.
  • При x>0x > 0 (где xx — независимая переменная) имеем y>0y'' > 0 (где yy'' — вторая производная). График выпуклый вниз (чаша), он лежит выше касательных.

Точка x=0x = 0 (где xx — абсцисса), в которой вторая производная равна нулю и меняет свой знак, является точкой перегиба. В этой точке руль, который мы крутили вправо, мгновенно проходит нейтральное положение и начинает крутиться влево. Касательная в этой точке протыкает график насквозь.

Главная иллюзия: почему вторая производная — это не кривизна

На этом этапе возникает сильный соблазн поставить знак равенства между модулем второй производной f(x)|f''(x)| (где f(x)|f''(x)| — модуль второй производной) и геометрической кривизной линии. Кажется логичным: чем больше вторая производная, тем быстрее меняется наклон, а значит, тем сильнее изогнута линия.

Это одно из самых опасных заблуждений в дифференциальной геометрии. Вторая производная не является объективной мерой кривизны.

Чтобы понять причину, нужно вспомнить, что именно дифференцирует вторая производная f(x)f''(x) (где f(x)f''(x) — вторая производная). Она измеряет скорость изменения тангенса угла наклона по отношению к горизонтальному смещению xx (где xx — независимая переменная).

f(x)=d(tanα)dxf''(x) = \frac{d(\tan \alpha)}{dx}

В этой формуле f(x)f''(x) — вторая производная, tanα\tan \alpha — тангенс угла наклона касательной, а xx — независимая переменная.

Но истинная геометрическая кривизна (насколько сильно изогнут кусок проволоки в пространстве) не должна зависеть от того, под каким углом мы на эту проволоку смотрим или как располагаем оси координат. Истинная кривизна — это скорость изменения самого угла α\alpha (где α\alpha — угол наклона касательной) по отношению к пройденному пути ss (где ss — длина дуги, то есть пройденный путь вдоль кривой).

Посмотрим, как связаны тангенс угла наклона tan(α)\tan(\alpha) (где tan(α)\tan(\alpha) — тангенс угла, а α\alpha — сам угол) и сам угол α\alpha. Из правил дифференцирования сложной функции:

d(tanα)dx=1cos2αdαdx\frac{d(\tan \alpha)}{dx} = \frac{1}{\cos^2 \alpha} \cdot \frac{d\alpha}{dx}

Здесь tanα\tan \alpha — тангенс угла наклона касательной, cosα\cos \alpha — косинус этого угла, xx — абсцисса, а dαdx\frac{d\alpha}{dx} — производная угла по координате xx.

Множитель 1cos2α\frac{1}{\cos^2 \alpha} (где cosα\cos \alpha — косинус угла наклона) играет роковую роль. Если касательная становится очень крутой (угол α\alpha, где α\alpha — угол наклона касательной, приближается к 9090^\circ или к значению π2\frac{\pi}{2}, где π\pi — математическая константа пи), косинус стремится к нулю, а дробь 1cos2α\frac{1}{\cos^2 \alpha} (где cosα\cos \alpha — косинус угла наклона) улетает в бесконечность.

Это означает, что для круто идущих вверх графиков вторая производная выдаст гигантские значения просто из-за математического свойства тангенса, даже если сама линия почти прямая!

Доказательство через окружность

Проверим этот парадокс на идеальном геометрическом объекте — окружности. Окружность имеет одинаковый изгиб в любой своей точке. Ее истинная кривизна должна быть константой.

Возьмем верхнюю полуокружность радиуса RR (где RR — радиус окружности) с центром в начале координат. Ее уравнение:

y=R2x2y = \sqrt{R^2 - x^2}

В этой формуле yy — ордината точки на верхней полуокружности, RR — радиус окружности, а xx — ее абсцисса.

Найдем первую производную (уклон):

y=xR2x2y' = \frac{-x}{\sqrt{R^2 - x^2}}

Здесь yy' — первая производная функции полуокружности, xx — абсцисса, а RR — радиус.

Теперь вычислим вторую производную, применив правило дифференцирования дроби:

y=1R2x2(x)xR2x2(R2x2)2y'' = \frac{-1 \cdot \sqrt{R^2 - x^2} - (-x) \cdot \frac{-x}{\sqrt{R^2 - x^2}}}{(\sqrt{R^2 - x^2})^2}

В этой формуле yy'' — вторая производная функции полуокружности, xx — абсцисса, а RR — радиус.

После приведения к общему знаменателю в числителе получим:

y=(R2x2)x2(R2x2)3/2=R2(R2x2)3/2y'' = \frac{-(R^2 - x^2) - x^2}{(R^2 - x^2)^{3/2}} = \frac{-R^2}{(R^2 - x^2)^{3/2}}

Здесь yy'' — упрощенное выражение для второй производной, xx — абсцисса, а RR — радиус окружности.

Проанализируем результат. Вторая производная yy'' (где yy'' — вторая производная) зависит от xx (где xx — абсцисса). В самой верхней точке окружности при x=0x = 0 (где xx — абсцисса) касательная горизонтальна. Здесь:

y=R2R3=1Ry'' = \frac{-R^2}{R^3} = -\frac{1}{R}

где yy'' — вторая производная, а RR — радиус окружности.

Но что происходит, когда мы сдвигаемся к краю окружности, например, когда xx стремится к RR (xRx \to R, где xx — абсцисса, а RR — радиус)? Знаменатель (R2x2)3/2(R^2 - x^2)^{3/2} (где RR — радиус, а xx — абсцисса) стремится к нулю. Следовательно, вторая производная стремится к минус бесконечности: yy'' \to -\infty (где yy'' — вторая производная).

Окружность изгибается везде одинаково, но ее вторая производная меняется от 1R-\frac{1}{R} до -\infty (где RR — радиус окружности). Этот пример неопровержимо доказывает: вторая производная искажает восприятие формы при больших углах наклона. Она хорошо описывает изгиб только на пологих участках (где косинус угла наклона касательной близок к единице, то есть cosα1\cos \alpha \approx 1, а пройденный путь ss примерно равен смещению по горизонтали xx, где ss — длина дуги, а xx — абсцисса). Для точного описания изгиба в любых условиях потребуется сконструировать новый инвариантный аппарат, свободный от привязки к осям координат.

Нулевая вторая производная и физика рывка

Вернемся к инженерной задаче, упомянутой в начале. Почему переход от прямой к дуге окружности вызывает проблемы?

На прямом участке пути уклон не меняется, поэтому вторая производная равна нулю (f(x)=0f''(x) = 0, где f(x)f''(x) — вторая производная). На участке круговой дуги вторая производная отлична от нуля (и, как мы выяснили, непостоянна относительно xx, но для пологих поворотов железной дороги она держится на уровне примерно 1R\approx -\frac{1}{R}, где RR — радиус дуги).

В физике первая производная координаты по времени — это скорость, вторая — ускорение, а третья — рывок (jerk). При движении вдоль трассы с постоянной скоростью поперечное (центробежное) ускорение пропорционально кривизне пути (которая при малых углах наклона близка ко второй производной f(x)f''(x), где f(x)f''(x) — вторая производная функции траектории). Если вторая производная f(x)f''(x) меняется скачком, то и центробежное ускорение меняется скачком. Это означает, что его производная по времени — поперечный рывок — в точке стыка стремится к бесконечности (математически выражается через дельта-функцию Дирака). Бесконечный рывок физически ощущается как жесткий удар, способный сбросить состав с рельсов или разрушить подвеску автомобиля.

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

Понимание геометрического смысла второй производной f(x)f''(x) (где f(x)f''(x) — вторая производная) как скорости изменения уклона — это мост между абстрактным дифференцированием и реальной геометрией кривых. Вторая производная безошибочно указывает, в какую сторону выгибается линия и где находятся точки перелома тенденций. Однако ее зависимость от крутизны самого графика ставит предел ее применимости как универсальной меры изгиба, открывая путь к более глубокому понятию — фундаментальной кривизне.

Исследование функций на выпуклость и вогнутость: критерии и интервалы

Тросы подвесного моста Золотые Ворота и каменные своды древнеримских акведуков выполняют одну и ту же задачу — удерживают колоссальный вес над пустотой. Однако с точки зрения математики они являются абсолютными противоположностями. Форма провисающего троса — это кривая, строго выпуклая вниз, где каждый элемент работает на растяжение. Форма каменной арки — кривая, строго выпуклая вверх, где материал испытывает только сжатие. Знак второй производной функции, описывающей эти конструкции, не просто определяет геометрический силуэт, он диктует инженерам выбор материала: сталь для положительной второй производной и камень или бетон — для отрицательной.

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

Достаточные условия выпуклости и вогнутости

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

В математическом анализе поведение функции на интервале определяется через достаточные условия. Пусть задана функция f(x)f(x), которая имеет первую и вторую производные на всем интервале (a,b)(a, b).

Если для любой точки xx из интервала (a,b)(a, b) выполняется неравенство f(x)>0f''(x) > 0, то график функции на этом интервале имеет строгую выпуклость вниз. Геометрически это означает, что любая хорда, соединяющая две точки графика на этом участке, будет лежать строго выше самой дуги кривой.

Если же для любой точки xx из интервала (a,b)(a, b) справедливо неравенство f(x)<0f''(x) < 0, то график функции строго выпуклый вверх. В этом случае любая хорда на данном отрезке окажется строго ниже дуги графика.

Здесь f(x)f(x) — значение функции в точке, xx — независимая переменная, а f(x)f''(x) — значение второй производной в этой же точке.

Важно понимать разницу между строгой и нестрогой выпуклостью. Если мы ослабим условие до f(x)0f''(x) \geq 0, мы получим определение просто выпуклой вниз функции. В таком случае график может содержать абсолютно плоские, прямолинейные участки, где вторая производная тождественно равна нулю на целом отрезке. Однако в дифференциальной геометрии кривых наибольший интерес представляют именно интервалы строгой выпуклости, где кривизна постоянно присутствует, пусть и в бесконечно малых значениях.

Алгоритм поиска интервалов: критические точки второго рода

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

Такие границы называются критическими точками второго рода. К ним относятся:

  1. Точки, в которых вторая производная равна нулю: f(x)=0f''(x) = 0.
  2. Точки, в которых вторая производная не существует (например, точки разрыва самой функции, вертикальные асимптоты или точки возврата, где график имеет острый "излом").

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

Влияние асимптот на смену выпуклости

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

Рассмотрим классический пример рациональной функции, описывающей сложное электростатическое поле или распределение скоростей в гидродинамике:

f(x)=xx21f(x) = \frac{x}{x^2 - 1}

Здесь f(x)f(x) — исследуемая функция, а xx — координата. Область определения этой функции прерывается в двух точках: x=1x = -1 и x=1x = 1. В этих точках знаменатель обращается в нуль, образуя вертикальные асимптоты.

Найдем первую производную, используя правило дифференцирования дроби:

f(x)=1(x21)x2x(x21)2=x21(x21)2f'(x) = \frac{1 \cdot (x^2 - 1) - x \cdot 2x}{(x^2 - 1)^2} = \frac{-x^2 - 1}{(x^2 - 1)^2}

В этой формуле числитель всегда отрицателен, а знаменатель всегда положителен (в области определения). Это означает, что f(x)<0f'(x) < 0 везде, то есть функция постоянно убывает. Но убывать можно по-разному: ускоряясь (выпуклость вверх) или замедляясь (выпуклость вниз). Чтобы выяснить это, найдем вторую производную:

f(x)=(2x)(x21)2(x21)2(x21)2x(x21)4f''(x) = \frac{(-2x)(x^2 - 1)^2 - (-x^2 - 1) \cdot 2(x^2 - 1) \cdot 2x}{(x^2 - 1)^4}

Сократив дробь на общий множитель (x21)(x^2 - 1), получим:

f(x)=2x(x21)+4x(x2+1)(x21)3=2x3+6x(x21)3=2x(x2+3)(x21)3f''(x) = \frac{-2x(x^2 - 1) + 4x(x^2 + 1)}{(x^2 - 1)^3} = \frac{2x^3 + 6x}{(x^2 - 1)^3} = \frac{2x(x^2 + 3)}{(x^2 - 1)^3}

Теперь найдем критические точки второго рода. Приравняем числитель к нулю: 2x(x2+3)=02x(x^2 + 3) = 0. Поскольку скобка (x2+3)(x^2 + 3) всегда строго больше нуля, единственным решением является x=0x = 0. Знаменатель обращается в нуль в точках x=1x = -1 и x=1x = 1.

Нанесем точки 1-1, 00 и 11 на числовую ось. Они разбивают область определения на четыре интервала. Проверим знак f(x)f''(x) на каждом из них:

  • На интервале (,1)(-\infty, -1) возьмем пробную точку x=2x = -2. Числитель: 2(2)(4+3)=28<02(-2)(4 + 3) = -28 < 0. Знаменатель: ((2)21)3=27>0((-2)^2 - 1)^3 = 27 > 0. Отношение отрицательное, f(x)<0f''(x) < 0. График выпуклый вверх.
  • На интервале (1,0)(-1, 0) возьмем x=0.5x = -0.5. Числитель отрицательный, знаменатель ((0.5)21)3((-0.5)^2 - 1)^3 тоже отрицательный. Отношение положительное, f(x)>0f''(x) > 0. График выпуклый вниз.
  • На интервале (0,1)(0, 1) возьмем x=0.5x = 0.5. Числитель положительный, знаменатель отрицательный. f(x)<0f''(x) < 0. График выпуклый вверх.
  • На интервале (1,)(1, \infty) возьмем x=2x = 2. Числитель и знаменатель положительны. f(x)>0f''(x) > 0. График выпуклый вниз.

Этот пример блестяще демонстрирует, что смена выпуклости происходит не только в точке x=0x = 0 (где график физически пересекает ось и плавно меняет кривизну), но и при переходе через вертикальные асимптоты x=1x = -1 и x=1x = 1, хотя самой функции в этих точках не существует. Кривая уходит в минус бесконечность, будучи выпуклой вверх, а возвращается из плюс бесконечности уже выпуклой вниз.

Нулевая вторая производная без смены выпуклости

Существует опасная ловушка при исследовании функций: равенство второй производной нулю в некоторой точке не гарантирует, что кривая меняет направление выпуклости. Условие f(x)=0f''(x) = 0 является необходимым для смены формы (если функция непрерывна), но не достаточным.

Рассмотрим функцию f(x)=x4f(x) = x^4. Ее первая производная: f(x)=4x3f'(x) = 4x^3. Вторая производная: f(x)=12x2f''(x) = 12x^2.

Найдем критические точки второго рода: 12x2=012x^2 = 0, откуда x=0x = 0. Однако, если мы проверим знаки на интервалах до и после нуля, мы увидим парадоксальную картину. При x<0x < 0 (например, x=1x = -1), f(1)=12(1)2=12>0f''(-1) = 12(-1)^2 = 12 > 0. При x>0x > 0 (например, x=1x = 1), f(1)=12(1)2=12>0f''(1) = 12(1)^2 = 12 > 0.

Вторая производная положительна по обе стороны от нуля. График функции f(x)=x4f(x) = x^4 строго выпуклый вниз на всей числовой оси, несмотря на то, что в точке x=0x = 0 его вторая производная обнуляется.

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

Глобальная выпуклость: экспонента и логарифм

Некоторые фундаментальные функции сохраняют характер своей выпуклости на всей области определения. Их форма задает базовые законы роста и затухания в природе.

Экспоненциальная функция f(x)=exf(x) = e^x. Ее производные уникальны: f(x)=exf'(x) = e^x и f(x)=exf''(x) = e^x. Поскольку число ee в любой степени строго больше нуля, f(x)>0f''(x) > 0 для любого действительного xx. Экспонента глобально выпукла вниз. Геометрически это означает, что скорость ее роста постоянно увеличивается. Какую бы касательную мы ни провели к графику exe^x, кривая всегда будет отрываться от нее и уходить вверх. Именно эта математическая особенность делает экспоненциальный рост таким "взрывным" — график неизбежно обгонит любую, даже самую крутую прямую линию.

Натуральный логарифм f(x)=ln(x)f(x) = \ln(x), определенный для x>0x > 0. Первая производная: f(x)=1xf'(x) = \frac{1}{x}. Вторая производная: f(x)=1x2f''(x) = -\frac{1}{x^2}. Для любого допустимого xx, квадрат x2>0x^2 > 0, следовательно, дробь с минусом всегда отрицательна: f(x)<0f''(x) < 0. Логарифмическая кривая глобально выпукла вверх. Ее рост постоянно замедляется, и кривая всегда будет прижиматься к земле, в конечном итоге оказываясь ниже любой своей касательной.

Инженерное приложение: эпюры изгибающих моментов

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

M(x)=EIy(x)M(x) = E \cdot I \cdot y''(x)

В этой формуле M(x)M(x) — изгибающий момент во внутреннем сечении балки на расстоянии xx от края, EE — модуль Юнга (характеризует жесткость самого материала, например, стали или бетона), II — момент инерции сечения (зависит от формы балки, например, двутавр или прямоугольник), а y(x)y''(x) — вторая производная функции прогиба балки.

Поскольку величины EE и II всегда строго положительны, знак изгибающего момента M(x)M(x) полностью совпадает со знаком второй производной y(x)y''(x). Это совпадение имеет критическое значение для безопасности конструкций.

Если на определенном интервале y(x)>0y''(x) > 0, балка выпукла вниз (провисает "чашей"). В этом состоянии нижние слои материала растягиваются, а верхние — сжимаются. Если y(x)<0y''(x) < 0, балка выпукла вверх (выгибается "куполом", что часто происходит над промежуточными опорами многопролетных мостов). В этом случае растягиваются верхние слои, а сжимаются нижние.

Бетон — материал, который прекрасно выдерживает колоссальные сжимающие нагрузки, но крайне хрупок при растяжении. Чтобы бетонная балка не треснула под собственным весом, в зону растяжения закладывают стальную арматуру. Таким образом, математическая задача поиска интервалов, где y(x)>0y''(x) > 0 и y(x)<0y''(x) < 0, превращается в прямой инженерный чертеж: на участках с положительной второй производной стальные стержни укладывают в нижнюю часть балки, а на участках с отрицательной — в верхнюю. Ошибка в определении знака второй производной приведет к тому, что арматура окажется в сжатой зоне, где она бесполезна, а неармированный бетон в зоне растяжения мгновенно разрушится.

Анализ второй производной позволяет не просто констатировать форму кривой в статике. Разбиение области определения на интервалы выпуклости и вогнутости дает инструмент для предсказания поведения функции. Мы видим, как кривая реагирует на приближение к асимптотам, как она может "обманывать" нулевой второй производной, сохраняя форму, и как строгие алгебраические неравенства материализуются в расположении стальных каркасов внутри бетонных монолитов.

Точки перегиба: необходимые и достаточные условия существования и алгоритмы поиска

Точки перегиба: необходимые и достаточные условия существования и алгоритмы поиска

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

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

Геометрическая анатомия перегиба

Строгое определение гласит: точка M(x0,f(x0))M(x_0, f(x_0)) называется точкой перегиба графика непрерывной функции f(x)f(x), если при переходе через абсциссу x0x_0 функция меняет направление выпуклости.

Из этого определения вытекает важнейшее геометрическое свойство. Мы знаем, что на участке выпуклости вниз график лежит выше своих касательных, а на участке выпуклости вверх — ниже. Если в точке x0x_0 выпуклость меняется, значит, слева от точки касания график должен находиться по одну сторону от касательной, а справа — по другую.

Единственный способ для непрерывной кривой удовлетворить этому условию — пересечь прямую касательной насквозь. В точке перегиба касательная перестает быть просто линией, к которой график «прижимается» с одной стороны. Она становится секущей, локально разрезая кривую на две части.

Необходимое условие: где искать перегиб

Поиск точек перегиба начинается с анализа второй производной f(x)f''(x). Если кривая в точке x0x_0 имеет перегиб, и в этой точке существует вторая производная, то она обязана быть равной нулю: f(x0)=0f''(x_0) = 0.

Логика этого утверждения опирается на теорему Больцано-Коши о промежуточном значении. Если функция меняет выпуклость, значит, вторая производная меняет знак (например, с плюса на минус). Если при этом f(x)f''(x) непрерывна, она не может перепрыгнуть из положительной зоны в отрицательную, не пересекая ось абсцисс.

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

Рассмотрим функцию f(x)=xxf(x) = x \cdot |x|. Раскрыв модуль, получаем кусочно-заданную функцию: x2x^2 при x0x \geq 0 и x2-x^2 при x<0x < 0. Первая производная f(x)=2xf'(x) = 2|x| существует везде, включая ноль (она равна нулю, касательная горизонтальна). Но вторая производная f(x)f''(x) равна 22 при x>0x > 0 и 2-2 при x<0x < 0. В самой точке x=0x = 0 происходит скачок, f(0)f''(0) не существует. Тем не менее, знак второй производной меняется с минуса на плюс, выпуклость меняется, и начало координат является полноценной точкой перегиба.

Поэтому кандидатами на точки перегиба (критическими точками второго рода) являются все значения xx, при которых f(x)=0f''(x) = 0 или f(x)f''(x) не существует.

Достаточные условия: как отсеять ложные срабатывания

Необходимое условие дает нам список подозреваемых. Но не каждая точка, где f(x)=0f''(x) = 0, является точкой перегиба. Чтобы подтвердить статус точки, требуется выполнить одно из достаточных условий.

Первый признак: смена знака

Самый универсальный метод — проверить знак f(x)f''(x) в окрестности критической точки. Если при переходе через x0x_0 вторая производная меняет знак, точка перегиба существует.

Алгоритм проверки (метод интервалов) прямолинеен:

  1. Найти область определения функции.
  2. Вычислить f(x)f''(x).
  3. Найти нули второй производной и точки ее несуществования.
  4. Разбить область определения этими точками на интервалы.
  5. Вычислить знак f(x)f''(x) внутри каждого интервала. Если знаки по бокам от критической точки разные — это перегиб. Если одинаковые — перегиба нет.

Второй признак: тест третьей производной

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

Если в точке x0x_0 вторая производная равна нулю (f(x0)=0f''(x_0) = 0), но третья производная отлична от нуля (f(x0)0f'''(x_0) \neq 0), то x0x_0 гарантированно является точкой перегиба.

Глубинный смысл этого правила раскрывается через разложение функции в ряд Тейлора в окрестности точки x0x_0:

f(x)f(x0)+f(x0)(xx0)+f(x0)2(xx0)2+f(x0)6(xx0)3f(x) \approx f(x_0) + f'(x_0)(x - x_0) + \frac{f''(x_0)}{2}(x - x_0)^2 + \frac{f'''(x_0)}{6}(x - x_0)^3

Первые два слагаемых f(x0)+f(x0)(xx0)f(x_0) + f'(x_0)(x - x_0) — это уравнение касательной линии к графику в точке x0x_0. Обозначим его как yкасy_{кас}. Разность между функцией и касательной показывает, как кривая отклоняется от прямой:

f(x)yкасf(x0)2(xx0)2+f(x0)6(xx0)3f(x) - y_{кас} \approx \frac{f''(x_0)}{2}(x - x_0)^2 + \frac{f'''(x_0)}{6}(x - x_0)^3

Если f(x0)0f''(x_0) \neq 0, то поведение кривой диктуется квадратичным членом (xx0)2(x - x_0)^2. Квадрат любого числа положителен, поэтому кривая отходит от касательной в одну и ту же сторону (вверх или вниз) независимо от того, двигаемся мы вправо (x>x0x > x_0) или влево (x<x0x < x_0). Это стандартная выпуклость.

Но если f(x0)=0f''(x_0) = 0, квадратичный член исчезает. Главным становится кубический член (xx0)3(x - x_0)^3.

Кубическая парабола меняет знак. При x>x0x > x_0 выражение (xx0)3(x - x_0)^3 положительно, а при x<x0x < x_0 — отрицательно. Это аналитически доказывает, что функция переходит с одной стороны касательной на другую. Если же и третья производная равна нулю, придется смотреть на четвертую (которая, будучи четной степенью, снова не даст перегиба, если она не равна нулю). Перегиб возникает только тогда, когда первая ненулевая производная после первой имеет нечетный порядок.

Аномальные случаи: перегиб с вертикальной касательной

Интуитивно мы представляем график в точке перегиба пологим. Однако математика допускает ситуации, когда в точке перегиба кривая становится абсолютно вертикальной.

Рассмотрим функцию f(x)=x3f(x) = \sqrt[3]{x}. Найдем первую производную:

f(x)=13x2/3=13x23f'(x) = \frac{1}{3}x^{-2/3} = \frac{1}{3\sqrt[3]{x^2}}

При x=0x = 0 знаменатель обращается в ноль, производная стремится к бесконечности. Это означает, что касательная совпадает с осью ординат (вертикальна), и функция в этой точке не дифференцируема в классическом смысле.

Теперь найдем вторую производную:

f(x)=29x5/3=29x53f''(x) = -\frac{2}{9}x^{-5/3} = -\frac{2}{9\sqrt[3]{x^5}}

Вторая производная при x=0x=0 также не существует. Применим первый достаточный признак и проверим знаки в окрестности нуля. При x<0x < 0 (например, x=1x = -1), значение x5x^5 отрицательно, корень из него отрицателен, и два минуса дают плюс: f(1)>0f''(-1) > 0 (выпуклость вниз). При x>0x > 0 (например, x=1x = 1), знаменатель положителен, перед дробью стоит минус: f(1)<0f''(1) < 0 (выпуклость вверх).

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

Физический и прикладной смысл точек перегиба

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

Кинематика: момент нулевого ускорения

Если функция s(t)s(t) описывает путь материальной точки в зависимости от времени, то первая производная s(t)s'(t) — это скорость v(t)v(t), а вторая производная s(t)s''(t) — ускорение a(t)a(t).

Точка перегиба графика пути — это момент времени, когда s(t)=0s''(t) = 0 и меняет знак. Физически это означает, что ускорение переходит от положительных значений к отрицательным (или наоборот). То есть объект перестает разгоняться и начинает тормозить. По теореме Ферма для первой производной, точка, где производная скорости (ускорение) равна нулю, является точкой экстремума скорости. Таким образом, точка перегиба на графике пути всегда соответствует моменту максимальной (или минимальной) скорости движения.

Динамика популяций: пик логистической кривой

В биологии, социологии и экономике рост систем в условиях ограниченных ресурсов описывается логистическим уравнением (моделью Ферхюльста). График численности популяции N(t)N(t) имеет характерную S-образную форму.

Аналитическое выражение логистической кривой:

N(t)=K1+er(tt0)N(t) = \frac{K}{1 + e^{-r(t-t_0)}}

где KK — емкость среды (максимально возможная численность), rr — коэффициент скорости роста.

На начальном этапе популяция мала, ресурсов в избытке, и рост идет по экспоненте (выпуклость вниз, вторая производная положительна — скорость роста увеличивается). Однако по мере приближения к пределу KK, нехватка еды или территории начинает замедлять рост. Скорость роста падает (выпуклость вверх, вторая производная отрицательна).

Точка перегиба находится точно посередине, при численности N=K/2N = K/2. Это критический момент в жизни популяции: именно здесь скорость прироста достигает своего абсолютного максимума, после чего система переходит в режим торможения. В эпидемиологии точка перегиба на графике числа зараженных означает прохождение пика заболеваемости — момент, когда новые случаи еще появляются, но их количество за день начинает снижаться.

Природные формы и инженерия

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

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

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

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

Фундаментальное понятие кривизны и радиуса кривизны как меры отклонения от прямой

Фундаментальное понятие кривизны и радиуса кривизны как меры отклонения от прямой

Водитель автомобиля, движущегося с постоянной скоростью по извилистой трассе, физически ощущает геометрию дороги через рулевое колесо. Чем круче поворот, тем сильнее приходится выкручивать руль. Если дорога прямая, руль находится в нейтральном положении. Если трасса переходит в идеальное круговое кольцо, водитель поворачивает руль на определенный угол и фиксирует его. В этот момент угол поворота рулевого колеса является точным механическим аналогом фундаментального математического понятия — кривизны траектории. Вторая производная, как мы выяснили ранее, зависит от выбора системы координат и искажается при больших углах наклона графика. Настоящая, «внутренняя» геометрия линии требует инструмента, который описывает форму кривой безотносительно осей XX и YY.

Отказ от декартовых координат: естественная параметризация

Чтобы измерить истинную изогнутость линии, необходимо отказаться от привязки к внешним осям координат. В декартовой системе мы оцениваем, как меняется функция yy при изменении аргумента xx. Но для самой кривой нет никаких xx и yy. Если нарисовать линию на листе бумаги, а затем повернуть лист, форма линии останется неизменной, хотя ее уравнение y=f(x)y = f(x) и значения всех производных изменятся кардинально.

Математика решает эту проблему переходом к внутренней системе отсчета, которая "живет" на самой кривой. Эта концепция называется естественной параметризацией.

Выберем на кривой произвольную точку отсчета M0M_0. Положение любой другой точки MM на этой линии можно однозначно задать одним числом ss — длиной дуги от M0M_0 до MM. Переменная ss становится естественным параметром. Движение вдоль кривой означает непрерывное изменение ss. Теперь любое геометрическое свойство линии можно изучать как функцию от пройденного пути ss, а не от абстрактной координаты xx.

Второе ключевое понятие — угол наклона касательной α\alpha. Проведем в точке MM касательную прямую и измерим угол α\alpha между этой касательной и неким фиксированным направлением (например, горизонталью). По мере того как точка движется по кривой, касательная вращается, и угол α\alpha меняется.

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

Средняя кривизна дуги

Рассмотрим участок кривой между точками M1M_1 и M2M_2. Пусть длина дуги между ними равна Δs\Delta s. В точке M1M_1 касательная образует угол α1\alpha_1, а в точке M2M_2 — угол α2\alpha_2.

Разность Δα=α2α1\Delta \alpha = \alpha_2 - \alpha_1 называется углом смежности. Этот угол показывает, на какую величину повернулась касательная при переходе от начала дуги к ее концу.

Средней кривизной дуги M1M2M_1M_2 называется отношение угла смежности к длине этой дуги:

Kcp=ΔαΔsK_{cp} = \left| \frac{\Delta \alpha}{\Delta s} \right|

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

Смысл средней кривизны предельно прост: это средняя скорость вращения касательной на единицу длины пути. Если Δα\Delta \alpha измеряется в радианах, а Δs\Delta s в метрах, то размерность кривизны — рад/м (или просто 1/м1/\text{м}). Значение Kcp=0.1K_{cp} = 0.1 означает, что в среднем на каждом метре пути касательная поворачивается на 0.10.1 радиана.

Строгое определение кривизны в точке

Средняя кривизна описывает участок целиком, но реальные кривые изгибаются неравномерно. Чтобы получить характеристику формы в конкретной точке MM, необходимо стянуть дугу в эту точку, устремив длину пути Δs\Delta s к нулю.

Кривизной кривой κ\kappa (каппа) в точке MM называется предел средней кривизны при Δs0\Delta s \to 0:

κ=limΔs0ΔαΔs=dαds\kappa = \lim_{\Delta s \to 0} \left| \frac{\Delta \alpha}{\Delta s} \right| = \left| \frac{d\alpha}{ds} \right|

Кривизна κ\kappa — это производная угла наклона касательной по длине дуги. В отличие от второй производной d2y/dx2d^2y/dx^2, которая показывает скорость изменения уклона относительно оси XX, истинная кривизна показывает скорость вращения касательной относительно реального пройденного пути по кривой. Это инвариантная величина: как бы мы ни вращали систему координат, длина дуги ss и взаимный угол поворота касательных Δα\Delta \alpha останутся неизменными, а значит, и кривизна κ\kappa в данной точке не изменится.

Идеальные эталоны: прямая и окружность

Чтобы лучше почувствовать работу формулы κ=dα/ds\kappa = |d\alpha/ds|, применим ее к двум фундаментальным геометрическим объектам, форма которых не меняется от точки к точке.

Прямая линия. У прямой линии касательная совпадает с самой линией во всех точках. При движении вдоль прямой угол наклона касательной остается константой: α=const\alpha = \text{const}. Следовательно, приращение угла на любом участке равно нулю (Δα=0\Delta \alpha = 0). Производная константы равна нулю, поэтому для прямой линии: κ=0\kappa = 0 Прямая — это линия нулевой кривизны.

Окружность. Рассмотрим окружность радиуса RR. Из элементарной геометрии известно, что длина дуги окружности ss прямо пропорциональна центральному углу, опирающемуся на эту дугу. Если центральный угол выражен в радианах, то s=Rθs = R \cdot \theta. При движении по окружности касательная всегда перпендикулярна радиусу. Это означает, что угол поворота касательной Δα\Delta \alpha в точности равен центральному углу Δθ\Delta \theta, на который повернулся радиус-вектор. Подставляя Δα\Delta \alpha вместо θ\theta в формулу длины дуги, получаем: s=RΔαs = R \cdot \Delta \alpha Отсюда выражаем отношение угла к дуге: Δαs=1R\frac{\Delta \alpha}{s} = \frac{1}{R} Поскольку это отношение постоянно на любом участке окружности, предел при Δs0\Delta s \to 0 равен самому отношению. Таким образом, кривизна окружности: κ=1R\kappa = \frac{1}{R}

Этот результат имеет колоссальное значение. Кривизна окружности постоянна во всех ее точках и обратно пропорциональна ее радиусу. Чем меньше радиус (тесный поворот), тем больше кривизна. Чем больше радиус, тем кривизна ближе к нулю. Земля имеет настолько большой радиус, что ее кривизна κ\kappa локально неотличима от нуля, из-за чего поверхность кажется нам плоской.

Радиус кривизны как мера локального сходства

Понимание того, что кривизна окружности равна 1/R1/R, позволяет ввести мощный геометрический инструмент для анализа любых, сколь угодно сложных кривых.

Если в некоторой точке MM произвольной кривой кривизна равна κ\kappaκ0\kappa \neq 0), мы можем поставить в соответствие этой точке величину RR, определяемую как:

R=1κR = \frac{1}{\kappa}

Эта величина называется радиусом кривизны кривой в точке MM.

Геометрический смысл радиуса кривизны заключается в следующем: это радиус такой идеальной окружности, которая «наилучшим образом» прилегает к нашей кривой в данной точке. Такая окружность называется соприкасающейся. Она не просто касается кривой (имеет общую касательную, как это делают бесконечно много окружностей). Она имеет с кривой общую касательную, общую нормаль и одинаковую кривизну в точке касания.

Представьте, что вы летите на самолете вдоль сложной извилистой траектории. В каждый момент времени вы можете зафиксировать штурвал. Самолет перестанет менять кривизну пути и полетит по идеальной окружности. Радиус этой окружности и есть радиус кривизны вашей изначальной траектории в той точке, где вы зафиксировали управление.

В точках, где кривая делает резкий изгиб (например, вершина параболы), кривизна κ\kappa велика, а радиус кривизны RR мал — соприкасающаяся окружность получается крошечной. На пологих участках кривизна κ\kappa стремится к нулю, а радиус кривизны RR уходит в бесконечность — соприкасающаяся окружность становится огромной, локально приближаясь к прямой линии. В точках перегиба, где выпуклость меняется на вогнутость, кривизна строго равна нулю, и радиус кривизны не существует (или формально равен бесконечности).

Физика изгиба: центростремительное ускорение

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

Из курса физики известна формула центростремительного ускорения для движения по окружности: an=v2Ra_n = \frac{v^2}{R} где vv — линейная скорость объекта, а RR — радиус окружности.

Используя соотношение R=1/κR = 1/\kappa, мы можем переписать эту формулу через кривизну траектории: an=v2κa_n = v^2 \cdot \kappa

Эта форма записи гораздо глубже отражает суть явления. Ускорение, которое испытывает тело (а значит, и сила, действующая на него согласно второму закону Ньютона F=manF = m a_n), прямо пропорционально квадрату скорости и локальной кривизне траектории.

Если пассажирский поезд движется с постоянной скоростью v=20v = 20 м/с, то боковая перегрузка, вдавливающая пассажиров в бок кресла, зависит исключительно от геометрии рельсов — от значения κ\kappa. Если путь прямой, κ=0\kappa = 0, боковое ускорение отсутствует. Если путь представляет собой дугу окружности радиусом 500500 метров, кривизна κ=1/500=0.002\kappa = 1/500 = 0.002 м1^{-1}. Боковое ускорение составит an=2020.002=0.8a_n = 20^2 \cdot 0.002 = 0.8 м/с2^2. Это постоянная и вполне комфортная величина.

Но что произойдет в точке стыковки прямого участка и круговой дуги? На прямой κ=0\kappa = 0. В следующую миллисекунду, когда поезд въезжает на дугу, кривизна мгновенно становится равной κ=0.002\kappa = 0.002. Происходит разрыв функции кривизны. Согласно формуле an=v2κa_n = v^2 \cdot \kappa, центростремительное ускорение также мгновенно прыгнет от 00 до 0.80.8 м/с2^2. В физике мгновенное изменение ускорения называется рывком (jerk). Пассажиры ощутят резкий удар вбок, а гребни колес поезда нанесут разрушительный удар по рельсам.

Инженерное решение: клотоида как кривая идеального перехода

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

Поскольку скорость поезда vv постоянна, для линейного роста ускорения ana_n необходимо, чтобы кривизна пути κ\kappa росла линейно по мере продвижения поезда. Математически это условие записывается так: κ=cs\kappa = c \cdot s где ss — пройденный путь (длина дуги от начала поворота), а cc — некий коэффициент пропорциональности.

Линия, кривизна которой изменяется строго пропорционально пройденному пути, называется клотоидой (или спиралью Корню).

В начале клотоиды (при s=0s = 0) ее кривизна κ=0\kappa = 0, что обеспечивает идеальную, гладкую стыковку с прямым участком пути. По мере движения поезда ss увеличивается, кривизна плавно растет, а радиус кривизны RR плавно уменьшается. Когда кривизна достигает значения 1/Rдуги1/R_{дуги}, клотоида заканчивается и идеально стыкуется с круговой дугой.

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

Фундаментальное определение кривизны через предел отношения угла смежности к длине дуги κ=dα/ds\kappa = |d\alpha/ds| дает ясное геометрическое и физическое понимание того, как изгибается линия в пространстве. Оно не зависит от того, как мы нарисуем оси координат. Однако на практике функции редко задаются через длину дуги и угол касательной. Чаще всего мы имеем дело с аналитическими уравнениями вида y=f(x)y = f(x) или параметрическими системами x(t),y(t)x(t), y(t). Для работы с такими функциями требуется перевести чистое геометрическое определение κ\kappa на язык производных декартовых координат, что позволит вычислять радиусы кривизны для любых графиков прямым дифференцированием.

Аналитический вывод формул кривизны для функций, заданных в явном виде

Аналитический вывод формул кривизны для функций, заданных в явном виде

Представьте себе инженера, проектирующего американские горки. Если он попытается рассчитать перегрузки, которые испытают пассажиры в мертвой петле, опираясь исключительно на вторую производную траектории yy'', составленный им проект окажется смертельно опасным. На крутых виражах, где полотно уходит вверх под углом в 4545^\circ или 6060^\circ, вторая производная может принимать огромные значения, в то время как реальная физическая кривизна трассы — а значит, и центростремительная сила — останется в пределах нормы. И наоборот: на пологих участках даже небольшое значение yy'' способно вызвать жесткий удар, если скорость движения высока.

Вторая производная не равна геометрической кривизне. Она является лишь «плоской» проекцией изгиба на вертикальную ось, привязанной к выбранной системе координат. Чтобы получить истинную, инвариантную меру изгиба кривой, нам необходимо преодолеть эту координатную зависимость. В этой статье мы шаг за шагом выведем аналитическую формулу кривизны для функций, заданных в явном виде y=f(x)y = f(x), и разберемся в глубоком геометрическом смысле каждого ее элемента.

Геометрический конфликт: почему вторая производная искажает изгиб

Напомним, что геометрическая кривизна κ\kappa в точке определяет скорость вращения касательной по мере движения вдоль самой кривой. Ее фундаментальное определение опирается на естественный параметр — длину дуги ss:

κ=dαds\kappa = \left| \frac{d\alpha}{ds} \right|

где:

  • κ\kappa — кривизна кривой в данной точке;
  • dαd\alpha — бесконечно малый (дифференциальный) угол поворота касательной;
  • dsds — дифференциал длины дуги кривой.

Когда функция задана в явном виде y=f(x)y = f(x), мы не знаем длину дуги ss напрямую. Нам доступна лишь независимая переменная xx, которая движется по горизонтали равномерным шагом dxdx.

Вторая производная yy'' по своему определению измеряет скорость изменения углового коэффициента касательной (то есть tanα\tan \alpha) по отношению к горизонтальному шагу dxdx:

y=d(tanα)dxy'' = \frac{d(\tan \alpha)}{dx}

В этом и кроется геометрический конфликт. Во-первых, yy'' измеряет изменение тангенса угла, а не самого угла α\alpha. Во-вторых, это изменение соотносится с шагом по горизонтали dxdx, а не с реальным путем dsds, пройденным по дуге кривой. Если кривая идет круто вверх, то даже при крошечном шаге dxdx расстояние dsds, преодолеваемое точкой вдоль кривой, оказывается значительно больше dxdx.

Чтобы получить формулу истинной кривизны, нам нужно связать дифференциалы dαd\alpha и dsds с дифференциалом независимой переменной dxdx.

Пошаговый аналитический вывод формулы кривизны

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

tanα=y\tan \alpha = y'

где yy' — первая производная функции y=f(x)y = f(x) по переменной xx. Из этого соотношения мы можем выразить сам угол α\alpha через арктангенс первой производной:

α=arctan(y)\alpha = \arctan(y')

Теперь наша задача — найти дифференциал угла dαd\alpha, то есть определить, как меняется этот угол при изменении аргумента xx на величину dxdx.

Шаг 1: Дифференцирование угла наклона

Применим правило дифференцирования сложной функции. Продифференцируем обе части равенства α=arctan(y)\alpha = \arctan(y') по переменной xx:

dαdx=ddx(arctan(y))\frac{d\alpha}{dx} = \frac{d}{dx} \big( \arctan(y') \big)

Вспомним, что производная арктангенса (arctan(u))=11+u2u\big(\arctan(u)\big)' = \frac{1}{1 + u^2} \cdot u'. В нашем случае роль промежуточного аргумента uu играет первая производная yy'. Следовательно:

dαdx=11+(y)2d(y)dx\frac{d\alpha}{dx} = \frac{1}{1 + (y')^2} \cdot \frac{d(y')}{dx}

Поскольку производная от первой производной по xx — это вторая производная yy'', получаем:

dαdx=y1+(y)2\frac{d\alpha}{dx} = \frac{y''}{1 + (y')^2}

Умножив обе части на дифференциал dxdx, мы находим точное аналитическое выражение для дифференциала угла наклона касательной:

dα=y1+(y)2dxd\alpha = \frac{y''}{1 + (y')^2} dx

где:

  • dαd\alpha — дифференциал угла наклона касательной;
  • yy' — первая производная функции;
  • yy'' — вторая производная функции;
  • dxdx — дифференциал независимой переменной.

Этот промежуточный результат чрезвычайно красив. Он показывает, что изменение угла наклона dαd\alpha прямо пропорционально второй производной yy'', но «гасится» знаменателем 1+(y)21 + (y')^2. Чем круче наклон кривой (чем больше величина yy'), тем медленнее меняется сам угол α\alpha при изменении тангенса.

Шаг 2: Дифференциал длины дуги

Теперь обратимся к знаменателю фундаментальной формулы кривизны — дифференциалу длины дуги dsds. Рассмотрим бесконечно малый прямоугольный треугольник, гипотенузой которого является элемент дуги dsds, а катетами — приращения по осям dxdx и dydy. По теореме Пифагора:

(ds)2=(dx)2+(dy)2(ds)^2 = (dx)^2 + (dy)^2

Вынесем (dx)2(dx)^2 за скобки в правой части равенства:

(ds)2=(1+(dydx)2)(dx)2(ds)^2 = \left( 1 + \left(\frac{dy}{dx}\right)^2 \right) (dx)^2

Вспомним, что отношение дифференциалов dydx\frac{dy}{dx} — это и есть первая производная yy'. Тогда выражение принимает вид:

(ds)2=(1+(y)2)(dx)2(ds)^2 = \big(1 + (y')^2\big) (dx)^2

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

ds=1+(y)2dxds = \sqrt{1 + (y')^2} dx

где:

  • dsds — дифференциал длины дуги;
  • yy' — первая производная функции;
  • dxdx — дифференциал независимой переменной.

Шаг 3: Объединение результатов

Теперь у нас есть все элементы, чтобы собрать воедино формулу кривизны. Подставим полученные выражения для dαd\alpha и dsds в исходное определение κ=dαds\kappa = \left| \frac{d\alpha}{ds} \right|:

κ=y1+(y)2dx1+(y)2dx\kappa = \left| \frac{\frac{y''}{1 + (y')^2} dx}{\sqrt{1 + (y')^2} dx} \right|

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

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

(1+(y)2)1+(y)2=(1+(y)2)1(1+(y)2)1/2=(1+(y)2)3/2\big(1 + (y')^2\big) \cdot \sqrt{1 + (y')^2} = \big(1 + (y')^2\big)^1 \cdot \big(1 + (y')^2\big)^{1/2} = \big(1 + (y')^2\big)^{3/2}

Таким образом, мы приходим к классической аналитической формуле кривизны плоской кривой, заданной явным уравнением y=f(x)y = f(x):

κ=y(1+(y)2)3/2\kappa = \frac{|y''|}{\big(1 + (y')^2\big)^{3/2}}

где:

  • κ\kappa — кривизна в данной точке;
  • yy' — первая производная функции в этой точке;
  • yy'' — вторая производная функции в этой точке.

Соответственно, радиус кривизны RR, будучи величиной, обратной кривизне, вычисляется по формуле:

R=(1+(y)2)3/2yR = \frac{\big(1 + (y')^2\big)^{3/2}}{|y''|}

где RR — радиус соприкасающейся окружности в исследуемой точке кривой.

Физический и геометрический смысл знаменателя

Почему в знаменателе формулы кривизны возникло именно выражение (1+(y)2)3/2\big(1 + (y')^2\big)^{3/2}? Какова его роль?

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

Случай 1: Экстремум функции (горизонтальная касательная)

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

y=0y' = 0

Подставим это значение в формулу кривизны:

κ=y(1+02)3/2=y\kappa = \frac{|y''|}{\big(1 + 0^2\big)^{3/2}} = |y''|

В точках экстремума, и только в них, вторая производная в точности равна геометрической кривизне кривой. Здесь нет координатного искажения, так как в окрестности этой точки направление движения вдоль кривой практически совпадает с направлением оси OxOx (то есть dsdxds \approx dx).

Случай 2: Крутой склон

Теперь представим точку, где кривая уходит круто вверх или вниз, например, под углом 4545^\circ. В этом случае первая производная y=1y' = 1 (или y=1y' = -1). Подставим это значение:

κ=y(1+12)3/2=y23/2=y8y2.83\kappa = \frac{|y''|}{\big(1 + 1^2\big)^{3/2}} = \frac{|y''|}{2^{3/2}} = \frac{|y''|}{\sqrt{8}} \approx \frac{|y''|}{2.83}

Посмотрите, что произошло: реальная геометрическая кривизна оказалась почти в три раза меньше, чем значение второй производной!

Это происходит потому, что при крутом наклоне один и тот же шаг по горизонтали dx=1dx = 1 заставляет нас пройти по дуге кривой гораздо большее расстояние ds=21.41ds = \sqrt{2} \approx 1.41. Изменение угла касательной «растягивается» по более длинному пути, из-за чего кривая изгибается более плавно, чем нам пытается показать вторая производная. Знаменатель (1+(y)2)3/2\big(1 + (y')^2\big)^{3/2} компенсирует это растяжение, возвращая нас к объективной геометрической реальности.

Подробный разбор практических примеров

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

Пример 1: Верификация формулы на окружности

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

y=R2x2y = \sqrt{R^2 - x^2}

где:

  • RR — радиус окружности;
  • xx — абсцисса, лежащая в интервале (R,R)(-R, R).

Найдем первую производную этой функции, используя правило дифференцирования сложной функции:

y=12R2x2(R2x2)=2x2R2x2=xR2x2y' = \frac{1}{2\sqrt{R^2 - x^2}} \cdot (R^2 - x^2)' = \frac{-2x}{2\sqrt{R^2 - x^2}} = -\frac{x}{\sqrt{R^2 - x^2}}

Теперь найдем вторую производную, применив правило дифференцирования частного (uv)=uvuvv2\left(\frac{u}{v}\right)' = \frac{u'v - uv'}{v^2}:

y=(x)R2x2x(R2x2)(R2x2)2y'' = -\frac{(x)' \cdot \sqrt{R^2 - x^2} - x \cdot \left(\sqrt{R^2 - x^2}\right)'}{(\sqrt{R^2 - x^2})^2}

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

y=1R2x2x(xR2x2)R2x2y'' = -\frac{1 \cdot \sqrt{R^2 - x^2} - x \cdot \left(-\frac{x}{\sqrt{R^2 - x^2}}\right)}{R^2 - x^2}

Приведем числитель к общему знаменателю R2x2\sqrt{R^2 - x^2}:

y=R2x2+x2R2x2R2x2=R2(R2x2)3/2y'' = -\frac{\frac{R^2 - x^2 + x^2}{\sqrt{R^2 - x^2}}}{R^2 - x^2} = -\frac{R^2}{(R^2 - x^2)^{3/2}}

Теперь подготовим знаменатель формулы кривизны. Сначала найдем выражение 1+(y)21 + (y')^2:

1+(y)2=1+(xR2x2)2=1+x2R2x2=R2x2+x2R2x2=R2R2x21 + (y')^2 = 1 + \left(-\frac{x}{\sqrt{R^2 - x^2}}\right)^2 = 1 + \frac{x^2}{R^2 - x^2} = \frac{R^2 - x^2 + x^2}{R^2 - x^2} = \frac{R^2}{R^2 - x^2}

Возведем полученное выражение в степень 3/23/2:

(1+(y)2)3/2=(R2R2x2)3/2=(R2)3/2(R2x2)3/2=R3(R2x2)3/2\big(1 + (y')^2\big)^{3/2} = \left(\frac{R^2}{R^2 - x^2}\right)^{3/2} = \frac{(R^2)^{3/2}}{(R^2 - x^2)^{3/2}} = \frac{R^3}{(R^2 - x^2)^{3/2}}

Подставим абсолютное значение второй производной y|y''| и знаменатель в общую формулу кривизны:

κ=y(1+(y)2)3/2=R2(R2x2)3/2R3(R2x2)3/2\kappa = \frac{|y''|}{\big(1 + (y')^2\big)^{3/2}} = \frac{\frac{R^2}{(R^2 - x^2)^{3/2}}}{\frac{R^3}{(R^2 - x^2)^{3/2}}}

Дроби в числителе и знаменателе имеют одинаковый знаменатель (R2x2)3/2(R^2 - x^2)^{3/2}, который благополучно сокращается. Остается:

κ=R2R3=1R\kappa = \frac{R^2}{R^3} = \frac{1}{R}

Мы получили потрясающий результат! Несмотря на то, что и первая производная yy', и вторая производная yy'' ведут себя крайне сложно и стремятся к бесконечности при приближении к краям окружности (x±Rx \to \pm R), их аналитическая комбинация в формуле кривизны дает абсолютно постоянную величину κ=1R\kappa = \frac{1}{R}. Это триумф математического анализа, подтверждающий строгую геометрическую инсистентность формулы.

Пример 2: Исследование кривизны квадратичной параболы

Рассмотрим стандартную параболу:

y=x2y = x^2

Вычислим ее производные:

y=2x,y=2y' = 2x, \quad y'' = 2

Подставим эти выражения в формулу кривизны:

κ(x)=2(1+4x2)3/2\kappa(x) = \frac{2}{(1 + 4x^2)^{3/2}}

Проанализируем, как меняется кривизна этой линии при движении вдоль оси OxOx:

  1. В вершине параболы (x=0x = 0) кривизна достигает своего максимального значения:

    κ(0)=2(1+0)3/2=2\kappa(0) = \frac{2}{(1 + 0)^{3/2}} = 2

    Соответственно, радиус соприкасающейся окружности в вершине равен R(0)=12=0.5R(0) = \frac{1}{2} = 0.5.

  2. При удалении от вершины (когда x±x \to \pm\infty) знаменатель формулы начинает стремительно расти. Например, в точке x=5x = 5:

    κ(5)=2(1+425)3/2=21013/2210150.00197\kappa(5) = \frac{2}{(1 + 4 \cdot 25)^{3/2}} = \frac{2}{101^{3/2}} \approx \frac{2}{1015} \approx 0.00197

    Кривизна упала практически до нуля (уменьшилась более чем в 1000 раз!).

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

Инженерные приложения и малые деформации

В строительной механике и теории упругости существует знаменитое уравнение Эйлера — Бернулли, описывающее изгиб упругой балки под действием нагрузки:

M(x)=EIκ(x)M(x) = E \cdot I \cdot \kappa(x)

где:

  • M(x)M(x) — изгибающий момент в сечении балки;
  • EE — модуль упругости материала (модуль Юнга);
  • II — момент инерции сечения балки;
  • κ(x)\kappa(x) — кривизна оси балки.

Если мы подставим в это уравнение точную формулу кривизны, то получим сложное нелинейное дифференциальное уравнение:

M(x)=EIy(1+(y)2)3/2M(x) = E \cdot I \cdot \frac{y''}{\big(1 + (y')^2\big)^{3/2}}

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

В реальных строительных конструкциях (мостах, перекрытиях зданий, несущих балках) прогибы под нагрузкой ничтожно малы по сравнению с длиной самой балки. Это значит, что угол наклона касательной α\alpha в любой точке практически равен нулю. Следовательно, первая производная yy' чрезвычайно мала:

y1    (y)20y' \ll 1 \implies (y')^2 \approx 0

В этом случае знаменатель формулы кривизны превращается в единицу:

(1+(y)2)3/2(1+0)3/2=1\big(1 + (y')^2\big)^{3/2} \approx (1 + 0)^{3/2} = 1

Благодаря этому допущению нелинейная кривизна заменяется простой второй производной:

κy\kappa \approx y''

И уравнение изгиба балки принимает классический линейный вид, знакомый любому студенту инженерного вуза:

y=M(x)EIy'' = \frac{M(x)}{E \cdot I}

Однако важно понимать границы применимости этого упрощения. Если мы проектируем гибкое крыло современного самолета, выполненное из композитных материалов, или длинную удочку, изгиб которых может быть огромным, пренебрегать квадратом первой производной (y)2(y')^2 категорически запрещено. В таких задачах инженерам приходится использовать полную, нелинейную формулу кривизны, которую мы сегодня вывели.

Сравнение кривизны различных функций

Для глубокого понимания геометрии изгиба полезно сравнить поведение кривизны различных математических семейств. Давайте сопоставим две кривые, которые часто путают из-за схожего внешнего вида: квадратичную параболу y=x2+1y = x^2 + 1 и цепную линию (график гиперболического косинуса) y=cosh(x)y = \cosh(x). Обе эти кривые проходят через точку (0,1)(0, 1), симметричны относительно оси ординат и направлены ветвями вверх.

Сравним их кривизну в нижней точке (при x=0x = 0):

Параметр Парабола y=x2+1y = x^2 + 1 Цепная линия y=cosh(x)y = \cosh(x)
Первая производная yy' 2x2x sinh(x)\sinh(x)
Вторая производная yy'' 22 cosh(x)\cosh(x)
Значение yy' при x=0x=0 00 00
Значение yy'' при x=0x=0 22 11
Кривизна κ\kappa при x=0x=0 22 11
Радиус кривизны RR при x=0x=0 0.50.5 1.01.0

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

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

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

Дифференциальная геометрия кривых, заданных параметрическими уравнениями

Дифференциальная геометрия кривых, заданных параметрическими уравнениями

Представьте, что вы проектируете «мертвую петлю» для американских горок. Пытаясь описать эту траекторию классическим уравнением вида y=f(x)y = f(x), вы мгновенно столкнетесь с непреодолимым математическим препятствием. В самой верхней точке петли, где вагонетка движется вертикально, касательная становится вертикальной, а производная yy' устремляется в бесконечность. Более того, в петле одной координате xx соответствуют сразу два значения yy — на подъеме и на спуске, что полностью нарушает само определение функции.

Чтобы описывать реальные траектории — от виражей гоночных трасс до движения планет — математики используют параметрическое представление. Вместо того чтобы связывать координаты xx и yy напрямую, мы выражаем каждую из них через независимый параметр tt (который в физике часто отождествляют со временем):

x=x(t),y=y(t)x = x(t), \quad y = y(t)

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

Дифференцирование параметрических функций: преодоление координатных сингулярностей

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

Пусть функции x(t)x(t) и y(t)y(t) дифференцируемы, причем x(t)0x'(t) \neq 0. Согласно правилу дифференцирования сложной функции, если yy является функцией от xx, а xx в свою очередь зависит от tt, то производная yy по tt выражается как:

dydt=dydxdxdt\frac{dy}{dt} = \frac{dy}{dx} \cdot \frac{dx}{dt}

Здесь dydt\frac{dy}{dt} — это производная y(t)y'(t), а dxdt\frac{dx}{dt} — это x(t)x'(t). Отсюда мы можем выразить первую производную yy по xx, которую обозначим как yxy'_x:

yx=dydx=y(t)x(t)y'_x = \frac{dy}{dx} = \frac{y'(t)}{x'(t)}

Эта формула показывает, что тангенс угла наклона касательной к параметрической кривой равен отношению скоростей изменения координат yy и xx.

Теперь сделаем следующий шаг и найдем вторую производную yxxy''_{xx}. По определению, вторая производная — это производная от первой производной по переменной xx:

yxx=ddx(yx)y''_{xx} = \frac{d}{dx} (y'_x)

Поскольку величина yxy'_x выражена через параметр tt, мы снова применяем правило дифференцирования сложной функции:

yxx=ddt(yx)dtdxy''_{xx} = \frac{d}{dt} (y'_x) \cdot \frac{dt}{dx}

Учитывая, что dtdx=1x(t)\frac{dt}{dx} = \frac{1}{x'(t)}, перепишем выражение:

yxx=ddt(y(t)x(t))1x(t)y''_{xx} = \frac{d}{dt} \left( \frac{y'(t)}{x'(t)} \right) \cdot \frac{1}{x'(t)}

Применим правило дифференцирования частного для выражения в скобках:

ddt(y(t)x(t))=y(t)x(t)y(t)x(t)(x(t))2\frac{d}{dt} \left( \frac{y'(t)}{x'(t)} \right) = \frac{y''(t)x'(t) - y'(t)x''(t)}{(x'(t))^2}

Здесь:

  • x(t)x'(t) и y(t)y'(t) — первые производные координат по параметру tt;
  • x(t)x''(t) и y(t)y''(t) — вторые производные координат по параметру tt.

Умножая полученную дробь на 1x(t)\frac{1}{x'(t)}, мы находим окончательную формулу для второй производной параметрической функции:

yxx=x(t)y(t)y(t)x(t)(x(t))3y''_{xx} = \frac{x'(t)y''(t) - y'(t)x''(t)}{(x'(t))^3}

Обратите внимание на знаменатель (x(t))3(x'(t))^3. Если точка движется так, что в некоторый момент времени ее горизонтальная скорость x(t)x'(t) становится равной нулю (что соответствует вертикальной касательной), то вторая производная yxxy''_{xx} обращается в бесконечность. Однако сама кривая в этой точке может быть абсолютно гладкой и иметь конечный, стабильный изгиб. Это доказывает, что вторая производная yxxy''_{xx} — это лишь проекционная характеристика, сильно зависящая от выбора осей координат, а не истинная геометрическая мера изгиба.

Аналитический вывод формулы параметрической кривизны

Чтобы получить истинную геометрическую кривизну κ\kappa, независимую от выбора системы координат, мы должны обратиться к формуле кривизны для явной функции y=f(x)y = f(x), которую мы вывели в предыдущих главах:

κ=yxx(1+(yx)2)3/2\kappa = \frac{|y''_{xx}|}{(1 + (y'_x)^2)^{3/2}}

Здесь:

  • yxy'_x — первая производная функции по xx;
  • yxxy''_{xx} — вторая производная функции по xx;
  • κ\kappa — геометрическая кривизна кривой.

Подставим в эту формулу наши параметрические выражения для yxy'_x и yxxy''_{xx}. Начнем со знаменателя:

1+(yx)2=1+(yx)2=(x)2+(y)2(x)21 + (y'_x)^2 = 1 + \left( \frac{y'}{x'} \right)^2 = \frac{(x')^2 + (y')^2}{(x')^2}

Возведем это выражение в степень 3/23/2:

(1+(yx)2)3/2=((x)2+(y)2(x)2)3/2=((x)2+(y)2)3/2x3(1 + (y'_x)^2)^{3/2} = \left( \frac{(x')^2 + (y')^2}{(x')^2} \right)^{3/2} = \frac{((x')^2 + (y')^2)^{3/2}}{|x'|^3}

Теперь подставим числитель yxx|y''_{xx}| и преобразованный знаменатель в общую формулу кривизны:

κ=xyyx(x)3((x)2+(y)2)3/2x3\kappa = \frac{\left| \frac{x'y'' - y'x''}{(x')^3} \right|}{\frac{((x')^2 + (y')^2)^{3/2}}{|x'|^3}}

Поскольку (x)3=x3|(x')^3| = |x'|^3, эти члены в числителе и знаменателе взаимно уничтожаются. Мы получаем элегантную, симметричную и фундаментальную формулу кривизны плоской параметрической кривой:

κ=x(t)y(t)y(t)x(t)((x(t))2+(y(t))2)3/2\kappa = \frac{|x'(t)y''(t) - y'(t)x''(t)|}{((x'(t))^2 + (y'(t))^2)^{3/2}}

Эта формула — шедевр дифференциальной геометрии. В ней больше нет сингулярностей, связанных с вертикальными касательными. Пока скорость движения точки отлична от нуля (то есть (x)2+(y)20(x')^2 + (y')^2 \neq 0), кривизна определена корректно в любой точке плоскости.

Кинематический мост: скорость, ускорение и векторное произведение

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

Пусть движение точки по плоскости задано радиус-вектором r(t)=(x(t),y(t))\vec{r}(t) = (x(t), y(t)). Тогда вектор скорости точки v(t)\vec{v}(t) — это первая производная радиус-вектора по времени:

v(t)=r(t)=(x(t),y(t))\vec{v}(t) = \vec{r}'(t) = (x'(t), y'(t))

Модуль вектора скорости (мгновенная скорость движения точки вдоль траектории) равен:

v=v=(x)2+(y)2v = |\vec{v}| = \sqrt{(x')^2 + (y')^2}

Вектор ускорения a(t)\vec{a}(t) — это вторая производная радиус-вектора:

a(t)=r(t)=(x(t),y(t))\vec{a}(t) = \vec{r}''(t) = (x''(t), y''(t))

Рассмотрим теперь числитель формулы кривизны: xyyx|x'y'' - y'x''|. В линейной алгебре эта величина представляет собой модуль псевдоскалярного (или косого) произведения двух векторов на плоскости — v\vec{v} и a\vec{a}. Если мы мысленно перенесем наши векторы в трехмерное пространство, добавив нулевую координату zz, то векторное произведение скорости и ускорения будет равно:

v×a=det(ijkxy0xy0)=(xyyx)k\vec{v} \times \vec{a} = \det \begin{pmatrix} \vec{i} & \vec{j} & \vec{k} \\ x' & y' & 0 \\ x'' & y'' & 0 \end{pmatrix} = (x'y'' - y'x'')\vec{k}

Модуль этого векторного произведения равен:

v×a=xyyx|\vec{v} \times \vec{a}| = |x'y'' - y'x''|

Подставив эти физические величины в формулу кривизны, мы получаем исключительно простую векторную запись:

κ=v×av3\kappa = \frac{|\vec{v} \times \vec{a}|}{v^3}

Здесь:

  • v\vec{v} — вектор скорости;
  • a\vec{a} — вектор ускорения;
  • vv — скалярная скорость (модуль вектора скорости);
  • κ\kappa — кривизна траектории.

Эта формула раскрывает глубокую связь между геометрией и физикой. Полное ускорение движущейся точки a\vec{a} всегда можно разложить на две составляющие: тангенциальное ускорение aτa_\tau (отвечающее за изменение величины скорости) и нормальное (центростремительное) ускорение ana_n (отвечающее за изменение направления движения):

a=aττ+ann\vec{a} = a_\tau \vec{\tau} + a_n \vec{n}

Где τ\vec{\tau} — единичный касательный вектор, а n\vec{n} — единичный вектор нормали. Из механики известно, что нормальное ускорение напрямую связано с кривизной траектории:

an=κv2    κ=anv2a_n = \kappa v^2 \implies \kappa = \frac{a_n}{v^2}

Поскольку тангенциальная составляющая ускорения коллинеарна скорости, векторное произведение v×a\vec{v} \times \vec{a} «замечает» только нормальную составляющую:

v×a=van|\vec{v} \times \vec{a}| = v \cdot a_n

Отсюда мгновенно следует, что an=v×ava_n = \frac{|\vec{v} \times \vec{a}|}{v}. Подставляя это в выражение для кривизны, мы вновь возвращаемся к формуле κ=v×av3\kappa = \frac{|\vec{v} \times \vec{a}|}{v^3}.

Таким образом, геометрическая кривизна — это мера «неэффективности» ускорения с точки зрения изменения скорости: чем больше ускорение отклоняется от направления движения, тем сильнее изгибается траектория.

Практикум: расчет кривизны для фундаментальных параметрических кривых

Чтобы закрепить понимание формулы, детально разберем два классических примера: окружность и циклоиду.

Пример 1. Параметрическая окружность

Проверим работоспособность формулы на объекте с известной геометрией. Зададим окружность радиуса RR с центром в начале координат параметрическими уравнениями:

x(t)=Rcost,y(t)=Rsintx(t) = R\cos t, \quad y(t) = R\sin t

Найдем первые и вторые производные по параметру tt:

x(t)=Rsint,y(t)=Rcostx'(t) = -R\sin t, \quad y'(t) = R\cos t

x(t)=Rcost,y(t)=Rsintx''(t) = -R\cos t, \quad y''(t) = -R\sin t

Вычислим числитель формулы кривизны (модуль косого произведения):

xyyx=(Rsint)(Rsint)(Rcost)(Rcost)=R2sin2t+R2cos2t=R2(sin2t+cos2t)=R2|x'y'' - y'x''| = |(-R\sin t)(-R\sin t) - (R\cos t)(-R\cos t)| = |R^2\sin^2 t + R^2\cos^2 t| = R^2(\sin^2 t + \cos^2 t) = R^2

Теперь вычислим знаменатель (куб скорости):

((x)2+(y)2)3/2=((Rsint)2+(Rcost)2)3/2=(R2sin2t+R2cos2t)3/2=(R2)3/2=R3((x')^2 + (y')^2)^{3/2} = ((-R\sin t)^2 + (R\cos t)^2)^{3/2} = (R^2\sin^2 t + R^2\cos^2 t)^{3/2} = (R^2)^{3/2} = R^3

Находим кривизну κ\kappa:

κ=R2R3=1R\kappa = \frac{R^2}{R^3} = \frac{1}{R}

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

Пример 2. Кривизна циклоиды

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

x(t)=a(tsint),y(t)=a(1cost)x(t) = a(t - \sin t), \quad y(t) = a(1 - \cos t)

Где aa — радиус катящегося колеса, а параметр tt — угол поворота колеса. Найдем производные:

x(t)=a(1cost),y(t)=asintx'(t) = a(1 - \cos t), \quad y'(t) = a\sin t

x(t)=asint,y(t)=acostx''(t) = a\sin t, \quad y''(t) = a\cos t

Вычислим числитель формулы кривизны:

xyyx=a(1cost)acostasintasint=a2(costcos2tsin2t)=a2(cost1)x'y'' - y'x'' = a(1 - \cos t) \cdot a\cos t - a\sin t \cdot a\sin t = a^2(\cos t - \cos^2 t - \sin^2 t) = a^2(\cos t - 1)

Так как нам нужен модуль этой величины, а cost1\cos t \leq 1, запишем:

xyyx=a2(1cost)|x'y'' - y'x''| = a^2(1 - \cos t)

Теперь найдем квадрат скорости (выражение под знаменателем):

(x)2+(y)2=a2(1cost)2+a2sin2t=a2(12cost+cos2t+sin2t)=a2(22cost)=2a2(1cost)(x')^2 + (y')^2 = a^2(1 - \cos t)^2 + a^2\sin^2 t = a^2(1 - 2\cos t + \cos^2 t + \sin^2 t) = a^2(2 - 2\cos t) = 2a^2(1 - \cos t)

Возведем полученное выражение в степень 3/23/2:

((x)2+(y)2)3/2=(2a2(1cost))3/2=22a3(1cost)3/2((x')^2 + (y')^2)^{3/2} = (2a^2(1 - \cos t))^{3/2} = 2\sqrt{2}a^3(1 - \cos t)^{3/2}

Теперь соберем формулу кривизны циклоиды:

κ=a2(1cost)22a3(1cost)3/2=122a1cost\kappa = \frac{a^2(1 - \cos t)}{2\sqrt{2}a^3(1 - \cos t)^{3/2}} = \frac{1}{2\sqrt{2}a\sqrt{1 - \cos t}}

Используем тригонометрическое тождество половинного угла: 1cost=2sin2(t/2)1 - \cos t = 2\sin^2(t/2). Для интервала t(0,2π)t \in (0, 2\pi) значение sin(t/2)>0\sin(t/2) > 0, поэтому 1cost=2sin(t/2)\sqrt{1 - \cos t} = \sqrt{2}\sin(t/2). Подставим это в формулу:

κ=122a2sin(t/2)=14asin(t/2)\kappa = \frac{1}{2\sqrt{2}a \cdot \sqrt{2}\sin(t/2)} = \frac{1}{4a\sin(t/2)}

Соответственно, радиус кривизны R=1/κR = 1/\kappa равен:

R=4asin(t/2)R = 4a\sin(t/2)

Этот результат поражает своей простотой и изяществом. Давайте проанализируем его физически:

  • В вершине циклоиды (t=πt = \pi), где колесо находится в наивысшей точке своего полета, sin(π/2)=1\sin(\pi/2) = 1. Кривизна минимальна и равна κmin=14a\kappa_{min} = \frac{1}{4a}, а радиус кривизны достигает максимума: Rmax=4aR_{max} = 4a. Это значит, что в верхней точке траектория циклоиды изгибается очень плавно — в четыре раза плавнее, чем само катящееся колесо!
  • В точках соприкосновения колеса с дорогой (t0t \to 0 или t2πt \to 2\pi), sin(t/2)0\sin(t/2) \to 0. Кривизна устремляется в бесконечность (κ\kappa \to \infty), а радиус кривизны падает до нуля (R0R \to 0). В этих точках циклоида образует острые «клювы» (точки возврата). В этот момент физическая скорость точки колеса падает до нуля, она мгновенно останавливается и резко меняет направление движения.

Сравнение кривых и инженерный смысл

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

Если состыковать прямой участок дороги (кривизна κ=0\kappa = 0) напрямую с дугой окружности (кривизна κ=1/R\kappa = 1/R), то в точке стыка кривизна претерпит мгновенный разрыв. Для пассажиров поезда это будет означать мгновенный боковой удар — бесконечный «рывок» (производная ускорения по времени).

Чтобы избежать этого, инженеры проектируют переходные кривые, задаваемые параметрически. Самым известным примером является клотоида (спираль Корню), в которой кривизна растет строго линейно с ростом пройденного пути ss:

κ(s)=Cs\kappa(s) = C \cdot s

Где CC — постоянный коэффициент масштаба. В параметрическом виде клотоида описывается через интегралы Френеля:

x(t)=0tcos(u2)du,y(t)=0tsin(u2)dux(t) = \int_0^t \cos(u^2) du, \quad y(t) = \int_0^t \sin(u^2) du

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

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

Соприкасающаяся окружность и определение координат центра кривизны

Соприкасающаяся окружность и определение координат центра кривизны

Представьте, что вы проектируете трассу для скоростного поезда. На прямом участке состав движется стабильно, но впереди его ждет поворот. Если соединить прямую линию с дугой окружности напрямую, в точке стыка пассажиры почувствуют мгновенный и сокрушительный боковой удар. Физически это объясняется тем, что центростремительное ускорение мгновенно прыгает от нуля до конечного значения. Чтобы избежать катастрофы, инженерам необходимо точно знать, как именно искривляется траектория в каждой микроскопической точке пути. Но как измерить это искривление геометрически? Касательная линия здесь не поможет — она лишь показывает направление движения, оставаясь абсолютно прямой. Нам нужен более совершенный инструмент аппроксимации — «идеальный круговой шаблон», который прижимается к кривой плотнее, чем любая прямая.

Геометрия наилучшего соприкосновения

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

Однако касательная совершенно «не чувствует» изгиб. Чтобы измерить кривизну, нам необходимо перейти к аппроксимации второго порядка. Простейшей нелинейной кривой с постоянным изгибом является окружность. Попробуем построить окружность, которая ближе всего подходит к нашей кривой в окрестности заданной точки.

Геометрически соприкасающуюся окружность можно определить через предельный переход. Возьмем на непрерывной кривой три близкие точки: M1M_1, MM и M2M_2. Через любые три точки, не лежащие на одной прямой, можно провести единственную окружность. Начнем стягивать точки M1M_1 и M2M_2 к точке MM вдоль кривой. В пределе, когда обе точки сольются с MM, секущая окружность займет строго определенное положение.

Соприкасающаяся окружность (или окружность кривизны) в точке MM — это предельное положение окружности, проходящей через три точки кривой, когда эти точки стремятся к точке MM.

Эта окружность обладает уникальными свойствами:

  1. В точке соприкосновения у нее и у исходной кривой совпадают координаты: yкривой=yокрy_{\text{кривой}} = y_{\text{окр}}.
  2. Совпадают их первые производные (общая касательная): yкривой=yокрy'_{\text{кривой}} = y'_{\text{окр}}.
  3. Совпадают их вторые производные (одинаковая скорость изменения наклона): yкривой=yокрy''_{\text{кривой}} = y''_{\text{окр}}.

Благодаря совпадению трех важнейших характеристик, соприкасающаяся окружность буквально «вжимается» в кривую, повторяя ее форму в бесконечно малой окрестности точки MM лучше, чем любая другая окружность. Радиус этой окружности RR называется радиусом кривизны, а ее центр CCцентром кривизны.

Вывод координат центра кривизны для явной функции

Наша задача — найти координаты центра кривизны C(ξ,η)C(\xi, \eta) для кривой, заданной уравнением y=f(x)y = f(x).

Пусть в точке M(x,y)M(x, y) кривая имеет радиус кривизны RR, а касательная к ней наклонена под углом α\alpha к оси абсцисс. Направление изгиба кривой указывает нам, в какую сторону направлена вогнутость. Центр соприкасающейся окружности всегда лежит на нормали к кривой (прямой, перпендикулярной касательной), проведенной в сторону вогнутости графика.

Рассмотрим геометрические связи. Вектор, направленный из точки M(x,y)M(x, y) в центр кривизны C(ξ,η)C(\xi, \eta), имеет длину RR и перпендикулярен касательной.

Угол наклона касательной α\alpha связан с первой производной фундаментальным соотношением:

tanα=y\tan \alpha = y'

Поскольку нормаль перпендикулярна касательной, угол ее наклона к оси OxOx равен α+π2\alpha + \frac{\pi}{2} (если кривая вогнута вверх, то есть y>0y'' > 0). Спроецируем отрезок MCMC, равный по длине радиусу кривизны RR, на координатные оси:

ξ=x+Rcos(α+π2)=xRsinα\xi = x + R \cos\left(\alpha + \frac{\pi}{2}\right) = x - R \sin \alpha

η=y+Rsin(α+π2)=y+Rcosα\eta = y + R \sin\left(\alpha + \frac{\pi}{2}\right) = y + R \cos \alpha

Теперь нам необходимо выразить тригонометрические функции sinα\sin \alpha и cosα\cos \alpha через производную yy'. Воспользуемся известными тригонометрическими тождествами:

cosα=11+tan2α=11+(y)2\cos \alpha = \frac{1}{\sqrt{1 + \tan^2 \alpha}} = \frac{1}{\sqrt{1 + (y')^2}}

sinα=tanα1+tan2α=y1+(y)2\sin \alpha = \frac{\tan \alpha}{\sqrt{1 + \tan^2 \alpha}} = \frac{y'}{\sqrt{1 + (y')^2}}

Вспомним формулу радиуса кривизны RR для функции, заданной в явном виде, которую мы получили ранее:

R=(1+(y)2)3/2yR = \frac{(1 + (y')^2)^{3/2}}{|y''|}

Подставим выражения для sinα\sin \alpha, cosα\cos \alpha и RR в уравнения для координат центра (ξ,η)(\xi, \eta). Рассмотрим случай, когда y>0y'' > 0 (график вогнут вверх, y=y|y''| = y''):

ξ=x(1+(y)2)3/2yy1+(y)2=xy(1+(y)2)y\xi = x - \frac{(1 + (y')^2)^{3/2}}{y''} \cdot \frac{y'}{\sqrt{1 + (y')^2}} = x - \frac{y' (1 + (y')^2)}{y''}

η=y+(1+(y)2)3/2y11+(y)2=y+1+(y)2y\eta = y + \frac{(1 + (y')^2)^{3/2}}{y''} \cdot \frac{1}{\sqrt{1 + (y')^2}} = y + \frac{1 + (y')^2}{y''}

Если график вогнут вниз (y<0y'' < 0), то центр кривизны лежит с противоположной стороны нормали, а y=y|y''| = -y''. Подстановка знаков в тригонометрические проекции и раскрытие модуля приводят к абсолютно тем же самым итоговым формулам! Знаменатель yy'' автоматически учитывает направление изгиба кривой благодаря своему знаку.

Таким образом, мы получили классические аналитические формулы для координат центра кривизны плоской кривой:

ξ=xy(1+(y)2)y\xi = x - \frac{y' (1 + (y')^2)}{y''}

η=y+1+(y)2y\eta = y + \frac{1 + (y')^2}{y''}

где:

  • x,yx, y — координаты исходной точки на кривой;
  • yy' — первая производная функции в этой точке;
  • yy'' — вторая производная функции в этой точке;
  • ξ\xi (кси) — абсцисса центра кривизны;
  • η\eta (эта) — ордината центра кривизны.

Параметрический случай: переход к инвариантной форме

На практике траектории движения физических тел или сложные инженерные профили редко удается описать простой явной функцией y=f(x)y = f(x). Гораздо чаще используется параметрическое представление:

{x=x(t)y=y(t)\begin{cases} x = x(t) \\ y = y(t) \end{cases}

где tt — вспомогательный параметр (например, время).

Чтобы найти координаты центра кривизны для параметрической кривой, мы можем воспользоваться правилами дифференцирования сложных функций. Выразим производные по xx через производные по параметру tt (обозначим их штрихами xx' и yy'):

yx=yxy'_x = \frac{y'}{x'}

yxx=yxxy(x)3y''_{xx} = \frac{y'' x' - x'' y'}{(x')^3}

Подставим эти соотношения в ранее выведенные формулы для ξ\xi и η\eta. Начнем с абсциссы ξ\xi:

ξ=xyx(1+(yx)2)yxxy(x)3\xi = x - \frac{\frac{y'}{x'} \left(1 + \left(\frac{y'}{x'}\right)^2\right)}{\frac{y'' x' - x'' y'}{(x')^3}}

Приведем выражение в числителе к общему знаменателю:

1+(yx)2=(x)2+(y)2(x)21 + \left(\frac{y'}{x'}\right)^2 = \frac{(x')^2 + (y')^2}{(x')^2}

Теперь подставим это обратно и сократим степени xx':

ξ=xyx(x)2+(y)2(x)2yxxy(x)3=xy((x)2+(y)2)xyyx\xi = x - \frac{\frac{y'}{x'} \cdot \frac{(x')^2 + (y')^2}{(x')^2}}{\frac{y'' x' - x'' y'}{(x')^3}} = x - \frac{y' \left((x')^2 + (y')^2\right)}{x' y'' - y' x''}

Проделаем аналогичную операцию для ординаты η\eta:

η=y+1+(yx)2yxxy(x)3=y+(x)2+(y)2(x)2yxxy(x)3=y+x((x)2+(y)2)xyyx\eta = y + \frac{1 + \left(\frac{y'}{x'}\right)^2}{\frac{y'' x' - x'' y'}{(x')^3}} = y + \frac{\frac{(x')^2 + (y')^2}{(x')^2}}{\frac{y'' x' - x'' y'}{(x')^3}} = y + \frac{x' \left((x')^2 + (y')^2\right)}{x' y'' - y' x''}

Мы получили изящные симметричные формулы для вычисления координат центра кривизны параметрической линии:

ξ=xy((x)2+(y)2)xyyx\xi = x - \frac{y' \left((x')^2 + (y')^2\right)}{x' y'' - y' x''}

η=y+x((x)2+(y)2)xyyx\eta = y + \frac{x' \left((x')^2 + (y')^2\right)}{x' y'' - y' x''}

где все производные берутся строго по параметру tt. Знаменатель этих дробей xyyxx' y'' - y' x'' представляет собой не что иное, как псевдоскалярное (косое) произведение векторов скорости v(x,y)\vec{v}(x', y') и ускорения a(x,y)\vec{a}(x'', y''). Оно напрямую характеризует интенсивность и направление поворота вектора скорости в пространстве.

Сравнение кривых по их кривизне

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

Семейство кривых Уравнение Кривизна κ\kappa Поведение центра кривизны C(ξ,η)C(\xi, \eta)
Прямая линия y=kx+by = kx + b κ=0\kappa = 0 Центр кривизны не существует (уходит в бесконечность)
Окружность (xa)2+(yb)2=R2(x-a)^2 + (y-b)^2 = R^2 κ=1R=const\kappa = \frac{1}{R} = \text{const} Центр кривизны неподвижен и совпадает с геометрическим центром окружности (a,b)(a, b)
Парабола y=ax2y = ax^2 κ=2a(1+4a2x2)3/2\kappa = \frac{2a}{(1 + 4a^2x^2)^{3/2}} Максимальна в вершине (x=0x=0), стремительно затухает до 00 при x±x \to \pm\infty. Центр непрерывно движется по сложной траектории.
Эллипс {x=acosty=bsint\begin{cases} x = a\cos t \\ y = b\sin t \end{cases} κ=ab(a2sin2t+b2cos2t)3/2\kappa = \frac{ab}{(a^2\sin^2 t + b^2\cos^2 t)^{3/2}} Периодически меняется. Максимальна на вершинах большой оси, минимальна на вершинах малой оси.

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

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

Подробный разбор практических примеров

Закрепим полученный математический аппарат на двух детальных примерах: для явного и параметрического задания кривых.

Пример 1. Расчет для параболы в произвольной точке

Найдем координаты центра кривизны и составим уравнение соприкасающейся окружности для параболы y=12x2y = \frac{1}{2}x^2 в точке M(2,2)M(2, 2).

Шаг 1. Нахождение производных. Вычислим первую и вторую производные нашей функции:

y=(12x2)=xy' = \left(\frac{1}{2}x^2\right)' = x

y=(x)=1y'' = (x)' = 1

В заданной точке x=2x = 2:

y(2)=2y'(2) = 2

y(2)=1y''(2) = 1

Шаг 2. Расчет радиуса кривизны. Применим формулу радиуса RR:

R=(1+(y)2)3/2y=(1+22)3/21=53/21=5511.18R = \frac{(1 + (y')^2)^{3/2}}{|y''|} = \frac{(1 + 2^2)^{3/2}}{1} = \frac{5^{3/2}}{1} = 5\sqrt{5} \approx 11.18

Шаг 3. Вычисление координат центра кривизны. Подставим значения в формулы для ξ\xi и η\eta:

ξ=xy(1+(y)2)y=22(1+22)1=2251=210=8\xi = x - \frac{y' (1 + (y')^2)}{y''} = 2 - \frac{2 \cdot (1 + 2^2)}{1} = 2 - \frac{2 \cdot 5}{1} = 2 - 10 = -8

η=y+1+(y)2y=2+1+221=2+51=7\eta = y + \frac{1 + (y')^2}{y''} = 2 + \frac{1 + 2^2}{1} = 2 + \frac{5}{1} = 7

Таким образом, центр соприкасающейся окружности находится в точке C(8,7)C(-8, 7).

Шаг 4. Составление уравнения соприкасающейся окружности. Каноническое уравнение окружности с центром (ξ,η)(\xi, \eta) и радиусом RR имеет вид:

(xξ)2+(yη)2=R2(x - \xi)^2 + (y - \eta)^2 = R^2

Подставим наши данные:

(x(8))2+(y7)2=(55)2(x - (-8))^2 + (y - 7)^2 = (5\sqrt{5})^2

(x+8)2+(y7)2=125(x + 8)^2 + (y - 7)^2 = 125

Проверим, лежит ли исходная точка M(2,2)M(2, 2) на этой окружности:

(2+8)2+(27)2=102+(5)2=100+25=125(2 + 8)^2 + (2 - 7)^2 = 10^2 + (-5)^2 = 100 + 25 = 125

Равенство выполняется идеально. Мы нашли окружность, которая в окрестности точки (2,2)(2, 2) сливается с параболой.


Пример 2. Расчет для вершины эллипса (параметрический вид)

Найдем координаты центра кривизны эллипса, заданного параметрически:

{x=acosty=bsint\begin{cases} x = a\cos t \\ y = b\sin t \end{cases}

в его правой вершине, где параметр t=0t = 0 (координаты точки M(a,0)M(a, 0)).

Шаг 1. Вычисление производных по параметру. Найдем первые и вторые производные функций x(t)x(t) и y(t)y(t) по переменной tt:

x=asintx' = -a\sin t

y=bcosty' = b\cos t

x=acostx'' = -a\cos t

y=bsinty'' = -b\sin t

Подставим значение параметра t=0t = 0:

x(0)=asin(0)=0x'(0) = -a\sin(0) = 0

y(0)=bcos(0)=by'(0) = b\cos(0) = b

x(0)=acos(0)=ax''(0) = -a\cos(0) = -a

y(0)=bsin(0)=0y''(0) = -b\sin(0) = 0

Шаг 2. Расчет знаменателя и вспомогательных сумм. Вычислим величину, стоящую в знаменателе формул координат центра кривизны:

D=xyyx=00b(a)=abD = x' y'' - y' x'' = 0 \cdot 0 - b \cdot (-a) = ab

Также найдем сумму квадратов первых производных:

(x)2+(y)2=02+b2=b2(x')^2 + (y')^2 = 0^2 + b^2 = b^2

Шаг 3. Вычисление координат центра кривизны. Подставим полученные значения в параметрические формулы:

ξ=xy((x)2+(y)2)xyyx=abb2ab=ab2a=a2b2a\xi = x - \frac{y' \left((x')^2 + (y')^2\right)}{x' y'' - y' x''} = a - \frac{b \cdot b^2}{ab} = a - \frac{b^2}{a} = \frac{a^2 - b^2}{a}

η=y+x((x)2+(y)2)xyyx=0+0b2ab=0\eta = y + \frac{x' \left((x')^2 + (y')^2\right)}{x' y'' - y' x''} = 0 + \frac{0 \cdot b^2}{ab} = 0

Центр кривизны эллипса в его правой вершине имеет координаты C(a2b2a,0)C\left(\frac{a^2 - b^2}{a}, 0\right).

Шаг 4. Геометрическая интерпретация. Радиус кривизны RR в этой точке равен расстоянию от точки M(a,0)M(a, 0) до центра кривизны C(a2b2a,0)C\left(\frac{a^2 - b^2}{a}, 0\right):

R=aξ=a(ab2a)=b2aR = a - \xi = a - \left(a - \frac{b^2}{a}\right) = \frac{b^2}{a}

Этот результат чрезвычайно важен для понимания геометрии эллипса. Если эллипс сильно вытянут (aba \gg b), то радиус кривизны в его вершине b2a\frac{b^2}{a} становится очень маленьким — кривая в этой области совершает чрезвычайно крутой, резкий разворот. Напротив, у вершин малой оси (при t=π2t = \frac{\pi}{2}) аналогичный расчет даст радиус кривизны R=a2bR = \frac{a^2}{b}, что соответствует очень пологому изгибу.

Особые случаи: когда формулы теряют силу

Математические формулы всегда имеют свои границы применимости. Внимательно взглянем на знаменатели в формулах координат центра кривизны:

  • Для явной функции в знаменателе стоит yy''.
  • Для параметрической кривой в знаменателе стоит xyyxx' y'' - y' x''.

Что происходит, когда эти величины обращаются в ноль?

Если y=0y'' = 0 в некоторой точке x0x_0, это означает, что в данной точке кривизна линии равна нулю (κ=0\kappa = 0). Геометрически в окрестности этой точки кривая ведет себя практически как прямая линия. Радиус кривизны стремится к бесконечности (RR \to \infty), а координаты центра кривизны ξ\xi и η\eta устремляются в бесконечность. Соприкасающаяся окружность в такой точке вырождается в прямую линию — касательную.

Типичным примером являются точки перегиба функций (например, x=0x = 0 для кубической параболы y=x3y = x^3). При переходе через такую точку кривизна меняет свой знак, а центр соприкасающейся окружности скачком перемещается с одной стороны кривой на другую, проходя через бесконечность.

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

Геометрическая теория эволют и эвольвент плоских кривых

Геометрическая теория эволют и эвольвент плоских кривых

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

Двойственная природа эволюты: геометрическое место центров и огибающая нормалей

В предыдущих разделах курса мы научились находить координаты центра соприкасающейся окружности (ξ,η)(\xi, \eta) для каждой точки плоской кривой. Если мы начнем непрерывно двигаться вдоль исходной кривой, центр соприкасающейся окружности также начнет перемещаться, описывая свою собственную траекторию.

Эволюта — это геометрическое место центров кривизны данной плоской кривой.

Если исходная кривая задана параметрически как r(t)=(x(t),y(t))\vec{r}(t) = (x(t), y(t)), то ее эволюта представляет собой новую кривую rE(t)=(ξ(t),η(t))\vec{r}_E(t) = (\xi(t), \eta(t)), координаты которой вычисляются по уже известным нам формулам. Однако у эволюты есть второе, не менее фундаментальное геометрическое определение, связывающее ее с понятием огибающей.

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

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

Докажем теорему, которая связывает эти два взгляда на эволюту.

Теорема об огибающей семейства нормалей

Эволюта плоской кривой является огибающей семейства ее нормалей.

Пусть исходная кривая задана в явном виде y=f(x)y = f(x). Уравнение нормали к этой кривой в некоторой точке с абсциссой uu имеет вид:

yf(u)=1f(u)(xu)y - f(u) = -\frac{1}{f'(u)} (x - u)

Здесь uu — абсцисса точки на исходной кривой, f(u)f(u) — ордината этой точки, f(u)f'(u) — значение первой производной функции в точке uu, а xx и yy — текущие координаты точек нормали.

Для удобства перенесем все члены в одну сторону и запишем уравнение семейства нормалей как функцию трех переменных F(x,y,u)=0F(x, y, u) = 0, где uu выступает в роли параметра семейства:

F(x,y,u)=(xu)+f(u)(yf(u))=0F(x, y, u) = (x - u) + f'(u)(y - f(u)) = 0

В этом уравнении F(x,y,u)F(x, y, u) — функция, задающая семейство нормалей, uu — параметр семейства (абсцисса точки касания), f(u)f(u) — значение исходной функции, f(u)f'(u) — ее первая производная, а xx и yy — координаты точек на нормали.

Согласно классической теории дифференциальной геометрии, для нахождения огибающей семейства кривых F(x,y,u)=0F(x, y, u) = 0 необходимо составить систему уравнений, состоящую из самого уравнения семейства и его частной производной по параметру uu:

{F(x,y,u)=0Fu(x,y,u)=0\begin{cases} F(x, y, u) = 0 \\ \frac{\partial F}{\partial u}(x, y, u) = 0 \end{cases}

Здесь F(x,y,u)F(x, y, u) — уравнение семейства нормалей, Fu(x,y,u)\frac{\partial F}{\partial u}(x, y, u) — частная производная этой функции по параметру uu, а xx и yy — координаты точек огибающей.

Вычислим частную производную функции FF по переменной uu, помня, что xx и yy при этом дифференцировании считаются константами:

Fu=u[xu+f(u)yf(u)f(u)]=1+f(u)yf(u)f(u)(f(u))2\frac{\partial F}{\partial u} = \frac{\partial}{\partial u} \left[ x - u + f'(u)y - f'(u)f(u) \right] = -1 + f''(u)y - f''(u)f(u) - (f'(u))^2

В данном выражении Fu\frac{\partial F}{\partial u} обозначает частную производную функции семейства по параметру uu, f(u)f'(u) и f(u)f''(u) — первую и вторую производные исходной функции соответственно, а xx и yy — координаты, рассматриваемые при дифференцировании как константы.

Сгруппируем слагаемые:

Fu=1(f(u))2+f(u)(yf(u))=0\frac{\partial F}{\partial u} = -1 - (f'(u))^2 + f''(u)(y - f(u)) = 0

Здесь f(u)f'(u) — первая производная исходной функции, f(u)f''(u) — ее вторая производная, uu — параметр семейства, а yy — ордината искомой точки огибающей.

Из этого уравнения мы можем выразить разность yf(u)y - f(u), предполагая, что вторая производная f(u)0f''(u) \neq 0:

yf(u)=1+(f(u))2f(u)y - f(u) = \frac{1 + (f'(u))^2}{f''(u)}

В этой формуле yy — ордината точки огибающей, f(u)f(u) — значение исходной функции в точке uu, а f(u)f'(u) и f(u)f''(u) — первая и вторая производные в этой точке.

Отсюда находим ординату точки огибающей, которую обозначим как η\eta:

η=f(u)+1+(f(u))2f(u)\eta = f(u) + \frac{1 + (f'(u))^2}{f''(u)}

Здесь η\eta — ордината центра кривизны (точки огибающей), f(u)f(u) — значение функции в точке uu, а f(u)f'(u) и f(u)f''(u) — первая и вторая производные исходной функции соответственно.

Теперь подставим полученное выражение для yf(u)y - f(u) в исходное уравнение нормали F(x,y,u)=0F(x, y, u) = 0, чтобы найти абсциссу точки огибающей, которую обозначим как ξ\xi:

(ξu)+f(u)(1+(f(u))2f(u))=0(\xi - u) + f'(u) \left( \frac{1 + (f'(u))^2}{f''(u)} \right) = 0

В данном уравнении ξ\xi — абсцисса центра кривизны, uu — абсцисса исходной точки, f(u)f'(u) и f(u)f''(u) — первая и вторая производные исходной функции соответственно.

ξ=uf(u)(1+(f(u))2)f(u)\xi = u - \frac{f'(u)(1 + (f'(u))^2)}{f''(u)}

Здесь ξ\xi — абсцисса центра кривизны, uu — абсцисса исходной точки, а f(u)f'(u) и f(u)f''(u) — первая и вторая производные функции в этой точке.

Сопоставив полученные выражения для ξ\xi и η\eta с формулами координат центра кривизны, мы видим идеальное совпадение.

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

Эвольвента и механическая модель развертывания нити

Если кривая BB является эволютой для кривой AA, то кривая AA называется эвольвентой (или разверткой) для кривой BB.

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

Переведем эту физическую модель на строгий язык векторного анализа. Пусть эволюта задана своим естественным параметрическим уравнением rE(s)\vec{r}_E(s), где ss — длина дуги, отсчитываемая от некоторой начальной точки.

В каждой точке эволюты определен единичный вектор касательной:

τE(s)=drEds\vec{\tau}_E(s) = \frac{d\vec{r}_E}{ds}

Здесь τE(s)\vec{\tau}_E(s) — единичный вектор касательной к эволюте, rE\vec{r}_E — радиус-вектор точки эволюты, а ss — длина дуги эволюты.

Пусть общая длина нити равна LL. Когда нить частично намотана на эволюту на участке от 00 до ss, длина ее свободного натянутого конца равна LsL - s. Так как свободный конец нити направлен по касательной к эволюте в точке rE(s)\vec{r}_E(s), векторное уравнение эвольвенты rI(s)\vec{r}_I(s) запишется в виде:

rI(s)=rE(s)+(Ls)τE(s)\vec{r}_I(s) = \vec{r}_E(s) + (L - s)\vec{\tau}_E(s)

В этом уравнении rI(s)\vec{r}_I(s) — радиус-вектор точки эвольвенты, rE(s)\vec{r}_E(s) — радиус-вектор точки эволюты, LL — полная длина нити, ss — длина дуги эволюты, а τE(s)\vec{\tau}_E(s) — единичный вектор касательной к эволюте.

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

drIds=drEds+dds[(Ls)τE(s)]\frac{d\vec{r}_I}{ds} = \frac{d\vec{r}_E}{ds} + \frac{d}{ds} \left[ (L - s)\vec{\tau}_E(s) \right]

Здесь drIds\frac{d\vec{r}_I}{ds} — производная радиус-вектора эвольвенты по длине дуги ss, drEds\frac{d\vec{r}_E}{ds} — производная радиус-вектора эволюты по длине дуги ss, LL — длина нити, а τE(s)\vec{\tau}_E(s) — единичный вектор касательной к эволюте.

Применим правило дифференцирования произведения:

drIds=τE(s)τE(s)+(Ls)dτEds\frac{d\vec{r}_I}{ds} = \vec{\tau}_E(s) - \vec{\tau}_E(s) + (L - s)\frac{d\vec{\tau}_E}{ds}

В этой формуле drIds\frac{d\vec{r}_I}{ds} — производная радиус-вектора эвольвенты по параметру ss, τE(s)\vec{\tau}_E(s) — единичный вектор касательной к эволюте, LL — длина нити, а dτEds\frac{d\vec{\tau}_E}{ds} — производная вектора касательной.

Первые два слагаемых взаимно уничтожаются. Согласно первой формуле Френе-Серре для плоских кривых, производная единичного вектора касательной по длине дуги пропорциональна единичному вектору нормали nE(s)\vec{n}_E(s) к эволюте:

drEds=τE(s)    dτEds=κE(s)nE(s)\frac{d\vec{r}_E}{ds} = \vec{\tau}_E(s) \implies \frac{d\vec{\tau}_E}{ds} = \kappa_E(s)\vec{n}_E(s)

Здесь rE\vec{r}_E — радиус-вектор эволюты, ss — длина ее дуги, τE(s)\vec{\tau}_E(s) — единичный вектор касательной, nE(s)\vec{n}_E(s) — единичный вектор нормали к эволюте, а κE(s)\kappa_E(s) — кривизна эволюты в данной точке.

Подставляя это выражение, получаем:

drIds=(Ls)κE(s)nE(s)\frac{d\vec{r}_I}{ds} = (L - s)\kappa_E(s)\vec{n}_E(s)

В данном уравнении drIds\frac{d\vec{r}_I}{ds} — производная радиус-вектора эвольвенты по длине дуги ss, LL — полная длина нити, κE(s)\kappa_E(s) — кривизна эволюты, а nE(s)\vec{n}_E(s) — единичный вектор нормали к ней.

Полученный результат чрезвычайно важен. Вектор drIds\frac{d\vec{r}_I}{ds} направлен по касательной к эвольвенте. Из формулы видно, что он коллинеарен вектору nE(s)\vec{n}_E(s) — нормали к эволюте.

Поскольку вектор нормали к эволюте nE\vec{n}_E перпендикулярен касательной к эволюте τE\vec{\tau}_E (которая совпадает с направлением натянутой нити), то и касательная к эвольвенте в любой точке перпендикулярна натянутой нити. Нить в каждый момент времени служит нормалью для строящейся эвольвенты.

Фундаментальные свойства пары «эволюта — эвольвента»

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

Свойство 1: Семейство эвольвент

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

Свойство 2: Теорема Гюйгенса о длине дуги эволюты

Пусть на некотором участке исходной кривой (эвольвенты) ее радиус кривизны R(s)R(s) изменяется монотонно. Тогда длина дуги эволюты между двумя точками C1C_1 и C2C_2 (которые являются центрами кривизны для точек M1M_1 и M2M_2 эвольвенты) равна абсолютному изменению радиуса кривизны эвольвенты на этом участке:

ΔsE=R(M2)R(M1)\Delta s_E = |R(M_2) - R(M_1)|

Здесь ΔsE\Delta s_E — длина дуги эволюты между двумя точками, а R(M1)R(M_1) и R(M2)R(M_2) — радиусы кривизны эвольвенты в соответствующих точках M1M_1 и M2M_2.

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

Свойство 3: Особые точки эволюты (точки возврата)

Если на исходной кривой есть точка, в которой кривизна достигает локального экстремума (максимума или минимума), то в соответствующей точке на эволюте возникает особенность — точка возврата (cusp).

Покажем это математически. Ранее мы выяснили, что вектор скорости движения вдоль эволюты равен:

drEds=dRdsτE(s)\frac{d\vec{r}_E}{ds} = \frac{dR}{ds} \vec{\tau}_E(s)

В этой формуле rE\vec{r}_E — радиус-вектор эволюты, ss — длина дуги исходной кривой, RR — радиус кривизны исходной кривой, а τE(s)\vec{\tau}_E(s) — единичный вектор касательной к эволюте.

Если в некоторой точке кривизна κ\kappa экстремальна, то и радиус кривизны RR имеет локальный экстремум. Следовательно, его производная по длине дуги обращается в нуль:

dRds=0\frac{dR}{ds} = 0

Здесь RR — радиус кривизны исходной кривой, а ss — длина ее дуги.

В этой точке вектор скорости движения по эволюте становится нулевым: drEds=0\frac{d\vec{r}_E}{ds} = \vec{0}. Точка на эволюте становится сингулярной. Двигаясь вдоль исходной кривой мимо ее вершины (точки экстремума кривизны), центр кривизны мгновенно останавливается, а затем начинает движение в обратном направлении под острым углом, формируя «клюв» или точку возврата первого рода.

Классические примеры расчетов

Пример 1: Эволюта параболы

Найдем уравнение эволюты для классической квадратичной параболы, заданной уравнением:

y=x2y = x^2

В этом уравнении yy — ордината квадратичной параболы, а xx — ее абсцисса.

Для этого нам понадобятся первая и вторая производные:

y=2x,y=2y' = 2x, \quad y'' = 2

Здесь yy' и yy'' — соответственно первая и вторая производные функции y=x2y = x^2 по переменной xx.

Воспользуемся формулами для координат центра кривизны (ξ,η)(\xi, \eta), которые мы подробно изучили ранее:

ξ=xy(1+(y)2)y\xi = x - \frac{y'(1 + (y')^2)}{y''}

В этой формуле ξ\xi — абсцисса центра кривизны, xx — абсцисса точки на кривой, а yy' и yy'' — первая и вторая производные функции в этой точке.

η=y+1+(y)2y\eta = y + \frac{1 + (y')^2}{y''}

Здесь η\eta — ордината центра кривизны, yy — ордината точки на кривой, а yy' и yy'' — первая и вторая производные функции.

Подставим наши производные в эти формулы:

ξ=x2x(1+4x2)2=xx(1+4x2)=xx4x3=4x3\xi = x - \frac{2x(1 + 4x^2)}{2} = x - x(1 + 4x^2) = x - x - 4x^3 = -4x^3

Здесь ξ\xi — абсцисса центра кривизны параболы, выраженная через ее абсциссу xx.

η=x2+1+4x22=x2+12+2x2=3x2+12\eta = x^2 + \frac{1 + 4x^2}{2} = x^2 + \frac{1}{2} + 2x^2 = 3x^2 + \frac{1}{2}

Здесь η\eta — ордината центра кривизны параболы, выраженная через ее абсциссу xx.

Мы получили параметрические уравнения эволюты, где роль параметра играет абсцисса xx исходной параболы:

{ξ=4x3η=3x2+12\begin{cases} \xi = -4x^3 \\ \eta = 3x^2 + \frac{1}{2} \end{cases}

В этой системе уравнений ξ\xi и η\eta — координаты центра кривизны параболы, а xx — ее абсцисса, выступающая в роли параметра.

Исключим параметр xx из этой системы, чтобы получить явную зависимость между координатами ξ\xi и η\eta. Из второго уравнения выразим x2x^2:

3x2=η12    x2=13(η12)3x^2 = \eta - \frac{1}{2} \implies x^2 = \frac{1}{3}\left(\eta - \frac{1}{2}\right)

Здесь xx — абсцисса точки исходной параболы, а η\eta — ордината центра кривизны.

Теперь возведем первое уравнение в квадрат, чтобы избавиться от знака минус и получить x6x^6:

ξ2=16x6=16(x2)3\xi^2 = 16x^6 = 16(x^2)^3

В этом уравнении ξ\xi — абсцисса центра кривизны параболы, а xx — абсцисса точки на исходной параболе.

Подставим выраженный ранее квадрат x2x^2 в это уравнение:

ξ2=16(13(η12))3=1627(η12)3\xi^2 = 16 \left( \frac{1}{3}\left(\eta - \frac{1}{2}\right) \right)^3 = \frac{16}{27} \left(\eta - \frac{1}{2}\right)^3

Здесь ξ\xi и η\eta — координаты центра кривизны параболы.

Умножим обе части на 2727:

27ξ2=16(η12)327\xi^2 = 16\left(\eta - \frac{1}{2}\right)^3

В этом уравнении ξ\xi и η\eta — координаты точек эволюты параболы в декартовой системе координат.

Полученное уравнение описывает полукубическую параболу (также известную как парабола Нейла).

Обратите внимание на точку (0,12)\left(0, \frac{1}{2}\right). При ξ=0\xi = 0 значение η=12\eta = \frac{1}{2}. В этой точке парабола Нейла имеет точку возврата («клюв»). Физически это объясняется тем, что в вершине исходной параболы (0,0)(0, 0) кривизна максимальна (κ=2\kappa = 2, где κ\kappa — кривизна; радиус кривизны минимален R=1/2R = 1/2, где RR — радиус кривизны). Именно здесь производная радиуса кривизны по длине дуги обращается в нуль, что и порождает сингулярность на эволюте.

Пример 2: Эвольвента окружности

Выведем параметрические уравнения эвольвенты окружности радиуса rr с центром в начале координат. Окружность задается параметрически:

{x(t)=rcosty(t)=rsint\begin{cases} x(t) = r\cos t \\ y(t) = r\sin t \end{cases}

В этой системе уравнений x(t)x(t) и y(t)y(t) — координаты точки на окружности, rr — радиус этой окружности, а tt — угловой параметр.

Представим нить, закрепленную в точке t=0t = 0 и намотанную на окружность. Начнем разматывать ее против часовой стрелки. При некотором значении угла tt точка касания нити с окружностью имеет координаты (rcost,rsint)(r\cos t, r\sin t).

Длина размотанного участка нити равна длине дуги окружности, соответствующей углу tt. Поскольку радиус равен rr, длина дуги составляет:

s=rts = r \cdot t

Здесь ss — длина дуги окружности, rr — ее радиус, а tt — угол в радианах.

Свободный конец нити направлен по касательной к окружности в точке касания в сторону размотки. Единичный вектор касательной к окружности в точке tt равен:

τ(t)=(sint,cost)\vec{\tau}(t) = (-\sin t, \cos t)

В этой формуле τ(t)\vec{\tau}(t) — единичный вектор касательной к окружности, а tt — угловой параметр.

Точка на конце нити (которая и чертит эвольвенту) смещена от точки касания вдоль вектора касательной на расстояние, равное длине размотанной нити rtr \cdot t. Запишем координаты (xI,yI)(x_I, y_I) этой точки:

xI(t)=x(t)+rtsint=rcost+rtsint=r(cost+tsint)x_I(t) = x(t) + rt\sin t = r\cos t + rt\sin t = r(\cos t + t\sin t)

Здесь xI(t)x_I(t) — абсцисса точки эвольвенты окружности, x(t)x(t) — абсцисса соответствующей точки окружности, rr — радиус окружности, а tt — угловой параметр.

yI(t)=y(t)rtcost=rsintrtcost=r(sinttcost)y_I(t) = y(t) - rt\cos t = r\sin t - rt\cos t = r(\sin t - t\cos t)

Здесь yI(t)y_I(t) — ордината точки эвольвенты окружности, y(t)y(t) — ордината соответствующей точки окружности, rr — радиус окружности, а tt — угловой параметр.

Мы получили классические уравнения эвольвенты окружности:

{xI(t)=r(cost+tsint)yI(t)=r(sinttcost)\begin{cases} x_I(t) = r(\cos t + t\sin t) \\ y_I(t) = r(\sin t - t\cos t) \end{cases}

В этой системе уравнений xI(t)x_I(t) и yI(t)y_I(t) — параметрические координаты точек эвольвенты окружности, rr — радиус базовой окружности, а tt — угловой параметр.

Инженерный триумф эвольвенты: теория зубчатых зацеплений

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

Общая нормаль к профилям зубьев в точке их соприкосновения должна всегда проходить через одну и ту же точку на линии центров валов — полюс зацепления.

Если это условие выполняется, отношение угловых скоростей шестерен (передаточное число) остается строго постоянным в любой микромомент вращения.

Эвольвента окружности идеально решает эту задачу. Профили зубьев обеих шестерен выполняются по эвольвенте. Базовыми окружностями (эволютами), с которых «сматывается» воображаемая нить, служат так называемые основные окружности шестерен.

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

Поскольку нормаль к эвольвенте всегда является касательной к ее базовой окружности, общая нормаль в точке контакта зубьев всегда совпадает с этой линией зацепления. Линия зацепления пересекает линию центров валов в фиксированной точке — полюсе зацепления.

Преимущества эвольвентного зацепления перед любыми другими (например, циклоидальными):

  1. Постоянство передаточного отношения: Полное отсутствие микроускорений и вибраций, что критично для высокоскоростных редукторов и точных приборов.
  2. Толерантность к монтажным погрешностям: Если в процессе сборки или из-за износа подшипников расстояние между центрами валов слегка изменится, линия зацепления просто немного повернется, но точка пересечения с линией центров останется неподвижной. Передаточное отношение не изменится, зацепление продолжит работать плавно.
  3. Технологичность производства: Инструмент для нарезания эвольвентных зубьев (червячная фреза) имеет простой прямолинейный профиль режущих кромок, что делает производство шестерен дешевым и высокоточным.

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

Прикладные аспекты кривизны в классической механике и инженерном проектировании

Прикладные аспекты кривизны в классической механике и инженерном проектировании

Представьте, что вы находитесь в поезде, который мчится по высокоскоростной магистрали со скоростью 300 км/ч300 \text{ км/ч}. Внезапно прямая траектория сменяется дугой окружности радиуса 800 м800 \text{ м}. Даже если поезд не сойдет с рельсов благодаря прочности колесных реборд, все пассажиры внутри вагона мгновенно ощутят сильнейший боковой удар: боковое ускорение скачкообразно вырастет от 00 до почти 8.7 м/с28.7 \text{ м/с}^2 (что близко к ускорению свободного падения). Этот сокрушительный толчок — не просто следствие высокой скорости, а результат математической сингулярности: разрыва функции геометрической кривизны в точке сопряжения прямой и окружности.

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


Динамика криволинейного движения и геометрия траекторий

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

Из кинематики известно, что вектор нормального ускорения всегда направлен по главной нормали к кривой (в сторону ее вогнутости), а его модуль выражается через скорость vv и геометрическую кривизу траектории κ\kappa:

an=v2κ=v2Ra_n = v^2 \kappa = \frac{v^2}{R}

где:

  • ana_n — нормальное ускорение движущейся точки,
  • vv — мгновенная скорость точки вдоль траектории,
  • κ\kappa — кривизна траектории в данной точке,
  • RR — радиус кривизны в данной точке (R=1/κR = 1/\kappa).

Согласно второму закону Ньютона, для удержания тела массой mm на траектории с кривизной κ\kappa внешние силы должны обеспечивать нормальную силу:

Fn=mv2κF_n = m v^2 \kappa

где FnF_n — центростремительная сила, удерживающая тело на криволинейной траектории. Если эта сила генерируется реакцией опоры (например, полотном дороги), то само тело по третьему закону Ньютона действует на опору с равной по модулю и противоположной по направлению центробежной силой.

Вертикальные кривые и эффект невесомости

Рассмотрим движение автомобиля по вертикальному профилю дороги, который на вершине холма аппроксимируется функцией y=f(x)y = f(x). В высшей точке холма касательная к профилю горизонтальна (y=0y' = 0), а значит, кривизна дороги равна абсолютному значению второй производной: κ=y\kappa = |y''|.

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

mgN=mv2κmg - N = m v^2 \kappa

где:

  • mm — масса автомобиля,
  • gg — ускорение свободного падения,
  • NN — сила нормальной реакции дороги,
  • vv — скорость автомобиля в вершине холма,
  • κ\kappa — кривизна дорожного профиля в вершине холма.

Выразим силу нормальной реакции, которая определяет кажущийся вес водителя и пассажиров:

N=m(gv2κ)N = m(g - v^2 \kappa)

Из этой формулы следует фундаментальный вывод: сила давления автомобиля на дорогу уменьшается с ростом скорости и с увеличением кривизны холма. Если скорость автомобиля достигнет критического значения vcritv_{\text{crit}}, при котором нормальная реакция обратится в нуль (N=0N = 0), автомобиль оторвется от дороги и перейдет в режим свободного полета (невесомости):

vcrit=gκ=gRv_{\text{crit}} = \sqrt{\frac{g}{\kappa}} = \sqrt{g R}

где RR — радиус кривизны дороги в вершине холма. Проектировщики автомобильных дорог обязаны рассчитывать вертикальные радиусы кривизны холмов так, чтобы даже при превышении разрешенной скорости движения сила NN оставалась существенно больше нуля, гарантируя сцепление шин с асфальтом.

Эволюция петель американских горок: от круга к клотоиде

Одним из самых ярких примеров торжества дифференциальной геометрии над слепой инженерной интуицией стала история создания вертикальных петель на американских горках. Первые петли, построенные в середине и конце XIX века (начиная с парижской Centrifugal Railway в 1846 году), были строго круговыми. Их кривизна была постоянной: κ=1/R\kappa = 1/R.

Чтобы вагонетка могла успешно преодолеть верхнюю точку круговой петли радиуса RR, сила тяжести в этой точке должна быть полностью компенсирована центростремительным ускорением. В верхней точке (h=2Rh = 2R) сила нормальной реакции NtopN_{\text{top}} и сила тяжести направлены вниз:

Ntop+mg=mvtop2R    Ntop=m(vtop2Rg)N_{\text{top}} + mg = \frac{m v_{\text{top}}^2}{R} \implies N_{\text{top}} = m\left(\frac{v_{\text{top}}^2}{R} - g\right)

Для предотвращения падения вагонетки необходимо, чтобы Ntop0N_{\text{top}} \geq 0, откуда минимальная скорость на вершине составляет vtop=gRv_{\text{top}} = \sqrt{gR}. Согласно закону сохранения механической энергии (пренебрегая трением), скорость вагонетки на входе в петлю v0v_0 (где высота h=0h = 0) связана со скоростью на вершине соотношением:

mv022=mvtop22+mg(2R)    v02=vtop2+4gR\frac{m v_0^2}{2} = \frac{m v_{\text{top}}^2}{2} + mg(2R) \implies v_0^2 = v_{\text{top}}^2 + 4gR

Подставляя минимальную скорость на вершине, получаем необходимую скорость на входе в круговую петлю:

v05gRv_0 \geq \sqrt{5gR}

Теперь рассчитаем нормальное ускорение an,0a_{n, 0}, которое испытывают пассажиры в самой нижней точке на входе в петлю при этой скорости:

an,0=v02R5ga_{n, 0} = \frac{v_0^2}{R} \geq 5g

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

N0=man,0+mg6mgN_0 = m a_{n, 0} + mg \geq 6mg

Пассажиры испытывали перегрузку не менее 6g6g (а на практике, из-за трения и необходимости иметь запас скорости, мгновенные перегрузки на круговых петлях вроде Flip Flap Railway достигали опасных для здоровья 12g12g14g14g), что приводило к травмам шейных позвонков и кратковременной потере зрения (блэкауту).

Решением проблемы стало использование петель переменной кривизны — так называемых клотоидных петель (в форме капли или перевернутой слезы). В такой петле кривизна κ(s)\kappa(s) непрерывно увеличивается по мере подъема вагонетки.

Параметр сравнения Круговая петля Клотоидная петля
Кривизна на входе Высокая (κ=1/R\kappa = 1/R) Близкая к нулю (κ0\kappa \to 0)
Кривизна на вершине Такая же (κ=1/R\kappa = 1/R) Максимальная (κ=κmax\kappa = \kappa_{\text{max}})
Входная скорость v0v_0 Очень высокая (требуется для преодоления вершины) Умеренная (за счет малого радиуса на вершине)
Максимальная перегрузка От 6g6g до 14g14g (опасна для здоровья) Комфортные 3g3g4g4g
Характер изменения сил Мгновенный удар на входе Плавное нарастание и спад

В клотоидной петле на входе, где скорость vv максимальна, кривизна κ\kappa минимальна, что удерживает перегрузку v2κv^2 \kappa в безопасных пределах. На вершине, где скорость падает, кривизна возрастает до максимума, обеспечивая необходимое центростремительное ускорение для удержания вагонетки на рельсах при гораздо меньшей кинетической энергии.


Математическое проектирование дорожных переходов и виражей

В дорожном и железнодорожном строительстве сопряжение прямых участков дорог с круговыми закруглениями — критически важная задача. Прямой участок имеет нулевую кривизну (κ=0\kappa = 0), а круговой — постоянную ненулевую кривизну (κ=1/R\kappa = 1/R).

Если соединить их напрямую, возникнет разрыв первого рода в функции кривизны. С точки зрения динамики автомобиля это означает мгновенный скачок боковой силы. Скорость изменения бокового ускорения называется рывком (jerk) JJ:

J=dandtJ = \frac{d a_n}{dt}

При постоянной скорости движения vv вдоль пути ss (так что ds=vdtds = v \, dt) выразим рывок через скорость изменения кривизны по длине дуги:

J=ddt(v2κ)=v3dκdsJ = \frac{d}{dt}(v^2 \kappa) = v^3 \frac{d\kappa}{ds}

где:

  • JJ — рывок (скорость изменения бокового ускорения),
  • vv — скорость движения транспортного средства,
  • κ\kappa — кривизна траектории как функция длины дуги ss.

При непосредственном сопряжении прямой и окружности величина dκ/dsd\kappa/ds стремится к бесконечности, порождая бесконечный теоретический рывок. На практике это приводит к резкому удару, быстрому износу колесных пар и рельсов, а для водителей автомобилей делает невозможным плавный поворот руля.

Клотоида как идеальная переходная кривая

Чтобы исключить скачок ускорения, между прямым и круговым участками укладывают переходную кривую, кривизна которой плавно и непрерывно нарастает от 00 до 1/R1/R. Идеальной переходной кривой является клотоида (в отечественной инженерной практике ее также называют спиралью Корню).

Закон изменения кривизны клотоиды по длине дуги ss строго линеен:

κ(s)=cs\kappa(s) = c \cdot s

где cc — постоянный коэффициент пропорциональности. Если общая длина переходной кривой равна LL, а на ее конце кривизна должна достигнуть значения кривизны кругового участка 1/R1/R, то коэффициент cc равен:

c=1RL    κ(s)=sRLc = \frac{1}{R L} \implies \kappa(s) = \frac{s}{R L}

где:

  • ss — расстояние, пройденное вдоль переходной кривой от ее начала (0sL0 \leq s \leq L),
  • RR — радиус сопрягаемой круговой дуги,
  • LL — полная длина переходной кривой.

При движении по клотоиде с постоянной скоростью vv боковое ускорение нарастает строго линейно во времени:

an(t)=v2κ(s(t))=v2(vt)RL=v3tRLa_n(t) = v^2 \kappa(s(t)) = \frac{v^2 (v t)}{R L} = \frac{v^3 t}{R L}

При этом рывок остается постоянным и конечным на всем протяжении переходного участка:

J=v3RLJ = \frac{v^3}{R L}

В инженерном проектировании длина переходного участка LL рассчитывается исходя из ограничений на максимально допустимый физиологический рывок (для железных дорог общего пользования обычно принимают Jmax0.30.5 м/с3J_{\text{max}} \approx 0.3 \text{–} 0.5 \text{ м/с}^3):

Lv3RJmaxL \geq \frac{v^3}{R J_{\text{max}}}

Геометрия виражей: компенсация центробежных сил

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

Чтобы минимизировать зависимость безопасности движения от трения, дорожное полотно на поворотах проектируют с наклоном внутрь — строят вираж (superelevation).

Рассмотрим поперечный профиль дороги, наклоненный к горизонту под углом θ\theta. На автомобиль действуют сила тяжести mgmg и сила нормальной реакции дороги NN. Наша цель — найти такой угол наклона θ\theta, при котором автомобиль может проходить поворот радиуса RR на расчетной скорости vv вообще без участия сил бокового трения (когда колеса не испытывают бокового усилия сдвига).

В отсутствие бокового трения запишем уравнения равновесия в проекциях на вертикальную и горизонтальную (радиальную) оси:

Ncosθ=mgN \cos\theta = mg

Nsinθ=mv2κN \sin\theta = m v^2 \kappa

где:

  • θ\theta — угол наклона виража к горизонту,
  • NN — сила нормальной реакции дороги,
  • mm — масса автомобиля,
  • gg — ускорение свободного падения,
  • vv — расчетная скорость движения,
  • κ\kappa — кривизна дороги в плане (κ=1/R\kappa = 1/R).

Разделив второе уравнение на первое, мы исключаем массу автомобиля и силу реакции опоры, получая классическую формулу виража:

tanθ=v2κg=v2gR\tan\theta = \frac{v^2 \kappa}{g} = \frac{v^2}{g R}

Так как вдоль переходной кривой (клотоиды) кривизна κ(s)\kappa(s) меняется линейно от нуля до максимального значения, угол наклона виража θ(s)\theta(s) также должен плавно увеличиваться вдоль пути. Для малых углов справедливо приближение tanθθ\tan\theta \approx \theta, что означает линейный закон изменения высоты внешнего рельса (или внешней кромки асфальта) относительно внутреннего:

θ(s)v2sgRL\theta(s) \approx \frac{v^2 s}{g R L}

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


Строительная механика: упругая линия балки при больших деформациях

В теории упругости изгиб стержней и балок описывается фундаментальным законом Эйлера — Бернулли: изгибающий момент MM в любом поперечном сечении балки прямо пропорционален изменению кривизны ее нейтральной оси под действием нагрузки:

M(x)=EIκ(x)M(x) = E I \kappa(x)

где:

  • M(x)M(x) — изгибающий момент в сечении с координатой xx,
  • EE — модуль упругости материала (модуль Юнга),
  • II — момент инерции поперечного сечения балки относительно нейтральной оси,
  • κ(x)\kappa(x) — кривизна изогнутой оси балки.

В подавляющем большинстве инженерных расчетов строительных конструкций (мостов, перекрытий зданий) деформации предполагаются крайне малыми. В этом случае угол наклона касательной к изогнутой оси балки близок к нулю (y1y' \ll 1), поэтому квадратом первой производной в формуле кривизны пренебрегают:

(y)20    κ=y(1+(y)2)3/2y(y')^2 \approx 0 \implies \kappa = \frac{y''}{(1 + (y')^2)^{3/2}} \approx y''

Это приводит к линейному дифференциальному уравнению изгиба:

EIy=M(x)E I y'' = M(x)

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

y(1+(y)2)3/2=M(x)EI\frac{y''}{(1 + (y')^2)^{3/2}} = \frac{M(x)}{E I}

Точный расчет изгиба консольной балки под действием силы

Проведем математический анализ изгиба тонкой консольной балки длиной LL, жестко защемленной на левом конце (x=0x=0) и нагруженной на свободном правом конце вертикальной силой PP.

При значительном изгибе правый конец балки смещается не только по вертикали на величину δy\delta_y, но и по горизонтали внутрь на величину δx\delta_x. Пусть фактическая координата нагруженного конца балки равна (xend,yend)(x_{\text{end}}, y_{\text{end}}). Изгибающий момент в произвольной точке балки с координатой xx равен:

M(x)=P(xendx)M(x) = P(x_{\text{end}} - x)

Тогда точное уравнение изгиба принимает вид:

EIy(1+(y)2)3/2=P(xendx)E I \frac{y''}{(1 + (y')^2)^{3/2}} = P(x_{\text{end}} - x)

Для решения этой нелинейной задачи перейдем к естественной параметризации кривой. Обозначим через ss длину дуги балки, отсчитываемую от точки защемления (0sL0 \leq s \leq L), а через θ(s)\theta(s) — угол наклона касательной к балке в точке ss относительно горизонтали. Тогда координаты точек балки выражаются через тригонометрические функции угла наклона:

dxds=cosθ(s),dyds=sinθ(s)\frac{dx}{ds} = \cos\theta(s), \quad \frac{dy}{ds} = \sin\theta(s)

Кривизна упругой линии в терминах угла наклона и длины дуги записывается как первая производная:

κ(s)=dθds\kappa(s) = \frac{d\theta}{ds}

Выразим изгибающий момент через интеграл по длине дуги. Координата x(s)x(s) произвольной точки балки равна:

x(s)=0scosθ(u)dux(s) = \int_{0}^{s} \cos\theta(u) \, du

Следовательно, координата нагруженного конца балки составляет xend=x(L)=0Lcosθ(u)dux_{\text{end}} = x(L) = \int_{0}^{L} \cos\theta(u) \, du. Изгибающий момент в точке ss равен:

M(s)=P(xendx(s))=PsLcosθ(u)duM(s) = P(x_{\text{end}} - x(s)) = P \int_{s}^{L} \cos\theta(u) \, du

Подставим это выражение в точный закон Эйлера — Бернулли:

EIdθds=PsLcosθ(u)duE I \frac{d\theta}{ds} = P \int_{s}^{L} \cos\theta(u) \, du

Чтобы избавиться от интеграла в правой части, продифференцируем обе части уравнения по переменной ss:

EId2θds2=Pcosθ(s)E I \frac{d^2\theta}{ds^2} = -P \cos\theta(s)

Разделим обе части на EIEI и перенесем все члены в левую часть:

d2θds2+PEIcosθ=0\frac{d^2\theta}{ds^2} + \frac{P}{E I} \cos\theta = 0

Мы получили нелинейное дифференциальное уравнение второго порядка, которое математически полностью эквивалентно уравнению колебаний математического маятника! Роль времени здесь играет длина дуги ss, а роль угла отклонения маятника — угол наклона касательной к балке θ\theta.

Для решения этого уравнения зададим граничные условия:

  1. На левом конце при s=0s = 0 балка жестко защемлена горизонтально: θ(0)=0\theta(0) = 0.
  2. На правом конце при s=Ls = L изгибающий момент равен нулю, следовательно, кривизна равна нулю: dθdss=L=0\left.\frac{d\theta}{ds}\right|_{s=L} = 0.

Интегрируя это уравнение один раз (умножив на dθ/dsd\theta/ds), получаем:

12(dθds)2+PEIsinθ=C\frac{1}{2}\left(\frac{d\theta}{ds}\right)^2 + \frac{P}{E I} \sin\theta = C

Используя граничное условие на свободном конце (θ(L)=θL\theta(L) = \theta_L, где θL\theta_L — пока неизвестный максимальный угол наклона на конце балки), найдем константу интегрирования CC:

C=PEIsinθL    dθds=2PEI(sinθLsinθ)C = \frac{P}{E I} \sin\theta_L \implies \frac{d\theta}{ds} = \sqrt{\frac{2P}{E I}(\sin\theta_L - \sin\theta)}

Разделяя переменные и интегрируя от 00 до LL, мы приходим к эллиптическому интегралу, связывающему нагрузку PP и максимальный угол наклона θL\theta_L:

0θLdθsinθLsinθ=2PEIL\int_{0}^{\theta_L} \frac{d\theta}{\sqrt{\sin\theta_L - \sin\theta}} = \sqrt{\frac{2P}{E I}} L

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

Сопоставление линейной и точной моделей

Сравним результаты точного расчета (Elastica) и классического линейного приближения для консольной балки. Введем безразмерный параметр нагрузки:

α=PL2EI\alpha = \frac{P L^2}{E I}

В рамках линейной теории максимальный прогиб на конце балки δlinear\delta_{\text{linear}} и максимальный угол наклона θlinear\theta_{\text{linear}} вычисляются по простым формулам:

δlinear=PL33EI=α3L\delta_{\text{linear}} = \frac{P L^3}{3 E I} = \frac{\alpha}{3} L

θlinear=PL22EI=α2\theta_{\text{linear}} = \frac{P L^2}{2 E I} = \frac{\alpha}{2}

Параметр нагрузки α\alpha Линейный прогиб δlinear/L\delta_{\text{linear}}/L Точный прогиб δexact/L\delta_{\text{exact}}/L Погрешность линейной модели
0.1 (малая нагрузка) 0.0330.033 0.0330.033 <0.1%< 0.1\%
1.0 (средняя нагрузка) 0.3330.333 0.3020.302 10.3%\approx 10.3\%
2.0 (высокая нагрузка) 0.6670.667 0.4930.493 35.3%\approx 35.3\%
5.0 (экстремальный изгиб) 1.6671.667 0.7140.714 133.5%\approx 133.5\% (физически бессмысленно)

При больших значениях параметра α\alpha линейная модель начинает стремительно терять адекватность. Например, при α=5\alpha = 5 линейная формула предсказывает вертикальный прогиб, превышающий длину самой балки в 1.671.67 раза, что физически невозможно для нерастяжимого стержня. Точное же решение показывает, что балка просто сильно сворачивается, ее вертикальный прогиб асимптотически стремится к длине балки LL, а горизонтальная проекция сокращается.

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


Замыкание мысли

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

Комплексный практикум по расчету и сравнению дифференциальных параметров кривых

Комплексный практикум по расчету и сравнению дифференциальных параметров кривых

Представьте себе, что вы проектируете трассу для высокоскоростного поезда, несущегося со скоростью 300 км/ч. Пассажиры внутри вагона не видят геометрию рельсов напрямую, но их тела мгновенно и чутко реагируют на малейшие изменения дифференциальных параметров пути. Почему визуально «гладкий» стык двух дуг может вызвать опасный боковой удар, и как математический анализ позволяет предотвратить аварии еще на этапе чертежей?

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


Методологический фундамент практикума

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

  1. Локальную кривизну κ\kappa, которая определяет мгновенную скорость изменения направления касательной относительно длины дуги.
  2. Радиус кривизны RR, задающий размер соприкасающейся окружности.
  3. Координаты центра кривизны (ξ,η)(\xi, \eta), описывающие траекторию мгновенного центра вращения (эволюту).

Вспомним базовые расчетные формулы, которые мы доказали в предыдущих главах. Для явной функции y=f(x)y = f(x) кривизна вычисляется как:

κ(x)=y(1+(y)2)3/2\kappa(x) = \frac{|y''|}{(1 + (y')^2)^{3/2}}

где:

  • yy' — первая производная функции по xx (тангенс угла наклона касательной);
  • yy'' — вторая производная функции по xx (скорость изменения наклона);
  • (1+(y)2)3/2(1 + (y')^2)^{3/2} — масштабный корректор, компенсирующий проекционные искажения при переходе от шага dxdx к шагу по дуге dsds.

Для параметрически заданной кривой x=x(t)x = x(t), y=y(t)y = y(t) формула принимает вид:

κ(t)=xyyx((x)2+(y)2)3/2\kappa(t) = \frac{|x' y'' - y' x''|}{((x')^2 + (y')^2)^{3/2}}

где штрихами обозначены производные по параметру tt.


Кейс 1. Битва в вершине: Парабола против цепной линии

Начнем с классического сравнения двух фундаментальных U-образных кривых: квадратичной параболы и цепной линии (которая физически описывает форму однородной нити, подвешенной за концы).

Пусть обе кривые имеют вершины в начале координат (0,0)(0,0) и одинаковую вторую производную в этой точке. Рассмотрим функции:

  • Парабола: f(x)=12x2f(x) = \frac{1}{2}x^2
  • Цепная линия: g(x)=cosh(x)1g(x) = \cosh(x) - 1

В точке x=0x = 0 обе функции дают y=0y = 0, их первые производные равны нулю (f(0)=g(0)=0f'(0) = g'(0) = 0), а вторые производные равны единице (f(0)=g(0)=1f''(0) = g''(0) = 1). Следовательно, в самой вершине их кривизна одинакова:

κf(0)=κg(0)=1\kappa_f(0) = \kappa_g(0) = 1

Соответственно, радиус соприкасающейся окружности для обеих кривых в вершине равен R=1R = 1.

Но что происходит, когда мы удаляемся от вершины? Как быстро эти кривые «выпрямляются», теряя свою кривизну? Проведем строгий асимптотический анализ при xx \to \infty.

Анализ параболы f(x)=12x2f(x) = \frac{1}{2}x^2

Найдем производные: f(x)=xf'(x) = x, f(x)=1f''(x) = 1. Подставим их в формулу кривизны:

κf(x)=1(1+x2)3/2\kappa_f(x) = \frac{1}{(1 + x^2)^{3/2}}

При больших значениях xx единица в знаменателе становится пренебрежимо малой по сравнению с x2x^2. Следовательно, мы можем записать асимптотическое поведение:

κf(x)1(x2)3/2=1x3\kappa_f(x) \approx \frac{1}{(x^2)^{3/2}} = \frac{1}{|x|^3}

Мы видим, что кривизна параболы убывает степенным (алгебраическим) образом — обратно пропорционально кубу расстояния от оси ординат.

Анализ цепной линии g(x)=cosh(x)1g(x) = \cosh(x) - 1

Напомним, что гиперболический косинус определяется через экспоненты:

cosh(x)=ex+ex2\cosh(x) = \frac{e^x + e^{-x}}{2}

Производные равны: g(x)=sinh(x)g'(x) = \sinh(x), g(x)=cosh(x)g''(x) = \cosh(x). Используя основное гиперболическое тождество cosh2(x)sinh2(x)=1\cosh^2(x) - \sinh^2(x) = 1, преобразуем знаменатель формулы кривизны:

κg(x)=cosh(x)(1+sinh2(x))3/2=cosh(x)(cosh2(x))3/2=1cosh2(x)\kappa_g(x) = \frac{\cosh(x)}{(1 + \sinh^2(x))^{3/2}} = \frac{\cosh(x)}{(\cosh^2(x))^{3/2}} = \frac{1}{\cosh^2(x)}

При xx \to \infty слагаемое exe^{-x} стремится к нулю, поэтому cosh(x)12ex\cosh(x) \approx \frac{1}{2}e^x. Подставим это в асимптотику:

κg(x)1(12ex)2=4e2x\kappa_g(x) \approx \frac{1}{(\frac{1}{2}e^x)^2} = 4e^{-2x}

Кривизна цепной линии убывает экспоненциально!

Сравнительный итог анализа

Несмотря на визуальное сходство графиков в окрестности вершины, их дифференциальное поведение на бесконечности принципиально различно.

Параметр Парабола y=12x2y = \frac{1}{2}x^2 Цепная линия y=cosh(x)1y = \cosh(x) - 1
Кривизна в вершине κ(0)\kappa(0) 11 11
Скорость затухания кривизны Степенная: $\sim \frac{1}{ x
Геометрический смысл Медленное выпрямление Стремительное, почти мгновенное выпрямление

Инсайт практикума: Вторая производная дает информацию лишь о локальном поведении кривой в бесконечно малой окрестности точки. Глобальное сравнение формы кривых требует обязательного учета первой производной, которая входит в знаменатель формулы кривизны в степени 3/23/2.


Кейс 2. Параметрический анализ эллипса и геометрия его эволюты

Перейдем к более сложному объекту — замкнутой параметрической кривой. Рассмотрим эллипс с полуосями a=4a = 4 (вдоль оси OxOx) и b=2b = 2 (вдоль оси OyOy):

{x(t)=4cos(t)y(t)=2sin(t)\begin{cases} x(t) = 4\cos(t) \\ y(t) = 2\sin(t) \end{cases}

где t[0,2π]t \in [0, 2\pi] — параметр.

Наша задача — рассчитать кривизну, радиус кривизны и координаты центра кривизны в двух характерных точках: в правой вершине (t=0t = 0, координаты (4,0)(4,0)) и в верхней вершине (t=π/2t = \pi/2, координаты (0,2)(0,2)), а также проанализировать форму эволюты этого эллипса.

Расчет производных

Вычислим первые две производные по параметру tt:

{x(t)=4sin(t)x(t)=4cos(t)и{y(t)=2cos(t)y(t)=2sin(t)\begin{cases} x'(t) = -4\sin(t) \\ x''(t) = -4\cos(t) \end{cases} \quad \text{и} \quad \begin{cases} y'(t) = 2\cos(t) \\ y''(t) = -2\sin(t) \end{cases}

Анализ в правой вершине (t=0t = 0)

Подставим t=0t = 0 в уравнения производных:

  • x(0)=0x'(0) = 0, x(0)=4x''(0) = -4
  • y(0)=2y'(0) = 2, y(0)=0y''(0) = 0

Вычислим кривизну κ(0)\kappa(0):

κ(0)=002(4)(02+22)3/2=88=1\kappa(0) = \frac{|0 \cdot 0 - 2 \cdot (-4)|}{(0^2 + 2^2)^{3/2}} = \frac{8}{8} = 1

Радиус кривизны в этой точке равен:

R(0)=1κ(0)=1R(0) = \frac{1}{\kappa(0)} = 1

Поскольку в точке (4,0)(4,0) касательная к эллипсу вертикальна, а сам эллипс вогнут влево (в сторону начала координат), центр соприкасающейся окружности должен лежать на оси OxOx левее самой точки на величину радиуса R=1R = 1.

Проверим это строгим расчетом координат центра кривизны (ξ,η)(\xi, \eta). Вспомним параметрические формулы для центра кривизны:

ξ=xy((x)2+(y)2)xyyx\xi = x - \frac{y' ((x')^2 + (y')^2)}{x' y'' - y' x''}

η=y+x((x)2+(y)2)xyyx\eta = y + \frac{x' ((x')^2 + (y')^2)}{x' y'' - y' x''}

Подставляя значения для t=0t = 0:

ξ(0)=42(0+4)8=41=3\xi(0) = 4 - \frac{2 \cdot (0 + 4)}{8} = 4 - 1 = 3

η(0)=0+0(0+4)8=0\eta(0) = 0 + \frac{0 \cdot (0 + 4)}{8} = 0

Центр кривизны находится в точке (3,0)(3,0).

Анализ в верхней вершине (t=π/2t = \pi/2)

Подставим t=π/2t = \pi/2 в уравнения производных:

  • x(π/2)=4x'(\pi/2) = -4, x(π/2)=0x''(\pi/2) = 0
  • y(π/2)=0y'(\pi/2) = 0, y(π/2)=2y''(\pi/2) = -2

Вычислим кривизну κ(π/2)\kappa(\pi/2):

κ(π/2)=4(2)00((4)2+02)3/2=864=18=0.125\kappa(\pi/2) = \frac{|-4 \cdot (-2) - 0 \cdot 0|}{((-4)^2 + 0^2)^{3/2}} = \frac{8}{64} = \frac{1}{8} = 0.125

Радиус кривизны в этой точке равен:

R(π/2)=1κ(π/2)=8R(\pi/2) = \frac{1}{\kappa(\pi/2)} = 8

В точке (0,2)(0,2) касательная горизонтальна, эллипс вогнут вниз. Центр соприкасающейся окружности должен лежать на оси OyOy ниже точки на величину радиуса R=8R = 8. Вычислим координаты центра кривизны:

ξ(π/2)=00(16+0)8=0\xi(\pi/2) = 0 - \frac{0 \cdot (16 + 0)}{8} = 0

η(π/2)=2+4(16+0)8=28=6\eta(\pi/2) = 2 + \frac{-4 \cdot (16 + 0)}{8} = 2 - 8 = -6

Центр кривизны находится в точке (0,6)(0,-6).

Геометрия эволюты эллипса

Если мы проведем аналогичные вычисления для всех значений параметра tt, то получим параметрические уравнения траектории центра кривизны (эволюты):

{ξ(t)=a2b2acos3(t)=3cos3(t)η(t)=a2b2bsin3(t)=6sin3(t)\begin{cases} \xi(t) = \frac{a^2 - b^2}{a} \cos^3(t) = 3\cos^3(t) \\ \eta(t) = -\frac{a^2 - b^2}{b} \sin^3(t) = -6\sin^3(t) \end{cases}

Это уравнение описывает растянутую вдоль оси ординат астроиду — звездообразную кривую с четырьмя острыми вершинами (точками возврата). Эти вершины соответствуют экстремумам кривизны исходного эллипса.


Кейс 3. Проектирование виражей: Клотоида против кубической параболы

Перейдем к важнейшей прикладной задаче дорожного и железнодорожного строительства — плавного сопряжения прямолинейного участка дороги с круговым закруглением.

Пусть нам нужно соединить прямую линию (где кривизна κ=0\kappa = 0) и дугу окружности радиуса R0R_0 (где кривизна постоянна и равна κ0=1/R0\kappa_0 = 1/R_0). Если соединить их напрямую, возникнет мгновенный скачок кривизны. Физически это означает, что боковое (центростремительное) ускорение an=v2κa_n = v^2 \kappa вырастет от нуля до конечного значения мгновенно. Возникнет бесконечный рывок (производная ускорения по времени), который приведет к сильному удару, разрушающему подвеску поезда и рельсы.

Для плавного изменения кривизны используют переходные кривые.

Идеальное решение: Клотоида (спираль Корню)

В клотоиде кривизна κ\kappa меняется строго линейно по мере движения вдоль кривой (пропорционально пройденной длине дуги ss):

κ(s)=Cs\kappa(s) = C \cdot s

где CC — постоянный коэффициент масштаба.

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

Инженерный компромисс: Кубическая парабола

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

y=kx3y = kx^3

где k>0k > 0 — малый коэффициент.

Давайте проведем строгий математический аудит этого компромисса и выясним, где именно и почему данная аппроксимация начинает давать опасные сбои.

Вычислим производные кубической параболы: y=3kx2y' = 3kx^2, y=6kxy'' = 6kx. Подставим их в явную формулу кривизны:

κ(x)=6kx(1+9k2x4)3/2\kappa(x) = \frac{6kx}{(1 + 9k^2 x^4)^{3/2}}

При малых значениях xx (в самом начале переходного участка), величина 9k2x49k^2 x^4 пренебрежимо мала по сравнению с единицей. Кроме того, длина дуги ss практически совпадает со своей проекцией на ось абсцисс (sxs \approx x). В этой области мы получаем:

κ(x)6kx6ks\kappa(x) \approx 6kx \approx 6ks

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

Точка отказа кубической параболы

Но что произойдет, если мы продолжим движение дальше по кубической параболе?

Знаменатель (1+9k2x4)3/2(1 + 9k^2 x^4)^{3/2} начнет стремительно расти и подавлять линейный рост числителя. В результате кривизна сначала замедлит свой рост, достигнет максимума, а затем начнет убывать, стремясь к нулю на бесконечности!

Найдем точку, в которой кривизна кубической параболы достигает своего теоретического максимума. Для этого исследуем функцию кривизны на экстремум. Найдем производную κ(x)\kappa(x) по xx (для простоты исследуем квадрат кривизны или воспользуемся правилом дифференцирования дроби):

dκdx=6k(1+9k2x4)3/26kx32(1+9k2x4)1/236k2x3(1+9k2x4)3\frac{d\kappa}{dx} = \frac{6k \cdot (1 + 9k^2 x^4)^{3/2} - 6kx \cdot \frac{3}{2}(1 + 9k^2 x^4)^{1/2} \cdot 36k^2 x^3}{(1 + 9k^2 x^4)^3}

Приравняем числитель к нулю и сократим на общие положительные множители 6k(1+9k2x4)1/26k(1 + 9k^2 x^4)^{1/2}:

(1+9k2x4)54k2x4=0(1 + 9k^2 x^4) - 54k^2 x^4 = 0

145k2x4=0    x4=145k21 - 45k^2 x^4 = 0 \implies x^4 = \frac{1}{45k^2}

Таким образом, точка максимума кривизны находится на расстоянии:

xmax=145k240.386kx_{\text{max}} = \frac{1}{\sqrt[4]{45k^2}} \approx \frac{0.386}{\sqrt{k}}

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


Сравнительный анализ переходных кривых

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

Свойство сравнения Идеальная клотоида Кубическая парабола y=kx3y = kx^3
Закон изменения кривизны Строго линейный по дуге: κ(s)=Cs\kappa(s) = C \cdot s Нелинейный: κ(x)=6kx(1+9k2x4)3/2\kappa(x) = \frac{6kx}{(1 + 9k^2 x^4)^{3/2}}
Поведение при малых xx Совпадает с параболой Практически линейно: κ6kx\kappa \approx 6kx
Глобальный предел кривизны Стремится к бесконечности Стремится к нулю при xx \to \infty
Физический эффект для водителя Плавный, монотонный поворот руля Необходимость контрруления на длинных дистанциях
Область применения Автомагистрали, скоростные ж/д пути Трамвайные пути, подъезды с малыми углами поворота

Заключение

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

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