Pyomo报错：未初始化NumericValue对象问题求助

我是Pyomo新手，找不到难度适配的示例，现在遇到变量初始化问题无法解决。前两个约束运行正常，但执行第三个约束PM_rule时出现错误。

原始代码 #

# Import of the pyomo module
from pyomo.environ import *
 
# Creation of a Concrete Model
model = ConcreteModel()

## Define sets ##
#  Sets

model.I = RangeSet(1, N, doc='Farms')
model.J = RangeSet(1, n_max, doc='Turbines')
model.T = RangeSet(1, T, doc='Time horizon')
model.T_minus_1 = RangeSet(1, T-1, doc='Time horizon up to T-1')
model.R = RangeSet(1, R, doc='Maintenance Crews')

## Define variables ##
#  Variables

model.Xpm_ijt = Var(model.I, model.J, model.T, domain=Binary, doc='Binary variable for shipment decision')
model.Xcm_ijt = Var(model.I, model.J, model.T, domain=Binary, doc='Binary variable for shipment decision')

model.y_rt = Var(model.R, model.T, domain=Binary, doc='Binary variable for shipment decision')
model.z_rit = Var(model.R, model.I, model.T, domain=Binary, doc='Binary variable for shipment decision')
model.f_ijt = Var(model.I, model.J, model.T, domain=Binary, doc='Binary variable for shipment decision')

model.a_ijt = Var(model.I, model.J, model.T,  domain=NonNegativeIntegers, bounds=(0, H*10), doc='Operating time of turbine j at farm i at time t')

#Constraints 
#Constraint 1

def init_a_ijt_rule(model, i, j, t):
    if t == 1 and i == 1:
        if j == 1:
            return 0
        elif j == 2:
            return 0
        elif j == 3:
            return 10
        elif j == 4:
            return 20
        elif j == 5:
            return 23
    elif t == 1 and i == 2:
        if j == 1:
            return 12
        elif j == 2:
            return 20
        elif j == 3:
            return 34
    else:
        return 0

model.constraint_a_ijt = Var(model.I, model.J, model.T, initialize=init_a_ijt_rule, doc='Operating time of turbine j at farm i at time t')

#Constraint 2

def operating_time_rule(model, i, j, t):
    if model.Xpm_ijt[i, j, t].value == 1 or model.Xcm_ijt[i, j, t].value == 1:
        return model.a_ijt[i, j, t+1] == 0
    else:
        return model.a_ijt[i, j, t+1] == model.a_ijt[i, j, t] + 1

model.constraint_operating_time = Constraint(model.I, model.J, model.T_minus_1, rule=operating_time_rule)

#Constraint 3

def PM_rule(model, i, j, t):
    if value(model.a_ijt[i, j, t]) <= (H - h-1) or value(model.a_ijt[i, j, t]) >= H:
        return model.Xpm_ijt[i, j, t] == 0
    else:
        return Constraint.Skip

model.constraint_Xpm = Constraint(model.I, model.J, model.T, rule=PM_rule)

错误信息 #

ERROR: evaluating object as numeric value: a_ijt[1,1,1] (object: <class 'pyomo.core.base.var._GeneralVarData'>) No value for uninitialized NumericValue object a_ijt[1,1,1]
ERROR: Rule failed when generating expression for Constraint constraint_Xpm with index (1, 1, 1): ValueError: No value for uninitialized NumericValue object a_ijt[1,1,1]
ERROR: Constructing component 'constraint_Xpm' from data=None failed: ValueError: No value for uninitialized NumericValue object a_ijt[1,1,1]

问题分析与修复 #

1. 变量初始化错误

你错误地创建了新变量model.constraint_a_ijt来存储初始化值，而没有给目标变量model.a_ijt赋值。正确的做法是直接在定义model.a_ijt时指定initialize参数：

# 修正model.a_ijt的定义，直接添加initialize参数
model.a_ijt = Var(model.I, model.J, model.T,  
                  domain=NonNegativeIntegers, 
                  bounds=(0, H*10), 
                  initialize=init_a_ijt_rule,  # 用初始化规则直接初始化变量
                  doc='Operating time of turbine j at farm i at time t')

# 删除错误的model.constraint_a_ijt定义

2. 约束规则中直接访问变量值的错误

Pyomo模型构建阶段，变量还未被求解器赋值，不能通过.value或value()函数获取变量值编写条件判断。约束必须是基于变量的代数表达式，可通过大M法或Indicator约束实现逻辑条件：

修正约束2（operating_time_rule）

使用大M法将逻辑条件转化为线性约束（M取大于变量最大可能值的数，比如H*10 +1）：

def operating_time_rule(model, i, j, t):
    M = H*10 + 1  # 足够大的常数，避免约束失效
    # 当Xpm或Xcm为1时，a_ijt[t+1]必须为0；否则a_ijt[t+1] = a_ijt[t]+1
    expr1 = model.a_ijt[i,j,t+1] <= M * (1 - model.Xpm_ijt[i,j,t])
    expr2 = model.a_ijt[i,j,t+1] <= M * (1 - model.Xcm_ijt[i,j,t])
    expr3 = model.a_ijt[i,j,t+1] >= (model.a_ijt[i,j,t] + 1) - M * (model.Xpm_ijt[i,j,t] + model.Xcm_ijt[i,j,t])
    expr4 = model.a_ijt[i,j,t+1] <= (model.a_ijt[i,j,t] + 1) + M * (model.Xpm_ijt[i,j,t] + model.Xcm_ijt[i,j,t])
    return expr1 & expr2 & expr3 & expr4

model.constraint_operating_time = Constraint(model.I, model.J, model.T_minus_1, rule=operating_time_rule)

修正约束3（PM_rule）

将基于变量值的条件转化为线性约束，确保当a_ijt[i,j,t]不在(H-h-1, H)区间时，Xpm_ijt[i,j,t]必须为0：

def PM_rule(model, i, j, t):
    L = H - h - 1
    # 当a_ijt <= L时，Xpm必须为0；当a_ijt >= H时，Xpm必须为0
    expr1 = model.Xpm_ijt[i,j,t] <= model.a_ijt[i,j,t] - L  # a<=L时右边<=0，Xpm只能为0
    expr2 = model.Xpm_ijt[i,j,t] <= H - model.a_ijt[i,j,t]  # a>=H时右边<=0，Xpm只能为0
    return expr1 & expr2

model.constraint_Xpm = Constraint(model.I, model.J, model.T, rule=PM_rule)

额外注意事项 #

• 确保代码中N、n_max、T、R、H、h这些参数已提前定义，否则模型无法构建。
• 大M法中M的取值要合理，既不能太大导致数值不稳定，也不能太小导致约束失效。

内容的提问来源于stack exchange，提问作者Khadija B