参数和计算分离。

data.py

@@ -0,0 +1,18 @@
+loop_end = 100000  # 最大循环次数
+# 架线时的状态
+# 取外过无风
+string_length = 9.2  # 串长 单位m
+string_g = 60 * 9.8  # 串重 单位N
+t_m = 15  # 导线架设时的气温。单位°C
+t_e = 20  # 架线时考虑初伸长的降温，取正值。单位°C
+alpha = 0.0000155  # 导线膨胀系数  1/°C
+elastic = 95900  # 弹性系数 N/mm2
+area = 154.48  # 导线面积 mm2
+lambda_m = 14.8129 / area  # 导线比载 N/(m.mm)
+# 取400m代表档距下
+sigma_m = 28517 / area  # 架线时，初伸长未释放前的最低点水平应力。单位N/mm2
+span_count = 3  # 几个档距
+# n个档距,n-1个直线塔
+h_array = [0, 0, 0]
+l_array = [400, 400, 400]
+t_array = [15, 15, 15]

main.py

@@ -2,6 +2,7 @@
 # 新版输电线路大手册 P328
 
 import math
+import data
 
 # h_i 悬点高差
 # l_i 悬点档距
@@ -82,24 +83,23 @@ def fun_sigma_i1(
 
 # 求解循环。
 def cal_loop():
-    loop_end = 100000  # 最大循环次数
+    loop_end = data.loop_end  # 最大循环次数
     # 架线时的状态
     # 取外过无风
-    string_length = 9.2  # 串长 单位m
-    string_g = 60 * 9.8  # 串重 单位N
-    t_m = 15  # 导线架设时的气温。单位°C
-    t_e = 20  # 架线时考虑初伸长的降温，取正值。单位°C
-    alpha = 0.0000155  # 导线膨胀系数  1/°C
-    elastic = 95900  # 弹性系数 N/mm2
-    area = 154.48  # 导线面积 mm2
-    lambda_m = 14.8129 / area  # 导线比载 N/(m.mm)
-    # 取400m代表档距下
-    sigma_m = 28517 / area  # 架线时，初伸长未释放前的最低点水平应力。单位N/mm2
-    span_count = 3  # 几个档距
+    string_length = data.string_length  # 串长 单位m
+    string_g = data.string_g  # 串重 单位N
+    t_m = data.t_m  # 导线架设时的气温。单位°C
+    t_e = data.t_e  # 架线时考虑初伸长的降温，取正值。单位°C
+    alpha = data.alpha  # 导线膨胀系数  1/°C
+    elastic = data.elastic  # 弹性系数 N/mm2
+    area = data.area  # 导线面积 mm2
+    lambda_m = data.lambda_m   # 导线比载 N/(m.mm)
+    sigma_m = data.sigma_m   # 架线时，初伸长未释放前的最低点水平应力。单位N/mm2
+    span_count = data.span_count  # 几个档距
     # n个档距,n-1个直线塔
-    h_array = [0, 0, 0]
-    l_array = [400, 400, 400]
-    t_array = [15, 15, 15]
+    h_array = data.h_array
+    l_array = data.l_array
+    t_array = data.t_array
     lambda_array = [lambda_m, lambda_m, lambda_m]
     loop_count = 1
     sigma_0 = sigma_m * 0.8