Note

Go to the end to download the full example code.

Shear#

Damage simulation for a plate subjected to shear.

Shear
  • $\phi$
  • TRI3: Ne = 2923, Nn = 1540
  • Shear
  • Boundary conditions
  • Shear
  • Summary
/home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/stable/examples/PhaseField/results/Shear2D/Test/Miehe_AT1_optimMesh
   1: 0.050µm, [0.00e+00;0.00e+00], 1:0.040s, tol=1.00e+00
   2: 0.100µm, [0.00e+00;0.00e+00], 1:0.030s, tol=7.50e-01  0.25% -> 12.05s
   3: 0.150µm, [0.00e+00;0.00e+00], 1:0.030s, tol=5.56e-01  0.50% -> 11.92s
   4: 0.200µm, [0.00e+00;0.00e+00], 1:0.030s, tol=4.38e-01  0.75% -> 11.81s
   5: 0.250µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.60e-01  1.00% -> 11.92s
   6: 0.300µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.06e-01  1.25% -> 11.83s
   7: 0.350µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.65e-01  1.50% -> 11.75s
   8: 0.400µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.34e-01  1.75% -> 11.86s
   9: 0.450µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.10e-01  2.00% -> 11.87s
  10: 0.500µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.90e-01  2.25% -> 11.67s
  11: 0.550µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.74e-01  2.50% -> 11.67s
  12: 0.600µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.60e-01  2.75% -> 11.79s
  13: 0.650µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.48e-01  3.00% -> 11.53s
  14: 0.700µm, [0.00e+00;0.00e+00], 1:0.031s, tol=1.38e-01  3.25% -> 11.81s
  15: 0.750µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.29e-01  3.50% -> 11.51s
  16: 0.800µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.21e-01  3.75% -> 11.57s
  17: 0.850µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.14e-01  4.00% -> 11.41s
  18: 0.900µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.08e-01  4.25% -> 11.51s
  19: 0.950µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.02e-01  4.50% -> 11.56s
  20: 1.000µm, [0.00e+00;0.00e+00], 1:0.030s, tol=9.75e-02  4.75% -> 11.39s
  21: 1.050µm, [0.00e+00;0.00e+00], 1:0.030s, tol=9.30e-02  5.00% -> 11.21s
  22: 1.100µm, [0.00e+00;0.00e+00], 1:0.030s, tol=8.88e-02  5.25% -> 11.35s
  23: 1.150µm, [0.00e+00;0.00e+00], 1:0.030s, tol=8.51e-02  5.50% -> 11.20s
  24: 1.200µm, [0.00e+00;0.00e+00], 1:0.030s, tol=8.16e-02  5.75% -> 11.15s
  25: 1.250µm, [0.00e+00;0.00e+00], 1:0.030s, tol=7.84e-02  6.00% -> 11.19s
  26: 1.300µm, [0.00e+00;0.00e+00], 1:0.030s, tol=7.54e-02  6.25% -> 11.16s
  27: 1.350µm, [0.00e+00;0.00e+00], 1:0.030s, tol=7.27e-02  6.50% -> 11.23s
  28: 1.400µm, [0.00e+00;0.00e+00], 1:0.037s, tol=7.02e-02  6.75% -> 13.98s
  29: 1.450µm, [0.00e+00;0.00e+00], 1:0.032s, tol=6.78e-02  7.00% -> 11.78s
  30: 1.500µm, [0.00e+00;0.00e+00], 1:0.031s, tol=6.56e-02  7.25% -> 11.60s
  31: 1.550µm, [0.00e+00;0.00e+00], 1:0.031s, tol=6.35e-02  7.50% -> 11.57s
  32: 1.600µm, [0.00e+00;0.00e+00], 1:0.031s, tol=6.15e-02  7.75% -> 11.41s
  33: 1.650µm, [0.00e+00;0.00e+00], 1:0.031s, tol=5.97e-02  8.00% -> 11.41s
  34: 1.700µm, [0.00e+00;0.00e+00], 1:0.031s, tol=5.80e-02  8.25% -> 11.39s
  35: 1.750µm, [0.00e+00;0.00e+00], 1:0.030s, tol=5.63e-02  8.50% -> 10.97s
  36: 1.800µm, [0.00e+00;0.00e+00], 1:0.030s, tol=5.48e-02  8.75% -> 11.06s
  37: 1.850µm, [0.00e+00;0.00e+00], 1:0.030s, tol=5.33e-02  9.00% -> 11.09s
  38: 1.900µm, [0.00e+00;0.00e+00], 1:0.031s, tol=5.19e-02  9.25% -> 11.33s
  39: 1.950µm, [0.00e+00;0.00e+00], 1:0.030s, tol=5.06e-02  9.50% -> 10.81s
  40: 2.000µm, [0.00e+00;0.00e+00], 1:0.030s, tol=4.94e-02  9.75% -> 10.78s
  41: 2.050µm, [0.00e+00;0.00e+00], 1:0.030s, tol=4.82e-02  10.00% -> 10.66s
  42: 2.100µm, [0.00e+00;0.00e+00], 1:0.030s, tol=4.71e-02  10.25% -> 10.67s
  43: 2.150µm, [0.00e+00;0.00e+00], 1:0.030s, tol=4.60e-02  10.50% -> 10.77s
  44: 2.200µm, [0.00e+00;0.00e+00], 1:0.030s, tol=4.49e-02  10.75% -> 10.65s
  45: 2.250µm, [0.00e+00;0.00e+00], 1:0.030s, tol=4.40e-02  11.00% -> 10.63s
  46: 2.300µm, [0.00e+00;0.00e+00], 1:0.030s, tol=4.30e-02  11.25% -> 10.60s
  47: 2.350µm, [0.00e+00;0.00e+00], 1:0.030s, tol=4.21e-02  11.50% -> 10.66s
  48: 2.400µm, [0.00e+00;0.00e+00], 1:0.030s, tol=4.12e-02  11.75% -> 10.59s
  49: 2.450µm, [0.00e+00;0.00e+00], 1:0.030s, tol=4.04e-02  12.00% -> 10.55s
  50: 2.500µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.96e-02  12.25% -> 10.61s
  51: 2.550µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.88e-02  12.50% -> 10.54s
  52: 2.600µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.81e-02  12.75% -> 10.37s
  53: 2.650µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.74e-02  13.00% -> 10.48s
  54: 2.700µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.67e-02  13.25% -> 10.44s
  55: 2.750µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.60e-02  13.50% -> 10.26s
  56: 2.800µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.54e-02  13.75% -> 10.34s
  57: 2.850µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.48e-02  14.00% -> 10.24s
  58: 2.900µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.42e-02  14.25% -> 10.13s
  59: 2.950µm, [0.00e+00;0.00e+00], 1:0.031s, tol=3.36e-02  14.50% -> 10.52s
  60: 3.000µm, [0.00e+00;0.00e+00], 1:0.031s, tol=3.31e-02  14.75% -> 10.51s
  61: 3.050µm, [0.00e+00;0.00e+00], 1:0.031s, tol=3.25e-02  15.00% -> 10.41s
  62: 3.100µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.20e-02  15.25% -> 10.10s
  63: 3.150µm, [0.00e+00;0.00e+00], 1:0.029s, tol=3.15e-02  15.50% -> 9.96s
  64: 3.200µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.10e-02  15.75% -> 10.24s
  65: 3.250µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.05e-02  16.00% -> 10.15s
  66: 3.300µm, [0.00e+00;0.00e+00], 1:0.030s, tol=3.01e-02  16.25% -> 10.04s
  67: 3.350µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.96e-02  16.50% -> 9.99s
  68: 3.400µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.92e-02  16.75% -> 9.86s
  69: 3.450µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.88e-02  17.00% -> 9.81s
  70: 3.500µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.84e-02  17.25% -> 10.04s
  71: 3.550µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.80e-02  17.50% -> 9.94s
  72: 3.600µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.76e-02  17.75% -> 9.95s
  73: 3.650µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.72e-02  18.00% -> 9.76s
  74: 3.700µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.68e-02  18.25% -> 9.73s
  75: 3.750µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.65e-02  18.50% -> 9.74s
  76: 3.800µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.61e-02  18.75% -> 9.60s
  77: 3.850µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.58e-02  19.00% -> 9.84s
  78: 3.900µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.55e-02  19.25% -> 9.56s
  79: 3.950µm, [0.00e+00;0.00e+00], 1:0.029s, tol=2.52e-02  19.50% -> 9.49s
  80: 4.000µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.48e-02  19.75% -> 9.53s
  81: 4.050µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.45e-02  20.00% -> 9.50s
  82: 4.100µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.42e-02  20.25% -> 9.43s
  83: 4.150µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.40e-02  20.50% -> 9.43s
  84: 4.200µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.37e-02  20.75% -> 9.38s
  85: 4.250µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.34e-02  21.00% -> 9.36s
  86: 4.300µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.31e-02  21.25% -> 9.46s
  87: 4.350µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.29e-02  21.50% -> 9.37s
  88: 4.400µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.26e-02  21.75% -> 9.35s
  89: 4.450µm, [0.00e+00;0.00e+00], 1:0.029s, tol=2.23e-02  22.00% -> 9.20s
  90: 4.500µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.21e-02  22.25% -> 9.24s
  91: 4.550µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.19e-02  22.50% -> 9.16s
  92: 4.600µm, [0.00e+00;0.00e+00], 1:0.031s, tol=2.16e-02  22.75% -> 9.43s
  93: 4.650µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.14e-02  23.00% -> 9.37s
  94: 4.700µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.12e-02  23.25% -> 9.30s
  95: 4.750µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.09e-02  23.50% -> 9.11s
  96: 4.800µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.07e-02  23.75% -> 9.05s
  97: 4.850µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.05e-02  24.00% -> 9.10s
  98: 4.900µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.03e-02  24.25% -> 9.07s
  99: 4.950µm, [0.00e+00;0.00e+00], 1:0.030s, tol=2.01e-02  24.50% -> 9.04s
 100: 5.000µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.99e-02  24.75% -> 9.10s
 101: 5.050µm, [0.00e+00;0.00e+00], 1:0.031s, tol=1.97e-02  25.00% -> 9.45s
 102: 5.100µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.95e-02  25.25% -> 9.02s
 103: 5.150µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.93e-02  25.50% -> 8.83s
 104: 5.200µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.91e-02  25.75% -> 8.98s
 105: 5.250µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.90e-02  26.00% -> 8.84s
 106: 5.300µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.88e-02  26.25% -> 8.82s
 107: 5.350µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.86e-02  26.50% -> 8.69s
 108: 5.400µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.84e-02  26.75% -> 8.66s
 109: 5.450µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.83e-02  27.00% -> 8.68s
 110: 5.500µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.81e-02  27.25% -> 8.66s
 111: 5.550µm, [0.00e+00;0.00e+00], 1:0.031s, tol=1.79e-02  27.50% -> 8.85s
 112: 5.600µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.78e-02  27.75% -> 8.73s
 113: 5.650µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.76e-02  28.00% -> 8.76s
 114: 5.700µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.75e-02  28.25% -> 8.53s
 115: 5.750µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.73e-02  28.50% -> 8.47s
 116: 5.800µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.72e-02  28.75% -> 8.41s
 117: 5.850µm, [0.00e+00;0.00e+00], 1:0.029s, tol=1.70e-02  29.00% -> 8.36s
 118: 5.900µm, [0.00e+00;0.00e+00], 1:0.029s, tol=1.69e-02  29.25% -> 8.34s
 119: 5.950µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.67e-02  29.50% -> 8.35s
 120: 6.000µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.66e-02  29.75% -> 8.31s
 121: 6.050µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.65e-02  30.00% -> 8.39s
 122: 6.100µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.63e-02  30.25% -> 8.34s
 123: 6.150µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.62e-02  30.50% -> 8.27s
 124: 6.200µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.61e-02  30.75% -> 8.29s
 125: 6.250µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.59e-02  31.00% -> 8.19s
 126: 6.300µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.58e-02  31.25% -> 8.16s
 127: 6.350µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.57e-02  31.50% -> 8.09s
 128: 6.400µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.56e-02  31.75% -> 8.07s
 129: 6.450µm, [0.00e+00;0.00e+00], 1:0.031s, tol=1.54e-02  32.00% -> 8.32s
 130: 6.500µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.53e-02  32.25% -> 8.07s
 131: 6.550µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.52e-02  32.50% -> 8.19s
 132: 6.600µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.51e-02  32.75% -> 8.06s
 133: 6.650µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.50e-02  33.00% -> 7.97s
 134: 6.700µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.49e-02  33.25% -> 8.03s
 135: 6.750µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.48e-02  33.50% -> 8.05s
 136: 6.800µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.47e-02  33.75% -> 7.94s
 137: 6.850µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.45e-02  34.00% -> 7.86s
 138: 6.900µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.44e-02  34.25% -> 7.91s
 139: 6.950µm, [0.00e+00;0.00e+00], 1:0.030s, tol=1.43e-02  34.50% -> 7.85s
 140: 7.000µm, [6.00e-07;1.59e-03], 1:0.030s, tol=1.42e-02  34.75% -> 7.91s
 141: 7.050µm, [2.22e-06;6.20e-03], 1:0.030s, tol=1.39e-02  35.00% -> 7.90s
 142: 7.100µm, [4.16e-06;1.23e-02], 1:0.030s, tol=1.38e-02  35.25% -> 7.76s
 143: 7.150µm, [6.13e-06;1.92e-02], 1:0.030s, tol=1.37e-02  35.50% -> 7.71s
 144: 7.200µm, [7.97e-06;2.65e-02], 1:0.030s, tol=1.35e-02  35.75% -> 7.73s
 145: 7.250µm, [9.65e-06;3.41e-02], 1:0.030s, tol=1.35e-02  36.00% -> 7.69s
 146: 7.300µm, [1.12e-05;4.18e-02], 1:0.031s, tol=1.34e-02  36.25% -> 7.83s
 147: 7.350µm, [1.25e-05;4.96e-02], 1:0.030s, tol=1.33e-02  36.50% -> 7.65s
 148: 7.400µm, [1.36e-05;5.75e-02], 1:0.030s, tol=1.32e-02  36.75% -> 7.63s
 149: 7.450µm, [1.46e-05;6.55e-02], 1:0.031s, tol=1.31e-02  37.00% -> 7.76s
 150: 7.500µm, [1.55e-05;7.35e-02], 1:0.031s, tol=1.30e-02  37.25% -> 7.71s
 151: 7.550µm, [1.62e-05;8.16e-02], 1:0.030s, tol=1.29e-02  37.50% -> 7.60s
 152: 7.600µm, [1.68e-05;8.98e-02], 1:0.030s, tol=1.28e-02  37.75% -> 7.54s
 153: 7.650µm, [1.72e-05;9.80e-02], 1:0.030s, tol=1.28e-02  38.00% -> 7.40s
 154: 7.700µm, [1.76e-05;1.06e-01], 1:0.030s, tol=1.27e-02  38.25% -> 7.42s
 155: 7.750µm, [1.79e-05;1.15e-01], 1:0.030s, tol=1.26e-02  38.50% -> 7.43s
 156: 7.800µm, [1.81e-05;1.23e-01], 1:0.030s, tol=1.25e-02  38.75% -> 7.42s
 157: 7.850µm, [1.82e-05;1.32e-01], 1:0.031s, tol=1.24e-02  39.00% -> 7.46s
 158: 7.900µm, [1.83e-05;1.40e-01], 1:0.030s, tol=1.24e-02  39.25% -> 7.30s
 159: 7.950µm, [1.83e-05;1.49e-01], 1:0.030s, tol=1.23e-02  39.50% -> 7.27s
 160: 8.000µm, [1.82e-05;1.58e-01], 1:0.030s, tol=1.22e-02  39.75% -> 7.21s
 161: 8.050µm, [1.81e-05;1.67e-01], 1:0.031s, tol=1.21e-02  40.00% -> 7.43s
 162: 8.100µm, [1.80e-05;1.76e-01], 1:0.030s, tol=1.21e-02  40.25% -> 7.27s
 163: 8.150µm, [1.78e-05;1.85e-01], 1:0.031s, tol=1.20e-02  40.50% -> 7.28s
 164: 8.200µm, [1.75e-05;1.94e-01], 1:0.031s, tol=1.19e-02  40.75% -> 7.44s
 165: 8.250µm, [1.73e-05;2.03e-01], 1:0.031s, tol=1.18e-02  41.00% -> 7.34s
 166: 8.300µm, [1.70e-05;2.12e-01], 1:0.030s, tol=1.18e-02  41.25% -> 7.04s
 167: 8.350µm, [1.67e-05;2.22e-01], 1:0.030s, tol=1.17e-02  41.50% -> 7.05s
 168: 8.400µm, [1.64e-05;2.31e-01], 1:0.030s, tol=1.16e-02  41.75% -> 7.00s
 169: 8.450µm, [1.60e-05;2.41e-01], 1:0.030s, tol=1.16e-02  42.00% -> 7.04s
 170: 8.500µm, [1.57e-05;2.51e-01], 1:0.030s, tol=1.15e-02  42.25% -> 6.99s
 171: 8.550µm, [1.53e-05;2.60e-01], 1:0.030s, tol=1.14e-02  42.50% -> 6.88s
 172: 8.600µm, [1.49e-05;2.70e-01], 1:0.030s, tol=1.14e-02  42.75% -> 6.94s
 173: 8.650µm, [1.45e-05;2.80e-01], 1:0.030s, tol=1.13e-02  43.00% -> 6.83s
 174: 8.700µm, [1.42e-05;2.90e-01], 1:0.030s, tol=1.13e-02  43.25% -> 6.82s
 175: 8.750µm, [1.38e-05;3.01e-01], 1:0.030s, tol=1.12e-02  43.50% -> 6.83s
 176: 8.800µm, [1.34e-05;3.11e-01], 1:0.031s, tol=1.11e-02  43.75% -> 6.97s
 177: 8.850µm, [1.30e-05;3.21e-01], 1:0.031s, tol=1.11e-02  44.00% -> 6.94s
 178: 8.900µm, [1.26e-05;3.32e-01], 1:0.031s, tol=1.10e-02  44.25% -> 6.89s
 179: 8.950µm, [1.22e-05;3.43e-01], 1:0.030s, tol=1.10e-02  44.50% -> 6.68s
 180: 9.000µm, [1.19e-05;3.56e-01], 1:0.030s, tol=1.09e-02  44.75% -> 6.68s
 181: 9.050µm, [1.15e-05;3.71e-01], 1:0.030s, tol=1.08e-02  45.00% -> 6.71s
 182: 9.100µm, [1.12e-05;3.86e-01], 1:0.031s, tol=1.07e-02  45.25% -> 6.69s
 183: 9.150µm, [1.09e-05;4.01e-01], 1:0.030s, tol=1.07e-02  45.50% -> 6.63s
 184: 9.200µm, [1.05e-05;4.17e-01], 1:0.030s, tol=1.06e-02  45.75% -> 6.51s
 185: 9.250µm, [1.02e-05;4.32e-01], 1:0.030s, tol=1.06e-02  46.00% -> 6.55s
 186: 9.300µm, [9.84e-06;4.47e-01], 1:0.030s, tol=1.06e-02  46.25% -> 6.52s
 187: 9.350µm, [9.49e-06;4.62e-01], 1:0.030s, tol=1.05e-02  46.50% -> 6.44s
 188: 9.400µm, [9.27e-06;4.77e-01], 1:0.030s, tol=1.04e-02  46.75% -> 6.49s
 189: 9.450µm, [9.10e-06;4.94e-01], 1:0.030s, tol=1.03e-02  47.00% -> 6.40s
 190: 9.500µm, [8.97e-06;5.14e-01], 1:0.030s, tol=1.03e-02  47.25% -> 6.33s
 191: 9.550µm, [8.88e-06;5.39e-01], 1:0.030s, tol=1.02e-02  47.50% -> 6.34s
 192: 9.600µm, [8.83e-06;5.71e-01], 1:0.030s, tol=1.03e-02  47.75% -> 6.30s
 193: 9.650µm, [8.78e-06;6.08e-01], 1:0.030s, tol=1.04e-02  48.00% -> 6.33s
 194: 9.700µm, [8.76e-06;6.49e-01], 1:0.030s, tol=1.07e-02  48.25% -> 6.25s
 195: 9.750µm, [8.75e-06;6.95e-01], 1:0.030s, tol=1.13e-02  48.50% -> 6.25s
 196: 9.800µm, [9.05e-06;7.43e-01], 1:0.031s, tol=1.23e-02  48.75% -> 6.27s
 197: 9.850µm, [9.39e-06;7.94e-01], 1:0.030s, tol=1.47e-02  49.00% -> 6.17s
 198: 9.900µm, [9.76e-06;8.46e-01], 1:0.031s, tol=1.92e-02  49.25% -> 6.26s
 199: 9.950µm, [9.84e-06;8.94e-01], 1:0.030s, tol=2.57e-02  49.50% -> 6.08s
 200: 10.000µm, [9.59e-06;9.28e-01], 1:0.030s, tol=3.11e-02  49.75% -> 6.12s
 201: 10.050µm, [9.14e-06;9.49e-01], 1:0.030s, tol=3.15e-02  50.00% -> 6.07s
 202: 10.100µm, [9.06e-06;9.58e-01], 1:0.031s, tol=2.98e-02  50.25% -> 6.07s
 203: 10.150µm, [9.00e-06;9.60e-01], 1:0.031s, tol=3.01e-02  50.50% -> 6.12s
 204: 10.200µm, [9.21e-06;9.60e-01], 1:0.030s, tol=3.37e-02  50.75% -> 5.93s
 205: 10.250µm, [9.27e-06;9.60e-01], 1:0.030s, tol=4.01e-02  51.00% -> 5.98s
 206: 10.300µm, [9.16e-06;9.71e-01], 1:0.030s, tol=4.63e-02  51.25% -> 5.91s
 207: 10.350µm, [8.87e-06;9.82e-01], 1:0.031s, tol=4.81e-02  51.50% -> 5.92s
 208: 10.400µm, [8.45e-06;9.89e-01], 1:0.031s, tol=4.32e-02  51.75% -> 5.94s
 209: 10.450µm, [8.11e-06;9.92e-01], 1:0.031s, tol=3.51e-02  52.00% -> 5.87s
 210: 10.500µm, [8.09e-06;9.94e-01], 1:0.030s, tol=3.08e-02  52.25% -> 5.80s
 211: 10.550µm, [8.04e-06;9.95e-01], 1:0.031s, tol=2.91e-02  52.50% -> 5.81s
 212: 10.600µm, [7.90e-06;9.95e-01], 1:0.030s, tol=2.80e-02  52.75% -> 5.73s
 213: 10.650µm, [7.80e-06;9.95e-01], 1:0.031s, tol=2.75e-02  53.00% -> 5.87s
 214: 10.700µm, [7.73e-06;9.94e-01], 1:0.032s, tol=2.85e-02  53.25% -> 5.93s
 215: 10.750µm, [7.62e-06;9.94e-01], 1:0.032s, tol=3.07e-02  53.50% -> 5.90s
 216: 10.800µm, [7.45e-06;9.93e-01], 1:0.031s, tol=3.36e-02  53.75% -> 5.71s
 217: 10.850µm, [7.23e-06;9.93e-01], 1:0.030s, tol=3.59e-02  54.00% -> 5.55s
 218: 10.900µm, [6.95e-06;9.92e-01], 1:0.030s, tol=3.61e-02  54.25% -> 5.55s
 219: 10.950µm, [6.71e-06;9.92e-01], 1:0.031s, tol=3.49e-02  54.50% -> 5.64s
 220: 11.000µm, [6.54e-06;9.92e-01], 1:0.031s, tol=3.27e-02  54.75% -> 5.53s
 221: 11.050µm, [6.47e-06;9.92e-01], 1:0.030s, tol=3.08e-02  55.00% -> 5.41s
 222: 11.100µm, [6.45e-06;9.92e-01], 1:0.030s, tol=3.00e-02  55.25% -> 5.39s
 223: 11.150µm, [6.36e-06;9.93e-01], 1:0.030s, tol=2.96e-02  55.50% -> 5.42s
 224: 11.200µm, [6.21e-06;9.93e-01], 1:0.030s, tol=2.93e-02  55.75% -> 5.31s
 225: 11.250µm, [6.16e-06;9.93e-01], 1:0.030s, tol=3.12e-02  56.00% -> 5.28s
 226: 11.300µm, [6.13e-06;9.93e-01], 1:0.030s, tol=3.48e-02  56.25% -> 5.33s
 227: 11.350µm, [6.03e-06;9.94e-01], 1:0.031s, tol=3.87e-02  56.50% -> 5.37s
 228: 11.400µm, [5.87e-06;9.94e-01], 1:0.031s, tol=4.11e-02  56.75% -> 5.40s
 229: 11.450µm, [5.71e-06;9.94e-01], 1:0.030s, tol=4.15e-02  57.00% -> 5.18s
 230: 11.500µm, [5.65e-06;9.95e-01], 1:0.030s, tol=4.01e-02  57.25% -> 5.17s
 231: 11.550µm, [5.70e-06;9.95e-01], 1:0.030s, tol=3.84e-02  57.50% -> 5.12s
 232: 11.600µm, [5.74e-06;9.95e-01], 1:0.030s, tol=3.71e-02  57.75% -> 5.13s
 233: 11.650µm, [5.66e-06;9.95e-01], 1:0.030s, tol=3.53e-02  58.00% -> 5.08s
 234: 11.700µm, [5.61e-06;9.95e-01], 1:0.030s, tol=3.48e-02  58.25% -> 5.03s
 235: 11.750µm, [5.60e-06;9.95e-01], 1:0.030s, tol=3.73e-02  58.50% -> 5.01s
 236: 11.800µm, [5.51e-06;9.96e-01], 1:0.030s, tol=4.00e-02  58.75% -> 4.99s
 237: 11.850µm, [5.41e-06;9.96e-01], 1:0.030s, tol=4.25e-02  59.00% -> 4.93s
 238: 11.900µm, [5.38e-06;9.96e-01], 1:0.030s, tol=4.26e-02  59.25% -> 4.95s
 239: 11.950µm, [5.51e-06;9.96e-01], 1:0.030s, tol=4.20e-02  59.50% -> 4.87s
 240: 12.000µm, [5.60e-06;9.96e-01], 1:0.030s, tol=4.03e-02  59.75% -> 4.86s
 241: 12.050µm, [5.53e-06;9.96e-01], 1:0.030s, tol=3.75e-02  60.00% -> 4.83s
 242: 12.100µm, [5.56e-06;9.96e-01], 1:0.030s, tol=3.75e-02  60.25% -> 4.84s
 243: 12.150µm, [5.52e-06;9.97e-01], 1:0.030s, tol=3.99e-02  60.50% -> 4.77s
 244: 12.200µm, [5.47e-06;9.97e-01], 1:0.030s, tol=4.29e-02  60.75% -> 4.76s
 245: 12.250µm, [5.50e-06;9.97e-01], 1:0.030s, tol=4.42e-02  61.00% -> 4.71s
 246: 12.300µm, [5.73e-06;9.97e-01], 1:0.030s, tol=4.43e-02  61.25% -> 4.69s
 247: 12.350µm, [5.83e-06;9.97e-01], 1:0.030s, tol=4.23e-02  61.50% -> 4.61s
 248: 12.400µm, [5.81e-06;9.97e-01], 1:0.030s, tol=3.95e-02  61.75% -> 4.58s
 249: 12.450µm, [5.86e-06;9.97e-01], 1:0.030s, tol=4.02e-02  62.00% -> 4.57s
 250: 12.500µm, [5.85e-06;9.97e-01], 1:0.031s, tol=4.32e-02  62.25% -> 4.64s
 251: 12.550µm, [5.91e-06;9.97e-01], 1:0.030s, tol=4.52e-02  62.50% -> 4.54s
 252: 12.600µm, [6.26e-06;9.97e-01], 1:0.030s, tol=4.57e-02  62.75% -> 4.49s
 253: 12.650µm, [6.37e-06;9.97e-01], 1:0.030s, tol=4.31e-02  63.00% -> 4.46s
 254: 12.700µm, [6.44e-06;9.98e-01], 1:0.030s, tol=4.07e-02  63.25% -> 4.45s
 255: 12.750µm, [6.46e-06;9.98e-01], 1:0.030s, tol=4.15e-02  63.50% -> 4.39s
 256: 12.800µm, [6.47e-06;9.98e-01], 1:0.031s, tol=4.43e-02  63.75% -> 4.43s
 257: 12.850µm, [6.91e-06;9.98e-01], 1:0.031s, tol=4.65e-02  64.00% -> 4.39s
 258: 12.900µm, [7.01e-06;9.98e-01], 1:0.030s, tol=4.45e-02  64.25% -> 4.31s
 259: 12.950µm, [6.99e-06;9.98e-01], 1:0.030s, tol=4.13e-02  64.50% -> 4.26s
 260: 13.000µm, [6.98e-06;9.98e-01], 1:0.031s, tol=4.14e-02  64.75% -> 4.31s
 261: 13.050µm, [6.91e-06;9.98e-01], 1:0.032s, tol=4.41e-02  65.00% -> 4.44s
 262: 13.100µm, [7.11e-06;9.98e-01], 1:0.031s, tol=4.64e-02  65.25% -> 4.32s
 263: 13.150µm, [7.16e-06;9.98e-01], 1:0.031s, tol=4.35e-02  65.50% -> 4.27s
 264: 13.200µm, [7.20e-06;9.98e-01], 1:0.031s, tol=4.08e-02  65.75% -> 4.19s
 265: 13.250µm, [7.21e-06;9.98e-01], 1:0.033s, tol=4.22e-02  66.00% -> 4.43s
 266: 13.300µm, [7.25e-06;9.98e-01], 1:0.031s, tol=4.48e-02  66.25% -> 4.19s
 267: 13.350µm, [7.50e-06;9.98e-01], 1:0.030s, tol=4.48e-02  66.50% -> 4.02s
 268: 13.400µm, [7.52e-06;9.98e-01], 1:0.034s, tol=4.09e-02  66.75% -> 4.55s
 269: 13.450µm, [7.59e-06;9.98e-01], 1:0.030s, tol=4.04e-02  67.00% -> 4.01s
 270: 13.500µm, [7.57e-06;9.98e-01], 1:0.030s, tol=4.28e-02  67.25% -> 3.94s
 271: 13.550µm, [7.89e-06;9.98e-01], 1:0.030s, tol=4.45e-02  67.50% -> 3.93s
 272: 13.600µm, [7.94e-06;9.98e-01], 1:0.030s, tol=4.04e-02  67.75% -> 3.90s
 273: 13.650µm, [8.07e-06;9.99e-01], 1:0.030s, tol=3.91e-02  68.00% -> 3.86s
 274: 13.700µm, [8.09e-06;9.99e-01], 1:0.030s, tol=4.15e-02  68.25% -> 3.81s
 275: 13.750µm, [8.47e-06;9.99e-01], 1:0.030s, tol=4.34e-02  68.50% -> 3.79s
 276: 13.800µm, [8.62e-06;9.99e-01], 1:0.030s, tol=3.95e-02  68.75% -> 3.81s
 277: 13.850µm, [8.76e-06;9.99e-01], 1:0.031s, tol=3.84e-02  69.00% -> 3.83s
 278: 13.900µm, [8.91e-06;9.99e-01], 1:0.031s, tol=4.11e-02  69.25% -> 3.85s
 279: 13.950µm, [9.37e-06;9.99e-01], 1:0.031s, tol=4.18e-02  69.50% -> 3.79s
 280: 14.000µm, [9.56e-06;9.99e-01], 1:0.031s, tol=3.76e-02  69.75% -> 3.80s
 281: 14.050µm, [9.70e-06;9.99e-01], 1:0.031s, tol=3.77e-02  70.00% -> 3.75s
 282: 14.100µm, [1.02e-05;9.99e-01], 1:0.032s, tol=4.09e-02  70.25% -> 3.80s
 283: 14.150µm, [1.05e-05;9.99e-01], 1:0.032s, tol=3.84e-02  70.50% -> 3.77s
 284: 14.200µm, [1.08e-05;9.99e-01], 1:0.032s, tol=3.55e-02  70.75% -> 3.69s
 285: 14.250µm, [1.10e-05;9.99e-01], 1:0.032s, tol=3.73e-02  71.00% -> 3.68s
 286: 14.300µm, [1.18e-05;9.99e-01], 1:0.032s, tol=3.94e-02  71.25% -> 3.68s
 287: 14.350µm, [1.21e-05;9.99e-01], 1:0.032s, tol=3.55e-02  71.50% -> 3.61s
 288: 14.400µm, [1.24e-05;9.99e-01], 1:0.031s, tol=3.51e-02  71.75% -> 3.51s
 289: 14.450µm, [1.32e-05;9.99e-01], 1:0.033s, tol=3.82e-02  72.00% -> 3.68s
 290: 14.500µm, [1.38e-05;9.99e-01], 1:0.032s, tol=3.54e-02  72.25% -> 3.51s
 291: 14.550µm, [1.41e-05;9.99e-01], 1:0.031s, tol=3.32e-02  72.50% -> 3.42s
 292: 14.600µm, [1.49e-05;9.99e-01], 1:0.032s, tol=3.57e-02  72.75% -> 3.44s
 293: 14.650µm, [1.57e-05;9.99e-01], 1:0.031s, tol=3.55e-02  73.00% -> 3.39s
 294: 14.700µm, [1.62e-05;9.99e-01], 1:0.032s, tol=3.20e-02  73.25% -> 3.41s
 295: 14.750µm, [1.68e-05;9.99e-01], 1:0.032s, tol=3.34e-02  73.50% -> 3.38s
 296: 14.800µm, [1.78e-05;9.99e-01], 1:0.032s, tol=3.49e-02  73.75% -> 3.32s
 297: 14.850µm, [1.83e-05;9.99e-01], 1:0.031s, tol=3.16e-02  74.00% -> 3.27s
 298: 14.900µm, [1.89e-05;9.99e-01], 1:0.031s, tol=3.20e-02  74.25% -> 3.21s
 299: 14.950µm, [1.95e-05;9.99e-01], 1:0.031s, tol=3.36e-02  74.50% -> 3.18s
 300: 15.000µm, [2.00e-05;9.99e-01], 1:0.031s, tol=3.07e-02  74.75% -> 3.12s
 301: 15.050µm, [2.04e-05;9.99e-01], 1:0.031s, tol=3.08e-02  75.00% -> 3.08s
 302: 15.100µm, [2.08e-05;9.99e-01], 1:0.031s, tol=3.21e-02  75.25% -> 3.06s
 303: 15.150µm, [2.11e-05;9.99e-01], 1:0.031s, tol=2.95e-02  75.50% -> 3.02s
 304: 15.200µm, [2.14e-05;9.99e-01], 1:0.031s, tol=2.96e-02  75.75% -> 3.00s
 305: 15.250µm, [2.18e-05;9.99e-01], 1:0.031s, tol=3.05e-02  76.00% -> 2.98s
 306: 15.300µm, [2.20e-05;9.99e-01], 1:0.031s, tol=2.84e-02  76.25% -> 2.92s
 307: 15.350µm, [2.23e-05;9.99e-01], 1:0.032s, tol=2.87e-02  76.50% -> 2.97s
 308: 15.400µm, [2.26e-05;9.99e-01], 1:0.032s, tol=2.87e-02  76.75% -> 2.95s
 309: 15.450µm, [2.28e-05;9.99e-01], 1:0.032s, tol=2.73e-02  77.00% -> 2.90s
 310: 15.500µm, [2.30e-05;9.99e-01], 1:0.032s, tol=2.78e-02  77.25% -> 2.88s
 311: 15.550µm, [2.32e-05;9.99e-01], 1:0.031s, tol=2.69e-02  77.50% -> 2.76s
 312: 15.600µm, [2.33e-05;9.99e-01], 1:0.030s, tol=2.61e-02  77.75% -> 2.69s
 313: 15.650µm, [2.34e-05;9.99e-01], 1:0.030s, tol=2.67e-02  78.00% -> 2.67s
 314: 15.700µm, [2.35e-05;9.99e-01], 1:0.030s, tol=2.52e-02  78.25% -> 2.62s
 315: 15.750µm, [2.36e-05;9.99e-01], 1:0.030s, tol=2.48e-02  78.50% -> 2.61s
 316: 15.800µm, [2.37e-05;9.99e-01], 1:0.031s, tol=2.54e-02  78.75% -> 2.59s
 317: 15.850µm, [2.37e-05;9.99e-01], 1:0.030s, tol=2.39e-02  79.00% -> 2.56s
 318: 15.900µm, [2.38e-05;9.99e-01], 1:0.031s, tol=2.32e-02  79.25% -> 2.54s
 319: 15.950µm, [2.39e-05;9.99e-01], 1:0.030s, tol=2.41e-02  79.50% -> 2.47s
 320: 16.000µm, [2.39e-05;9.99e-01], 1:0.030s, tol=2.23e-02  79.75% -> 2.43s
 321: 16.050µm, [2.39e-05;9.99e-01], 1:0.030s, tol=2.20e-02  80.00% -> 2.39s
 322: 16.100µm, [2.40e-05;9.99e-01], 1:0.030s, tol=2.28e-02  80.25% -> 2.41s
 323: 16.150µm, [2.40e-05;9.99e-01], 1:0.030s, tol=2.10e-02  80.50% -> 2.36s
 324: 16.200µm, [2.40e-05;9.99e-01], 1:0.031s, tol=2.09e-02  80.75% -> 2.40s
 325: 16.250µm, [2.40e-05;9.99e-01], 1:0.031s, tol=2.12e-02  81.00% -> 2.32s
 326: 16.300µm, [2.40e-05;9.99e-01], 1:0.031s, tol=1.93e-02  81.25% -> 2.30s
 327: 16.350µm, [2.40e-05;9.99e-01], 1:0.031s, tol=1.93e-02  81.50% -> 2.30s
 328: 16.400µm, [2.41e-05;9.99e-01], 1:0.030s, tol=1.98e-02  81.75% -> 2.20s
 329: 16.450µm, [2.41e-05;9.99e-01], 1:0.030s, tol=1.81e-02  82.00% -> 2.16s
 330: 16.500µm, [2.41e-05;9.99e-01], 1:0.030s, tol=1.66e-02  82.25% -> 2.15s
 331: 16.550µm, [2.41e-05;9.99e-01], 1:0.030s, tol=1.73e-02  82.50% -> 2.12s
 332: 16.600µm, [2.41e-05;9.99e-01], 1:0.031s, tol=1.72e-02  82.75% -> 2.14s
 333: 16.650µm, [2.41e-05;9.99e-01], 1:0.031s, tol=1.66e-02  83.00% -> 2.09s
 334: 16.700µm, [2.41e-05;9.99e-01], 1:0.030s, tol=1.63e-02  83.25% -> 2.02s
 335: 16.750µm, [2.41e-05;9.99e-01], 1:0.032s, tol=1.71e-02  83.50% -> 2.13s
 336: 16.800µm, [2.41e-05;9.99e-01], 1:0.031s, tol=1.88e-02  83.75% -> 1.99s
 337: 16.850µm, [2.41e-05;9.99e-01], 1:0.030s, tol=2.21e-02  84.00% -> 1.92s
 338: 16.900µm, [2.42e-05;9.99e-01], 1:0.030s, tol=2.61e-02  84.25% -> 1.90s
 339: 16.950µm, [2.42e-05;9.99e-01], 1:0.030s, tol=3.05e-02  84.50% -> 1.86s
 340: 17.000µm, [2.42e-05;1.00e+00], 1:0.030s, tol=3.18e-02  84.75% -> 1.85s
 341: 17.050µm, [2.42e-05;1.00e+00], 1:0.030s, tol=3.25e-02  85.00% -> 1.81s
 342: 17.100µm, [2.42e-05;1.00e+00], 1:0.031s, tol=3.29e-02  85.25% -> 1.83s
 343: 17.150µm, [2.43e-05;1.00e+00], 1:0.031s, tol=2.48e-02  85.50% -> 1.80s
 344: 17.200µm, [2.43e-05;1.01e+00], 1:0.031s, tol=1.33e-02  85.75% -> 1.74s
 345: 17.250µm, [2.43e-05;1.01e+00], 1:0.030s, tol=8.31e-03  86.00% -> 1.68s
 346: 17.300µm, [2.43e-05;1.01e+00], 1:0.030s, tol=6.99e-03  86.25% -> 1.66s
 347: 17.350µm, [2.43e-05;1.01e+00], 1:0.031s, tol=6.52e-03  86.50% -> 1.67s
 348: 17.400µm, [2.43e-05;1.01e+00], 1:0.031s, tol=6.25e-03  86.75% -> 1.62s
 349: 17.450µm, [2.43e-05;1.01e+00], 1:0.030s, tol=6.08e-03  87.00% -> 1.58s
 350: 17.500µm, [2.43e-05;1.01e+00], 1:0.030s, tol=5.96e-03  87.25% -> 1.53s
 351: 17.550µm, [2.43e-05;1.01e+00], 1:0.031s, tol=5.88e-03  87.50% -> 1.55s
 352: 17.600µm, [2.43e-05;1.01e+00], 1:0.031s, tol=5.83e-03  87.75% -> 1.52s
 353: 17.650µm, [2.43e-05;1.01e+00], 1:0.030s, tol=5.79e-03  88.00% -> 1.46s
 354: 17.700µm, [2.43e-05;1.01e+00], 1:0.031s, tol=5.76e-03  88.25% -> 1.46s
 355: 17.750µm, [2.43e-05;1.01e+00], 1:0.031s, tol=5.74e-03  88.50% -> 1.43s
 356: 17.800µm, [2.43e-05;1.01e+00], 1:0.032s, tol=5.71e-03  88.75% -> 1.42s
 357: 17.850µm, [2.43e-05;1.01e+00], 1:0.031s, tol=5.69e-03  89.00% -> 1.36s
 358: 17.900µm, [2.43e-05;1.01e+00], 1:0.031s, tol=5.67e-03  89.25% -> 1.31s
 359: 17.950µm, [2.43e-05;1.01e+00], 1:0.031s, tol=5.65e-03  89.50% -> 1.29s
 360: 18.000µm, [2.43e-05;1.01e+00], 1:0.031s, tol=5.64e-03  89.75% -> 1.26s
 361: 18.050µm, [2.43e-05;1.01e+00], 1:0.030s, tol=5.62e-03  90.00% -> 1.20s
 362: 18.100µm, [2.43e-05;1.01e+00], 1:0.030s, tol=5.60e-03  90.25% -> 1.18s
 363: 18.150µm, [2.43e-05;1.01e+00], 1:0.031s, tol=5.59e-03  90.50% -> 1.16s
 364: 18.200µm, [2.43e-05;1.01e+00], 1:0.031s, tol=5.58e-03  90.75% -> 1.13s
 365: 18.250µm, [2.43e-05;1.01e+00], 1:0.030s, tol=5.57e-03  91.00% -> 1.10s
 366: 18.300µm, [2.43e-05;1.01e+00], 1:0.031s, tol=5.55e-03  91.25% -> 1.08s
 367: 18.350µm, [2.43e-05;1.01e+00], 1:0.030s, tol=5.54e-03  91.50% -> 1.03s
 368: 18.400µm, [2.43e-05;1.01e+00], 1:0.030s, tol=5.52e-03  91.75% -> 990.23ms
