Модель Имитации Долговременных Адаптаций проф. В.Н. Селуянова
Девять ОДУ, три режима работы мышц, прогноз адаптации на любой горизонт от 1 до 365 дней — с учётом рациона и режима питания.
1. Введение
МИДА (Модель Имитации Долговременных Адаптаций) — компьютерная математическая модель долговременной адаптации организма к физической нагрузке, разработанная проф. В.Н. Селуяновым и адаптированная для экосистемы Samson-fit.
Модель — научное ядро сервиса прогноза адаптации: она преобразует тренировочный план и данные рациона в траектории структурной адаптации — массу миофибрилл (mf) и массу митохондрий (mx).
Горизонт прогноза: по умолчанию 90 дней. Пользователь задаёт любой горизонт от 1 до 365 дней с точностью до одного дня; Оптимизатор работает с горизонтами 30–365 дней.
Поддерживаются все три режима работы мышц: изометрический (нативный для модели), динамический (концентрический) и эксцентрический.
2. Теоретические основы
2.1. Принцип суперкомпенсации
Тренировочная нагрузка разрушает часть клеточных структур (миофибриллы, митохондрии). Восстановление не только возвращает исходный уровень, но может превысить его — суперкомпенсация. Повторение нагрузки на пике суперкомпенсации обеспечивает структурный рост; слишком ранняя или слишком поздняя нагрузка накапливает усталость вместо роста.
2.2. Принцип дозирования
МИДА формализует это правило: стимул роста зависит от доли интенсивности выше стрессового порога (ПАНО для митохондрий, около 30% ПМВ для миофибрилл), от продолжительности сессии и от интервала восстановления между сессиями. Слишком высокая интенсивность без достаточного восстановления разрушает структуры.
2.3. Характеристики нагрузки упражнения (формулы (2.7)–(2.11))
Исходная модель Селуянова рассматривает влияние упражнения как простую функцию Prel × t × n. Однако одно и то же упражнение, выполняемое с разной техникой, может давать значительно разные тренировочные эффекты. Реальные упражнения обладают биомеханическими характеристиками, которые существенно модифицируют тренировочный стимул.
Например, приседания со штангой и жим лёжа при одинаковых Prel, t и n по-разному влияют на одни и те же группы мышц: приседания задействуют больше мышечных волокон за счёт большей амплитуды движения, в то время как жим лёжа локализует нагрузку на грудных и трицепсах. Учёт этих различий требует введения поправочных коэффициентов.
Расширенная функция влияния упражнения:
где Kbiomech,j,i — биомеханический поправочный коэффициент упражнения j на систему i (зависит от углов в суставах, плеч силы, рычажных отношений). Диапазон значений: от 0 до 1.
Kфаза,j,i(φ) — поправочный коэффициент фазы на фазе φ упражнения j для системы i. Он учитывает тот факт, что разные мышцы вносят разный вклад на разных фазах упражнения.
Kтемп,j — поправочный коэффициент темпа (зависит от длительности эксцентрической, изометрической и концентрической фаз).
Kампл,j — поправочный коэффициент амплитуды/диапазона движения.
Биомеханический коэффициент Kбиомех
Биомеханический коэффициент определяется отношением момента силы мышцы в данном упражнении к максимально возможному моменту и угловой зависимостью:
где M_j,i — момент силы мышцы i в упражнении j, M_i,макс — максимально возможный момент для мышцы i, theta_j,i — рабочий угол в суставе, theta_i,опт — оптимальный угол для мышцы i (угол, при котором мышца развивает максимальную силу).
Функция cos(θj,i − θi,opt) учитывает тот факт, что сила мышцы зависит от её длины, которая, в свою очередь, определяется углом в суставе. При оптимальном угле (коэффициент cos = 1) мышца развивает максимальную силу; отклонение от оптимума снижает силу пропорционально cos(δθ).
Фазовый коэффициент Kфаза
где Cj,i(φ) — относительный вклад мышцы i на фазе φ упражнения j (из анатомических таблиц), Cj,i,макс — максимальный вклад по всем фазам упражнения.
Практическое значение: во время приседания квадрицепсы доминируют в фазе подъёма (разгибание колена), в то время как ягодицы больше задействованы в фазе опускания (разгибание бедра). Kphase позволяет дифференцировать вклад мышц на каждой фазе упражнения.
Коэффициент темпа Kтемп
где tecc — длительность эксцентрической фазы (с), tiso — длительность изометрической фазы (с), tcon — длительность концентрической фазы (с), ttot = tecc + tiso + tcon — общее время повторения, αtemp ≈ 1.3 — коэффициент эксцентрического напряжения, βtemp ≈ 1.5 — коэффициент изометрического напряжения.
Эксцентрическая фаза создаёт большее мышечное напряжение (микроповреждения), в то время как изометрическая фаза создаёт больший метаболический стресс (окклюзия). Коэффициент Ktempo учитывает эти различия: пауэрлифтинг (3-0-X) даёт Ktempo ≈ 1,30, гипертрофия (2-1-1) даёт Ktempo ≈ 1,50.
Коэффициент амплитуды Kампл
где ROMj — фактический диапазон движения в упражнении j (градусы), ROMмакс — максимально возможный диапазон для данных суставов.
Упражнения с большей амплитудой задействуют больше мышечных волокон на большей длине, создавая более полный стимул. Полные приседания (ROM ≈ 140°) дают Kampl = 0,93, в то время как жим лёжа (ROM ≈ 60°) даёт только Kampl = 0,67.
3. Математическая структура
Состояние модели — система из девяти обыкновенных дифференциальных уравнений.
| Переменная | Смысл | Масштаб |
|---|---|---|
| mg | масса эндокринной железы (MG/MGmax) | медленные — обновление раз в сутки |
| hg | концентрация гормона в железе (HG/HGmax) | быстрые — шаг 1 с |
| hb | концентрация гормона в крови (HB/HBmax) | быстрые — шаг 1 с |
| hm | концентрация гормона в мышце (HM/HMmax) | быстрые — шаг 1 с |
| mf | масса миофибрилл (MF/MFmax) | медленные — обновление раз в сутки |
| mx | масса митохондрий (MX/MXmax) | медленные — обновление раз в сутки |
| hplus | концентрация ионов водорода (H+) | быстрые — шаг 1 с |
3.1. Три режима работы мышц
- Изометрический (нативный): выше 30% ПМВ сосуды пережаты (окклюзия) — доставка гормона в мышцу (hm) ограничена, H+ накапливается быстрее и выводится медленнее; при расслаблении реактивная гиперемия усиливает вынос метаболитов и доставку гормонов. Метаболический стресс становится самостоятельным анаболическим стимулом.
- Динамический (концентрический): классические коэффициенты модели; усталость внутри подхода снижает фактическую интенсивность последующих повторений.
- Эксцентрический: интенсивность может превышать 100% ПМВ (коэффициент нагрузки 1.30, ограничение около 180% ПМВ); повышенное микроповреждение даёт усиленный стресс-стимул (множитель 1.3).
3.2. Ночной отдых
Во время ночного сна и дальнейшего отдыха (16 часов суммарно) быстрые переменные экспоненциально приближаются к равновесию; суточные итоги закрываются однократным обновлением медленных переменных.
4. Параметры модели
4.1. Базовые коэффициенты модели
| Коэффициент | Значение | Описание |
|---|---|---|
| k1 | 0.50 | Синтез гормона в железе во время упражнения (1/ч) |
| k1,rest | 0.15 | Синтез гормона в покое (1/ч) |
| k2 | 5.00 | Секреция гормона из железы в кровь (1/ч) |
| k3 | 0.12 | Элиминация гормона из крови (печень, 1/ч) |
| k4 | 3.00 | Поглощение гормона мышцей из крови (динамика, 1/ч) |
| k6 | 0.05 | Элиминация гормона из мышцы (1/ч) |
| kgs | 0.060 | Синтез эндокринной железы (1/день) |
| kgd | 0.040 | Деградация эндокринной железы (1/день) |
| kgi | 4.00 | Ингибирование деградации железы гормонами |
| kms | 0.045 | Синтез миофибрилл (1/день) |
| kmd | 0.0020 | Деградация миофибрилл (1/день) |
| kmi | 3.00 | Множитель стрессового стимула для MF |
| kxs | 0.055 | Синтез митохондрий (1/день) |
| kx0 | 0.010 | Базовый синтез митохондрий в покое (1/день) |
| kxd | 0.020 | Деградация митохондрий (1/день) |
| kxh | 0.20 | Повреждение митохондрий H+ |
| kh | 6.00 | Скорость накопления H+ (1/ч) |
| khr | 8.00 | Скорость удаления H+ (1/ч) |
| n | 2.0 | Коэффициент Хилла (нелинейность) |
| Iпри | 0.30 | Интенсивность анаэробного порога |
| stress_threshold | 0.50 | Порог интенсивности для мышечного стресса |
В изометрическом режиме дополнительно изменяются коэффициенты (факторы окклюзии), в эксцентрическом применяется коэффициент нагрузки выше ПМВ. Точные значения находятся в модуле движка и подлежат пер-пользовательской калибровке.
4.2. Личные параметры
Начальное состояние и коэффициенты масштабируются массой тела, возрастной группой, полом, типом телосложения и стажем тренировок; личные измерения ПМВ уточняют нормировку интенсивности.
5. Метод расчёта
5.1. Алгоритм расчётного дня
- Активности плана интегрируются шагом 1 с: шаги упражнений обновляют быстрые переменные; накапливаются интегральные дозы гормонов, H+, стресса и аэробного стимула.
- Защита интенсивности: рацион-модификатор k-Energy дня масштабирует фактическую интенсивность сессий, не поднимая её выше 100% ПМВ.
- Дневной отдых между концом тренировки и ночным сном моделируется шагами по 1 минуте.
- Ночной сон и дальнейший отдых (16 часов суммарно) закрываются аналитически; скорость выведения H+ умножается на k-Recovery дня.
- Медленные переменные обновляются однократно: эндокринная железа, миофибриллы (масштабирование k-Anabolic) и митохондрии (×0.3 рационного отклонения).
5.2. Результаты симуляции
Выход включает итоговые массы (mf/mx/mg), посуточные траектории со значениями модификаторов дня, покомпонентную (по мышцам) симуляцию каждого упражнения, пик H+ и интегральные дозы дня. Результаты питают страницу тренировки, Оптимизатор, ИИ-Тренера и отчёты точности.
6. Питание и прогноз адаптации
С этапа интеграции питания прогноз МИДА учитывает не только тренировочную нагрузку, но и рацион и режим питания. Система из девяти уравнений Селуянова остаётся без изменений — рацион входит через модификаторы дня: k-Anabolic (белок на кг массы, распределение по приёмам, энергобаланс) масштабирует рост миофибрилл, k-Energy (энергобаланс, оценка гликогена) корректирует фактическую интенсивность сессий дня, k-Recovery (вечерний белок, достаточность энергии) ускоряет выведение H+ в 16-часовом ночном блоке.
Доступны два источника рациона: рекомендуемые цели фазы (по умолчанию) или план питания пользователя на горизонт прогноза (1–365 дней). Сильный дефицит энергии не усиливает стимул — оптимизатор снижает объём (число подходов), сохраняя интенсивность, а ИИ-Тренер создаёт корректировку рациона при расхождении факта с целями фазы.
6.1. Модификаторы дня
k-Anabolic определяется белком на кг массы тела, распределением белка по приёмам (30–40 г на приём) и неотрицательным энергобалансом; k-Energy — оценкой гликогена (доля углеводов от цели); k-Recovery — вечерним белком и достаточностью энергии.
6.2. Источники рациона
Цели фазы (по умолчанию): прогноз считается при рекомендуемом рационе фазы — все коэффициенты дня равны 1.0.
План питания: пользователь задаёт посуточный/понедельный рацион на горизонт прогноза (1–365 дней). Дни ниже безопасного минимума 1200 ккал заменяются целями фазы; план короче горизонта циклится. Редактор предзаполняет план из целей фазы («Создать из целей фазы») и проверяет JSON на лету.
6.3. Замыкание контура питания
Оптимизатор: при сегодняшнем энергобалансе −500 ккал/сут или ниже объём оптимизированных планов снижается (×0.80, при −1000 ккал/сут ×0.60) за счёт сокращения числа подходов — интенсивность сохраняется.
ИИ-Тренер: при расхождении факта с целями фазы три дня подряд (белок ниже 70% или калории ниже 85% цели) создаётся корректировка рациона с конкретными цифрами.
6.4. Калибровка
Верифицированные записи симуляций MF (с включённым учётом рациона) дают пер-пользовательский масштаб k-Anabolic: среднее отношение факт/прогноз сохраняется в границах 0.85–1.15 и хранится в профиле питания.
7. Оптимизатор и ИИ-Тренер
Оптимизатор подбирает недельные тренировочные планы: генерация кандидатов → surrogate-фильтрация → точный расчёт МИДА (с рационным профилем и периодизацией) → целевая функция по цели → топ-5 планов. Режимы: строгий, лояльный, максимальный; горизонты 30–365 дней; применяются пер-пользовательская калибровка ОДУ и личные ПМВ.
ИИ-Тренер связывает прогнозы с ежедневной практикой: автоматические корректировки по светофору здоровья, плато, перетренированности и рациону; эскалация к физическому тренеру; комбо-режим.
8. Верификация и точность
Каждый прогноз можно верифицировать: платформа сравнивает предсказанные массы с реальными измерениями (сила хвата, индекс BOLT, состав тела, HRV) и фиксирует ошибку. Верифицированные записи питают пер-пользовательскую калибровку и обучающие датасеты ИИ; динамика точности публикуется на странице точности.
9. Ограничения модели
- Модель — упрощение: она не заменяет консультацию врача и не является индивидуальным предписанием.
- Влияние рациона моделируется модификаторами дня с консервативной калибровкой; индивидуальный метаболический ответ может отличаться.
- Качество прогноза зависит от качества данных: честный дневник, точные веса и интенсивности.
- Рационные коэффициенты k-Energy/k-Recovery калибруются по мере накопления статистики верификации (v2).
10. Заключение
МИДА в Samson-fit сочетает научную школу проф. В.Н. Селуянова с программной верификацией: три режима сокращения, персональный горизонт прогноза от 1 до 365 дней, рацион как вход первого класса и автоматический контур корректировок через Оптимизатор и ИИ-Тренер.
Ссылки
- Селуянов В.Н. Модель имитации долгосрочной адаптации (МИДА). — М.: РГАФК, 1995.
- Селуянов В.Н. Подготовка бегуна на средние дистанции. — М.: СпортАкадемПресс, 2001.
- Селуянов В.Н., Сарсания С.К., Сарсания К.С. Ферментативная миопатия и миофибриллярная гипертрофия // Теория и практика физической культуры. — 2004. — № 6.
- Вирю А.А. Гормональные механизмы адаптации к физическим нагрузкам. — Л.: Наука, 1981.
- Меерсон Ф.З. Адаптация, стресс и профилактика. — М.: Наука, 1981.
- МакАрдл У.Д., Кэтч Ф.И., Кэтч В.Л. Физиология спорта: энергетика, питание и работоспособность человека. — 7-е изд. — Lippincott Williams & Wilkins, 2010.
- Зациорский В.М., Крамер У.Дж. Наука и практика силовой тренировки. — 2-е изд. — Human Kinetics, 2006.
- Крамер У.Дж., Ратамесс Н.А. Гормональные реакции и адаптации к силовым упражнениям и тренировкам // Sports Medicine. — 2005. — Т. 35(4). — С. 339–361.
- Шёнфельд Б.Дж. Механизмы мышечной гипертрофии и их применение к силовым тренировкам // Journal of Strength and Conditioning Research. — 2010. — Т. 24(10). — С. 2857–2872.
- Такарада Ю. и др. Влияние силовых упражнений в сочетании с умеренной сосудистой окклюзией на мышечную функцию человека // Journal of Applied Physiology. — 2000. — Т. 88(6). — С. 2097–2106.
Приложение А. Сводная таблица постоянных времени
| Переменная | Процесс | Постоянная времени | Т½ (дни) |
|---|---|---|---|
| mg | Синтез железы | 1/kgs = 16.7 д. | 11.6 |
| mg | Деградация железы | 1/kgd = 25.0 д. | 17.3 |
| mf | Синтез миофибрилл | 1/kms = 22.2 д. | 15.4 |
| mf | Деградация MF | 1/kmd = 500 д. | 346 |
| mx | Синтез митохондрий | 1/kxs = 18.2 д. | 12.6 |
| mx | Деградация MX | 1/kxd = 50 д. | 34.7 |
| hg | Синтез гормона (упражнение) | 1/k1= 2 ч | 1.4 |
| hg | Синтез гормона (покой) | 1/k1,rest = 6.7 ч | 4.6 |
| hb | Элиминация из крови | 1/k3= 8,3 ч | 5.8 |
| hm | Поглощение мышцей | 1/k4= 0,33 ч | 0.23 |
| hm | Элиминация из мышцы | 1/k6= 20 ч | 14 |
| H+ | Накопление | 1/kh= 0,17 ч | 0.12 |
| H+ | Удаление (динамика) | 1/khr = 0.125 ч | 0.087 |
| H+ | Удаление (окклюзия) | 1/khr,iso = 0.5 ч | 0.35 |
Приложение B. Псевдокод алгоритма дня
Упрощённый псевдокод текущего движка (apps/mida/engine.py):
for day in 1..num_days:
nd = nutrition[(day - 1) % len(nutrition)] # рационные коэффициенты дня
for activity in day_plan:
I_eff = clip(activity.intensity * day_factor * nd.k_energy, MVC_limit)
step(state, I_eff, dt=1s) # быстрые переменные (hg/hb/hm/H+)
accumulate integrals (hb, H+, stress, aero)
rest_steps(state, daytime_rest_hours, dt=1min)
rest_decay(state, 16h, recovery_mult=nd.k_recovery) # ночь — аналитически
slow_update(state, doses, k_anabolic=nd.k_anabolic) # mg, mf, mx раз в сутки
Приложение В. Биомеханические коэффициенты для основных упражнений
Данное приложение содержит расширенную таблицу биомеханических коэффициентов для основных силовых упражнений. Значения Kbiomech рассчитываются по формуле (2.8)... Значения Ktempo рассчитываются по формуле (2.10) для стандартного темпа «гипертрофия (2-1-1)», если не указано иное.
Таблица 4: Биомеханические коэффициенты для основных упражнений
| Упражнение | Целевая мышца | Kбиомех | θопт (°) | ROM (°) | Kампл |
|---|---|---|---|---|---|
| Приседания со штангой | Квадрицепс | 0.92 | 90 | 140 | 0.93 |
| Приседания со штангой | Ягодицы | 0.85 | 90 | 140 | 0.93 |
| Жим штанги лежа | Большая грудная | 0.88 | 90 | 60 | 0.67 |
| Жим штанги лежа | Трицепс | 0.72 | 90 | 60 | 0.67 |
| Становая тяга | Разгибатели спины | 0.90 | 45 | 90 | 0.75 |
| Становая тяга | Ягодицы | 0.78 | 45 | 90 | 0.75 |
| Подтягивания | Широчайшая | 0.90 | 120 | 120 | 0.80 |
| Отжимания | Большая грудная | 0.75 | 90 | 90 | 0.60 |
| Выпады | Квадрицепс | 0.80 | 90 | 100 | 0.67 |
Таблица 5: Коэффициенты темпа
| Режим | tecc(с) | tiso(с) | tcon(с) | Kтемп |
|---|---|---|---|---|
| Пауэрлифтинг (3-0-X) | 3 | 0 | X (быстро) | 1.30 |
| Гипертрофия (2-1-1) | 2 | 1 | 1 | 1.50 |
| Скорость (X-0-1) | X (быстро) | 0 | 1 | 0.70 |
| Изометрия (0-5-0) | 0 | 5 | 0 | 1.50 |
Таблица 6: Фазовые коэффициенты для приседания
| Фаза упражнения | Угол (°) | Kфаза | Доля времени |
|---|---|---|---|
| Опускание (эксцентрика) | 0–50 | 0.60 | 35% |
| Нижняя точка (присед) | 50–90 | 0.85 | 10% |
| Подъём (концентрика) | 90–100 | 1.00 | 45% |
| Верхняя точка (фиксация) | 100–140 | 0.75 | 10% |
Таблица 8: Расширенная таблица биомеханических коэффициентов
| Упражнение | Целевая мышца | Kбиомех | θопт (°) | ROM (°) | Kампл | Kтемп |
|---|---|---|---|---|---|---|
| Приседания со штангой | Квадрицепс | 0.92 | 90 | 140 | 0.93 | 1.30 |
| Приседания со штангой | Ягодицы | 0.85 | 90 | 140 | 0.93 | 1.30 |
| Приседания со штангой | Разгибатели спины | 0.70 | 45 | 140 | 0.93 | 1.30 |
| Жим штанги лежа | Большая грудная | 0.88 | 90 | 60 | 0.67 | 1.30 |
| Жим штанги лежа | Трицепс | 0.72 | 90 | 60 | 0.67 | 1.30 |
| Жим штанги лежа | Дельтовидная (пер.) | 0.60 | 90 | 60 | 0.67 | 1.30 |
| Становая тяга | Разгибатели спины | 0.90 | 45 | 90 | 0.75 | 1.30 |
| Становая тяга | Ягодицы | 0.78 | 45 | 90 | 0.75 | 1.30 |
| Становая тяга | Квадрицепс | 0.65 | 90 | 90 | 0.75 | 1.30 |
| Подтягивания | Широчайшая | 0.90 | 120 | 120 | 0.80 | 1.50 |
| Подтягивания | Бицепс | 0.55 | 90 | 120 | 0.80 | 1.50 |
| Отжимания | Большая грудная | 0.75 | 90 | 90 | 0.60 | 1.50 |
| Отжимания | Трицепс | 0.65 | 90 | 90 | 0.60 | 1.50 |
| Выпады | Квадрицепс | 0.80 | 90 | 100 | 0.67 | 1.30 |
| Выпады | Ягодицы | 0.75 | 90 | 100 | 0.67 | 1.30 |
| Жим ногами | Квадрицепс | 0.80 | 90 | 80 | 0.53 | 1.30 |
| Жим ногами | Ягодицы | 0.65 | 90 | 80 | 0.53 | 1.30 |
| Подъёмы гантелей | Дельтовидная (ср.) | 0.70 | 90 | 90 | 0.60 | 1.50 |
| Сгибания рук с гантелями | Бицепс | 0.85 | 90 | 140 | 0.93 | 1.50 |
| Разгибания ног | Квадрицепс | 0.88 | 90 | 90 | 0.60 | 1.30 |
| Тяга в наклоне | Широчайшая | 0.82 | 120 | 80 | 0.53 | 1.30 |
| Тяга в наклоне | Разгибатели спины | 0.78 | 45 | 80 | 0.53 | 1.30 |
Таблица 9: Коэффициенты темпа для различных режимов
Параметры расчёта: αtemp = 1,3 (коэффициент эксцентрического напряжения), βtemp = 1,5 (коэффициент изометрического напряжения).
| Режим | tecc(с) | tiso(с) | tcon(с) | ttot(с) | Kтемп |
|---|---|---|---|---|---|
| Пауэрлифтинг (3-0-X) | 3.0 | 0.0 | 0.5 | 3.5 | 1.30 |
| Гипертрофия (2-1-1) | 2.0 | 1.0 | 1.0 | 4.0 | 1.50 |
| Скорость (X-0-1) | 0.5 | 0.0 | 1.0 | 1.5 | 0.70 |
| Изометрия (0-5-0) | 0.0 | 5.0 | 0.0 | 5.0 | 1.50 |
| Супер-медленно (4-2-4) | 4.0 | 2.0 | 4.0 | 10.0 | 1.52 |
| Взрывной (X-0-X) | 0.5 | 0.0 | 0.5 | 1.0 | 0.65 |
| Контроль (2-0-2) | 2.0 | 0.0 | 2.0 | 4.0 | 1.15 |
Примечания к таблицам
- Kбиомех рассчитывается по формуле (2.8) с использованием анатомических данных о моментах силы мышц при различных углах в суставах.
- Kампл рассчитывается по формуле (2.11). Максимальные значения ROM взяты из анатомических справочников: колено — 150°, плечо — 180°, тазобедренный — 150°, локтевой — 150°.
- Kтемп приведён для стандартного темпа «гипертрофия (2-1-1)», если не указано иное. Для расчёта Ktempo для других темпов используйте формулу (2.10).
- Значения Kbiomech для вспомогательных мышц (антагонистов, стабилизаторов) обычно ниже, чем для целевых мышц, и могут быть уточнены с использованием данных ЭМГ.
- «X» означает взрывное (максимально быстрое) выполнение фазы, tcon ≈ 0,5 с для практических расчётов.