Merge branch 'main' of github.com:lxmota/Norma.jl

sandialabs · Dec 11, 2024 · dbdb327 · dbdb327
2 parents 260e3bb + 5000e2f
commit dbdb327
Show file tree

Hide file tree

Showing 52 changed files with 450 additions and 221 deletions.
diff --git a/examples/ahead/nonoverlap/clamped/dynamic/clamped-1.yaml b/examples/ahead/nonoverlap/clamped/dynamic/clamped-1.yaml
@@ -7,10 +7,9 @@ model:
     blocks:
       coarse: hyperelastic
     hyperelastic:
-      #model: neohookean
       model: neohookean
       elastic modulus: 1.0e+09
-      Poisson's ratio: 0.25
+      Poisson's ratio: 0.0
       density: 1000.0
 time integrator:
   type: Newmark
@@ -38,12 +37,11 @@ boundary conditions:
     - node set: nsz-
       component: z
       function: "0.0"
-  Schwarz coupling:
+  Schwarz nonoverlap:
     - side set: ssz+
       source: clamped-2.yaml
       source block: fine
       source side set: ssz-
-      coupling type: nonoverlap 
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/nonoverlap/clamped/dynamic/clamped-2.yaml b/examples/ahead/nonoverlap/clamped/dynamic/clamped-2.yaml
@@ -7,10 +7,9 @@ model:
     blocks:
       fine: hyperelastic
     hyperelastic:
-            #model: neohookean
       model: neohookean
       elastic modulus: 1.0e+09
-      Poisson's ratio: 0.25
+      Poisson's ratio: 0.0
       density: 1000.0
 time integrator:
   type: Newmark
@@ -38,12 +37,11 @@ boundary conditions:
     - node set: nsz+
       component: z
       function: "0.0"
-  Schwarz coupling:
+  Schwarz nonoverlap:
     - side set: ssz-
       source: clamped-1.yaml
       source block: coarse
       source side set: ssz+
-      coupling type: nonoverlap 
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/nonoverlap/cuboid/cuboids-hex.jou b/examples/ahead/nonoverlap/cuboid/cuboids-hex.jou
@@ -1,8 +1,28 @@
 ${side = 1.0}
 ${height = 1.5}
 ${overlap_vol_frac = 0.0}
-${h_fine = 1.00}
-${h_coarse = 1.00}
+${h_fine = 0.25}
+${h_coarse = 0.50}
+
+${area = side * side}
+${volume = area * height}
+${overlap_volume = overlap_vol_frac * volume}
+${overlap_length = overlap_volume / area}
+${domain_height = (height + overlap_length) / 2.0}
+${offset = abs(domain_height - overlap_length) / 2.0}
+
+create brick x {side} y {side} z {domain_height}
+move volume 1 z {-offset} include_merged
+#volume 1 scheme tetmesh
+volume 1 size {h_fine}
+mesh volume 1
+block 1 volume 1
+#block 1 element type tetra4
+block 1 name "fine"
+nodeset 1 surface 4
+nodeset 1 name "nsx-"
+nodeset 2 surface 6
+nodeset 2 name "nsx+"
 
 ${area = side * side}
 ${volume = area * height}

diff --git a/examples/ahead/nonoverlap/cuboid/quasistatic/cuboid-1.yaml b/examples/ahead/nonoverlap/cuboid/quasistatic/cuboid-1.yaml
@@ -24,12 +24,11 @@ boundary conditions:
     - node set: nsz-
       component: z
       function: "0.0"
-  Schwarz coupling:
+  Schwarz nonoverlap:
     - side set: ssz+
       source: cuboid-2.yaml
       source block: coarse
       source side set: ssz-
-      coupling type: nonoverlap 
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/nonoverlap/cuboid/quasistatic/cuboid-2.yaml b/examples/ahead/nonoverlap/cuboid/quasistatic/cuboid-2.yaml
@@ -24,12 +24,11 @@ boundary conditions:
     - node set: nsz+
       component: z
       function: "1.0 * t"
-  Schwarz coupling:
+  Schwarz nonoverlap:
     - side set: ssz-
       source: cuboid-1.yaml
       source block: fine
       source side set: ssz+
-      coupling type: nonoverlap 
 solver:
 solver:
   type: Hessian minimizer

diff --git a/examples/ahead/nonoverlap/laser-weld/laser-weld.jou b/examples/ahead/nonoverlap/laser-weld/laser-weld.jou
@@ -268,9 +268,6 @@ refine surface 23 26 numsplit 1 depth 1
 
 mesh volume 1 2 5
 
-nodeset 1 surface 30 37 26 8 20 41 47
-nodeset 1 name "surface_symmetry"
-
 nodeset 2 surface 16
 nodeset 2 name "surface_negative_y"
 
@@ -291,9 +288,6 @@ volume all size {sizeElementFine}
 
 mesh volume all 
 
