# MINLP written by GAMS Convert at 02/17/22 17:19:49 # # Equation counts # Total E G L N X C B # 62 36 12 14 0 0 0 0 # # Variable counts # x b i s1s s2s sc si # Total cont binary integer sos1 sos2 scont sint # 57 39 18 0 0 0 0 0 # FX 0 # # Nonzero counts # Total const NL # 212 58 154 # # Reformulation has removed 1 variable and 1 equation from pyomo.environ import * model = m = ConcreteModel() m.b1 = Var(within=Binary, bounds=(0,1), initialize=0) m.b2 = Var(within=Binary, bounds=(0,1), initialize=1) m.b3 = Var(within=Binary, bounds=(0,1), initialize=0) m.b4 = Var(within=Binary, bounds=(0,1), initialize=0) m.b5 = Var(within=Binary, bounds=(0,1), initialize=1) m.b6 = Var(within=Binary, bounds=(0,1), initialize=0) m.b7 = Var(within=Binary, bounds=(0,1), initialize=1) m.b8 = Var(within=Binary, bounds=(0,1), initialize=0) m.b9 = Var(within=Binary, bounds=(0,1), initialize=0) m.b10 = Var(within=Binary, bounds=(0,1), initialize=0) m.b11 = Var(within=Binary, bounds=(0,1), initialize=1) m.b12 = Var(within=Binary, bounds=(0,1), initialize=0) m.b13 = Var(within=Binary, bounds=(0,1), initialize=1) m.b14 = Var(within=Binary, bounds=(0,1), initialize=0) m.b15 = Var(within=Binary, bounds=(0,1), initialize=0) m.b16 = Var(within=Binary, bounds=(0,1), initialize=1) m.b17 = Var(within=Binary, bounds=(0,1), initialize=0) m.b18 = Var(within=Binary, bounds=(0,1), initialize=0) m.x19 = Var(within=Reals, bounds=(0.526315789473684,1.05263157894737) , initialize=1) m.x20 = Var(within=Reals, bounds=(0.961538461538462,2.11538461538462) , initialize=1) m.x21 = Var(within=Reals, bounds=(0.2,1), initialize=0.3) m.x22 = Var(within=Reals, bounds=(0,0.8), initialize=0.8) m.x23 = Var(within=Reals, bounds=(6,13), initialize=13) m.x24 = Var(within=Reals, bounds=(6,13), initialize=10.6604143883993) m.x25 = Var(within=Reals, bounds=(0.26,0.35), initialize=0.3273636888081) m.x26 = Var(within=Reals, bounds=(0,90), initialize=27.4760701692312) m.x27 = Var(within=Reals, bounds=(4.9,5.5), initialize=5.30016785671899) m.x28 = Var(within=Reals, bounds=(55,63), initialize=58.4294759668997) m.x29 = Var(within=Reals, bounds=(0,None), initialize=0) m.x30 = Var(within=Reals, bounds=(0,200), initialize=82.6048888012009) m.x31 = Var(within=Reals, bounds=(0.296392099803303,0.404171045186323) , initialize=0.332831162843986) m.x32 = Var(within=Reals, bounds=(0.134723681728774,0.229030258938916) , initialize=0.159266373947579) m.x33 = Var(within=Reals, bounds=(0,26), initialize=6.62084546202556) m.x34 = Var(within=Reals, bounds=(0,1), initialize=0.556542937336767) m.x35 = Var(within=Reals, bounds=(0,34.1), initialize=23.6919439882642) m.x36 = Var(within=Reals, bounds=(0,1), initialize=0.913247121927004) m.x37 = Var(within=Reals, bounds=(0.526315789473684,1.05263157894737) , initialize=0.8) m.x38 = Var(within=Reals, bounds=(0.961538461538462,2.11538461538462) , initialize=0.961538461538462) m.x39 = Var(within=Reals, bounds=(0.2,1), initialize=0.4) m.x40 = Var(within=Reals, bounds=(0,0.8), initialize=0.8) m.x41 = Var(within=Reals, bounds=(6,13), initialize=13) m.x42 = Var(within=Reals, bounds=(6,13), initialize=12.0499480711157) m.x43 = Var(within=Reals, bounds=(0.26,0.35), initialize=0.31546674627608) m.x44 = Var(within=Reals, bounds=(0,90), initialize=36.408786197323) m.x45 = Var(within=Reals, bounds=(4.9,5.5), initialize=5.26747485493678) m.x46 = Var(within=Reals, bounds=(55,63), initialize=59.9333303564984) m.x47 = Var(within=Reals, bounds=(0,None), initialize=0) m.x48 = Var(within=Reals, bounds=(0,200), initialize=58.2740173254705) m.x49 = Var(within=Reals, bounds=(0.296392099803303,0.404171045186323) , initialize=0.334395741702399) m.x50 = Var(within=Reals, bounds=(0.134723681728774,0.229030258938916) , initialize=0.165538447790303) m.x51 = Var(within=Reals, bounds=(0,26), initialize=9.8912291834487) m.x52 = Var(within=Reals, bounds=(0,1), initialize=0.659560560750654) m.x53 = Var(within=Reals, bounds=(0,34.1), initialize=17.380013231931) m.x54 = Var(within=Reals, bounds=(0,1), initialize=0.821759148100612) m.x55 = Var(within=Reals, bounds=(0,200), initialize=30) m.x56 = Var(within=Reals, bounds=(0,34.1), initialize=7.53758887227811) m.x57 = Var(within=Reals, bounds=(0,1), initialize=0.588460387570554) m.obj = Objective(sense=minimize, expr= m.x29 + m.x47) m.e1 = Constraint(expr= -21.6 * m.x21 * m.x22 * m.x27 + m.x26 == 0) m.e2 = Constraint(expr= -0.007852585706 * m.x19**3 - 0.154288922601 * m.x20**3 - 0.352933730854 * m.x21**3 + 0.004816150342 * m.x19**2 * m.x20 + 0.00547943134 * m.x20**2 * m.x19 + 0.02533808214 * m.x19**2 * m.x21 - 0.00021201136 * m.x21**2 * m.x19 + 0.057118497613 * m.x20**2 * m.x21 + 0.042739509965 * m.x21**2 * m.x20 + 0.01583097252 * m.x19 * m.x20 * m.x21 + 0.001028174658 * m.x19**2 + 0.805369774847 * m.x20**2 + 0.655580550098 * m.x21**2 - 0.057270405947 * m.x19 * m.x20 - 0.07973036236 * m.x19 * m.x21 - 0.342091579946 * m.x20 * m.x21 + 0.191345333621 * m.x19 - 1.188971392024 * m.x20 + 0.346682012779 * m.x21 + m.x27 == 4.960068215723) m.e3 = Constraint(expr= -2.21406746341 * m.x19**3 - 1.086659693552 * m.x20**3 - 5.577874978662 * m.x21**3 + 0.815241697738 * m.x19**2 * m.x20 - 0.509578110533 * m.x20**2 * m.x19 - 1.561758113326 * m.x19**2 * m.x21 - 2.212321055022 * m.x21**2 * m.x19 + 0.612567680918 * m.x20**2 * m.x21 - 0.254008083604 * m.x21**2 * m.x20 + 0.159429747244 * m.x19 * m.x20 * m.x21 + 8.905599398536 * m.x19**2 + 6.095001164559 * m.x20**2 + 15.207539664993 * m.x21**2 - 0.089172114876 * m.x19 * m.x20 + 3.273526677614 * m.x19 * m.x21 - 2.498376358946 * m.x20 * m.x21 - 2.621894664006 * m.x19 - 9.284846067558 * m.x20 - 5.837143728557 * m.x21 + m.x28 == 57.679680208231) m.e4 = Constraint(expr= -1.456640469666 * m.x19**3 + 5.495718264905 * m.x20**3 + 28.456261951645 * m.x21**3 - 0.912917970314 * m.x19**2 * m.x20 + 0.88119920631 * m.x20**2 * m.x19 + 1.049763024383 * m.x19**2 * m.x21 + 0.308107344863 * m.x21**2 * m.x19 - 2.043536297441 * m.x20**2 * m.x21 - 15.609611231641 * m.x21**2 * m.x20 - 0.336486837518 * m.x19 * m.x20 * m.x21 + 4.634160849469 * m.x19**2 - 31.478262635483 * m.x20**2 - 34.016843490037 * m.x21**2 - 1.153148892739 * m.x19 * m.x20 - 1.168601192983 * m.x19 * m.x21 + 32.056936006397 * m.x20 * m.x21 - 3.405095041238 * m.x19 + 54.472915571467 * m.x20 - 9.56987912824 * m.x21 + m.x23 == 44.230616625681) m.e5 = Constraint(expr= -3.334445194766 * m.x19**3 + 2.265666208775 * m.x20**3 + 20.256566414583 * m.x21**3 - 0.413782262402 * m.x19**2 * m.x20 + 3.523622273943 * m.x20**2 * m.x19 + 0.285910055687 * m.x19**2 * m.x21 + 10.110726634622 * m.x21**2 * m.x19 - 1.95072196814 * m.x20**2 * m.x21 - 10.308512463418 * m.x21**2 * m.x20 - 5.808426325827 * m.x19 * m.x20 * m.x21 + 6.932398033967 * m.x19**2 - 15.80019426934 * m.x20**2 - 39.197963873266 * m.x21**2 - 7.900706395772 * m.x19 * m.x20 - 6.58186092156 * m.x19 * m.x21 + 30.119438887106 * m.x20 * m.x21 + 6.733798415788 * m.x19 + 26.385308892431 * m.x20 + 4.098268423019 * m.x21 + m.x24 == 32.102172356117) m.e6 = Constraint(expr= 0.194075741585 * m.x19**2 + 0.004843420334 * m.x20**2 - 0.04736686635 * m.x21**2 - 9.4029979e-05 * m.x19 * m.x20 - 0.011329651793 * m.x19 * m.x21 + 0.001017352942 * m.x20 * m.x21 - 0.382275988592 * m.x19 - 0.019484588535 * m.x20 + 0.077357069039 * m.x21 + m.x25 == 0.140278656706) m.e7 = Constraint(expr= m.x23 <= 11) m.e8 = Constraint(expr= m.x24 <= 11) m.e9 = Constraint(expr= m.x25 >= 0.32) m.e10 = Constraint(expr= -21.6 * m.x39 * m.x40 * m.x45 + m.x44 == 0) m.e11 = Constraint(expr= -0.007852585706 * m.x37**3 - 0.154288922601 * m.x38**3 - 0.352933730854 * m.x39**3 + 0.004816150342 * m.x37**2 * m.x38 + 0.00547943134 * m.x38**2 * m.x37 + 0.02533808214 * m.x37**2 * m.x39 - 0.00021201136 * m.x39**2 * m.x37 + 0.057118497613 * m.x38**2 * m.x39 + 0.042739509965 * m.x39**2 * m.x38 + 0.01583097252 * m.x37 * m.x38 * m.x39 + 0.001028174658 * m.x37**2 + 0.805369774847 * m.x38**2 + 0.655580550098 * m.x39**2 - 0.057270405947 * m.x37 * m.x38 - 0.07973036236 * m.x37 * m.x39 - 0.342091579946 * m.x38 * m.x39 + 0.191345333621 * m.x37 - 1.188971392024 * m.x38 + 0.346682012779 * m.x39 + m.x45 == 4.960068215723) m.e12 = Constraint(expr= -2.21406746341 * m.x37**3 - 1.086659693552 * m.x38**3 - 5.577874978662 * m.x39**3 + 0.815241697738 * m.x37**2 * m.x38 - 0.509578110533 * m.x38**2 * m.x37 - 1.561758113326 * m.x37**2 * m.x39 - 2.212321055022 * m.x39**2 * m.x37 + 0.612567680918 * m.x38**2 * m.x39 - 0.254008083604 * m.x39**2 * m.x38 + 0.159429747244 * m.x37 * m.x38 * m.x39 + 8.905599398536 * m.x37**2 + 6.095001164559 * m.x38**2 + 15.207539664993 * m.x39**2 - 0.089172114876 * m.x37 * m.x38 + 3.273526677614 * m.x37 * m.x39 - 2.498376358946 * m.x38 * m.x39 - 2.621894664006 * m.x37 - 9.284846067558 * m.x38 - 5.837143728557 * m.x39 + m.x46 == 57.679680208231) m.e13 = Constraint(expr= -1.456640469666 * m.x37**3 + 5.495718264905 * m.x38**3 + 28.456261951645 * m.x39**3 - 0.912917970314 * m.x37**2 * m.x38 + 0.88119920631 * m.x38**2 * m.x37 + 1.049763024383 * m.x37**2 * m.x39 + 0.308107344863 * m.x39**2 * m.x37 - 2.043536297441 * m.x38**2 * m.x39 - 15.609611231641 * m.x39**2 * m.x38 - 0.336486837518 * m.x37 * m.x38 * m.x39 + 4.634160849469 * m.x37**2 - 31.478262635483 * m.x38**2 - 34.016843490037 * m.x39**2 - 1.153148892739 * m.x37 * m.x38 - 1.168601192983 * m.x37 * m.x39 + 32.056936006397 * m.x38 * m.x39 - 3.405095041238 * m.x37 + 54.472915571467 * m.x38 - 9.56987912824 * m.x39 + m.x41 == 44.230616625681) m.e14 = Constraint(expr= -3.334445194766 * m.x37**3 + 2.265666208775 * m.x38**3 + 20.256566414583 * m.x39**3 - 0.413782262402 * m.x37**2 * m.x38 + 3.523622273943 * m.x38**2 * m.x37 + 0.285910055687 * m.x37**2 * m.x39 + 10.110726634622 * m.x39**2 * m.x37 - 1.95072196814 * m.x38**2 * m.x39 - 10.308512463418 * m.x39**2 * m.x38 - 5.808426325827 * m.x37 * m.x38 * m.x39 + 6.932398033967 * m.x37**2 - 15.80019426934 * m.x38**2 - 39.197963873266 * m.x39**2 - 7.900706395772 * m.x37 * m.x38 - 6.58186092156 * m.x37 * m.x39 + 30.119438887106 * m.x38 * m.x39 + 6.733798415788 * m.x37 + 26.385308892431 * m.x38 + 4.098268423019 * m.x39 + m.x42 == 32.102172356117) m.e15 = Constraint(expr= 0.194075741585 * m.x37**2 + 0.004843420334 * m.x38**2 - 0.04736686635 * m.x39**2 - 9.4029979e-05 * m.x37 * m.x38 - 0.011329651793 * m.x37 * m.x39 + 0.001017352942 * m.x38 * m.x39 - 0.382275988592 * m.x37 - 0.019484588535 * m.x38 + 0.077357069039 * m.x39 + m.x43 == 0.140278656706) m.e16 = Constraint(expr= m.x41 <= 11) m.e17 = Constraint(expr= m.x42 <= 11) m.e18 = Constraint(expr= m.x43 >= 0.32) m.e19 = Constraint(expr= exp(-0.029595 * m.x26) * (33.7894914681534 + m.x26) + m.x33 == 33.7894914681534) m.e20 = Constraint(expr= exp(-0.029595 * m.x26) + m.x34 == 1) m.e21 = Constraint(expr= -0.134723681728774 * (0.010073140669 * m.x19**2 + 0.011394190823 * m.x20**2 + 0.052910213683 * m.x21**2 + 0.000159410872 * m.x19 * m.x20 + 0.008036404292 * m.x19 * m.x21 - 0.003423392047 * m.x20 * m.x21 + 0.097124049148 * m.x19 + 0.03829180344 * m.x20 + 0.370440556286 * m.x21) + m.x31 == 0.29587368369345) m.e22 = Constraint(expr= -0.134723681728774 * (0.46598008632 * m.x19**2 - 0.00797004615 * m.x20**2 - 0.01779288613 * m.x21**2 - 0.01429434551 * m.x19 * m.x20 - 0.03832188467 * m.x19 * m.x21 + 0.00970510229 * m.x20 * m.x21 - 0.88981702163 * m.x19 + 0.07730602595 * m.x20 + 0.39988032723 * m.x21) + m.x32 == 0.194162178290626) m.e23 = Constraint(expr= -m.x26 * (2715.7894736842 / m.x27 + 5187 * m.x31 - 5187 * m.x32) / (4320 * m.x21 - 5187 * m.x32) + m.x30 == 0) m.e24 = Constraint(expr= exp(-0.029595 * m.x30) * (33.7894914681534 + m.x30) + m.x35 == 33.7894914681534) m.e25 = Constraint(expr= exp(-0.029595 * m.x30) + m.x36 == 1) m.e26 = Constraint(expr= exp(-0.029595 * m.x44) * (33.7894914681534 + m.x44) + m.x51 == 33.7894914681534) m.e27 = Constraint(expr= exp(-0.029595 * m.x44) + m.x52 == 1) m.e28 = Constraint(expr= -0.134723681728774 * (0.010073140669 * m.x37**2 + 0.011394190823 * m.x38**2 + 0.052910213683 * m.x39**2 + 0.000159410872 * m.x37 * m.x38 + 0.008036404292 * m.x37 * m.x39 - 0.003423392047 * m.x38 * m.x39 + 0.097124049148 * m.x37 + 0.03829180344 * m.x38 + 0.370440556286 * m.x39) + m.x49 == 0.29587368369345) m.e29 = Constraint(expr= -0.134723681728774 * (0.46598008632 * m.x37**2 - 0.00797004615 * m.x38**2 - 0.01779288613 * m.x39**2 - 0.01429434551 * m.x37 * m.x38 - 0.03832188467 * m.x37 * m.x39 + 0.00970510229 * m.x38 * m.x39 - 0.88981702163 * m.x37 + 0.07730602595 * m.x38 + 0.39988032723 * m.x39) + m.x50 == 0.194162178290626) m.e30 = Constraint(expr= -m.x44 * (2715.7894736842 / m.x45 + 5187 * m.x49 - 5187 * m.x50) / (4320 * m.x39 - 5187 * m.x50) + m.x48 == 0) m.e31 = Constraint(expr= exp(-0.029595 * m.x48) * (33.7894914681534 + m.x48) + m.x53 == 33.7894914681534) m.e32 = Constraint(expr= exp(-0.029595 * m.x48) + m.x54 == 1) m.e33 = Constraint(expr= m.b1 + m.b2 + m.b3 == 1) m.e34 = Constraint(expr= m.b1 * m.x26 <= 0) m.e35 = Constraint(expr= m.b2 * m.x26 >= 0) m.e36 = Constraint(expr= m.b2 * (m.x26 - m.x55) <= 0) m.e37 = Constraint(expr= m.b3 * (m.x26 - m.x55) >= 0) m.e38 = Constraint(expr= m.b4 + m.b5 + m.b6 == 1) m.e39 = Constraint(expr= m.b4 * m.x30 <= 0) m.e40 = Constraint(expr= m.b5 * m.x30 >= 0) m.e41 = Constraint(expr= m.b5 * (m.x30 - m.x55) <= 0) m.e42 = Constraint(expr= m.b6 * (m.x30 - m.x55) >= 0) m.e43 = Constraint(expr= m.b3 + m.b5 <= 1) m.e44 = Constraint(expr= m.b4 * (-1 + m.b1) >= 0) m.e45 = Constraint(expr= m.b7 + m.b8 + m.b9 == 1) m.e46 = Constraint(expr= m.b7 * (m.x44 - m.x55) <= 0) m.e47 = Constraint(expr= m.b8 * (m.x44 - m.x55) >= 0) m.e48 = Constraint(expr= m.b8 * (-200 + m.x44) <= 0) m.e49 = Constraint(expr= m.b9 * (-200 + m.x44) >= 0) m.e50 = Constraint(expr= m.b10 + m.b11 + m.b12 == 1) m.e51 = Constraint(expr= m.b10 * (m.x48 - m.x55) <= 0) m.e52 = Constraint(expr= m.b11 * (m.x48 - m.x55) >= 0) m.e53 = Constraint(expr= m.b11 * (-200 + m.x48) <= 0) m.e54 = Constraint(expr= m.b12 * (-200 + m.x48) >= 0) m.e55 = Constraint(expr= m.b9 + m.b11 <= 1) m.e56 = Constraint(expr= m.b10 * (-1 + m.b7) >= 0) m.e57 = Constraint(expr= -0.412152560083595 * m.b15 * m.x56 - 0.247360857459789 * m.b14 * m.x56 - m.b13 * (m.b1 * (11.34 * m.x56 / m.x28 + m.x26 * ( 0.854659090909091 / m.x27 - 11.34 / m.x28) * m.x57) + 0.0566666666666667 * m.x56 + m.b2 * (0.854659090909091 * m.x33 / m.x27 + (-11.34 * m.x33 + 11.34 * m.x56) / m.x28 + m.x26 * (0.854659090909091 / m.x27 - 11.34 / m.x28) * ( -m.x34 + m.x57)) + 0.854659090909091 * m.x56 / m.x27 * m.b3 + 0.01728 * m.b4 * m.x21 * m.x56 + m.b1 * m.b5 * (0.6 * (0.03458 * m.x32 * m.x35 + m.x26 * (0.0181052631578947 / m.x27 + 0.03458 * m.x31 - 0.03458 * m.x32) * m.x36) + 0.01728 * m.x21 * (-m.x35 + m.x56)) + m.b2 * m.b5 * (0.6 * (m.x33 * (0.0181052631578947 / m.x27 + 0.03458 * m.x31) + 0.03458 * m.x32 * ( -m.x33 + m.x35) + m.x26 * (0.0181052631578947 / m.x27 + 0.03458 * m.x31 - 0.03458 * m.x32) * (-m.x34 + m.x36)) + 0.01728 * m.x21 * (-m.x35 + m.x56)) + 0.6 * m.b1 * m.b6 * (0.03458 * m.x32 * m.x56 + m.x26 * ( 0.0181052631578947 / m.x27 + 0.03458 * m.x31 - 0.03458 * m.x32) * m.x57) + 0.6 * m.b2 * m.b6 * (m.x33 * (0.0181052631578947 / m.x27 + 0.03458 * m.x31) + 0.03458 * m.x32 * (-m.x33 + m.x56) + m.x26 * (0.0181052631578947 / m.x27 + 0.03458 * m.x31 - 0.03458 * m.x32) * (-m.x34 + m.x57)) + 0.6 * m.b3 * m.b6 * (0.0181052631578947 / m.x27 + 0.03458 * m.x31) * m.x56) + m.x29 == 0) m.e58 = Constraint(expr= m.b13 + m.b14 + m.b15 == 1) m.e59 = Constraint(expr= -m.b18 * (13.6674288800091 - 0.412152560083595 * m.x56) - m.b17 * (8.20275610163388 - 0.247360857459789 * m.x56) - m.b16 * (m.b7 * ((376.046780997472 - 11.34 * m.x56) / m.x46 + m.x44 * (0.854659090909091 / m.x45 - 11.34 / m.x46) * (0.997312113279821 - m.x57)) - 0.0566666666666667 * m.x56 + m.b8 * ((0.854659090909091 * m.x51 - 0.854659090909091 * m.x56) / m.x45 + (376.046780997472 - 11.34 * m.x51) / m.x46 + m.x44 * (0.854659090909091 / m.x45 - 11.34 / m.x46) * ( 0.997312113279821 - m.x52)) + m.b9 * (28.341428570246 - 0.854659090909091 * m.x56) / m.x45 + 0.01728 * m.b10 * m.x39 * (33.1610917987189 - m.x56) + m.b7 * m.b11 * (0.6 * (0.03458 * m.x50 * (m.x53 - m.x56) + m.x44 * ( 0.0181052631578947 / m.x45 + 0.03458 * m.x49 - 0.03458 * m.x50) * (m.x54 - m.x57)) + 0.01728 * m.x39 * (33.1610917987189 - m.x53)) + m.b8 * m.b11 * ( 0.6 * ((0.0181052631578947 / m.x45 + 0.03458 * m.x49) * (m.x51 - m.x56) + 0.03458 * m.x50 * (-m.x51 + m.x53) + m.x44 * (0.0181052631578947 / m.x45 + 0.03458 * m.x49 - 0.03458 * m.x50) * (-m.x52 + m.x54)) + 0.01728 * m.x39 * (33.1610917987189 - m.x53)) + 0.6 * m.b7 * m.b12 * (0.03458 * m.x50 * ( 33.1610917987189 - m.x56) + m.x44 * (0.0181052631578947 / m.x45 + 0.03458 * m.x49 - 0.03458 * m.x50) * (0.997312113279821 - m.x57)) + 0.6 * m.b8 * m.b12 * ((0.0181052631578947 / m.x45 + 0.03458 * m.x49) * (m.x51 - m.x56) + 0.03458 * m.x50 * (33.1610917987189 - m.x51) + m.x44 * ( 0.0181052631578947 / m.x45 + 0.03458 * m.x49 - 0.03458 * m.x50) * ( 0.997312113279821 - m.x52)) + 0.6 * m.b9 * m.b12 * (0.0181052631578947 / m.x45 + 0.03458 * m.x49) * (33.1610917987189 - m.x56) + 1.87912853526074) + m.x47 == 0) m.e60 = Constraint(expr= m.b16 + m.b17 + m.b18 == 1) m.e61 = Constraint(expr= exp(-0.029595 * m.x55) * (33.7894914681534 + m.x55) + m.x56 == 33.7894914681534) m.e62 = Constraint(expr= exp(-0.029595 * m.x55) + m.x57 == 1)