弹性波波场模拟matlab

弹性波波场模拟matlab程序SoundSim_ElasticEngine2D

function [ CatchAline ] = RunElastic2D_lambwave（ input_args ）
% 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 kevinrudd@SoundSim.com
% -------------------------------------------------------------------------------------
% Lamb Wave Example （Guided Elastic Wave）-
% An 5 cycle 0.68 MHz excitation on a 3.2mm thick Aluminum Plate.
% As the guided waves propagate down the plate The first symmetric （So）
% and anti-symmetric Modes （Ao） separate and are easily distinguishable.
%
%
%                                                                                         Aprox. Material Values From
%                                                                                Stress Waves in Solids （H. Kolsky） Dover 1963
% Material Properties                                                                   Aluminum  Steel  Copper  Glass
p.cl = 6320;                         % Longitudinal Wave Speed （m/s）                     6320      5940   4560    5800
p.cs = 3100;                         % Sheer Wave Speed （m/s）                            3100      3220   2250    3350
p.den = 2700;                        % Material density （kg/m^3）                         2700      7800   8900    2500 
p.platethickness = 0.0032;           % plate thickness （meters）  （in x2-direction）
p.platelength    = 0.30;             % plate 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 plate
p.tthickness = 0.001;                % transducer diameter  （meters） - this is 2D
p.tfreq = 680000;                    % 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 is 
% -- can only be placed on the top of the plate
p.catchtpos = 0.25;                  % 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 = 40;                          % thickness of absorbing boundary layer （in simulation units ds） Only on left and right
p.plotevery = 3;                     % update the plot every  timesteps
p.SimulationTime =.00012;            % siulation runtime （seconds）

p.plotmode = ‘x2_velocity‘;          % plot mode （‘abs_velocity‘ = absolute velocity）;
                                     %           （‘x1_velocity‘  = velocity in the x1 position （left and right））
                 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        3164  2006-07-31 16:20  RunElastic2D_lambwave.m
     文件        3230  2006-07-31 16:21  RunElastic2D_steelblock.m
     文件      136491  2006-07-31 16:19  SoundSim_2DElastic_Readme.pdf
     文件       13287  2006-07-31 16:19  SoundSim_ElasticEngine2D.p