2 files changed, 152 insertions, 0 deletions
diff --git a/code/fmincon/Continuous through circulation dryer.sce b/code/fmincon/Continuous through circulation dryer.sce
new file mode 100644
index 0000000..bbf38ab
--- /dev/null
+++ b/code/fmincon/Continuous through circulation dryer.sce
@@ -0,0 +1,66 @@
+//Reference: Rein Luus and Taina . I. Jaakola, Optimization of Nonlinear Functions Subject to Equality Constraints. Judicious Use of Elementary Calculus and Random Numbers,Ind. Eng. Chem. Process Des. Develop., Vol. 12, No. 3, 1973
+
+// Optimization of a continuous through circulation dyer
+// Maximize the production rate by considering the fluid velocity (X1) and bed length (X2) as the design parameters subjected to power constraint and the moisture content distribution constraint
+
+//======================================================================
+// Copyright (C) 2018 - IIT Bombay - FOSSEE
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+// Author: Remya Kommadath
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+//======================================================================
+clc; 
+// Objective function
+function f = ObjectiveFunction(X)
+    f = -(0.0064*X(1)*(1 - exp(-0.184*(X(1)^0.3)*X(2))));
+endfunction
+// Non linear equality constraints
+function [C,Ceq] = NLConstraints(X)
+    Ceq(1,1) = (3000 + X(1))*(X(1)^2)*X(2) - 1.2D13;
+    Ceq(1,2) = exp(0.184*X(1)^0.3*X(2)) - ln(4.1);
+    C = [];
+endfunction
+// Bounds of the problem
+lb = [0 0];
+ub = [];
+// Initial guess
+x0 = [1000 20];
+// Calling the solver
+[xopt,fopt,exitflag,output,lambda] = fmincon(ObjectiveFunction,x0,[],[],[],[],lb,[],NLConstraints)
+
+// Result representation
+select exitflag
+case 0
+    disp("Optimal Solution Found")
+    disp(x0,"Initial guess given to the solver")
+    disp(xopt',"The optimum solution obtained")
+    disp(fopt,"The optimum value of the objective function")
+case 1
+    disp(" Maximum Number of Iterations Exceeded. Output may not be optimal")
+    disp(x0,"Initial guess given to the solver")
+    disp(xopt',"The solution obtained")
+    disp(fopt,"The value of the objective function")
+case 2
+    disp("Maximum amount of CPU Time exceeded. Output may not be optimal")
+    disp(x0,"Initial guess given to the solver")
+    disp(xopt',"The solution obtained")
+    disp(fopt,"The value of the objective function")
+case 3
+    disp("Stop at Tiny Step")
+    disp(x0,"Initial guess given to the solver")
+    disp(xopt',"The solution obtained")
+    disp(fopt,"The value of the objective function")
+case 4
+    disp("Solved To Acceptable Level")
+    disp(x0,"Initial guess given to the solver")
+    disp(xopt',"The solution obtained")
+    disp(fopt,"The value of the objective function")
+case 5
+    disp("Converged to a point of local infeasibility")
+end
+
diff --git a/code/fmincon/Welded Beam Design.sce b/code/fmincon/Welded Beam Design.sce
new file mode 100644
index 0000000..b520f31
--- /dev/null
+++ b/code/fmincon/Welded Beam Design.sce
@@ -0,0 +1,86 @@
+//Reference: K. M. Ragsdell and D. T. Phillips,"Optimal Design of a Class of Welded Structures Using Geometric Programming",ASME Journal of Engineering for Industry,Vol.98, pp 1021-1025, 1976
+//A welded beam is designed for minimum cost subject to constraints on shear stress,bending stress in the beam,buckling load on the bar,end deflectionof the beam and the bound constraints. There are four design variables such as Weld thickness(inch), weld length (inch), bar thickness (inch), bar breadth(inch).
+//=====================================================================
+// Copyright (C) 2018 - IIT Bombay - FOSSEE
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+// Author: Remya Kommadath
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+//=====================================================================
+clc; 
+// Objective function
+function f = ObjectiveFunction(X)
+    f = 1.10471*X(1)^2*X(2) + 0.04811*X(3)*X(4)*(14+X(2));
+endfunction
+// Non linear equality and inequality constraints
+function [C,Ceq] = NLconstraints(X)
+    P = 6000; L = 14; E = 3*10^7; G = 12*10^6; tauMax = 13600;
+    sigmaMax = 30000; deltaMax = 0.25;
+
+    M = P*(L+(X(2)/2));
+    R = sqrt((X(2)^2/4) + ((X(1) + X(3))/2)^2);
+    J = 2*0.7071068*X(1)*X(2)*((X(2)^2/12)+((X(1)+X(3))/2)^2);
+    sigma = (6*P*L)/(X(3)^2*X(4));
+    delta = (4*P*L^3)/(E*X(3)^3*X(4));
+    Pc1 = (4.013*(sqrt(E*G*(X(3)^2*X(4)^6)/36)))/(L^2);
+    Pc2 = 1-(X(3)/(2*L))*sqrt(E/(4*G));
+    Pc = Pc1*Pc2;
+    tauPrime = P/(sqrt(2)*X(1)*X(2));
+    tauDprime = (M*R)/J;
+    tau = sqrt(tauPrime^2 + 2*tauPrime*tauDprime*(X(2)/(2*R))+tauDprime^2);
+
+    C(1) = tau - tauMax;
+    C(2) = sigma - sigmaMax;
+    C(3) = 0.10471*X(1)^2 + 0.04811*X(3)*X(4)*(14+X(2)) - 5;
+    C(4) = delta - deltaMax;
+    C(5) = P - Pc;
+    C = C';
+    Ceq = [];
+endfunction
+// Linear inequality constraints
+A = [1 0 0 -1;-1 0 0 0];
+b = [0 -0.125]';
+// Bounds of the problem
+lb = [0.1 0.1 0.1 0.1];
+ub = [2 10 10 2];
+// Initial guess
+x0 = rand(1,4).*(ub-lb);
+// Design parameters of the problem
+designParameters = {'Weld thickness(inch)','weld length (inch)','bar thickness (inch)','bar breadth(inch)'}
+inGuess = [designParameters; string(x0)]
+disp(inGuess,"Initial guess given to the solver")
+input("Press enter to proceed: ")
+// Calling the solver function
+[xopt,fopt,exitflag,output,lambda] = fmincon(ObjectiveFunction,x0,A,b,[],[],lb,ub,NLconstraints)
+
+// Result representation
+clc; 
+optSol = [designParameters; string(xopt')]
+select exitflag
+case 0
+    disp(optSol,"The optimum solution obtained")
+    disp(fopt,"The minimum cost for the weldment assembly is")
+case 1
+    disp(" Maximum Number of Iterations Exceeded. Output may not be optimal")
+    disp(optSol,"The solution obtained")
+    disp(fopt,"The minimum cost for the weldment assembly is")
+case 2
+    disp("Maximum amount of CPU Time exceeded. Output may not be optimal")
+    disp(optSol,"The solution obtained")
+    disp(fopt,"The minimum cost for the weldment assembly is")
+case 3
+    disp("Stop at Tiny Step")
+    disp(optSol,"The solution obtained")
+    disp(fopt,"The minimum cost for the weldment assembly is")
+case 4
+    disp("Solved To Acceptable Level")
+    disp(optSol,"The solution obtained")
+    disp(fopt,"The minimum cost for the weldment assembly is")
+case 5
+    disp("Converged to a point of local infeasibility")
+end
+disp(output)