Issues with Tabular Data Input Displacement

[Lucansys]

Good morning!

I am trying to automate my Ansys workflow integrating the PyAnsys Geometry with parametric dimension values and PyMechanical with different input displacement scenarios that are tabular data from bearing tested values.

I tried putting these values inside Globals as .txt files:

And tried to read them like this:

But every time I initialize and try to run this:

This happens:

The values from the tabular data are not inserted inside Ansys, what should I do?

[dipinknair]

@Lucansys Could you attach your entire script with necessary geometry files here so I can use them to debug ?

[Lucansys]

Of couse! Sorry for the late response.

########## 00_Global.py ##########
geometryPath_P = r"C:\Users\Lucas\Desktop\Projetos-PYTHON\pymech\INICIAÇÃO CIENTÍFICA\Dados\Geometrias\PyAnsysGeoExport\Design1.scdocx"
save_file = r'C:\Users\Lucas\Desktop\Projetos-PYTHON\pymech\INICIAÇÃO CIENTÍFICA\Dados\Salvos\Teste1.mechdb'
result_dir = r"C:\Users\Lucas\Desktop\Projetos-PYTHON\pymech\INICIAÇÃO CIENTÍFICA\Resultados"
freq_saudavel_str = r"C:\Users\Lucas\Desktop\Projetos-PYTHON\pymech\INICIAÇÃO CIENTÍFICA\Dados\Input Displacement\freq_saudavel.txt"
amp_saudavel_str = r"C:\Users\Lucas\Desktop\Projetos-PYTHON\pymech\INICIAÇÃO CIENTÍFICA\Dados\Input Displacement\amp_saudavel.txt"
fase_saudavel_str = r"C:\Users\Lucas\Desktop\Projetos-PYTHON\pymech\INICIAÇÃO CIENTÍFICA\Dados\Input Displacement\fase_saudavel.txt"


########## 01_AccessPoints.py ##########
proj = ExtAPI.DataModel.Project
model = proj.Model
model.AddHarmonicResponseAnalysis()
geo = model.Geometry
mesh = model.Mesh
connection = model.Connections
materials = model.Materials
analysis = model.Analyses[0]
solution = analysis.Solution



########## 02_Unidades.py ##########
ExtAPI.Application.ActiveUnitSystem = MechanicalUnitSystem.StandardMKS
ExtAPI.Application.ActiveAngleUnit = AngleUnitType.Radian



########## 03_GeometryImport.py ##########
geomImport = Model.GeometryImportGroup.AddGeometryImport()
# Format of Geometry file
geomImport_format = Ansys.Mechanical.DataModel.Enums.GeometryImportPreference.Format.Automatic

geomImport_preferences = Ansys.ACT.Mechanical.Utilities.GeometryImportPreferences()
geomImport_preferences.ProcessSolids = True
geomImport_preferences.ProcessSurfaces = False
geomImport_preferences.ProcessLines = False

geomImport_preferences.ProcessNamedSelections = True
geomImport_preferences.NamedSelectionKey = "NSel"

geomImport_preferences.ProcessMaterialProperties = True
geomImport_preferences.MixedImportResolution = GeometryImportPreference.MixedImportResolution.Solid

geomImport_preferences.ProcessCoordinateSystems = True
geomImport_preferences.CoordinateSystemKey = "CSys"
geomImport.Import(geometryPath_P, geomImport_format, geomImport_preferences)



########## 04_Mesh.py ##########
mesh = Model.Mesh
mesh.ElementSize = Quantity("0.001 [m]")
mesh.GenerateMesh()



########## 05_NamedSelections.py ##########
selinfo = ExtAPI.SelectionManager.CreateSelectionInfo(SelectionTypeEnum.GeometryEntities)
ns_faceDisplacement = Model.AddNamedSelection()
ns_faceDisplacement.Name = "ns_faceDisplacement"
selinfo.Ids = [10]
ns_faceDisplacement.Location = selinfo
ns_faceFixedSupport = Model.AddNamedSelection()
ns_faceFixedSupport.Name = "ns_faceFixedSupport"
selinfo.Ids = [11, 9]
ns_faceFixedSupport.Location = selinfo
ns_faceOutput = Model.AddNamedSelection()
ns_faceOutput.Name = "ns_faceOutput"
selinfo.Ids = [8]
ns_faceOutput.Location = selinfo



