Фізика · Динаміка рідин
·
⏱ ≈ 13 хв читання
·
Останнє оновлення: 23 червня 2026 р.
Балістика та аеродинамічний опір — від гарматних ядер до надзвукових снарядів
Ідеальна парабола Галілея — лише початок. Реальні снаряди зазнають нелінійного опору, зміни режиму поблизу числа Маха 1, сили Магнуса та відхилення Коріоліса. Точкова балістична модель із залежним від Маха коефіцієнтом опору пояснює, чому снайпери вносять поправки на вітер і температуру, а ракети йдуть пологими траєкторіями.
1. Сила опору та коефіцієнт опору
Сила опору (у напрямку, протилежному швидкості): F_D = ½ · ρ · v² · C_D ·
A ρ = густина повітря (кг/м³), v = швидкість (м/с), C_D = коефіцієнт опору
(безрозмірний), A = опорна площа лобового перерізу (м²) — зазвичай πd²/4
Сповільнення через опір: a_drag = F_D / m = (ρ · v² · C_D · A) / (2m)
Число Рейнольдса визначає режим течії: Re = ρvd/μ де μ = динамічна
в'язкість ≈ 1.81×10⁻⁵ Па·с C_D сфери: Re < 1000: C_D ≈ 24/Re +
6/(1+√Re) + 0.4 (Шиллер-Науман) Re = 1000–2×10⁵: C_D ≈ 0.44
(домінує інерція) Re > 2×10⁵: C_D ≈ 0.10 (турбулентний примежовий шар, криза опору
— напр. ямки на м'ячі для гольфу)
2. Залежний від Маха C_D та трансзвукове зростання опору
Число Маха: M = v / a_sound де a = √(γRT) ≈ 340 м/с за 15°C Режим C_D
(обтічний снаряд): M < 0.8 (дозвуковий): C_D ≈ 0.15–0.25 M
0.8–1.2 (трансзвуковий): C_D піковий на 0.4–0.6 (хвильовий опір) M > 1.2
(надзвуковий): C_D ≈ 0.35–0.55 і спадає з M Коефіцієнт опору залежно від
Маха — емпіричний кусковий (модель кулі G7): M < 0.8: C_D = 0.18 M
0.8–1.0: C_D = 0.18 + 1.2(M−0.8)² (швидке зростання поблизу M=1) M 1.0–1.2: C_D
= 0.42 − 0.5(M−1.0) M > 1.2: C_D = 0.32 − 0.03(M−1.2) Джерело
хвильового опору: ударна хвиля формується на носі та основі тіла Коефіцієнт хвильового опору ≈
4/(γM²) · sin²θ_shock (кут конуса θ) → Гостріший ніс → слабша ударна хвиля →
менший хвильовий опір → конструкції з конічним хвостовиком (boat-tail)
3. Балістичний коефіцієнт та поперечна щільність
Балістичний коефіцієнт (BC): BC = m / (C_D · A) [кг/м²] Вище BC →
менше сповільнення опором за тієї самої початкової швидкості Важчий, вужчий, більш
аеродинамічний снаряд → вище BC Поперечна щільність (SD): SD = m / d²
(де d = калібр у дюймах, m у гранах) SD пов'язує масу з
перерізом; кулі з плоским носом мають нижче BC, ніж загострені (spitzer)
конструкції з тією ж SD Еталонні моделі G1 та G7: G1: еталонний снаряд із
плоским дном і тупим носом (історична артилерія) G7: загострений з конічним хвостовиком для великих
дальностей (сучасні гвинтівкові кулі з високим BC) BC_G7 ≈ 0.3–0.65 для
типових далекобійних куль BC_G1 ≈ 0.45–0.90 для тих самих куль
(інша еталонна крива) Час польоту / збереження швидкості: v(t)
≈ v₀ · exp(−(ρ·A·C_D)/(2m) · v₀t) (наближення сталого CD)
4. Стандартна атмосфера ISA
Тропосфера ISA (0–11 км): T = 288.15 − 6.5·h/1000 [K], h у м
Тиск: P(h) = 101325 · (T/288.15)^5.2561 [Па] Густина: ρ(h) = P(h)
/ (287.05 · T) [кг/м³] ρ(0) = 1.225 кг/м³ (рівень моря) ρ(5500м) ≈ 0.736
кг/м³ (≈ 60% від рівня моря) Швидкість звуку: a = √(γRT), γ=1.4, R=287
Дж/(кг·К) a(0°C) = 331.3 м/с, a(15°C) = 340.3 м/с, a(−50°C) = 299 м/с
Вплив великої висоти на дальність: нижче ρ → менше опору На висоті 1500 м
(5000 фт): ρ ≈ 1.05 кг/м³ (−14%) → куля довше зберігає швидкість →
~3% більша дальність Стрільці на великій висоті мусять коригувати
точки влучання вгору
5. Рівняння точкової траєкторії
Вектор стану: [x, y, z, vx, vy, vz] Компоненти прискорення: a_x =
−(F_D/m)·vx/v + F_M_x/m + F_Cor_x/m a_y = −(F_D/m)·vy/v + F_M_y/m +
F_Cor_y/m − g a_z = −(F_D/m)·vz/v + F_M_z/m + F_Cor_z/m F_D =
½ρv²C_D(M)A (напрямок: протилежний v) g = 9.80665 м/с² (змінюється ±0.05%
з широтою) Інтегрування RK4: k₁ = f(t, s) k₂ = f(t+h/2, s+h/2·k₁)
k₃ = f(t+h/2, s+h/2·k₂) k₄ = f(t+h, s+h·k₃) s(t+h) = s(t) +
h·(k₁+2k₂+2k₃+k₄)/6 Крок часу: dt = 0.001–0.01 с (0.001 с для
надзвукових снарядів)
6. Деривація, Магнус та Коріоліс
Ефект Магнуса (обертовий снаряд): F_M = C_M · (ρ/2) · v · A · (ω ×
v̂) Гвинтівка з правим нарізом: вісь обертання ≈ вісь швидкості → гіроскопічна
прецесія спричиняє повільне рискання → права деривація на дальності Наближення
деривації: SD(дальність) ≈ 1.25 · (1 + 4/stability) · twist_inches⁻¹ ·
range² [дюйми/100 ярдів²] Відхилення Коріоліса (по кожній осі): F_Cor = −2m(Ω ×
v) Ω = 7.292×10⁻⁵ рад/с (обертання Землі) Ефект на 1000 м
(горизонтально): ≈ 4–8 см праворуч у Пн. півкулі на середніх широтах Вертикальний
ефект Етвеша: +v_E·2Ω·sin(lat) — впливає на підвищення на будь-якій дальності Аеродинамічний
стрибок: початкове рискання кулі від пориву вітру → бічне відхилення На 1000 м
бічний вітер 10 миль/год → 1–2 MOA (≈28–55 см) убік
7. JavaScript-симулятор балістичної траєкторії
// Точкова балістична траєкторія з атмосферою ISA + моделлю Mach-CD
function isaAtmosphere(h) {
const T = Math.max(216.65, 288.15 - 0.0065*h);
const P = 101325 * Math.pow(T/288.15, 5.2561);
const rho = P / (287.05 * T);
const a = Math.sqrt(1.4 * 287.05 * T);
return {rho, a}; // густина кг/м³, швидкість звуку м/с
}
function dragCoeff(Mach) {
if (Mach < 0.8) return 0.18;
if (Mach < 1.0) return 0.18 + 1.2*(Mach-0.8)**2;
if (Mach < 1.2) return 0.42 - 0.5*(Mach-1.0);
return Math.max(0.15, 0.32 - 0.03*(Mach-1.2));
}
function ballisticTrajectory({
v0, angleDeg, mass, diameter, dt = 0.002
}) {
const angleRad = angleDeg * Math.PI / 180;
const A = Math.PI * (diameter/2)**2;
const g = 9.80665;
let x=0, y=0;
let vx = v0*Math.cos(angleRad), vy = v0*Math.sin(angleRad);
const path = [{x, y}];
while (y >= 0) {
const v = Math.hypot(vx, vy);
const {rho, a} = isaAtmosphere(y);
const Mach = v / a;
const Cd = dragCoeff(Mach);
const Fd = 0.5 * rho * v*v * Cd * A;
// Підкрок RK4 (тут спрощений метод Ейлера для стислості)
const ax = -(Fd/mass) * (vx/v);
const ay = -(Fd/mass) * (vy/v) - g;
vx += ax * dt;
vy += ay * dt;
x += vx * dt;
y += vy * dt;
if (path.length % 50 === 0) path.push({x, y});
}
return path;
}
// Приклад: 7.62×51 NATO на 800 м/с, кут підвищення 30°
const traj = ballisticTrajectory({
v0: 800, angleDeg: 30, mass: 0.0104, diameter: 0.00762
});
console.log(`Range: ${traj[traj.length-1].x.toFixed(0)} m`);
// Вакуум 30° → 32.7 км; з опором → значно менше
8. Оптимізація дальності та реальні поправки
Оптимальний кут запуску у вакуумі: θ_opt = 45° (макс. дальність = v₀²/g) З
опором: θ_opt падає нижче 45° Снаряд із малим опором (M=3): θ_opt ≈ 40°
Великий опір (мушкетна куля): θ_opt ≈ 30–35° Причина: за великих кутів зайвий
час польоту → більше роботи опору Зведення артилерійських поправок:
Початкова швидкість: −5 м/с → поправка дальності −3% Температура −10°C: повітря
щільніше → −1% дальності Боковий вітер 5 м/с: бічне відхилення ≈ 0.1 мрад ×
дальність Висота +1000 м: +0.5% дальності (нижча густина) Коріоліс на 45° пн. ш.:
дрейф праворуч ~4 см на 1000 м (стрілець у Пн. півкулі) Високоточна гвинтівка: Сучасне
ПЗ зовнішньої балістики (Applied Ballistics, Hornady 4DOF)
чисельно інтегрує повні рівняння 6DOF, зокрема гіроскопічну
стабільність, кут спокою при рисканні, Коріоліс/Етвеш по всіх осях, щоб передбачити
точку влучання з точністю до сантиметра на дальностях 1–4 км.