小灯泡的动态电阻曲线理论推导

小型白炽灯珠伏安特性理论分析

• 重新测试白炽灯珠的伏安特性

 

01 动态电阻

---

一、前言

  昨天测试了这种小型白炽灯的伏安特性。 可以看到它具有比较特殊的曲线形式。  特别是， 通过计算电流变化与电压变化的比值， 所得到它的动态电阻， 可以看到具有特殊的峰值分布。  只要温度上升和下降足够慢， 得到的 伏安曲线在电压上升和下降过程是一致的。 下面希望通过对小灯泡热耗散理论建模推导出 它对应的伏安特性曲线。

二、理论分析

1、基本假设

  （1） 假设灯丝的电阻与温度成线性关系： $R_a=a\cdot \left(t-t_0\right)+b$ ;
  （2） 假设灯丝热辐射与温度成四次方关系： $P_r=c\cdot {\left(t-t_0\right)}^4$
  （3） 假设灯丝热传导耗散功率与温度成正比： $P_t=d\cdot \left(t-t_0\right)$

首先，给出一些基本假设。 假设电灯泡的钨丝电阻与温度呈线性关系。  灯丝热辐射功率与温度呈四次方关系； 钨丝通过热传导耗散功率与温度呈正比。  根据所测量得到的伏安曲线可以计算出对应的温度。 再根据输入电功率等于辐射功率加上热传导功率， 便可以得到能量平衡方程。 将方程中右边的系数进行简化。  最终得到钨丝在功率守恒情况下的伏安曲线守恒方程。

2、电压电流关系

  在电压 $U$ 作用下， 对应的电流为 $I$ 。 此时输入电功率为： $U\cdot I$， 对应的电阻为 $U/I$ 。 假设对应的 温度为 $t$ ， 根据基本假设 （1）， 可知： $$\frac{U}{I}=a\cdot \left(t-t_0\right)+b$$

  所以： $t=\frac{1}{a}\left(\frac{U}{I}-b\right)+t_0$ ， 根据假设 （2）、（3）， 输入电功率等于耗散功率与辐射功率之和：

$$U\cdot I=\frac{c}{a^4}{\left(\frac{U}{I}-b\right)}^4+\frac{d}{a}\left(\frac{U}{I}-b\right)$$

  将 $c/a^4,d/a$ 分别以 $C,D$ 替换， 那么有：

$$U\cdot I=C\cdot {\left(\frac{U}{I}-b\right)}^4+D\cdot \left(\frac{U}{I}-b\right)$$

3、数据拟合

  使用前面所测量得到的数据， 来拟合钨丝伏安特性能量守恒方程中的系数。 最终确定下其中的 b,C,D 三个参数的数值。  通过这个数值， 可以计算出不同测量电压下的 钨丝辐射功率。 可以看到与输入的电功率曲线基本上是吻合的。

▲ 图1.2.1 实际测试数据

vdim=[0.0159,0.0325,0.0487,0.0656,0.0822,0.0997,0.1172,0.1358,0.1551,0.1746,0.1957,0.2172,0.2405,0.2642,0.2898,0.3154,0.3427,0.3699,0.3984,0.4271,0.4550,0.4838,0.5118,0.5406,0.5683,0.5969,0.6245,0.6529,0.6804,0.7087,0.7370,0.7645,0.7927,0.8200,0.8482,0.8755,0.9037,0.9310,0.9592,0.9865,1.0147,1.0429,1.0702,1.0984,1.1258,1.1540,1.1813,1.2096,1.2369,1.2652,1.2926,1.3209,1.3492,1.3766,1.4049,1.4323,1.4607,1.4882,1.5166,1.5441,1.5726,1.6001,1.6286,1.6571,1.6847,1.7133,1.7409,1.7695,1.7972,1.8258,1.8536,1.8822,1.9100,1.9388,1.9675,1.9953,2.0240,2.0519,2.0806,2.1084,2.1372,2.1650,2.1937,2.2225,2.2503,2.2791,2.3069,2.3356,2.3633,2.3922,2.4202,2.4490,2.4769,2.5057,2.5345,2.5624,2.5912,2.6192,2.6482]
idim=[0.0118,0.0244,0.0365,0.0489,0.0607,0.0725,0.0836,0.0945,0.1047,0.1139,0.1226,0.1300,0.1366,0.1420,0.1467,0.1504,0.1534,0.1557,0.1578,0.1596,0.1612,0.1628,0.1644,0.1662,0.1680,0.1699,0.1719,0.1739,0.1760,0.1781,0.1803,0.1824,0.1847,0.1869,0.1892,0.1914,0.1936,0.1959,0.1982,0.2004,0.2027,0.2050,0.2072,0.2095,0.2116,0.2139,0.2161,0.2183,0.2205,0.2227,0.2248,0.2270,0.2291,0.2312,0.2334,0.2354,0.2375,0.2395,0.2416,0.2436,0.2457,0.2476,0.2497,0.2517,0.2536,0.2556,0.2574,0.2594,0.2612,0.2631,0.2650,0.2668,0.2686,0.2705,0.2723,0.2740,0.2758,0.2776,0.2794,0.2811,0.2829,0.2846,0.2864,0.2881,0.2898,0.2916,0.2933,0.2950,0.2967,0.2984,0.3001,0.3018,0.3035,0.3052,0.3069,0.3085,0.3102,0.3118,0.3135]

