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

Movie_func 0/20
Movie_func 1/20 (5.00 %) 2.20 s
Movie_func 2/20 (10.00 %) 2.13 s
Movie_func 3/20 (15.00 %) 2.01 s
Movie_func 4/20 (20.00 %) 1.86 s
Movie_func 5/20 (25.00 %) 1.73 s
Movie_func 6/20 (30.00 %) 1.63 s
Movie_func 7/20 (35.00 %) 1.48 s
Movie_func 8/20 (40.00 %) 1.38 s
Movie_func 9/20 (45.00 %) 1.28 s
Movie_func 10/20 (50.00 %) 1.18 s
Movie_func 11/20 (55.00 %) 1.05 s
Movie_func 12/20 (60.00 %) 919.61 ms
Movie_func 13/20 (65.00 %) 818.03 ms
Movie_func 14/20 (70.00 %) 688.12 ms
Movie_func 15/20 (75.00 %) 571.72 ms
Movie_func 16/20 (80.00 %) 459.39 ms
Movie_func 17/20 (85.00 %) 351.81 ms
Movie_func 18/20 (90.00 %) 228.41 ms
Movie_func 19/20 (95.00 %) 114.55 ms
Movie_func 20/20 (100.00 %) 0.00 µs
PostProcessing : 103.170 ms
Display : 483.521 ms
Mesh : 48.319 ms
Boundary Conditions : 4.857 ms
Matrix : 2.280 s
Split : 1.277 s
Solver : 5.101 s
Resolution phase field : 9.565 s
PyVista_Interface : 3.984 s


 14 from EasyFEA import (
 15     Display,
 16     Folder,
 17     Models,
 18     plt,
 19     np,
 20     Tic,
 21     ElemType,
 22     Mesh,
 23     Simulations,
 24     Paraview,
 25     PyVista,
 26 )
 27 from EasyFEA.Geoms import Point, Points, Domain, Line, Contour
 28
 29 import multiprocessing
 30
 31 # Display.Clear()
 32
 33 useParallel = False
 34 nProcs = 4  # number of processes in parallel
 35
 36 # ----------------------------------------------
 37 # Configuration
 38 # ----------------------------------------------
 39 dim = 2
 40
 41 # simu options
 42 doSimu = True
 43 meshTest = True
 44 openCrack = True
 45 optimMesh = True
 46
 47 # outputs
 48 plotResult = True
 49 showResult = True
 50 plotMesh = False
 51 plotEnergy = False
 52 saveParaview = False
 53 makeMovie = True
 54
 55 # phasefield
 56 maxIter = 1000
 57
 58 tolConv = 1e-0  # 1e-1, 1e-2, 1e-3
 59
 60 pfmSolver = Models.PhaseField.SolverType.History
 61
 62 # splits = ["Bourdin","Amor","Miehe","Stress"] # Splits Isotropes
 63 # splits = ["He","AnisotStrain","AnisotStress","Zhang"] # Splits Anisotropes sans bourdin
 64 # splits = ["Bourdin","Amor","Miehe","Stress","He","AnisotStrain","AnisotStress","Zhang"]
 65 splits = ["Miehe"]
 66
 67 # regus = ["AT1", "AT2"]
 68 regus = ["AT1"]  # "AT1", "AT2"
 69
 70 # ----------------------------------------------
 71 # Mesh
 72 # ----------------------------------------------
 73 L = 1e-3
 74 # m
 75 l0 = 1e-5
 76 thickness = 1 if dim == 2 else 0.1 / 1000
 77
 78
 79 def DoMesh(split: str) -> Mesh:
 80     # meshSize
 81     clC = l0 * 2 if meshTest else l0 / 2
 82     if optimMesh:
 83         # a coarser mesh can be used outside the refined zone
 84         clD = clC * 4
 85         # refines the mesh in the area where the crack will propagate
 86         gap = L * 0.05
 87         h = L if split == "Bourdin" else L / 2 + gap
 88         refineDomain = Domain(Point(L / 2 - gap, 0), Point(L, h, thickness), clC)
 89     else:
 90         clD = clC
 91         refineDomain = None
 92
 93     # geom
 94     pt1 = Point()
 95     pt2 = Point(L)
 96     pt3 = Point(L, L)
 97     pt4 = Point(0, L)
 98     contour = Points([pt1, pt2, pt3, pt4], clD)
 99
