run_nlbm.m

Benchmark Control Problems for Seismically Excited Nonlinear Building

%-----------------------------------------------------------------------
%|*********************************************************************|
%|*      Simulation Control Program for Nonlinear Benchmark Study     *|
%|*                                                                   *|
%|*                  University of Notre Dame                         *|
%|*                       November, 1999                              *|
%|*                                                                   *|
%|*               Coded by      Y.Ohtori                              *|
%|*                             R.E.Christenson                       *|
%|*               Supervised by B.F.Spencer, Jr.  	              *|
%|*********************************************************************|
%-----------------------------------------------------------------------

clear

%---------------------------------------
% --- Define Building for Simulation ---
%---------------------------------------
 No_bld     = 20;			% 20 Story Building

%---------------------------------------------
% --- Build Benchmark Model for Simulation ---
%---------------------------------------------
 sprintf(' ---> BUILDING MODEL FOR SIMULATION')
 Bld_NLBM
 
%----------------------------------------------
% --- Simulate El Centro Earthquake Records ---
%----------------------------------------------
 OPTIONS = simset('solver','ode5','FixedStep',dt_cal);
 isim = 0;	
 tf       = 100;            % Duration (sec)
 for intensity=0.5:0.5:1.5
  isim = isim+1;

  sprintf([' ---> SIMULATION: ElCentro ' num2str(intensity,2)])
  sim('Sim_NLBM',[0 tf],OPTIONS,[])
  load Mem_damage.out -ascii

  save(['BLD' num2str(No_bld) '_' num2str(isim) '.mat'],...
        'No_bld','yf','ye','Mem_damage','tf')
 end

%----------------------------------------------
% --- Simulate Hachinohe Earthquake Records ---
%----------------------------------------------
 delete_line('Sim_NLBM','Ground Accel/1','intensity/1')
 add_line   ('Sim_NLBM','Ground Accel/2','intensity/1')
 tf       = 100;            % Duration (sec)
 for intensity=0.5:0.5:1.5
  isim = isim+1;

  sprintf([' ---> SIMULATION: Hachinohe ' num2str(intensity,2)])
  sim('Sim_NLBM',[0 tf],OPTIONS,[])
  load Mem_damage.out -ascii

  save(['BLD' num2str(No_bld) '_' num2str(isim) '.mat'],...
        'No_bld','yf','ye','Mem_damage','tf')
 end

%-----------------------------------------------
% --- Simulate Northridge Earthquake Records ---
%-----------------------------------------------
 delete_line('Sim_NLBM','Ground Accel/2','intensity/1')
 add_line   ('Sim_NLBM','Ground Accel/3','intensity/1')
 tf       = 100;            % Duration (sec)
 for intensity=0.5:0.5:1.0
  isim = isim+1;

  sprintf([' ---> SIMULATION: Northridge ' num2str(intensity,2)])
  sim('Sim_NLBM',[0 tf],OPTIONS,[])
  load Mem_damage.out -ascii

  save(['BLD' num2str(No_bld) '_' num2str(isim) '.mat'],...
        'No_bld','yf','ye','Mem_damage','tf')
 end

%-----------------------------------------
% --- Simulate Kobe Earthquake Records ---
%-----------------------------------------
 delete_line('Sim_NLBM','Ground Accel/3','intensity/1')
 add_line   ('Sim_NLBM','Ground Accel/4','intensity/1')
 tf       = 180;            % Duration (sec)
 for intensity=0.5:0.5:1.0
  isim = isim+1;

  sprintf([' ---> SIMULATION: Kobe ' num2str(intensity,2)])
  sim('Sim_NLBM',[0 tf],OPTIONS,[])
  load Mem_damage.out -ascii

  save(['BLD' num2str(No_bld) '_' num2str(isim) '.mat'],...
        'No_bld','yf','ye','Mem_damage','tf')
end

%------------------------------------------------------
%--- RETURN SIMULATION TO ORIGINAL CONFIGURATION ------
%------------------------------------------------------
 delete_line('Sim_NLBM','Ground Accel/4','intensity/1')
 add_line   ('Sim_NLBM','Ground Accel/1','intensity/1')