使用sympy.solve触发IndexError: list index out of range错误求助

我尝试推导关于B的表达式，但运行代码时触发错误：IndexError: list index out of range。除最后一行代码外，其余部分均运行正常。此前一直使用该方法，本次却无法成功。相关代码如下：

#Tut 10 Question 3 Attempt Almost got there...

import sympy as sym  # 原代码缺少导入语句，需补充

x = sym.Symbol("x")
L,a,P = sym.Symbol("L"),sym.Symbol("a"),sym.Symbol("P")
A,B,C,D,E,F = sym.Symbol("A"),sym.Symbol("B"),sym.Symbol("C"),sym.Symbol("D"),sym.Symbol("E"),sym.Symbol("F")

#LHS
print("LHS")
lV = (L-a)*(P) #Shear force remains constant throughout 0 < x < a

lM = sym.integrate(lV,x) + A
lM = lM.subs(A,0) 
display(lM)

ltheta = sym.integrate(lM,x) + B
display(ltheta)

lv = sym.integrate(ltheta,x) + C
lv = lv.subs(C,0)
display(lv)

#RHS
print("RHS")
rV = (a)*(P) #Shear force remains constant throughout 0 < x < a  # 此处注释疑似错误，应为a < x < L？

rM = sym.integrate(rV,x) + D
rM = rM.subs(D,0) 
display(rM)

rtheta = sym.integrate(rM,x) + E
display(rtheta)

rv = sym.integrate(rtheta,x) + F
rv = rv.subs(F,0)
display(rv)

print("combined")
CombinedTheta = rtheta - ltheta
CombinedTheta = CombinedTheta.subs(x,a)   ####  Look at the Hints of the question thethaLHS = thetaRHS at x =a
display(CombinedTheta)

Combinedv = rv - lv
Combinedv = Combinedv.subs(x,a) ####  Look at the Hints of the question vLHS = vRHS at x =a
display(Combinedv)

#now we will solve E i.t.o B
print("solving for E")
factorE = sym.solve(CombinedTheta,E)[0]
display(factorE)

#now we will solve B by plugging factorE into E of Combinedv
print("solving for B")
Combinedv = Combinedv.subs(E,factorE)
Equation = Combinedv
display(Equation)
FactorB = sym.solve(Equation,B)[0]

错误原因分析 #

IndexError的直接原因是sym.solve(Equation,B)返回了空列表，说明当前方程无法解出关于B的有效表达式——要么方程是恒等式（B被完全消去），要么是矛盾式，或者推导过程中存在错误。

排查方向 #

1. 检查方程内容：在调用solve前打印Equation，确认是否包含变量B：

   print(Equation)
   

   如果输出中没有B，说明前面的推导中B已被消去，需回溯检查。
2. 剪力定义错误：代码中RHS的剪力注释范围明显矛盾，且剪力的正负号、取值是否符合题目中的梁受力情况？错误的剪力会导致后续积分结果全部出错。
3. 边界条件误用：提前将积分常数A、C、D、F设为0是否合理？部分边界条件需要联立求解，而非直接代入0，这可能导致方程失去求解B的条件。
4. 积分过程验证：逐段检查弯矩、转角、位移的积分结果是否正确，比如LHS的弯矩积分是否应为(L-a)*P*x + A，代入A=0后是否符合预期。

内容的提问来源于stack exchange，提问作者Benro Snyman