Pyomo定义约束触发无效表达式返回布尔False报错咨询

报错原因 #

这个错误的核心诱因是生成器表达式的条件筛选位置书写错误，导致约束表达式在Python原生层面被提前计算为普通布尔值，没有生成Pyomo可识别的符号约束对象：

• 原代码里的if i!=j被写在了外层遍历r的生成器末尾，但j是内层遍历索引集modelo.j的循环变量，外层生成器执行时无法访问到内层的j，会直接读取当前运行环境外层作用域已定义的j变量。
• 当外层作用域的j取值恰好等于当前迭代的i时，外层生成器的筛选条件i!=j恒为False，外层sum无法遍历到任何r，两层嵌套求和的结果直接是Python原生整数0，此时0 == 1会直接返回原生布尔值False，不符合Pyomo约束的返回要求，就会触发该报错。

修正方法 #

把if i!=j的筛选条件移动到内层遍历j的生成器末尾，让条件正确作用于j的筛选，保证整个求和过程生成Pyomo符号表达式，不会提前计算为原生数值：

def rest1(modelo, i):
    return sum(sum(modelo.M[r, i, j] for j in modelo.j if i != j) for r in modelo.r) == 1

如果要提升代码鲁棒性，避免极端场景（比如某个i对应的j集合只有i自身，导致求和结果恒为0）下再次触发同类报错，可以增加常数判断逻辑，按照Pyomo规范返回不可行标记或跳过冗余约束：

def rest1(modelo, i):
    total = sum(sum(modelo.M[r, i, j] for j in modelo.j if i != j) for r in modelo.r)
    if total.is_constant():
        # 求和结果为固定常数，不满足等于1的要求时标记为不可行，满足则跳过该约束
        return pyEnv.Constraint.Infeasible if total() != 1 else pyEnv.Constraint.Skip
    return total == 1

修改后再执行原有的约束注册代码即可正常生成约束。

内容的提问来源于stack exchange，提问作者abner costa