$$b=1.192$$$$C=1.854\times 10^{-4}$$$$D=4.62\times 10^{-2}$$

▲ 图1.2.2 绘制UI以及拟合 UI

二、参数计算

  根据前面拟合过程， 对于 温度 $t_0$ 无法进行直接拟合。 下面假设实验室温度 $t_0=28$ 摄氏度。

  参数 $a$ 也无法直接估计， 根据 金属钨的温度系数 $r=4.5\times 10^{-3}$ , 可以知道 ：

$$a=b\cdot r=5.364\times 10^{-3}$$

  那么
$$c=C\cdot a^4=1.5348\times 10^{-13}$$ $$d=D\cdot a=2.478\times 10^{-4}$$

三、拟合曲线

  前面钨丝伏安特性所满足的功率守恒公式中的参数通过实际测量数据拟合求出，  利用这个方程，  可以针对每一个电压 U 求解出对应的 电流 I。 利用 sympy 函数中的 nsolve 数值方程求解工具， 计算实验中每一个电压对应电流。  这样便可以获得模型对应的预测电流。 实际计算结果中， 除了在电压非常小的时候计算出的电流出现负值。  其它电压对应方程求解的电流与实际电流还是非常吻合的。 只是在伏安曲线拐点处， 还是出现了比较大的差异。

▲ 图1.4.1 求解方程搜得到的拟合曲线

#!/usr/local/bin/python
# -*- coding: gbk -*-
#============================================================
# TEST2.PY                     -- by Dr. ZhuoQing 2024-12-01
#
# Note:
#============================================================
from headm import *
from sympy                  import symbols,simplify,expand,print_latex
from sympy                  import *
#------------------------------------------------------------
vdim, idim = tspload('measure', 'vdim', 'idim')
b = 1.192
C = 1.854e-4
D = 4.62e-2
#------------------------------------------------------------
x = Symbol('x')
ifit = []
for U,i in zip(vdim, idim):
    ii = nsolve(C*(U/x-b)**4+D*(U/x-b)-U*x, x, i)
    if ii < 0: ii = i
    printff(U,i,ii)
    ifit.append(ii)
plt.plot(vdim, idim, lw=3, label='Origin')
plt.plot(vdim, ifit, lw=3, label='Fit')
plt.xlabel("Voltage(V)")
plt.ylabel("Current(A)")
plt.grid(True)
plt.legend(loc="upper left")
plt.tight_layout()
plt.show()
#------------------------------------------------------------
result = 'result'
#------------------------------------------------------------
#------------------------------------------------------------
mstr = latex(result)
printf(mstr)
_=tspexecutepythoncmd("msg2latex")
#------------------------------------------------------------
#        END OF FILE : TEST2.PY
#============================================================

 

※ 总  结 ※

---

  本文对于小型白炽灯珠的伏安特性进行了理论分析。 主要根据测量所得的伏安特性数据，  根据基本的物理假设。 利用输入电功率与热传导和热辐射功率平衡，  得到了灯珠的伏安特性平衡方程。 利用测量曲线获得方程中的系数。  通过该方程， 可以求解出电压对应的电流。 可以看到与实际测量的电流变化趋势是一致的。

---

■ 相关文献链接:

• 重新测试白炽灯珠的伏安特性
• 金属钨的温度系数

● 相关图表链接:

• 图1.2.1 实际测试数据
• 图1.2.2 绘制UI以及拟合 UI
• 图1.4.1 求解方程搜得到的拟合曲线