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,