本文整理汇总了Python中pyomo.core.ConcreteModel.P方法的典型用法代码示例。如果您正苦于以下问题：Python ConcreteModel.P方法的具体用法？Python ConcreteModel.P怎么用？Python ConcreteModel.P使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类pyomo.core.ConcreteModel的用法示例。
在下文中一共展示了ConcreteModel.P方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: build_model
# 需要导入模块: from pyomo.core import ConcreteModel [as 别名]
# 或者: from pyomo.core.ConcreteModel import P [as 别名]
def build_model(use_mccormick=False):
"""Build the GDP model."""
m = ConcreteModel()
m.F = Var(bounds=(0, 8), doc="Flow into reactor")
m.X = Var(bounds=(0, 1), doc="Reactor conversion")
m.d = Param(initialize=2, doc="Max product demand")
m.c = Param([1, 2, 'I', 'II'], doc="Costs", initialize={
1: 2, # Value of product
2: 0.2, # Cost of raw material
'I': 2.5, # Cost of reactor I
'II': 1.5 # Cost of reactor II
})
m.alpha = Param(['I', 'II'], doc="Reactor coefficient",
initialize={'I': -8, 'II': -10})
m.beta = Param(['I', 'II'], doc="Reactor coefficient",
initialize={'I': 9, 'II': 15})
m.X_LB = Param(['I', 'II'], doc="Reactor conversion lower bound",
initialize={'I': 0.2, 'II': 0.7})
m.X_UB = Param(['I', 'II'], doc="Reactor conversion upper bound",
initialize={'I': 0.95, 'II': 0.99})
m.C_rxn = Var(bounds=(1.5, 2.5), doc="Cost of reactor")
m.reactor_choice = Disjunction(expr=[
# Disjunct 1: Choose reactor I
[m.F == m.alpha['I'] * m.X + m.beta['I'],
m.X_LB['I'] <= m.X,
m.X <= m.X_UB['I'],
m.C_rxn == m.c['I']],
# Disjunct 2: Choose reactor II
[m.F == m.alpha['II'] * m.X + m.beta['II'],
m.X_LB['II'] <= m.X,
m.X <= m.X_UB['II'],
m.C_rxn == m.c['II']]
], xor=True)
if use_mccormick:
m.P = Var(bounds=(0, 8), doc="McCormick approximation of F*X")
m.mccormick_1 = Constraint(
expr=m.P <= m.F.lb * m.X + m.F * m.X.ub - m.F.lb * m.X.ub,
doc="McCormick overestimator")
m.mccormick_2 = Constraint(
expr=m.P <= m.F.ub * m.X + m.F * m.X.lb - m.F.ub * m.X.lb,
doc="McCormick underestimator")
m.max_demand = Constraint(expr=m.P <= m.d, doc="product demand")
m.profit = Objective(
expr=m.c[1] * m.P - m.c[2] * m.F - m.C_rxn, sense=maximize)
else:
m.max_demand = Constraint(expr=m.F * m.X <= m.d, doc="product demand")
m.profit = Objective(
expr=m.c[1] * m.F * m.X - m.c[2] * m.F - m.C_rxn, sense=maximize)
return m