基于受力计算3D悬索长度的代码偏差问题求助

我编写了build_3D_line函数，用于基于受力参数计算3D空间中悬索的长度，输入参数包括：

• 水平恒力H
• 单位长度垂直载荷q_verti（如重力）
• 单位长度横向载荷q_trans（如风载）
• 支座水平间距displacement
• 支座高度差dheight
• 计算精度accuracy

测试参数 #

H = 8969.779999999999
q_verti = 8.29
q_trans = 0
displacement = 233.40000000002328
dheight = 1.379999999999999
accuracy = 10**5

问题现象 #

函数返回的悬索长度为233.84650993989436，但预期结果应为233.85685400689795。通过返回结果反推的水平力与给定的H值不符，反推公式参考F. Kiessling等人著作（ISBN 3-540-00297-9）：

horizontal_force = (L2EA) / (L1*(1 + eps_t*(T2 - T1))) + H1 - E*A

其中用到的材料参数：

• E=127530000000（杨氏模量）
• A=9.327e-05（截面积）
• eps_t=1.7e-5（热伸长系数）
• T1=T2=10（温度）

求助需求 #

我无法定位代码中的错误，希望有人能指出build_3D_line函数的问题所在，或提供基于给定输入计算3D悬索长度的其他可靠方法。

内容的提问来源于stack exchange，提问作者Tim Moore