# 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='--')