Чому велосипед стійкий у русі? Модель Уіппла

1. Міф про гіроскоп

Поширене пояснення стверджує, що колеса, що обертаються, протидіють перекиданню так само, як гіроскоп протидіє зміні орієнтації. Експерименти спростовують це як основну причину: велосипеди з додатковими протиобертовими колесами (що скасовують сумарний кутовий момент) залишаються самостійкими, а велосипеди зі зворотним виносом переднього колеса стають нестійкими за будь-якої швидкості навіть зі звичайними гіроскопічними колесами. Гіроскопічні ефекти дають додатковий момент керма, але домінувальний механізм — геометричний зв'язок, описаний повною моделлю Уіппла.

2. Модель Уіппла — чотири тверді тіла

Еталонна модель велосипеда Уіппла розглядає машину як чотири з'єднані тверді тіла, ігноруючи деформацію шин, згин рами й керування вершником (некерована система «без рук»):

Два ступені свободи описують малі відхилення від рівноважного стану прямолінійного руху у вертикальному положенні: кут нахилу φ (крен задньої рами) і кут повороту керма δ (обертання вилки відносно рами).

3. Лінеаризовані рівняння руху

Лінеаризація навколо рівноважного стану (вертикальний стан, стала швидкість v уперед) дає компактне матричне рівняння другого порядку:

4. Власні моди: weave, capsize, wobble

Підставляючи q = q₀·e^(λt), рівняння перетворюються на задачу на власні значення, залежну від швидкості. З типових параметрів велосипеда виникають три характерні моди:

5. Винос, кут нахилу вилки та геометрія керма

Механічний винос (trail) — горизонтальна відстань між точкою контакту осі керма із землею й контактною плямою переднього колеса — це найвпливовіший геометричний параметр:

6. Діапазон самостійкості

Дослідження на основі моделі Уіппла (Meijaard, Papadopoulos, Ruina, Schwab, 2007) підтвердили за допомогою фізичного «некерованого» велосипеда, що самостійкість реальна й кількісно відповідає передбаченням лінеаризованих власних значень:

• Нижче v_weave (~4 м/с): мода weave нестійка — класичне похитування, яке відчуває початківець на швидкості пішохода, вимагаючи активної корекції керма.
• Між v_weave і v_capsize: і weave, і capsize загасають — велосипед самостійкий і відновлюється після малих збурень без рук на кермі.
• Вище v_capsize (якщо вона існує): мода capsize знову стає нестійкою, хоча й дуже повільно (велика стала часу), тож на звичайних швидкостях їзди це рідко помітно.

7. JavaScript-перевірка лінійної стійкості

// Обчислення власних значень моделі Уіппла на заданій швидкості (умовні матриці)
// M, C1, K0, K2 — матриці 2x2 з еталонних параметрів велосипеда
function stateMatrix(M, C1, K0, K2, v, g = 9.81) {
  // Побудова матриці стану 4x4 для x = [phi, delta, phiDot, deltaDot]
  const Minv = invert2x2(M);
  const K = addScaled(scale2x2(K0, g), K2, v * v); // g*K0 + v^2*K2
  const C = scale2x2(C1, v);
  const negMinvK = scale2x2(mul2x2(Minv, K), -1);
  const negMinvC = scale2x2(mul2x2(Minv, C), -1);
  // A = [[0, I], [-Minv*K, -Minv*C]]  (блочна форма 4x4)
  return assembleBlock(negMinvK, negMinvC);
}

function isStable(A) {
  const eigenvalues = eig4(A); // чисельний розв'язувач власних значень
  return eigenvalues.every(lambda => lambda.re < 0);
}

// Перебір швидкості для пошуку вікна самостійкості
function findStableRange(M, C1, K0, K2, vMin = 0, vMax = 10, step = 0.05) {
  const stableSpeeds = [];
  for (let v = vMin; v <= vMax; v += step) {
    const A = stateMatrix(M, C1, K0, K2, v);
    if (isStable(A)) stableSpeeds.push(v);
  }
  return { min: stableSpeeds[0], max: stableSpeeds[stableSpeeds.length-1] };
}

8. Інженерні застосування