在Pyomo约束中定义循环/周期性边界条件

2024-05-17 • 问答

我正在尝试对利用循环条件的Pyomo模型定义约束。下面是我认为它应该工作的方式（GAMS的循环语法）。

from __future__ import division
from pyomo.environ import *

model = ConcreteModel()

## define sets
model.t                 = Set(initialize = [i for i in range(8760)]) 

## define variables
model.ESS_SOC           = Var(model.t,domain = NonNegativeReals) # battery state of charge
model.ESS_c             = Var(model.t,domain = NonNegativeReals) # battery charging
model.ESS_d             = Var(model.t,domain = NonNegativeReals) # battery discharging

## skip obj for this example

## define constraints
#SOC constraint
model.SOC_const = ConstraintList()
for i in model.t: 
    model.SOC_const.add( model.ESS_SOC[i] == model.ESS_SOC[i--1] + model.ESS_c[i] - model.ESS_d[i] )

但是当我运行上面的示例时，我收到以下错误消息：

KeyError: "Index '8760' is not valid for indexed component 'ESS_SOC'"

在给出model.t定义的情况下，我同意错误，但是错误让我相信它几乎正在执行我想要的操作，即：

model.ESS_SOC[0] == model.ESS_SOC[8759] + model.ESS_c[0] - model.ESS_d[0]
model.ESS_SOC[1] == model.ESS_SOC[0] + model.ESS_c[1] - model.ESS_d[1]
...
model.ESS_SOC[8759] == model.ESS_SOC[8758] + model.ESS_c[8759] - model.ESS_d[8759]

有没有一种方法来定义约束，这就是我所得到的？

from __future__ import division from pyomo.environ import * model = ConcreteModel() ## define sets model.t = Set(initialize = [i for i in range(8760)],ordered=True) ## define variables model.ESS_SOC = Var(model.t,domain = NonNegativeReals) # battery state of charge model.ESS_c = Var(model.t,domain = NonNegativeReals) # battery charging model.ESS_d = Var(model.t,domain = NonNegativeReals) # battery discharging ## skip obj for this example ## define constraints #SOC constraint def _SOC_const(m,i): if i == m.t.first(): return model.ESS_SOC[i] == model.ESS_SOC[m.t.last()] + model.ESS_c[i] - model.ESS_d[i] return model.ESS_SOC[i] == model.ESS_SOC[i-1] + model.ESS_c[i] - model.ESS_d[i] model.SOC_const = Constraint(model.t,rule=_SOC_const)

在Pyomo约束中定义循环/周期性边界条件

ildar_kh 回答：在Pyomo约束中定义循环/周期性边界条件

大家都在问