http://www.k-wave.org/forum/topic/simulation-of-waves-on-a-metallic-plate-generated-by-force-from-an-acuator Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 17:27:20 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/simulation-of-waves-on-a-metallic-plate-generated-by-force-from-an-acuator#post-6230 Tue, 12 Dec 2017 22:20:51 +0000 Bradley Treeby 6230@http://www.k-wave.org/forum/ <p>Hi jrjulio18,</p> <p>I'd suggest using a velocity source. This is converted to a force source which represents the input of body forces per unit mass, and added to the momentum equation.</p> <p>Brad </p> http://www.k-wave.org/forum/topic/simulation-of-waves-on-a-metallic-plate-generated-by-force-from-an-acuator#post-6206 Thu, 09 Nov 2017 01:00:06 +0000 jrjulio18 6206@http://www.k-wave.org/forum/ <p>What kind of source should I generate in the simulation if i want to simulate a force hitting a plate? </p> http://www.k-wave.org/forum/topic/simulation-of-waves-on-a-metallic-plate-generated-by-force-from-an-acuator#post-6189 Sun, 05 Nov 2017 04:23:19 +0000 Fei Li 6189@http://www.k-wave.org/forum/ <p>k-Wave doesn't work with displacement directly, but you can record the particle velocity by setting sensor.record = {'u'} (see the Recording The Particle Velocity Example). You could then numerically integrate the velocity to calculate the displacement.</p> <p>---<a href="http://www.k-wave.org/forum/topic/displacement-of-a-viscous-medium#post-6006">'displacement in the viscous medium'</a> answered by Treeby </p> http://www.k-wave.org/forum/topic/simulation-of-waves-on-a-metallic-plate-generated-by-force-from-an-acuator#post-6188 Fri, 03 Nov 2017 23:02:33 +0000 jrjulio18 6188@http://www.k-wave.org/forum/ <p>Hi all,<br /> I am working on a project where we are trying to create a realistic simulation of the propagation of shear waves on a plate and to graph the displacement of the plate at multiple points. Currently I am using the source.s_mask to simulate a time varying pressure being applied to the plate at an small area and reading data using sensor.mask at different points on the grid. The simulation used is pstdElastic2D(kgrid, medium,source,sensor).Is the data recorded by the sensor Stress or Strain? and if is stress how can I manipulate that recorded data to a displacement data. Im sharing my code so maybe you guys can see it and run it. Have a good day.</p> <p>-Julio </p> <p>clear all ; close all;<br /> %% Create the Computational Grid %%<br /> dx = 0.2e-1; % [m]<br /> dy = 0.2e-1; % [m]<br /> Nx = round(2/dx)+1; % [grid points]<br /> Ny = round(2/dy)+1; % [grid points]<br /> kgrid = makeGrid(Nx, dx, Ny, dy);</p> <p>%% Medium Properties %%<br /> % define the compressional sound speed [m/s]<br /> medium.sound_speed_compression = 6420*ones(Nx, Ny);<br /> % define the shear sound speed [m/s]<br /> medium.sound_speed_shear = 3040*ones(Nx, Ny);<br /> % define the mass density [kg/m^3]<br /> medium.density = 2700*ones(Nx, Ny); </p> <p>%% Delta Sorce %%<br /> t = [linspace(-1,1+.39325,101)]; % Time Vector<br /> w = 2; % Triangle Width<br /> x = tripuls(t,w); % Sampled aperiodic triangle</p> <p>%% Square Source %%<br /> y = square(pi*t);</p> <p>%% Source (Presure) %%<br /> source.s_mask = zeros(Nx, Ny); %Creates my matrix for the Source<br /> source.s_mask(round(Nx/8),round(Ny/2)) =1;<br /> %source.s_mask(7*round(Nx/8),round(Ny/2)) =1;<br /> %source.s_mask(round(Nx/2),round(Ny/8)) =1;<br /> %source.s_mask(round(Nx/2),7*round(Ny/8)) =1;<br /> source.sxy =1*y; %Square and Delta Source </p> <p>%% Define my Sensors in the Grid %%<br /> sensor.mask = zeros(Nx, Ny);<br /> sensor.mask(round(Nx/8),round(Ny/2))=1 ; % Sensor at the 1/8 in the X-direction<br /> sensor.mask(round(Nx/4),round(Ny/2))=1 ; % Sensor at the 1/4 in the X-direction<br /> sensor.mask(round(Nx/2),round(Ny/2)) = 1 ;% Sensor at the center<br /> %sensor.mask(3*round(Nx/4),round(Ny/2))=1 ; % Sensor at the 3/4 in the X-direction </p> <p>%% run the simulation<br /> sensor_data = pstdElastic2D(kgrid, medium, source,sensor);<br /> [r, c] = size(sensor_data);<br /> %% Plot the Sensor Simulated Data %%<br /> for x = 1:4<br /> figure(1+r)<br /> hold on<br /> plot(1:c,sensor_data(x,1:c));<br /> xlabel('Time (usec)');ylabel('Displacement (um)');<br /> pause(.25)<br /> grid on;<br /> end</p> <p>-Julio </p>