RunElastic2D_SteelBlock( )

function [ CatchAline ] = RunElastic2D_SteelBlock(  )
% SoundSim - 2D Elastic Simulation - For Academic Use Only
% -------------------------------------------------------------------------------------
% Limited Feature Beta Release of a SoundSim (www.SoundSim.com)
% Please Send all Questions, Suggestions, and Comments to [e-mail]
% No Error Checking - May become unstable if incorrect parameters are used
% -------------------------------------------------------------------------------------
% Steel Block Example -
% A 0.8 MHz Impulse on a top of a 7cm by 7cm steel block.
% The longitudinal, sheer, and Raleigh waves are visible
% The mode conversions are also visible as these waves 
% reflect from the edges of the block


%                                                                                        Aprox. Material Values From
%                                                                               Stress Waves in Solids (H. Kolsky) Dover 1963
% Material Properties                                                                   Aluminum  Steel  Copper  Glass 
p.cl = 5940;                         % Longitudinal Wave Speed (m/s)                     6320      5940   4560    5800
p.cs = 3220;                         % Sheer Wave Speed (m/s)                            3100      3220   2250    3350
p.den = 7800;                        % Material density (kg/m^3)                         2700      7800   8900    2500     
p.platethickness = 0.07;             % block thickness (meters)  (in x2-direction)
p.platelength    = 0.07;             % block length (meters)     (in x1-direction)

% Driving Transducer Properties
% -- These parameters adust the transducer position, size, frequency, and drive function length
% -- of the driving transducer.  The transducer is on the top of the plate.  
p.tpos = 0.025;                      % transducer position  (meters) - from left side of the block
p.tthickness = 0.001;                % transducer diameter  (meters) - this is 2D
p.tfreq = 800000;                    % frequency of transducer (Hertz) 
p.tpulselength = (.5)*(1/(p.tfreq)); % transducer pulse length (seconds)

% Catch Transducer Properties
% -- These parameters adjust the catch transducer location and size.  Again, this transducer  
% -- can only be placed on the top of the block
p.catchtpos = 0.05;                  % catch transducer position  (meters) - from left side of the plate
p.catchtthickness = 0.001;           % catch transducer diameter  (meters) - this is 2D

% Other Params
p.abc = 0;                           % thickness of absorbing boundary layer (in simulation units ds) Only on left and right
p.plotevery = 3;                     % update the plot every <plotevery> timesteps
p.SimulationTime =.000013;           % siulation runtime (seconds)

p.plotmode = 'abs_velocity';         % plot mode ('abs_velocity' = absolute velocity);
                                     %           ('x1_velocity'  = velocity in the x1 position (left and right))
                                     %           ('x2_velocity'  = velocity in the x2 position (up and down))
                             
CatchAline=SoundSim_ElasticEngine2D( p );  % Run the Simulation