bartruss3_0.gms:
Reference:
- Ferris, M C, and Tin-Loi, F, On the Solution of a Minimum Weight Elastoplastic Problem Involving Displacement and Complementarity Constraints. Comp. Meth. in Appl. Mech. and Engng 174 (1999), 107-120.
Point:
<![CDATA[
* MPEC written by GAMS Convert at 11/06/01 17:01:56
*
* Equation counts
* Total E G L N X
* 29 23 6 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 36 36 0 0 0 0 0 0
* FX 6 6 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 96 58 38 0
*
* Solve m using MPEC minimizing objvar;
Variables objvar,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18
,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35
,x36;
Positive Variables x12,x13,x16,x17,x20,x21,x28,x29,x30,x31,x32,x33,x34,x35,x36;
Equations e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13,e14,e15,e16,e17,e18,e19
,e20,e21,e22,e23,e24,e25,e26,e27,e28,e29;
e1.. objvar - 500*x28 - 400*x29 - 500*x30 =E= 0;
e2.. - x28 + x29 =E= 0;
e3.. - x28 + x30 =E= 0;
e4.. x2 - 40*x28 =E= 0;
e5.. x3 - 50*x29 =E= 0;
e6.. x4 - 40*x30 =E= 0;
e7.. x5 - 50*x28 =E= 0;
e8.. x6 - 50*x28 =E= 0;
e9.. x7 - 50*x29 =E= 0;
e10.. x8 - 50*x29 =E= 0;
e11.. x9 - 50*x30 =E= 0;
e12.. x10 - 50*x30 =E= 0;
e13.. x11 - 5*x28 =E= 0;
e14.. x14 - 5*x28 =E= 0;
e15.. x15 - 6.25*x29 =E= 0;
e16.. x18 - 6.25*x29 =E= 0;
e17.. x19 - 5*x30 =E= 0;
e18.. x22 - 5*x30 =E= 0;
e19.. x2*(0.6*x26 + 0.8*x27) - x2*(x31 - x32) - x23 =E= 0;
e20.. x3*x27 - x3*(x33 - x34) - x24 =E= 0;
e21.. x4*(0.8*x27 - 0.6*x26) - x4*(x35 - x36) - x25 =E= 0;
e22.. 0.6*x23 - 0.6*x25 =E= 400;
e23.. 0.8*x23 + x24 + 0.8*x25 =E= 600;
e24.. x11*x31 + x12*x32 + x5 - x23 =G= 0;
e25.. x13*x31 + x14*x32 + x6 + x23 =G= 0;
e26.. x15*x33 + x16*x34 + x7 - x24 =G= 0;
e27.. x17*x33 + x18*x34 + x8 + x24 =G= 0;
e28.. x19*x35 + x20*x36 + x9 - x25 =G= 0;
e29.. x21*x35 + x22*x36 + x10 + x25 =G= 0;
* set non default bounds
x12.fx = 0;
x13.fx = 0;
x16.fx = 0;
x17.fx = 0;
x20.fx = 0;
x21.fx = 0;
x26.lo = -4; x26.up = 4;
x27.lo = -4; x27.up = 4;
* set non default levels
* set non default marginals
Model m / e1,e2,e3,e4.x2,e5.x3,e6.x4,e7.x5,e8.x6,e9.x7,e10.x8,e11.x9,e12.x10
,e13.x11,e14.x14,e15.x15,e16.x18,e17.x19,e18.x22,e19.x23,e20.x24
,e21.x25,e22.x26,e23.x27,e24.x31,e25.x32,e26.x33,e27.x34,e28.x35
,e29.x36 /;
m.limrow=0; m.limcol=0;
$if NOT '%gams.u1%' == '' $include '%gams.u1%'
Solve m using MPEC minimizing objvar;
]]></plaintext></body>