.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/PhaseField/Shear.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code. .. rst-class:: sphx-glr-example-title .. _sphx_glr_examples_PhaseField_Shear.py: Shear ===== Damage simulation for a plate subjected to shear. .. GENERATED FROM PYTHON SOURCE LINES 13-367 .. image-sg:: /examples/PhaseField/images/sphx_glr_Shear_001.gif :alt: Shear :srcset: /examples/PhaseField/images/sphx_glr_Shear_001.gif :class: sphx-glr-single-img .. rst-class:: sphx-glr-horizontal * .. image-sg:: /examples/PhaseField/images/sphx_glr_Shear_002.png :alt: $\phi$ :srcset: /examples/PhaseField/images/sphx_glr_Shear_002.png :class: sphx-glr-multi-img * .. image-sg:: /examples/PhaseField/images/sphx_glr_Shear_003.png :alt: TRI3: Ne = 2923, Nn = 1540 :srcset: /examples/PhaseField/images/sphx_glr_Shear_003.png :class: sphx-glr-multi-img * .. image-sg:: /examples/PhaseField/images/sphx_glr_Shear_004.png :alt: Shear :srcset: /examples/PhaseField/images/sphx_glr_Shear_004.png :class: sphx-glr-multi-img * .. image-sg:: /examples/PhaseField/images/sphx_glr_Shear_005.png :alt: Boundary conditions :srcset: /examples/PhaseField/images/sphx_glr_Shear_005.png :class: sphx-glr-multi-img * .. image-sg:: /examples/PhaseField/images/sphx_glr_Shear_006.png :alt: Shear :srcset: /examples/PhaseField/images/sphx_glr_Shear_006.png :class: sphx-glr-multi-img * .. image-sg:: /examples/PhaseField/images/sphx_glr_Shear_007.png :alt: Summary :srcset: /examples/PhaseField/images/sphx_glr_Shear_007.png :class: sphx-glr-multi-img .. rst-class:: sphx-glr-script-out .. code-block:: none /home/docs/checkouts/readthedocs.org/user_builds/easyfea/envs/v1.5.5/lib/python3.11/site-packages/results/PhaseField/Shear_2D/Test/ElasIsot_Miehe_AT1_DP_optimMesh 1 : 0.050 µm, [0.00e+00; 0.00e+00], 1:0.151 s, tol=1.00e+00 2 : 0.100 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=1.00e+00 0.25 % -> 14.87 s 3 : 0.150 µm, [0.00e+00; 0.00e+00], 1:0.152 s, tol=5.56e-01 0.50 % -> 1.01 m 4 : 0.200 µm, [0.00e+00; 0.00e+00], 1:0.164 s, tol=4.38e-01 0.75 % -> 1.09 m 5 : 0.250 µm, [0.00e+00; 0.00e+00], 1:0.161 s, tol=3.60e-01 1.00 % -> 1.06 m 6 : 0.300 µm, [0.00e+00; 0.00e+00], 1:0.158 s, tol=3.06e-01 1.25 % -> 1.04 m 7 : 0.350 µm, [0.00e+00; 0.00e+00], 1:0.158 s, tol=2.65e-01 1.50 % -> 1.04 m 8 : 0.400 µm, [0.00e+00; 0.00e+00], 1:0.153 s, tol=2.34e-01 1.75 % -> 1.00 m 9 : 0.450 µm, [0.00e+00; 0.00e+00], 1:0.141 s, tol=2.10e-01 2.00 % -> 55.37 s 10 : 0.500 µm, [0.00e+00; 0.00e+00], 1:0.151 s, tol=1.90e-01 2.25 % -> 59.22 s 11 : 0.550 µm, [0.00e+00; 0.00e+00], 1:0.159 s, tol=1.74e-01 2.50 % -> 1.03 m 12 : 0.600 µm, [0.00e+00; 0.00e+00], 1:0.162 s, tol=1.60e-01 2.75 % -> 1.05 m 13 : 0.650 µm, [0.00e+00; 0.00e+00], 1:0.164 s, tol=1.48e-01 3.00 % -> 1.06 m 14 : 0.700 µm, [0.00e+00; 0.00e+00], 1:0.156 s, tol=1.38e-01 3.25 % -> 1.01 m 15 : 0.750 µm, [0.00e+00; 0.00e+00], 1:0.075 s, tol=1.29e-01 3.50 % -> 28.89 s 16 : 0.800 µm, [0.00e+00; 0.00e+00], 1:0.152 s, tol=1.21e-01 3.75 % -> 58.35 s 17 : 0.850 µm, [0.00e+00; 0.00e+00], 1:0.153 s, tol=1.14e-01 4.00 % -> 58.86 s 18 : 0.900 µm, [0.00e+00; 0.00e+00], 1:0.148 s, tol=1.08e-01 4.25 % -> 56.75 s 19 : 0.950 µm, [0.00e+00; 0.00e+00], 1:0.160 s, tol=1.02e-01 4.50 % -> 1.02 m 20 : 1.000 µm, [0.00e+00; 0.00e+00], 1:0.156 s, tol=9.75e-02 4.75 % -> 59.51 s 21 : 1.050 µm, [0.00e+00; 0.00e+00], 1:0.149 s, tol=9.30e-02 5.00 % -> 56.71 s 22 : 1.100 µm, [0.00e+00; 0.00e+00], 1:0.152 s, tol=8.88e-02 5.25 % -> 57.70 s 23 : 1.150 µm, [0.00e+00; 0.00e+00], 1:0.161 s, tol=8.51e-02 5.50 % -> 1.02 m 24 : 1.200 µm, [0.00e+00; 0.00e+00], 1:0.152 s, tol=8.16e-02 5.75 % -> 57.32 s 25 : 1.250 µm, [0.00e+00; 0.00e+00], 1:0.158 s, tol=7.84e-02 6.00 % -> 59.51 s 26 : 1.300 µm, [0.00e+00; 0.00e+00], 1:0.148 s, tol=7.54e-02 6.25 % -> 55.57 s 27 : 1.350 µm, [0.00e+00; 0.00e+00], 1:0.162 s, tol=7.27e-02 6.50 % -> 1.01 m 28 : 1.400 µm, [0.00e+00; 0.00e+00], 1:0.142 s, tol=7.02e-02 6.75 % -> 52.88 s 29 : 1.450 µm, [0.00e+00; 0.00e+00], 1:0.165 s, tol=6.78e-02 7.00 % -> 1.02 m 30 : 1.500 µm, [0.00e+00; 0.00e+00], 1:0.162 s, tol=6.56e-02 7.25 % -> 1.00 m 31 : 1.550 µm, [0.00e+00; 0.00e+00], 1:0.156 s, tol=6.35e-02 7.50 % -> 57.82 s 32 : 1.600 µm, [0.00e+00; 0.00e+00], 1:0.160 s, tol=6.15e-02 7.75 % -> 59.11 s 33 : 1.650 µm, [0.00e+00; 0.00e+00], 1:0.160 s, tol=5.97e-02 8.00 % -> 58.97 s 34 : 1.700 µm, [0.00e+00; 0.00e+00], 1:0.176 s, tol=5.80e-02 8.25 % -> 1.08 m 35 : 1.750 µm, [0.00e+00; 0.00e+00], 1:0.161 s, tol=5.63e-02 8.50 % -> 59.05 s 36 : 1.800 µm, [0.00e+00; 0.00e+00], 1:0.160 s, tol=5.48e-02 8.75 % -> 58.23 s 37 : 1.850 µm, [0.00e+00; 0.00e+00], 1:0.159 s, tol=5.33e-02 9.00 % -> 57.84 s 38 : 1.900 µm, [0.00e+00; 0.00e+00], 1:0.157 s, tol=5.19e-02 9.25 % -> 57.07 s 39 : 1.950 µm, [0.00e+00; 0.00e+00], 1:0.162 s, tol=5.06e-02 9.50 % -> 58.72 s 40 : 2.000 µm, [0.00e+00; 0.00e+00], 1:0.158 s, tol=4.94e-02 9.75 % -> 57.13 s 41 : 2.050 µm, [0.00e+00; 0.00e+00], 1:0.063 s, tol=4.82e-02 10.00 % -> 22.64 s 42 : 2.100 µm, [0.00e+00; 0.00e+00], 1:0.111 s, tol=4.71e-02 10.25 % -> 39.71 s 43 : 2.150 µm, [0.00e+00; 0.00e+00], 1:0.157 s, tol=4.60e-02 10.50 % -> 56.26 s 44 : 2.200 µm, [0.00e+00; 0.00e+00], 1:0.163 s, tol=4.49e-02 10.75 % -> 58.26 s 45 : 2.250 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=4.40e-02 11.00 % -> 55.25 s 46 : 2.300 µm, [0.00e+00; 0.00e+00], 1:0.157 s, tol=4.30e-02 11.25 % -> 55.81 s 47 : 2.350 µm, [0.00e+00; 0.00e+00], 1:0.164 s, tol=4.21e-02 11.50 % -> 58.13 s 48 : 2.400 µm, [0.00e+00; 0.00e+00], 1:0.162 s, tol=4.12e-02 11.75 % -> 57.26 s 49 : 2.450 µm, [0.00e+00; 0.00e+00], 1:0.164 s, tol=4.04e-02 12.00 % -> 57.82 s 50 : 2.500 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=3.96e-02 12.25 % -> 54.50 s 51 : 2.550 µm, [0.00e+00; 0.00e+00], 1:0.161 s, tol=3.88e-02 12.50 % -> 56.40 s 52 : 2.600 µm, [0.00e+00; 0.00e+00], 1:0.157 s, tol=3.81e-02 12.75 % -> 54.88 s 53 : 2.650 µm, [0.00e+00; 0.00e+00], 1:0.149 s, tol=3.74e-02 13.00 % -> 51.91 s 54 : 2.700 µm, [0.00e+00; 0.00e+00], 1:0.147 s, tol=3.67e-02 13.25 % -> 51.08 s 55 : 2.750 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=3.60e-02 13.50 % -> 53.71 s 56 : 2.800 µm, [0.00e+00; 0.00e+00], 1:0.154 s, tol=3.54e-02 13.75 % -> 53.21 s 57 : 2.850 µm, [0.00e+00; 0.00e+00], 1:0.150 s, tol=3.48e-02 14.00 % -> 51.68 s 58 : 2.900 µm, [0.00e+00; 0.00e+00], 1:0.154 s, tol=3.42e-02 14.25 % -> 52.92 s 59 : 2.950 µm, [0.00e+00; 0.00e+00], 1:0.159 s, tol=3.36e-02 14.50 % -> 54.41 s 60 : 3.000 µm, [0.00e+00; 0.00e+00], 1:0.158 s, tol=3.31e-02 14.75 % -> 53.95 s 61 : 3.050 µm, [0.00e+00; 0.00e+00], 1:0.165 s, tol=3.25e-02 15.00 % -> 56.19 s 62 : 3.100 µm, [0.00e+00; 0.00e+00], 1:0.151 s, tol=3.20e-02 15.25 % -> 51.35 s 63 : 3.150 µm, [0.00e+00; 0.00e+00], 1:0.162 s, tol=3.15e-02 15.50 % -> 54.72 s 64 : 3.200 µm, [0.00e+00; 0.00e+00], 1:0.160 s, tol=3.10e-02 15.75 % -> 53.99 s 65 : 3.250 µm, [0.00e+00; 0.00e+00], 1:0.152 s, tol=3.05e-02 16.00 % -> 51.15 s 66 : 3.300 µm, [0.00e+00; 0.00e+00], 1:0.153 s, tol=3.01e-02 16.25 % -> 51.34 s 67 : 3.350 µm, [0.00e+00; 0.00e+00], 1:0.098 s, tol=2.96e-02 16.50 % -> 32.66 s 68 : 3.400 µm, [0.00e+00; 0.00e+00], 1:0.158 s, tol=2.92e-02 16.75 % -> 52.45 s 69 : 3.450 µm, [0.00e+00; 0.00e+00], 1:0.158 s, tol=2.88e-02 17.00 % -> 52.54 s 70 : 3.500 µm, [0.00e+00; 0.00e+00], 1:0.161 s, tol=2.84e-02 17.25 % -> 53.39 s 71 : 3.550 µm, [0.00e+00; 0.00e+00], 1:0.157 s, tol=2.80e-02 17.50 % -> 51.88 s 72 : 3.600 µm, [0.00e+00; 0.00e+00], 1:0.159 s, tol=2.76e-02 17.75 % -> 52.36 s 73 : 3.650 µm, [0.00e+00; 0.00e+00], 1:0.159 s, tol=2.72e-02 18.00 % -> 52.23 s 74 : 3.700 µm, [0.00e+00; 0.00e+00], 1:0.160 s, tol=2.68e-02 18.25 % -> 52.37 s 75 : 3.750 µm, [0.00e+00; 0.00e+00], 1:0.150 s, tol=2.65e-02 18.50 % -> 48.98 s 76 : 3.800 µm, [0.00e+00; 0.00e+00], 1:0.163 s, tol=2.61e-02 18.75 % -> 53.05 s 77 : 3.850 µm, [0.00e+00; 0.00e+00], 1:0.160 s, tol=2.58e-02 19.00 % -> 51.92 s 78 : 3.900 µm, [0.00e+00; 0.00e+00], 1:0.159 s, tol=2.55e-02 19.25 % -> 51.42 s 79 : 3.950 µm, [0.00e+00; 0.00e+00], 1:0.158 s, tol=2.52e-02 19.50 % -> 50.96 s 80 : 4.000 µm, [0.00e+00; 0.00e+00], 1:0.089 s, tol=2.48e-02 19.75 % -> 28.44 s 81 : 4.050 µm, [0.00e+00; 0.00e+00], 1:0.156 s, tol=2.45e-02 20.00 % -> 49.85 s 82 : 4.100 µm, [0.00e+00; 0.00e+00], 1:0.161 s, tol=2.42e-02 20.25 % -> 51.42 s 83 : 4.150 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=2.40e-02 20.50 % -> 49.35 s 84 : 4.200 µm, [0.00e+00; 0.00e+00], 1:0.154 s, tol=2.37e-02 20.75 % -> 48.88 s 85 : 4.250 µm, [0.00e+00; 0.00e+00], 1:0.157 s, tol=2.34e-02 21.00 % -> 49.68 s 86 : 4.300 µm, [0.00e+00; 0.00e+00], 1:0.158 s, tol=2.31e-02 21.25 % -> 49.83 s 87 : 4.350 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=2.29e-02 21.50 % -> 48.77 s 88 : 4.400 µm, [0.00e+00; 0.00e+00], 1:0.160 s, tol=2.26e-02 21.75 % -> 50.13 s 89 : 4.450 µm, [0.00e+00; 0.00e+00], 1:0.159 s, tol=2.23e-02 22.00 % -> 49.66 s 90 : 4.500 µm, [0.00e+00; 0.00e+00], 1:0.152 s, tol=2.21e-02 22.25 % -> 47.34 s 91 : 4.550 µm, [0.00e+00; 0.00e+00], 1:0.160 s, tol=2.19e-02 22.50 % -> 49.67 s 92 : 4.600 µm, [0.00e+00; 0.00e+00], 1:0.160 s, tol=2.16e-02 22.75 % -> 49.53 s 93 : 4.650 µm, [0.00e+00; 0.00e+00], 1:0.099 s, tol=2.14e-02 23.00 % -> 30.61 s 94 : 4.700 µm, [0.00e+00; 0.00e+00], 1:0.033 s, tol=2.12e-02 23.25 % -> 9.99 s 95 : 4.750 µm, [0.00e+00; 0.00e+00], 1:0.035 s, tol=2.09e-02 23.50 % -> 10.79 s 96 : 4.800 µm, [0.00e+00; 0.00e+00], 1:0.035 s, tol=2.07e-02 23.75 % -> 10.76 s 97 : 4.850 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=2.05e-02 24.00 % -> 47.05 s 98 : 4.900 µm, [0.00e+00; 0.00e+00], 1:0.157 s, tol=2.03e-02 24.25 % -> 47.58 s 99 : 4.950 µm, [0.00e+00; 0.00e+00], 1:0.158 s, tol=2.01e-02 24.50 % -> 47.80 s 100 : 5.000 µm, [0.00e+00; 0.00e+00], 1:0.156 s, tol=1.99e-02 24.75 % -> 47.00 s 101 : 5.050 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=1.97e-02 25.00 % -> 46.56 s 102 : 5.100 µm, [0.00e+00; 0.00e+00], 1:0.152 s, tol=1.95e-02 25.25 % -> 45.53 s 103 : 5.150 µm, [0.00e+00; 0.00e+00], 1:0.161 s, tol=1.93e-02 25.50 % -> 48.05 s 104 : 5.200 µm, [0.00e+00; 0.00e+00], 1:0.153 s, tol=1.91e-02 25.75 % -> 45.48 s 105 : 5.250 µm, [0.00e+00; 0.00e+00], 1:0.160 s, tol=1.90e-02 26.00 % -> 47.44 s 106 : 5.300 µm, [0.00e+00; 0.00e+00], 1:0.161 s, tol=1.88e-02 26.25 % -> 47.56 s 107 : 5.350 µm, [0.00e+00; 0.00e+00], 1:0.148 s, tol=1.86e-02 26.50 % -> 43.55 s 108 : 5.400 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=1.84e-02 26.75 % -> 45.51 s 109 : 5.450 µm, [0.00e+00; 0.00e+00], 1:0.124 s, tol=1.83e-02 27.00 % -> 36.23 s 110 : 5.500 µm, [0.00e+00; 0.00e+00], 1:0.158 s, tol=1.81e-02 27.25 % -> 46.05 s 111 : 5.550 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=1.79e-02 27.50 % -> 44.99 s 112 : 5.600 µm, [0.00e+00; 0.00e+00], 1:0.153 s, tol=1.78e-02 27.75 % -> 44.30 s 113 : 5.650 µm, [0.00e+00; 0.00e+00], 1:0.162 s, tol=1.76e-02 28.00 % -> 46.71 s 114 : 5.700 µm, [0.00e+00; 0.00e+00], 1:0.166 s, tol=1.75e-02 28.25 % -> 47.61 s 115 : 5.750 µm, [0.00e+00; 0.00e+00], 1:0.167 s, tol=1.73e-02 28.50 % -> 47.89 s 116 : 5.800 µm, [0.00e+00; 0.00e+00], 1:0.154 s, tol=1.72e-02 28.75 % -> 43.94 s 117 : 5.850 µm, [0.00e+00; 0.00e+00], 1:0.162 s, tol=1.70e-02 29.00 % -> 46.05 s 118 : 5.900 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=1.69e-02 29.25 % -> 43.93 s 119 : 5.950 µm, [0.00e+00; 0.00e+00], 1:0.156 s, tol=1.67e-02 29.50 % -> 44.03 s 120 : 6.000 µm, [0.00e+00; 0.00e+00], 1:0.157 s, tol=1.66e-02 29.75 % -> 44.17 s 121 : 6.050 µm, [0.00e+00; 0.00e+00], 1:0.138 s, tol=1.65e-02 30.00 % -> 38.65 s 122 : 6.100 µm, [0.00e+00; 0.00e+00], 1:0.152 s, tol=1.63e-02 30.25 % -> 42.50 s 123 : 6.150 µm, [0.00e+00; 0.00e+00], 1:0.157 s, tol=1.62e-02 30.50 % -> 43.70 s 124 : 6.200 µm, [0.00e+00; 0.00e+00], 1:0.157 s, tol=1.61e-02 30.75 % -> 43.55 s 125 : 6.250 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=1.59e-02 31.00 % -> 42.83 s 126 : 6.300 µm, [0.00e+00; 0.00e+00], 1:0.148 s, tol=1.58e-02 31.25 % -> 40.75 s 127 : 6.350 µm, [0.00e+00; 0.00e+00], 1:0.157 s, tol=1.57e-02 31.50 % -> 43.09 s 128 : 6.400 µm, [0.00e+00; 0.00e+00], 1:0.161 s, tol=1.56e-02 31.75 % -> 43.99 s 129 : 6.450 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=1.54e-02 32.00 % -> 42.21 s 130 : 6.500 µm, [0.00e+00; 0.00e+00], 1:0.157 s, tol=1.53e-02 32.25 % -> 42.65 s 131 : 6.550 µm, [0.00e+00; 0.00e+00], 1:0.164 s, tol=1.52e-02 32.50 % -> 44.27 s 132 : 6.600 µm, [0.00e+00; 0.00e+00], 1:0.155 s, tol=1.51e-02 32.75 % -> 41.74 s 133 : 6.650 µm, [0.00e+00; 0.00e+00], 1:0.158 s, tol=1.50e-02 33.00 % -> 42.40 s 134 : 6.700 µm, [0.00e+00; 0.00e+00], 1:0.110 s, tol=1.49e-02 33.25 % -> 29.50 s 135 : 6.750 µm, [0.00e+00; 0.00e+00], 1:0.027 s, tol=1.48e-02 33.50 % -> 7.29 s 136 : 6.800 µm, [0.00e+00; 0.00e+00], 1:0.026 s, tol=1.47e-02 33.75 % -> 6.95 s 137 : 6.850 µm, [0.00e+00; 0.00e+00], 1:0.161 s, tol=1.45e-02 34.00 % -> 42.41 s 138 : 6.900 µm, [0.00e+00; 0.00e+00], 1:0.148 s, tol=1.44e-02 34.25 % -> 38.97 s 139 : 6.950 µm, [0.00e+00; 0.00e+00], 1:0.152 s, tol=1.43e-02 34.50 % -> 39.88 s 140 : 7.000 µm, [6.00e-07; 1.59e-03], 1:0.164 s, tol=1.42e-02 34.75 % -> 42.87 s 141 : 7.050 µm, [2.22e-06; 6.20e-03], 1:0.162 s, tol=1.39e-02 35.00 % -> 42.19 s 142 : 7.100 µm, [4.16e-06; 1.23e-02], 1:0.163 s, tol=1.38e-02 35.25 % -> 42.27 s 143 : 7.150 µm, [6.13e-06; 1.92e-02], 1:0.161 s, tol=1.37e-02 35.50 % -> 41.57 s 144 : 7.200 µm, [7.97e-06; 2.65e-02], 1:0.176 s, tol=1.35e-02 35.75 % -> 45.29 s 145 : 7.250 µm, [9.65e-06; 3.41e-02], 1:0.166 s, tol=1.35e-02 36.00 % -> 42.55 s 146 : 7.300 µm, [1.12e-05; 4.18e-02], 1:0.167 s, tol=1.34e-02 36.25 % -> 42.61 s 147 : 7.350 µm, [1.25e-05; 4.96e-02], 1:0.158 s, tol=1.33e-02 36.50 % -> 40.16 s 148 : 7.400 µm, [1.36e-05; 5.75e-02], 1:0.157 s, tol=1.32e-02 36.75 % -> 39.67 s 149 : 7.450 µm, [1.46e-05; 6.55e-02], 1:0.032 s, tol=1.31e-02 37.00 % -> 8.13 s 150 : 7.500 µm, [1.55e-05; 7.35e-02], 1:0.035 s, tol=1.30e-02 37.25 % -> 8.83 s 151 : 7.550 µm, [1.62e-05; 8.16e-02], 1:0.161 s, tol=1.29e-02 37.50 % -> 40.14 s 152 : 7.600 µm, [1.68e-05; 8.98e-02], 1:0.162 s, tol=1.28e-02 37.75 % -> 40.38 s 153 : 7.650 µm, [1.72e-05; 9.80e-02], 1:0.167 s, tol=1.28e-02 38.00 % -> 41.48 s 154 : 7.700 µm, [1.76e-05; 1.06e-01], 1:0.166 s, tol=1.27e-02 38.25 % -> 41.04 s 155 : 7.750 µm, [1.79e-05; 1.15e-01], 1:0.168 s, tol=1.26e-02 38.50 % -> 41.37 s 156 : 7.800 µm, [1.81e-05; 1.23e-01], 1:0.167 s, tol=1.25e-02 38.75 % -> 40.97 s 157 : 7.850 µm, [1.82e-05; 1.32e-01], 1:0.167 s, tol=1.24e-02 39.00 % -> 40.80 s 158 : 7.900 µm, [1.83e-05; 1.40e-01], 1:0.065 s, tol=1.24e-02 39.25 % -> 15.76 s 159 : 7.950 µm, [1.83e-05; 1.49e-01], 1:0.035 s, tol=1.23e-02 39.50 % -> 8.53 s 160 : 8.000 µm, [1.82e-05; 1.58e-01], 1:0.117 s, tol=1.22e-02 39.75 % -> 28.08 s 161 : 8.050 µm, [1.81e-05; 1.67e-01], 1:0.165 s, tol=1.21e-02 40.00 % -> 39.65 s 162 : 8.100 µm, [1.80e-05; 1.76e-01], 1:0.167 s, tol=1.21e-02 40.25 % -> 39.95 s 163 : 8.150 µm, [1.78e-05; 1.85e-01], 1:0.156 s, tol=1.20e-02 40.50 % -> 37.18 s 164 : 8.200 µm, [1.75e-05; 1.94e-01], 1:0.157 s, tol=1.19e-02 40.75 % -> 37.27 s 165 : 8.250 µm, [1.73e-05; 2.03e-01], 1:0.162 s, tol=1.18e-02 41.00 % -> 38.26 s 166 : 8.300 µm, [1.70e-05; 2.12e-01], 1:0.165 s, tol=1.18e-02 41.25 % -> 38.81 s 167 : 8.350 µm, [1.67e-05; 2.22e-01], 1:0.164 s, tol=1.17e-02 41.50 % -> 38.43 s 168 : 8.400 µm, [1.64e-05; 2.31e-01], 1:0.165 s, tol=1.16e-02 41.75 % -> 38.47 s 169 : 8.450 µm, [1.60e-05; 2.41e-01], 1:0.161 s, tol=1.16e-02 42.00 % -> 37.41 s 170 : 8.500 µm, [1.57e-05; 2.51e-01], 1:0.166 s, tol=1.15e-02 42.25 % -> 38.39 s 171 : 8.550 µm, [1.53e-05; 2.60e-01], 1:0.162 s, tol=1.14e-02 42.50 % -> 37.31 s 172 : 8.600 µm, [1.49e-05; 2.70e-01], 1:0.163 s, tol=1.14e-02 42.75 % -> 37.36 s 173 : 8.650 µm, [1.45e-05; 2.80e-01], 1:0.165 s, tol=1.13e-02 43.00 % -> 37.68 s 174 : 8.700 µm, [1.42e-05; 2.90e-01], 1:0.160 s, tol=1.13e-02 43.25 % -> 36.35 s 175 : 8.750 µm, [1.38e-05; 3.01e-01], 1:0.163 s, tol=1.12e-02 43.50 % -> 36.90 s 176 : 8.800 µm, [1.34e-05; 3.11e-01], 1:0.100 s, tol=1.11e-02 43.75 % -> 22.46 s 177 : 8.850 µm, [1.30e-05; 3.21e-01], 1:0.169 s, tol=1.11e-02 44.00 % -> 37.96 s 178 : 8.900 µm, [1.26e-05; 3.32e-01], 1:0.156 s, tol=1.10e-02 44.25 % -> 34.83 s 179 : 8.950 µm, [1.22e-05; 3.43e-01], 1:0.161 s, tol=1.10e-02 44.50 % -> 35.78 s 180 : 9.000 µm, [1.19e-05; 3.56e-01], 1:0.154 s, tol=1.09e-02 44.75 % -> 34.10 s 181 : 9.050 µm, [1.15e-05; 3.71e-01], 1:0.172 s, tol=1.08e-02 45.00 % -> 37.87 s 182 : 9.100 µm, [1.12e-05; 3.86e-01], 1:0.168 s, tol=1.07e-02 45.25 % -> 36.83 s 183 : 9.150 µm, [1.09e-05; 4.01e-01], 1:0.165 s, tol=1.07e-02 45.50 % -> 36.01 s 184 : 9.200 µm, [1.05e-05; 4.17e-01], 1:0.167 s, tol=1.06e-02 45.75 % -> 36.28 s 185 : 9.250 µm, [1.02e-05; 4.32e-01], 1:0.164 s, tol=1.06e-02 46.00 % -> 35.47 s 186 : 9.300 µm, [9.84e-06; 4.47e-01], 1:0.159 s, tol=1.06e-02 46.25 % -> 34.21 s 187 : 9.350 µm, [9.49e-06; 4.62e-01], 1:0.167 s, tol=1.05e-02 46.50 % -> 35.78 s 188 : 9.400 µm, [9.27e-06; 4.77e-01], 1:0.161 s, tol=1.04e-02 46.75 % -> 34.33 s 189 : 9.450 µm, [9.10e-06; 4.94e-01], 1:0.134 s, tol=1.03e-02 47.00 % -> 28.45 s 190 : 9.500 µm, [8.97e-06; 5.14e-01], 1:0.158 s, tol=1.03e-02 47.25 % -> 33.38 s 191 : 9.550 µm, [8.88e-06; 5.39e-01], 1:0.168 s, tol=1.02e-02 47.50 % -> 35.31 s 192 : 9.600 µm, [8.83e-06; 5.71e-01], 1:0.169 s, tol=1.03e-02 47.75 % -> 35.35 s 193 : 9.650 µm, [8.78e-06; 6.08e-01], 1:0.157 s, tol=1.04e-02 48.00 % -> 32.72 s 194 : 9.700 µm, [8.76e-06; 6.49e-01], 1:0.167 s, tol=1.07e-02 48.25 % -> 34.59 s 195 : 9.750 µm, [8.75e-06; 6.95e-01], 1:0.164 s, tol=1.13e-02 48.50 % -> 33.83 s 196 : 9.800 µm, [9.05e-06; 7.43e-01], 1:0.167 s, tol=1.23e-02 48.75 % -> 34.27 s 197 : 9.850 µm, [9.39e-06; 7.94e-01], 1:0.169 s, tol=1.47e-02 49.00 % -> 34.50 s 198 : 9.900 µm, [9.76e-06; 8.46e-01], 1:0.157 s, tol=1.92e-02 49.25 % -> 31.92 s 199 : 9.950 µm, [9.84e-06; 8.94e-01], 1:0.165 s, tol=2.57e-02 49.50 % -> 33.35 s 200 : 10.000 µm, [9.59e-06; 9.28e-01], 1:0.159 s, tol=3.11e-02 49.75 % -> 31.99 s 201 : 10.050 µm, [9.14e-06; 9.49e-01], 1:0.083 s, tol=3.15e-02 50.00 % -> 16.60 s 202 : 10.100 µm, [9.06e-06; 9.58e-01], 1:0.029 s, tol=2.98e-02 50.25 % -> 5.68 s 203 : 10.150 µm, [9.00e-06; 9.60e-01], 1:0.033 s, tol=3.01e-02 50.50 % -> 6.48 s 204 : 10.200 µm, [9.21e-06; 9.60e-01], 1:0.032 s, tol=3.37e-02 50.75 % -> 6.36 s 205 : 10.250 µm, [9.27e-06; 9.60e-01], 1:0.161 s, tol=4.01e-02 51.00 % -> 31.50 s 206 : 10.300 µm, [9.16e-06; 9.71e-01], 1:0.164 s, tol=4.63e-02 51.25 % -> 32.00 s 207 : 10.350 µm, [8.87e-06; 9.82e-01], 1:0.160 s, tol=4.81e-02 51.50 % -> 31.08 s 208 : 10.400 µm, [8.45e-06; 9.89e-01], 1:0.154 s, tol=4.32e-02 51.75 % -> 29.75 s 209 : 10.450 µm, [8.11e-06; 9.92e-01], 1:0.164 s, tol=3.51e-02 52.00 % -> 31.52 s 210 : 10.500 µm, [8.09e-06; 9.94e-01], 1:0.167 s, tol=3.08e-02 52.25 % -> 31.94 s 211 : 10.550 µm, [8.04e-06; 9.95e-01], 1:0.152 s, tol=2.91e-02 52.50 % -> 28.90 s 212 : 10.600 µm, [7.90e-06; 9.95e-01], 1:0.165 s, tol=2.80e-02 52.75 % -> 31.22 s 213 : 10.650 µm, [7.80e-06; 9.95e-01], 1:0.164 s, tol=2.75e-02 53.00 % -> 30.86 s 214 : 10.700 µm, [7.73e-06; 9.94e-01], 1:0.165 s, tol=2.85e-02 53.25 % -> 30.88 s 215 : 10.750 µm, [7.62e-06; 9.94e-01], 1:0.168 s, tol=3.07e-02 53.50 % -> 31.29 s 216 : 10.800 µm, [7.45e-06; 9.93e-01], 1:0.111 s, tol=3.36e-02 53.75 % -> 20.51 s 217 : 10.850 µm, [7.23e-06; 9.93e-01], 1:0.166 s, tol=3.59e-02 54.00 % -> 30.46 s 218 : 10.900 µm, [6.95e-06; 9.92e-01], 1:0.167 s, tol=3.61e-02 54.25 % -> 30.59 s 219 : 10.950 µm, [6.71e-06; 9.92e-01], 1:0.163 s, tol=3.49e-02 54.50 % -> 29.69 s 220 : 11.000 µm, [6.54e-06; 9.92e-01], 1:0.156 s, tol=3.27e-02 54.75 % -> 28.27 s 221 : 11.050 µm, [6.47e-06; 9.92e-01], 1:0.155 s, tol=3.08e-02 55.00 % -> 27.92 s 222 : 11.100 µm, [6.45e-06; 9.92e-01], 1:0.169 s, tol=3.00e-02 55.25 % -> 30.29 s 223 : 11.150 µm, [6.36e-06; 9.93e-01], 1:0.163 s, tol=2.96e-02 55.50 % -> 29.08 s 224 : 11.200 µm, [6.21e-06; 9.93e-01], 1:0.173 s, tol=2.93e-02 55.75 % -> 30.62 s 225 : 11.250 µm, [6.16e-06; 9.93e-01], 1:0.095 s, tol=3.12e-02 56.00 % -> 16.73 s 226 : 11.300 µm, [6.13e-06; 9.93e-01], 1:0.166 s, tol=3.48e-02 56.25 % -> 29.12 s 227 : 11.350 µm, [6.03e-06; 9.94e-01], 1:0.162 s, tol=3.87e-02 56.50 % -> 28.20 s 228 : 11.400 µm, [5.87e-06; 9.94e-01], 1:0.167 s, tol=4.11e-02 56.75 % -> 28.94 s 229 : 11.450 µm, [5.71e-06; 9.94e-01], 1:0.089 s, tol=4.15e-02 57.00 % -> 15.32 s 230 : 11.500 µm, [5.65e-06; 9.95e-01], 1:0.160 s, tol=4.01e-02 57.25 % -> 27.40 s 231 : 11.550 µm, [5.70e-06; 9.95e-01], 1:0.158 s, tol=3.84e-02 57.50 % -> 26.88 s 232 : 11.600 µm, [5.74e-06; 9.95e-01], 1:0.169 s, tol=3.71e-02 57.75 % -> 28.59 s 233 : 11.650 µm, [5.66e-06; 9.95e-01], 1:0.163 s, tol=3.53e-02 58.00 % -> 27.42 s 234 : 11.700 µm, [5.61e-06; 9.95e-01], 1:0.160 s, tol=3.48e-02 58.25 % -> 26.74 s 235 : 11.750 µm, [5.60e-06; 9.95e-01], 1:0.161 s, tol=3.73e-02 58.50 % -> 26.75 s 236 : 11.800 µm, [5.51e-06; 9.96e-01], 1:0.168 s, tol=4.00e-02 58.75 % -> 27.75 s 237 : 11.850 µm, [5.41e-06; 9.96e-01], 1:0.166 s, tol=4.25e-02 59.00 % -> 27.25 s 238 : 11.900 µm, [5.38e-06; 9.96e-01], 1:0.167 s, tol=4.26e-02 59.25 % -> 27.25 s 239 : 11.950 µm, [5.51e-06; 9.96e-01], 1:0.166 s, tol=4.20e-02 59.50 % -> 26.89 s 240 : 12.000 µm, [5.60e-06; 9.96e-01], 1:0.166 s, tol=4.03e-02 59.75 % -> 26.78 s 241 : 12.050 µm, [5.53e-06; 9.96e-01], 1:0.168 s, tol=3.75e-02 60.00 % -> 26.91 s 242 : 12.100 µm, [5.56e-06; 9.96e-01], 1:0.123 s, tol=3.75e-02 60.25 % -> 19.59 s 243 : 12.150 µm, [5.52e-06; 9.97e-01], 1:0.165 s, tol=3.99e-02 60.50 % -> 26.09 s 244 : 12.200 µm, [5.47e-06; 9.97e-01], 1:0.168 s, tol=4.29e-02 60.75 % -> 26.39 s 245 : 12.250 µm, [5.50e-06; 9.97e-01], 1:0.165 s, tol=4.42e-02 61.00 % -> 25.78 s 246 : 12.300 µm, [5.73e-06; 9.97e-01], 1:0.168 s, tol=4.43e-02 61.25 % -> 26.05 s 247 : 12.350 µm, [5.83e-06; 9.97e-01], 1:0.151 s, tol=4.23e-02 61.50 % -> 23.29 s 248 : 12.400 µm, [5.81e-06; 9.97e-01], 1:0.163 s, tol=3.95e-02 61.75 % -> 24.97 s 249 : 12.450 µm, [5.86e-06; 9.97e-01], 1:0.161 s, tol=4.02e-02 62.00 % -> 24.51 s 250 : 12.500 µm, [5.85e-06; 9.97e-01], 1:0.115 s, tol=4.32e-02 62.25 % -> 17.42 s 251 : 12.550 µm, [5.91e-06; 9.97e-01], 1:0.025 s, tol=4.52e-02 62.50 % -> 3.82 s 252 : 12.600 µm, [6.26e-06; 9.97e-01], 1:0.026 s, tol=4.57e-02 62.75 % -> 3.82 s 253 : 12.650 µm, [6.37e-06; 9.97e-01], 1:0.026 s, tol=4.31e-02 63.00 % -> 3.80 s 254 : 12.700 µm, [6.44e-06; 9.98e-01], 1:0.030 s, tol=4.07e-02 63.25 % -> 4.34 s 255 : 12.750 µm, [6.46e-06; 9.98e-01], 1:0.155 s, tol=4.15e-02 63.50 % -> 22.60 s 256 : 12.800 µm, [6.47e-06; 9.98e-01], 1:0.158 s, tol=4.43e-02 63.75 % -> 22.91 s 257 : 12.850 µm, [6.91e-06; 9.98e-01], 1:0.164 s, tol=4.65e-02 64.00 % -> 23.64 s 258 : 12.900 µm, [7.01e-06; 9.98e-01], 1:0.055 s, tol=4.45e-02 64.25 % -> 7.81 s 259 : 12.950 µm, [6.99e-06; 9.98e-01], 1:0.026 s, tol=4.13e-02 64.50 % -> 3.74 s 260 : 13.000 µm, [6.98e-06; 9.98e-01], 1:0.033 s, tol=4.14e-02 64.75 % -> 4.63 s 261 : 13.050 µm, [6.91e-06; 9.98e-01], 1:0.128 s, tol=4.41e-02 65.00 % -> 17.93 s 262 : 13.100 µm, [7.11e-06; 9.98e-01], 1:0.085 s, tol=4.64e-02 65.25 % -> 11.76 s 263 : 13.150 µm, [7.16e-06; 9.98e-01], 1:0.036 s, tol=4.35e-02 65.50 % -> 5.01 s 264 : 13.200 µm, [7.20e-06; 9.98e-01], 1:0.120 s, tol=4.08e-02 65.75 % -> 16.39 s 265 : 13.250 µm, [7.21e-06; 9.98e-01], 1:0.169 s, tol=4.22e-02 66.00 % -> 23.01 s 266 : 13.300 µm, [7.25e-06; 9.98e-01], 1:0.165 s, tol=4.48e-02 66.25 % -> 22.31 s 267 : 13.350 µm, [7.50e-06; 9.98e-01], 1:0.167 s, tol=4.48e-02 66.50 % -> 22.38 s 268 : 13.400 µm, [7.52e-06; 9.98e-01], 1:0.165 s, tol=4.09e-02 66.75 % -> 21.97 s 269 : 13.450 µm, [7.59e-06; 9.98e-01], 1:0.159 s, tol=4.04e-02 67.00 % -> 21.03 s 270 : 13.500 µm, [7.57e-06; 9.98e-01], 1:0.167 s, tol=4.28e-02 67.25 % -> 21.88 s 271 : 13.550 µm, [7.89e-06; 9.98e-01], 1:0.162 s, tol=4.45e-02 67.50 % -> 21.08 s 272 : 13.600 µm, [7.94e-06; 9.98e-01], 1:0.167 s, tol=4.04e-02 67.75 % -> 21.57 s 273 : 13.650 µm, [8.07e-06; 9.99e-01], 1:0.164 s, tol=3.91e-02 68.00 % -> 21.02 s 274 : 13.700 µm, [8.09e-06; 9.99e-01], 1:0.155 s, tol=4.15e-02 68.25 % -> 19.71 s 275 : 13.750 µm, [8.47e-06; 9.99e-01], 1:0.158 s, tol=4.34e-02 68.50 % -> 19.94 s 276 : 13.800 µm, [8.62e-06; 9.99e-01], 1:0.160 s, tol=3.95e-02 68.75 % -> 20.00 s 277 : 13.850 µm, [8.76e-06; 9.99e-01], 1:0.169 s, tol=3.84e-02 69.00 % -> 20.99 s 278 : 13.900 µm, [8.91e-06; 9.99e-01], 1:0.157 s, tol=4.11e-02 69.25 % -> 19.33 s 279 : 13.950 µm, [9.37e-06; 9.99e-01], 1:0.167 s, tol=4.18e-02 69.50 % -> 20.41 s 280 : 14.000 µm, [9.56e-06; 9.99e-01], 1:0.166 s, tol=3.76e-02 69.75 % -> 20.08 s 281 : 14.050 µm, [9.70e-06; 9.99e-01], 1:0.174 s, tol=3.77e-02 70.00 % -> 20.89 s 282 : 14.100 µm, [1.02e-05; 9.99e-01], 1:0.165 s, tol=4.09e-02 70.25 % -> 19.67 s 283 : 14.150 µm, [1.05e-05; 9.99e-01], 1:0.159 s, tol=3.84e-02 70.50 % -> 18.78 s 284 : 14.200 µm, [1.08e-05; 9.99e-01], 1:0.170 s, tol=3.55e-02 70.75 % -> 19.91 s 285 : 14.250 µm, [1.10e-05; 9.99e-01], 1:0.163 s, tol=3.73e-02 71.00 % -> 18.93 s 286 : 14.300 µm, [1.18e-05; 9.99e-01], 1:0.087 s, tol=3.94e-02 71.25 % -> 10.06 s 287 : 14.350 µm, [1.21e-05; 9.99e-01], 1:0.036 s, tol=3.55e-02 71.50 % -> 4.16 s 288 : 14.400 µm, [1.24e-05; 9.99e-01], 1:0.165 s, tol=3.51e-02 71.75 % -> 18.60 s 289 : 14.450 µm, [1.32e-05; 9.99e-01], 1:0.169 s, tol=3.82e-02 72.00 % -> 18.94 s 290 : 14.500 µm, [1.38e-05; 9.99e-01], 1:0.171 s, tol=3.54e-02 72.25 % -> 19.01 s 291 : 14.550 µm, [1.41e-05; 9.99e-01], 1:0.158 s, tol=3.32e-02 72.50 % -> 17.40 s 292 : 14.600 µm, [1.49e-05; 9.99e-01], 1:0.167 s, tol=3.57e-02 72.75 % -> 18.22 s 293 : 14.650 µm, [1.57e-05; 9.99e-01], 1:0.166 s, tol=3.55e-02 73.00 % -> 17.95 s 294 : 14.700 µm, [1.62e-05; 9.99e-01], 1:0.168 s, tol=3.20e-02 73.25 % -> 18.00 s 295 : 14.750 µm, [1.68e-05; 9.99e-01], 1:0.167 s, tol=3.34e-02 73.50 % -> 17.71 s 296 : 14.800 µm, [1.78e-05; 9.99e-01], 1:0.169 s, tol=3.49e-02 73.75 % -> 17.76 s 297 : 14.850 µm, [1.83e-05; 9.99e-01], 1:0.162 s, tol=3.16e-02 74.00 % -> 16.87 s 298 : 14.900 µm, [1.89e-05; 9.99e-01], 1:0.171 s, tol=3.20e-02 74.25 % -> 17.63 s 299 : 14.950 µm, [1.95e-05; 9.99e-01], 1:0.122 s, tol=3.36e-02 74.50 % -> 12.40 s 300 : 15.000 µm, [2.00e-05; 9.99e-01], 1:0.036 s, tol=3.07e-02 74.75 % -> 3.68 s 301 : 15.050 µm, [2.04e-05; 9.99e-01], 1:0.130 s, tol=3.08e-02 75.00 % -> 12.98 s 302 : 15.100 µm, [2.08e-05; 9.99e-01], 1:0.167 s, tol=3.21e-02 75.25 % -> 16.53 s 303 : 15.150 µm, [2.11e-05; 9.99e-01], 1:0.170 s, tol=2.95e-02 75.50 % -> 16.67 s 304 : 15.200 µm, [2.14e-05; 9.99e-01], 1:0.157 s, tol=2.96e-02 75.75 % -> 15.26 s 305 : 15.250 µm, [2.18e-05; 9.99e-01], 1:0.155 s, tol=3.05e-02 76.00 % -> 14.88 s 306 : 15.300 µm, [2.20e-05; 9.99e-01], 1:0.166 s, tol=2.84e-02 76.25 % -> 15.79 s 307 : 15.350 µm, [2.23e-05; 9.99e-01], 1:0.166 s, tol=2.87e-02 76.50 % -> 15.62 s 308 : 15.400 µm, [2.26e-05; 9.99e-01], 1:0.145 s, tol=2.87e-02 76.75 % -> 13.46 s 309 : 15.450 µm, [2.28e-05; 9.99e-01], 1:0.027 s, tol=2.73e-02 77.00 % -> 2.48 s 310 : 15.500 µm, [2.30e-05; 9.99e-01], 1:0.031 s, tol=2.78e-02 77.25 % -> 2.82 s 311 : 15.550 µm, [2.32e-05; 9.99e-01], 1:0.163 s, tol=2.69e-02 77.50 % -> 14.67 s 312 : 15.600 µm, [2.33e-05; 9.99e-01], 1:0.163 s, tol=2.61e-02 77.75 % -> 14.52 s 313 : 15.650 µm, [2.34e-05; 9.99e-01], 1:0.166 s, tol=2.67e-02 78.00 % -> 14.62 s 314 : 15.700 µm, [2.35e-05; 9.99e-01], 1:0.148 s, tol=2.52e-02 78.25 % -> 12.84 s 315 : 15.750 µm, [2.36e-05; 9.99e-01], 1:0.027 s, tol=2.48e-02 78.50 % -> 2.34 s 316 : 15.800 µm, [2.37e-05; 9.99e-01], 1:0.029 s, tol=2.54e-02 78.75 % -> 2.44 s 317 : 15.850 µm, [2.37e-05; 9.99e-01], 1:0.182 s, tol=2.39e-02 79.00 % -> 15.32 s 318 : 15.900 µm, [2.38e-05; 9.99e-01], 1:0.160 s, tol=2.32e-02 79.25 % -> 13.28 s 319 : 15.950 µm, [2.39e-05; 9.99e-01], 1:0.155 s, tol=2.41e-02 79.50 % -> 12.73 s 320 : 16.000 µm, [2.39e-05; 9.99e-01], 1:0.163 s, tol=2.23e-02 79.75 % -> 13.21 s 321 : 16.050 µm, [2.39e-05; 9.99e-01], 1:0.165 s, tol=2.20e-02 80.00 % -> 13.20 s 322 : 16.100 µm, [2.40e-05; 9.99e-01], 1:0.161 s, tol=2.28e-02 80.25 % -> 12.74 s 323 : 16.150 µm, [2.40e-05; 9.99e-01], 1:0.167 s, tol=2.10e-02 80.50 % -> 13.04 s 324 : 16.200 µm, [2.40e-05; 9.99e-01], 1:0.173 s, tol=2.09e-02 80.75 % -> 13.32 s 325 : 16.250 µm, [2.40e-05; 9.99e-01], 1:0.162 s, tol=2.12e-02 81.00 % -> 12.33 s 326 : 16.300 µm, [2.40e-05; 9.99e-01], 1:0.161 s, tol=1.93e-02 81.25 % -> 12.09 s 327 : 16.350 µm, [2.40e-05; 9.99e-01], 1:0.172 s, tol=1.93e-02 81.50 % -> 12.74 s 328 : 16.400 µm, [2.41e-05; 9.99e-01], 1:0.116 s, tol=1.98e-02 81.75 % -> 8.47 s 329 : 16.450 µm, [2.41e-05; 9.99e-01], 1:0.156 s, tol=1.81e-02 82.00 % -> 11.26 s 330 : 16.500 µm, [2.41e-05; 9.99e-01], 1:0.173 s, tol=1.66e-02 82.25 % -> 12.30 s 331 : 16.550 µm, [2.41e-05; 9.99e-01], 1:0.164 s, tol=1.73e-02 82.50 % -> 11.48 s 332 : 16.600 µm, [2.41e-05; 9.99e-01], 1:0.176 s, tol=1.72e-02 82.75 % -> 12.15 s 333 : 16.650 µm, [2.41e-05; 9.99e-01], 1:0.167 s, tol=1.66e-02 83.00 % -> 11.37 s 334 : 16.700 µm, [2.41e-05; 9.99e-01], 1:0.163 s, tol=1.63e-02 83.25 % -> 10.93 s 335 : 16.750 µm, [2.41e-05; 9.99e-01], 1:0.163 s, tol=1.71e-02 83.50 % -> 10.77 s 336 : 16.800 µm, [2.41e-05; 9.99e-01], 1:0.170 s, tol=1.88e-02 83.75 % -> 11.06 s 337 : 16.850 µm, [2.41e-05; 9.99e-01], 1:0.165 s, tol=2.21e-02 84.00 % -> 10.57 s 338 : 16.900 µm, [2.42e-05; 9.99e-01], 1:0.163 s, tol=2.61e-02 84.25 % -> 10.28 s 339 : 16.950 µm, [2.42e-05; 9.99e-01], 1:0.172 s, tol=3.05e-02 84.50 % -> 10.67 s 340 : 17.000 µm, [2.42e-05; 1.00e+00], 1:0.159 s, tol=3.18e-02 84.75 % -> 9.69 s 341 : 17.050 µm, [2.42e-05; 1.00e+00], 1:0.157 s, tol=3.25e-02 85.00 % -> 9.39 s 342 : 17.100 µm, [2.42e-05; 1.00e+00], 1:0.175 s, tol=3.29e-02 85.25 % -> 10.34 s 343 : 17.150 µm, [2.43e-05; 1.00e+00], 1:0.167 s, tol=2.48e-02 85.50 % -> 9.70 s 344 : 17.200 µm, [2.43e-05; 1.01e+00], 1:0.157 s, tol=1.33e-02 85.75 % -> 8.96 s 345 : 17.250 µm, [2.43e-05; 1.01e+00], 1:0.164 s, tol=8.31e-03 86.00 % -> 9.19 s 346 : 17.300 µm, [2.43e-05; 1.01e+00], 1:0.159 s, tol=6.99e-03 86.25 % -> 8.76 s 347 : 17.350 µm, [2.43e-05; 1.01e+00], 1:0.157 s, tol=6.52e-03 86.50 % -> 8.48 s 348 : 17.400 µm, [2.43e-05; 1.01e+00], 1:0.156 s, tol=6.25e-03 86.75 % -> 8.28 s 349 : 17.450 µm, [2.43e-05; 1.01e+00], 1:0.163 s, tol=6.08e-03 87.00 % -> 8.49 s 350 : 17.500 µm, [2.43e-05; 1.01e+00], 1:0.165 s, tol=5.96e-03 87.25 % -> 8.42 s 351 : 17.550 µm, [2.43e-05; 1.01e+00], 1:0.174 s, tol=5.88e-03 87.50 % -> 8.70 s 352 : 17.600 µm, [2.43e-05; 1.01e+00], 1:0.159 s, tol=5.83e-03 87.75 % -> 7.80 s 353 : 17.650 µm, [2.43e-05; 1.01e+00], 1:0.149 s, tol=5.79e-03 88.00 % -> 7.16 s 354 : 17.700 µm, [2.43e-05; 1.01e+00], 1:0.168 s, tol=5.76e-03 88.25 % -> 7.90 s 355 : 17.750 µm, [2.43e-05; 1.01e+00], 1:0.168 s, tol=5.74e-03 88.50 % -> 7.74 s 356 : 17.800 µm, [2.43e-05; 1.01e+00], 1:0.155 s, tol=5.71e-03 88.75 % -> 6.98 s 357 : 17.850 µm, [2.43e-05; 1.01e+00], 1:0.171 s, tol=5.69e-03 89.00 % -> 7.53 s 358 : 17.900 µm, [2.43e-05; 1.01e+00], 1:0.167 s, tol=5.67e-03 89.25 % -> 7.19 s 359 : 17.950 µm, [2.43e-05; 1.01e+00], 1:0.175 s, tol=5.65e-03 89.50 % -> 7.36 s 360 : 18.000 µm, [2.43e-05; 1.01e+00], 1:0.168 s, tol=5.64e-03 89.75 % -> 6.89 s 361 : 18.050 µm, [2.43e-05; 1.01e+00], 1:0.153 s, tol=5.62e-03 90.00 % -> 6.13 s 362 : 18.100 µm, [2.43e-05; 1.01e+00], 1:0.165 s, tol=5.60e-03 90.25 % -> 6.44 s 363 : 18.150 µm, [2.43e-05; 1.01e+00], 1:0.167 s, tol=5.59e-03 90.50 % -> 6.35 s 364 : 18.200 µm, [2.43e-05; 1.01e+00], 1:0.168 s, tol=5.58e-03 90.75 % -> 6.22 s 365 : 18.250 µm, [2.43e-05; 1.01e+00], 1:0.159 s, tol=5.57e-03 91.00 % -> 5.73 s 366 : 18.300 µm, [2.43e-05; 1.01e+00], 1:0.167 s, tol=5.55e-03 91.25 % -> 5.85 s 367 : 18.350 µm, [2.43e-05; 1.01e+00], 1:0.164 s, tol=5.54e-03 91.50 % -> 5.58 s 368 : 18.400 µm, [2.43e-05; 1.01e+00], 1:0.172 s, tol=5.52e-03 91.75 % -> 5.68 s Saved: /results/PhaseField/Shear_2D/Test/ElasIsot_Miehe_AT1_DP_optimMesh/force-displacement.pickle Saved: /results/PhaseField/Shear_2D/Test/ElasIsot_Miehe_AT1_DP_optimMesh/simulation.pickle Saved: /results/PhaseField/Shear_2D/Test/ElasIsot_Miehe_AT1_DP_optimMesh/summary.txt Generate movie 0/20 Generate movie 1/20 (5.00 %) 2.67 s Generate movie 2/20 (10.00 %) 2.53 s Generate movie 3/20 (15.00 %) 2.39 s Generate movie 4/20 (20.00 %) 2.24 s Generate movie 5/20 (25.00 %) 2.10 s Generate movie 6/20 (30.00 %) 1.95 s Generate movie 7/20 (35.00 %) 1.82 s Generate movie 8/20 (40.00 %) 1.68 s Generate movie 9/20 (45.00 %) 1.54 s Generate movie 10/20 (50.00 %) 1.40 s Generate movie 11/20 (55.00 %) 1.26 s Generate movie 12/20 (60.00 %) 1.12 s Generate movie 13/20 (65.00 %) 978.24 ms Generate movie 14/20 (70.00 %) 843.93 ms Generate movie 15/20 (75.00 %) 693.27 ms Generate movie 16/20 (80.00 %) 559.36 ms Generate movie 17/20 (85.00 %) 418.34 ms Generate movie 18/20 (90.00 %) 278.31 ms Generate movie 19/20 (95.00 %) 139.72 ms Generate movie 20/20 (100.00 %) 0.00 µs PostProcessing : 103.351 ms Display : 486.569 ms Mesh : 48.083 ms Boundary Conditions : 5.236 ms Matrix : 8.270 s Split : 2.370 s Solver : 47.404 s Resolution phase field : 54.538 s PyVista_Interface : 4.650 s | .. code-block:: Python :lineno-start: 14 from EasyFEA import ( Display, Folder, Models, plt, np, Tic, ElemType, Mesh, Simulations, Paraview, PyVista, ) from EasyFEA.Geoms import Point, Points, Domain, Line, Contour import multiprocessing # Display.Clear() useParallel = False nProcs = 4 # number of processes in parallel # ---------------------------------------------- # Configuration # ---------------------------------------------- dim = 2 # simu options doSimu = True meshTest = True openCrack = True optimMesh = True # outputs plotResult = True showResult = True plotMesh = False plotEnergy = False saveParaview = False makeMovie = True # phasefield maxIter = 1000 tolConv = 1e-0 # 1e-1, 1e-2, 1e-3 pfmSolver = Models.PhaseField.SolverType.History # splits = ["Bourdin","Amor","Miehe","Stress"] # Splits Isotropes # splits = ["He","AnisotStrain","AnisotStress","Zhang"] # Splits Anisotropes sans bourdin # splits = ["Bourdin","Amor","Miehe","Stress","He","AnisotStrain","AnisotStress","Zhang"] splits = ["Miehe"] # regus = ["AT1", "AT2"] regus = ["AT1"] # "AT1", "AT2" # ---------------------------------------------- # Mesh # ---------------------------------------------- L = 1e-3 # m l0 = 1e-5 thickness = 1 if dim == 2 else 0.1 / 1000 def DoMesh(split: str) -> Mesh: # meshSize clC = l0 * 2 if meshTest else l0 / 2 if optimMesh: # a coarser mesh can be used outside the refined zone clD = clC * 4 # refines the mesh in the area where the crack will propagate gap = L * 0.05 h = L if split == "Bourdin" else L / 2 + gap refineDomain = Domain(Point(L / 2 - gap, 0), Point(L, h, thickness), clC) else: clD = clC refineDomain = None # geom pt1 = Point() pt2 = Point(L) pt3 = Point(L, L) pt4 = Point(0, L) contour = Points([pt1, pt2, pt3, pt4], clD) if dim == 2: ptC1 = Point(0, L / 2, isOpen=openCrack) ptC2 = Point(L / 2, L / 2) cracks = [Line(ptC1, ptC2, clC, isOpen=openCrack)] if dim == 3: meshSize = clD if optimMesh else clC ptC1 = Point(0, L / 2, 0, isOpen=openCrack) ptC2 = Point(L / 2, L / 2, 0) ptC3 = Point(L / 2, L / 2, thickness) ptC4 = Point(0, L / 2, thickness, isOpen=openCrack) l1 = Line(ptC1, ptC2, meshSize, openCrack) l2 = Line(ptC2, ptC3, meshSize, False) l3 = Line(ptC3, ptC4, meshSize, openCrack) l4 = Line(ptC4, ptC1, meshSize, openCrack) cracks = [Contour([l1, l2, l3, l4], isOpen=openCrack)] # folder = Folder.Join("",results=True) # ax = Display.Init_Axes() # contour.Get_Contour().Plot(ax, color='k', plotPoints=False) # cracks[0].Plot(ax, color='k', plotPoints=False) # # if refineDomain != None: # # refineDomain.Plot(ax, color='k', plotPoints=False) # ax.axis('off') # Display.Save_fig(folder,"sample",True) if dim == 2: mesh = contour.Mesh_2D([], ElemType.TRI3, cracks, [refineDomain]) elif dim == 3: mesh = contour.Mesh_Extrude( [], [0, 0, thickness], [4], ElemType.PRISM6, cracks, [refineDomain], additionalLines=[l1], ) return mesh # ---------------------------------------------- # Simu # ---------------------------------------------- def DoSimu(split: str, regu: str): # Builds the path to the folder based on the problem data folder_save = Folder.PhaseField_Folder( f"Shear_{dim}D", "ElasIsot", split, regu, "DP", tolConv, pfmSolver, meshTest, optimMesh, not openCrack, ) Display.MyPrint(folder_save, "green") if doSimu: mesh = DoMesh(split) # Nodes recovery nodes_crack = mesh.Nodes_Conditions(lambda x, y, z: (y == L / 2) & (x <= L / 2)) nodes_upper = mesh.Nodes_Conditions(lambda x, y, z: y == L) nodes_lower = mesh.Nodes_Conditions(lambda x, y, z: y == 0) nodes_left = mesh.Nodes_Conditions(lambda x, y, z: (x == 0) & (y > 0) & (y < L)) nodes_right = mesh.Nodes_Conditions( lambda x, y, z: (x == L) & (y > 0) & (y < L) ) # Builds edge nodes nodes_edges = [] for nodes in [nodes_lower, nodes_right, nodes_upper]: nodes_edges.extend(nodes) # ---------------------------------------------- # Material # ---------------------------------------------- material = Models.ElasIsot( dim, E=210e9, v=0.3, planeStress=False, thickness=thickness ) Gc = 2.7e3 # J/m2 pfm = Models.PhaseField(material, split, regu, Gc, l0, pfmSolver) # ---------------------------------------------- # Boundary conditions # ---------------------------------------------- u_inc = 5e-8 if meshTest else 1e-8 N = 400 if meshTest else 2000 loadings = np.linspace(u_inc, u_inc * N, N, endpoint=True) config = f""" u_inc = {u_inc:.1e} N = {N} for iter, dep in enumerate(loadings): loadings = np.linspace(u_inc, u_inc*N, N, endpoint=True) if not openCrack: simu.add_dirichlet(nodes_crack, [1], ["d"], problemType="damage") simu.add_dirichlet(nodes_left, [0], ["y"]) simu.add_dirichlet(nodes_right, [0], ["y"]) simu.add_dirichlet(nodes_upper, [dep,0], ["x","y"]) simu.add_dirichlet(nodes_lower, [0]*dim, simu.Get_dofs()) """ def Loading(dep): """Boundary conditions""" simu.Bc_Init() if not openCrack: simu.add_dirichlet(nodes_crack, [1], ["d"], problemType="damage") simu.add_dirichlet(nodes_left, [0], ["y"]) simu.add_dirichlet(nodes_right, [0], ["y"]) simu.add_dirichlet(nodes_upper, [dep, 0], ["x", "y"]) simu.add_dirichlet(nodes_lower, [0] * dim, simu.Get_unknowns()) # ---------------------------------------------- # Simulation # ---------------------------------------------- simu = Simulations.PhaseFieldSimu(mesh, pfm, verbosity=False) simu.Results_Set_Bc_Summary(config) dofsX_upper = simu.Bc_dofs_nodes(nodes_upper, ["x"]) # INIT N = len(loadings) nDetect = 0 displacement = [] force = [] for iter, dep in enumerate(loadings): # apply new boundary conditions Loading(dep) # solve and save iter u, _, Kglob, converg = simu.Solve(tolConv, maxIter, convOption=2) simu.Save_Iter() # print iter solution simu.Results_Set_Iteration_Summary(iter, dep * 1e6, "µm", iter / N, True) # If the solver has not converged, stop the simulation. if not converg: break # resulting force on upper edge f = np.sum(Kglob[dofsX_upper, :] @ u) displacement.append(dep) force.append(f) # Detect damaged edges if np.any(simu.damage[nodes_edges] >= 1): nDetect += 1 if nDetect == 10: break # ---------------------------------------------- # Saving # ---------------------------------------------- force = np.asarray(force) displacement = np.asarray(displacement) print() Simulations.Save_pickle( (force, displacement), folder_save, "force-displacement" ) simu.Save(folder_save) else: simu: Simulations.PhaseFieldSimu = Simulations.Load_Simu(folder_save) force, displacement = Simulations.Load_pickle(folder_save, "force-displacement") # ---------------------------------------------- # Results # --------------------------------------------- if plotResult: ax = Display.Plot_Result( simu, "damage", nodeValues=True, plotMesh=False, folder=folder_save, filename="damage", ncolors=25, ) Display.Plot_Mesh(simu) Display.Plot_Iter_Summary(simu, folder_save, None, None) Display.Plot_BoundaryConditions(simu) Display.Plot_Force_Displacement( force * 1e-6, displacement * 1e6, "ud [µm]", "f [kN/mm]", folder_save ) # ax = Display.Plot_Result(simu, "damage", 1.5, ncolors=21, clim=(0,0.9)) # ax.axis('off'); ax.set_title("") # Display.Save_fig(folder, "deform damage") if plotMesh: # PyVista.Plot_Mesh(simu.mesh).show() ax = Display.Plot_Mesh(simu.mesh, lw=0.3, facecolors="white") ax.axis("off") ax.set_title("") Display.Save_fig(folder_save, "mesh", transparent=True) if saveParaview: Paraview.Save_simu(simu, folder_save, 400) if makeMovie: simu.Set_Iter(-1) nodes_upper = simu.mesh.Nodes_Conditions(lambda x, y, z: y == L) depMax = simu.Result("displacement_norm")[nodes_upper].max() deformFactor = L * 0.05 / depMax iterations = np.arange(0, simu.Niter, simu.Niter // 20) def Func(plotter, iter): simu.Set_Iter(iterations[iter]) grid = PyVista._pvMesh(simu, "damage", deformFactor) tresh = grid.threshold((0, 0.8)) PyVista.Plot( tresh, "damage", deformFactor, plotMesh=True, plotter=plotter, clim=(0, 1), ) PyVista.Movie_func(Func, iterations.size, folder_save, "damage.gif") if plotEnergy: Display.Plot_Energy(simu, N=400, folder=folder_save) Tic.Resume() if doSimu: Tic.Plot_History(folder_save, False) if showResult: plt.show() # plt.close("all") Tic.Clear() if __name__ == "__main__": # generates configs Splits = [] Regus = [] for split in splits.copy(): for regu in regus.copy(): Splits.append(split) Regus.append(regu) if useParallel: items = [(split, regu) for split, regu in zip(Splits, Regus)] with multiprocessing.Pool(nProcs) as pool: for result in pool.starmap(DoSimu, items): pass else: [DoSimu(split, regu) for split, regu in zip(Splits, Regus)] .. rst-class:: sphx-glr-timing **Total running time of the script:** (1 minutes 0.288 seconds) .. _sphx_glr_download_examples_PhaseField_Shear.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: Shear.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: Shear.py ` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: Shear.zip ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_