Harmonic acoustic analysis#

This example examines a harmonic acoustic analysis that uses surface velocity to determine the steady-state response of a structure and the surrounding fluid medium to loads and excitations that vary sinusoidally with time.

Import necessary libraries#

import os

from PIL import Image
import ansys.mechanical.core as mech
from ansys.mechanical.core.examples import delete_downloads, download_file
from matplotlib import image as mpimg
from matplotlib import pyplot as plt
from matplotlib.animation import FuncAnimation

Embed mechanical and set global variables

app = mech.App(version=241)
globals().update(mech.global_variables(app, True))
print(app)

cwd = os.path.join(os.getcwd(), "out")


def display_image(image_name):
    plt.figure(figsize=(16, 9))
    plt.imshow(mpimg.imread(os.path.join(cwd, image_name)))
    plt.xticks([])
    plt.yticks([])
    plt.axis("off")
    plt.show()
Ansys Mechanical [Ansys Mechanical Enterprise]
Product Version:241
Software build date: 11/27/2023 10:24:20

Configure graphics for image export#

ExtAPI.Graphics.Camera.SetSpecificViewOrientation(ViewOrientationType.Iso)
image_export_format = GraphicsImageExportFormat.PNG
settings_720p = Ansys.Mechanical.Graphics.GraphicsImageExportSettings()
settings_720p.Resolution = GraphicsResolutionType.EnhancedResolution
settings_720p.Background = GraphicsBackgroundType.White
settings_720p.Width = 1280
settings_720p.Height = 720
settings_720p.CurrentGraphicsDisplay = False
ExtAPI.Graphics.Camera.Rotate(180, CameraAxisType.ScreenY)

Download geometry and materials files#

geometry_path = download_file("C_GEOMETRY.agdb", "pymechanical", "embedding")
mat_path = download_file("Air-material.xml", "pymechanical", "embedding")

Import the geometry#

geometry_import = Model.GeometryImportGroup.AddGeometryImport()
geometry_import_format = (
    Ansys.Mechanical.DataModel.Enums.GeometryImportPreference.Format.Automatic
)
geometry_import_preferences = Ansys.ACT.Mechanical.Utilities.GeometryImportPreferences()
geometry_import_preferences.ProcessNamedSelections = True
geometry_import.Import(
    geometry_path, geometry_import_format, geometry_import_preferences
)


GEOM = Model.Geometry

solid1 = GEOM.Children[0]
solid2 = GEOM.Children[1]
solid3 = GEOM.Children[2]
solid4 = GEOM.Children[3]
solid5 = GEOM.Children[4]
solid6 = GEOM.Children[5]
solid7 = GEOM.Children[6]
solid8 = GEOM.Children[7]
solid9 = GEOM.Children[8]
solid10 = GEOM.Children[9]
solid11 = GEOM.Children[10]

solid1.Suppressed = True
solid2.Suppressed = True
solid3.Suppressed = True
solid4.Suppressed = True
solid5.Suppressed = True
solid7.Suppressed = True
solid10.Suppressed = True
solid11.Suppressed = True


ExtAPI.Graphics.Camera.SetFit()
ExtAPI.Graphics.ExportImage(
    os.path.join(cwd, "geometry.png"), image_export_format, settings_720p
)
display_image("geometry.png")
harmonic acoustics

Store all Variables necessary for analysis#

MESH = Model.Mesh
NS = Model.NamedSelections
CONN = Model.Connections
CS = Model.CoordinateSystems
MAT = Model.Materials

Setup the Analysis#

Add harmonic acoustics and unit system

Model.AddHarmonicAcousticAnalysis()
ExtAPI.Application.ActiveUnitSystem = (
    Ansys.ACT.Interfaces.Common.MechanicalUnitSystem.StandardMKS
)

Import and assign materials

MAT.Import(mat_path)
solid6.Material = "Air"
solid8.Material = "Air"
solid9.Material = "Air"

Create coordinate system

LCS1 = CS.AddCoordinateSystem()
LCS1.OriginX = Quantity("0 [mm]")
LCS1.OriginY = Quantity("0 [mm]")
LCS1.OriginZ = Quantity("0 [mm]")
LCS1.PrimaryAxisDefineBy = CoordinateSystemAlignmentType.GlobalZ

Generate mesh

MESH.ElementSize = Quantity("200 [mm]")
MESH.GenerateMesh()

Create named selections#

SF_Velo = Model.AddNamedSelection()
SF_Velo.ScopingMethod = GeometryDefineByType.Worksheet
SF_Velo.Name = "SF_Velo"
GEN_CRT1 = SF_Velo.GenerationCriteria
CRT1 = Ansys.ACT.Automation.Mechanical.NamedSelectionCriterion()
CRT1.Active = True
CRT1.Action = SelectionActionType.Add
CRT1.EntityType = SelectionType.GeoFace
CRT1.Criterion = SelectionCriterionType.Size
CRT1.Operator = SelectionOperatorType.Equal
CRT1.Value = Quantity("3e6 [mm^2]")
GEN_CRT1.Add(CRT1)
CRT2 = Ansys.ACT.Automation.Mechanical.NamedSelectionCriterion()
CRT2.Active = True
CRT2.Action = SelectionActionType.Filter
CRT2.EntityType = SelectionType.GeoFace
CRT2.Criterion = SelectionCriterionType.LocationZ
CRT2.Operator = SelectionOperatorType.Equal
CRT2.Value = Quantity("15000 [mm]")
GEN_CRT1.Add(CRT2)
SF_Velo.Activate()
SF_Velo.Generate()

ABS_Face = Model.AddNamedSelection()
ABS_Face.ScopingMethod = GeometryDefineByType.Worksheet
ABS_Face.Name = "ABS_Face"
GEN_CRT2 = ABS_Face.GenerationCriteria
CRT1 = Ansys.ACT.Automation.Mechanical.NamedSelectionCriterion()
CRT1.Active = True
CRT1.Action = SelectionActionType.Add
CRT1.EntityType = SelectionType.GeoFace
CRT1.Criterion = SelectionCriterionType.Size
CRT1.Operator = SelectionOperatorType.Equal
CRT1.Value = Quantity("1.5e6 [mm^2]")
GEN_CRT2.Add(CRT1)
CRT2 = Ansys.ACT.Automation.Mechanical.NamedSelectionCriterion()
CRT2.Active = True
CRT2.Action = SelectionActionType.Filter
CRT2.EntityType = SelectionType.GeoFace
CRT2.Criterion = SelectionCriterionType.LocationY
CRT2.Operator = SelectionOperatorType.Equal
CRT2.Value = Quantity("500 [mm]")
GEN_CRT2.Add(CRT2)
ABS_Face.Activate()
ABS_Face.Generate()

PRES_Face = Model.AddNamedSelection()
PRES_Face.ScopingMethod = GeometryDefineByType.Worksheet
PRES_Face.Name = "PRES_Face"
GEN_CRT3 = PRES_Face.GenerationCriteria
CRT1 = Ansys.ACT.Automation.Mechanical.NamedSelectionCriterion()
CRT1.Active = True
CRT1.Action = SelectionActionType.Add
CRT1.EntityType = SelectionType.GeoFace
CRT1.Criterion = SelectionCriterionType.Size
CRT1.Operator = SelectionOperatorType.Equal
CRT1.Value = Quantity("1.5e6 [mm^2]")
GEN_CRT3.Add(CRT1)
CRT2 = Ansys.ACT.Automation.Mechanical.NamedSelectionCriterion()
CRT2.Active = True
CRT2.Action = SelectionActionType.Filter
CRT2.EntityType = SelectionType.GeoFace
CRT2.Criterion = SelectionCriterionType.LocationY
CRT2.Operator = SelectionOperatorType.Equal
CRT2.Value = Quantity("4500 [mm]")
GEN_CRT3.Add(CRT2)
PRES_Face.Activate()
PRES_Face.Generate()

