У меня есть следующее ограничение:
Код: Выделить всё
def grid_power_calculation_with_balance(model, t):
if t == 0:
return model.P_grid[t] == 0
return model.P_grid[t] == model.PCh[t] - model.PDch[t] + \
sum(model.gamma[i, t] * (cooling_capacity(value(model.T[i, t - 1])) / COP(value(model.T[i, t - 1])) * model.delta[i, t] + P_aux[i]) for i in I) - model.P_used_renew[t]
model.grid_power_calculation = Constraint(T,rule=grid_power_calculation_with_balance)
Код: Выделить всё
model.delta = Var(I, T, within=Binary) # Decision variable: 1 if cooling, 0 if heating
model.P_used_renew = Var(T, within=NonNegativeReals)
def renew_usage_limit(model, t):
return model.P_used_renew[t]
Подробнее здесь: [url]https://stackoverflow.com/questions/79309902/pyomo-not-respecting-constraint[/url]
Мобильная версия