참조

Neculai Andrei, 슬롯 나라 기술을 사용한 비선형 최적화 애플리케이션, 스프링거 최적화 및 그 애플리케이션, 모델속도(5.2) 장기계공학의 응용, 2013

메인파일 : speed.gms

$onText
감속기의 무게를 최소화합니다.
감속기의 무게는 제약에 따라 최소화되어야 합니다.
기어 톱니의 굽힘 응력, 표면 응력, 가로 편향
샤프트의 응력과 샤프트의 응력.

Datseris, P., 휴리스틱을 통한 감속기의 무게 최소화
및 분해 기술. 메커니즘과 기계이론, vol.17, 1982,
255-262쪽.

Aguirre, A.H., Munoz Zavala, A.E., Villa Diharce, E., Botello Rionada, S.,
COPSO: PSO 알고리즘을 통한 제한된 최적화. Comunicacion Tecnica
아니요 I-07-04/22-02-2007. 멕시코 수학 연구 센터(CIMAT).
$offText

변수 x1, x2, x3, x4, x5, x6, x7, obj;

방정식 g1, g2, g3, g4, g5, g6, g7, g8, g9, g10, g11, g;

* 최소화할 목적 함수:

g.. obj =e= 0.7854*x1*POWER(x2,2)*(3.3333*POWER(x3,2) + 14.9334*x3 -
            43.0934) - 1.508*x1*(파워(x6,2)+파워(X7,2)) +
            7.4777*(파워(x6,3)+파워(x7,3))+
            0.7854*(x4*전력(x6,2)+x5*전력(x7,2));

* 제약:

g1.. 27/(x1*전력(x2,2)*x3) - 1 =l= 0;

g2.. 397.5/(x1*전력(x2,2)*전력(x3,2)) - 1 =l= 0;

g3.. (1.93*POWER(x4,3))/(x2*x3*POWER(x6,4)) - 1 =l= 0;

g4.. (1.93*POWER(x5,3))/(x2*x3*POWER(x7,4)) - 1 =l= 0;

g5.. (sqrt(POWER((745*x4)/(x2*x3),2)+16900000))/(110*POWER(x6,3)) - 1 =l= 0;

g6.. (sqrt(POWER((745*x5)/(x2*x3),2)+15750000))/(85*POWER(x7,3)) - 1 =l= 0;

g7..(x2*x3)/40 - 1 =l= 0;

g8..(5*x2)/x1 - 1 =l= 0;

g9.. x1/(12*x2) - 1 =l= 0;

g10.. (1.5*x6 + 1.9)/x4 - 1 =l= 0;

g11.. (1.1*x7 + 1.9)/x5 - 1 =l= 0;

* 변수의 한계

x1.lo = 2.6;     x1.up = 3.6;
x2.lo = 0.7;     x2.up = 0.8;
x3.lo = 17;     x3.up = 28;
x4.lo = 7.3;     x4.up = 8.3;
x5.lo = 7.8;     x5.up = 8.3;
x6.lo = 2.9;     x6.up = 3.9;
x7.lo = 5.0;     x7.up = 5.5;

모델 속도 /all/;

nlp를 사용하여 obj를 최소화하는 속도를 해결합니다.
* 종료 속도