ACOUSTIC_Region = Model.AddNamedSelection()
ACOUSTIC_Region.ScopingMethod = GeometryDefineByType.Worksheet
ACOUSTIC_Region.Name = "ACOUSTIC_Region"
GEN_CRT4 = ACOUSTIC_Region.GenerationCriteria
CRT1 = Ansys.ACT.Automation.Mechanical.NamedSelectionCriterion()
CRT1.Active = True
CRT1.Action = SelectionActionType.Add
CRT1.EntityType = SelectionType.GeoBody
CRT1.Criterion = SelectionCriterionType.Type
CRT1.Operator = SelectionOperatorType.Equal
CRT1.Value = 8
GEN_CRT4.Add(CRT1)
ACOUSTIC_Region.Activate()
ACOUSTIC_Region.Generate()

Analysis settings#

ANALYSIS_SETTINGS = Model.Analyses[0].AnalysisSettings
ANALYSIS_SETTINGS.RangeMaximum = Quantity("100 [Hz]")
ANALYSIS_SETTINGS.SolutionIntervals = 50
ANALYSIS_SETTINGS.CalculateVelocity = True
ANALYSIS_SETTINGS.CalculateEnergy = True
ANALYSIS_SETTINGS.CalculateVolumeEnergy = True

Boundary conditions and load#

HARM_ACOUST = Model.Analyses[0]

Acoustic region

Acoustic_region = [x for x in HARM_ACOUST.Children if x.Name == "Acoustics Region"][0]
Acoustic_region.Location = ACOUSTIC_Region

Surface velocity

SURF_VEL = HARM_ACOUST.AddAcousticSurfaceVelocity()
SURF_VEL.Location = SF_Velo
SURF_VEL.Magnitude.Output.DiscreteValues = [Quantity("5000 [mm s-1]")]

Acoustic pressure

ACOUST_PRES = HARM_ACOUST.AddAcousticPressure()
ACOUST_PRES.Location = PRES_Face
ACOUST_PRES.Magnitude = Quantity("1.5e-7 [MPa]")

Acoustic absoption surface

ABSORP_SURF = HARM_ACOUST.AddAcousticAbsorptionSurface()
ABSORP_SURF.Location = ABS_Face
ABSORP_SURF.AbsorptionCoefficient.Output.DiscreteValues = [Quantity("0.02")]

HARM_ACOUST.Activate()
ExtAPI.Graphics.Camera.SetFit()
ExtAPI.Graphics.ExportImage(
    os.path.join(cwd, "bounday_conditions.png"), image_export_format, settings_720p
)
display_image("bounday_conditions.png")
harmonic acoustics

Add results#

SOLN = Model.Analyses[0].Solution

Acoustic pressure

ACOUST_PRES_RES1 = SOLN.AddAcousticPressureResult()
ACOUST_PRES_RES1.By = SetDriverStyle.ResultSet
ACOUST_PRES_RES1.SetNumber = 25

Acoustic velocity - total and directional

ACOUST_TOT_VEL1 = SOLN.AddAcousticTotalVelocityResult()
ACOUST_TOT_VEL1.Frequency = Quantity("50 [Hz]")

ACOUST_DIR_VEL1 = SOLN.AddAcousticDirectionalVelocityResult()
ACOUST_DIR_VEL1.Frequency = Quantity("50 [Hz]")
ACOUST_DIR_VEL1.CoordinateSystem = LCS1

ACOUST_DIR_VEL2 = SOLN.AddAcousticDirectionalVelocityResult()
ACOUST_DIR_VEL2.NormalOrientation = NormalOrientationType.ZAxis
ACOUST_DIR_VEL2.By = SetDriverStyle.ResultSet
ACOUST_DIR_VEL2.SetNumber = 25

Acoustic sound pressure and frequency bands

ACOUST_SPL = SOLN.AddAcousticSoundPressureLevel()
ACOUST_SPL.Frequency = Quantity("50 [Hz]")

ACOUST_A_SPL = SOLN.AddAcousticAWeightedSoundPressureLevel()
ACOUST_A_SPL.Frequency = Quantity("50 [Hz]")

ACOUST_FRQ_BAND_SPL = SOLN.AddAcousticFrequencyBandSPL()

A_FREQ_BAND_SPL = SOLN.AddAcousticFrequencyBandAWeightedSPL()

Z_VELO_RESP = SOLN.AddAcousticVelocityFrequencyResponse()
Z_VELO_RESP.NormalOrientation = NormalOrientationType.ZAxis
Z_VELO_RESP.Location = PRES_Face
Z_VELO_RESP.NormalOrientation = NormalOrientationType.ZAxis

Acoustic kinetic and potentional energy frequency response

KE_RESP = SOLN.AddAcousticKineticEnergyFrequencyResponse()
KE_RESP.Location = ABS_Face
KE_display = KE_RESP.TimeHistoryDisplay

PE_RESP = SOLN.AddAcousticPotentialEnergyFrequencyResponse()
PE_RESP.Location = ABS_Face
PE_display = PE_RESP.TimeHistoryDisplay

Acoustic total and directional velocity

ACOUST_TOT_VEL2 = SOLN.AddAcousticTotalVelocityResult()
ACOUST_TOT_VEL2.Location = PRES_Face
ACOUST_TOT_VEL2.Frequency = Quantity("30 [Hz]")
ACOUST_TOT_VEL2.Amplitude = True

ACOUST_DIR_VEL3 = SOLN.AddAcousticDirectionalVelocityResult()
ACOUST_DIR_VEL3.NormalOrientation = NormalOrientationType.ZAxis
ACOUST_DIR_VEL3.Location = PRES_Face
ACOUST_DIR_VEL3.Frequency = Quantity("10 [Hz]")
ACOUST_DIR_VEL3.Amplitude = True

ACOUST_KE = SOLN.AddAcousticKineticEnergy()
ACOUST_KE.Location = ABS_Face
ACOUST_KE.Frequency = Quantity("68 [Hz]")
ACOUST_KE.Amplitude = True

ACOUST_PE = SOLN.AddAcousticPotentialEnergy()
ACOUST_PE.Location = ABS_Face
ACOUST_PE.Frequency = Quantity("10 [Hz]")
ACOUST_PE.Amplitude = True

Solve#

SOLN.Solve(True)

Messages#

Messages = ExtAPI.Application.Messages
if Messages:
    for message in Messages:
        print(f"[{message.Severity}] {message.DisplayString}")
else:
    print("No [Info]/[Warning]/[Error] Messages")
[Warning] One or more Contact Regions are connected to Acoustic Physics Regions. To improve the solution, the contact settings may have been overwritten internally.  Refer to the notes for Acoustics Analysis in the Help System for more details.

Postprocessing#

Total acoustic pressure#

Tree.Activate([ACOUST_PRES_RES1])
ExtAPI.Graphics.ExportImage(
    os.path.join(cwd, "acou_pressure.png"), image_export_format, settings_720p
)
display_image("acou_pressure.png")
harmonic acoustics

Total acoustic velocity#

Tree.Activate([ACOUST_PRES_RES1])
ExtAPI.Graphics.ExportImage(
    os.path.join(cwd, "totalvelocity.png"), image_export_format, settings_720p
)
display_image("totalvelocity.png")
harmonic acoustics

Sound pressure level#

Tree.Activate([ACOUST_SPL])
ExtAPI.Graphics.ExportImage(
    os.path.join(cwd, "sound_pressure.png"), image_export_format, settings_720p
)
display_image("sound_pressure.png")
harmonic acoustics

