def thermal_properties(U, r, R, A, rho, t0, du):
"""
计算导热系数和比热容
参数:
U (float): 加热器的电压
r (float): 加热器的电阻
R (float): 样品的厚度
A (float): 样品的横截面积
rho (float): 样品的密度
t0 (float): 样品的初始温度
du (float): 加热面和中心面的温度差
返回:
k (float): 样品的导热系数
c (float): 样品的比热容
"""
# 计算加热功率
q_c = U ** 2 / (2 * r)
# 计算导热系数
k= q_c * R / du
# 计算比热容
c = q_c /(rho*R*(alpha))
return k, c
# 测试数据
U = 10
r = 2
R = 0.1
A = 0.01
rho = 2700
t0 = 20
du = 5
t = 60
alpha=1
k, c = thermal_properties(U, r, R, A, rho, t0, du)
print(f"导热系数 k = {k:.2f} W/(m·K)")
print(f"比热容 c = {c:.2f} J/(kg·K)")