-nodeset 1 surface 30 37 26 8 20 41 47
-nodeset 1 name "surface_symmetry"
-
 nodeset 2 surface 16
 nodeset 2 name "surface_negative_y"
 
@@ -320,15 +314,11 @@ mesh volume 1 2 4 5 6
 
 reset nodeset
 
-nodeset 1 surface 30 37 26 8 20 41 47 
+#nodeset 2 surface 29
+#nodeset 2 name "nodeset_negative_y"
 
-nodeset 1 name "nodeset_symmetry"
-
-nodeset 2 surface 29
-nodeset 2 name "nodeset_negative_y"
-
-nodeset 3 surface 39
-nodeset 3 name "nodeset_positive_y"
+#nodeset 3 surface 39
+#nodeset 3 name "nodeset_positive_y"
 
 #sideset 2 surface 29
 #sideset 2 name "sideset_negative_y"
@@ -352,9 +342,6 @@ mesh volume 3 6
 
 reset nodeset
 
-nodeset 1 surface 30 37
-
-nodeset 1 name "nodeset_symmetry"
 
 nodeset 2 surface 16
 nodeset 2 name "nodeset_negative_y"
@@ -380,10 +367,6 @@ mesh volume 4 7
 
 reset nodeset
 
-nodeset 1 surface 41 47
-
-nodeset 1 name "nodeset_symmetry"
-
 nodeset 2 surface 17
 nodeset 2 name "nodeset_negative_y"
 
@@ -423,7 +406,7 @@ reset
 import mesh geometry "gauge.g" 
 sideset 1 surface 17 
 sideset 1 name "sideset_negative_y" 
-sideset 2 surface 2 
+sideset 2 surface 1 
 sideset 2 name "sideset_positive_y" 
 block all element type HEX8
 export mesh "gauge.g" block all overwrite
diff --git a/examples/ahead/nonoverlap/laser-weld/quasistatic/gauge.yaml b/examples/ahead/nonoverlap/laser-weld/quasistatic/gauge.yaml
@@ -16,19 +16,17 @@ model:
 time integrator:
   type: quasi static
 boundary conditions:
-  Schwarz coupling:
+  Schwarz nonoverlap:
     - side set: sideset_negative_y
       source: holder-0.yaml
       source block: holder_0
       source side set: sideset_positive_y
       search tolerance: 5.0e-02
-      coupling type: nonoverlap
     - side set: sideset_positive_y
       source: holder-1.yaml
       source block: holder_1
       source side set: sideset_negative_y
       search tolerance: 5.0e-02
-      coupling type: nonoverlap
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/nonoverlap/laser-weld/quasistatic/holder-0.yaml b/examples/ahead/nonoverlap/laser-weld/quasistatic/holder-0.yaml
@@ -24,13 +24,12 @@ boundary conditions:
     - node set: nodeset_negative_y
       component: z
       function: "0.0"
-  Schwarz coupling:
+  Schwarz nonoverlap:
     - side set: sideset_positive_y
       source: gauge.yaml
       source block: region_weld_0
       source side set: sideset_negative_y
       search tolerance: 5.0e-02
-      coupling type: nonoverlap
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/nonoverlap/laser-weld/quasistatic/holder-1.yaml b/examples/ahead/nonoverlap/laser-weld/quasistatic/holder-1.yaml
@@ -24,13 +24,12 @@ boundary conditions:
     - node set: nodeset_positive_y
       component: z
       function: "0.0"
-  Schwarz coupling:
+  Schwarz nonoverlap:
     - side set: sideset_negative_y
       source: gauge.yaml
       source block: region_weld_1
       source side set: sideset_positive_y
       search tolerance: 5.0e-02
-      coupling type: nonoverlap
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/nonoverlap/notched-cylinder/dynamic/notched-cylinder-1.yaml b/examples/ahead/nonoverlap/notched-cylinder/dynamic/notched-cylinder-1.yaml
@@ -26,13 +26,12 @@ boundary conditions:
     - node set: -Z_bottom
       component: z
       function: "0.0"
-  Schwarz coupling:
+  Schwarz nonoverlap:
     - side set: +Z_bottomSS
       source: notched-cylinder-2.yaml
       source block: coarse
       source side set: -Z_topSS
       search tolerance: 5.0e-02
-      coupling type: overlap
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/nonoverlap/notched-cylinder/dynamic/notched-cylinder-2.yaml b/examples/ahead/nonoverlap/notched-cylinder/dynamic/notched-cylinder-2.yaml
@@ -26,12 +26,11 @@ boundary conditions:
     - node set: +Z_top
       component: z
       function: "0.0064 * t"
-  Schwarz coupling:
+  Schwarz nonoverlap:
     - side set: -Z_topSS
       source: notched-cylinder-1.yaml
       source side set: +Z_bottomSS
       source block: fine
