Added text to nodal results plot

Author: akaszynski

doc/conf.py

@@ -36,9 +36,6 @@
 templates_path = ['_templates']
 
 # The suffix(es) of source filenames.
-# You can specify multiple suffix as a list of string:
-#
-# source_suffix = ['.rst', '.md']
 source_suffix = '.rst'
 
 # The master toctree document.
@@ -54,9 +51,10 @@
 # built documents.
 #
 # The short X.Y version.
-version = u'0.16'
+version = u'0.17'
+
 # The full version, including alpha/beta/rc tags.
-release = u'0.16'
+release = u'0.17.1'
 
 # The language for content autogenerated by Sphinx. Refer to documentation
 # for a list of supported languages.

doc/loading_results.rst

@@ -29,14 +29,14 @@ accesses the results through a python object `result` which you can create with
     import pyansys
     result = pyansys.ResultReader('file.rst')
     
-Upon initialization this result file contains several properties to include the
-time values from the analysis, node numbering, element numbering, etc.
+Upon initialization the ``ResultReader`` object contains several properties to 
+include the time values from the analysis, node numbering, element numbering, etc.
 
 
-Result File Properties
-~~~~~~~~~~~~~~~~~~~~~~
+ResultReader Properties
+~~~~~~~~~~~~~~~~~~~~~~~
 
-The properties of the result file can be listed the below code.  At the moment,
+The properties of the ``ResultReader`` are listed below.  At the moment,
 the property listing is limited to only the number of results in the file.
 
 .. code:: python
@@ -63,7 +63,7 @@ Geometry
 ~~~~~~~~
 
 The geometry of the model can be accessed directly from the dictionary by 
-accessing
+accessing:
 
 .. code:: python
 
@@ -103,9 +103,9 @@ degrees of freedom.
     result.PlotNodalResult(0, label='Normalized')
 
 Stress can be obtained as well using the below code.  The nodal stress is 
-computed in the same manner that ANSYS uses by to determine the stress at each
-node by averaging the stress evaluated at that node for all attached elements.
-For now, only component stresses can be displayed.
+computed in the same manner as ANSYS by averaging the stress evaluated at
+that node for all attached elements.  For now, only component stresses can 
+be displayed.
 
 .. code:: python
     

pyansys/_version.py

@@ -1,5 +1,5 @@
 # major, minor, patch
-version_info = 0, 17, 1
+version_info = 0, 17, 2
 
 # Nice string for the version
 __version__ = '.'.join(map(str, version_info))

pyansys/binary_reader.py

@@ -180,6 +180,9 @@ def LoadKM(self, as_sparse=True, utri=True):
         m_diag = m_block[3]
         m_data_diag = m_block[4]
 
+        self.m_block = m_block
+        self.k_block = k_block
+        
         # assemble data
         if utri:
             # stiffness matrix
@@ -203,6 +206,13 @@ def LoadKM(self, as_sparse=True, utri=True):
             mcol = np.hstack((m_block[1], m_block[0], m_diag)) # col and diag
             mdata= np.hstack((m_block[2], m_block[2], m_data_diag)) # data and diag
 
+        # store data for later reference
+        self.krow  = krow
+        self.kcol  = kcol
+        self.kdata = kdata
+        self.mrow  = mrow
+        self.mcol  = mcol
+        self.mdata = mdata
 
         # number of dimentions
         ndim = nref.size
@@ -301,7 +311,7 @@ def __init__(self, filename, logger=False, load_geometry=True):
         if load_geometry:
             self.StoreGeometry()
 
-        if self.resultheader['nSector'] and load_geometry:
+        if self.resultheader['nSector'] > 1 and load_geometry:
             self.iscyclic = True
 
             # Add cyclic properties
@@ -608,11 +618,22 @@ def PlotNodalResult(self, rnum, comp='norm', as_abs=False, label=''):
             d = np.abs(d)
 
         # Generate plot
-        text = 'Result {:d} at {:f}'.format(rnum + 1, self.tvalues[rnum])
+#        text = 'Result {:d} at {:f}'.format(rnum + 1, self.tvalues[rnum])
+        
+        # setup text
+        ls_table = self.resultheader['ls_table']
+        freq = self.GetTimeValues()[rnum]
+        text  = 'Cumulative Index: {:3d}\n'.format(ls_table[rnum, 2])
+        text += 'Loadstep:         {:3d}\n'.format(ls_table[rnum, 0])
+        text += 'Substep:          {:3d}\n'.format(ls_table[rnum, 1])
+#        text += 'Harmonic Index:   {:3d}\n'.format(hindex)
+        text += 'Frequency:      {:10.4f} Hz'.format(freq)
+
+        
         plobj = vtkInterface.PlotClass()
         plobj.AddMesh(self.uGrid, no_copy=True, scalars=d, stitle=stitle, 
                       flipscalars=True, interpolatebeforemap=True)
-        plobj.AddText(text)
+        plobj.AddText(text, fontsize=20)
         cpos = plobj.Plot();  # store camera position
         del plobj