Saved mesh data in:
/home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/stable/examples/PhaseField/results/Shear2D/Test/Miehe_AT1_optimMesh/Meshes
Saved simulation and summary in:
/home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/stable/examples/PhaseField/results/Shear2D/Test/Miehe_AT1_optimMesh

Generate movie 01/21 (4.76 %) 2.92 s
Generate movie 02/21 (9.52 %) 2.61 s
Generate movie 03/21 (14.29 %) 2.47 s
Generate movie 04/21 (19.05 %) 2.33 s
Generate movie 05/21 (23.81 %) 2.23 s
Generate movie 06/21 (28.57 %) 2.08 s
Generate movie 07/21 (33.33 %) 1.92 s
Generate movie 08/21 (38.10 %) 1.79 s
Generate movie 09/21 (42.86 %) 1.64 s
Generate movie 10/21 (47.62 %) 1.51 s
Generate movie 11/21 (52.38 %) 1.37 s
Generate movie 12/21 (57.14 %) 1.25 s
Generate movie 13/21 (61.90 %) 1.10 s
Generate movie 14/21 (66.67 %) 967.50 ms
Generate movie 15/21 (71.43 %) 825.28 ms
Generate movie 16/21 (76.19 %) 682.78 ms
Generate movie 17/21 (80.95 %) 575.21 ms
Generate movie 18/21 (85.71 %) 406.58 ms
Generate movie 19/21 (90.48 %) 275.01 ms
Generate movie 20/21 (95.24 %) 137.58 ms
Generate movie 21/21 (100.00 %) 0.00 µs
Mesh: 55.156 ms
Boundary Conditions: 6.541 ms
Matrix: 6.566 s
Split: 1.623 s
Solver: 5.165 s
PostProcessing: 334.306 ms
Resolution phase field: 11.193 s
Display: 583.565 ms
PyVista_Interface: 4.629 s


 14 import matplotlib.pyplot as plt
 15 import numpy as np
 16
 17 from EasyFEA import (
 18     Display,
 19     Folder,
 20     Models,
 21     Tic,
 22     ElemType,
 23     Simulations,
 24     Paraview,
 25     PyVista,
 26 )
 27 from EasyFEA.Geoms import Point, Points, Domain, Line, Contour
 28
 29 if __name__ == "__main__":
 30
 31     Display.Clear()
 32
 33     # ----------------------------------------------
 34     # Configuration
 35     # ----------------------------------------------
 36     dim = 2
 37
 38     # simu options
 39     doSimu = True
 40     meshTest = True
 41     openCrack = True
 42     optimMesh = True
 43
 44     # outputs
 45     folder = Folder.Results_Dir() + f"{dim}D"
 46     plotMesh = False
 47     plotEnergy = False
 48     makeParaview = False
 49     makeMovie = True
 50
 51     # phasefield
 52     maxIter = 1000
 53
 54     tolConv = 1e-0  # 1e-1, 1e-2, 1e-3
 55
 56     pfmSolver = Models.PhaseField.SolverType.History
 57
 58     # Available splits: Bourdin, Amor, Miehe, Stress (isotropic)
 59     #                   He, AnisotStrain, AnisotStress, Zhang (anisotropic)
 60     split = Models.PhaseField.SplitType.Miehe
 61
 62     # Available regus: AT1, AT2
 63     regu = Models.PhaseField.ReguType.AT1
 64
 65     # ----------------------------------------------
 66     # Geometry
 67     # ----------------------------------------------
 68     L = 1e-3  # m
 69     l0 = 1e-5  # m
 70     thickness = 1 if dim == 2 else 0.1 / 1000
 71
 72     # ----------------------------------------------
 73     # Material
 74     # ----------------------------------------------
 75     E = 210e9  # Pa
 76     v = 0.3
 77     Gc = 2.7e3  # J/m2
 78
 79     folder_save = Simulations.PhaseField.Folder(
 80         folder,
 81         "",
 82         split,
 83         regu,
 84         "",
 85         tolConv,
 86         pfmSolver,
 87         meshTest,
 88         optimMesh,
 89         not openCrack,
 90     )
 91     Display.MyPrint(folder_save, "green", end="\n")
 92
 93     if doSimu:
 94         # ----------------------------------------------
 95         # Mesh
 96         # ----------------------------------------------
 97         clC = l0 * 2 if meshTest else l0 / 2
 98         if optimMesh:
 99             clD = clC * 4
