How do I get real and imaginary stress tensor for a harmonic analysis using DPF-Core?

I am using DPF-Core (PyAnsys) via an IDE (Integrated Development Environment eg. Visual Studio Code) and using PyDPF to retrieve Data from an RST file of a Harmonic analysis. The RST file contains a named-selection "NODE" (a set of nodes). How do I retrieve the stress-tensor (real and imaginary part) for each node at each frequency?

Find more posts tagged with

PyDPF

Accepted answers

Ayush Kumar

from ansys.dpf import core as dpf
from ansys.dpf.core import operators as ops

model = dpf.Model(r"\Path\to\my\file.rst")
print(model)

tf = model.metadata.time_freq_support
print('Number of solution sets', tf.n_sets)

stress = ops.result.stress(data_sources=model)
freq_ids = list(range(1, model.metadata.time_freq_support.n_sets+1))

print(model.metadata.available_named_selections)

ns_operator = dpf.operators.scoping.on_named_selection()
ns_operator.inputs.named_selection_name.connect('NODE')
ns_operator.inputs.data_sources.connect(model)
stress.inputs.mesh_scoping.connect(ns_operator.outputs.mesh_scoping)
stress.inputs.time_scoping.connect(freq_ids)
stress.inputs.requested_location.connect("Nodal") 

fields = stress.outputs.fields_container.get_data()

for freq_id in freq_ids:
    print("REAL - %s" % fields.get_field_by_time_id(freq_id).data)
    print("IMAG - %s" % fields.get_imaginary_field(freq_id).data)

Ayush Kumar

In case you want to retrieve this data in Mechanical, mech_dpf must be used, the following code will give the same output in Mechanical Scripting Console.

NOTE: Please set Mechanical analysis system index in line number 6. This code assumes the first analysis to be Harmonic and retrieves data from that analysis.

import mech_dpf
mech_dpf.setExtAPI(ExtAPI)
import Ans.DataProcessing as dpf
dpf.DataProcessingCore.LoadLibrary("math",r"C:\Program Files\ANSYS Inc\v212\aisol\bin\winx64\Ans.Dpf.Math.dll")

ds = mech_dpf.GetDataSources(0)
timeProvider = dpf.operators.metadata.time_freq_provider()
timeProvider.inputs.data_sources.Connect(ds)
freqs = timeProvider.outputs.time_freq_support.GetData()
numSets = freqs.NumberSets
ids = [index + 1 for index in range(numSets)]

freq_scoping = dpf.Scoping()
freq_scoping.Ids = ids

ns_op = dpf.Operator("scoping_provider_by_ns")
ns_op.Connect(1, "NODE")
ns_op.Connect(4, ds)

#Stress Y
op_sy = dpf.Operator("S")
op_sy.Connect(0, freq_scoping)
op_sy.Connect(1, ns_op)
op_sy.Connect(4, ds)
op_sy.Connect(9, "Nodal")

for id in ids:
    print "REAL - %s " % op_sy.GetOutputAsFieldsContainer(0).GetFieldByTimeId(id).Data
    print "IMAG - %s " % op_sy.GetOutputAsFieldsContainer(0).GetImaginaryField(id).Data

All comments

Ayush Kumar

from ansys.dpf import core as dpf
from ansys.dpf.core import operators as ops

model = dpf.Model(r"\Path\to\my\file.rst")
print(model)

tf = model.metadata.time_freq_support
print('Number of solution sets', tf.n_sets)

stress = ops.result.stress(data_sources=model)
freq_ids = list(range(1, model.metadata.time_freq_support.n_sets+1))

print(model.metadata.available_named_selections)

ns_operator = dpf.operators.scoping.on_named_selection()
ns_operator.inputs.named_selection_name.connect('NODE')
ns_operator.inputs.data_sources.connect(model)
stress.inputs.mesh_scoping.connect(ns_operator.outputs.mesh_scoping)
stress.inputs.time_scoping.connect(freq_ids)
stress.inputs.requested_location.connect("Nodal") 

fields = stress.outputs.fields_container.get_data()

for freq_id in freq_ids:
    print("REAL - %s" % fields.get_field_by_time_id(freq_id).data)
    print("IMAG - %s" % fields.get_imaginary_field(freq_id).data)

Ayush Kumar

In case you want to retrieve this data in Mechanical, mech_dpf must be used, the following code will give the same output in Mechanical Scripting Console.

NOTE: Please set Mechanical analysis system index in line number 6. This code assumes the first analysis to be Harmonic and retrieves data from that analysis.

import mech_dpf
mech_dpf.setExtAPI(ExtAPI)
import Ans.DataProcessing as dpf
dpf.DataProcessingCore.LoadLibrary("math",r"C:\Program Files\ANSYS Inc\v212\aisol\bin\winx64\Ans.Dpf.Math.dll")

ds = mech_dpf.GetDataSources(0)
timeProvider = dpf.operators.metadata.time_freq_provider()
timeProvider.inputs.data_sources.Connect(ds)
freqs = timeProvider.outputs.time_freq_support.GetData()
numSets = freqs.NumberSets
ids = [index + 1 for index in range(numSets)]

freq_scoping = dpf.Scoping()
freq_scoping.Ids = ids

ns_op = dpf.Operator("scoping_provider_by_ns")
ns_op.Connect(1, "NODE")
ns_op.Connect(4, ds)

#Stress Y
op_sy = dpf.Operator("S")
op_sy.Connect(0, freq_scoping)
op_sy.Connect(1, ns_op)
op_sy.Connect(4, ds)
op_sy.Connect(9, "Nodal")

for id in ids:
    print "REAL - %s " % op_sy.GetOutputAsFieldsContainer(0).GetFieldByTimeId(id).Data
    print "IMAG - %s " % op_sy.GetOutputAsFieldsContainer(0).GetImaginaryField(id).Data