root/OptimizingMATLABCode/Truss/paramSweep.m @ 10
| 10 | anderm8 | %% Truss Deflection - Parameter Sweep
|
|
% This is a parameter sweep study of the effect of the number of elements
|
|||
% and element cross sectional area on the displacement at the tip of a
|
|||
% cantilevered truss.
|
|||
%
|
|||
%
|
|||
% Copyright 2015 The MathWorks, Inc.
|
|||
%
|
|||
%% Initialize Problem
|
|||
% Number of elements and cross sectional area of each element
|
|||
nNum = 8;
|
|||
aNum = 8;
|
|||
% Vectors of number of elements and cross section areas to sweep
|
|||
nVals = (1:nNum)+10; % number of segments, start with 11
|
|||
aVals = linspace(100, 200, aNum); % cross sectional area
|
|||
% Grid of all combinations
|
|||
[nGrid, aGrid] = meshgrid(nVals, aVals);
|
|||
% Peak value results matrix
|
|||
peakVals = nan(nNum,aNum);
|
|||
%% Parameter Sweep
|
|||
t0 = tic;
|
|||
for ii = 1:numel(aGrid)
|
|||
% Solve ODE
|
|||
Y = trussCantilever(nGrid(ii),aGrid(ii));
|
|||
% Determine peak deflection in Y direction at the tip
|
|||
peakVals(ii) = max(Y(:,end));
|
|||
end
|
|||
toc(t0)
|
|||
%% Visualize results
|
|||
% Show the parameter sweep grid results
|
|||
visualizeParamSweep(nVals, aVals, peakVals);
|
|||
%% Animate Truss Deflection
|
|||
% Animate the truss over the simulation time at the largest deflection
|
|||
% Which combination yielded biggest deflection?
|
|||
[~, idx] = max(peakVals(:));
|
|||
% Get the full results at this location and animate it
|
|||
[Yr,bars,L,N,H] = trussCantilever(nGrid(idx),aGrid(idx));
|
|||
plotTruss(Yr,bars,L,N,H);
|