Total velocity on pressure surface#

Tree.Activate([ACOUST_TOT_VEL2])
ExtAPI.Graphics.ExportImage(
    os.path.join(cwd, "totalvelocity_pressure.png"), image_export_format, settings_720p
)
display_image("totalvelocity_pressure.png")
harmonic acoustics

Kinetic energy on absorption face#

Tree.Activate([ACOUST_KE])
ExtAPI.Graphics.ExportImage(
    os.path.join(cwd, "kineticenergy.png"), image_export_format, settings_720p
)
display_image("kineticenergy.png")
harmonic acoustics

Total acoustic pressure animation#

animation_export_format = (
    Ansys.Mechanical.DataModel.Enums.GraphicsAnimationExportFormat.GIF
)
settings_720p = Ansys.Mechanical.Graphics.AnimationExportSettings()
settings_720p.Width = 1280
settings_720p.Height = 720

ACOUST_PRES_RES1.ExportAnimation(
    os.path.join(cwd, "press.gif"), animation_export_format, settings_720p
)
gif = Image.open(os.path.join(cwd, "press.gif"))
fig, ax = plt.subplots(figsize=(16, 9))
ax.axis("off")
img = ax.imshow(gif.convert("RGBA"))


def update(frame):
    gif.seek(frame)
    img.set_array(gif.convert("RGBA"))
    return [img]


ani = FuncAnimation(
    fig, update, frames=range(gif.n_frames), interval=200, repeat=True, blit=True
)
plt.show()

Display output file from solve#

def write_file_contents_to_console(path):
    """Write file contents to console."""
    with open(path, "rt") as file:
        for line in file:
            print(line, end="")


