# 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])
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='--')
line2 = ax1.lines[-1]
ax1.legend((line1, line2), ('Lamb module', 'Dispersion'))
# Plot group velocity using the Lamb class.
fig2, ax2 = alum.plot_group_velocity(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).
ax2.get_legend().remove()
# Plot group velocity obtained by Dispersion.
line1 = ax2.lines[0]
for mode in df_vg.columns[::2]:
ax2.plot(df_vg[mode]*1e3, df_vg[mode+1]*1e3,
color = 'orange',
linestyle='--')