Adapt .dat, fix INFO and EXPORT mode

dragos-ana · dragos-ana · commit 0481a4a0af69 · 2024-11-30T21:29:43.000+01:00
- The reading and writing of .dat files was adapted to the new 4C .dat
file structure
- Added new default .dat file from the current 4C repo
- Fixed issue regarding INFO_MODE and EXPORT_MODE, which were not
displayed after clicking on the respective buttons
- Added EXPORT_STATUS to signal to the user that the export was
successful (or not!)
diff --git a/src/fourc_webviewer/gui_utils.py b/src/fourc_webviewer/gui_utils.py
@@ -62,11 +62,29 @@ def _bottom_sheet_info():
 
 def _bottom_sheet_export(server_controller):
     with vuetify.VBottomSheet(v_model=("EXPORT_MODE",), inset=True):
-        with vuetify.VCard(classes="text-center", height=150, title="Export"):
+        with vuetify.VCard(classes="text-center", height=250, title="Export"):
             with vuetify.VCardText():
                 vuetify.VTextField(
                     label="Export .dat file", v_model=("EXPORT_DAT_PATH",)
                 )
+                vuetify.VAlert(
+                   title="Click on <save> to export the changed .dat file under the entered path",
+                   type="info", 
+                   v_if  = ("EXPORT_STATUS == EXPORT_STATUS_POSSIB[0]", ),
+                   classes = "h-50"
+                )
+                vuetify.VAlert(
+                   title="Your file was exported correctly!",
+                   type="success",
+                   v_if  = ("EXPORT_STATUS == EXPORT_STATUS_POSSIB[1]", ),
+                   classes = "h-50"
+                )
+                vuetify.VAlert(
+                   title="There was a problem while trying to export! Check the entered path and the modified .dat file settings!",
+                   type="error",
+                   v_if  = ("EXPORT_STATUS == EXPORT_STATUS_POSSIB[2]", ),
+                   classes = "h-50"
+                )
         vuetify.VBtn(text="SAVE", color="primary", click=server_controller)
 
 
diff --git a/src/fourc_webviewer/input_file_utils/read_dat_file.py b/src/fourc_webviewer/input_file_utils/read_dat_file.py
@@ -306,19 +306,6 @@ def read_dat_file(dat_file_path):
     )  # category variable -> used further down below for the sectioning of the file content
     while l < len(file_lines):
         match file_lines[l][0]:  # sectioning based on first character
-            case "=":  # get file title
-                curr_line_ind = l + 1
-                while curr_line_ind < len(file_lines):
-                    if file_lines[curr_line_ind][0] == "=":
-                        l = curr_line_ind + 1
-                        break
-                    else:
-                        file_title.append(file_lines[curr_line_ind])
-                        curr_line_ind += 1
-                        if curr_line_ind == len(file_lines):
-                            raise Exception(
-                                "The entire file was read without sectioning a file title"
-                            )
             case (
                 "-"
             ):  # Symbol indicates a new category -> The use of "-" as a bullet point in the file description is discussed below
@@ -520,6 +507,9 @@ def read_dat_file(dat_file_path):
     # take all lines from the geometry_categories and save them in a geometry_lines list
     geometry_lines = [l for l in file_lines[geometry_categories_indices[0] :]]
 