solve_path = HARM_ACOUST.WorkingDir
solve_out_path = os.path.join(solve_path, "solve.out")
if solve_out_path:
    write_file_contents_to_console(solve_out_path)
 Ansys Mechanical Enterprise


 *------------------------------------------------------------------*
 |                                                                  |
 |   W E L C O M E   T O   T H E   A N S Y S (R)  P R O G R A M     |
 |                                                                  |
 *------------------------------------------------------------------*




 ***************************************************************
 *         ANSYS MAPDL 2024 R1          LEGAL NOTICES          *
 ***************************************************************
 *                                                             *
 * Copyright 1971-2024 Ansys, Inc.  All rights reserved.       *
 * Unauthorized use, distribution or duplication is            *
 * prohibited.                                                 *
 *                                                             *
 * Ansys is a registered trademark of Ansys, Inc. or its       *
 * subsidiaries in the United States or other countries.       *
 * See the Ansys, Inc. online documentation or the Ansys, Inc. *
 * documentation CD or online help for the complete Legal      *
 * Notice.                                                     *
 *                                                             *
 ***************************************************************
 *                                                             *
 * THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION       *
 * INCLUDE TRADE SECRETS AND CONFIDENTIAL AND PROPRIETARY      *
 * PRODUCTS OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS.    *
 * The software products and documentation are furnished by    *
 * Ansys, Inc. or its subsidiaries under a software license    *
 * agreement that contains provisions concerning               *
 * non-disclosure, copying, length and nature of use,          *
 * compliance with exporting laws, warranties, disclaimers,    *
 * limitations of liability, and remedies, and other           *
 * provisions.  The software products and documentation may be *
 * used, disclosed, transferred, or copied only in accordance  *
 * with the terms and conditions of that software license      *
 * agreement.                                                  *
 *                                                             *
 * Ansys, Inc. is a UL registered                              *
 * ISO 9001:2015 company.                                      *
 *                                                             *
 ***************************************************************
 *                                                             *
 * This product is subject to U.S. laws governing export and   *
 * re-export.                                                  *
 *                                                             *
 * For U.S. Government users, except as specifically granted   *
 * by the Ansys, Inc. software license agreement, the use,     *
 * duplication, or disclosure by the United States Government  *
 * is subject to restrictions stated in the Ansys, Inc.        *
 * software license agreement and FAR 12.212 (for non-DOD      *
 * licenses).                                                  *
 *                                                             *
 ***************************************************************

 2024 R1

 Point Releases and Patches installed:

 Ansys, Inc. License Manager 2024 R1
 Structures 2024 R1
 LS-DYNA 2024 R1
 Mechanical Products 2024 R1


          *****  MAPDL COMMAND LINE ARGUMENTS  *****
  BATCH MODE REQUESTED (-b)    = NOLIST
  INPUT FILE COPY MODE (-c)    = COPY
  DISTRIBUTED MEMORY PARALLEL REQUESTED
       4 PARALLEL PROCESSES REQUESTED WITH SINGLE THREAD PER PROCESS
    TOTAL OF     4 CORES REQUESTED
  INPUT FILE NAME              = /github/home/.mw/Application Data/Ansys/v241/AnsysMech1B1F/Project_Mech_Files/HarmonicAcoustics/dummy.dat
  OUTPUT FILE NAME             = /github/home/.mw/Application Data/Ansys/v241/AnsysMech1B1F/Project_Mech_Files/HarmonicAcoustics/solve.out
  START-UP FILE MODE           = NOREAD
  STOP FILE MODE               = NOREAD

 RELEASE= 2024 R1              BUILD= 24.1      UP20231106   VERSION=LINUX x64
 CURRENT JOBNAME=file0  07:41:07  MAY 06, 2024 CP=      0.229


 PARAMETER _DS_PROGRESS =     999.0000000

 /INPUT FILE= ds.dat  LINE=       0



 *** NOTE ***                            CP =       0.296   TIME= 07:41:07
 The /CONFIG,NOELDB command is not valid in a distributed memory
 parallel solution.  Command is ignored.

 *GET  _WALLSTRT  FROM  ACTI  ITEM=TIME WALL  VALUE=  7.68527778

 TITLE=
 --Harmonic Acoustics


 SET PARAMETER DIMENSIONS ON  _WB_PROJECTSCRATCH_DIR
  TYPE=STRI  DIMENSIONS=      248        1        1

 PARAMETER _WB_PROJECTSCRATCH_DIR(1) = /github/home/.mw/Application Data/Ansys/v241/AnsysMech1B1F/Project_Mech_Files/HarmonicAcoustics/

 SET PARAMETER DIMENSIONS ON  _WB_SOLVERFILES_DIR
  TYPE=STRI  DIMENSIONS=      248        1        1

 PARAMETER _WB_SOLVERFILES_DIR(1) = /github/home/.mw/Application Data/Ansys/v241/AnsysMech1B1F/Project_Mech_Files/HarmonicAcoustics/

 SET PARAMETER DIMENSIONS ON  _WB_USERFILES_DIR
  TYPE=STRI  DIMENSIONS=      248        1        1

 PARAMETER _WB_USERFILES_DIR(1) = /github/home/.mw/Application Data/Ansys/v241/AnsysMech1B1F/Project_Mech_Files/UserFiles/
 --- Data in consistent MKS units. See Solving Units in the help system for more

 MKS UNITS SPECIFIED FOR INTERNAL
  LENGTH        (l)  = METER (M)
  MASS          (M)  = KILOGRAM (KG)
  TIME          (t)  = SECOND (SEC)
  TEMPERATURE   (T)  = CELSIUS (C)
  TOFFSET            = 273.0
  CHARGE        (Q)  = COULOMB
  FORCE         (f)  = NEWTON (N) (KG-M/SEC2)
  HEAT               = JOULE (N-M)

  PRESSURE           = PASCAL (NEWTON/M**2)
  ENERGY        (W)  = JOULE (N-M)
  POWER         (P)  = WATT (N-M/SEC)
  CURRENT       (i)  = AMPERE (COULOMBS/SEC)
  CAPACITANCE   (C)  = FARAD
  INDUCTANCE    (L)  = HENRY
  MAGNETIC FLUX      = WEBER
  RESISTANCE    (R)  = OHM
  ELECTRIC POTENTIAL = VOLT

 INPUT  UNITS ARE ALSO SET TO MKS

 *** MAPDL - ENGINEERING ANALYSIS SYSTEM  RELEASE 2024 R1          24.1     ***
 Ansys Mechanical Enterprise
 00000000  VERSION=LINUX x64     07:41:07  MAY 06, 2024 CP=      0.300

 --Harmonic Acoustics



          ***** MAPDL ANALYSIS DEFINITION (PREP7) *****
 *********** Nodes for the whole assembly ***********
 *********** Elements for Body 1 "Solid" ***********
 *********** Elements for Body 2 "Solid" ***********
 *********** Elements for Body 3 "Solid" ***********
 *********** Send User Defined Coordinate System(s) ***********
 *********** Set Reference Temperature ***********
 *********** Send Materials ***********
 *********** Create Contact "Contact Region 6" ***********
             Real Constant Set For Above Contact Is 5 & 4
 *********** Create Contact "Contact Region 9" ***********
             Real Constant Set For Above Contact Is 7 & 6
 *********** Send Named Selection as Node Component ***********
 *********** Send Named Selection as Node Component ***********
 *********** Send Named Selection as Node Component ***********
 *********** Send Named Selection as Element Component ***********
 *********** Create Acoustic Pressure ***********
 *********** Create Absorption Surface ***********
 *********** Begin Program Controlled Acoustic Far-field Radiation Surface *****
 *********** End Program Controlled Acoustic Far-field Radiation Surface *******
 *********** Create Surface Velocity ***********


 ***** ROUTINE COMPLETED *****  CP =         0.438


 --- Number of total nodes = 9687
 --- Number of contact elements = 320
 --- Number of spring elements = 0
 --- Number of bearing elements = 0
 --- Number of solid elements = 1840
 --- Number of condensed parts = 0
 --- Number of total elements = 2160

 *GET  _WALLBSOL  FROM  ACTI  ITEM=TIME WALL  VALUE=  7.68527778

 *****  MAPDL SOLUTION ROUTINE  *****


 PERFORM A HARMONIC ANALYSIS
  THIS WILL BE A NEW ANALYSIS


 PERFORM A FULL HARMONIC RESPONSE ANALYSIS

 THERMAL STRAINS ARE NOT INCLUDED IN THE LOAD VECTOR.

 STEP BOUNDARY CONDITION KEY= 1

 HARMONIC FREQUENCY RANGE - END=  100.00      BEGIN=  0.0000

 USE    50 SUBSTEP(S) THIS LOAD STEP FOR ALL  DEGREES OF FREEDOM

 *** NOTE ***                            CP =       0.443   TIME= 07:41:07
 The acoustic solver is set to the total field formulation.

 ERASE THE CURRENT DATABASE OUTPUT CONTROL TABLE.


 WRITE ALL  ITEMS TO THE DATABASE WITH A FREQUENCY OF NONE
   FOR ALL APPLICABLE ENTITIES

 WRITE NSOL ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL
   FOR ALL APPLICABLE ENTITIES

 WRITE EANG ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL
   FOR ALL APPLICABLE ENTITIES

 WRITE ETMP ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL
   FOR ALL APPLICABLE ENTITIES

 WRITE VENG ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL
   FOR ALL APPLICABLE ENTITIES

 PRINTOUT RESUMED BY /GOP

 WRITE MISC ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL
   FOR THE ENTITIES DEFINED BY COMPONENT _ELMISC

 WRITE FGRA ITEMS TO THE DATABASE WITH A FREQUENCY OF ALL
   FOR ALL APPLICABLE ENTITIES

 *GET  ANSINTER_  FROM  ACTI  ITEM=INT        VALUE=  0.00000000

 *IF  ANSINTER_                         ( =   0.00000     )  NE
      0                                 ( =   0.00000     )  THEN

 *ENDIF

 *** NOTE ***                            CP =       0.498   TIME= 07:41:07
 The automatic domain decomposition logic has selected the FREQ domain
 decomposition method with 1 processes per frequency solution.

 *****  MAPDL SOLVE    COMMAND  *****



   D I S T R I B U T E D   D O M A I N   D E C O M P O S E R

 ...Number of frequency solutions:  50
 ...Decompose to 4 frequency domains (with 1 processes per domain)


 *** WARNING ***                         CP =       0.671   TIME= 07:41:07
 Element shape checking is currently inactive.  Issue SHPP,ON or
 SHPP,WARN to reactivate, if desired.

 *** NOTE ***                            CP =       0.684   TIME= 07:41:07
 The model data was checked and warning messages were found.
  Please review output or errors file ( /github/home/.mw/Application
 Data/Ansys/v241/AnsysMech1B1F/Project_Mech_Files/HarmonicAcoustics/file
 e0.err ) for these warning messages.

 *** MAPDL - ENGINEERING ANALYSIS SYSTEM  RELEASE 2024 R1          24.1     ***
 Ansys Mechanical Enterprise
 00000000  VERSION=LINUX x64     07:41:07  MAY 06, 2024 CP=      0.685

 --Harmonic Acoustics



                       S O L U T I O N   O P T I O N S

   PROBLEM DIMENSIONALITY. . . . . . . . . . . . .3-D
   DEGREES OF FREEDOM. . . . . . PRES
   ANALYSIS TYPE . . . . . . . . . . . . . . . . .HARMONIC
      SOLUTION METHOD. . . . . . . . . . . . . . .AUTO
   OFFSET TEMPERATURE FROM ABSOLUTE ZERO . . . . .  273.15
   THERMAL EXPANSION . . . . . . . . . . . . . . .OFF
   COMPLEX DISPLACEMENT PRINT OPTION . . . . . . .REAL AND IMAGINARY
   GLOBALLY ASSEMBLED MATRIX . . . . . . . . . . .SYMMETRIC

 *** NOTE ***                            CP =       0.695   TIME= 07:41:07
 The Solution Control Option is only valid for single field structural,
 single field thermal, single field diffusion, coupled
 thermal-diffusion analyses and coupled-field analyses with structural
 degrees of freedom.  The SOLCONTROL,ON command (if present) has been
 de-activated.

 *** NOTE ***                            CP =       0.719   TIME= 07:41:07
 A target element 3397 is attached to a fluid element 1836.  You may
 consider to flip the contact and target surfaces.

 *** NOTE ***                            CP =       1.525   TIME= 07:41:07
 Symmetric Deformable- deformable contact pair identified by real
 constant set 4 and contact element type 4 has been set up.  The
 companion pair has real constant set ID 5.  Both pairs should have the
 same behavior.
 *WARNING*: The contact pairs have similar mesh patterns.  MAPDL will
 keep the current pair and deactivate its companion pair.
 Pure pore fluid contact is activated.
 Contact algorithm: MPC based approach

 *** NOTE ***                            CP =       1.525   TIME= 07:41:07
 Contact related postprocess items (ETABLE, pressure ...) are not
 available.
 Contact detection at: nodal point (Dual shape function based)
 MPC will be built internally to handle bonded contact.
 Average contact surface length               0.19070
 Average contact pair depth                   0.17708
 Average target surface length                0.19375
 Default pinball region factor PINB           0.25000
 The resulting pinball region                 0.44271E-01
 Initial penetration/gap is excluded.
 Bonded contact (always) is defined.

 *** NOTE ***                            CP =       1.525   TIME= 07:41:07
 Max.  Initial penetration 2.852608047E-15 was detected between contact
 element 3333 and target element 3413.
 ****************************************


 *** WARNING ***                         CP =       1.525   TIME= 07:41:07
 Element shape checking is currently inactive.  Issue SHPP,ON or
 SHPP,WARN to reactivate, if desired.

 *** WARNING ***                         CP =       1.525   TIME= 07:41:07
 Element shape checking is currently inactive.  Issue SHPP,ON or
 SHPP,WARN to reactivate, if desired.

 *** WARNING ***                         CP =       1.525   TIME= 07:41:07
 Element shape checking is currently inactive.  Issue SHPP,ON or
 SHPP,WARN to reactivate, if desired.

 *** NOTE ***                            CP =       1.525   TIME= 07:41:07
 Symmetric Deformable- deformable contact pair identified by real
 constant set 5 and contact element type 4 has been set up.  The
 companion pair has real constant set ID 4.  Both pairs should have the
 same behavior.
 MAPDL will deactivate the current pair and keep its companion pair,
 resulting in asymmetric contact.
 Pure pore fluid contact is activated.
 Contact algorithm: MPC based approach

 *** NOTE ***                            CP =       1.525   TIME= 07:41:07
 Contact related postprocess items (ETABLE, pressure ...) are not
 available.
 Contact detection at: nodal point (Dual shape function based)
 MPC will be built internally to handle bonded contact.
 Average contact surface length               0.19070
 Average contact pair depth                   0.20000
 Average target surface length                0.19375
 Default pinball region factor PINB           0.25000
 The resulting pinball region                 0.50000E-01
 Initial penetration/gap is excluded.
 Bonded contact (always) is defined.

 *** NOTE ***                            CP =       1.525   TIME= 07:41:07
 Max.  Initial penetration 2.404977522E-15 was detected between contact
 element 3381 and target element 3285.
 ****************************************


 *** NOTE ***                            CP =       1.525   TIME= 07:41:07
 Symmetric Deformable- deformable contact pair identified by real
 constant set 6 and contact element type 6 has been set up.  The
 companion pair has real constant set ID 7.  Both pairs should have the
 same behavior.
 *WARNING*: The contact pairs have similar mesh patterns.  MAPDL will
 keep the current pair and deactivate its companion pair.
 Pure pore fluid contact is activated.
 Contact algorithm: MPC based approach

 *** NOTE ***                            CP =       1.525   TIME= 07:41:07
 Contact related postprocess items (ETABLE, pressure ...) are not
 available.
 Contact detection at: nodal point (Dual shape function based)
 MPC will be built internally to handle bonded contact.
 Average contact surface length               0.19070
 Average contact pair depth                   0.17708
 Average target surface length                0.19375
 Default pinball region factor PINB           0.25000
 The resulting pinball region                 0.44271E-01
 Initial penetration/gap is excluded.
 Bonded contact (always) is defined.

 *** NOTE ***                            CP =       1.525   TIME= 07:41:07
 Max.  Initial penetration 1.025113101E-15 was detected between contact
 element 3502 and target element 3587.
 ****************************************


 *** NOTE ***                            CP =       1.525   TIME= 07:41:07
 Symmetric Deformable- deformable contact pair identified by real
 constant set 7 and contact element type 6 has been set up.  The
 companion pair has real constant set ID 6.  Both pairs should have the
 same behavior.
 MAPDL will deactivate the current pair and keep its companion pair,
 resulting in asymmetric contact.
 Pure pore fluid contact is activated.
 Contact algorithm: MPC based approach

 *** NOTE ***                            CP =       1.525   TIME= 07:41:07
 Contact related postprocess items (ETABLE, pressure ...) are not
 available.
 Contact detection at: nodal point (Dual shape function based)
 MPC will be built internally to handle bonded contact.
 Average contact surface length               0.19070
 Average contact pair depth                   0.20000
 Average target surface length                0.19375
 Default pinball region factor PINB           0.25000
 The resulting pinball region                 0.50000E-01
 Initial penetration/gap is excluded.
 Bonded contact (always) is defined.

 *** NOTE ***                            CP =       1.525   TIME= 07:41:07
 Max.  Initial penetration 1.419491659E-15 was detected between contact
 element 3525 and target element 3452.
 ****************************************




                      L O A D   S T E P   O P T I O N S

   LOAD STEP NUMBER. . . . . . . . . . . . . . . .     1
   FREQUENCY RANGE . . . . . . . . . . . . . . . .  0.0000     TO   100.00
   NUMBER OF SUBSTEPS. . . . . . . . . . . . . . .    50
   STEP CHANGE BOUNDARY CONDITIONS . . . . . . . .   YES
   PRINT OUTPUT CONTROLS . . . . . . . . . . . . .NO PRINTOUT
   DATABASE OUTPUT CONTROLS
      ITEM     FREQUENCY   COMPONENT
       ALL       NONE
      NSOL        ALL
      EANG        ALL
      ETMP        ALL
      VENG        ALL
      MISC        ALL       _ELMISC
      FGRA        ALL

   AUTO SELECTION OF VT FOR FREQUENCY SWEEP. . . . NO

 *** NOTE ***                            CP =       2.461   TIME= 07:41:07
 Symmetric Deformable- deformable contact pair identified by real
 constant set 4 and contact element type 4 has been set up.  The
 companion pair has real constant set ID 5.  Both pairs should have the
 same behavior.
 MAPDL will keep the current pair and deactivate its companion pair,
 resulting in asymmetric contact.
 Pure pore fluid contact is activated.
 Contact algorithm: MPC based approach

 *** NOTE ***                            CP =       2.461   TIME= 07:41:07
 Contact related postprocess items (ETABLE, pressure ...) are not
 available.
 Contact detection at: nodal point (Dual shape function based)
 MPC will be built internally to handle bonded contact.
 Average contact surface length               0.19070
 Average contact pair depth                   0.17708
 Average target surface length                0.19375
 Default pinball region factor PINB           0.25000
 The resulting pinball region                 0.44271E-01
 Initial penetration/gap is excluded.
 Bonded contact (always) is defined.

 *** NOTE ***                            CP =       2.461   TIME= 07:41:07
 Max.  Initial penetration 2.852608047E-15 was detected between contact
 element 3333 and target element 3413.
 ****************************************


 *** NOTE ***                            CP =       2.461   TIME= 07:41:07
 Symmetric Deformable- deformable contact pair identified by real
 constant set 5 and contact element type 4 has been set up.  The
 companion pair has real constant set ID 4.  Both pairs should have the
 same behavior.
 MAPDL will deactivate the current pair and keep its companion pair,
 resulting in asymmetric contact.
 Pure pore fluid contact is activated.
 Contact algorithm: MPC based approach

 *** NOTE ***                            CP =       2.461   TIME= 07:41:07
 Contact related postprocess items (ETABLE, pressure ...) are not
 available.
 Contact detection at: nodal point (Dual shape function based)
 MPC will be built internally to handle bonded contact.
 Average contact surface length               0.19070
 Average contact pair depth                   0.20000
 Average target surface length                0.19375
 Default pinball region factor PINB           0.25000
 The resulting pinball region                 0.50000E-01
 Initial penetration/gap is excluded.
 Bonded contact (always) is defined.

 *** NOTE ***                            CP =       2.461   TIME= 07:41:07
 Max.  Initial penetration 2.404977522E-15 was detected between contact
 element 3381 and target element 3285.
 ****************************************


 *** NOTE ***                            CP =       2.461   TIME= 07:41:07
 Symmetric Deformable- deformable contact pair identified by real
 constant set 6 and contact element type 6 has been set up.  The
 companion pair has real constant set ID 7.  Both pairs should have the
 same behavior.
 MAPDL will keep the current pair and deactivate its companion pair,
 resulting in asymmetric contact.
 Pure pore fluid contact is activated.
 Contact algorithm: MPC based approach

 *** NOTE ***                            CP =       2.461   TIME= 07:41:07
 Contact related postprocess items (ETABLE, pressure ...) are not
 available.
 Contact detection at: nodal point (Dual shape function based)
 MPC will be built internally to handle bonded contact.
 Average contact surface length               0.19070
 Average contact pair depth                   0.17708
 Average target surface length                0.19375
 Default pinball region factor PINB           0.25000
 The resulting pinball region                 0.44271E-01
 Initial penetration/gap is excluded.
 Bonded contact (always) is defined.

 *** NOTE ***                            CP =       2.461   TIME= 07:41:07
 Max.  Initial penetration 1.025113101E-15 was detected between contact
 element 3502 and target element 3587.
 ****************************************


 *** NOTE ***                            CP =       2.461   TIME= 07:41:07
 Symmetric Deformable- deformable contact pair identified by real
 constant set 7 and contact element type 6 has been set up.  The
 companion pair has real constant set ID 6.  Both pairs should have the
 same behavior.
 MAPDL will deactivate the current pair and keep its companion pair,
 resulting in asymmetric contact.
 Pure pore fluid contact is activated.
 Contact algorithm: MPC based approach

 *** NOTE ***                            CP =       2.461   TIME= 07:41:07
 Contact related postprocess items (ETABLE, pressure ...) are not
 available.
 Contact detection at: nodal point (Dual shape function based)
 MPC will be built internally to handle bonded contact.
 Average contact surface length               0.19070
 Average contact pair depth                   0.20000
 Average target surface length                0.19375
 Default pinball region factor PINB           0.25000
 The resulting pinball region                 0.50000E-01
 Initial penetration/gap is excluded.
 Bonded contact (always) is defined.

 *** NOTE ***                            CP =       2.461   TIME= 07:41:07
 Max.  Initial penetration 1.419491659E-15 was detected between contact
 element 3525 and target element 3452.
 ****************************************




            **** CENTER OF MASS, MASS, AND MASS MOMENTS OF INERTIA ****

  CALCULATIONS ASSUME ELEMENT MASS AT ELEMENT CENTROID

  TOTAL MASS =   16.537

                           MOM. OF INERTIA         MOM. OF INERTIA
  CENTER OF MASS            ABOUT ORIGIN        ABOUT CENTER OF MASS

  XC =  0.75000          IXX =    3329.          IXX =    61.45
  YC =   2.5000          IYY =    3189.          IYY =    15.40
  ZC =   13.833          IZZ =    164.8          IZZ =    52.15
                         IXY =   -31.01          IXY =   0.2558E-12
                         IYZ =   -571.9          IYZ =   0.6594E-11
                         IZX =   -171.6          IZX =   0.2018E-11


  *** MASS SUMMARY BY ELEMENT TYPE ***

  TYPE      MASS
     1   6.43125
     2   6.43125
     3   3.67500

 Range of element maximum matrix coefficients in global coordinates
 Maximum = 0.195402284 at element 1153.
 Minimum = 0.18606905 at element 1433.

   *** ELEMENT MATRIX FORMULATION TIMES
     TYPE    NUMBER   ENAME      TOTAL CP  AVE CP

        1       720  FLUID220      0.060   0.000084
        2       720  FLUID220      0.062   0.000087
        3       400  FLUID220      0.035   0.000088
        4        80  CONTA174      0.010   0.000129
        5        80  TARGE170      0.000   0.000002
        6        80  CONTA174      0.011   0.000133
        7        80  TARGE170      0.000   0.000003
 Time at end of element matrix formulation CP = 2.76435089.

 SPARSE MATRIX DIRECT SOLVER.
  Number of equations =        9246,    Maximum wavefront =     81
  Memory allocated for solver              =    52.753 MB
  Memory required for in-core solution     =    50.547 MB
  Memory required for out-of-core solution =    22.878 MB

 *** NOTE ***                            CP =       2.897   TIME= 07:41:08
 The Sparse Matrix Solver is currently running in the in-core memory
 mode.  This memory mode uses the most amount of memory in order to
 avoid using the hard drive as much as possible, which most often
 results in the fastest solution time.  This mode is recommended if
 enough physical memory is present to accommodate all of the solver
 data.
 Sparse solver maximum pivot= 0.772343337 at node 1569 PRES.
 Sparse solver minimum pivot= 0.104583107 at node 646 PRES.
 Sparse solver minimum pivot in absolute value= 0.104583107 at node 646
 PRES.

   *** ELEMENT RESULT CALCULATION TIMES
     TYPE    NUMBER   ENAME      TOTAL CP  AVE CP

        1       720  FLUID220      0.052   0.000072
        2       720  FLUID220      0.053   0.000074
        3       400  FLUID220      0.029   0.000073
        4        80  CONTA174      0.002   0.000024
        6        80  CONTA174      0.002   0.000023

   *** NODAL LOAD CALCULATION TIMES
     TYPE    NUMBER   ENAME      TOTAL CP  AVE CP

        1       720  FLUID220      0.004   0.000005
        2       720  FLUID220      0.004   0.000005
        3       400  FLUID220      0.002   0.000005
        4        80  CONTA174      0.000   0.000003
        6        80  CONTA174      0.000   0.000003
 *** LOAD STEP     1   SUBSTEP     1  COMPLETED.  FREQUENCY=   2.00000
 *** LOAD STEP     1   SUBSTEP     2  COMPLETED.  FREQUENCY=   4.00000
 *** LOAD STEP     1   SUBSTEP     3  COMPLETED.  FREQUENCY=   6.00000
 *** LOAD STEP     1   SUBSTEP     4  COMPLETED.  FREQUENCY=   8.00000
 *** LOAD STEP     1   SUBSTEP     5  COMPLETED.  FREQUENCY=   10.0000
 *** LOAD STEP     1   SUBSTEP     6  COMPLETED.  FREQUENCY=   12.0000
 *** LOAD STEP     1   SUBSTEP     7  COMPLETED.  FREQUENCY=   14.0000
 *** LOAD STEP     1   SUBSTEP     8  COMPLETED.  FREQUENCY=   16.0000
 *** LOAD STEP     1   SUBSTEP     9  COMPLETED.  FREQUENCY=   18.0000
 *** LOAD STEP     1   SUBSTEP    10  COMPLETED.  FREQUENCY=   20.0000
 *** LOAD STEP     1   SUBSTEP    11  COMPLETED.  FREQUENCY=   22.0000
 *** LOAD STEP     1   SUBSTEP    12  COMPLETED.  FREQUENCY=   24.0000
 *** LOAD STEP     1   SUBSTEP    13  COMPLETED.  FREQUENCY=   26.0000
 *** LOAD STEP     1   SUBSTEP    14  COMPLETED.  FREQUENCY=   28.0000
 *** LOAD STEP     1   SUBSTEP    15  COMPLETED.  FREQUENCY=   30.0000
 *** LOAD STEP     1   SUBSTEP    16  COMPLETED.  FREQUENCY=   32.0000
 *** LOAD STEP     1   SUBSTEP    17  COMPLETED.  FREQUENCY=   34.0000
 *** LOAD STEP     1   SUBSTEP    18  COMPLETED.  FREQUENCY=   36.0000
 *** LOAD STEP     1   SUBSTEP    19  COMPLETED.  FREQUENCY=   38.0000
 *** LOAD STEP     1   SUBSTEP    20  COMPLETED.  FREQUENCY=   40.0000
 *** LOAD STEP     1   SUBSTEP    21  COMPLETED.  FREQUENCY=   42.0000
 *** LOAD STEP     1   SUBSTEP    22  COMPLETED.  FREQUENCY=   44.0000
 *** LOAD STEP     1   SUBSTEP    23  COMPLETED.  FREQUENCY=   46.0000
 *** LOAD STEP     1   SUBSTEP    24  COMPLETED.  FREQUENCY=   48.0000
 *** LOAD STEP     1   SUBSTEP    25  COMPLETED.  FREQUENCY=   50.0000
 *** LOAD STEP     1   SUBSTEP    26  COMPLETED.  FREQUENCY=   52.0000
 *** LOAD STEP     1   SUBSTEP    27  COMPLETED.  FREQUENCY=   54.0000
 *** LOAD STEP     1   SUBSTEP    28  COMPLETED.  FREQUENCY=   56.0000
 *** LOAD STEP     1   SUBSTEP    29  COMPLETED.  FREQUENCY=   58.0000
 *** LOAD STEP     1   SUBSTEP    30  COMPLETED.  FREQUENCY=   60.0000
 *** LOAD STEP     1   SUBSTEP    31  COMPLETED.  FREQUENCY=   62.0000
 *** LOAD STEP     1   SUBSTEP    32  COMPLETED.  FREQUENCY=   64.0000
 *** LOAD STEP     1   SUBSTEP    33  COMPLETED.  FREQUENCY=   66.0000
 *** LOAD STEP     1   SUBSTEP    34  COMPLETED.  FREQUENCY=   68.0000
 *** LOAD STEP     1   SUBSTEP    35  COMPLETED.  FREQUENCY=   70.0000
 *** LOAD STEP     1   SUBSTEP    36  COMPLETED.  FREQUENCY=   72.0000
 *** LOAD STEP     1   SUBSTEP    37  COMPLETED.  FREQUENCY=   74.0000
 *** LOAD STEP     1   SUBSTEP    38  COMPLETED.  FREQUENCY=   76.0000
 *** LOAD STEP     1   SUBSTEP    39  COMPLETED.  FREQUENCY=   78.0000
 *** LOAD STEP     1   SUBSTEP    40  COMPLETED.  FREQUENCY=   80.0000
 *** LOAD STEP     1   SUBSTEP    41  COMPLETED.  FREQUENCY=   82.0000
 *** LOAD STEP     1   SUBSTEP    42  COMPLETED.  FREQUENCY=   84.0000
 *** LOAD STEP     1   SUBSTEP    43  COMPLETED.  FREQUENCY=   86.0000
 *** LOAD STEP     1   SUBSTEP    44  COMPLETED.  FREQUENCY=   88.0000
 *** LOAD STEP     1   SUBSTEP    45  COMPLETED.  FREQUENCY=   90.0000
 *** LOAD STEP     1   SUBSTEP    46  COMPLETED.  FREQUENCY=   92.0000
 *** LOAD STEP     1   SUBSTEP    47  COMPLETED.  FREQUENCY=   94.0000
 *** LOAD STEP     1   SUBSTEP    48  COMPLETED.  FREQUENCY=   96.0000
 *** LOAD STEP     1   SUBSTEP    49  COMPLETED.  FREQUENCY=   98.0000
 *** LOAD STEP     1   SUBSTEP    50  COMPLETED.  FREQUENCY=   100.000


 *** MAPDL BINARY FILE STATISTICS
  BUFFER SIZE USED= 16384
        9.938 MB WRITTEN ON ELEMENT MATRIX FILE: file0.emat
        4.125 MB WRITTEN ON ELEMENT SAVED DATA FILE: file0.esav
        6.188 MB WRITTEN ON ASSEMBLED MATRIX FILE: file0.full
       20.875 MB WRITTEN ON RESULTS FILE: file0.rst
 *************** Write FE CONNECTORS *********

 WRITE OUT CONSTRAINT EQUATIONS TO FILE= file.ce

 FINISH SOLUTION PROCESSING


 ***** ROUTINE COMPLETED *****  CP =         9.813



 PRINTOUT RESUMED BY /GOP

 *GET  _WALLASOL  FROM  ACTI  ITEM=TIME WALL  VALUE=  7.68750000

 *** MAPDL - ENGINEERING ANALYSIS SYSTEM  RELEASE 2024 R1          24.1     ***
 Ansys Mechanical Enterprise
 00000000  VERSION=LINUX x64     07:41:15  MAY 06, 2024 CP=      9.816

 --Harmonic Acoustics



          ***** MAPDL RESULTS INTERPRETATION (POST1) *****

 *** NOTE ***                            CP =       9.816   TIME= 07:41:15
 Reading results into the database (SET command) will update the current
 displacement and force boundary conditions in the database with the
 values from the results file for that load set.  Note that any
 subsequent solutions will use these values unless action is taken to
 either SAVE the current values or not overwrite them (/EXIT,NOSAVE).

 Set Encoding of XML File to:ISO-8859-1

 Set Output of XML File to:
     PARM,     ,     ,     ,     ,     ,     ,     ,     ,     ,     ,     ,
         ,     ,     ,     ,     ,     ,     ,

 DATABASE WRITTEN ON FILE  parm.xml

 EXIT THE MAPDL POST1 DATABASE PROCESSOR


 ***** ROUTINE COMPLETED *****  CP =         9.818



 PRINTOUT RESUMED BY /GOP

 *GET  _WALLDONE  FROM  ACTI  ITEM=TIME WALL  VALUE=  7.68750000

 PARAMETER _PREPTIME =     0.000000000

 PARAMETER _SOLVTIME =     8.000000000

 PARAMETER _POSTTIME =     0.000000000

 PARAMETER _TOTALTIM =     8.000000000

 *GET  _DLBRATIO  FROM  ACTI  ITEM=SOLU DLBR  VALUE=  0.00000000

 *GET  _COMBTIME  FROM  ACTI  ITEM=SOLU COMB  VALUE= 0.121714840

 *GET  _SSMODE   FROM  ACTI  ITEM=SOLU SSMM  VALUE=  2.00000000

 *GET  _NDOFS    FROM  ACTI  ITEM=SOLU NDOF  VALUE=  9246.00000

 /FCLEAN COMMAND REMOVING ALL LOCAL FILES
 --- Total number of nodes = 9687
 --- Total number of elements = 2160
 --- Element load balance ratio = 0
 --- Time to combine distributed files = 0.12171484
 --- Sparse memory mode = 2
 --- Number of DOF = 9246

 EXIT MAPDL WITHOUT SAVING DATABASE


 NUMBER OF WARNING MESSAGES ENCOUNTERED=          4
 NUMBER OF ERROR   MESSAGES ENCOUNTERED=          0

