# 计算直线塔不平衡张力
# 新版输电线路大手册 P328

import math

# 架线时的状态
t_m = 100  # 导线架设时的气温。单位°C
sigma_m = 10  # 架线时，初伸长未释放前的最低点水平应力。单位N/mm2
lambda_m = 10  # 导线比载 N/(m.mm)
t_e = 25  # 架线时考虑初伸长的降温，取正值。单位°C


# h_i 悬点高差
# l_i 悬点档距
# alpha 导线膨胀系数  1/°C
# E 弹性系数 N/mm2
# lambda_i 计算不平衡张力时导线比载 N/(m.mm)
# sigma_i 计算不平衡张力时最低点水平应力 单位N/mm2
# t_i 计算不平衡张力时导线温度 单位°C
# lambda_m 导线架线时时导线比载 N/(m.mm)
# sigma_m 导线架线时时最低点水平应力 单位N/mm2
# t_m 导线架线时时导线温度 单位°C
def delta_li(h_i, l_i, lambda_i, alpha, E, t_i, sigma_i, lambda_m, t_m, sigma_m):
    beta_i = math.atan(h_i / l_i)
    t = (
        l_i
        / ((math.cos(beta_i) ** 2) * (1 + (lambda_i * l_i / sigma_i) ** 2 / 8))
        * (
            (l_i * math.cos(beta_i)) ** 2
            / 24
            * ((lambda_m / sigma_m) ** 2 - (lambda_i / sigma_i) ** 2)
            + ((sigma_i - sigma_m) / E / math.cos(beta_i))
            + alpha * (t_i + t_e - t_m)
        )
    )


print("Finished.")