Note

Go to the end to download the full example code.

Tension#

Damage simulation for a plate subjected to tension.

Tension
  • $\phi$
  • TRI3: Ne = 1518, Nn = 819
  • Tension
  • Boundary conditions
  • Tension
  • Summary
/home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/v1.7.1/examples/PhaseField/results/Tension2D/Test/Elas_Isot_Bourdin_AT1_DP_optimMesh   1 : 0.000 µm, [0.00e+00; 0.00e+00], 1:0.035 s, tol=0.00e+00
   2 : 0.100 µm, [0.00e+00; 0.00e+00], 1:0.023 s, tol=0.00e+00
   3 : 0.200 µm, [0.00e+00; 0.00e+00], 1:0.044 s, tol=0.00e+00
   4 : 0.300 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
   5 : 0.400 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
   6 : 0.500 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
   7 : 0.600 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
   8 : 0.700 µm, [0.00e+00; 0.00e+00], 1:0.052 s, tol=0.00e+00
   9 : 0.800 µm, [0.00e+00; 0.00e+00], 1:0.052 s, tol=0.00e+00
  10 : 0.900 µm, [0.00e+00; 0.00e+00], 1:0.058 s, tol=0.00e+00
  11 : 1.000 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  12 : 1.100 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
  13 : 1.200 µm, [0.00e+00; 0.00e+00], 1:0.054 s, tol=0.00e+00
  14 : 1.300 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  15 : 1.400 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  16 : 1.500 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  17 : 1.600 µm, [0.00e+00; 0.00e+00], 1:0.059 s, tol=0.00e+00
  18 : 1.700 µm, [0.00e+00; 0.00e+00], 1:0.054 s, tol=0.00e+00
  19 : 1.800 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
  20 : 1.900 µm, [0.00e+00; 0.00e+00], 1:0.046 s, tol=0.00e+00
  21 : 2.000 µm, [0.00e+00; 0.00e+00], 1:0.023 s, tol=0.00e+00
  22 : 2.100 µm, [0.00e+00; 0.00e+00], 1:0.042 s, tol=0.00e+00
  23 : 2.200 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  24 : 2.300 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  25 : 2.400 µm, [0.00e+00; 0.00e+00], 1:0.058 s, tol=0.00e+00
  26 : 2.500 µm, [0.00e+00; 0.00e+00], 1:0.058 s, tol=0.00e+00
  27 : 2.600 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  28 : 2.700 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
  29 : 2.800 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
  30 : 2.900 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  31 : 3.000 µm, [0.00e+00; 0.00e+00], 1:0.057 s, tol=0.00e+00
  32 : 3.100 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  33 : 3.200 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  34 : 3.300 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  35 : 3.400 µm, [0.00e+00; 0.00e+00], 1:0.057 s, tol=0.00e+00
  36 : 3.500 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
  37 : 3.600 µm, [0.00e+00; 0.00e+00], 1:0.059 s, tol=0.00e+00
  38 : 3.700 µm, [0.00e+00; 0.00e+00], 1:0.057 s, tol=0.00e+00
  39 : 3.800 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
  40 : 3.900 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  41 : 4.000 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  42 : 4.100 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  43 : 4.110 µm, [0.00e+00; 0.00e+00], 1:0.058 s, tol=0.00e+00
  44 : 4.120 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
  45 : 4.130 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
  46 : 4.140 µm, [0.00e+00; 0.00e+00], 1:0.054 s, tol=0.00e+00
  47 : 4.150 µm, [0.00e+00; 0.00e+00], 1:0.059 s, tol=0.00e+00
  48 : 4.160 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
  49 : 4.170 µm, [0.00e+00; 0.00e+00], 1:0.052 s, tol=0.00e+00
  50 : 4.180 µm, [0.00e+00; 0.00e+00], 1:0.058 s, tol=0.00e+00
  51 : 4.190 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  52 : 4.200 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  53 : 4.210 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  54 : 4.220 µm, [0.00e+00; 0.00e+00], 1:0.057 s, tol=0.00e+00
  55 : 4.230 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  56 : 4.240 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
  57 : 4.250 µm, [0.00e+00; 0.00e+00], 1:0.024 s, tol=0.00e+00
  58 : 4.260 µm, [0.00e+00; 0.00e+00], 1:0.023 s, tol=0.00e+00
  59 : 4.270 µm, [0.00e+00; 0.00e+00], 1:0.024 s, tol=0.00e+00
  60 : 4.280 µm, [0.00e+00; 0.00e+00], 1:0.051 s, tol=0.00e+00
  61 : 4.290 µm, [0.00e+00; 0.00e+00], 1:0.057 s, tol=0.00e+00
  62 : 4.300 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
  63 : 4.310 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
  64 : 4.320 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  65 : 4.330 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  66 : 4.340 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  67 : 4.350 µm, [0.00e+00; 0.00e+00], 1:0.057 s, tol=0.00e+00
  68 : 4.360 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
  69 : 4.370 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
  70 : 4.380 µm, [0.00e+00; 0.00e+00], 1:0.058 s, tol=0.00e+00
  71 : 4.390 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  72 : 4.400 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
  73 : 4.410 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
  74 : 4.420 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
  75 : 4.430 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  76 : 4.440 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  77 : 4.450 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
  78 : 4.460 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
  79 : 4.470 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  80 : 4.480 µm, [0.00e+00; 0.00e+00], 1:0.058 s, tol=0.00e+00
  81 : 4.490 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
  82 : 4.500 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
  83 : 4.510 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  84 : 4.520 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  85 : 4.530 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  86 : 4.540 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  87 : 4.550 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  88 : 4.560 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  89 : 4.570 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  90 : 4.580 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  91 : 4.590 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  92 : 4.600 µm, [0.00e+00; 0.00e+00], 1:0.058 s, tol=0.00e+00
  93 : 4.610 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
  94 : 4.620 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=0.00e+00
  95 : 4.630 µm, [0.00e+00; 0.00e+00], 1:0.023 s, tol=0.00e+00
  96 : 4.640 µm, [0.00e+00; 0.00e+00], 1:0.050 s, tol=0.00e+00
  97 : 4.650 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
  98 : 4.660 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
  99 : 4.670 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 100 : 4.680 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 101 : 4.690 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 102 : 4.700 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 103 : 4.710 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 104 : 4.720 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 105 : 4.730 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 106 : 4.740 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 107 : 4.750 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 108 : 4.760 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 109 : 4.770 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 110 : 4.780 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 111 : 4.790 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 112 : 4.800 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 113 : 4.810 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 114 : 4.820 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 115 : 4.830 µm, [0.00e+00; 0.00e+00], 1:0.054 s, tol=0.00e+00
 116 : 4.840 µm, [0.00e+00; 0.00e+00], 1:0.057 s, tol=0.00e+00
 117 : 4.850 µm, [0.00e+00; 0.00e+00], 1:0.059 s, tol=0.00e+00
 118 : 4.860 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 119 : 4.870 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 120 : 4.880 µm, [0.00e+00; 0.00e+00], 1:0.071 s, tol=0.00e+00
 121 : 4.890 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 122 : 4.900 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 123 : 4.910 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 124 : 4.920 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
 125 : 4.930 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 126 : 4.940 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 127 : 4.950 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 128 : 4.960 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
 129 : 4.970 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 130 : 4.980 µm, [0.00e+00; 0.00e+00], 1:0.052 s, tol=0.00e+00
 131 : 4.990 µm, [0.00e+00; 0.00e+00], 1:0.068 s, tol=0.00e+00
 132 : 5.000 µm, [0.00e+00; 0.00e+00], 1:0.058 s, tol=0.00e+00
 133 : 5.010 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 134 : 5.020 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 135 : 5.030 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 136 : 5.040 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 137 : 5.050 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 138 : 5.060 µm, [0.00e+00; 0.00e+00], 1:0.054 s, tol=0.00e+00
 139 : 5.070 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 140 : 5.080 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 141 : 5.090 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 142 : 5.100 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
 143 : 5.110 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 144 : 5.120 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 145 : 5.130 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 146 : 5.140 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 147 : 5.150 µm, [0.00e+00; 0.00e+00], 1:0.057 s, tol=0.00e+00
 148 : 5.160 µm, [0.00e+00; 0.00e+00], 1:0.043 s, tol=0.00e+00
 149 : 5.170 µm, [0.00e+00; 0.00e+00], 1:0.027 s, tol=0.00e+00
 150 : 5.180 µm, [0.00e+00; 0.00e+00], 1:0.025 s, tol=0.00e+00
 151 : 5.190 µm, [0.00e+00; 0.00e+00], 1:0.023 s, tol=0.00e+00
 152 : 5.200 µm, [0.00e+00; 0.00e+00], 1:0.030 s, tol=0.00e+00
 153 : 5.210 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 154 : 5.220 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 155 : 5.230 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
 156 : 5.240 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 157 : 5.250 µm, [0.00e+00; 0.00e+00], 1:0.033 s, tol=0.00e+00
 158 : 5.260 µm, [0.00e+00; 0.00e+00], 1:0.034 s, tol=0.00e+00
 159 : 5.270 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 160 : 5.280 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 161 : 5.290 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 162 : 5.300 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 163 : 5.310 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 164 : 5.320 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 165 : 5.330 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 166 : 5.340 µm, [0.00e+00; 0.00e+00], 1:0.057 s, tol=0.00e+00
 167 : 5.350 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 168 : 5.360 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
 169 : 5.370 µm, [0.00e+00; 0.00e+00], 1:0.028 s, tol=0.00e+00
 170 : 5.380 µm, [0.00e+00; 0.00e+00], 1:0.035 s, tol=0.00e+00
 171 : 5.390 µm, [0.00e+00; 0.00e+00], 1:0.052 s, tol=0.00e+00
 172 : 5.400 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 173 : 5.410 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 174 : 5.420 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 175 : 5.430 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 176 : 5.440 µm, [0.00e+00; 0.00e+00], 1:0.056 s, tol=0.00e+00
 177 : 5.450 µm, [0.00e+00; 0.00e+00], 1:0.053 s, tol=0.00e+00
 178 : 5.460 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 179 : 5.470 µm, [0.00e+00; 0.00e+00], 1:0.058 s, tol=0.00e+00
 180 : 5.480 µm, [0.00e+00; 0.00e+00], 1:0.055 s, tol=0.00e+00
 181 : 5.490 µm, [3.22e-11; 2.30e-04], 1:0.053 s, tol=1.00e+00
 182 : 5.500 µm, [1.10e-10; 8.12e-04], 1:0.055 s, tol=9.20e-01
 183 : 5.510 µm, [1.95e-10; 1.49e-03], 1:0.059 s, tol=7.05e-01
 184 : 5.520 µm, [2.78e-10; 2.20e-03], 1:0.059 s, tol=5.41e-01
 185 : 5.530 µm, [3.56e-10; 2.92e-03], 1:0.060 s, tol=4.32e-01
 186 : 5.540 µm, [4.29e-10; 3.63e-03], 1:0.057 s, tol=3.57e-01
 187 : 5.550 µm, [4.97e-10; 4.35e-03], 1:0.056 s, tol=3.04e-01
 188 : 5.560 µm, [5.60e-10; 5.07e-03], 1:0.053 s, tol=2.65e-01
 189 : 5.570 µm, [6.18e-10; 5.79e-03], 1:0.057 s, tol=2.34e-01
 190 : 5.580 µm, [6.71e-10; 6.50e-03], 1:0.055 s, tol=2.10e-01
 191 : 5.590 µm, [7.20e-10; 7.22e-03], 1:0.060 s, tol=1.90e-01
 192 : 5.600 µm, [7.65e-10; 7.93e-03], 1:0.058 s, tol=1.74e-01
 193 : 5.610 µm, [8.06e-10; 8.65e-03], 1:0.060 s, tol=1.60e-01
 194 : 5.620 µm, [8.43e-10; 9.37e-03], 1:0.060 s, tol=1.49e-01
 195 : 5.630 µm, [8.77e-10; 1.01e-02], 1:0.056 s, tol=1.39e-01
 196 : 5.640 µm, [9.08e-10; 1.08e-02], 1:0.059 s, tol=1.30e-01
 197 : 5.650 µm, [9.35e-10; 1.15e-02], 1:0.053 s, tol=1.22e-01
 198 : 5.660 µm, [9.60e-10; 1.22e-02], 1:0.056 s, tol=1.15e-01
 199 : 5.670 µm, [9.82e-10; 1.29e-02], 1:0.061 s, tol=1.09e-01
 200 : 5.680 µm, [1.00e-09; 1.36e-02], 1:0.058 s, tol=1.04e-01
 201 : 5.690 µm, [1.02e-09; 1.44e-02], 1:0.063 s, tol=9.86e-02
 202 : 5.700 µm, [1.03e-09; 1.51e-02], 1:0.058 s, tol=9.41e-02
 203 : 5.710 µm, [1.04e-09; 1.58e-02], 1:0.061 s, tol=9.00e-02
 204 : 5.720 µm, [1.05e-09; 1.65e-02], 1:0.037 s, tol=8.62e-02
 205 : 5.730 µm, [1.06e-09; 1.72e-02], 1:0.023 s, tol=8.28e-02
 206 : 5.740 µm, [1.07e-09; 1.79e-02], 1:0.045 s, tol=7.96e-02
 207 : 5.750 µm, [1.07e-09; 1.86e-02], 1:0.060 s, tol=7.67e-02
 208 : 5.760 µm, [1.08e-09; 1.93e-02], 1:0.058 s, tol=7.40e-02
 209 : 5.770 µm, [1.08e-09; 2.00e-02], 1:0.057 s, tol=7.14e-02
 210 : 5.780 µm, [1.08e-09; 2.08e-02], 1:0.057 s, tol=6.91e-02
 211 : 5.790 µm, [1.08e-09; 2.15e-02], 1:0.061 s, tol=6.68e-02
 212 : 5.800 µm, [1.07e-09; 2.22e-02], 1:0.060 s, tol=6.48e-02
 213 : 5.810 µm, [1.07e-09; 2.29e-02], 1:0.057 s, tol=6.28e-02
 214 : 5.820 µm, [1.06e-09; 2.36e-02], 1:0.059 s, tol=6.10e-02
 215 : 5.830 µm, [1.06e-09; 2.43e-02], 1:0.058 s, tol=5.93e-02
 216 : 5.840 µm, [1.05e-09; 2.50e-02], 1:0.057 s, tol=5.77e-02
 217 : 5.850 µm, [1.04e-09; 2.57e-02], 1:0.059 s, tol=5.61e-02
 218 : 5.860 µm, [1.04e-09; 2.64e-02], 1:0.058 s, tol=5.47e-02
 219 : 5.870 µm, [1.03e-09; 2.71e-02], 1:0.058 s, tol=5.33e-02
 220 : 5.880 µm, [1.02e-09; 2.78e-02], 1:0.061 s, tol=5.20e-02
 221 : 5.890 µm, [1.01e-09; 2.85e-02], 1:0.056 s, tol=5.07e-02
 222 : 5.900 µm, [9.96e-10; 2.92e-02], 1:0.061 s, tol=4.95e-02
 223 : 5.910 µm, [9.84e-10; 2.99e-02], 1:0.061 s, tol=4.84e-02
 224 : 5.920 µm, [9.72e-10; 3.06e-02], 1:0.053 s, tol=4.73e-02
 225 : 5.930 µm, [1.00e-09; 3.14e-02], 1:0.056 s, tol=5.78e-02
 226 : 5.940 µm, [1.05e-09; 3.23e-02], 1:0.059 s, tol=6.26e-02
 227 : 5.950 µm, [1.10e-09; 3.32e-02], 1:0.060 s, tol=6.34e-02
 228 : 5.960 µm, [1.14e-09; 3.41e-02], 1:0.055 s, tol=6.27e-02
 229 : 5.970 µm, [1.18e-09; 3.49e-02], 1:0.053 s, tol=6.15e-02
 230 : 5.980 µm, [1.22e-09; 3.58e-02], 1:0.055 s, tol=6.01e-02
 231 : 5.990 µm, [1.25e-09; 3.67e-02], 1:0.055 s, tol=5.88e-02
 232 : 6.000 µm, [1.28e-09; 3.76e-02], 1:0.060 s, tol=5.75e-02
 233 : 6.010 µm, [1.30e-09; 3.85e-02], 1:0.053 s, tol=5.62e-02
 234 : 6.020 µm, [1.33e-09; 3.94e-02], 1:0.059 s, tol=5.49e-02
 235 : 6.030 µm, [1.34e-09; 4.03e-02], 1:0.055 s, tol=5.37e-02
 236 : 6.040 µm, [1.36e-09; 4.11e-02], 1:0.055 s, tol=5.25e-02
 237 : 6.050 µm, [1.37e-09; 4.20e-02], 1:0.059 s, tol=5.14e-02
 238 : 6.060 µm, [1.38e-09; 4.29e-02], 1:0.059 s, tol=5.03e-02
 239 : 6.070 µm, [1.39e-09; 4.38e-02], 1:0.040 s, tol=4.92e-02
 240 : 6.080 µm, [1.40e-09; 4.46e-02], 1:0.035 s, tol=4.82e-02
 241 : 6.090 µm, [1.40e-09; 4.55e-02], 1:0.060 s, tol=4.72e-02
 242 : 6.100 µm, [1.40e-09; 4.64e-02], 1:0.059 s, tol=4.63e-02
 243 : 6.110 µm, [1.41e-09; 4.73e-02], 1:0.059 s, tol=4.54e-02
 244 : 6.120 µm, [1.40e-09; 4.81e-02], 1:0.064 s, tol=4.45e-02
 245 : 6.130 µm, [1.40e-09; 4.90e-02], 1:0.061 s, tol=4.36e-02
 246 : 6.140 µm, [1.40e-09; 4.98e-02], 1:0.061 s, tol=4.28e-02
 247 : 6.150 µm, [1.39e-09; 5.07e-02], 1:0.059 s, tol=4.20e-02
 248 : 6.160 µm, [1.38e-09; 5.16e-02], 1:0.058 s, tol=4.12e-02
 249 : 6.170 µm, [1.37e-09; 5.24e-02], 1:0.065 s, tol=4.05e-02
 250 : 6.180 µm, [1.36e-09; 5.33e-02], 1:0.057 s, tol=3.97e-02
 251 : 6.190 µm, [1.35e-09; 5.41e-02], 1:0.057 s, tol=3.90e-02
 252 : 6.200 µm, [1.34e-09; 5.50e-02], 1:0.057 s, tol=3.84e-02
 253 : 6.210 µm, [1.33e-09; 5.58e-02], 1:0.059 s, tol=3.77e-02
 254 : 6.220 µm, [1.31e-09; 5.67e-02], 1:0.057 s, tol=3.71e-02
 255 : 6.230 µm, [1.30e-09; 5.75e-02], 1:0.056 s, tol=3.64e-02
 256 : 6.240 µm, [1.28e-09; 5.84e-02], 1:0.061 s, tol=3.58e-02
 257 : 6.250 µm, [1.27e-09; 5.94e-02], 1:0.061 s, tol=3.53e-02
 258 : 6.260 µm, [1.25e-09; 6.04e-02], 1:0.057 s, tol=3.47e-02
 259 : 6.270 µm, [1.23e-09; 6.15e-02], 1:0.057 s, tol=3.42e-02
 260 : 6.280 µm, [1.22e-09; 6.25e-02], 1:0.057 s, tol=3.36e-02
 261 : 6.290 µm, [1.20e-09; 6.36e-02], 1:0.061 s, tol=3.31e-02
 262 : 6.300 µm, [1.18e-09; 6.46e-02], 1:0.063 s, tol=3.26e-02
 263 : 6.310 µm, [1.16e-09; 6.56e-02], 1:0.055 s, tol=3.21e-02
 264 : 6.320 µm, [1.14e-09; 6.66e-02], 1:0.055 s, tol=3.17e-02
 265 : 6.330 µm, [1.12e-09; 6.77e-02], 1:0.059 s, tol=3.12e-02
 266 : 6.340 µm, [1.10e-09; 6.87e-02], 1:0.057 s, tol=3.08e-02
 267 : 6.350 µm, [1.08e-09; 6.97e-02], 1:0.057 s, tol=3.03e-02
 268 : 6.360 µm, [1.06e-09; 7.07e-02], 1:0.057 s, tol=2.99e-02
 269 : 6.370 µm, [1.05e-09; 7.17e-02], 1:0.060 s, tol=2.95e-02
 270 : 6.380 µm, [1.03e-09; 7.28e-02], 1:0.054 s, tol=2.91e-02
 271 : 6.390 µm, [1.01e-09; 7.38e-02], 1:0.061 s, tol=2.87e-02
 272 : 6.400 µm, [9.86e-10; 7.48e-02], 1:0.057 s, tol=2.83e-02
 273 : 6.410 µm, [9.66e-10; 7.58e-02], 1:0.053 s, tol=2.80e-02
 274 : 6.420 µm, [9.46e-10; 7.68e-02], 1:0.055 s, tol=2.76e-02
 275 : 6.430 µm, [9.26e-10; 7.78e-02], 1:0.058 s, tol=2.73e-02
 276 : 6.440 µm, [9.07e-10; 7.88e-02], 1:0.060 s, tol=2.69e-02
 277 : 6.450 µm, [8.87e-10; 7.98e-02], 1:0.055 s, tol=2.66e-02
 278 : 6.460 µm, [8.68e-10; 8.08e-02], 1:0.062 s, tol=2.63e-02
 279 : 6.470 µm, [8.48e-10; 8.18e-02], 1:0.057 s, tol=2.59e-02
 280 : 6.480 µm, [8.29e-10; 8.28e-02], 1:0.061 s, tol=2.56e-02
 281 : 6.490 µm, [8.10e-10; 8.38e-02], 1:0.057 s, tol=2.53e-02
 282 : 6.500 µm, [7.91e-10; 8.48e-02], 1:0.031 s, tol=2.50e-02
 283 : 6.510 µm, [7.73e-10; 8.58e-02], 1:0.022 s, tol=2.47e-02
 284 : 6.520 µm, [7.54e-10; 8.67e-02], 1:0.021 s, tol=2.44e-02
 285 : 6.530 µm, [7.36e-10; 8.77e-02], 1:0.021 s, tol=2.42e-02
 286 : 6.540 µm, [7.18e-10; 8.87e-02], 1:0.022 s, tol=2.39e-02
 287 : 6.550 µm, [7.00e-10; 8.97e-02], 1:0.024 s, tol=2.36e-02
 288 : 6.560 µm, [6.82e-10; 9.07e-02], 1:0.024 s, tol=2.34e-02
 289 : 6.570 µm, [6.65e-10; 9.16e-02], 1:0.052 s, tol=2.31e-02
 290 : 6.580 µm, [6.48e-10; 9.26e-02], 1:0.058 s, tol=2.29e-02
 291 : 6.590 µm, [6.31e-10; 9.36e-02], 1:0.061 s, tol=2.26e-02
 292 : 6.600 µm, [6.14e-10; 9.46e-02], 1:0.061 s, tol=2.24e-02
 293 : 6.610 µm, [5.98e-10; 9.55e-02], 1:0.061 s, tol=2.22e-02
 294 : 6.620 µm, [5.82e-10; 9.65e-02], 1:0.053 s, tol=2.19e-02
 295 : 6.630 µm, [5.66e-10; 9.74e-02], 1:0.061 s, tol=2.17e-02
 296 : 6.640 µm, [5.50e-10; 9.84e-02], 1:0.057 s, tol=2.15e-02
 297 : 6.650 µm, [5.36e-10; 9.95e-02], 1:0.056 s, tol=2.40e-02
 298 : 6.660 µm, [5.22e-10; 1.01e-01], 1:0.061 s, tol=2.80e-02
 299 : 6.670 µm, [5.09e-10; 1.02e-01], 1:0.053 s, tol=3.00e-02
 300 : 6.680 µm, [4.96e-10; 1.04e-01], 1:0.055 s, tol=3.07e-02
 301 : 6.690 µm, [4.83e-10; 1.05e-01], 1:0.059 s, tol=3.09e-02
 302 : 6.700 µm, [4.72e-10; 1.07e-01], 1:0.059 s, tol=3.30e-02
 303 : 6.710 µm, [4.64e-10; 1.09e-01], 1:0.055 s, tol=4.01e-02
 304 : 6.720 µm, [4.57e-10; 1.12e-01], 1:0.057 s, tol=4.52e-02
 305 : 6.730 µm, [4.52e-10; 1.15e-01], 1:0.060 s, tol=4.81e-02
 306 : 6.740 µm, [4.47e-10; 1.18e-01], 1:0.057 s, tol=4.96e-02
 307 : 6.750 µm, [4.42e-10; 1.21e-01], 1:0.053 s, tol=5.02e-02
 308 : 6.760 µm, [4.37e-10; 1.24e-01], 1:0.064 s, tol=5.01e-02
 309 : 6.770 µm, [4.32e-10; 1.27e-01], 1:0.057 s, tol=4.97e-02
 310 : 6.780 µm, [4.26e-10; 1.30e-01], 1:0.060 s, tol=4.91e-02
 311 : 6.790 µm, [4.20e-10; 1.33e-01], 1:0.066 s, tol=4.84e-02
 312 : 6.800 µm, [4.14e-10; 1.37e-01], 1:0.025 s, tol=4.76e-02
 313 : 6.810 µm, [4.08e-10; 1.40e-01], 1:0.039 s, tol=4.68e-02
 314 : 6.820 µm, [4.01e-10; 1.43e-01], 1:0.060 s, tol=4.60e-02
 315 : 6.830 µm, [3.94e-10; 1.46e-01], 1:0.060 s, tol=4.52e-02
 316 : 6.840 µm, [3.87e-10; 1.49e-01], 1:0.057 s, tol=4.44e-02
 317 : 6.850 µm, [3.80e-10; 1.52e-01], 1:0.061 s, tol=4.36e-02
 318 : 6.860 µm, [3.73e-10; 1.56e-01], 1:0.058 s, tol=4.99e-02
 319 : 6.870 µm, [3.67e-10; 1.60e-01], 1:0.059 s, tol=5.47e-02
 320 : 6.880 µm, [3.60e-10; 1.65e-01], 1:0.057 s, tol=5.79e-02
 321 : 6.890 µm, [3.55e-10; 1.70e-01], 1:0.056 s, tol=6.00e-02
 322 : 6.900 µm, [3.49e-10; 1.75e-01], 1:0.056 s, tol=6.12e-02
 323 : 6.910 µm, [3.43e-10; 1.80e-01], 1:0.053 s, tol=6.17e-02
 324 : 6.920 µm, [3.40e-10; 1.86e-01], 1:0.057 s, tol=6.60e-02
 325 : 6.930 µm, [3.40e-10; 1.92e-01], 1:0.061 s, tol=7.21e-02
 326 : 6.940 µm, [3.42e-10; 2.00e-01], 1:0.058 s, tol=7.72e-02
 327 : 6.950 µm, [3.45e-10; 2.08e-01], 1:0.060 s, tol=8.11e-02
 328 : 6.960 µm, [3.49e-10; 2.16e-01], 1:0.061 s, tol=8.41e-02
 329 : 6.970 µm, [3.53e-10; 2.25e-01], 1:0.057 s, tol=8.62e-02
 330 : 6.980 µm, [3.57e-10; 2.35e-01], 1:0.056 s, tol=8.76e-02
 331 : 6.990 µm, [3.61e-10; 2.45e-01], 1:0.059 s, tol=8.85e-02
 332 : 7.000 µm, [3.65e-10; 2.56e-01], 1:0.061 s, tol=8.89e-02
 333 : 7.010 µm, [3.69e-10; 2.66e-01], 1:0.061 s, tol=8.89e-02
 334 : 7.020 µm, [3.73e-10; 2.77e-01], 1:0.064 s, tol=8.85e-02
 335 : 7.030 µm, [3.76e-10; 2.89e-01], 1:0.064 s, tol=8.80e-02
 336 : 7.040 µm, [3.78e-10; 3.00e-01], 1:0.057 s, tol=8.73e-02
 337 : 7.050 µm, [3.81e-10; 3.13e-01], 1:0.057 s, tol=9.31e-02
 338 : 7.060 µm, [3.86e-10; 3.28e-01], 1:0.060 s, tol=1.04e-01
 339 : 7.070 µm, [3.92e-10; 3.46e-01], 1:0.061 s, tol=1.17e-01
 340 : 7.080 µm, [4.01e-10; 3.66e-01], 1:0.061 s, tol=1.30e-01
 341 : 7.090 µm, [4.11e-10; 3.90e-01], 1:0.057 s, tol=1.42e-01
 342 : 7.100 µm, [4.24e-10; 4.16e-01], 1:0.061 s, tol=1.54e-01
 343 : 7.110 µm, [4.39e-10; 4.46e-01], 1:0.063 s, tol=1.66e-01
 344 : 7.120 µm, [4.56e-10; 4.80e-01], 1:0.060 s, tol=1.76e-01
 345 : 7.130 µm, [4.75e-10; 5.17e-01], 1:0.064 s, tol=1.86e-01
 346 : 7.140 µm, [4.95e-10; 5.58e-01], 1:0.064 s, tol=1.95e-01
 347 : 7.150 µm, [5.16e-10; 6.02e-01], 1:0.059 s, tol=2.09e-01
 348 : 7.160 µm, [5.37e-10; 6.51e-01], 1:0.051 s, tol=2.26e-01
 349 : 7.170 µm, [5.55e-10; 7.03e-01], 1:0.029 s, tol=2.44e-01
 350 : 7.180 µm, [5.67e-10; 7.57e-01], 1:0.050 s, tol=2.73e-01
 351 : 7.190 µm, [5.65e-10; 8.13e-01], 1:0.057 s, tol=3.30e-01
 352 : 7.200 µm, [5.61e-10; 8.73e-01], 1:0.057 s, tol=4.26e-01
 353 : 7.210 µm, [5.62e-10; 9.33e-01], 1:0.061 s, tol=5.09e-01
 354 : 7.220 µm, [5.68e-10; 9.75e-01], 1:0.062 s, tol=4.89e-01
 355 : 7.230 µm, [5.73e-10; 9.90e-01], 1:0.057 s, tol=3.64e-01
 356 : 7.240 µm, [5.81e-10; 9.90e-01], 1:0.062 s, tol=2.73e-01
 357 : 7.250 µm, [5.90e-10; 9.95e-01], 1:0.061 s, tol=2.23e-01
 358 : 7.260 µm, [5.97e-10; 9.99e-01], 1:0.061 s, tol=1.88e-01
 359 : 7.270 µm, [6.03e-10; 1.00e+00], 1:0.061 s, tol=1.78e-01
 360 : 7.280 µm, [6.10e-10; 1.00e+00], 1:0.062 s, tol=1.62e-01
 361 : 7.290 µm, [6.14e-10; 1.00e+00], 1:0.061 s, tol=1.52e-01
 362 : 7.300 µm, [6.17e-10; 1.00e+00], 1:0.061 s, tol=1.42e-01
 363 : 7.310 µm, [6.23e-10; 1.00e+00], 1:0.061 s, tol=1.18e-01
 364 : 7.320 µm, [6.25e-10; 1.00e+00], 1:0.062 s, tol=1.11e-01
 365 : 7.330 µm, [6.28e-10; 1.00e+00], 1:0.057 s, tol=1.11e-01
 366 : 7.340 µm, [6.34e-10; 1.00e+00], 1:0.057 s, tol=9.32e-02
 367 : 7.350 µm, [6.36e-10; 1.00e+00], 1:0.060 s, tol=8.10e-02
 368 : 7.360 µm, [6.37e-10; 1.00e+00], 1:0.061 s, tol=7.99e-02
 369 : 7.370 µm, [6.42e-10; 1.00e+00], 1:0.057 s, tol=8.51e-02
 370 : 7.380 µm, [6.42e-10; 1.01e+00], 1:0.057 s, tol=9.51e-02
 371 : 7.390 µm, [6.43e-10; 1.01e+00], 1:0.057 s, tol=8.91e-02
 372 : 7.400 µm, [6.46e-10; 1.01e+00], 1:0.057 s, tol=7.11e-02
 373 : 7.410 µm, [6.42e-10; 1.01e+00], 1:0.059 s, tol=6.24e-02
 374 : 7.420 µm, [6.37e-10; 1.01e+00], 1:0.057 s, tol=6.64e-02
 375 : 7.430 µm, [6.41e-10; 1.01e+00], 1:0.061 s, tol=7.39e-02
 376 : 7.440 µm, [6.36e-10; 1.01e+00], 1:0.061 s, tol=7.67e-02
 377 : 7.450 µm, [6.32e-10; 1.01e+00], 1:0.032 s, tol=7.30e-02
 378 : 7.460 µm, [6.32e-10; 1.01e+00], 1:0.025 s, tol=6.50e-02
 379 : 7.470 µm, [6.13e-10; 1.01e+00], 1:0.045 s, tol=6.81e-02
 380 : 7.480 µm, [6.04e-10; 1.01e+00], 1:0.060 s, tol=6.60e-02
 381 : 7.490 µm, [6.01e-10; 1.00e+00], 1:0.057 s, tol=5.40e-02
 382 : 7.500 µm, [5.77e-10; 1.00e+00], 1:0.057 s, tol=5.04e-02
 383 : 7.510 µm, [5.76e-10; 1.00e+00], 1:0.070 s, tol=5.08e-02
 384 : 7.520 µm, [5.59e-10; 1.00e+00], 1:0.063 s, tol=5.55e-02
 385 : 7.530 µm, [5.38e-10; 1.01e+00], 1:0.066 s, tol=5.31e-02
 386 : 7.540 µm, [5.25e-10; 1.01e+00], 1:0.063 s, tol=4.30e-02
 387 : 7.550 µm, [4.73e-10; 1.00e+00], 1:0.057 s, tol=3.96e-02
 388 : 7.560 µm, [4.40e-10; 1.00e+00], 1:0.057 s, tol=3.94e-02
 389 : 7.570 µm, [4.04e-10; 1.00e+00], 1:0.061 s, tol=4.37e-02
 390 : 7.580 µm, [3.44e-10; 1.01e+00], 1:0.061 s, tol=4.78e-02
 391 : 7.590 µm, [3.14e-10; 1.01e+00], 1:0.062 s, tol=4.37e-02
 392 : 7.600 µm, [2.39e-10; 1.01e+00], 1:0.061 s, tol=4.01e-02
 393 : 7.610 µm, [1.57e-10; 1.01e+00], 1:0.062 s, tol=4.05e-02
 394 : 7.620 µm, [1.23e-10; 1.01e+00], 1:0.057 s, tol=4.04e-02
 395 : 7.630 µm, [2.97e-11; 1.01e+00], 1:0.057 s, tol=4.07e-02
 396 : 7.640 µm, [-1.22e-12; 1.01e+00], 1:0.061 s, tol=3.94e-02
 397 : 7.650 µm, [-4.92e-11; 1.01e+00], 1:0.061 s, tol=3.89e-02
 398 : 7.660 µm, [-1.63e-10; 1.01e+00], 1:0.061 s, tol=3.78e-02
 399 : 7.670 µm, [-2.93e-10; 1.01e+00], 1:0.057 s, tol=3.26e-02
 400 : 7.680 µm, [-3.82e-10; 1.00e+00], 1:0.063 s, tol=2.67e-02
 401 : 7.690 µm, [-5.54e-10; 1.00e+00], 1:0.057 s, tol=2.73e-02
 402 : 7.700 µm, [-7.68e-10; 1.00e+00], 1:0.061 s, tol=3.29e-02
 403 : 7.710 µm, [-9.49e-10; 1.00e+00], 1:0.062 s, tol=3.54e-02
 404 : 7.720 µm, [-1.13e-09; 1.00e+00], 1:0.061 s, tol=3.55e-02
 405 : 7.730 µm, [-1.37e-09; 1.00e+00], 1:0.061 s, tol=3.75e-02
 406 : 7.740 µm, [-1.59e-09; 1.00e+00], 1:0.057 s, tol=3.31e-02
 407 : 7.750 µm, [-1.81e-09; 1.00e+00], 1:0.061 s, tol=3.03e-02
 408 : 7.760 µm, [-1.97e-09; 1.01e+00], 1:0.057 s, tol=3.18e-02
 409 : 7.770 µm, [-2.27e-09; 1.01e+00], 1:0.057 s, tol=3.65e-02
 410 : 7.780 µm, [-2.61e-09; 1.01e+00], 1:0.061 s, tol=3.77e-02
 411 : 7.790 µm, [-2.88e-09; 1.01e+00], 1:0.063 s, tol=3.59e-02
 412 : 7.800 µm, [-3.30e-09; 1.01e+00], 1:0.061 s, tol=3.50e-02
 413 : 7.810 µm, [-3.79e-09; 1.01e+00], 1:0.061 s, tol=2.97e-02
 414 : 7.820 µm, [-4.37e-09; 1.01e+00], 1:0.057 s, tol=2.78e-02
 415 : 7.830 µm, [-4.88e-09; 1.01e+00], 1:0.061 s, tol=3.12e-02
 416 : 7.840 µm, [-5.54e-09; 1.01e+00], 1:0.057 s, tol=3.46e-02
 417 : 7.850 µm, [-6.16e-09; 1.01e+00], 1:0.061 s, tol=3.47e-02
 418 : 7.860 µm, [-6.64e-09; 1.01e+00], 1:0.061 s, tol=3.61e-02
 419 : 7.870 µm, [-7.02e-09; 1.01e+00], 1:0.057 s, tol=3.92e-02
 420 : 7.880 µm, [-7.38e-09; 1.00e+00], 1:0.057 s, tol=4.47e-02
 421 : 7.890 µm, [-7.83e-09; 1.00e+00], 1:0.053 s, tol=5.55e-02
 422 : 7.900 µm, [-8.12e-09; 1.00e+00], 1:0.057 s, tol=9.23e-02
 423 : 7.910 µm, [-8.56e-09; 1.00e+00], 1:0.061 s, tol=2.06e-01
 424 : 7.920 µm, [-8.56e-09; 1.00e+00], 1:0.057 s, tol=5.10e-02
 425 : 7.930 µm, [-8.56e-09; 1.00e+00], 1:0.055 s, tol=8.99e-03
 426 : 7.940 µm, [-8.56e-09; 1.00e+00], 1:0.062 s, tol=1.90e-03
 427 : 7.950 µm, [-8.56e-09; 1.00e+00], 1:0.061 s, tol=1.26e-03
 428 : 7.960 µm, [-8.56e-09; 1.00e+00], 1:0.057 s, tol=1.30e-03
 429 : 7.970 µm, [-8.56e-09; 1.00e+00], 1:0.058 s, tol=1.38e-03
 430 : 7.980 µm, [-8.56e-09; 1.00e+00], 1:0.061 s, tol=1.44e-03
 431 : 7.990 µm, [-8.56e-09; 1.00e+00], 1:0.061 s, tol=1.51e-03
 432 : 8.000 µm, [-8.56e-09; 1.00e+00], 1:0.061 s, tol=1.55e-03
 433 : 8.010 µm, [-8.56e-09; 1.00e+00], 1:0.061 s, tol=1.56e-03
 434 : 8.020 µm, [-8.56e-09; 1.00e+00], 1:0.061 s, tol=1.56e-03