+--------------------- M A P D L   S T A T I S T I C S ------------------------+

Release: 2024 R1            Build: 24.1       Update: UP20231106   Platform: LINUX x64
Date Run: 05/06/2024   Time: 07:41     Process ID: 6362
Operating System: Ubuntu 20.04.6 LTS

Processor Model: AMD EPYC 7763 64-Core Processor

Compiler: Intel(R) Fortran Compiler Classic Version 2021.9  (Build: 20230302)
          Intel(R) C/C++ Compiler Classic Version 2021.9  (Build: 20230302)
          Intel(R) Math Kernel Library Version 2020.0.0 Product Build 20191122
          BLAS Library supplied by AMD BLIS

Number of machines requested            :    1
Total number of cores available         :    8
Number of physical cores available      :    4
Number of processes requested           :    4
Number of threads per process requested :    1
Total number of cores requested         :    4 (Distributed Memory Parallel)
MPI Type: INTELMPI
MPI Version: Intel(R) MPI Library 2021.10 for Linux* OS


GPU Acceleration: Not Requested

Job Name: file0
Input File: dummy.dat

  Core                Machine Name   Working Directory
 -----------------------------------------------------
     0                f8cd179a5ae7   /github/home/.mw/Application Data/Ansys/v241/AnsysMech1B1F/Project_Mech_Files/HarmonicAcoustics
     1                f8cd179a5ae7   /github/home/.mw/Application Data/Ansys/v241/AnsysMech1B1F/Project_Mech_Files/HarmonicAcoustics
     2                f8cd179a5ae7   /github/home/.mw/Application Data/Ansys/v241/AnsysMech1B1F/Project_Mech_Files/HarmonicAcoustics
     3                f8cd179a5ae7   /github/home/.mw/Application Data/Ansys/v241/AnsysMech1B1F/Project_Mech_Files/HarmonicAcoustics

