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.7.0/examples/PhaseField/results/Shear2D/Test/ElasIsot_Miehe_AT1_DP_optimMesh   1 : 0.050 µm, [0.00e+00; 0.00e+00], 1:0.073 s, tol=1.00e+00
   2 : 0.100 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=1.00e+00  0.25 % -> 24.13 s
   3 : 0.150 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=5.56e-01  0.50 % -> 24.38 s
   4 : 0.200 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=4.38e-01  0.75 % -> 24.31 s
   5 : 0.250 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=3.60e-01  1.00 % -> 23.96 s
   6 : 0.300 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=3.06e-01  1.25 % -> 24.47 s
   7 : 0.350 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=2.65e-01  1.50 % -> 24.26 s
   8 : 0.400 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.34e-01  1.75 % -> 23.92 s
   9 : 0.450 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.10e-01  2.00 % -> 24.03 s
  10 : 0.500 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.90e-01  2.25 % -> 23.89 s
  11 : 0.550 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=1.74e-01  2.50 % -> 23.32 s
  12 : 0.600 µm, [0.00e+00; 0.00e+00], 1:0.059 s, tol=1.60e-01  2.75 % -> 23.04 s
  13 : 0.650 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=1.48e-01  3.00 % -> 23.26 s
  14 : 0.700 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.38e-01  3.25 % -> 23.72 s
  15 : 0.750 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=1.29e-01  3.50 % -> 23.86 s
  16 : 0.800 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=1.21e-01  3.75 % -> 23.78 s
  17 : 0.850 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.14e-01  4.00 % -> 23.53 s
  18 : 0.900 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.08e-01  4.25 % -> 23.27 s
  19 : 0.950 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=1.02e-01  4.50 % -> 22.94 s
  20 : 1.000 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=9.75e-02  4.75 % -> 23.09 s
  21 : 1.050 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=9.30e-02  5.00 % -> 22.79 s
  22 : 1.100 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=8.88e-02  5.25 % -> 22.63 s
  23 : 1.150 µm, [0.00e+00; 0.00e+00], 1:0.059 s, tol=8.51e-02  5.50 % -> 22.26 s
  24 : 1.200 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=8.16e-02  5.75 % -> 22.69 s
  25 : 1.250 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=7.84e-02  6.00 % -> 22.90 s
  26 : 1.300 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=7.54e-02  6.25 % -> 22.92 s
  27 : 1.350 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=7.27e-02  6.50 % -> 22.84 s
  28 : 1.400 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=7.02e-02  6.75 % -> 22.91 s
  29 : 1.450 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=6.78e-02  7.00 % -> 22.93 s
  30 : 1.500 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=6.56e-02  7.25 % -> 22.90 s
  31 : 1.550 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=6.35e-02  7.50 % -> 22.75 s
  32 : 1.600 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=6.15e-02  7.75 % -> 23.04 s
  33 : 1.650 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=5.97e-02  8.00 % -> 22.46 s
  34 : 1.700 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=5.80e-02  8.25 % -> 22.58 s
  35 : 1.750 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=5.63e-02  8.50 % -> 22.43 s
  36 : 1.800 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=5.48e-02  8.75 % -> 22.62 s
  37 : 1.850 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=5.33e-02  9.00 % -> 22.59 s
  38 : 1.900 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=5.19e-02  9.25 % -> 22.51 s
  39 : 1.950 µm, [0.00e+00; 0.00e+00], 1:0.063 s, tol=5.06e-02  9.50 % -> 22.66 s
  40 : 2.000 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=4.94e-02  9.75 % -> 22.17 s
  41 : 2.050 µm, [0.00e+00; 0.00e+00], 1:0.063 s, tol=4.82e-02  10.00 % -> 22.70 s
  42 : 2.100 µm, [0.00e+00; 0.00e+00], 1:0.063 s, tol=4.71e-02  10.25 % -> 22.51 s
  43 : 2.150 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=4.60e-02  10.50 % -> 22.21 s
  44 : 2.200 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=4.49e-02  10.75 % -> 21.84 s
  45 : 2.250 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=4.40e-02  11.00 % -> 21.82 s
  46 : 2.300 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=4.30e-02  11.25 % -> 22.10 s
  47 : 2.350 µm, [0.00e+00; 0.00e+00], 1:0.063 s, tol=4.21e-02  11.50 % -> 22.21 s
  48 : 2.400 µm, [0.00e+00; 0.00e+00], 1:0.063 s, tol=4.12e-02  11.75 % -> 22.07 s
  49 : 2.450 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=4.04e-02  12.00 % -> 21.95 s
  50 : 2.500 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=3.96e-02  12.25 % -> 21.62 s
  51 : 2.550 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=3.88e-02  12.50 % -> 21.51 s
  52 : 2.600 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=3.81e-02  12.75 % -> 21.71 s
  53 : 2.650 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=3.74e-02  13.00 % -> 21.66 s
  54 : 2.700 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=3.67e-02  13.25 % -> 21.41 s
  55 : 2.750 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=3.60e-02  13.50 % -> 21.21 s
  56 : 2.800 µm, [0.00e+00; 0.00e+00], 1:0.063 s, tol=3.54e-02  13.75 % -> 21.86 s
  57 : 2.850 µm, [0.00e+00; 0.00e+00], 1:0.063 s, tol=3.48e-02  14.00 % -> 21.62 s
  58 : 2.900 µm, [0.00e+00; 0.00e+00], 1:0.066 s, tol=3.42e-02  14.25 % -> 22.48 s
  59 : 2.950 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=3.36e-02  14.50 % -> 20.90 s
  60 : 3.000 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=3.31e-02  14.75 % -> 21.06 s
  61 : 3.050 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=3.25e-02  15.00 % -> 21.12 s
  62 : 3.100 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=3.20e-02  15.25 % -> 20.64 s
  63 : 3.150 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=3.15e-02  15.50 % -> 20.66 s
  64 : 3.200 µm, [0.00e+00; 0.00e+00], 1:0.063 s, tol=3.10e-02  15.75 % -> 21.13 s
  65 : 3.250 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=3.05e-02  16.00 % -> 20.60 s
  66 : 3.300 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=3.01e-02  16.25 % -> 20.43 s
  67 : 3.350 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.96e-02  16.50 % -> 20.52 s
  68 : 3.400 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=2.92e-02  16.75 % -> 20.54 s
  69 : 3.450 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.88e-02  17.00 % -> 20.27 s
  70 : 3.500 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.84e-02  17.25 % -> 20.09 s
  71 : 3.550 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.80e-02  17.50 % -> 20.07 s
  72 : 3.600 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.76e-02  17.75 % -> 20.08 s
  73 : 3.650 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=2.72e-02  18.00 % -> 20.34 s
  74 : 3.700 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.68e-02  18.25 % -> 19.96 s
  75 : 3.750 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.65e-02  18.50 % -> 19.46 s
  76 : 3.800 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.61e-02  18.75 % -> 19.49 s
  77 : 3.850 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.58e-02  19.00 % -> 19.55 s
  78 : 3.900 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.55e-02  19.25 % -> 19.82 s
  79 : 3.950 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.52e-02  19.50 % -> 19.70 s
  80 : 4.000 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=2.48e-02  19.75 % -> 19.81 s
  81 : 4.050 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.45e-02  20.00 % -> 19.41 s
  82 : 4.100 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.42e-02  20.25 % -> 19.50 s
  83 : 4.150 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.40e-02  20.50 % -> 19.07 s
  84 : 4.200 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.37e-02  20.75 % -> 19.19 s
  85 : 4.250 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.34e-02  21.00 % -> 19.09 s
  86 : 4.300 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.31e-02  21.25 % -> 19.06 s
  87 : 4.350 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.29e-02  21.50 % -> 19.05 s
  88 : 4.400 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.26e-02  21.75 % -> 18.98 s
  89 : 4.450 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.23e-02  22.00 % -> 18.78 s
  90 : 4.500 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.21e-02  22.25 % -> 18.62 s
  91 : 4.550 µm, [0.00e+00; 0.00e+00], 1:0.059 s, tol=2.19e-02  22.50 % -> 18.41 s
  92 : 4.600 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.16e-02  22.75 % -> 18.56 s
  93 : 4.650 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.14e-02  23.00 % -> 18.50 s
  94 : 4.700 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.12e-02  23.25 % -> 18.54 s
  95 : 4.750 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.09e-02  23.50 % -> 18.61 s
  96 : 4.800 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=2.07e-02  23.75 % -> 18.59 s
  97 : 4.850 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.05e-02  24.00 % -> 18.33 s
  98 : 4.900 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.03e-02  24.25 % -> 18.15 s
  99 : 4.950 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=2.01e-02  24.50 % -> 18.04 s
 100 : 5.000 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=1.99e-02  24.75 % -> 18.03 s
 101 : 5.050 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=1.97e-02  25.00 % -> 17.96 s
 102 : 5.100 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=1.95e-02  25.25 % -> 18.40 s
 103 : 5.150 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.93e-02  25.50 % -> 18.30 s
 104 : 5.200 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=1.91e-02  25.75 % -> 18.30 s
 105 : 5.250 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=1.90e-02  26.00 % -> 18.41 s
 106 : 5.300 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=1.88e-02  26.25 % -> 18.30 s
 107 : 5.350 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=1.86e-02  26.50 % -> 18.09 s
 108 : 5.400 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.84e-02  26.75 % -> 17.89 s
 109 : 5.450 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.83e-02  27.00 % -> 17.86 s
 110 : 5.500 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.81e-02  27.25 % -> 17.82 s
 111 : 5.550 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=1.79e-02  27.50 % -> 17.93 s
 112 : 5.600 µm, [0.00e+00; 0.00e+00], 1:0.063 s, tol=1.78e-02  27.75 % -> 18.06 s
 113 : 5.650 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=1.76e-02  28.00 % -> 17.76 s
 114 : 5.700 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.75e-02  28.25 % -> 17.49 s
 115 : 5.750 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.73e-02  28.50 % -> 17.34 s
 116 : 5.800 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.72e-02  28.75 % -> 17.42 s
 117 : 5.850 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=1.70e-02  29.00 % -> 16.97 s
 118 : 5.900 µm, [0.00e+00; 0.00e+00], 1:0.065 s, tol=1.69e-02  29.25 % -> 18.37 s
 119 : 5.950 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.67e-02  29.50 % -> 17.07 s
 120 : 6.000 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.66e-02  29.75 % -> 17.13 s
 121 : 6.050 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.65e-02  30.00 % -> 17.02 s
 122 : 6.100 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=1.63e-02  30.25 % -> 16.77 s
 123 : 6.150 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=1.62e-02  30.50 % -> 17.10 s
 124 : 6.200 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=1.61e-02  30.75 % -> 17.05 s
 125 : 6.250 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.59e-02  31.00 % -> 16.91 s
 126 : 6.300 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.58e-02  31.25 % -> 16.81 s
 127 : 6.350 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.57e-02  31.50 % -> 16.60 s
 128 : 6.400 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.56e-02  31.75 % -> 16.60 s
 129 : 6.450 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=1.54e-02  32.00 % -> 16.36 s
 130 : 6.500 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=1.53e-02  32.25 % -> 16.26 s
 131 : 6.550 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.52e-02  32.50 % -> 16.51 s
 132 : 6.600 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.51e-02  32.75 % -> 16.31 s
 133 : 6.650 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.50e-02  33.00 % -> 16.46 s
 134 : 6.700 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.49e-02  33.25 % -> 16.31 s
 135 : 6.750 µm, [0.00e+00; 0.00e+00], 1:0.060 s, tol=1.48e-02  33.50 % -> 15.98 s
 136 : 6.800 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.47e-02  33.75 % -> 16.20 s
 137 : 6.850 µm, [0.00e+00; 0.00e+00], 1:0.061 s, tol=1.45e-02  34.00 % -> 16.22 s
 138 : 6.900 µm, [0.00e+00; 0.00e+00], 1:0.063 s, tol=1.44e-02  34.25 % -> 16.52 s
 139 : 6.950 µm, [0.00e+00; 0.00e+00], 1:0.062 s, tol=1.43e-02  34.50 % -> 16.33 s
 140 : 7.000 µm, [6.00e-07; 1.59e-03], 1:0.062 s, tol=1.42e-02  34.75 % -> 16.07 s
 141 : 7.050 µm, [2.22e-06; 6.20e-03], 1:0.060 s, tol=1.39e-02  35.00 % -> 15.67 s
 142 : 7.100 µm, [4.16e-06; 1.23e-02], 1:0.061 s, tol=1.38e-02  35.25 % -> 15.81 s
 143 : 7.150 µm, [6.13e-06; 1.92e-02], 1:0.063 s, tol=1.37e-02  35.50 % -> 16.16 s
 144 : 7.200 µm, [7.97e-06; 2.65e-02], 1:0.062 s, tol=1.35e-02  35.75 % -> 15.84 s
 145 : 7.250 µm, [9.65e-06; 3.41e-02], 1:0.062 s, tol=1.35e-02  36.00 % -> 15.88 s
 146 : 7.300 µm, [1.12e-05; 4.18e-02], 1:0.062 s, tol=1.34e-02  36.25 % -> 15.79 s
 147 : 7.350 µm, [1.25e-05; 4.96e-02], 1:0.062 s, tol=1.33e-02  36.50 % -> 15.75 s
 148 : 7.400 µm, [1.36e-05; 5.75e-02], 1:0.062 s, tol=1.32e-02  36.75 % -> 15.79 s
 149 : 7.450 µm, [1.46e-05; 6.55e-02], 1:0.061 s, tol=1.31e-02  37.00 % -> 15.42 s
 150 : 7.500 µm, [1.55e-05; 7.35e-02], 1:0.060 s, tol=1.30e-02  37.25 % -> 15.15 s
 151 : 7.550 µm, [1.62e-05; 8.16e-02], 1:0.060 s, tol=1.29e-02  37.50 % -> 15.12 s
 152 : 7.600 µm, [1.68e-05; 8.98e-02], 1:0.060 s, tol=1.28e-02  37.75 % -> 15.02 s
 153 : 7.650 µm, [1.72e-05; 9.80e-02], 1:0.061 s, tol=1.28e-02  38.00 % -> 15.06 s
 154 : 7.700 µm, [1.76e-05; 1.06e-01], 1:0.061 s, tol=1.27e-02  38.25 % -> 14.95 s
 155 : 7.750 µm, [1.79e-05; 1.15e-01], 1:0.060 s, tol=1.26e-02  38.50 % -> 14.76 s
 156 : 7.800 µm, [1.81e-05; 1.23e-01], 1:0.061 s, tol=1.25e-02  38.75 % -> 14.99 s
 157 : 7.850 µm, [1.82e-05; 1.32e-01], 1:0.060 s, tol=1.24e-02  39.00 % -> 14.65 s
 158 : 7.900 µm, [1.83e-05; 1.40e-01], 1:0.061 s, tol=1.24e-02  39.25 % -> 14.75 s
 159 : 7.950 µm, [1.83e-05; 1.49e-01], 1:0.061 s, tol=1.23e-02  39.50 % -> 14.80 s
 160 : 8.000 µm, [1.82e-05; 1.58e-01], 1:0.061 s, tol=1.22e-02  39.75 % -> 14.71 s
 161 : 8.050 µm, [1.81e-05; 1.67e-01], 1:0.060 s, tol=1.21e-02  40.00 % -> 14.43 s
 162 : 8.100 µm, [1.80e-05; 1.76e-01], 1:0.060 s, tol=1.21e-02  40.25 % -> 14.45 s
 163 : 8.150 µm, [1.78e-05; 1.85e-01], 1:0.061 s, tol=1.20e-02  40.50 % -> 14.43 s
 164 : 8.200 µm, [1.75e-05; 1.94e-01], 1:0.060 s, tol=1.19e-02  40.75 % -> 14.33 s
 165 : 8.250 µm, [1.73e-05; 2.03e-01], 1:0.060 s, tol=1.18e-02  41.00 % -> 14.12 s
 166 : 8.300 µm, [1.70e-05; 2.12e-01], 1:0.061 s, tol=1.18e-02  41.25 % -> 14.38 s
 167 : 8.350 µm, [1.67e-05; 2.22e-01], 1:0.061 s, tol=1.17e-02  41.50 % -> 14.30 s
 168 : 8.400 µm, [1.64e-05; 2.31e-01], 1:0.061 s, tol=1.16e-02  41.75 % -> 14.32 s
 169 : 8.450 µm, [1.60e-05; 2.41e-01], 1:0.061 s, tol=1.16e-02  42.00 % -> 14.14 s
 170 : 8.500 µm, [1.57e-05; 2.51e-01], 1:0.062 s, tol=1.15e-02  42.25 % -> 14.27 s
 171 : 8.550 µm, [1.53e-05; 2.60e-01], 1:0.061 s, tol=1.14e-02  42.50 % -> 14.13 s
 172 : 8.600 µm, [1.49e-05; 2.70e-01], 1:0.060 s, tol=1.14e-02  42.75 % -> 13.80 s
 173 : 8.650 µm, [1.45e-05; 2.80e-01], 1:0.061 s, tol=1.13e-02  43.00 % -> 13.82 s
 174 : 8.700 µm, [1.42e-05; 2.90e-01], 1:0.060 s, tol=1.13e-02  43.25 % -> 13.72 s
 175 : 8.750 µm, [1.38e-05; 3.01e-01], 1:0.061 s, tol=1.12e-02  43.50 % -> 13.86 s
 176 : 8.800 µm, [1.34e-05; 3.11e-01], 1:0.062 s, tol=1.11e-02  43.75 % -> 13.99 s
 177 : 8.850 µm, [1.30e-05; 3.21e-01], 1:0.062 s, tol=1.11e-02  44.00 % -> 13.91 s
 178 : 8.900 µm, [1.26e-05; 3.32e-01], 1:0.061 s, tol=1.10e-02  44.25 % -> 13.61 s
 179 : 8.950 µm, [1.22e-05; 3.43e-01], 1:0.062 s, tol=1.10e-02  44.50 % -> 13.66 s
 180 : 9.000 µm, [1.19e-05; 3.56e-01], 1:0.062 s, tol=1.09e-02  44.75 % -> 13.70 s
 181 : 9.050 µm, [1.15e-05; 3.71e-01], 1:0.062 s, tol=1.08e-02  45.00 % -> 13.65 s
 182 : 9.100 µm, [1.12e-05; 3.86e-01], 1:0.062 s, tol=1.07e-02  45.25 % -> 13.55 s
 183 : 9.150 µm, [1.09e-05; 4.01e-01], 1:0.062 s, tol=1.07e-02  45.50 % -> 13.49 s
 184 : 9.200 µm, [1.05e-05; 4.17e-01], 1:0.062 s, tol=1.06e-02  45.75 % -> 13.39 s
 185 : 9.250 µm, [1.02e-05; 4.32e-01], 1:0.062 s, tol=1.06e-02  46.00 % -> 13.33 s
 186 : 9.300 µm, [9.84e-06; 4.47e-01], 1:0.062 s, tol=1.06e-02  46.25 % -> 13.39 s
 187 : 9.350 µm, [9.49e-06; 4.62e-01], 1:0.062 s, tol=1.05e-02  46.50 % -> 13.20 s
 188 : 9.400 µm, [9.27e-06; 4.77e-01], 1:0.062 s, tol=1.04e-02  46.75 % -> 13.22 s
 189 : 9.450 µm, [9.10e-06; 4.94e-01], 1:0.062 s, tol=1.03e-02  47.00 % -> 13.07 s
 190 : 9.500 µm, [8.97e-06; 5.14e-01], 1:0.062 s, tol=1.03e-02  47.25 % -> 12.99 s
 191 : 9.550 µm, [8.88e-06; 5.39e-01], 1:0.063 s, tol=1.02e-02  47.50 % -> 13.13 s
 192 : 9.600 µm, [8.83e-06; 5.71e-01], 1:0.062 s, tol=1.03e-02  47.75 % -> 13.03 s
 193 : 9.650 µm, [8.78e-06; 6.08e-01], 1:0.062 s, tol=1.04e-02  48.00 % -> 12.94 s
 194 : 9.700 µm, [8.76e-06; 6.49e-01], 1:0.062 s, tol=1.07e-02  48.25 % -> 12.77 s
 195 : 9.750 µm, [8.75e-06; 6.95e-01], 1:0.062 s, tol=1.13e-02  48.50 % -> 12.70 s
 196 : 9.800 µm, [9.05e-06; 7.43e-01], 1:0.064 s, tol=1.23e-02  48.75 % -> 13.02 s
 197 : 9.850 µm, [9.39e-06; 7.94e-01], 1:0.063 s, tol=1.47e-02  49.00 % -> 12.83 s
 198 : 9.900 µm, [9.76e-06; 8.46e-01], 1:0.063 s, tol=1.92e-02  49.25 % -> 12.77 s
 199 : 9.950 µm, [9.84e-06; 8.94e-01], 1:0.064 s, tol=2.57e-02  49.50 % -> 12.85 s
 200 : 10.000 µm, [9.59e-06; 9.28e-01], 1:0.062 s, tol=3.11e-02  49.75 % -> 12.53 s
 201 : 10.050 µm, [9.14e-06; 9.49e-01], 1:0.063 s, tol=3.15e-02  50.00 % -> 12.60 s
 202 : 10.100 µm, [9.06e-06; 9.58e-01], 1:0.063 s, tol=2.98e-02  50.25 % -> 12.50 s
 203 : 10.150 µm, [9.00e-06; 9.60e-01], 1:0.063 s, tol=3.01e-02  50.50 % -> 12.40 s
 204 : 10.200 µm, [9.21e-06; 9.60e-01], 1:0.062 s, tol=3.37e-02  50.75 % -> 12.30 s
 205 : 10.250 µm, [9.27e-06; 9.60e-01], 1:0.062 s, tol=4.01e-02  51.00 % -> 12.10 s
 206 : 10.300 µm, [9.16e-06; 9.71e-01], 1:0.063 s, tol=4.63e-02  51.25 % -> 12.35 s
 207 : 10.350 µm, [8.87e-06; 9.82e-01], 1:0.064 s, tol=4.81e-02  51.50 % -> 12.48 s
 208 : 10.400 µm, [8.45e-06; 9.89e-01], 1:0.063 s, tol=4.32e-02  51.75 % -> 12.07 s
 209 : 10.450 µm, [8.11e-06; 9.92e-01], 1:0.062 s, tol=3.51e-02  52.00 % -> 11.87 s
 210 : 10.500 µm, [8.09e-06; 9.94e-01], 1:0.063 s, tol=3.08e-02  52.25 % -> 12.05 s
 211 : 10.550 µm, [8.04e-06; 9.95e-01], 1:0.061 s, tol=2.91e-02  52.50 % -> 11.65 s
 212 : 10.600 µm, [7.90e-06; 9.95e-01], 1:0.062 s, tol=2.80e-02  52.75 % -> 11.76 s
 213 : 10.650 µm, [7.80e-06; 9.95e-01], 1:0.061 s, tol=2.75e-02  53.00 % -> 11.48 s
 214 : 10.700 µm, [7.73e-06; 9.94e-01], 1:0.062 s, tol=2.85e-02  53.25 % -> 11.54 s
 215 : 10.750 µm, [7.62e-06; 9.94e-01], 1:0.061 s, tol=3.07e-02  53.50 % -> 11.43 s
 216 : 10.800 µm, [7.45e-06; 9.93e-01], 1:0.061 s, tol=3.36e-02  53.75 % -> 11.35 s
 217 : 10.850 µm, [7.23e-06; 9.93e-01], 1:0.062 s, tol=3.59e-02  54.00 % -> 11.49 s
 218 : 10.900 µm, [6.95e-06; 9.92e-01], 1:0.061 s, tol=3.61e-02  54.25 % -> 11.25 s
 219 : 10.950 µm, [6.71e-06; 9.92e-01], 1:0.061 s, tol=3.49e-02  54.50 % -> 11.16 s
 220 : 11.000 µm, [6.54e-06; 9.92e-01], 1:0.061 s, tol=3.27e-02  54.75 % -> 10.97 s
 221 : 11.050 µm, [6.47e-06; 9.92e-01], 1:0.061 s, tol=3.08e-02  55.00 % -> 10.97 s
 222 : 11.100 µm, [6.45e-06; 9.92e-01], 1:0.062 s, tol=3.00e-02  55.25 % -> 11.06 s
 223 : 11.150 µm, [6.36e-06; 9.93e-01], 1:0.061 s, tol=2.96e-02  55.50 % -> 10.91 s
 224 : 11.200 µm, [6.21e-06; 9.93e-01], 1:0.061 s, tol=2.93e-02  55.75 % -> 10.86 s
 225 : 11.250 µm, [6.16e-06; 9.93e-01], 1:0.061 s, tol=3.12e-02  56.00 % -> 10.78 s
 226 : 11.300 µm, [6.13e-06; 9.93e-01], 1:0.062 s, tol=3.48e-02  56.25 % -> 10.81 s
 227 : 11.350 µm, [6.03e-06; 9.94e-01], 1:0.061 s, tol=3.87e-02  56.50 % -> 10.68 s
 228 : 11.400 µm, [5.87e-06; 9.94e-01], 1:0.062 s, tol=4.11e-02  56.75 % -> 10.69 s
 229 : 11.450 µm, [5.71e-06; 9.94e-01], 1:0.063 s, tol=4.15e-02  57.00 % -> 10.83 s
 230 : 11.500 µm, [5.65e-06; 9.95e-01], 1:0.062 s, tol=4.01e-02  57.25 % -> 10.68 s
 231 : 11.550 µm, [5.70e-06; 9.95e-01], 1:0.063 s, tol=3.84e-02  57.50 % -> 10.79 s
 232 : 11.600 µm, [5.74e-06; 9.95e-01], 1:0.063 s, tol=3.71e-02  57.75 % -> 10.71 s
 233 : 11.650 µm, [5.66e-06; 9.95e-01], 1:0.063 s, tol=3.53e-02  58.00 % -> 10.60 s
 234 : 11.700 µm, [5.61e-06; 9.95e-01], 1:0.064 s, tol=3.48e-02  58.25 % -> 10.67 s
 235 : 11.750 µm, [5.60e-06; 9.95e-01], 1:0.062 s, tol=3.73e-02  58.50 % -> 10.36 s
 236 : 11.800 µm, [5.51e-06; 9.96e-01], 1:0.062 s, tol=4.00e-02  58.75 % -> 10.17 s
 237 : 11.850 µm, [5.41e-06; 9.96e-01], 1:0.062 s, tol=4.25e-02  59.00 % -> 10.19 s
 238 : 11.900 µm, [5.38e-06; 9.96e-01], 1:0.063 s, tol=4.26e-02  59.25 % -> 10.33 s
 239 : 11.950 µm, [5.51e-06; 9.96e-01], 1:0.064 s, tol=4.20e-02  59.50 % -> 10.39 s
 240 : 12.000 µm, [5.60e-06; 9.96e-01], 1:0.063 s, tol=4.03e-02  59.75 % -> 10.22 s
 241 : 12.050 µm, [5.53e-06; 9.96e-01], 1:0.064 s, tol=3.75e-02  60.00 % -> 10.27 s
 242 : 12.100 µm, [5.56e-06; 9.96e-01], 1:0.063 s, tol=3.75e-02  60.25 % -> 10.00 s
 243 : 12.150 µm, [5.52e-06; 9.97e-01], 1:0.063 s, tol=3.99e-02  60.50 % -> 9.92 s
 244 : 12.200 µm, [5.47e-06; 9.97e-01], 1:0.062 s, tol=4.29e-02  60.75 % -> 9.77 s
 245 : 12.250 µm, [5.50e-06; 9.97e-01], 1:0.063 s, tol=4.42e-02  61.00 % -> 9.77 s
 246 : 12.300 µm, [5.73e-06; 9.97e-01], 1:0.062 s, tol=4.43e-02  61.25 % -> 9.68 s
 247 : 12.350 µm, [5.83e-06; 9.97e-01], 1:0.063 s, tol=4.23e-02  61.50 % -> 9.72 s
 248 : 12.400 µm, [5.81e-06; 9.97e-01], 1:0.063 s, tol=3.95e-02  61.75 % -> 9.64 s
 249 : 12.450 µm, [5.86e-06; 9.97e-01], 1:0.062 s, tol=4.02e-02  62.00 % -> 9.49 s
 250 : 12.500 µm, [5.85e-06; 9.97e-01], 1:0.062 s, tol=4.32e-02  62.25 % -> 9.42 s
 251 : 12.550 µm, [5.91e-06; 9.97e-01], 1:0.062 s, tol=4.52e-02  62.50 % -> 9.31 s
 252 : 12.600 µm, [6.26e-06; 9.97e-01], 1:0.062 s, tol=4.57e-02  62.75 % -> 9.24 s
 253 : 12.650 µm, [6.37e-06; 9.97e-01], 1:0.063 s, tol=4.31e-02  63.00 % -> 9.28 s
 254 : 12.700 µm, [6.44e-06; 9.98e-01], 1:0.063 s, tol=4.07e-02  63.25 % -> 9.20 s
 255 : 12.750 µm, [6.46e-06; 9.98e-01], 1:0.063 s, tol=4.15e-02  63.50 % -> 9.25 s
 256 : 12.800 µm, [6.47e-06; 9.98e-01], 1:0.062 s, tol=4.43e-02  63.75 % -> 9.03 s
 257 : 12.850 µm, [6.91e-06; 9.98e-01], 1:0.063 s, tol=4.65e-02  64.00 % -> 9.01 s
 258 : 12.900 µm, [7.01e-06; 9.98e-01], 1:0.063 s, tol=4.45e-02  64.25 % -> 8.98 s
 259 : 12.950 µm, [6.99e-06; 9.98e-01], 1:0.062 s, tol=4.13e-02  64.50 % -> 8.83 s
 260 : 13.000 µm, [6.98e-06; 9.98e-01], 1:0.061 s, tol=4.14e-02  64.75 % -> 8.64 s
 261 : 13.050 µm, [6.91e-06; 9.98e-01], 1:0.062 s, tol=4.41e-02  65.00 % -> 8.65 s
 262 : 13.100 µm, [7.11e-06; 9.98e-01], 1:0.062 s, tol=4.64e-02  65.25 % -> 8.67 s
 263 : 13.150 µm, [7.16e-06; 9.98e-01], 1:0.063 s, tol=4.35e-02  65.50 % -> 8.63 s
 264 : 13.200 µm, [7.20e-06; 9.98e-01], 1:0.064 s, tol=4.08e-02  65.75 % -> 8.73 s
 265 : 13.250 µm, [7.21e-06; 9.98e-01], 1:0.063 s, tol=4.22e-02  66.00 % -> 8.51 s
 266 : 13.300 µm, [7.25e-06; 9.98e-01], 1:0.063 s, tol=4.48e-02  66.25 % -> 8.53 s
 267 : 13.350 µm, [7.50e-06; 9.98e-01], 1:0.063 s, tol=4.48e-02  66.50 % -> 8.44 s
 268 : 13.400 µm, [7.52e-06; 9.98e-01], 1:0.062 s, tol=4.09e-02  66.75 % -> 8.31 s
 269 : 13.450 µm, [7.59e-06; 9.98e-01], 1:0.062 s, tol=4.04e-02  67.00 % -> 8.19 s
 270 : 13.500 µm, [7.57e-06; 9.98e-01], 1:0.064 s, tol=4.28e-02  67.25 % -> 8.40 s
 271 : 13.550 µm, [7.89e-06; 9.98e-01], 1:0.063 s, tol=4.45e-02  67.50 % -> 8.13 s
 272 : 13.600 µm, [7.94e-06; 9.98e-01], 1:0.065 s, tol=4.04e-02  67.75 % -> 8.37 s
 273 : 13.650 µm, [8.07e-06; 9.99e-01], 1:0.065 s, tol=3.91e-02  68.00 % -> 8.33 s
 274 : 13.700 µm, [8.09e-06; 9.99e-01], 1:0.062 s, tol=4.15e-02  68.25 % -> 7.87 s
 275 : 13.750 µm, [8.47e-06; 9.99e-01], 1:0.062 s, tol=4.34e-02  68.50 % -> 7.79 s
 276 : 13.800 µm, [8.62e-06; 9.99e-01], 1:0.061 s, tol=3.95e-02  68.75 % -> 7.64 s
 277 : 13.850 µm, [8.76e-06; 9.99e-01], 1:0.061 s, tol=3.84e-02  69.00 % -> 7.62 s
 278 : 13.900 µm, [8.91e-06; 9.99e-01], 1:0.061 s, tol=4.11e-02  69.25 % -> 7.55 s
 279 : 13.950 µm, [9.37e-06; 9.99e-01], 1:0.063 s, tol=4.18e-02  69.50 % -> 7.63 s
 280 : 14.000 µm, [9.56e-06; 9.99e-01], 1:0.062 s, tol=3.76e-02  69.75 % -> 7.51 s
 281 : 14.050 µm, [9.70e-06; 9.99e-01], 1:0.062 s, tol=3.77e-02  70.00 % -> 7.47 s
 282 : 14.100 µm, [1.02e-05; 9.99e-01], 1:0.062 s, tol=4.09e-02  70.25 % -> 7.35 s
 283 : 14.150 µm, [1.05e-05; 9.99e-01], 1:0.062 s, tol=3.84e-02  70.50 % -> 7.34 s
 284 : 14.200 µm, [1.08e-05; 9.99e-01], 1:0.061 s, tol=3.55e-02  70.75 % -> 7.18 s
 285 : 14.250 µm, [1.10e-05; 9.99e-01], 1:0.062 s, tol=3.73e-02  71.00 % -> 7.16 s
 286 : 14.300 µm, [1.18e-05; 9.99e-01], 1:0.062 s, tol=3.94e-02  71.25 % -> 7.18 s
 287 : 14.350 µm, [1.21e-05; 9.99e-01], 1:0.064 s, tol=3.55e-02  71.50 % -> 7.25 s
 288 : 14.400 µm, [1.24e-05; 9.99e-01], 1:0.062 s, tol=3.51e-02  71.75 % -> 7.04 s
 289 : 14.450 µm, [1.32e-05; 9.99e-01], 1:0.063 s, tol=3.82e-02  72.00 % -> 7.00 s
 290 : 14.500 µm, [1.38e-05; 9.99e-01], 1:0.062 s, tol=3.54e-02  72.25 % -> 6.89 s
 291 : 14.550 µm, [1.41e-05; 9.99e-01], 1:0.063 s, tol=3.32e-02  72.50 % -> 6.94 s
 292 : 14.600 µm, [1.49e-05; 9.99e-01], 1:0.062 s, tol=3.57e-02  72.75 % -> 6.77 s
 293 : 14.650 µm, [1.57e-05; 9.99e-01], 1:0.062 s, tol=3.55e-02  73.00 % -> 6.64 s
 294 : 14.700 µm, [1.62e-05; 9.99e-01], 1:0.061 s, tol=3.20e-02  73.25 % -> 6.50 s
 295 : 14.750 µm, [1.68e-05; 9.99e-01], 1:0.061 s, tol=3.34e-02  73.50 % -> 6.50 s
 296 : 14.800 µm, [1.78e-05; 9.99e-01], 1:0.064 s, tol=3.49e-02  73.75 % -> 6.67 s
 297 : 14.850 µm, [1.83e-05; 9.99e-01], 1:0.082 s, tol=3.16e-02  74.00 % -> 8.55 s
 298 : 14.900 µm, [1.89e-05; 9.99e-01], 1:0.062 s, tol=3.20e-02  74.25 % -> 6.36 s
 299 : 14.950 µm, [1.95e-05; 9.99e-01], 1:0.061 s, tol=3.36e-02  74.50 % -> 6.26 s
 300 : 15.000 µm, [2.00e-05; 9.99e-01], 1:0.061 s, tol=3.07e-02  74.75 % -> 6.20 s
 301 : 15.050 µm, [2.04e-05; 9.99e-01], 1:0.062 s, tol=3.08e-02  75.00 % -> 6.19 s
 302 : 15.100 µm, [2.08e-05; 9.99e-01], 1:0.063 s, tol=3.21e-02  75.25 % -> 6.19 s
 303 : 15.150 µm, [2.11e-05; 9.99e-01], 1:0.063 s, tol=2.95e-02  75.50 % -> 6.13 s
 304 : 15.200 µm, [2.14e-05; 9.99e-01], 1:0.062 s, tol=2.96e-02  75.75 % -> 5.97 s
 305 : 15.250 µm, [2.18e-05; 9.99e-01], 1:0.061 s, tol=3.05e-02  76.00 % -> 5.89 s
 306 : 15.300 µm, [2.20e-05; 9.99e-01], 1:0.061 s, tol=2.84e-02  76.25 % -> 5.83 s
 307 : 15.350 µm, [2.23e-05; 9.99e-01], 1:0.061 s, tol=2.87e-02  76.50 % -> 5.73 s
 308 : 15.400 µm, [2.26e-05; 9.99e-01], 1:0.061 s, tol=2.87e-02  76.75 % -> 5.68 s
 309 : 15.450 µm, [2.28e-05; 9.99e-01], 1:0.061 s, tol=2.73e-02  77.00 % -> 5.59 s
 310 : 15.500 µm, [2.30e-05; 9.99e-01], 1:0.061 s, tol=2.78e-02  77.25 % -> 5.57 s
 311 : 15.550 µm, [2.32e-05; 9.99e-01], 1:0.062 s, tol=2.69e-02  77.50 % -> 5.55 s
 312 : 15.600 µm, [2.33e-05; 9.99e-01], 1:0.061 s, tol=2.61e-02  77.75 % -> 5.47 s
 313 : 15.650 µm, [2.34e-05; 9.99e-01], 1:0.061 s, tol=2.67e-02  78.00 % -> 5.37 s
 314 : 15.700 µm, [2.35e-05; 9.99e-01], 1:0.061 s, tol=2.52e-02  78.25 % -> 5.32 s
 315 : 15.750 µm, [2.36e-05; 9.99e-01], 1:0.061 s, tol=2.48e-02  78.50 % -> 5.24 s
 316 : 15.800 µm, [2.37e-05; 9.99e-01], 1:0.062 s, tol=2.54e-02  78.75 % -> 5.31 s
 317 : 15.850 µm, [2.37e-05; 9.99e-01], 1:0.061 s, tol=2.39e-02  79.00 % -> 5.14 s
 318 : 15.900 µm, [2.38e-05; 9.99e-01], 1:0.062 s, tol=2.32e-02  79.25 % -> 5.11 s
 319 : 15.950 µm, [2.39e-05; 9.99e-01], 1:0.061 s, tol=2.41e-02  79.50 % -> 5.00 s
 320 : 16.000 µm, [2.39e-05; 9.99e-01], 1:0.061 s, tol=2.23e-02  79.75 % -> 4.95 s
 321 : 16.050 µm, [2.39e-05; 9.99e-01], 1:0.061 s, tol=2.20e-02  80.00 % -> 4.87 s
 322 : 16.100 µm, [2.40e-05; 9.99e-01], 1:0.061 s, tol=2.28e-02  80.25 % -> 4.83 s
 323 : 16.150 µm, [2.40e-05; 9.99e-01], 1:0.060 s, tol=2.10e-02  80.50 % -> 4.70 s
 324 : 16.200 µm, [2.40e-05; 9.99e-01], 1:0.061 s, tol=2.09e-02  80.75 % -> 4.69 s
 325 : 16.250 µm, [2.40e-05; 9.99e-01], 1:0.061 s, tol=2.12e-02  81.00 % -> 4.67 s
 326 : 16.300 µm, [2.40e-05; 9.99e-01], 1:0.061 s, tol=1.93e-02  81.25 % -> 4.59 s
 327 : 16.350 µm, [2.40e-05; 9.99e-01], 1:0.061 s, tol=1.93e-02  81.50 % -> 4.55 s
 328 : 16.400 µm, [2.41e-05; 9.99e-01], 1:0.062 s, tol=1.98e-02  81.75 % -> 4.52 s
 329 : 16.450 µm, [2.41e-05; 9.99e-01], 1:0.062 s, tol=1.81e-02  82.00 % -> 4.45 s
 330 : 16.500 µm, [2.41e-05; 9.99e-01], 1:0.062 s, tol=1.66e-02  82.25 % -> 4.38 s
 331 : 16.550 µm, [2.41e-05; 9.99e-01], 1:0.061 s, tol=1.73e-02  82.50 % -> 4.30 s
 332 : 16.600 µm, [2.41e-05; 9.99e-01], 1:0.061 s, tol=1.72e-02  82.75 % -> 4.24 s
 333 : 16.650 µm, [2.41e-05; 9.99e-01], 1:0.064 s, tol=1.66e-02  83.00 % -> 4.36 s
 334 : 16.700 µm, [2.41e-05; 9.99e-01], 1:0.063 s, tol=1.63e-02  83.25 % -> 4.23 s
 335 : 16.750 µm, [2.41e-05; 9.99e-01], 1:0.062 s, tol=1.71e-02  83.50 % -> 4.11 s
 336 : 16.800 µm, [2.41e-05; 9.99e-01], 1:0.062 s, tol=1.88e-02  83.75 % -> 4.00 s
 337 : 16.850 µm, [2.41e-05; 9.99e-01], 1:0.062 s, tol=2.21e-02  84.00 % -> 3.96 s
 338 : 16.900 µm, [2.42e-05; 9.99e-01], 1:0.062 s, tol=2.61e-02  84.25 % -> 3.88 s
 339 : 16.950 µm, [2.42e-05; 9.99e-01], 1:0.061 s, tol=3.05e-02  84.50 % -> 3.78 s
 340 : 17.000 µm, [2.42e-05; 1.00e+00], 1:0.061 s, tol=3.18e-02  84.75 % -> 3.74 s
 341 : 17.050 µm, [2.42e-05; 1.00e+00], 1:0.061 s, tol=3.25e-02  85.00 % -> 3.65 s
 342 : 17.100 µm, [2.42e-05; 1.00e+00], 1:0.062 s, tol=3.29e-02  85.25 % -> 3.69 s
 343 : 17.150 µm, [2.43e-05; 1.00e+00], 1:0.063 s, tol=2.48e-02  85.50 % -> 3.66 s
 344 : 17.200 µm, [2.43e-05; 1.01e+00], 1:0.063 s, tol=1.33e-02  85.75 % -> 3.59 s
 345 : 17.250 µm, [2.43e-05; 1.01e+00], 1:0.062 s, tol=8.31e-03  86.00 % -> 3.47 s
 346 : 17.300 µm, [2.43e-05; 1.01e+00], 1:0.062 s, tol=6.99e-03  86.25 % -> 3.41 s
 347 : 17.350 µm, [2.43e-05; 1.01e+00], 1:0.063 s, tol=6.52e-03  86.50 % -> 3.39 s
 348 : 17.400 µm, [2.43e-05; 1.01e+00], 1:0.063 s, tol=6.25e-03  86.75 % -> 3.32 s
 349 : 17.450 µm, [2.43e-05; 1.01e+00], 1:0.064 s, tol=6.08e-03  87.00 % -> 3.32 s
 350 : 17.500 µm, [2.43e-05; 1.01e+00], 1:0.063 s, tol=5.96e-03  87.25 % -> 3.22 s
 351 : 17.550 µm, [2.43e-05; 1.01e+00], 1:0.063 s, tol=5.88e-03  87.50 % -> 3.14 s
 352 : 17.600 µm, [2.43e-05; 1.01e+00], 1:0.063 s, tol=5.83e-03  87.75 % -> 3.08 s
 353 : 17.650 µm, [2.43e-05; 1.01e+00], 1:0.062 s, tol=5.79e-03  88.00 % -> 2.95 s
 354 : 17.700 µm, [2.43e-05; 1.01e+00], 1:0.062 s, tol=5.76e-03  88.25 % -> 2.92 s
 355 : 17.750 µm, [2.43e-05; 1.01e+00], 1:0.062 s, tol=5.74e-03  88.50 % -> 2.87 s
 356 : 17.800 µm, [2.43e-05; 1.01e+00], 1:0.061 s, tol=5.71e-03  88.75 % -> 2.74 s
 357 : 17.850 µm, [2.43e-05; 1.01e+00], 1:0.061 s, tol=5.69e-03  89.00 % -> 2.70 s
 358 : 17.900 µm, [2.43e-05; 1.01e+00], 1:0.062 s, tol=5.67e-03  89.25 % -> 2.65 s
 359 : 17.950 µm, [2.43e-05; 1.01e+00], 1:0.063 s, tol=5.65e-03  89.50 % -> 2.65 s
 360 : 18.000 µm, [2.43e-05; 1.01e+00], 1:0.063 s, tol=5.64e-03  89.75 % -> 2.58 s
 361 : 18.050 µm, [2.43e-05; 1.01e+00], 1:0.062 s, tol=5.62e-03  90.00 % -> 2.47 s
 362 : 18.100 µm, [2.43e-05; 1.01e+00], 1:0.062 s, tol=5.60e-03  90.25 % -> 2.41 s
 363 : 18.150 µm, [2.43e-05; 1.01e+00], 1:0.063 s, tol=5.59e-03  90.50 % -> 2.38 s
 364 : 18.200 µm, [2.43e-05; 1.01e+00], 1:0.062 s, tol=5.58e-03  90.75 % -> 2.30 s
 365 : 18.250 µm, [2.43e-05; 1.01e+00], 1:0.062 s, tol=5.57e-03  91.00 % -> 2.24 s
 366 : 18.300 µm, [2.43e-05; 1.01e+00], 1:0.062 s, tol=5.55e-03  91.25 % -> 2.18 s
 367 : 18.350 µm, [2.43e-05; 1.01e+00], 1:0.062 s, tol=5.54e-03  91.50 % -> 2.10 s
 368 : 18.400 µm, [2.43e-05; 1.01e+00], 1:0.063 s, tol=5.52e-03  91.75 % -> 2.09 s
Saved:
/home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/v1.7.0/examples/PhaseField/results/Shear2D/Test/ElasIsot_Miehe_AT1_DP_optimMesh/force-displacement.pickle
Saved:
/home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/v1.7.0/examples/PhaseField/results/Shear2D/Test/ElasIsot_Miehe_AT1_DP_optimMesh/simulation.pickle
Saved:
/home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/v1.7.0/examples/PhaseField/results/Shear2D/Test/ElasIsot_Miehe_AT1_DP_optimMesh/summary.txt

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

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

Gallery generated by Sphinx-Gallery