Merge pull request #106 from yangbai90/master

add new damask notebook for different analysis
pyiron · Mar 8, 2022 · bbb66ce · bbb66ce
2 parents 0d5ff80 + 8109e41
commit bbb66ce
Show file tree

Hide file tree

Showing 5 changed files with 1,043 additions and 2 deletions.
diff --git a/notebooks/damask_tutorial_elastoplasticity-isotropic.ipynb b/notebooks/damask_tutorial_elastoplasticity-isotropic.ipynb
diff --git a/notebooks/damask_tutorial_elastoplasticity-powerlaw.ipynb b/notebooks/damask_tutorial_elastoplasticity-powerlaw.ipynb
diff --git a/notebooks/damask_tutorial_linearelastic.ipynb b/notebooks/damask_tutorial_linearelastic.ipynb
diff --git a/pyiron_continuum/damask/damaskjob.py b/pyiron_continuum/damask/damaskjob.py
@@ -110,6 +110,16 @@ def _load_results(self, file_name="damask_loading.hdf5"):
         self.output.strain = self.average_spatio_temporal_tensors('epsilon_V^0.0(F)')
         self.output.stress_von_Mises = self.average_spatio_temporal_tensors('sigma_vM')
         self.output.strain_von_Mises = self.average_spatio_temporal_tensors('epsilon_V^0.0(F)_vM')
+
+    def writeresults2vtk(self):
+        """
+            save results to vtk files
+        """
+        cwd = os.getcwd() # get the current dir
+        os.chdir(self.working_directory) # cd into the working dir
+        result=self._results
+        result.export_VTK()
+        os.chdir(cwd) # cd back to the notebook dir
 
     def temporal_spatial_shape(self, name):
         property_dict = self._results.get(name)

diff --git a/pyiron_continuum/damask/factory.py b/pyiron_continuum/damask/factory.py
@@ -173,7 +173,7 @@ def homogenization(method, parameters):
         return {method: parameters}
 
     @staticmethod
-    def phase(composition, lattice, output_list, elasticity, plasticity):
+    def phase(composition, lattice, output_list, elasticity, plasticity=None):
         """
         Returns a dictionary describing the phases for damask.
         Args:
@@ -196,7 +196,12 @@ def phase(composition, lattice, output_list, elasticity, plasticity):
                                     type='phenopowerlaw', xi_0_sl=[31e6],
                                     xi_inf_sl=[63e6])
         """
-        return {composition: {'lattice': lattice,
+        if plasticity == None:
+            return {composition: {'lattice': lattice,
+                              'mechanical': {'output': output_list,
+                                             'elastic': elasticity}}}
+        else:
+            return {composition: {'lattice': lattice,
                               'mechanical': {'output': output_list,
                                              'elastic': elasticity,
                                              'plasticity': plasticity}}}