def base_rotor_example():
"""Internal routine that create an example of a rotor, to be used in
the associated crack problems as a prerequisite.
This function returns an instance of a 6 DoF rotor, with a number of
components attached. As this is not the focus of the example here, but
only a requisite, see the example in "rotor assembly" for additional
information on the rotor object.
Returns
-------
rotor : ross.Rotor Object
An instance of a flexible 6 DoF rotor object.
Examples
--------
>>> rotor = base_rotor_example()
>>> rotor.Ip
0.015118294226367068
"""
steel2 = ross.Material(name="Steel", rho=7850, E=2.17e11, G_s=81.2e9)
# Rotor with 6 DoFs, with internal damping, with 10 shaft elements, 2 disks and 2 bearings.
i_d = 0
o_d = 0.019
n = 33
# fmt: off
L = np.array(
[0 , 25, 64, 104, 124, 143, 175, 207, 239, 271,
303, 335, 345, 355, 380, 408, 436, 466, 496, 526,
556, 586, 614, 647, 657, 667, 702, 737, 772, 807,
842, 862, 881, 914]
)/ 1000
# fmt: on
L = [L[i] - L[i - 1] for i in range(1, len(L))]
shaft_elem = [
ross.ShaftElement6DoF(
material=steel2,
L=l,
idl=i_d,
odl=o_d,
idr=i_d,
odr=o_d,
alpha=8.0501,
beta=1.0e-5,
rotary_inertia=True,
shear_effects=True,
)
for l in L
]
Id = 0.003844540885417
Ip = 0.007513248437500
disk0 = ross.DiskElement6DoF(n=12, m=2.6375, Id=Id, Ip=Ip)
disk1 = ross.DiskElement6DoF(n=24, m=2.6375, Id=Id, Ip=Ip)
kxx1 = 4.40e5
kyy1 = 4.6114e5
kzz = 0
cxx1 = 27.4
cyy1 = 2.505
czz = 0
kxx2 = 9.50e5
kyy2 = 1.09e8
cxx2 = 50.4
cyy2 = 100.4553
bearing0 = ross.BearingElement6DoF(
n=4, kxx=kxx1, kyy=kyy1, cxx=cxx1, cyy=cyy1, kzz=kzz, czz=czz
)
bearing1 = ross.BearingElement6DoF(
n=31, kxx=kxx2, kyy=kyy2, cxx=cxx2, cyy=cyy2, kzz=kzz, czz=czz
)
rotor = ross.Rotor(shaft_elem, [disk0, disk1], [bearing0, bearing1])
return rotor
import os
from pathlib import Path
from tempfile import tempdir
import numpy as np
import pytest
from numpy.testing import assert_allclose, assert_almost_equal
import ross as rs
from ross.defects.misalignment import MisalignmentFlex
from ross.units import Q_
steel2 = rs.Material(name="Steel", rho=7850, E=2.17e11, G_s=81.2e9)
# Rotor with 6 DoFs, with internal damping, with 10 shaft elements, 2 disks and 2 bearings.
i_d = 0
o_d = 0.019
n = 33
# fmt: off
L = np.array([
0, 25, 64, 104, 124, 143, 175, 207, 239, 271, 303, 335, 345, 355, 380, 408,
436, 466, 496, 526, 556, 586, 614, 647, 657, 667, 702, 737, 772, 807, 842,
862, 881, 914
]) / 1000
# fmt: on
L = [L[i] - L[i - 1] for i in range(1, len(L))]
shaft_elem = [
rs.ShaftElement6DoF(
import os
from pathlib import Path
from tempfile import tempdir
import numpy as np
import pytest
from numpy.testing import assert_allclose, assert_almost_equal
import ross as rs
from ross.defects.misalignment import MisalignmentFlex
from ross.units import Q_
steel2 = rs.Material(name="Steel",
rho=7850,
E=2.17e11,
Poisson=0.2992610837438423)
# Rotor with 6 DoFs, with internal damping, with 10 shaft elements, 2 disks and 2 bearings.
i_d = 0
o_d = 0.019
n = 33
# fmt: off
L = np.array([
0, 25, 64, 104, 124, 143, 175, 207, 239, 271, 303, 335, 345, 355, 380, 408,
436, 466, 496, 526, 556, 586, 614, 647, 657, 667, 702, 737, 772, 807, 842,
862, 881, 914
]) / 1000
# fmt: on
L = [L[i] - L[i - 1] for i in range(1, len(L))]
import os
from pathlib import Path
import ross as rs
import numpy as np
# MATERIAL
#
# E - Módulo de young
# G_s - Módulo de cisalhamento (shear)
# Poisson - Coeficiente de Poisson
# rho - Densidade do material
ABNT1020 = rs.Material(name="ABNT1020", rho=7850, E=2.05e11, Poisson=0.29)
ASTMA36 = rs.Material(name="ASTMA36", rho=7850, E=2e11, G_s=7.93e10)
# ===========================================
# EIXO
#
L = [0.3, 0.3, 0.052] # comprimento
i_d = 0 # diâmetro interno
o_d = 0.0127 # diâmetro externo
N_shaft = 3 # quantidade de elementos
shaft_elements = [
rs.ShaftElement(
n=i,
L=L[i],
idl=i_d,
odl=o_d,
material=ABNT1020,