-      coupling type: overlap
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/nonoverlap/notched-cylinder/quasistatic/notched-cylinder-1.yaml b/examples/ahead/nonoverlap/notched-cylinder/quasistatic/notched-cylinder-1.yaml
@@ -24,13 +24,12 @@ boundary conditions:
     - node set: -Z_bottom
       component: z
       function: "0.0"
-  Schwarz coupling:
+  Schwarz nonoverlap:
     - side set: +Z_bottomSS
       source: notched-cylinder-2.yaml
       source block: coarse
       source side set: -Z_topSS
       search tolerance: 5.0e-02
-      coupling type: nonoverlap
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/nonoverlap/notched-cylinder/quasistatic/notched-cylinder-2.yaml b/examples/ahead/nonoverlap/notched-cylinder/quasistatic/notched-cylinder-2.yaml
@@ -24,12 +24,11 @@ boundary conditions:
     - node set: +Z_top
       component: z
       function: "0.0064 * t"
-  Schwarz coupling:
+  Schwarz nonoverlap:
     - side set: -Z_topSS
       source: notched-cylinder-1.yaml
       source side set: +Z_bottomSS
       source block: fine
-      coupling type: nonoverlap
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/overlap/clamped/dynamic/clamped-1.yaml b/examples/ahead/overlap/clamped/dynamic/clamped-1.yaml
@@ -37,12 +37,11 @@ boundary conditions:
     - node set: nsz-
       component: z
       function: "0.0"
-  Schwarz coupling:
+  Schwarz overlap:
     - side set: ssz+
       source: clamped-2.yaml
       source block: fine
       source side set: ssz-
-      coupling type: overlap 
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/overlap/clamped/dynamic/clamped-2.yaml b/examples/ahead/overlap/clamped/dynamic/clamped-2.yaml
@@ -37,12 +37,11 @@ boundary conditions:
     - node set: nsz+
       component: z
       function: "0.0"
-  Schwarz coupling:
+  Schwarz overlap:
     - side set: ssz-
       source: clamped-1.yaml
       source block: coarse
       source side set: ssz+
-      coupling type: overlap 
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/overlap/cuboid/dynamic/cuboid-1.yaml b/examples/ahead/overlap/cuboid/dynamic/cuboid-1.yaml
@@ -26,12 +26,11 @@ boundary conditions:
     - node set: nsz-
       component: z
       function: "0.0"
-  Schwarz coupling:
+  Schwarz overlap:
     - side set: ssz+
       source: cuboid-2.yaml
       source block: coarse
       source side set: ssz-
-      coupling type: overlap 
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/overlap/cuboid/dynamic/cuboid-2.yaml b/examples/ahead/overlap/cuboid/dynamic/cuboid-2.yaml
@@ -26,12 +26,11 @@ boundary conditions:
     - node set: nsz+
       component: z
       function: "1.0 * t"
-  Schwarz coupling:
+  Schwarz overlap:
     - side set: ssz-
       source: cuboid-1.yaml
       source block: fine
       source side set: ssz+
-      coupling type: overlap 
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/overlap/cuboid/quasistatic/cuboid-1.yaml b/examples/ahead/overlap/cuboid/quasistatic/cuboid-1.yaml
@@ -24,12 +24,11 @@ boundary conditions:
     - node set: nsz-
       component: z
       function: "0.0"
-  Schwarz coupling:
+  Schwarz overlap:
     - side set: ssz+
       source: cuboid-2.yaml
       source block: coarse
       source side set: ssz-
-      coupling type: overlap 
 solver:
   type: Hessian minimizer
   step: full Newton

diff --git a/examples/ahead/overlap/cuboid/quasistatic/cuboid-2.yaml b/examples/ahead/overlap/cuboid/quasistatic/cuboid-2.yaml
@@ -24,12 +24,11 @@ boundary conditions:
     - node set: nsz+
       component: z
       function: "1.0 * t"
-  Schwarz coupling:
+  Schwarz overlap:
     - side set: ssz-
       source: cuboid-1.yaml
       source block: fine
       source side set: ssz+
-      coupling type: overlap 
 solver:
 solver:
   type: Hessian minimizer

diff --git a/examples/ahead/overlap/laser-weld/dynamic/gauge.yaml b/examples/ahead/overlap/laser-weld/dynamic/gauge.yaml
@@ -22,17 +22,15 @@ boundary conditions:
     - node set: nodeset_symmetry
       component: x
       function: "0.0"
-  Schwarz coupling:
+  Schwarz overlap:
     - side set: sideset_negative_y
       source: holder-0.yaml
       source block: holder_0
       source side set: sideset_positive_y
-      coupling type: overlap
     - side set: sideset_positive_y
       source: holder-1.yaml
       source block: holder_1
       source side set: sideset_negative_y
-      coupling type: overlap
 solver:
   type: Hessian minimizer
   step: full Newton