def rotor_example2(n_el=48):
"""
This function instantiate a overhung rotor similar to the example 5.9.9, page 218 (Dynamics of rotating machine,
FRISSWELL).
The following functions test_example2_w_equals0rpm() and test_example2_w_equals4000rpm() are the test functions of
this example.
P.S.: Overhung rotor.
:param n_el: number of shaft elements.
:return: Rotor object.
"""
shaft_elm = []
for i in range(n_el):
shaft_elm.append(
rs.ShaftElement(L=1.5 / n_el, material=steel, n=i, i_d=0,
o_d=0.05))
return rs.Rotor(
shaft_elm,
[
rs.DiskElement.from_geometry(
n=n_el, material=steel, width=0.07, i_d=0.05, o_d=0.35)
],
[
rs.BearingElement(n=0, kxx=10e6, kyy=10e6, cxx=0, cyy=0),
rs.BearingElement(
n=int((n_el / 1.5)), kxx=10e6, kyy=10e6, cxx=0, cyy=0),
],
)
def rotor_example1(n_el=48):
"""
This function instantiate a rotor similar to the example 5.9.1, page 206 (Dynamics of rotating machine, FRISSWELL)
The following functions test_example1_w_equals0rpm() and test_example1_w_equals4000rpm() are the test functions of
this example.
P.S.:Isotropic bearings.
:param w: speed of rotation.
:param n_el: number of shaft elements.
:return: Rotor object.
"""
shaft_elm = []
for i in range(n_el):
shaft_elm.append(
rs.ShaftElement(L=1.5 / n_el, material=steel, n=i, i_d=0,
o_d=0.05))
disk0 = rs.DiskElement.from_geometry(n=(n_el / 1.5) * 0.5,
material=steel,
width=0.07,
i_d=0.05,
o_d=0.28)
disk1 = rs.DiskElement.from_geometry(n=(n_el / 1.5),
material=steel,
width=0.07,
i_d=0.05,
o_d=0.35)
bearing0 = rs.BearingElement(n=0, kxx=1e6, kyy=1e6, cxx=0, cyy=0)
bearing1 = rs.BearingElement(n=n_el, kxx=1e6, kyy=1e6, cxx=0, cyy=0)
return rs.Rotor(shaft_elm, [disk0, disk1], [bearing0, bearing1])
def rotor_example7(w=0, n_el=48):
"""
This function instantiate a rotor similar to the example 5.9.10, page 219 (Dynamics of rotating machine, FRISSWELL)
The following functions test_example7_w_equals0rpm() and test_example7_w_equals4000rpm() are the test functions of
this example.
P.S.: Tapered shaft and damped bearings with anisotropic properties.
:param w: speed of rotation.
:param n_el: number of shaft elements.
:return: Rotor object.
"""
shaft_elm = []
for i in range(n_el):
shaft_elm.append(
rs.ShaftElement(L=1.5 / n_el,
material=steel,
n=i,
i_d=0,
o_d=0.025 + 0.015 * (i / n_el)))
disk0 = rs.DiskElement.from_geometry(n=(n_el / 2),
material=steel,
width=0.07,
i_d=0.05,
o_d=0.28)
bearing0 = rs.BearingElement(n=0, kxx=1e7, kyy=1e7, cxx=1e3, cyy=1e3)
bearing1 = rs.BearingElement(n=n_el, kxx=1e7, kyy=1e7, cxx=1e3, cyy=1e3)
return rs.Rotor(shaft_elm, [disk0], [bearing0, bearing1], w=w)