<jupyter_start><jupyter_text>Example 4 - Anisotropic Bearings.====In this example, we use the rotor seen in Example 5.9.2 from 'Dynamics of Rotating Machinery' by <NAME>, <NAME>, <NAME> & <NAME>, published by Cambridge University Press, 2010.Both bearings have a stiffness of 1 MN/m in the x direction and 0.8 MN/m in they direction. Calculate the eigenvalues and mode shapes at 0 and 4,000 rev/minand plot the natural frequency map for rotational speeds up to 4,500 rev/min.<jupyter_code>from bokeh.io import output_notebook, show
import ross as rs
import numpy as np
output_notebook()
#Classic Instantiation of the rotor
shaft_elements = []
bearing_seal_elements = []
disk_elements = []
Steel = rs.steel

for i in range(6):
    shaft_elements.append(rs.ShaftElement(L=0.25, material=Steel, n=i, i_d=0, o_d=0.05))
    
disk_elements.append(rs.DiskElement.from_geometry(n=2,
                                                  material=Steel, 
                                                  width=0.07,
                                                  i_d=0.05, 
                                                  o_d=0.28
                                                 )
                    )

disk_elements.append(rs.DiskElement.from_geometry(n=4,
                                                  material=Steel, 
                                                  width=0.07,
                                                  i_d=0.05, 
                                                  o_d=0.35
                                                 )
                    )
bearing_seal_elements.append(rs.BearingElement(n=0, kxx=1e6, kyy=.8e6, cxx=0, cyy=0))
bearing_seal_elements.append(rs.BearingElement(n=6, kxx=1e6, kyy=.8e6, cxx=0, cyy=0))

rotor592c = rs.Rotor(shaft_elements=shaft_elements,
                     bearing_seal_elements=bearing_seal_elements,
                     disk_elements=disk_elements,n_eigen = 12)

rotor592c.plot_rotor(plot_type='bokeh')
#From_section class method instantiation.
bearing_seal_elements = []
disk_elements = []
shaft_length_data = 3*[0.5]
i_d = 3*[0]
o_d = 3*[0.05]

disk_elements.append(rs.DiskElement.from_geometry(n=1,
                                                  material=Steel, 
                                                  width=0.07,
                                                  i_d=0.05, 
                                                  o_d=0.28
                                                 )
                    )

disk_elements.append(rs.DiskElement.from_geometry(n=2,
                                                  material=Steel, 
                                                  width=0.07,
                                                  i_d=0.05, 
                                                  o_d=0.35
                                                 )
                    )
bearing_seal_elements.append(rs.BearingElement(n=0, kxx=1e6, kyy=1e6, cxx=0, cyy=0))
bearing_seal_elements.append(rs.BearingElement(n=3, kxx=1e6, kyy=1e6, cxx=0, cyy=0))

rotor592fs = rs.Rotor.from_section(brg_seal_data=bearing_seal_elements,
                                   disk_data=disk_elements,leng_data=shaft_length_data,
                                   i_ds_data=i_d,o_ds_data=o_d                                  )
rotor592fs.plot_rotor(plot_type='bokeh')

#Obtaining results (wn is in rad/s)

print('Normal Instantiation =', rotor592c.wn)
print('\n')
print('From Section Instantiation =', rotor592fs.wn)
#Obtaining results for w=4000RPM (wn is in rad/s)

rotor592c.w=4000*np.pi/30

print('Normal Instantiation =', rotor592c.wn)<jupyter_output>Normal Instantiation = [ 82.32547398  86.86369902 239.64228361 287.24958074 583.48782916
 806.8872843 ]<jupyter_text>- Campbell Diagram<jupyter_code>campbell = rotor592c.run_campbell(np.linspace(0,4000*np.pi/30,100))
show(campbell.plot(plot_type='bokeh'))<jupyter_output><empty_output><jupyter_text>- Mode Shapes<jupyter_code>mode_shapes = rotor592c.run_mode_shapes()

for i in np.arange(0,5.1,1):
    mode_shapes.plot(mode=int(i))<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.09271037181996086, 'fraction_numerical': 0.0547945205479452, 'mean_word_length': 1.5749370277078085, 'pattern_counts': {'":': 0, '<': 16, '<?xml version=': 0, '>': 16, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 42, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '3', 'dir': 'jupyter-structured-clean-dedup', 'id': '27189', 'license': 'MIT', 'n_tokens_mistral': 1360, 'n_tokens_neox': 1204, 'n_words': 211, 'path': 'docs/examples/example_05_09_02.ipynb', 'repo_name': 'micanica/ross', 'size': '4100'}