# Example: A 10 mm aluminum plate.
alum = Lamb(thickness=10,
nmodes_sym=5,
nmodes_antisym=5,
fd_max=10000,
vp_max=15000,
c_L=c_L,
c_S=c_S,
c_R=c_R,
material='Aluminum')
# Plot phase velocity, group velocity and wavenumber.
alum.plot_phase_velocity()
alum.plot_group_velocity()
alum.plot_wave_number()
# Plot wave structure (displacement profiles across thickness) for A0
# and S0 modes at different fd values.
alum.plot_wave_structure(mode='A0',
nrows=3,
ncols=2,
fd=[500, 1000, 1500, 2000, 2500, 3000])
alum.plot_wave_structure(mode='S0',
nrows=4,
ncols=2,
fd=[500, 1000, 1500, 2000, 2500, 3000, 3500, 4000])
header=None)
# Create an instance of the same material using the Lamb class.
alum = Lamb(thickness=1,
nmodes_sym=5,
nmodes_antisym=5,
fd_max=10000,
vp_max=15000,
c_L=6420,
c_S=3040)
# Plot phase velocity using the Lamb class.
fig1, ax1 = alum.plot_phase_velocity(material_velocities=False,
cutoff_frequencies=False,
sym_style={'color' : 'black'},
antisym_style={'color' : 'black'})
# Remove the legend that labels Symmetric and Antisymmetric modes
# (we are interested in labeling only Lamb module and Dispersion).
ax1.get_legend().remove()
# Plot phase velocity obtained by Dispersion.
line1 = ax1.lines[0]
for mode in df_vp.columns[::2]:
ax1.plot(df_vp[mode]*1e3, df_vp[mode+1]*1e3,
color = 'orange',
linestyle='--')
