%file name = tre_redun_opt.m
%redundant truss optimization
%main program
%Prof. Reyolando Brasil - May 2021
clc
clear
%problem data
Prob_data(1)=1;%gravity load P, KN
Prob_data(2)=1000;%density, kg/m³
Prob_data(3)=10;%alowble stress, KN/m²
Prob_data(4)=100;%Young's Modulus, KN/m²
%options
options=optimset ('LargeScale','off','TolCon',1e-8,'TolX',1e-8);
%Lower and upper bounds of design variables
Lb=[0 0];Ub=[1 1];
%initial design
x0=[0.1 0.1];
%optimization function call
[x,FunVal,ExitFlag,Output]=...
    fmincon('tre_redun_obj',x0,[],[],[],[],Lb,Ub,'tre_redun_con',options,Prob_data)
%