%main program
%file name = port_wood_opt.m
%optimization of redundant wood planar portal frame
%Prof. Reyolando Brasil - May 2021
clc
clear
%portal frame data
h=3;%clumn height, m
L=6;%beam span, m
b=0.075;%section b dimension, m
P=15;%design load at bem mid span, KN
%options
options=optimset ('LargeScale','off','TolCon',1e-8,'TolX',1e-8);
%Variables lower and upper bounds
Lb=[0.15 0.15];Ub=[0.4 0.4];
%initial dimensions
x0=[0.15 0.15];
%optimization function call
[x,FunVal,ExitFlag,Output]=...
    fmincon('port_wood_obj',x0,[],[],[],[],Lb,Ub,'port_wood_con',options,h,L,b,P)
%