Кейс: медный провод нагревается под воздействием тока — проанализируйте, как изменение температуры влияет на сопротивление, механические свойства и долговечность, и какие расчёты нужны для оценки теплового режима в реальном устройстве
Кратко и по делу — влияние температуры и необходимые расчёты. 1) Сопротивление - Температура влияет почти линейно при умеренных TTT: R(T)=R0 [1+α (T−T0)],
R(T)=R_0\,[1+\alpha\,(T-T_0)], R(T)=R0[1+α(T−T0)],
где для меди α≈3.9×10−3 K−1\alpha\approx 3.9\times10^{-3}\ \mathrm{K^{-1}}α≈3.9×10−3K−1 (около 20∘C20^\circ\mathrm{C}20∘C). Удельное сопротивление при 20∘C20^\circ\mathrm{C}20∘C: ρ20≈1.68×10−8 Ω⋅m\rho_{20}\approx1.68\times10^{-8}\ \Omega\cdot\mathrm{m}ρ20≈1.68×10−8Ω⋅m. - Важный момент: при заданном токе III мощность выделяется как P=I2R(T).
P=I^2R(T). P=I2R(T).
Для источника напряжения VVV мощность P=V2/R(T)P=V^2/R(T)P=V2/R(T) — знак обратной связи меняется (при росте TTT и RRR мощность уменьшается). - На высоких температурах зависимость нестрого-линейна; при приближении к плавлению (Tm≈1084.6∘CT_m\approx1084.6^\circ\mathrm{C}Tm≈1084.6∘C) свойства резко меняются. 2) Механические свойства и долговечность - Модуль упругости EEE снижается с температурой, прочность и предел текучести заметно падают. Типично уже при T>0.3…0.4 TmT>0.3\ldots0.4\ T_mT>0.3…0.4Tm (в К) начинается существенный спад механической прочности и активируется ползучесть. - Для меди: рекристаллизация и отжиг происходят при сотнях градусов (характерно изменение твердости и пластичности при \( \sim200\mbox{–}400^\circ\mathrm{C}\)). - Ползучесть (creep) становится существенной при T≳0.4 TmT\gtrsim 0.4\,T_mT≳0.4Tm (в К). Для оценки долговечности применяют законы ползучести (например, степенные или Аррениус-зависимости скорости ползучести). - Повреждение от термоциклирования (усталость из‑за циклов нагрева/охлаждения) — зависит от амплитуды и числа циклов; используют fatigue-кривые и критерии накопленной пластической деформации. - Коррозия/окисление при высокой температуре и деградация изоляции (если провод в изоляции) ограничивают срок службы. 3) Тепловой режим — какие расчёты нужны - Определить исходные параметры: геометрию провода (L,dL, dL,d), ток/напряжение (III или VVV), материалные свойства как функции TTT: ρ(T), k(T), cp(T), ρmass\rho(T),\ k(T),\ c_p(T),\ \rho_{mass}ρ(T),k(T),cp(T),ρmass, эмиссивность ϵ\epsilonϵ, термическое сопротивление контактов. - Баланс мощности (стационарный): - Для длинного тонкого провода с конвекцией и излучением в окружение: I2R(T)=hAs (Ts−T∞)+ϵσAs (Ts4−T∞4),
I^2R(T)=hA_s\,(T_s-T_\infty)+\epsilon\sigma A_s\,(T_s^4-T_\infty^4), I2R(T)=hAs(Ts−T∞)+ϵσAs(Ts4−T∞4),
где hhh — коэффициент теплоотдачи, AsA_sAs — площадь поверхности провода, TsT_sTs — температура поверхности, T∞T_\inftyT∞ — температура окружения. - В простейшем приближении для равновесия можно пренебречь излучением при невысоких TTT. - Транзиент (если нагрев быстрый): уравнение теплоёмкости ρcpVdTdt=I2R(T)−теплопотери.
\rho c_p V\frac{dT}{dt}=I^2R(T)-\text{теплопотери}. ρcpVdtdT=I2R(T)−теплопотери.
В приближении «одного узла» характерная постоянная времени τ=ρcpVhAs.
\tau=\frac{\rho c_p V}{hA_s}. τ=hAsρcpV.
Решение: T(t)=T∞+(T(0)−T∞)e−t/τ+I2RhAs(1−e−t/τ)T(t)=T_\infty+(T(0)-T_\infty)\mathrm{e}^{-t/\tau}+\frac{I^2R}{hA_s}(1-\mathrm{e}^{-t/\tau})T(t)=T∞+(T(0)−T∞)e−t/τ+hAsI2R(1−e−t/τ) (при слабой зависимости R(T)R(T)R(T)). - Простая моделировка поперечного распределения температуры требует решения уравнения теплопроводности ∇⋅(k∇T)+q=ρcp∂T/∂t\nabla\cdot(k\nabla T)+q= \rho c_p \partial T/\partial t∇⋅(k∇T)+q=ρcp∂T/∂t. - Оценка уместности одноузловой модели: число Био Bi=hLc/kBi=hL_c/kBi=hLc/k. Для цилиндра Lc≈r/2L_c\approx r/2Lc≈r/2. Если Bi<0.1Bi<0.1Bi<0.1, можно считать температуру однородной по сечению. - Для переменного тока учитывать скин‑эффект: глубина скина δ=2ρωμ,
\delta=\sqrt{\frac{2\rho}{\omega\mu}}, δ=ωμ2ρ,
что влияет на эффективное сопротивление при высоких частот. 4) Что сравнивать и проверять - Максимальную температуру TmaxT_\text{max}Tmax с лимитами: рекристаллизация, рабочая температура изоляции, температурный класс материала, температуры начала существенного ползучести. - Термальные градиенты → механические напряжения (термонапряжения), особенно в закреплённых проводах. - Долговечность: вычислить скорость ползучести или время до достижения критической деформации (используются эмпирические/модельные зависимости, часто вида ε˙=Aσnexp(−Q/(RT))\dot\varepsilon = A\sigma^n\exp(-Q/(RT))ε˙=Aσnexp(−Q/(RT))). - При циклическом нагреве оценить число циклов до разрушения по соответствующим методикам усталости. 5) Практические данные для меди (примерные числа) - ρ20≈1.68×10−8 Ω⋅m\rho_{20}\approx1.68\times10^{-8}\ \Omega\cdot\mathrm{m}ρ20≈1.68×10−8Ω⋅m, α≈3.9×10−3 K−1\alpha\approx3.9\times10^{-3}\ \mathrm{K^{-1}}α≈3.9×10−3K−1. - Плотность ρmass≈8960 kg/m3\rho_{\text{mass}}\approx8960\ \mathrm{kg/m^3}ρmass≈8960kg/m3, cp≈385 J/(kg⋅K)c_p\approx385\ \mathrm{J/(kg\cdot K)}cp≈385J/(kg⋅K), теплопроводность \(k\approx 380\mbox{–}400\ \mathrm{W/(m\cdot K)}\) при 20∘C20^\circ\mathrm{C}20∘C. - Плавление: Tm≈1084.6∘CT_m\approx1084.6^\circ\mathrm{C}Tm≈1084.6∘C. 6) Рекомендуемая последовательность расчётов для практического случая - Снять геометрию и режим (I или V, частота). - Получить начальные R при T0T_0T0 и формулу R(T)R(T)R(T). - Оценить стационарную TTT решая баланс мощности (включая конвекцию/излучение/теплопроводность в креплениях). - Выполнить транзиентный расчёт (если требуются пуски/циклы) или численную модель (FEM) при значительных градиентах. - На основе T(t)T(t)T(t) оценить механические свойства, время до ползучести и усталости. Если нужно, могу: 1) рассчитать пример (с вашими I,d,LI, d, LI,d,L, условиями охлаждения), 2) составить упрощённую модель в Excel/FEM.
1) Сопротивление
- Температура влияет почти линейно при умеренных TTT:
R(T)=R0 [1+α (T−T0)], R(T)=R_0\,[1+\alpha\,(T-T_0)],
R(T)=R0 [1+α(T−T0 )], где для меди α≈3.9×10−3 K−1\alpha\approx 3.9\times10^{-3}\ \mathrm{K^{-1}}α≈3.9×10−3 K−1 (около 20∘C20^\circ\mathrm{C}20∘C). Удельное сопротивление при 20∘C20^\circ\mathrm{C}20∘C: ρ20≈1.68×10−8 Ω⋅m\rho_{20}\approx1.68\times10^{-8}\ \Omega\cdot\mathrm{m}ρ20 ≈1.68×10−8 Ω⋅m.
- Важный момент: при заданном токе III мощность выделяется как
P=I2R(T). P=I^2R(T).
P=I2R(T). Для источника напряжения VVV мощность P=V2/R(T)P=V^2/R(T)P=V2/R(T) — знак обратной связи меняется (при росте TTT и RRR мощность уменьшается).
- На высоких температурах зависимость нестрого-линейна; при приближении к плавлению (Tm≈1084.6∘CT_m\approx1084.6^\circ\mathrm{C}Tm ≈1084.6∘C) свойства резко меняются.
2) Механические свойства и долговечность
- Модуль упругости EEE снижается с температурой, прочность и предел текучести заметно падают. Типично уже при T>0.3…0.4 TmT>0.3\ldots0.4\ T_mT>0.3…0.4 Tm (в К) начинается существенный спад механической прочности и активируется ползучесть.
- Для меди: рекристаллизация и отжиг происходят при сотнях градусов (характерно изменение твердости и пластичности при \( \sim200\mbox{–}400^\circ\mathrm{C}\)).
- Ползучесть (creep) становится существенной при T≳0.4 TmT\gtrsim 0.4\,T_mT≳0.4Tm (в К). Для оценки долговечности применяют законы ползучести (например, степенные или Аррениус-зависимости скорости ползучести).
- Повреждение от термоциклирования (усталость из‑за циклов нагрева/охлаждения) — зависит от амплитуды и числа циклов; используют fatigue-кривые и критерии накопленной пластической деформации.
- Коррозия/окисление при высокой температуре и деградация изоляции (если провод в изоляции) ограничивают срок службы.
3) Тепловой режим — какие расчёты нужны
- Определить исходные параметры: геометрию провода (L,dL, dL,d), ток/напряжение (III или VVV), материалные свойства как функции TTT: ρ(T), k(T), cp(T), ρmass\rho(T),\ k(T),\ c_p(T),\ \rho_{mass}ρ(T), k(T), cp (T), ρmass , эмиссивность ϵ\epsilonϵ, термическое сопротивление контактов.
- Баланс мощности (стационарный):
- Для длинного тонкого провода с конвекцией и излучением в окружение:
I2R(T)=hAs (Ts−T∞)+ϵσAs (Ts4−T∞4), I^2R(T)=hA_s\,(T_s-T_\infty)+\epsilon\sigma A_s\,(T_s^4-T_\infty^4),
I2R(T)=hAs (Ts −T∞ )+ϵσAs (Ts4 −T∞4 ), где hhh — коэффициент теплоотдачи, AsA_sAs — площадь поверхности провода, TsT_sTs — температура поверхности, T∞T_\inftyT∞ — температура окружения.
- В простейшем приближении для равновесия можно пренебречь излучением при невысоких TTT.
- Транзиент (если нагрев быстрый): уравнение теплоёмкости
ρcpVdTdt=I2R(T)−теплопотери. \rho c_p V\frac{dT}{dt}=I^2R(T)-\text{теплопотери}.
ρcp VdtdT =I2R(T)−теплопотери. В приближении «одного узла» характерная постоянная времени
τ=ρcpVhAs. \tau=\frac{\rho c_p V}{hA_s}.
τ=hAs ρcp V . Решение: T(t)=T∞+(T(0)−T∞)e−t/τ+I2RhAs(1−e−t/τ)T(t)=T_\infty+(T(0)-T_\infty)\mathrm{e}^{-t/\tau}+\frac{I^2R}{hA_s}(1-\mathrm{e}^{-t/\tau})T(t)=T∞ +(T(0)−T∞ )e−t/τ+hAs I2R (1−e−t/τ) (при слабой зависимости R(T)R(T)R(T)).
- Простая моделировка поперечного распределения температуры требует решения уравнения теплопроводности ∇⋅(k∇T)+q=ρcp∂T/∂t\nabla\cdot(k\nabla T)+q= \rho c_p \partial T/\partial t∇⋅(k∇T)+q=ρcp ∂T/∂t.
- Оценка уместности одноузловой модели: число Био Bi=hLc/kBi=hL_c/kBi=hLc /k. Для цилиндра Lc≈r/2L_c\approx r/2Lc ≈r/2. Если Bi<0.1Bi<0.1Bi<0.1, можно считать температуру однородной по сечению.
- Для переменного тока учитывать скин‑эффект: глубина скина
δ=2ρωμ, \delta=\sqrt{\frac{2\rho}{\omega\mu}},
δ=ωμ2ρ , что влияет на эффективное сопротивление при высоких частот.
4) Что сравнивать и проверять
- Максимальную температуру TmaxT_\text{max}Tmax с лимитами: рекристаллизация, рабочая температура изоляции, температурный класс материала, температуры начала существенного ползучести.
- Термальные градиенты → механические напряжения (термонапряжения), особенно в закреплённых проводах.
- Долговечность: вычислить скорость ползучести или время до достижения критической деформации (используются эмпирические/модельные зависимости, часто вида ε˙=Aσnexp(−Q/(RT))\dot\varepsilon = A\sigma^n\exp(-Q/(RT))ε˙=Aσnexp(−Q/(RT))).
- При циклическом нагреве оценить число циклов до разрушения по соответствующим методикам усталости.
5) Практические данные для меди (примерные числа)
- ρ20≈1.68×10−8 Ω⋅m\rho_{20}\approx1.68\times10^{-8}\ \Omega\cdot\mathrm{m}ρ20 ≈1.68×10−8 Ω⋅m, α≈3.9×10−3 K−1\alpha\approx3.9\times10^{-3}\ \mathrm{K^{-1}}α≈3.9×10−3 K−1.
- Плотность ρmass≈8960 kg/m3\rho_{\text{mass}}\approx8960\ \mathrm{kg/m^3}ρmass ≈8960 kg/m3, cp≈385 J/(kg⋅K)c_p\approx385\ \mathrm{J/(kg\cdot K)}cp ≈385 J/(kg⋅K), теплопроводность \(k\approx 380\mbox{–}400\ \mathrm{W/(m\cdot K)}\) при 20∘C20^\circ\mathrm{C}20∘C.
- Плавление: Tm≈1084.6∘CT_m\approx1084.6^\circ\mathrm{C}Tm ≈1084.6∘C.
6) Рекомендуемая последовательность расчётов для практического случая
- Снять геометрию и режим (I или V, частота).
- Получить начальные R при T0T_0T0 и формулу R(T)R(T)R(T).
- Оценить стационарную TTT решая баланс мощности (включая конвекцию/излучение/теплопроводность в креплениях).
- Выполнить транзиентный расчёт (если требуются пуски/циклы) или численную модель (FEM) при значительных градиентах.
- На основе T(t)T(t)T(t) оценить механические свойства, время до ползучести и усталости.
Если нужно, могу: 1) рассчитать пример (с вашими I,d,LI, d, LI,d,L, условиями охлаждения), 2) составить упрощённую модель в Excel/FEM.