pyOMA Jupyter Notebook#

Run the cells below with Shift+Enter

%matplotlib widget

import matplotlib
import sys
import os
from pathlib import Path
import numpy as np
import matplotlib.pyplot as plt

please modify the path below to point to the code directory and run the cell

sys.path.append('/home/sima9999/git/pyOMA')

from pyOMA.core.PreProcessingTools import PreProcessSignals, GeometryProcessor, SignalPlot
from pyOMA.core.StabilDiagram import StabilCalc, StabilCluster, StabilPlot
from pyOMA.core.PostProcessingTools import MergePoSER
from pyOMA.core.PlotMSH import ModeShapePlot

from pyOMA.GUI.JupyterGUI import PlotMSHWeb, StabilGUIWeb, SnappingCursor, OutputWidgetHandler, ConfigGUIWeb

please modify the paths and setting below to point to the project and setup directories and run the cell

project_dir = Path('/home/sima9999/git/pyOMA/tests/files/')
nodes_file = project_dir / 'grid.txt'
lines_file = project_dir / 'lines.txt'
parent_child_file = ''

meas_name = 'measurement_1'
setup_dir = project_dir / meas_name
result_dir = setup_dir

setup_info_file = setup_dir / 'setup_info.txt'
meas_file = setup_dir / (meas_name + '.npy')
oma_conf_file = project_dir / 'ssi_config.txt'
chan_dofs_file = setup_dir / "channel_dofs.txt"

skip_existing=True
save_results=False

PreProcessSignals.load_measurement_file = np.load

display(ConfigGUIWeb(locals()))

geometry_data = GeometryProcessor.load_geometry(
    nodes_file=nodes_file,
    lines_file=lines_file,
    parent_childs_file=parent_child_file)

if not os.path.exists(setup_dir / 
        'prep_signals.npz') or not skip_existing:
    prep_signals = PreProcessSignals.init_from_config(
        conf_file=setup_info_file,
        meas_file=meas_file,
        chan_dofs_file=chan_dofs_file)

    if save_results:
        prep_signals.save_state(result_dir / 'prep_signals.npz')
else:
    prep_signals = PreProcessSignals.load_state(result_dir / 'prep_signals.npz')

prep_signals.decimate_signals(3)
prep_signals.decimate_signals(3)

INFO:pyOMA.core.PreProcessingTools:Loading results from  /home/sima9999/git/pyOMA/tests/files/measurement_1/prep_signals.npz
INFO:pyOMA.core.PreProcessingTools:Decimating signals by factor 3 to a sampling rate of 85.33333333333333 Hz
INFO:pyOMA.core.PreProcessingTools:Filtering signals in the band: None .. 34.13333333333333 with a 8 order cheby1 filter.
INFO:pyOMA.core.PreProcessingTools:Decimating signals by factor 3 to a sampling rate of 28.444444444444443 Hz
INFO:pyOMA.core.PreProcessingTools:Filtering signals in the band: None .. 11.377777777777776 with a 8 order cheby1 filter.

if not os.path.exists(result_dir /
        'modal_data.npz') or not skip_existing:

    modal_data = method.init_from_config(oma_conf_file, prep_signals)

    if save_results:
        
        modal_data.save_state(result_dir / 'modal_data.npz')
else:
    modal_data = method.load_state(
        result_dir / 'modal_data.npz', prep_signals)

INFO:pyOMA.core.SSICovRef:Loading results from  /home/sima9999/git/pyOMA/tests/files/measurement_1/modal_data.npz
INFO:pyOMA.core.SSICovRef:Covariance Matrices Built
INFO:pyOMA.core.SSICovRef:Modal Parameters Computed

if os.path.exists(result_dir 
                  / 'stabil_data.npz') and skip_existing:
    stabil_calc = StabilCalc.load_state(
        result_dir / 'stabil_data.npz', modal_data)
else:
    stabil_calc = StabilCalc(modal_data)
stabil_plot = StabilPlot(stabil_calc)

INFO:pyOMA.core.StabilDiagram:Now loading previous results from  /home/sima9999/git/pyOMA/tests/files/measurement_1/stabil_data.npz

vbox, snap_cursor =StabilGUIWeb(stabil_plot)
display(vbox)

INFO:pyOMA.core.StabilDiagram:Checking stabilisation criteria...
....................................................................................................

if save_results:
    stabil_calc.save_state(result_dir / 'stabil_data.npz')

mode_shape_plot = ModeShapePlot(amplitude=10,
                                prep_signals=prep_signals,
                                stabil_calc=stabil_calc,
                                geometry_data=geometry_data,
                                modal_data=modal_data)

display(PlotMSHWeb(mode_shape_plot))

'done'

'done'

This Page