Моделирование цикла Карно в Excel: от теории к практическим расчётам

Введение: зачем моделировать цикл Карно в электронных таблицах?

Цикл Карно — фундаментальная концепция термодинамики, описывающая идеальный тепловой двигатель с максимально возможным КПД. Для инженеров, студентов технических специальностей и исследователей умение визуализировать этот цикл в Microsoft Excel или Google Sheets становится ценным навыком. Электронные таблицы позволяют не только рассчитать ключевые параметры (давление, объём, температуру), но и построить PV-диаграмму — графическое представление процесса, которое наглядно демонстрирует работу цикла.

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

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

Подготовка структуры таблицы: какие данные нужны для расчётов

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

• 🔹 Констант: универсальная газовая постоянная (R = 8.314 Дж/(моль·К)), показатель адиабаты (γ = Cp/Cv, для двухатомного газа ≈1.4).
• 🔹 Входных параметров: температуры нагревателя (T₁) и холодильника (T₂), начальные давление (P₁) и объём (V₁), количество вещества (n).
• 🔹 Промежуточных расчётов: давление и объём в ключевых точках цикла (A→B→C→D→A).
• 🔹 Результатов: работа на каждом этапе, суммарная работа цикла (A_cycle), КПД (η).

Пример организации данных:

⚠️ Внимание: Все температуры в формулах должны быть в кельвинах. Если ваши исходные данные в градусах Цельсия, добавьте столбец для конвертации: =C2+273.15, где C2 — ячейка с °C.

📊 Какой газ вы чаще всего моделируете в термодинамических задачах?

Воздух

Азот

Гелий

Водяной пар

Другой

Расчёт параметров в ключевых точках цикла

Цикл Карно состоит из четырёх процессов, каждый из которых связывает две точки на PV-диаграмме. Для точек A→B→C→D→A необходимо рассчитать давление и объём, используя уравнения состояния идеального газа и законы термодинамики.

1. Точка A (начало изотермического расширения)

Здесь известны P₁, V₁, T₁. Используем уравнение Менделеева-Клапейрона:

=B8  ' V₁ (м³)
=(n*R*B5)/B8  ' P₁ = nRT₁/V₁ (Па)

2. Точка B (конец изотермического расширения)

Объём V₂ выбираем произвольно (например, в 2 раза больше V₁). Давление рассчитываем из изотермического условия P₁V₁ = P₂V₂:

=B8*2  ' V₂ = 2V₁
=(D2*B8)/D3  ' P₂ = P₁V₁/V₂

3. Точка C (конец адиабатического расширения)

Используем уравнение адиабаты P₂V₂^γ = P₃V₃^γ и связь температур T₂/T₁ = (V₂/V₃)^(γ-1). Решаем систему уравнений:

=D3*(B5/B6)^(1/(B3-1))  ' V₃
=D4*B6/D3  ' P₃ = nRT₂/V₃

Точка A: P₁V₁ = nRT₁

Точка B: P₂V₂ = P₁V₁ (изотерма)

Точка C: P₃V₃^γ = P₂V₂^γ (адиабата)

Точка D: T₄ = T₁ (замыкание цикла)-->

Автоматизация расчётов с помощью формул Excel

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

=n*R*B5*LN(D3/B8)  ' Работа A→B = nRT₁ ln(V₂/V₁)

Для адиабатического процесса B→C работа равна изменению внутренней энергии:

=n*B2*(B6-B5)/(B3-1)  ' Работа B→C = nCv(T₂-T₁)

⚠️ Внимание: При копировании формул Excel автоматически сдвигает ссылки на ячейки. Чтобы зафиксировать адрес (например, для константы R), используйте абсолютную ссылку: $B$2.

Построение PV-диаграммы цикла Карно

Визуализация — критически важный этап. PV-диаграмма должна отображать все четыре процесса цикла с чётким разделением на изотермы и адиабаты. Алгоритм:

1. Выделите столбцы с данными P и V для точек A, B, C, D, A (замкните цикл).
2. Перейдите на вкладку Вставка → Вставить график → Точечная с гладкими кривыми.
3. Добавьте подписи осей: горизонтальная — Объём (V), м³, вертикальная — Давление (P), Па.
4. Отформатируйте линии: изотермы (A→B и C→D) — красным, адиабаты (B→C и D→A) — синим.

⚠️ Внимание: Если график получился «рваным» (не замкнулся), проверьте последовательность точек в исходных данных. Последняя точка должна совпадать с первой (A), чтобы кривая замкнулась.

Как добавить вторую ось для температуры?

Чтобы отобразить температуру на PV-диаграмме, создайте вспомогательный столбец с температурами для каждой точки (T₁, T₁, T₂, T₂, T₁). Затем:

1. Кликните правой кнопкой по графику → Выбрать данные.

2. Добавьте новый ряд для температуры.

3. Кликните правой кнопкой по новому ряду → Формат ряда данных → выберите Вторичная ось.

4. Настройте вторую ось (справа) с подписью Температура (T), К.

Расчёт КПД цикла и проверка на ошибки

