''' t0 = time() # PREPARING MESH speed_of_sound = 331.2 # speed of sound at 0ºC density = 1.204 air = Fluid('air', density, speed_of_sound) steel = Material('Steel', 7860, young_modulus=210e9, poisson_ratio=0.3) offset = [0.005, 0.005] element_type = 'pipe_1' cross_section = CrossSection(0.05, 0.008, offset_y=offset[0], offset_z=offset[1], poisson_ratio=steel.poisson_ratio) cross_section.update_properties() preprocessor = Preprocessor() preprocessor.generate('examples/iges_files/tube_2.iges', 0.01) # preprocessor.set_acoustic_pressure_bc_by_node([50], 1e5 + 0j) # preprocessor.set_specific_impedance_bc_by_node(1086, speed_of_sound * density+ 0j) preprocessor.set_prescribed_dofs_bc_by_node([1136, 1236], np.zeros(6) + 0j) preprocessor.set_element_type_by_element('all', element_type) preprocessor.set_fluid_by_element('all', air) preprocessor.set_material_by_element('all', steel) preprocessor.set_cross_section_by_element('all', cross_section) preprocessor.set_structural_load_bc_by_node([50], np.array([1, 1, 1, 0, 0, 0]) + 0j) f_max = 200
from pulse.preprocessing.cross_section import CrossSection from pulse.preprocessing.material import Material from pulse.preprocessing.fluid import Fluid from pulse.preprocessing.preprocessor import Preprocessor from pulse.processing.assembly_acoustic import AssemblyAcoustic from pulse.processing.solution_acoustic import SolutionAcoustic from pulse.postprocessing.plot_acoustic_data import get_acoustic_frf # Fluid setup speed_of_sound = 343.21 density = 1.2041 air = Fluid('air', density, speed_of_sound) steel = Material('Steel', 7860, young_modulus=210e9, poisson_ratio=0.3) # Tube setup cross_section = CrossSection(0.05, 0.008, offset_y=0, offset_z=0) cross_section.update_properties() cross_section_expansion = CrossSection(0.288, 0.008, offset_y=0, offset_z=0) cross_section_expansion.update_properties() cross_section_branch = CrossSection(0.025, 0.004, offset_y=0, offset_z=0) cross_section_branch.update_properties() # Preprocessor init preprocessor = Preprocessor() preprocessor.generate('examples/iges_files/tube_2.iges', 0.01) preprocessor.set_material_by_element('all', steel) preprocessor.set_acoustic_pressure_bc_by_node(50, 1 + 0j) unflanged = True flanged = False if unflanged: preprocessor.set_radiation_impedance_bc_by_node(1086 , 1)
# Tube setup section_parameters = { "outer_diameter": 0.04859, "thickness": 0.003, "offset_y": 0, "offset_z": 0, "insulation_thickness": 0, "insulation_density": 0 } pipe_section_info = { "section_type_label": "Pipe section", "section_parameters": section_parameters } cross_section1 = CrossSection(pipe_section_info=pipe_section_info) cross_section1.update_properties() section_parameters = { "outer_diameter": 0.04044, "thickness": 0.003, "offset_y": 0, "offset_z": 0, "insulation_thickness": 0, "insulation_density": 0 } pipe_section_info = { "section_type_label": "Pipe section", "section_parameters": section_parameters } cross_section2 = CrossSection(pipe_section_info=pipe_section_info)