commit f38c73c001faf5c2416d0f3e6f8c9b330e20026a
Author: facat <facat@facat.com>
Date:   Fri Jun 26 22:12:16 2020 +0800

    第一次提交。

diff --git a/.gitignore b/.gitignore
new file mode 100644
index 0000000..f5e96db
--- /dev/null
+++ b/.gitignore
@@ -0,0 +1 @@
+venv
\ No newline at end of file
diff --git a/main.py b/main.py
new file mode 100644
index 0000000..fed635b
--- /dev/null
+++ b/main.py
@@ -0,0 +1,38 @@
+# 计算直线塔不平衡张力
+# 新版输电线路大手册 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.")