Latency time from master to core     1 =    2.216 microseconds
Latency time from master to core     2 =    2.179 microseconds
Latency time from master to core     3 =    2.159 microseconds

Communication speed from master to core     1 =  9134.04 MB/sec
Communication speed from master to core     2 = 12991.90 MB/sec
Communication speed from master to core     3 = 15433.93 MB/sec

Total CPU time for main thread                    :        8.3 seconds
Total CPU time summed for all threads             :       10.1 seconds

Elapsed time spent obtaining a license            :        0.3 seconds
Elapsed time spent pre-processing model (/PREP7)  :        0.0 seconds
Elapsed time spent solution - preprocessing       :        0.3 seconds
Elapsed time spent computing solution             :        7.5 seconds
Elapsed time spent solution - postprocessing      :        0.1 seconds
Elapsed time spent post-processing model (/POST1) :        0.0 seconds

Equation solver used                              :            Sparse (symmetric)
Equation solver computational rate                :        6.2 Gflops

Sum of disk space used on all processes           :      218.9 MB

Sum of memory used on all processes               :      612.0 MB
Sum of memory allocated on all processes          :     3246.0 MB
Physical memory available                         :         31 GB
Total amount of I/O written to disk               :        1.2 GB
Total amount of I/O read from disk                :        2.1 GB

