How to sum up HGENRATE NMISC 39 using Mechanical scripting?

Question

How to sum up HGENRATE NMISC 39 (SOLID278/Thermal) using Mechanical scripting?

Erik Kostson · Accepted Answer

The below script does that for all steps and elements:

APDL way:

SET,LAST,LAST ! set last result

ETABLE,mynmisc39, NMISC, 39 ! get nmisc39

SSUM ! sum

*GET,my_sum,SSUM,,ITEM,mynmisc39 ! get sum as my_sum

Mech. Script way

analysis = ExtAPI.DataModel.AnalysisList[0]

# Get result reader

reader = analysis.GetResultsData()

# Get global mesh data (contains element IDs)

mesh = analysis.MeshData     # or MeshDataByName("Global") if preferred

elemIds = mesh.ElementIds    # list/array of element IDs

nElems = len(elemIds)        # total element count

# Example: retrieve stress result "S" as an elemental result

datasets = reader.ListTimeFreq

#Loop over results

for i in range(len(datasets)):

    reader.CurrentTimeFreq = reader.ListTimeFreq[i] #Here is where time is selected

    rs = reader.GetResult("NMISC39")

    sValues = rs.GetElementValues(elemIds, False)

    totw=0

# Loop over elements with index controlled by nElems

    for e in range(nElems):

        eid = elemIds[e]

        totw=totw+sValues[e]

    print(str(totw) + " , " + str(i))