########## 06_Config.py ##########

analysis_settings = analysis.AnalysisSettings

analysis_settings.RangeMaximum = Quantity('10000 [Hz]')
analysis_settings.SolutionIntervals = 100 

analysis_settings.SolutionMethod = Ansys.Mechanical.DataModel.Enums.HarmonicMethod.ProgramControlled



########## 07_CondContorno.py ##########
fixedsupport = analysis.AddDisplacement()
fixedsupport.Location = ns_faceFixedSupport # Displacement em Y e Z = 0

fixedsupport.YComponent.Inputs[0].DiscreteValues = [Quantity('0 [sec]')]
fixedsupport.YComponent.Output.DiscreteValues = [Quantity('0 [m]')]
fixedsupport.ZComponent.Inputs[0].DiscreteValues = [Quantity('0 [Hz]')]
fixedsupport.ZComponent.Output.DiscreteValues = [Quantity('0 [m]')]

displacement = analysis.AddDisplacement()
displacement.Location = ns_faceDisplacement

displacement.XComponent.Inputs[0].DiscreteValues = freq_saudavel_str
displacement.YComponent.Inputs[0].DiscreteValues = [Quantity('0 [Hz]')] # Em Y = 0

displacement.XComponent.Output.DiscreteValues = amp_saudavel_str
displacement.YComponent.Output.DiscreteValues = [Quantity('0 [m]')] # Em Y = 0

displacement.XPhaseAngle.Output.DiscreteValues = fase_saudavel_str
displacement.YPhaseAngle.Output.DiscreteValues = [Quantity('0 [rad]')] # Em Y = 0



########## 08_Solve.py ##########
solution.Solve()
solution.GetResults()



########## 09_Post.py ##########
dirDef_Saida = solution.AddDirectionalDeformation()
dirDef_Saida.Location = ns_faceOutput

freqResp_Saida = solution.AddDeformationFrequencyResponse()
freqResp_Saida.Location = ns_faceOutput
freqResp_Saida.NormalOrientation = NormalOrientationType.XAxis

freqResp_Ent = solution.AddDeformationFrequencyResponse()
freqResp_Ent.Location = ns_faceDisplacement
freqResp_Ent.NormalOrientation = NormalOrientationType.XAxis

solution.EvaluateAllResults()

proj.Save(save_file)

The files are on this drive link:
https://drive.google.com/drive/folders/1u-8gMXezglSwOyCj3S8XzZ6AlfhCisC0?usp=sharing

[mikerife]

Hi @Lucansys
You can't do this in Mechanical using the UI (pointing directly to a external txt file) so scripting will not work either. Instead use the Python file open and read the file. Then apply the value. Something like this:

f = open(r"C:\Users\Lucas\Desktop\Projetos-PYTHON\pymech\INICIAÇÃO CIENTÍFICA\Dados\Input Displacement\freq_saudavel.txt")
freq_saudavel_str = f.read()
f.close()

displacement.XComponent.Inputs[0].DiscreteValues = [{freq_saudavel_str}]

Mike

[Lucansys]

I will try! Thank you for the response!

[Lucansys]

I tried something like this and didn't work. Am I doing something wrong?

[dipinknair]

@Lucansys Do you see the same issue without PyMechanical ? Just launching Mechanical and running your script inside scripting window ?

[Lucansys]

This is happening.

[Lucansys]

@dipinknair @mikerife Thank you so much for your help, but I managed to do what I was looking in a different way and it worked!

Now I'm having another problem on exporting the Directional Deformation values as txt or csv and the Frequency Responses Tabular Data as csv.

[dipinknair]

@Lucansys The right place to ask question related to Mechanical scripting is https://discuss.ansys.com/categories/structures.