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

Generate movie 01/21 (4.76 %) 2.07 s
Generate movie 02/21 (9.52 %) 1.68 s
Generate movie 03/21 (14.29 %) 1.58 s
Generate movie 04/21 (19.05 %) 1.48 s
Generate movie 05/21 (23.81 %) 1.43 s
Generate movie 06/21 (28.57 %) 1.34 s
Generate movie 07/21 (33.33 %) 1.20 s
Generate movie 08/21 (38.10 %) 1.12 s
Generate movie 09/21 (42.86 %) 1.03 s
Generate movie 10/21 (47.62 %) 956.09 ms
Generate movie 11/21 (52.38 %) 859.83 ms
Generate movie 12/21 (57.14 %) 777.95 ms
Generate movie 13/21 (61.90 %) 684.49 ms
Generate movie 14/21 (66.67 %) 598.05 ms
Generate movie 15/21 (71.43 %) 517.03 ms
Generate movie 16/21 (76.19 %) 431.13 ms
Generate movie 17/21 (80.95 %) 347.56 ms
Generate movie 18/21 (85.71 %) 268.68 ms
Generate movie 19/21 (90.48 %) 175.04 ms
Generate movie 20/21 (95.24 %) 85.91 ms
Generate movie 21/21 (100.00 %) 0.00 µs
Mesh : 45.805 ms
Boundary Conditions : 4.706 ms
Matrix : 5.263 s
Split : 1.429 s
Solver : 8.931 s
PostProcessing : 96.072 ms
Resolution phase field : 13.494 s
Display : 423.666 ms
PyVista_Interface : 2.889 s


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

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

Gallery generated by Sphinx-Gallery