Use more coarse configs for examples

Author: matnoel

examples/Contact/Contact2.py

@@ -29,7 +29,7 @@
     thickness = R / 3  # out-of-plane extent (3D)
     meshSize = R / 20 if dim == 2 else R / 8
 
-    N = 30  # load steps
+    N = 10  # load steps
     delta = 0.1
     penalty = 1e7  # contact stiffness (larger -> less penetration)
 

examples/Contact/Contact3.py

@@ -33,7 +33,7 @@
     thickness = e / 2  # out-of-plane extent (3D)
     mS = t / 5 if dim == 2 else t
 
-    N = 20  # load steps
+    N = 10  # load steps
     delta = 2 * t
     penalty = 1e6  # contact stiffness
 

examples/DIC/DIC1.py

@@ -66,7 +66,7 @@ def Plot(simu: Elastic, result: str, img=None, title="", plotMesh=True):
 
     x0, y0 = 35, 25
     x1, y1 = 423, 814
-    meshSize = (x1 - x0) / 20
+    meshSize = (x1 - x0) / 10
     contour = Domain((x0, y0), (x1, y1), meshSize)
     contour.Plot(ax)
 

examples/Hyperelasticity/AutoDiffPotential.py

@@ -154,7 +154,3 @@ def chi(Ii):
             rel(d2Wde(state, *aux), shipped.Compute_d2Wde(state)),
         )
         print(f"{type(shipped).__name__:24s}" + "".join(f"{e:12.2e}" for e in errors))
-
-    print(
-        "\nHolzapfel-Ogden carries a fibre direction per Gauss point, as the cardiac mesh does."
-    )

examples/Inelasticity/RelaxationPlate.py

@@ -49,7 +49,7 @@
     Models.InElastic.ViscoElastic.Maxwell(g=0.30, tau=1.0),
     Models.InElastic.ViscoElastic.Maxwell(g=0.20, tau=10.0),
 ]
-dt, nStep = 0.5, 40
+dt, nStep = 0.5, 10
 
 
 def Relaxation(n: int) -> float:

examples/LinearizedElasticity/Elas10.py

@@ -48,7 +48,7 @@
 
     # Time parameters
     tMax = 1e-6
-    Nt = 20
+    Nt = 10
     dt = tMax / Nt
 
     # Load parameters

examples/LinearizedElasticity/Elas9.py

@@ -45,7 +45,7 @@
 
     # time
     Tmax = 0.5
-    N = 50
+    N = 25
     dt = Tmax / N
     time = -dt
 

examples/LinearizedElasticity/Homog4.py

@@ -197,7 +197,7 @@ def Get_pairedNodes(
     # ----------------------------------------------
     # Mesh
     # ----------------------------------------------
-    meshSize = 1 / 5
+    meshSize = 1 / 3
 
     p0 = (-1 / 2, -1 / 2)
     p1 = (1 / 2, -1 / 2)
@@ -212,7 +212,7 @@ def Get_pairedNodes(
     inclusion = Circle((0, 0), 2 * r, meshSize, isFilled=True)
     # contour.Plot_Geoms([contour, inclusion])
 
-    elemType = ElemType.PRISM15
+    elemType = ElemType.PRISM18
     mesh = contour.Mesh_Extrude([inclusion], [0, 0, 1], [1 / meshSize], elemType)
     mesh.Translate(*-mesh.center)  # center mesh on 0,0,0
 

examples/LinearizedElasticity/Homog5.py

@@ -141,7 +141,7 @@
     U_e[..., 4] = u13_e
     U_e[..., 5] = u12_e
 
-    matrixType = "mass"
+    matrixType = "rigi"
     wJ_e_pg = mesh.groupElem.Get_weightedJacobian_e_pg(matrixType)
     B_e_pg = mesh.groupElem.Get_B_e_pg(matrixType)
 

examples/LinearizedElasticity/MeshOptim3.py

@@ -60,9 +60,7 @@
     surfLoad = P / h1 / b  # N/mm2
 
     # criteria
-    threshold = (
-        1 / 100 if dim == 2 else 0.04
-    )  # Target error for the optimization process
+    threshold = 0.04  # Target error for the optimization process
     iterMax = 20  # Maximum number of iterations
     coef = 1 / 10  # Scaling coefficient for the optimization process