100     if dim == 2:
101         ptC1 = Point(0, L / 2, isOpen=openCrack)
102         ptC2 = Point(L / 2, L / 2)
103         cracks = [Line(ptC1, ptC2, clC, isOpen=openCrack)]
104     if dim == 3:
105         meshSize = clD if optimMesh else clC
106         ptC1 = Point(0, L / 2, 0, isOpen=openCrack)
107         ptC2 = Point(L / 2, L / 2, 0)
108         ptC3 = Point(L / 2, L / 2, thickness)
109         ptC4 = Point(0, L / 2, thickness, isOpen=openCrack)
110         l1 = Line(ptC1, ptC2, meshSize, openCrack)
111         l2 = Line(ptC2, ptC3, meshSize, False)
112         l3 = Line(ptC3, ptC4, meshSize, openCrack)
113         l4 = Line(ptC4, ptC1, meshSize, openCrack)
114         cracks = [Contour([l1, l2, l3, l4], isOpen=openCrack)]
115
116     # folder = Folder.Join("",results=True)
117     # ax = Display.Init_Axes()
118     # contour.Get_Contour().Plot(ax, color='k', plotPoints=False)
119     # cracks[0].Plot(ax, color='k', plotPoints=False)
120     # # if refineDomain != None:
121     # #     refineDomain.Plot(ax, color='k', plotPoints=False)
122     # ax.axis('off')
123     # Display.Save_fig(folder,"sample",True)
124
125     if dim == 2:
126         mesh = contour.Mesh_2D([], ElemType.TRI3, cracks, [refineDomain])
127     elif dim == 3:
128         mesh = contour.Mesh_Extrude(
129             [],
130             [0, 0, thickness],
131             [4],
132             ElemType.PRISM6,
133             cracks,
134             [refineDomain],
135             additionalLines=[l1],
136         )
137
138     return mesh
139
140
141 # ----------------------------------------------
142 # Simu
143 # ----------------------------------------------
144 def DoSimu(split: str, regu: str):
145     # Builds the path to the folder based on the problem data
146     folder_save = Folder.PhaseField_Folder(
147         f"Shear_{dim}D",
148         "ElasIsot",
149         split,
150         regu,
151         "DP",
152         tolConv,
153         pfmSolver,
154         meshTest,
155         optimMesh,
156         not openCrack,
157     )
158
159     Display.MyPrint(folder_save, "green")
160
161     if doSimu:
162         mesh = DoMesh(split)
163
164         # Nodes recovery
165         nodes_crack = mesh.Nodes_Conditions(lambda x, y, z: (y == L / 2) & (x <= L / 2))
166         nodes_upper = mesh.Nodes_Conditions(lambda x, y, z: y == L)
167         nodes_lower = mesh.Nodes_Conditions(lambda x, y, z: y == 0)
168         nodes_left = mesh.Nodes_Conditions(lambda x, y, z: (x == 0) & (y > 0) & (y < L))
169         nodes_right = mesh.Nodes_Conditions(
170             lambda x, y, z: (x == L) & (y > 0) & (y < L)
171         )
172
173         # Builds edge nodes
174         nodes_edges = []
175         for nodes in [nodes_lower, nodes_right, nodes_upper]:
176             nodes_edges.extend(nodes)
177
178         # ----------------------------------------------
179         # Material
180         # ----------------------------------------------
181         material = Models.ElasIsot(
182             dim, E=210e9, v=0.3, planeStress=False, thickness=thickness
183         )
184         Gc = 2.7e3  # J/m2
185         pfm = Models.PhaseField(material, split, regu, Gc, l0, pfmSolver)
186
187         # ----------------------------------------------
188         # Boundary conditions
189         # ----------------------------------------------
190         u_inc = 5e-8 if meshTest else 1e-8
191         N = 400 if meshTest else 2000
192
193         loadings = np.linspace(u_inc, u_inc * N, N, endpoint=True)
194
195         config = f"""
196         u_inc = {u_inc:.1e}
197         N = {N}
198
199         for iter, dep in enumerate(loadings):
200
201         loadings = np.linspace(u_inc, u_inc*N, N, endpoint=True)
202
203         if not openCrack:
204             simu.add_dirichlet(nodes_crack, [1], ["d"], problemType="damage")
205         simu.add_dirichlet(nodes_left, [0], ["y"])
206         simu.add_dirichlet(nodes_right, [0], ["y"])
207         simu.add_dirichlet(nodes_upper, [dep,0], ["x","y"])
208         simu.add_dirichlet(nodes_lower, [0]*dim, simu.Get_dofs())
209         """
210
211         def Loading(dep):
212             """Boundary conditions"""
213             simu.Bc_Init()
214             if not openCrack:
215                 simu.add_dirichlet(nodes_crack, [1], ["d"], problemType="damage")
216             simu.add_dirichlet(nodes_left, [0], ["y"])
217             simu.add_dirichlet(nodes_right, [0], ["y"])
218             simu.add_dirichlet(nodes_upper, [dep, 0], ["x", "y"])
219             simu.add_dirichlet(nodes_lower, [0] * dim, simu.Get_unknowns())
220
221         # ----------------------------------------------
222         # Simulation
223         # ----------------------------------------------
224         simu = Simulations.PhaseFieldSimu(mesh, pfm, verbosity=False)
225         simu.Results_Set_Bc_Summary(config)
226
227         dofsX_upper = simu.Bc_dofs_nodes(nodes_upper, ["x"])
228
229         # INIT
230         N = len(loadings)
231         nDetect = 0
232         displacement = []
233         force = []
234         for iter, dep in enumerate(loadings):
235             # apply new boundary conditions
236             Loading(dep)
237
238             # solve and save iter
239             u, _, Kglob, converg = simu.Solve(tolConv, maxIter, convOption=2)
240             simu.Save_Iter()
241
242             # print iter solution
243             simu.Results_Set_Iteration_Summary(iter, dep * 1e6, "µm", iter / N, True)
244
245             # If the solver has not converged, stop the simulation.
246             if not converg:
247                 break
248
249             # resulting force on upper edge
250             f = np.sum(Kglob[dofsX_upper, :] @ u)
251
252             displacement.append(dep)
253             force.append(f)
254
255             # Detect damaged edges
256             if np.any(simu.damage[nodes_edges] >= 1):
257                 nDetect += 1
258                 if nDetect == 10:
259                     break
260
261         # ----------------------------------------------
262         # Saving
263         # ----------------------------------------------
264         force = np.asarray(force)
265         displacement = np.asarray(displacement)
266         print()
267         Simulations.Save_pickle(
268             (force, displacement), folder_save, "force-displacement"
269         )
270         simu.Save(folder_save)
271
272     else:
273         simu: Simulations.PhaseFieldSimu = Simulations.Load_Simu(folder_save)
274         force, displacement = Simulations.Load_pickle(folder_save, "force-displacement")
275
276     # ----------------------------------------------
277     # Results
278     # ---------------------------------------------
279     if plotResult:
280         ax = Display.Plot_Result(
281             simu,
282             "damage",
283             nodeValues=True,
284             plotMesh=False,
285             folder=folder_save,
286             filename="damage",
287             ncolors=25,
288         )
289         Display.Plot_Mesh(simu)
290         Display.Plot_Iter_Summary(simu, folder_save, None, None)
291         Display.Plot_BoundaryConditions(simu)
292         Display.Plot_Force_Displacement(
293             force * 1e-6, displacement * 1e6, "ud [µm]", "f [kN/mm]", folder_save
294         )
295
296         # ax = Display.Plot_Result(simu, "damage", 1.5, ncolors=21, clim=(0,0.9))
297         # ax.axis('off'); ax.set_title("")
298         # Display.Save_fig(folder, "deform damage")
299
300     if plotMesh:
301         # PyVista.Plot_Mesh(simu.mesh).show()
302         ax = Display.Plot_Mesh(simu.mesh, lw=0.3, facecolors="white")
303         ax.axis("off")
304         ax.set_title("")
305         Display.Save_fig(folder_save, "mesh", transparent=True)
306
307     if saveParaview:
308         Paraview.Save_simu(simu, folder_save, 400)
309
310     if makeMovie:
311         simu.Set_Iter(-1)
312         nodes_upper = simu.mesh.Nodes_Conditions(lambda x, y, z: y == L)
313         depMax = simu.Result("displacement_norm")[nodes_upper].max()
314         deformFactor = L * 0.05 / depMax
315
316         iterations = np.arange(0, simu.Niter, simu.Niter // 20)
317
318         def Func(plotter, iter):
319             simu.Set_Iter(iterations[iter])
320
321             grid = PyVista._pyVistaMesh(simu, "damage", deformFactor)
322
323             tresh = grid.threshold((0, 0.8))
324
325             PyVista.Plot(
326                 tresh,
327                 "damage",
328                 deformFactor,
329                 plotMesh=True,
330                 plotter=plotter,
331                 clim=(0, 1),
332             )
333
334         PyVista.Movie_func(Func, iterations.size, folder_save, "damage.gif")
335
336     if plotEnergy:
337         Display.Plot_Energy(simu, N=400, folder=folder_save)
338
339     Tic.Resume()
340
341     if doSimu:
342         Tic.Plot_History(folder_save, False)
343
344     if showResult:
345         plt.show()
346
347     # plt.close("all")
348     Tic.Clear()
349
350
351 if __name__ == "__main__":
352     # generates configs
353     Splits = []
354     Regus = []
355     for split in splits.copy():
356         for regu in regus.copy():
357             Splits.append(split)
358             Regus.append(regu)
359
360     if useParallel:
361         items = [(split, regu) for split, regu in zip(Splits, Regus)]
362         with multiprocessing.Pool(nProcs) as pool:
363             for result in pool.starmap(DoSimu, items):
364                 pass
365     else:
366         [DoSimu(split, regu) for split, regu in zip(Splits, Regus)]

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

Gallery generated by Sphinx-Gallery