Contents

TWO STORY STEEL FRAME, PUSH-OVER ANALYSIS WITH FIXED GRAVITY LOADS AND LATERAL FORCES UNDER LOAD CONTROL

%  =========================================================================================
%  FEDEASLab - Release 2.6, July 2004
%  Matlab Finite Elements for Design, Evaluation and Analysis of Structures
%  Copyright(c) 1998-2004. The Regents of the University of California. All Rights Reserved.
%  Created by Professor Filip C. Filippou (filippou@ce.berkeley.edu)
%  Department of Civil and Environmental Engineering, UC Berkeley
%  =========================================================================================

Initialization: clear memory and define global variables

% all units in kip and inches
CleanStart

Create output file

IOW = Create_File (mfilename);

Create Model

Model_TwoStoryFrm
% echo input data of structural model to output file (optional)
Print_Model (Model,'Push-over analysis of two story steel frame');

Element properties

SimpleNLElemData
% print element properties (optional)
Structure ('data',Model,ElemData);

Loading (distributed loads only)

% define loading
for el=5:6 ElemData{el}.w = [0;-0.50]; end
for el=7:8 ElemData{el}.w = [0;-0.35]; end

GravLoading = Create_Loading (Model);

Incremental analysis for distributed element loading (single load step)

% initialize state
State = Initialize_State(Model,ElemData);
% initialize solution strategy parameters
SolStrat = Initialize_SolStrat;
% specify initial load increment (even though it is the same as the default value and could be omitted)
SolStrat.IncrStrat.Dlam0 = 1;
% initialize analysis sequence
[State SolStrat] = Initialize(Model,ElemData,GravLoading,State,SolStrat);
% apply load in one increment
[State SolStrat] = Increment(Model,ElemData,GravLoading,State,SolStrat);
% perform equilibrium iterations (we assume that convergence will occur!)
[State SolStrat] = Iterate  (Model,ElemData,GravLoading,State,SolStrat);
% update State
State = Update_State(Model,ElemData,State);
% determine resisting force vector
State   = Structure ('forc',Model,ElemData,State);
% set plot counter and store results for post-processing
k = 1;
Post(k) = Structure ('post',Model,ElemData,State);

 Iteration: 1  Relative work increment = 1.466268e-029
 Number of iterations = 1

2. Loading in sequence: horizontal forces

% specify nodal forces
Pe(2,1) =  20;
Pe(3,1) =  40;
Pe(5,1) =  20;
Pe(6,1) =  40;
LatLoading = Create_Loading (Model,Pe);

Incremental analysis for horizontal force pattern (load control is on)

%  (gravity forces are left on by not initializing State!)
% specify initial load increment and turn load control on (default value is 'no')
SolStrat.IncrStrat.Dlam0 = 0.40;
SolStrat.IncrStrat.LoadCtrl = 'yes';
SolStrat.IterStrat.LoadCtrl = 'yes';
% specify number of load steps
nostep = 20;
% initialize analysis sequence
[State SolStrat] = Initialize(Model,ElemData,LatLoading,State,SolStrat);

tic;
% for specified number of steps, Increment, Iterate and Update_State (we assume again convergence!)
for j=1:nostep
   [State SolStrat] = Increment(Model,ElemData,LatLoading,State,SolStrat);
   [State SolStrat] = Iterate  (Model,ElemData,LatLoading,State,SolStrat);
   State = Update_State(Model,ElemData,State);
   k = k+1;
   Post(k) = Structure ('post',Model,ElemData,State);
   % print results to output file
   Structure('prin',Model,ElemData,State);
end
toc;

 Iteration: 1  Relative work increment = 4.015228e-029
 Number of iterations = 1

 Iteration: 1  Relative work increment = 3.866660e-029
 Number of iterations = 1

 Iteration: 1  Relative work increment = 3.831656e-029
 Number of iterations = 1

 Iteration: 1  Relative work increment = 3.931251e-029
 Number of iterations = 1

 Iteration: 1  Relative work increment = 2.510484e-002
 Iteration: 2  Relative work increment = 1.303911e-030
 Number of iterations = 2

 Iteration: 1  Relative work increment = 1.986888e-003
 Iteration: 2  Relative work increment = 1.292956e-030
 Number of iterations = 2

 Iteration: 1  Relative work increment = 7.605427e-029
 Number of iterations = 1

 Iteration: 1  Relative work increment = 7.431033e-029
 Number of iterations = 1

 Iteration: 1  Relative work increment = 4.447378e-003
 Iteration: 2  Relative work increment = 6.083769e-030
 Number of iterations = 2

 Iteration: 1  Relative work increment = 3.078007e-002
 Iteration: 2  Relative work increment = 3.087329e-029
 Number of iterations = 2

 Iteration: 1  Relative work increment = 8.192344e-029
 Number of iterations = 1

 Iteration: 1  Relative work increment = 9.350141e-029
 Number of iterations = 1

 Iteration: 1  Relative work increment = 3.948878e-003
 Iteration: 2  Relative work increment = 1.325691e-024
 Number of iterations = 2

 Iteration: 1  Relative work increment = 1.892417e-024
 Number of iterations = 1

 Iteration: 1  Relative work increment = 2.075135e-024
 Number of iterations = 1

 Iteration: 1  Relative work increment = 1.494838e-024
 Number of iterations = 1

 Iteration: 1  Relative work increment = 1.669817e-024
 Number of iterations = 1

 Iteration: 1  Relative work increment = 2.278438e-024
 Number of iterations = 1

 Iteration: 1  Relative work increment = 1.870416e-024
 Number of iterations = 1

 Iteration: 1  Relative work increment = 2.176312e-024
 Number of iterations = 1
Elapsed time is 1.156000 seconds.

Post-processing

% extract displacements from Post
np = length(Post);
x = zeros(np,1);
y = zeros(np,1);
pltDOF = Model.DOF(1,6);
supDOF = [Model.DOF(1,1) Model.DOF(1,4)];
for k=1:np
   x(k) = Post(k).U(pltDOF);
   y(k) = -sum(Post(k).Pr(supDOF))/sum(Pe);
end

% plot force displacement relation in new window
fig = Create_Window(0.70,0.70);
ph1 = plot(x,y,'s-');
set (ph1,'MarkerSize',4,'MarkerFaceColor','b');
grid('on');
xlabel ('Horizontal roof displacement');
ylabel ('Load factor {\lambda}');
axis ([0 10 0 3]);
title ('Load factor-displacement');

% plot moment distribution at end of push-over
Create_Window(0.70,0.70);
title('Moment distribution and plastic hinge locations near incipient collapse');
Plot_Model(Model);
Plot_ForcDistr (Model,ElemData,Post(end),'Mz');
% show plastic hinge locations
Plot_PlasticHingeswPost(Model,Post(end));

% plot deformed shape of structure and plastic hinge locations
Create_Window(0.70,0.70);
Plot_Model(Model);
Structure('defo',Model,ElemData,State);
Plot_PlasticHingeswPost(Model,Post(end),State.U);
Resize_Plot;

% close output file
fclose(IOW);

Published with MATLAB® 7.0