+    # add file name as a file title
+    file_title.append(Path(dat_file_path).name)
+
     return_dict = {
         "file_title": file_title,
         "file_description": file_description,
diff --git a/src/fourc_webviewer/input_file_utils/write_dat_file.py b/src/fourc_webviewer/input_file_utils/write_dat_file.py
@@ -38,23 +38,14 @@ def write_dat_file(
     #       geometry_lines: [<lines with geometry specific info>]; Example: ["-----------------------------------------------DNODE-NODE TOPOLOGY", "NODE    36 DNODE 1", "NODE    163 DNODE 1", ..., "-----------------------------------------------DLINE-NODE TOPOLOGY", "NODE    36 DLINE 1", "NODE    105 DLINE 1",...]
 
     # initialize lines_list, containing all the file lines
-    lines_list = ["//", "//"]
+    lines_list = []
 
     # number of characters for "===" and "---" lines
     num_of_chars = 100
 
     # line length for the category lines with <name     value>
     line_length = 60
 
-    # append title section
-    lines_list.append("=" * num_of_chars)
-    lines_list.append(
-        " " * math.floor((num_of_chars - len(title)) / 2)
-        + title
-        + " " * math.ceil((num_of_chars - len(title)) / 2)
-    )
-    lines_list.append("=" * num_of_chars)
-
     # append description section
     lines_list.append("-" * (num_of_chars - len("TITLE")) + "TITLE")
     lines_list.extend(description.split("\n"))
diff --git a/src/fourc_webviewer/run_webserver.py b/src/fourc_webviewer/run_webserver.py
@@ -25,6 +25,7 @@
 from vtkmodules.vtkCommonDataModel import vtkDataObject
 import fourc_webviewer.vturender as vtu
 from fourc_webviewer_default_files import DEFAULT_INPUT_FILE
+import time
 
 # ------------------------------------------------------------------------------#
 #                              COMMON SERVER SETUP                              #
@@ -431,6 +432,15 @@ def STATE_initialization(temp_dir, dat_path, vtu_path, dat_file_content):
     STATE.EDIT_MODE_POSSIB = ["VIEW MODE", "EDIT MODE"]
     STATE.EDIT_MODE = STATE.EDIT_MODE_POSSIB[0]
 
+    # initialize info mode value: False (bottom sheet with infos is not displayed until "INFO" button is pressed, and INFO_MODE is then set to True)
+    STATE.INFO_MODE = False
+
+    # initialize export mode value: False (bottom sheet with export settings is not displayed until "EXPORT" button is pressed, and EXPORT_MODE is then set to True)
+    STATE.EXPORT_MODE = False
+
+    # initialize the export status and its possible choices
+    STATE.EXPORT_STATUS_POSSIB = ["INFO", "SUCCESS", "ERROR"] # INFO: button was not yet clicked, SUCCESS: export was successful, ERROR: there was an error after trying to export
+    STATE.EXPORT_STATUS = "INFO"
 
 # ------------------------------------------------------------------------------#
 #                               STATE CHANGES                                  #
@@ -450,6 +460,12 @@ def change_dat_path(DAT_PATH, **kwargs):
     STATE.RENDER_COUNT = 1
 
 
+# after changing the export .dat file path
+@STATE.change("EXPORT_DAT_PATH")
+def change_export_dat_path(EXPORT_DAT_PATH, **kwargs):
+    # set export status to neutral ("INFO")
+    STATE.EXPORT_STATUS = STATE.EXPORT_STATUS_POSSIB[0]
+
 # after selecting another main category
 @STATE.change("SELECTED_MAIN_CATEGORY")
 def change_selected_main_category(SELECTED_MAIN_CATEGORY, **kwargs):
@@ -650,6 +666,12 @@ def change_materials_modif_attr(MATERIALS_MODIF_ATTR, **kwargs):
     STATE.dirty("MATERIALS")
 
 
+# after modifying the export mode (either after clicking on EXPORT or after canceling the export bottom sheet)
+@STATE.change("EXPORT_MODE")
+def change_export_mode(EXPORT_MODE, **kwargs):
+    STATE.EXPORT_STATUS = STATE.EXPORT_STATUS_POSSIB[0] # set export status to neutral ("INFO")
+
+
 # ------------------------------------------------------------------------------#
 #                            CLICK EVENT FUNCTIONS                             #
 # ------------------------------------------------------------------------------#
@@ -662,7 +684,6 @@ def click_info_button():
 def click_export_button():
     STATE.EXPORT_MODE = not STATE.EXPORT_MODE
 
-
 # convert provided dat file
 def click_convert_button():
     # global dat and vtu file paths
@@ -691,19 +712,23 @@ def click_convert_button():
 
 # export dat file
 def click_save_button():
-    write_dat_file(
-        STATE.TITLE,
-        STATE.DESCRIPTION,
-        STATE.CATEGORIES,
-        STATE.CATEGORY_ITEMS,
-        STATE.MATERIALS,
-        STATE.CLONING_MATERIAL_MAP,
-        STATE.FUNCT,
-        STATE.COND_GENERAL_TYPES,
-        STATE.COND_ENTITY_LIST,
-        STATE.COND_CONTEXT_LIST,
-        STATE.COND_TYPE_LIST,
-        STATE.RESULT_DESCRIPTION,
-        STATE.GEOMETRY_LINES,
-        STATE.EXPORT_DAT_PATH,
-    )
+    try:
+        write_dat_file(
+            STATE.TITLE,
+            STATE.DESCRIPTION,
+            STATE.CATEGORIES,
+            STATE.CATEGORY_ITEMS,
+            STATE.MATERIALS,
+            STATE.CLONING_MATERIAL_MAP,
+            STATE.FUNCT,
+            STATE.COND_GENERAL_TYPES,
+            STATE.COND_ENTITY_LIST,
+            STATE.COND_CONTEXT_LIST,
+            STATE.COND_TYPE_LIST,
+            STATE.RESULT_DESCRIPTION,
+            STATE.GEOMETRY_LINES,
+            STATE.EXPORT_DAT_PATH,
+        )
+        STATE.EXPORT_STATUS = STATE.EXPORT_STATUS_POSSIB[1]
+    except:
+        STATE.EXPORT_STATUS = STATE.EXPORT_STATUS_POSSIB[2]
diff --git a/src/fourc_webviewer_default_files/default_file.dat b/src/fourc_webviewer_default_files/default_file.dat
@@ -1,73 +1,104 @@
-==================================================================
-                  General Data File BACI
-==================================================================
 -------------------------------------------------------------TITLE
-Test: Apply constant displacement rate up to a certain displacement and hold this for a number of seconds.
-      It should be tested whether the implemented model MAT_InelasticDefgradUnifiedViscoplastic is able to 
-      model stress relaxation. 
+Simple test case based on 1hex8 element with isotropic elastic-viscoplastic
+behavior in a finite strain setting. The following conditions are assumed:
+- element is stretched at a constant strain rate
+- element can freely contract in lateral dimensions
+- viscoplasticity law: flow rule and hardening behavior modeled with the
+Reformulated Johnson-Cook law, as shown in:
+   Mareau: A thermodynamically consistent formulation of the Johnson-Cook
+           model, Mechanics of Materials 143, 2020
+- time integration of internal variables using logarithmic substepping, as shown in:
+   Ana: Continuum Modeling and Calibration of Viscoplasticity in the Context
+	of the Lithium Anode in Solid-State Batteries
+	(Master's Thesis, Chair of Computational Mechanics, TU Munich, 2024)
 ------------------------------------------------------PROBLEM SIZE
-//ELEMENTS    1
-//NODES       8
 DIM           3
---------------------------------------------------------PROBLEM TYP
-PROBLEMTYP                      Structure
------------------------------------------------------------------IO
+-------------------------------------------------------PROBLEM TYPE
+PROBLEMTYPE                      Structure
+----------------------------------------------------------------IO
 STRUCT_STRESS                   cauchy
 STRUCT_STRAIN                   log
-------------------------------------IO/RUNTIME VTK OUTPUT/STRUCTURE
+-----------------------------------IO/RUNTIME VTK OUTPUT/STRUCTURE
 OUTPUT_STRUCTURE                yes
 DISPLACEMENT                    yes
 STRESS_STRAIN                   yes
-----------------------------------------------IO/RUNTIME VTK OUTPUT
+---------------------------------------------IO/RUNTIME VTK OUTPUT
 OUTPUT_DATA_FORMAT              binary
 INTERVAL_STEPS                  1
--------------------------------------------------STRUCTURAL DYNAMIC
+------------------------------------------------STRUCTURAL DYNAMIC
 INT_STRATEGY                    Standard
-DYNAMICTYP                      Statics
+DYNAMICTYPE                      Statics
 RESULTSEVRY                     1
 RESTARTEVRY                     30
-TIMESTEP                        1.25e-1
-NUMSTEP                         1000000
-MAXTIME                         1.25e1
+TIMESTEP                        1.00e-6
+NUMSTEP                         100000000
+MAXTIME                         1.00e-4
 LINEAR_SOLVER                   1
 TOLRES                          1e-8
-TOLDISP                         1e-10
------------------------------------------------------------SOLVER 1
+TOLDISP                         1e-8
+NEGLECTINERTIA			           Yes
+----------------------------------------------------------SOLVER 1
 SOLVER                          UMFPACK
-----------------------------------------------------------MATERIALS
-MAT 1 MAT_MultiplicativeSplitDefgradElastHyper NUMMATEL 1 MATIDSEL 2 NUMFACINEL 1 INELDEFGRADFACIDS 3 DENS 5.34e-4
-MAT 2 ELAST_CoupNeoHooke YOUNG 7.810e3 NUE 0.38
-MAT 3 MAT_InelasticDefgradGenericIsotropicViscoplastic VP_FR_MATID 4
-MAT 4 MAT_ViscoplasticFlowRuleReformulatedJohnsonCook SIM_TEMPERATURE 293.0 MELT_TEMPERATURE 1793.0 REF_TEMPERATURE 293.0 TEMPERATURE_EXP 1.03 STRAIN_RATE_PREFAC 0.02 STRAIN_RATE_EXP_FAC 5 INIT_YIELD_STRENGTH 1.0 ISOTROP_HARDEN_PREFAC 10.0 ISOTROP_HARDEN_EXP 2.0
--------------------------------------------------------------FUNCT1
+---------------------------------------------------------MATERIALS
+MAT 1 MAT_MultiplicativeSplitDefgradElastHyper NUMMATEL 1 MATIDSEL 2 NUMFACINEL 1 INELDEFGRADFACIDS 3 DENS 7830.0
+MAT 2 ELAST_CoupNeoHooke YOUNG 200.0e3 NUE 0.29
+MAT 3 MAT_InelasticDefgradTransvIsotropElastViscoplast VISCOPLAST_LAW_ID 4 FIBER_READER_ID 5 YIELD_COND_A 0.0 YIELD_COND_B 0.0 YIELD_COND_F 0.0 ANISOTROPY isotrop LOG_SUBSTEP True MAX_HALVE_NUM_SUBSTEP 10
+MAT 4 MAT_ViscoplasticLawReformulatedJohnsonCook STRAIN_RATE_PREFAC 1.0e0 STRAIN_RATE_EXP_FAC 0.014 INIT_YIELD_STRENGTH 792.0e-0 ISOTROP_HARDEN_PREFAC 510.0e-0 ISOTROP_HARDEN_EXP 0.26
+MAT 5 ELAST_CoupTransverselyIsotropic ALPHA 1.0 BETA 1.0 GAMMA 1.0 ANGLE 0 STR_TENS_ID 100
+MAT 100 ELAST_StructuralTensor STRATEGY Standard
+------------------------------------------------------------FUNCT1
 COMPONENT 0 SYMBOLIC_FUNCTION_OF_SPACE_TIME 1.2*x
 COMPONENT 1 SYMBOLIC_FUNCTION_OF_SPACE_TIME -0.080830788307857
 ------------------------------------------------RESULT DESCRIPTION
-STRUCTURE DIS structure NODE 4 QUANTITY dispy VALUE -1.31962249833408790e-01 TOLERANCE 9.1e-10
-STRUCTURE DIS structure NODE 5 QUANTITY dispx VALUE -1.31962249833408513e-01 TOLERANCE 9.1e-10
-STRUCTURE DIS structure NODE 8 QUANTITY dispx VALUE -9.06271329067036280e-02 TOLERANCE 9.1e-10
-STRUCTURE DIS structure NODE 8 QUANTITY dispy VALUE -9.06271329067068060e-02 TOLERANCE 9.1e-10
-STRUCTURE DIS structure NODE 8 QUANTITY dispz VALUE 2.84025416687741394e-01 TOLERANCE 2.8e-09
-STRUCTURE DIS structure NODE 6 QUANTITY dispx VALUE -9.06271329067021153e-02 TOLERANCE 9.1e-10
-STRUCTURE DIS structure NODE 4 QUANTITY stress_zz VALUE  1.83879098311578976e-01 TOLERANCE 1.8e-07
-STRUCTURE DIS structure NODE 8 QUANTITY stress_zz VALUE  1.83879098311589884e-01 TOLERANCE 1.8e-07
-------------------------------------------------DESIGN DESCRIPTION
-NDPOINT 1
-NDLINE  0
-NDSURF  2
-NDVOL   0
+STRUCTURE DIS structure NODE 4 QUANTITY dispy VALUE -3.96619353942104843e-03 TOLERANCE 4.0e-11
+STRUCTURE DIS structure NODE 5 QUANTITY dispx VALUE -3.96619353942213090e-03 TOLERANCE 4.0e-11
+STRUCTURE DIS structure NODE 8 QUANTITY dispx VALUE -3.96619353942162696e-03 TOLERANCE 4.0e-11
+STRUCTURE DIS structure NODE 8 QUANTITY dispy VALUE -3.96619353942309367e-03 TOLERANCE 4.0e-11
+STRUCTURE DIS structure NODE 8 QUANTITY dispz VALUE  1.00501670841679491e-02 TOLERANCE 1.0e-10
+STRUCTURE DIS structure NODE 6 QUANTITY dispx VALUE -3.96619353942174319e-03 TOLERANCE 4.0e-11
+STRUCTURE DIS structure NODE 4 QUANTITY stress_zz VALUE  9.78458600965294636e+02 TOLERANCE 9.8e-04
+STRUCTURE DIS structure NODE 8 QUANTITY stress_zz VALUE  9.78458600965331243e+02 TOLERANCE 9.8e-04
 ------------------------------------DESIGN POINT DIRICH CONDITIONS
-DPOINT 1
+DPOINT 2
 // bottom_vertex_full
 E 1 - NUMDOF 3 ONOFF 1 1 1 VAL 0.0 0.0 0.0 FUNCT none none none
+E 2 - NUMDOF 3 ONOFF 1 1 1 VAL 0.0 0.0 1.0 FUNCT none none 1
+-------------------------------------DESIGN LINE DIRICH CONDITIONS
+DLINE 5
+// top_back
+E 1 - NUMDOF 3 ONOFF 1 0 1 VAL 0.0 0.0 1.0 FUNCT none none 1
+// top_left
+E 2 - NUMDOF 3 ONOFF 0 1 1 VAL 0.0 0.0 1.0 FUNCT none none 1
+// top_bottom
+E 3 - NUMDOF 3 ONOFF 1 1 0 VAL 0.0 0.0 0.0 FUNCT none none none
+// bottom_back
+E 4 - NUMDOF 3 ONOFF 1 0 1 VAL 0.0 0.0 0.0 FUNCT none none none
+// bottom_left
+E 5 - NUMDOF 3 ONOFF 0 1 1 VAL 0.0 0.0 0.0 FUNCT none none none
 -------------------------------------DESIGN SURF DIRICH CONDITIONS
-DSURF  2
+DSURF 4
 // z_top
 E 1 - NUMDOF 3 ONOFF 0 0 1 VAL 0.0 0.0 1.0 FUNCT none none 1
 // z_bottom
 E 2 - NUMDOF 3 ONOFF 0 0 1 VAL 0.0 0.0 0.0 FUNCT none none none
+// z_back
+E 3 - NUMDOF 3 ONOFF 1 0 0 VAL 0.0 0.0 0.0 FUNCT none none none
+// z_left
+E 4 - NUMDOF 3 ONOFF 0 1 0 VAL 0.0 0.0 0.0 FUNCT none none none
 -----------------------------------------------DNODE-NODE TOPOLOGY
 NODE    2 DNODE 1
+NODE    1 DNODE 2
+-----------------------------------------------DLINE-NODE TOPOLOGY
+NODE    1 DLINE 1
+NODE    4 DLINE 1
+NODE    1 DLINE 2
+NODE    5 DLINE 2
+NODE    1 DLINE 3
+NODE    2 DLINE 3
+NODE    2 DLINE 4
+NODE    3 DLINE 4
+NODE    2 DLINE 5
+NODE    6 DLINE 5
 -----------------------------------------------DSURF-NODE TOPOLOGY
 NODE    1 DSURFACE 1
 NODE    4 DSURFACE 1
@@ -77,6 +108,14 @@ NODE    2 DSURFACE 2
 NODE    3 DSURFACE 2
 NODE    6 DSURFACE 2
 NODE    7 DSURFACE 2
+NODE    1 DSURFACE 3
+NODE    2 DSURFACE 3
+NODE    3 DSURFACE 3
+NODE    4 DSURFACE 3
+NODE    1 DSURFACE 4
+NODE    2 DSURFACE 4
+NODE    5 DSURFACE 4
+NODE    6 DSURFACE 4
 -------------------------------------------------------NODE COORDS
 NODE         1 COORD -5.0000000000000000e-01 -5.0000000000000000e-01  5.0000000000000000e-01
 NODE         2 COORD -5.0000000000000000e-01 -5.0000000000000000e-01 -5.0000000000000000e-01
@@ -87,6 +126,4 @@ NODE         6 COORD  5.0000000000000000e-01 -5.0000000000000000e-01 -5.00000000
 NODE         7 COORD  5.0000000000000000e-01  5.0000000000000000e-01 -5.0000000000000000e-01
 NODE         8 COORD  5.0000000000000000e-01  5.0000000000000000e-01  5.0000000000000000e-01
 ------------------------------------------------STRUCTURE ELEMENTS
-   1 SOLIDH8 HEX8  1 2 3 4 5 6 7 8    MAT 1 KINEM nonlinear EAS none
----------------------------------------------------------------END
-// END
+1 SOLID HEX8  1 2 3 4 5 6 7 8  MAT 1 KINEM nonlinear FIBER1 0 0 1.0
