Pyomo约束报错：布尔值True替代Pyomo对象的问题求助

问题背景 #

在基于Pyomo实现的课程规划模型中，读取Excel数据后生成约束max_kurs_pro_lehrer_in_periode时触发报错：

ERROR: Rule failed when generating expression for Constraint max_kurs_pro_lehrer_in_periode with index ('AMU', 'CW40-Mon'): ValueError: Invalid constraint expression. The constraint expression resolved to a trivial Boolean (True) instead of a Pyomo object. Please modify your rule to return Constraint.Feasible instead of True.

报错核心原因：当某教师可教授的课程类型在规划周期内无需求时，约束中的求和项迭代器为空，sum()返回0，此时表达式0 <= lcap[l,t]会直接计算为布尔值True（无论lcap[l,t]是1还是0），而非Pyomo要求的表达式对象，从而触发错误。

解决方案 #

方法1：在约束规则内判断求和项是否为空 #

修改约束规则函数，先检查求和项是否存在，若为空则返回Constraint.Feasible（表示该约束恒成立，无需添加），否则返回正常的约束表达式：

from pyomo.environ import Constraint

def max_kurs_pro_lehrer_in_periode_rule(model, l, t):
    # 构建求和项列表，用于判断是否为空
    sum_terms = [model.gamma[l, k, s, t] 
                 for m in ML[l] if m not in D_gesamt 
                 for k in KM[m] 
                 for s in SM[m]]
    if not sum_terms:
        # 无课程可分配，约束恒成立，返回可行标识
        return Constraint.Feasible
    # 有课程项时返回正常约束
    return sum(sum_terms) <= lcap[l, t]

model.max_kurs_pro_lehrer_in_periode = Constraint(L, T, rule=max_kurs_pro_lehrer_in_periode_rule)

方法2：提前筛选有效索引对 #

先筛选出存在可分配课程的教师-时段对，仅对这些索引生成约束，从根源避免空求和的情况：

# 筛选存在可分配课程的教师-时段对
valid_lt_pairs = []
for l in L:
    for t in T:
        sum_terms = [model.gamma[l, k, s, t] 
                     for m in ML[l] if m not in D_gesamt 
                     for k in KM[m] 
                     for s in SM[m]]
        if sum_terms:
            valid_lt_pairs.append((l, t))

# 仅对有效索引对生成约束
model.max_kurs_pro_lehrer_in_periode = Constraint(valid_lt_pairs, rule=max_kurs_pro_lehrer_in_periode_rule)

原理说明 #

Pyomo要求约束规则的返回值必须是Pyomo表达式对象（如不等式、等式），而当求和项为空时，0 <= lcap[l,t]会被Python直接计算为布尔值，而非Pyomo表达式。使用Constraint.Feasible可以明确告诉Pyomo该约束无需添加，因为它恒成立，从而避免报错。

内容的提问来源于stack exchange，提问作者Zeineb Kyrkianidou-Toumi