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

Generate movie 01/21 (4.76 %) 2.09 s
Generate movie 02/21 (9.52 %) 1.67 s
Generate movie 03/21 (14.29 %) 1.55 s
Generate movie 04/21 (19.05 %) 1.47 s
Generate movie 05/21 (23.81 %) 1.41 s
Generate movie 06/21 (28.57 %) 1.29 s
Generate movie 07/21 (33.33 %) 1.21 s
Generate movie 08/21 (38.10 %) 1.14 s
Generate movie 09/21 (42.86 %) 1.05 s
Generate movie 10/21 (47.62 %) 987.96 ms
Generate movie 11/21 (52.38 %) 872.30 ms
Generate movie 12/21 (57.14 %) 800.42 ms
Generate movie 13/21 (61.90 %) 705.62 ms
Generate movie 14/21 (66.67 %) 607.35 ms
Generate movie 15/21 (71.43 %) 516.65 ms
Generate movie 16/21 (76.19 %) 434.23 ms
Generate movie 17/21 (80.95 %) 345.66 ms
Generate movie 18/21 (85.71 %) 260.99 ms
Generate movie 19/21 (90.48 %) 177.74 ms
Generate movie 20/21 (95.24 %) 91.87 ms
Generate movie 21/21 (100.00 %) 0.00 µs
Mesh : 45.504 ms
Boundary Conditions : 4.525 ms
Matrix : 5.019 s
Split : 1.410 s
Solver : 8.830 s
PostProcessing : 247.411 ms
Resolution phase field : 13.350 s
Display : 432.959 ms
PyVista_Interface : 2.913 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         E = {E:.2e} Pa;  v = {v};  Gc = {Gc:.2e} J/m2;  l0 = {l0:.2e} m
169
170         u_inc = {u_inc:.1e};  N = {N}
171         loadings = np.linspace(u_inc, u_inc*N, N, endpoint=True)
172
173         for iter, dep in enumerate(loadings):
174
175         if not openCrack:
176             simu.add_dirichlet(nodes_crack, [1], ["d"], problemType="damage")
177         simu.add_dirichlet(nodes_left, [0], ["y"])
178         simu.add_dirichlet(nodes_right, [0], ["y"])
179         simu.add_dirichlet(nodes_upper, [dep, 0], ["x", "y"])
180         simu.add_dirichlet(nodes_lower, [0]*dim, simu.Get_unknowns())
181         """
182
183         def Loading(dep):
184             simu.Bc_Init()
185             if not openCrack:
186                 simu.add_dirichlet(nodes_crack, [1], ["d"], problemType="damage")
187             simu.add_dirichlet(nodes_left, [0], ["y"])
188             simu.add_dirichlet(nodes_right, [0], ["y"])
189             simu.add_dirichlet(nodes_upper, [dep, 0], ["x", "y"])
190             simu.add_dirichlet(nodes_lower, [0] * dim, simu.Get_unknowns())
191
192         # ----------------------------------------------
193         # Simulation
194         # ----------------------------------------------
195         simu = Simulations.PhaseField(mesh, pfm, verbosity=False)
196         simu.Results_Set_Bc_Summary(config)
197
198         dofsX_upper = simu.Bc_dofs_nodes(nodes_upper, ["x"])
199
200         N = len(loadings)
201         nDetect = 0
202         list_dep = []
203         list_f = []
204         for iter, dep in enumerate(loadings):
205             Loading(dep)
206
207             u, _, converg = simu.Solve(tolConv, maxIter, convOption=2)
208             simu.Save_Iter()
209
210             simu.Results_Set_Iteration_Summary(iter, dep * 1e6, "µm", iter / N, True)
211
212             if not converg:
213                 break
214
215             f = np.sum(simu.Calc_Reaction(dofsX_upper, problemType="elastic"))
216             list_dep.append(dep)
217             list_f.append(f)
218
219             if simu.Detect_Damage(nodes_edges, 1):
220                 nDetect += 1
221                 if nDetect == 10:
222                     break
223
224         # ----------------------------------------------
225         # Saving
226         # ----------------------------------------------
227         print()
228         Simulations.Save_pickle((list_f, list_dep), folder_save, "force-displacement")
229         simu.Save(folder_save)
230
231     else:
232         simu: Simulations.PhaseField = Simulations.Load_Simu(folder_save)
233         list_f, list_dep = Simulations.Load_pickle(folder_save, "force-displacement")
234
235     # ----------------------------------------------
236     # Results
237     # ----------------------------------------------
238     Display.Plot_Result(
239         simu,
240         "damage",
241         nodeValues=True,
242         plotMesh=False,
243         folder=folder_save,
244         filename="damage",
245         ncolors=25,
246     )
247     Display.Plot_Mesh(simu)
248     Display.Plot_Iter_Summary(simu, folder_save, None, None)
249     Display.Plot_BoundaryConditions(simu)
250
251     ax = Display.Init_Axes()
252     ax.plot(np.abs(list_dep) * 1e6, np.abs(list_f) * 1e-6, c="blue")
253     ax.set_xlabel("ud [µm]")
254     ax.set_ylabel("f [kN/mm]")
255     ax.grid()
256     Display.Save_fig(folder_save, "force-displacement")
257
258     if plotMesh:
259         ax = Display.Plot_Mesh(simu.mesh, lw=0.3, facecolors="white")
260         ax.axis("off")
261         ax.set_title("")
262         Display.Save_fig(folder_save, "mesh", transparent=True)
263
264     if makeParaview:
265         Paraview.Save_simu(simu, folder_save, 400)
266
267     if makeMovie:
268         simu.Set_Iter(-1)
269         deformFactor = L * 0.05 / simu.Result("displacement_norm").max()
270
271         iterations = np.arange(0, simu.Niter, simu.Niter // 20)
272
273         def Func(plotter, iter):
274             simu.Set_Iter(iterations[iter])
275             thresh = PyVista._pvMesh(simu, "damage", deformFactor).threshold((0, 0.8))
276             PyVista.Plot(thresh, "damage", plotMesh=True, plotter=plotter, clim=(0, 1))
277
278         PyVista.Movie_func(Func, iterations.size, folder_save, "damage.gif")
279
280     if plotEnergy:
281         Display.Plot_Energy(simu, N=400, folder=folder_save)
282
283     Tic.Resume()
284
285     if doSimu:
286         Tic.Plot_History(folder_save, False)
287
288     plt.show()

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

Gallery generated by Sphinx-Gallery