如何使用 Pyomo 重新制定和解决双层优化问题？

我想我已经使用 pyomo.bilevel 正确构建了我的模型和子模型。但是，我一直在尝试实际求解模型。当我在代码末尾使用以下几行时，它解决了“模型”也就是上层问题，同时完全忽略了下层问题：

opt = pyo.solverFactory('scip')

opt.solve(model)

我已经看到我需要使用类似的东西：

xfrm = TransformationFactory('gdp.bigm')

xfrm.apply_to(model)

但我不知道正确的实现方法，而且我找不到解决问题的代码行的任何示例双层问题。我想避免使用命令行来求解模型，而是希望在我的代码末尾有一些东西。非常感谢提前！