+------------------ E N D   M A P D L   S T A T I S T I C S -------------------+


 *-----------------------------------------------------------------------------*
 |                                                                             |
 |                               RUN COMPLETED                                 |
 |                                                                             |
 |-----------------------------------------------------------------------------|
 |                                                                             |
 |  Ansys MAPDL 2024 R1         Build 24.1         UP20231106    LINUX x64     |
 |                                                                             |
 |-----------------------------------------------------------------------------|
 |                                                                             |
 |  Database Requested(-db)     1024 MB     Scratch Memory Requested   1024 MB |
 |  Max Database Used(Master)      5 MB     Max Scratch Used(Master)    148 MB |
 |  Max Database Used(Workers)     5 MB     Max Scratch Used(Workers)   148 MB |
 |  Sum Database Used(All)        20 MB     Sum Scratch Used(All)       592 MB |
 |                                                                             |
 |-----------------------------------------------------------------------------|
 |                                                                             |
 |        CP Time      (sec) =         10.084       Time  =  07:41:15          |
 |        Elapsed Time (sec) =         10.000       Date  =  05/06/2024        |
 |                                                                             |
 *-----------------------------------------------------------------------------*

Project tree#

def print_tree(node, indentation=""):
    if hasattr(node, "Suppressed") and node.Suppressed is True:
        print(f"{indentation}├── {node.Name} (Suppressed)")
    else:
        print(f"{indentation}├── {node.Name}")

    if (
        hasattr(node, "Children")
        and node.Children is not None
        and node.Children.Count > 0
    ):
        for child in node.Children:
            print_tree(child, indentation + "|  ")


