Adding temperature dependent properties (arrays).

ansys/mapdl/reader/rst.py

@@ -381,9 +381,32 @@ def _load_materials(self):
         def read_mat_data(ptr):
             """reads double precision material data given an offset from ptrMAT"""
             arr, sz = self.read_record(self._geometry_header["ptrMAT"] + ptr, True)
-            arr = arr[arr != 0]
-            if arr.size == 1:
-                return arr[0]
+            is_single_value = arr[arr != 0].size == 1
+
+            if is_single_value:
+                return arr[arr != 0][0]
+            else:
+                # We need to account for 0 at the temp and properties.
+                temps_ = arr[:sz]
+                values_ = arr[-1 : -1 * (sz + 1) : -1]
+
+                prop_ = np.vstack((temps_, values_))
+                for each in range(prop_.shape[1] - 1, -1, -1):
+                    # checking that two records are empty starting from the end.
+                    if (
+                        temps_[each] == 0
+                        and values_[each] == 0
+                        and temps_[each - 1] == 0
+                        and values_[each - 1] == 0
+                    ):
+                        continue  # they are zero, so we keep going.
+                    else:
+                        # found a non-zero report
+                        break
+
+                prop_ = prop_[:, :each]
+                prop_[1, :] = prop_[1, ::-1]
+                return prop_
 
         mat_table = self.read_record(self._geometry_header["ptrMAT"])
         if mat_table[0] != -101:  # pragma: no cover
@@ -423,9 +446,7 @@ def read_mat_data(ptr):
                         self._cfile._seekg(fs)
                         material[key] = np.fromstring(self._cfile._read(8))[0]
                     else:
-                        for i in range(-5, 10):
-                            material[key] = read_mat_data(ptr)
-                            # print(read_mat_data(ptr))
+                        material[key] = read_mat_data(ptr)
 
             # documentation for this is hard to follow, but it boils down to
             # the fact that "The record includes MP pointers followed by TB

tests/test_rst.py

@@ -636,3 +636,18 @@ def test_materials_v150():
         ans_mat = rst.materials[mat_type]
         for key, value in known_material.items():
             assert pytest.approx(ans_mat[key]) == known_material[key]
+
+
+def test_temperature_dependent_properties():
+    path_rth = os.path.join(testfiles_path, "temp_dependent_results.rth")
+    rst = pymapdl_reader.read_binary(path_rth)
+    mats = rst.materials
+    kxx = mats[1]["KXX"]
+    expected_value = np.array(
+        [[-100.0, 0.0, 100.0, 200.0], [114.0, 144.0, 165.0, 175.0]]
+    )
+
+    assert kxx is not None
+    assert isinstance(kxx, np.ndarray)
+    assert kxx.size > 0
+    assert np.allclose(kxx, expected_value)