100             gap = L * 0.05
101             h = L if split == Models.PhaseField.SplitType.Bourdin else L / 2 + gap
102             refineDomain = Domain(Point(L / 2 - gap, 0), Point(L, h, thickness), clC)
103         else:
104             clD = clC
105             refineDomain = None
106
107         pt1 = Point()
108         pt2 = Point(L)
109         pt3 = Point(L, L)
110         pt4 = Point(0, L)
111         contour = Points([pt1, pt2, pt3, pt4], clD)
112
113         if dim == 2:
114             ptC1 = Point(0, L / 2, isOpen=openCrack)
115             ptC2 = Point(L / 2, L / 2)
116             cracks = [Line(ptC1, ptC2, clC, isOpen=openCrack)]
117         elif dim == 3:
118             meshSize = clD if optimMesh else clC
119             ptC1 = Point(0, L / 2, 0, isOpen=openCrack)
120             ptC2 = Point(L / 2, L / 2, 0)
121             ptC3 = Point(L / 2, L / 2, thickness)
122             ptC4 = Point(0, L / 2, thickness, isOpen=openCrack)
123             l1 = Line(ptC1, ptC2, meshSize, openCrack)
124             l2 = Line(ptC2, ptC3, meshSize, False)
125             l3 = Line(ptC3, ptC4, meshSize, openCrack)
126             l4 = Line(ptC4, ptC1, meshSize, openCrack)
127             cracks = [Contour([l1, l2, l3, l4], isOpen=openCrack)]
128
129         if dim == 2:
130             mesh = contour.Mesh_2D([], ElemType.TRI3, cracks, [refineDomain])
131         elif dim == 3:
132             mesh = contour.Mesh_Extrude(
133                 [],
134                 [0, 0, thickness],
135                 [4],
136                 ElemType.PRISM6,
137                 cracks,
138                 [refineDomain],
139                 additionalLines=[l1],
140             )
141
142         # Nodes
143         nodes_crack = mesh.Nodes_Conditions(lambda x, y, z: (y == L / 2) & (x <= L / 2))
144         nodes_upper = mesh.Nodes_Conditions(lambda x, y, z: y == L)
145         nodes_lower = mesh.Nodes_Conditions(lambda x, y, z: y == 0)
146         nodes_left = mesh.Nodes_Conditions(lambda x, y, z: (x == 0) & (y > 0) & (y < L))
147         nodes_right = mesh.Nodes_Conditions(
148             lambda x, y, z: (x == L) & (y > 0) & (y < L)
149         )
150         nodes_edges = np.concatenate([nodes_lower, nodes_right, nodes_upper])
151
152         # ----------------------------------------------
153         # Material
154         # ----------------------------------------------
155         material = Models.Elastic.Isotropic(
156             dim, E=E, v=v, planeStress=False, thickness=thickness
157         )
158         pfm = Models.PhaseField(material, split, regu, Gc, l0, pfmSolver)
159
160         # ----------------------------------------------
161         # Boundary conditions
162         # ----------------------------------------------
163         u_inc = 5e-8 if meshTest else 1e-8
164         N = 400 if meshTest else 2000
165         loadings = np.linspace(u_inc, u_inc * N, N, endpoint=True)
166
167         config = f"""
168         u_inc = {u_inc:.1e};  N = {N}
169         loadings = np.linspace(u_inc, u_inc*N, N, endpoint=True)
170
171         for iter, dep in enumerate(loadings):
172
173             if not openCrack:
174                 simu.add_dirichlet(nodes_crack, [1], ["d"], problemType="damage")
175             simu.add_dirichlet(nodes_left, [0], ["y"])
176             simu.add_dirichlet(nodes_right, [0], ["y"])
177             simu.add_dirichlet(nodes_upper, [dep, 0], ["x", "y"])
178             simu.add_dirichlet(nodes_lower, [0]*dim, simu.Get_unknowns())
179         """
180
181         def Loading(dep):
182             simu.Bc_Init()
183             if not openCrack:
184                 simu.add_dirichlet(nodes_crack, [1], ["d"], problemType="damage")
185             simu.add_dirichlet(nodes_left, [0], ["y"])
186             simu.add_dirichlet(nodes_right, [0], ["y"])
187             simu.add_dirichlet(nodes_upper, [dep, 0], ["x", "y"])
188             simu.add_dirichlet(nodes_lower, [0] * dim, simu.Get_unknowns())
189
190         # ----------------------------------------------
191         # Simulation
192         # ----------------------------------------------
193         simu = Simulations.PhaseField(mesh, pfm, folder=folder_save)
194         simu.Results_Set_Bc_Summary(config)
195
196         dofsX_upper = simu.Bc_dofs_nodes(nodes_upper, ["x"])
197
198         N = len(loadings)
199         nDetect = 0
200         list_dep = []
201         list_f = []
202         for iter, dep in enumerate(loadings):
203             Loading(dep)
204
205             u, _, converg = simu.Solve(tolConv, maxIter, convOption=2)
206             simu.Save_Iter()
207
208             simu.Results_Set_Iteration_Summary(iter, dep * 1e6, "µm", iter / N, True)
209
210             if not converg:
211                 break
212
213             f = np.sum(simu.Calc_Reaction(dofsX_upper, problemType="elastic"))
214             list_dep.append(dep)
215             list_f.append(f)
216
217             if simu.Detect_Damage(nodes_edges, 1):
218                 nDetect += 1
219                 if nDetect == 10:
220                     break
221
222         # ----------------------------------------------
223         # Saving
224         # ----------------------------------------------
225         print()
226         Simulations.Save_pickle((list_f, list_dep), folder_save, "force-displacement")
227         simu.Save(folder_save)
228
229     else:
230         simu: Simulations.PhaseField = Simulations.Load_Simu(folder_save)
231         list_f, list_dep = Simulations.Load_pickle(folder_save, "force-displacement")
232
233     # ----------------------------------------------
234     # Results
235     # ----------------------------------------------
236     Display.Plot_Result(
237         simu,
238         "damage",
239         nodeValues=True,
240         plotMesh=False,
241         folder=folder_save,
242         filename="damage",
243         ncolors=25,
244     )
245     Display.Plot_Mesh(simu)
246     Display.Plot_Iter_Summary(simu, folder_save, None, None)
247     Display.Plot_BoundaryConditions(simu)
248
249     ax = Display.Init_Axes()
250     ax.plot(np.abs(list_dep) * 1e6, np.abs(list_f) * 1e-6, c="blue")
251     ax.set_xlabel("ud [µm]")
252     ax.set_ylabel("f [kN/mm]")
253     ax.grid()
254     Display.Save_fig(folder_save, "force-displacement")
255
256     if plotMesh:
257         ax = Display.Plot_Mesh(simu.mesh, lw=0.3, facecolors="white")
258         ax.axis("off")
259         ax.set_title("")
260         Display.Save_fig(folder_save, "mesh", transparent=True)
261
262     if makeParaview:
263         Paraview.Save_simu(simu, folder_save, 400)
264
265     if makeMovie:
266         simu.Set_Iter(-1)
267         deformFactor = L * 0.05 / simu.Result("displacement_norm").max()
268
269         iterations = np.arange(0, simu.Niter, simu.Niter // 20)
270
271         def Func(plotter, iter):
272             simu.Set_Iter(iterations[iter])
273             thresh = PyVista._pvMesh(simu, "damage", deformFactor).threshold((0, 0.8))
274             PyVista.Plot(thresh, "damage", plotMesh=True, plotter=plotter, clim=(0, 1))
275
276         PyVista.Movie_func(Func, iterations.size, folder_save, "damage.gif")
277
278     if plotEnergy:
279         Display.Plot_Energy(simu, N=400, folder=folder_save)
280
281     Tic.Resume()
282
283     if doSimu:
284         Tic.Plot_History(folder_save, False)
285
286     plt.show()

Total running time of the script: (0 minutes 18.583 seconds)

Gallery generated by Sphinx-Gallery