root_node = DataModel.Project
print_tree(root_node)
├── Project
|  ├── Model
|  |  ├── Geometry Imports
|  |  |  ├── Geometry Import
|  |  ├── Geometry
|  |  |  ├── Solid (Suppressed)
|  |  |  |  ├── Solid (Suppressed)
|  |  |  ├── Solid (Suppressed)
|  |  |  |  ├── Solid (Suppressed)
|  |  |  ├── Solid (Suppressed)
|  |  |  |  ├── Solid (Suppressed)
|  |  |  ├── Solid (Suppressed)
|  |  |  |  ├── Solid (Suppressed)
|  |  |  ├── Solid (Suppressed)
|  |  |  |  ├── Solid (Suppressed)
|  |  |  ├── Solid
|  |  |  |  ├── Solid
|  |  |  ├── Solid (Suppressed)
|  |  |  |  ├── Solid (Suppressed)
|  |  |  ├── Solid
|  |  |  |  ├── Solid
|  |  |  ├── Solid
|  |  |  |  ├── Solid
|  |  |  ├── Solid (Suppressed)
|  |  |  |  ├── Solid (Suppressed)
|  |  |  ├── Solid (Suppressed)
|  |  |  |  ├── Solid (Suppressed)
|  |  ├── Materials
|  |  |  ├── Structural Steel
|  |  |  ├── Air
|  |  ├── Coordinate Systems
|  |  |  ├── Global Coordinate System
|  |  |  ├── Coordinate System
|  |  ├── Remote Points
|  |  ├── Connections
|  |  |  ├── Contacts
|  |  |  |  ├── Contact Region
|  |  |  |  ├── Contact Region 2
|  |  |  |  ├── Contact Region 3
|  |  |  |  ├── Contact Region 4
|  |  |  |  ├── Contact Region 5
|  |  |  |  ├── Contact Region 6
|  |  |  |  ├── Contact Region 7
|  |  |  |  ├── Contact Region 8
|  |  |  |  ├── Contact Region 9
|  |  |  |  ├── Contact Region 10
|  |  ├── Mesh Workflows
|  |  ├── Mesh
|  |  ├── Named Selections
|  |  |  ├── SF_Velo
|  |  |  ├── ABS_Face
|  |  |  ├── PRES_Face
|  |  |  ├── ACOUSTIC_Region
|  |  ├── Harmonic Acoustics
|  |  |  ├── Pre-Stress/Modal (None)
|  |  |  ├── Analysis Settings
|  |  |  ├── Acoustics Region
|  |  |  ├── Surface Velocity
|  |  |  ├── Pressure
|  |  |  ├── Absorption Surface
|  |  |  ├── Solution
|  |  |  |  ├── Solution Information
|  |  |  |  ├── Acoustic Pressure
|  |  |  |  ├── Acoustic Total Velocity
|  |  |  |  ├── Acoustic Directional Velocity
|  |  |  |  ├── Acoustic Directional Velocity 2
|  |  |  |  ├── Sound Pressure Level
|  |  |  |  ├── A-Weighted Sound Pressure Level
|  |  |  |  ├── Frequency Band SPL
|  |  |  |  ├── A-Weighted Frequency Band SPL
|  |  |  |  ├── Acoustic Total Velocity 2
|  |  |  |  ├── Acoustic Directional Velocity 3
|  |  |  |  ├── Kinetic Energy
|  |  |  |  ├── Potential Energy
|  |  |  |  ├── Frequency Response
|  |  |  |  ├── Frequency Response 2
|  |  |  |  ├── Frequency Response 3

Cleanup#

Save project

app.save(os.path.join(cwd, "harmnonic_acoustics.mechdat"))
app.new()

# delete example file
delete_downloads()
True

Total running time of the script: (0 minutes 27.669 seconds)

Gallery generated by Sphinx-Gallery