КПД цикла Карно определяется только температурами нагревателя и холодильника: η = 1 − T₂/T₁. В Excel это реализуется простой формулой:

=1-B6/B5  ' КПД = 1 - T₂/T₁

Однако для проверки корректности модели сравните этот результат с КПД, рассчитанным через работы:

=A_cycle/Q_in  ' A_cycle — суммарная работа, Q_in — теплота, полученная от нагревателя

Разница между двумя значениями КПД не должна превышать 0.1%. Если расхождение больше — ищите ошибку в расчётах промежуточных точек.

Параметр	Формула	Пример значения
Суммарная работа A_cycle	=A_AB + A_BC + A_CD + A_DA	1256.6 Дж
Теплота Q_in	=n*R*T₁*LN(V₂/V₁)	4188.8 Дж
КПД по температурам	=1-T₂/T₁	0.4 (40%)
КПД по работам	=A_cycle/Q_in	0.3999 (39.99%)

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

Чтобы сделать модель интерактивной, добавьте элементы управления:

• 🎚️ Ползунки для температур (T₁ и T₂): Разработчик → Вставить → Ползунок (Form Control). Свяжите ползунок с ячейкой, содержащей температуру.
• 📊 Условное форматирование для КПД: если η > 0.5, закрасьте ячейку зелёным, если η < 0.3 — красным.
• 🔄 Кнопка сброса: назначьте макрос, который вернёт все параметры к начальным значениям.

Пример настройки ползунка:

1. Активируйте вкладку Разработчик (если её нет: Файл → Параметры → Настройка ленты).
2. Вставьте ползунок и кликните по нему правой кнопкой → Формат объекта.
3. В поле Связать с ячейкой укажите адрес, где хранится T₁ (например, $B$5).
4. Задайте диапазон значений (например, минимальное 300, максимальное 1000).

Частые ошибки и как их избежать

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

• 🔥 Неверные единицы измерения: Убедитесь, что давление в паскалях, объём в кубических метрах, температура в кельвинах. Например, если подставить объём в литрах, формулы дадут некорректный результат.
• 🔄 Незамкнутый цикл: Проверьте, что последняя точка (D→A) возвращается к исходным P₁ и V₁. Используйте формулу =И(AND(P_A=P_D; V_A=V_D)) для автоматической проверки.
• 📉 Отрицательная работа: На участках сжатия (C→D и D→A) работа должна быть отрицательной. Если получается положительное значение, поменяйте знак в формуле или проверьте направление процесса.

⚠️ Внимание: Если при изменении температур нагревателя и холодильника КПД цикла не меняется — проверьте абсолютные ссылки на ячейки в формуле =1-B6/B5. Возможно, Excel автоматически сдвинул адреса при копировании.

FAQ: Ответы на частые вопросы

Можно ли использовать эту модель для реальных газов, а не идеальных?

Нет, приведённые формулы основаны на уравнении состояния идеального газа PV = nRT. Для реальных газов (например, водяного пара) требуется учитывать поправки ван-дер-Ваальса или использовать таблицы термодинамических свойств. Однако структуру таблицы в Excel можно адаптировать, заменив уравнения на соответствующие зависимости для реального газа.

Как добавить в модель расчёт энтропии?

Энтропия изменяется только на изотермических участках (A→B и C→D). Для расчёта используйте формулы:

=n*R*LN(V₂/V₁)  ' ΔS_AB (изотерма T₁)
=-n*R*LN(V₃/V₄)  ' ΔS_CD (изотерма T₂)

На адиабатических участках энтропия постоянна (ΔS = 0). Суммарное изменение энтропии за цикл должно быть равно нулю.

Почему график получается не гладким, а ломаным?

Excel соединяет точки на графике прямыми линиями. Чтобы сгладить кривые:

1. Кликните правой кнопкой по линии графика → Формат ряда данных.
2. Выберите Сглаженная линия (опция доступна для точечных графиков).
3. Для более точного отображения увеличьте количество промежуточных точек между A→B, B→C и т.д. (используйте вспомогательные столбцы с шагом 0.1 по объёму).

Как экспортировать данные из Excel в программу для численных расчётов (например, MATLAB)?

Сохраните лист Excel в формате .csv (Файл → Сохранить как → CSV (разделители — запятые)). В MATLAB импортируйте данные с помощью:

data = readmatrix('carnot_cycle.csv');
T1 = data(1,2); V1 = data(2,2); % и т.д.

Убедитесь, что в CSV-файле нет пустых строк или объединённых ячеек — это вызовет ошибки при чтении.

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

Да. Создайте выпадающий список (Данные → Проверка данных → Список) с вариантами газов (например, "Гелий", "Азот", "Воздух"). В соседних столбцах храните соответствующие значения γ и Cv:

|   A   |     B     |   C   |   D   |
|-------|-----------|-------|-------|
| Газ   | γ         | Cv    | M     |
| Гелий | 1.667     | 12.5  | 4     |
| Азот  | 1.4       | 20.8  | 28    |

Затем используйте функцию ВПР для подстановки параметров:

=VLOOKUP(B1; A2:D3; 2; FALSE)  ' Подставляет γ для выбранного газа