Saved:
/home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/v1.7.1/examples/PhaseField/results/Tension2D/Test/Elas_Isot_Bourdin_AT1_DP_optimMesh/force-displacement.pickle
Saved:
/home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/v1.7.1/examples/PhaseField/results/Tension2D/Test/Elas_Isot_Bourdin_AT1_DP_optimMesh/simulation.pickle
Saved:
/home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/v1.7.1/examples/PhaseField/results/Tension2D/Test/Elas_Isot_Bourdin_AT1_DP_optimMesh/summary.txt

Generate movie 01/21 (4.76 %) 4.46 s
Generate movie 02/21 (9.52 %) 3.58 s
Generate movie 03/21 (14.29 %) 3.46 s
Generate movie 04/21 (19.05 %) 3.28 s
Generate movie 05/21 (23.81 %) 3.00 s
Generate movie 06/21 (28.57 %) 2.84 s
Generate movie 07/21 (33.33 %) 2.65 s
Generate movie 08/21 (38.10 %) 2.45 s
Generate movie 09/21 (42.86 %) 2.28 s
Generate movie 10/21 (47.62 %) 2.13 s
Generate movie 11/21 (52.38 %) 1.90 s
Generate movie 12/21 (57.14 %) 1.71 s
Generate movie 13/21 (61.90 %) 1.53 s
Generate movie 14/21 (66.67 %) 1.34 s
Generate movie 15/21 (71.43 %) 1.15 s
Generate movie 16/21 (76.19 %) 952.26 ms
Generate movie 17/21 (80.95 %) 759.61 ms
Generate movie 18/21 (85.71 %) 560.64 ms
Generate movie 19/21 (90.48 %) 381.11 ms
Generate movie 20/21 (95.24 %) 191.45 ms
Generate movie 21/21 (100.00 %) 0.00 µs
Mesh : 48.226 ms
Boundary Conditions : 8.029 ms
Matrix : 10.934 s
Split : 48.551 ms
PostProcessing : 133.003 ms
Solver : 14.233 s
Resolution phase field : 23.907 s
Display : 838.299 ms
PyVista_Interface : 6.400 s


 12 import matplotlib.pyplot as plt
 13 import numpy as np
 14
 15 from EasyFEA import (
 16     Display,
 17     Folder,
 18     Models,
 19     Tic,
 20     ElemType,
 21     Mesh,
 22     Simulations,
 23     PyVista,
 24     Paraview,
 25 )
 26 from EasyFEA.Geoms import Point, Points, Domain, Line, Contour
 27
 28 import multiprocessing
 29
 30 # Display.Clear()
 31
 32 useParallel = False
 33 nProcs = 4  # number of processes in parallel
 34
 35 # ----------------------------------------------
 36 # Configurations
 37 # ----------------------------------------------
 38 dim = 2
 39
 40 # simu options
 41 doSimu = True
 42 meshTest = True
 43 openCrack = True
 44 optimMesh = True
 45
 46 # outputs
 47 folder = Folder.Results_Dir() + f"{dim}D"
 48 plotMesh = False
 49 plotResult = True
 50 showResult = True
 51 plotEnergy = False
 52 saveParaview = False
 53 makeMovie = True
 54
 55 # material
 56 materialType = "Elas_Isot"  #  "Elas_Isot", "ElasAnisot"
 57
 58 # phasefield
 59 maxIter = 1000
 60 tolConv = 1e-0  # 1e-1, 1e-2, 1e-3
 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 = ["Bourdin"]
 67
 68 regus = ["AT1"]  # "AT1", "AT2"
 69 # regus = ["AT1", "AT2"]
 70
 71 thetas = [-70, -80, -90]  # [-0, -10, -20, -30, -45, -60]
 72 theta = -0  # default value
 73
 74 # ----------------------------------------------
 75 # Mesh
 76 # ----------------------------------------------
 77 L = 1e-3
 78 # m
 79 l0 = 8.5e-6 if materialType == "ElasAnisot" else 1e-5
 80 thickness = 1 if dim == 2 else 0.1 / 1000
 81
 82
 83 def DoMesh(materialType: str = "Elas_Isot") -> Mesh:
 84     # meshSize
 85     clC = l0 * 2 if meshTest else l0 / 2
 86     if optimMesh:
 87         # a coarser mesh can be used outside the refined zone
 88         clD = clC * 4
 89         # refines the mesh in the area where the crack will propagate
 90         gap = L * 0.05
 91         if materialType == "Elas_Isot":
 92             refineDomain = Domain(
 93                 Point(L / 2 - gap, L / 2 - gap), Point(L, L / 2 + gap, thickness), clC
 94             )
 95         else:
 96             refineDomain = Domain(
 97                 Point(L / 2 - gap, L / 2 - gap), Point(L, L * 0.8, thickness), clC
 98             )
 99     else:
100         clD = clC
101         refineDomain = None
102
103     # geom
104     pt1 = Point()
105     pt2 = Point(L)
106     pt3 = Point(L, L)
107     pt4 = Point(0, L)
108     contour = Points([pt1, pt2, pt3, pt4], clD)
109
110     if dim == 2:
111         ptC1 = Point(0, L / 2, isOpen=openCrack)
112         ptC2 = Point(L / 2, L / 2)
113         cracks = [Line(ptC1, ptC2, clC, isOpen=openCrack)]
114     if dim == 3:
115         meshSize = clD if optimMesh else clC
116         ptC1 = Point(0, L / 2, 0, isOpen=openCrack)
117         ptC2 = Point(L / 2, L / 2, 0)
118         ptC3 = Point(L / 2, L / 2, thickness)
119         ptC4 = Point(0, L / 2, thickness, isOpen=openCrack)
120         l1 = Line(ptC1, ptC2, meshSize, openCrack)
121         l2 = Line(ptC2, ptC3, meshSize, False)
122         l3 = Line(ptC3, ptC4, meshSize, openCrack)
123         l4 = Line(ptC4, ptC1, meshSize, openCrack)
124         cracks = [Contour([l1, l2, l3, l4])]
125
126     if dim == 2:
127         mesh = contour.Mesh_2D([], ElemType.TRI3, cracks, [refineDomain])
128     elif dim == 3:
129         mesh = contour.Mesh_Extrude(
130             [], [0, 0, thickness], [3], ElemType.TETRA4, cracks, [refineDomain]
131         )
132
133     return mesh
134
135
136 # ----------------------------------------------
137 # Simu
138 # ----------------------------------------------
139
140
141 def DoSimu(split: str, regu: str):
142     # Builds the path to the folder based on the problem data
143     folder_save = Simulations.PhaseField.Folder(
144         folder,
145         materialType,
146         split,
147         regu,
148         "DP",
149         tolConv,
150         pfmSolver,
151         meshTest,
152         optimMesh,
153         not openCrack,
154         theta=theta,
155     )
156
157     Display.MyPrint(folder_save, "green")
158
159     if doSimu:
160         mesh = DoMesh(materialType)
161
162         # Nodes recovery
163         nodes_upper = mesh.Nodes_Conditions(lambda x, y, z: y == L)
164         nodes_lower = mesh.Nodes_Conditions(lambda x, y, z: y == 0)
165         nodes_right = mesh.Nodes_Conditions(
166             lambda x, y, z: (x == L) & (y > 0) & (y < L)
167         )
168         nodes_crack = mesh.Nodes_Conditions(lambda x, y, z: (y == L / 2) & (x <= L / 2))
169         if openCrack:
170             nodes_detect = mesh.nodes.copy()
171         else:
172             nodes_detect = np.array(list(set(mesh.nodes) - set(nodes_crack)))
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         if materialType == "Elas_Isot":
183             material = Models.Elastic.Isotropic(
184                 dim, E=210e9, v=0.3, planeStress=False, thickness=thickness
185             )
186             Gc = 2.7e3  # J/m2
187         elif materialType == "ElasAnisot":
188             if dim == 2:
189                 c11 = 65
190                 c22 = 260
191                 c33 = 30
192                 c12 = 20
193                 C_voigt = np.array([[c11, c12, 0], [c12, c22, 0], [0, 0, c33]]) * 1e9
194
195                 theta_rad = theta * np.pi / 180
196                 axis1 = np.array([np.cos(theta_rad), np.sin(theta_rad), 0])
197                 axis2 = np.array([-np.sin(theta_rad), np.cos(theta_rad), 0])
198
199                 material = Models.ElasAnisot(
200                     dim,
201                     C=C_voigt,
202                     useVoigtNotation=True,
203                     axis1=axis1,
204                     axis2=axis2,
205                     planeStress=False,
206                     thickness=thickness,
207                 )
208                 Gc = 1e3  # J/m2
209             else:
210                 raise Exception("Not implemented in 3D")
211
212         pfm = Models.PhaseField(material, split, regu, Gc=Gc, l0=l0, solver=pfmSolver)
213
214         # ----------------------------------------------
215         # Boundary conditions
216         # ----------------------------------------------
217         if materialType == "Elas_Isot":
218             # load < threshold
219             uinc0 = 1e-7 if meshTest else 1e-8
220             N0 = 40 if meshTest else 400
221             dep0 = uinc0 * N0
222
223             # load >= threshold
224             uinc1 = 1e-8 if meshTest else 1e-9
225             N1 = 400 if meshTest else 4000
226             dep1 = dep0 + uinc1 * N1
227
228             threshold = uinc0 * N0
229
230             config = f"""
231             while True: # simu until break
232
233             uinc0 = {uinc0:.1e};  N0 = {N0};  dep0 = uinc0*N0 = {dep0}
234             uinc1 = {uinc1:.1e};  N0 = {N1};  dep1 = dep0 + uinc1*N1 = {dep1}
235
236             threshold = uinc0*N0 = {threshold}
237
238             dep += uinc0 if dep < threshold else uinc1
239
240             if not openCrack:
241                 simu.add_dirichlet(nodes_crack, [1], ["d"], problemType="damage")
242             if dim == 2:
243                 simu.add_dirichlet(nodes_upper, [0,dep], ["x","y"])
244             elif dim == 3:
245                 simu.add_dirichlet(nodes_upper, [0,dep,0], ["x","y","z"])
246             simu.add_dirichlet(nodes_lower, [0],["y"])
247             """
248
249         else:
250             # load < threshold
251             uinc0 = 12e-8 if meshTest else 6e-8
252             # load >= threshold
253             uinc1 = 4e-8 if meshTest else 2e-8
254
255             threshold = 0.6
256
257             config = f"""
258             while True: # simu until break
259
260             uinc0 = {uinc0:.1e} (simu.damage.max() < {threshold})
261             uinc1 = {uinc1:.1e}
262
263             if not openCrack:
264                 simu.add_dirichlet(nodes_crack, [1], ["d"], problemType="damage")
265             if dim == 2:
266                 simu.add_dirichlet(nodes_upper, [0,dep], ["x","y"])
267             elif dim == 3:
268                 simu.add_dirichlet(nodes_upper, [0,dep,0], ["x","y","z"])
269             simu.add_dirichlet(nodes_lower, [0],["y"])
270             """
271
272         def Loading(dep):
273             """Boundary conditions"""
274
275             simu.Bc_Init()
276
277             if not openCrack:
278                 simu.add_dirichlet(nodes_crack, [1], ["d"], problemType="damage")
279
280             if dim == 2:
281                 simu.add_dirichlet(nodes_upper, [0, dep], ["x", "y"])
282             elif dim == 3:
283                 simu.add_dirichlet(nodes_upper, [0, dep, 0], ["x", "y", "z"])
284             simu.add_dirichlet(nodes_lower, [0], ["y"])
285
286         # ----------------------------------------------
287         # Simulation
288         # ----------------------------------------------
289         simu = Simulations.PhaseField(mesh, pfm, verbosity=False)
290         simu.Results_Set_Bc_Summary(config)
291
292         dofsY_upper = simu.Bc_dofs_nodes(nodes_upper, ["y"])
293
294         # INIT
295         nDetect = 0
296         displacement = []
297         force = []
298
299         dep = -uinc0
300         iter = -1
301         while True:  # simu until break
302             iter += 1
303             if materialType == "Elas_Isot":
304                 dep += uinc0 if dep < threshold else uinc1
305             else:
306                 if np.max(simu.damage[nodes_detect]) < threshold:
307                     dep += uinc0
308                 else:
309                     dep += uinc1
310
311             # apply new boundary conditions
312             Loading(dep)
313
314             # solve and save iter
315             u, _, Ku, converg = simu.Solve(tolConv, maxIter, convOption=1)
316             simu.Save_Iter()
317
318             # print iter solution
319             simu.Results_Set_Iteration_Summary(iter, dep * 1e6, "µm", 0, True)
320
321             # If the solver has not converged, stop the simulation.
322             if not converg:
323                 break
324
325             # resulting force on upper edge
326             f = np.sum(Ku[dofsY_upper, :] @ u)
327
328             displacement.append(dep)
329             force.append(f)
330
331             # check for damaged edges
332             if np.any(simu.damage[nodes_edges] >= 1):
333                 nDetect += 1
334                 if nDetect == 10:
335                     # If the edge has been touched 10 times, stop the simulation
336                     break
337
338         # ----------------------------------------------
339         # Saving
340         # ----------------------------------------------
341         force = np.asarray(force)
342         displacement = np.asarray(displacement)
343         print()
344         Simulations.Save_pickle(
345             (force, displacement), folder_save, "force-displacement"
346         )
347         simu.Save(folder_save)
348
349     else:
350         # ----------------------------------------------
351         # Loading
352         # ---------------------------------------------
353         simu: Simulations.PhaseField = Simulations.Load_Simu(folder_save)
354         force, displacement = Simulations.Load_pickle(folder_save, "force-displacement")
355
356     # ----------------------------------------------
357     # Results
358     # ---------------------------------------------
359     if plotResult:
360         Display.Plot_Result(
361             simu,
362             "damage",
363             nodeValues=True,
364             plotMesh=False,
365             folder=folder_save,
366             filename="damage",
367         )
368         Display.Plot_Mesh(simu)
369         Display.Plot_Iter_Summary(simu, folder_save, None, None)
370         Display.Plot_BoundaryConditions(simu)
371         Display.Plot_Force_Displacement(
372             force * 1e-6, displacement * 1e6, "ud [µm]", "f [kN/mm]", folder_save
373         )
374
375     if plotMesh:
376         Display.Plot_Mesh(simu.mesh)
377
378     if plotEnergy:
379         Display.Plot_Energy(simu, N=400, folder=folder_save)
380
381     if saveParaview:
382         Paraview.Save_simu(simu, folder_save, 400)
383
384     if makeMovie:
385         simu.Set_Iter(-1)
386         nodes_upper = simu.mesh.Nodes_Conditions(lambda x, y, z: y == L)
387         depMax = simu.Result("displacement_norm")[nodes_upper].max()
388         deformFactor = L * 0.05 / depMax
389
390         iterations = np.arange(0, simu.Niter, simu.Niter // 20)
391
392         def Func(plotter, iter):
393             simu.Set_Iter(iterations[iter])
394
395             grid = PyVista._pvMesh(simu, "damage", deformFactor)
396
397             tresh = grid.threshold((0, 0.8))
398
399             PyVista.Plot(
400                 tresh,
401                 "damage",
402                 deformFactor,
403                 plotMesh=True,
404                 plotter=plotter,
405                 clim=(0, 1),
406             )
407
408         PyVista.Movie_func(Func, iterations.size, folder_save, "damage.gif")
409
410     Tic.Resume()
411
412     if doSimu:
413         Tic.Plot_History(folder_save, False)
414
415     if showResult:
416         plt.show()
417
418     # Display.plt.close("all")
419     Tic.Clear()
420
421
422 if __name__ == "__main__":
423     # generates configs
424     Splits = []
425     Regus = []
426     for split in splits.copy():
427         for regu in regus.copy():
428             Splits.append(split)
429             Regus.append(regu)
430
431     if useParallel:
432         items = [(split, regu) for split, regu in zip(Splits, Regus)]
433         with multiprocessing.Pool(nProcs) as pool:
434             for result in pool.starmap(DoSimu, items):
435                 pass
436     else:
437         [DoSimu(split, regu) for split, regu in zip(Splits, Regus)]

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

Gallery generated by Sphinx-Gallery