.. 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 14-374 .. 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/checkouts/v1.7.3/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.42 s 3 : 0.150 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=5.56e-01 0.50 % -> 14.34 s 4 : 0.200 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=4.38e-01 0.75 % -> 14.40 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.21 s 7 : 0.350 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.65e-01 1.50 % -> 14.38 s 8 : 0.400 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.34e-01 1.75 % -> 14.33 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.15 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.09 s 13 : 0.650 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.48e-01 3.00 % -> 13.97 s 14 : 0.700 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.38e-01 3.25 % -> 13.96 s 15 : 0.750 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=1.29e-01 3.50 % -> 14.11 s 16 : 0.800 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.21e-01 3.75 % -> 13.94 s 17 : 0.850 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.14e-01 4.00 % -> 13.91 s 18 : 0.900 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.08e-01 4.25 % -> 13.85 s 19 : 0.950 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.02e-01 4.50 % -> 13.78 s 20 : 1.000 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=9.75e-02 4.75 % -> 13.72 s 21 : 1.050 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=9.30e-02 5.00 % -> 13.73 s 22 : 1.100 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=8.88e-02 5.25 % -> 13.61 s 23 : 1.150 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=8.51e-02 5.50 % -> 13.66 s 24 : 1.200 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=8.16e-02 5.75 % -> 13.68 s 25 : 1.250 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=7.84e-02 6.00 % -> 13.57 s 26 : 1.300 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=7.54e-02 6.25 % -> 13.58 s 27 : 1.350 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=7.27e-02 6.50 % -> 13.55 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.37 s 30 : 1.500 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=6.56e-02 7.25 % -> 13.36 s 31 : 1.550 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=6.35e-02 7.50 % -> 13.38 s 32 : 1.600 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=6.15e-02 7.75 % -> 13.36 s 33 : 1.650 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=5.97e-02 8.00 % -> 13.44 s 34 : 1.700 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=5.80e-02 8.25 % -> 13.54 s 35 : 1.750 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=5.63e-02 8.50 % -> 13.45 s 36 : 1.800 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=5.48e-02 8.75 % -> 13.54 s 37 : 1.850 µm, [0.00e+00; 0.00e+00], 1:0.038 s, tol=5.33e-02 9.00 % -> 13.92 s 38 : 1.900 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=5.19e-02 9.25 % -> 13.53 s 39 : 1.950 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=5.06e-02 9.50 % -> 13.31 s 40 : 2.000 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=4.94e-02 9.75 % -> 13.11 s 41 : 2.050 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=4.82e-02 10.00 % -> 13.12 s 42 : 2.100 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=4.71e-02 10.25 % -> 13.16 s 43 : 2.150 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=4.60e-02 10.50 % -> 13.10 s 44 : 2.200 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=4.49e-02 10.75 % -> 13.09 s 45 : 2.250 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=4.40e-02 11.00 % -> 13.17 s 46 : 2.300 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=4.30e-02 11.25 % -> 12.87 s 47 : 2.350 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=4.21e-02 11.50 % -> 12.82 s 48 : 2.400 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=4.12e-02 11.75 % -> 12.91 s 49 : 2.450 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=4.04e-02 12.00 % -> 12.71 s 50 : 2.500 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=3.96e-02 12.25 % -> 12.74 s 51 : 2.550 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=3.88e-02 12.50 % -> 12.74 s 52 : 2.600 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=3.81e-02 12.75 % -> 12.62 s 53 : 2.650 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=3.74e-02 13.00 % -> 12.61 s 54 : 2.700 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=3.67e-02 13.25 % -> 12.53 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.038 s, tol=3.54e-02 13.75 % -> 13.28 s 57 : 2.850 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=3.48e-02 14.00 % -> 12.56 s 58 : 2.900 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=3.42e-02 14.25 % -> 12.42 s 59 : 2.950 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=3.36e-02 14.50 % -> 12.33 s 60 : 3.000 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=3.31e-02 14.75 % -> 12.35 s 61 : 3.050 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=3.25e-02 15.00 % -> 12.50 s 62 : 3.100 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=3.20e-02 15.25 % -> 12.38 s 63 : 3.150 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=3.15e-02 15.50 % -> 12.24 s 64 : 3.200 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=3.10e-02 15.75 % -> 12.15 s 65 : 3.250 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=3.05e-02 16.00 % -> 12.22 s 66 : 3.300 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=3.01e-02 16.25 % -> 12.18 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.036 s, tol=2.92e-02 16.75 % -> 12.07 s 69 : 3.450 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=2.88e-02 17.00 % -> 12.27 s 70 : 3.500 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.84e-02 17.25 % -> 12.01 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.95 s 73 : 3.650 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.72e-02 18.00 % -> 11.86 s 74 : 3.700 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.68e-02 18.25 % -> 11.84 s 75 : 3.750 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.65e-02 18.50 % -> 11.83 s 76 : 3.800 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=2.61e-02 18.75 % -> 12.05 s 77 : 3.850 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=2.58e-02 19.00 % -> 12.02 s 78 : 3.900 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=2.55e-02 19.25 % -> 11.87 s 79 : 3.950 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=2.52e-02 19.50 % -> 11.83 s 80 : 4.000 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=2.48e-02 19.75 % -> 11.87 s 81 : 4.050 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=2.45e-02 20.00 % -> 11.76 s 82 : 4.100 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=2.42e-02 20.25 % -> 11.78 s 83 : 4.150 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.40e-02 20.50 % -> 11.59 s 84 : 4.200 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.37e-02 20.75 % -> 11.46 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.46 s 87 : 4.350 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=2.29e-02 21.50 % -> 11.49 s 88 : 4.400 µm, [0.00e+00; 0.00e+00], 1:0.039 s, tol=2.26e-02 21.75 % -> 12.24 s 89 : 4.450 µm, [0.00e+00; 0.00e+00], 1:0.039 s, tol=2.23e-02 22.00 % -> 12.11 s 90 : 4.500 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=2.21e-02 22.25 % -> 11.45 s 91 : 4.550 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.19e-02 22.50 % -> 11.27 s 92 : 4.600 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.16e-02 22.75 % -> 11.19 s 93 : 4.650 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.14e-02 23.00 % -> 11.15 s 94 : 4.700 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.12e-02 23.25 % -> 11.15 s 95 : 4.750 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.09e-02 23.50 % -> 11.17 s 96 : 4.800 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.07e-02 23.75 % -> 11.10 s 97 : 4.850 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=2.05e-02 24.00 % -> 11.03 s 98 : 4.900 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=2.03e-02 24.25 % -> 11.08 s 99 : 4.950 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=2.01e-02 24.50 % -> 11.13 s 100 : 5.000 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.99e-02 24.75 % -> 10.95 s 101 : 5.050 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.97e-02 25.00 % -> 10.89 s 102 : 5.100 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.95e-02 25.25 % -> 10.83 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.81 s 105 : 5.250 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=1.90e-02 26.00 % -> 10.91 s 106 : 5.300 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=1.88e-02 26.25 % -> 10.78 s 107 : 5.350 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=1.86e-02 26.50 % -> 10.90 s 108 : 5.400 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.84e-02 26.75 % -> 10.67 s 109 : 5.450 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.83e-02 27.00 % -> 10.57 s 110 : 5.500 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.81e-02 27.25 % -> 10.51 s 111 : 5.550 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.79e-02 27.50 % -> 10.48 s 112 : 5.600 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.78e-02 27.75 % -> 10.46 s 113 : 5.650 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.76e-02 28.00 % -> 10.41 s 114 : 5.700 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.75e-02 28.25 % -> 10.39 s 115 : 5.750 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.73e-02 28.50 % -> 10.43 s 116 : 5.800 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=1.72e-02 28.75 % -> 10.46 s 117 : 5.850 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.70e-02 29.00 % -> 10.33 s 118 : 5.900 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.69e-02 29.25 % -> 10.22 s 119 : 5.950 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.67e-02 29.50 % -> 10.23 s 120 : 6.000 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.66e-02 29.75 % -> 10.22 s 121 : 6.050 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.65e-02 30.00 % -> 10.16 s 122 : 6.100 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=1.63e-02 30.25 % -> 10.21 s 123 : 6.150 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=1.62e-02 30.50 % -> 10.16 s 124 : 6.200 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.61e-02 30.75 % -> 10.08 s 125 : 6.250 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=1.59e-02 31.00 % -> 10.15 s 126 : 6.300 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.58e-02 31.25 % -> 10.02 s 127 : 6.350 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.57e-02 31.50 % -> 9.94 s 128 : 6.400 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=1.56e-02 31.75 % -> 10.00 s 129 : 6.450 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.54e-02 32.00 % -> 9.89 s 130 : 6.500 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.53e-02 32.25 % -> 9.86 s 131 : 6.550 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=1.52e-02 32.50 % -> 9.92 s 132 : 6.600 µm, [0.00e+00; 0.00e+00], 1:0.037 s, tol=1.51e-02 32.75 % -> 9.87 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.67 s 135 : 6.750 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.48e-02 33.50 % -> 9.68 s 136 : 6.800 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.47e-02 33.75 % -> 9.66 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.036 s, tol=1.44e-02 34.25 % -> 9.52 s 139 : 6.950 µm, [0.00e+00; 0.00e+00], 1:0.036 s, tol=1.43e-02 34.50 % -> 9.56 s 140 : 7.000 µm, [6.00e-07; 1.59e-03], 1:0.037 s, tol=1.42e-02 34.75 % -> 9.61 s 141 : 7.050 µm, [2.22e-06; 6.20e-03], 1:0.037 s, tol=1.39e-02 35.00 % -> 9.51 s 142 : 7.100 µm, [4.16e-06; 1.23e-02], 1:0.038 s, tol=1.38e-02 35.25 % -> 9.74 s 143 : 7.150 µm, [6.13e-06; 1.92e-02], 1:0.037 s, tol=1.37e-02 35.50 % -> 9.54 s 144 : 7.200 µm, [7.97e-06; 2.65e-02], 1:0.037 s, tol=1.35e-02 35.75 % -> 9.39 s 145 : 7.250 µm, [9.65e-06; 3.41e-02], 1:0.036 s, tol=1.35e-02 36.00 % -> 9.34 s 146 : 7.300 µm, [1.12e-05; 4.18e-02], 1:0.037 s, tol=1.34e-02 36.25 % -> 9.44 s 147 : 7.350 µm, [1.25e-05; 4.96e-02], 1:0.037 s, tol=1.33e-02 36.50 % -> 9.28 s 148 : 7.400 µm, [1.36e-05; 5.75e-02], 1:0.036 s, tol=1.32e-02 36.75 % -> 9.21 s 149 : 7.450 µm, [1.46e-05; 6.55e-02], 1:0.037 s, tol=1.31e-02 37.00 % -> 9.38 s 150 : 7.500 µm, [1.55e-05; 7.35e-02], 1:0.037 s, tol=1.30e-02 37.25 % -> 9.21 s 151 : 7.550 µm, [1.62e-05; 8.16e-02], 1:0.037 s, tol=1.29e-02 37.50 % -> 9.15 s 152 : 7.600 µm, [1.68e-05; 8.98e-02], 1:0.037 s, tol=1.28e-02 37.75 % -> 9.12 s 153 : 7.650 µm, [1.72e-05; 9.80e-02], 1:0.037 s, tol=1.28e-02 38.00 % -> 9.13 s 154 : 7.700 µm, [1.76e-05; 1.06e-01], 1:0.037 s, tol=1.27e-02 38.25 % -> 9.03 s 155 : 7.750 µm, [1.79e-05; 1.15e-01], 1:0.037 s, tol=1.26e-02 38.50 % -> 8.99 s 156 : 7.800 µm, [1.81e-05; 1.23e-01], 1:0.036 s, tol=1.25e-02 38.75 % -> 8.92 s 157 : 7.850 µm, [1.82e-05; 1.32e-01], 1:0.037 s, tol=1.24e-02 39.00 % -> 8.96 s 158 : 7.900 µm, [1.83e-05; 1.40e-01], 1:0.037 s, tol=1.24e-02 39.25 % -> 8.89 s 159 : 7.950 µm, [1.83e-05; 1.49e-01], 1:0.037 s, tol=1.23e-02 39.50 % -> 8.90 s 160 : 8.000 µm, [1.82e-05; 1.58e-01], 1:0.037 s, tol=1.22e-02 39.75 % -> 8.82 s 161 : 8.050 µm, [1.81e-05; 1.67e-01], 1:0.036 s, tol=1.21e-02 40.00 % -> 8.74 s 162 : 8.100 µm, [1.80e-05; 1.76e-01], 1:0.037 s, tol=1.21e-02 40.25 % -> 8.80 s 163 : 8.150 µm, [1.78e-05; 1.85e-01], 1:0.037 s, tol=1.20e-02 40.50 % -> 8.71 s 164 : 8.200 µm, [1.75e-05; 1.94e-01], 1:0.036 s, tol=1.19e-02 40.75 % -> 8.61 s 165 : 8.250 µm, [1.73e-05; 2.03e-01], 1:0.036 s, tol=1.18e-02 41.00 % -> 8.54 s 166 : 8.300 µm, [1.70e-05; 2.12e-01], 1:0.037 s, tol=1.18e-02 41.25 % -> 8.61 s 167 : 8.350 µm, [1.67e-05; 2.22e-01], 1:0.036 s, tol=1.17e-02 41.50 % -> 8.52 s 168 : 8.400 µm, [1.64e-05; 2.31e-01], 1:0.037 s, tol=1.16e-02 41.75 % -> 8.51 s 169 : 8.450 µm, [1.60e-05; 2.41e-01], 1:0.037 s, tol=1.16e-02 42.00 % -> 8.53 s 170 : 8.500 µm, [1.57e-05; 2.51e-01], 1:0.036 s, tol=1.15e-02 42.25 % -> 8.42 s 171 : 8.550 µm, [1.53e-05; 2.60e-01], 1:0.036 s, tol=1.14e-02 42.50 % -> 8.34 s 172 : 8.600 µm, [1.49e-05; 2.70e-01], 1:0.037 s, tol=1.14e-02 42.75 % -> 8.36 s 173 : 8.650 µm, [1.45e-05; 2.80e-01], 1:0.037 s, tol=1.13e-02 43.00 % -> 8.33 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.32 s 176 : 8.800 µm, [1.34e-05; 3.11e-01], 1:0.037 s, tol=1.11e-02 43.75 % -> 8.41 s 177 : 8.850 µm, [1.30e-05; 3.21e-01], 1:0.037 s, tol=1.11e-02 44.00 % -> 8.33 s 178 : 8.900 µm, [1.26e-05; 3.32e-01], 1:0.037 s, tol=1.10e-02 44.25 % -> 8.30 s 179 : 8.950 µm, [1.22e-05; 3.43e-01], 1:0.037 s, tol=1.10e-02 44.50 % -> 8.19 s 180 : 9.000 µm, [1.19e-05; 3.56e-01], 1:0.037 s, tol=1.09e-02 44.75 % -> 8.18 s 181 : 9.050 µm, [1.15e-05; 3.71e-01], 1:0.037 s, tol=1.08e-02 45.00 % -> 8.09 s 182 : 9.100 µm, [1.12e-05; 3.86e-01], 1:0.037 s, tol=1.07e-02 45.25 % -> 8.08 s 183 : 9.150 µm, [1.09e-05; 4.01e-01], 1:0.037 s, tol=1.07e-02 45.50 % -> 8.16 s 184 : 9.200 µm, [1.05e-05; 4.17e-01], 1:0.038 s, tol=1.06e-02 45.75 % -> 8.19 s 185 : 9.250 µm, [1.02e-05; 4.32e-01], 1:0.039 s, tol=1.06e-02 46.00 % -> 8.34 s 186 : 9.300 µm, [9.84e-06; 4.47e-01], 1:0.037 s, tol=1.06e-02 46.25 % -> 7.98 s 187 : 9.350 µm, [9.49e-06; 4.62e-01], 1:0.038 s, tol=1.05e-02 46.50 % -> 8.09 s 188 : 9.400 µm, [9.27e-06; 4.77e-01], 1:0.038 s, tol=1.04e-02 46.75 % -> 8.12 s 189 : 9.450 µm, [9.10e-06; 4.94e-01], 1:0.038 s, tol=1.03e-02 47.00 % -> 8.13 s 190 : 9.500 µm, [8.97e-06; 5.14e-01], 1:0.037 s, tol=1.03e-02 47.25 % -> 7.90 s 191 : 9.550 µm, [8.88e-06; 5.39e-01], 1:0.037 s, tol=1.02e-02 47.50 % -> 7.81 s 192 : 9.600 µm, [8.83e-06; 5.71e-01], 1:0.037 s, tol=1.03e-02 47.75 % -> 7.77 s 193 : 9.650 µm, [8.78e-06; 6.08e-01], 1:0.037 s, tol=1.04e-02 48.00 % -> 7.76 s 194 : 9.700 µm, [8.76e-06; 6.49e-01], 1:0.038 s, tol=1.07e-02 48.25 % -> 7.78 s 195 : 9.750 µm, [8.75e-06; 6.95e-01], 1:0.038 s, tol=1.13e-02 48.50 % -> 7.91 s 196 : 9.800 µm, [9.05e-06; 7.43e-01], 1:0.038 s, tol=1.23e-02 48.75 % -> 7.82 s 197 : 9.850 µm, [9.39e-06; 7.94e-01], 1:0.038 s, tol=1.47e-02 49.00 % -> 7.82 s 198 : 9.900 µm, [9.76e-06; 8.46e-01], 1:0.038 s, tol=1.92e-02 49.25 % -> 7.70 s 199 : 9.950 µm, [9.84e-06; 8.94e-01], 1:0.038 s, tol=2.57e-02 49.50 % -> 7.70 s 200 : 10.000 µm, [9.59e-06; 9.28e-01], 1:0.038 s, tol=3.11e-02 49.75 % -> 7.64 s 201 : 10.050 µm, [9.14e-06; 9.49e-01], 1:0.038 s, tol=3.15e-02 50.00 % -> 7.58 s 202 : 10.100 µm, [9.06e-06; 9.58e-01], 1:0.039 s, tol=2.98e-02 50.25 % -> 7.84 s 203 : 10.150 µm, [9.00e-06; 9.60e-01], 1:0.038 s, tol=3.01e-02 50.50 % -> 7.54 s 204 : 10.200 µm, [9.21e-06; 9.60e-01], 1:0.037 s, tol=3.37e-02 50.75 % -> 7.38 s 205 : 10.250 µm, [9.27e-06; 9.60e-01], 1:0.038 s, tol=4.01e-02 51.00 % -> 7.41 s 206 : 10.300 µm, [9.16e-06; 9.71e-01], 1:0.037 s, tol=4.63e-02 51.25 % -> 7.30 s 207 : 10.350 µm, [8.87e-06; 9.82e-01], 1:0.037 s, tol=4.81e-02 51.50 % -> 7.18 s 208 : 10.400 µm, [8.45e-06; 9.89e-01], 1:0.037 s, tol=4.32e-02 51.75 % -> 7.18 s 209 : 10.450 µm, [8.11e-06; 9.92e-01], 1:0.037 s, tol=3.51e-02 52.00 % -> 7.16 s 210 : 10.500 µm, [8.09e-06; 9.94e-01], 1:0.038 s, tol=3.08e-02 52.25 % -> 7.21 s 211 : 10.550 µm, [8.04e-06; 9.95e-01], 1:0.038 s, tol=2.91e-02 52.50 % -> 7.21 s 212 : 10.600 µm, [7.90e-06; 9.95e-01], 1:0.037 s, tol=2.80e-02 52.75 % -> 7.06 s 213 : 10.650 µm, [7.80e-06; 9.95e-01], 1:0.038 s, tol=2.75e-02 53.00 % -> 7.20 s 214 : 10.700 µm, [7.73e-06; 9.94e-01], 1:0.039 s, tol=2.85e-02 53.25 % -> 7.20 s 215 : 10.750 µm, [7.62e-06; 9.94e-01], 1:0.038 s, tol=3.07e-02 53.50 % -> 7.11 s 216 : 10.800 µm, [7.45e-06; 9.93e-01], 1:0.038 s, tol=3.36e-02 53.75 % -> 7.04 s 217 : 10.850 µm, [7.23e-06; 9.93e-01], 1:0.038 s, tol=3.59e-02 54.00 % -> 6.93 s 218 : 10.900 µm, [6.95e-06; 9.92e-01], 1:0.037 s, tol=3.61e-02 54.25 % -> 6.81 s 219 : 10.950 µm, [6.71e-06; 9.92e-01], 1:0.037 s, tol=3.49e-02 54.50 % -> 6.78 s 220 : 11.000 µm, [6.54e-06; 9.92e-01], 1:0.038 s, tol=3.27e-02 54.75 % -> 6.83 s 221 : 11.050 µm, [6.47e-06; 9.92e-01], 1:0.038 s, tol=3.08e-02 55.00 % -> 6.77 s 222 : 11.100 µm, [6.45e-06; 9.92e-01], 1:0.037 s, tol=3.00e-02 55.25 % -> 6.71 s 223 : 11.150 µm, [6.36e-06; 9.93e-01], 1:0.037 s, tol=2.96e-02 55.50 % -> 6.60 s 224 : 11.200 µm, [6.21e-06; 9.93e-01], 1:0.038 s, tol=2.93e-02 55.75 % -> 6.74 s 225 : 11.250 µm, [6.16e-06; 9.93e-01], 1:0.037 s, tol=3.12e-02 56.00 % -> 6.60 s 226 : 11.300 µm, [6.13e-06; 9.93e-01], 1:0.038 s, tol=3.48e-02 56.25 % -> 6.58 s 227 : 11.350 µm, [6.03e-06; 9.94e-01], 1:0.038 s, tol=3.87e-02 56.50 % -> 6.66 s 228 : 11.400 µm, [5.87e-06; 9.94e-01], 1:0.038 s, tol=4.11e-02 56.75 % -> 6.51 s 229 : 11.450 µm, [5.71e-06; 9.94e-01], 1:0.037 s, tol=4.15e-02 57.00 % -> 6.42 s 230 : 11.500 µm, [5.65e-06; 9.95e-01], 1:0.037 s, tol=4.01e-02 57.25 % -> 6.41 s 231 : 11.550 µm, [5.70e-06; 9.95e-01], 1:0.038 s, tol=3.84e-02 57.50 % -> 6.41 s 232 : 11.600 µm, [5.74e-06; 9.95e-01], 1:0.038 s, tol=3.71e-02 57.75 % -> 6.35 s 233 : 11.650 µm, [5.66e-06; 9.95e-01], 1:0.038 s, tol=3.53e-02 58.00 % -> 6.36 s 234 : 11.700 µm, [5.61e-06; 9.95e-01], 1:0.038 s, tol=3.48e-02 58.25 % -> 6.36 s 235 : 11.750 µm, [5.60e-06; 9.95e-01], 1:0.038 s, tol=3.73e-02 58.50 % -> 6.30 s 236 : 11.800 µm, [5.51e-06; 9.96e-01], 1:0.039 s, tol=4.00e-02 58.75 % -> 6.41 s 237 : 11.850 µm, [5.41e-06; 9.96e-01], 1:0.039 s, tol=4.25e-02 59.00 % -> 6.39 s 238 : 11.900 µm, [5.38e-06; 9.96e-01], 1:0.039 s, tol=4.26e-02 59.25 % -> 6.29 s 239 : 11.950 µm, [5.51e-06; 9.96e-01], 1:0.039 s, tol=4.20e-02 59.50 % -> 6.29 s 240 : 12.000 µm, [5.60e-06; 9.96e-01], 1:0.038 s, tol=4.03e-02 59.75 % -> 6.17 s 241 : 12.050 µm, [5.53e-06; 9.96e-01], 1:0.039 s, tol=3.75e-02 60.00 % -> 6.20 s 242 : 12.100 µm, [5.56e-06; 9.96e-01], 1:0.038 s, tol=3.75e-02 60.25 % -> 6.11 s 243 : 12.150 µm, [5.52e-06; 9.97e-01], 1:0.038 s, tol=3.99e-02 60.50 % -> 5.98 s 244 : 12.200 µm, [5.47e-06; 9.97e-01], 1:0.037 s, tol=4.29e-02 60.75 % -> 5.88 s 245 : 12.250 µm, [5.50e-06; 9.97e-01], 1:0.038 s, tol=4.42e-02 61.00 % -> 5.99 s 246 : 12.300 µm, [5.73e-06; 9.97e-01], 1:0.038 s, tol=4.43e-02 61.25 % -> 5.89 s 247 : 12.350 µm, [5.83e-06; 9.97e-01], 1:0.039 s, tol=4.23e-02 61.50 % -> 5.99 s 248 : 12.400 µm, [5.81e-06; 9.97e-01], 1:0.038 s, tol=3.95e-02 61.75 % -> 5.81 s 249 : 12.450 µm, [5.86e-06; 9.97e-01], 1:0.038 s, tol=4.02e-02 62.00 % -> 5.71 s 250 : 12.500 µm, [5.85e-06; 9.97e-01], 1:0.038 s, tol=4.32e-02 62.25 % -> 5.76 s 251 : 12.550 µm, [5.91e-06; 9.97e-01], 1:0.038 s, tol=4.52e-02 62.50 % -> 5.70 s 252 : 12.600 µm, [6.26e-06; 9.97e-01], 1:0.038 s, tol=4.57e-02 62.75 % -> 5.67 s 253 : 12.650 µm, [6.37e-06; 9.97e-01], 1:0.038 s, tol=4.31e-02 63.00 % -> 5.62 s 254 : 12.700 µm, [6.44e-06; 9.98e-01], 1:0.038 s, tol=4.07e-02 63.25 % -> 5.62 s 255 : 12.750 µm, [6.46e-06; 9.98e-01], 1:0.038 s, tol=4.15e-02 63.50 % -> 5.54 s 256 : 12.800 µm, [6.47e-06; 9.98e-01], 1:0.037 s, tol=4.43e-02 63.75 % -> 5.39 s 257 : 12.850 µm, [6.91e-06; 9.98e-01], 1:0.038 s, tol=4.65e-02 64.00 % -> 5.41 s 258 : 12.900 µm, [7.01e-06; 9.98e-01], 1:0.037 s, tol=4.45e-02 64.25 % -> 5.34 s 259 : 12.950 µm, [6.99e-06; 9.98e-01], 1:0.037 s, tol=4.13e-02 64.50 % -> 5.27 s 260 : 13.000 µm, [6.98e-06; 9.98e-01], 1:0.037 s, tol=4.14e-02 64.75 % -> 5.23 s 261 : 13.050 µm, [6.91e-06; 9.98e-01], 1:0.037 s, tol=4.41e-02 65.00 % -> 5.22 s 262 : 13.100 µm, [7.11e-06; 9.98e-01], 1:0.037 s, tol=4.64e-02 65.25 % -> 5.16 s 263 : 13.150 µm, [7.16e-06; 9.98e-01], 1:0.038 s, tol=4.35e-02 65.50 % -> 5.18 s 264 : 13.200 µm, [7.20e-06; 9.98e-01], 1:0.037 s, tol=4.08e-02 65.75 % -> 5.12 s 265 : 13.250 µm, [7.21e-06; 9.98e-01], 1:0.038 s, tol=4.22e-02 66.00 % -> 5.19 s 266 : 13.300 µm, [7.25e-06; 9.98e-01], 1:0.038 s, tol=4.48e-02 66.25 % -> 5.07 s 267 : 13.350 µm, [7.50e-06; 9.98e-01], 1:0.038 s, tol=4.48e-02 66.50 % -> 5.05 s 268 : 13.400 µm, [7.52e-06; 9.98e-01], 1:0.038 s, tol=4.09e-02 66.75 % -> 5.06 s 269 : 13.450 µm, [7.59e-06; 9.98e-01], 1:0.040 s, tol=4.04e-02 67.00 % -> 5.25 s 270 : 13.500 µm, [7.57e-06; 9.98e-01], 1:0.039 s, tol=4.28e-02 67.25 % -> 5.13 s 271 : 13.550 µm, [7.89e-06; 9.98e-01], 1:0.040 s, tol=4.45e-02 67.50 % -> 5.16 s 272 : 13.600 µm, [7.94e-06; 9.98e-01], 1:0.038 s, tol=4.04e-02 67.75 % -> 4.94 s 273 : 13.650 µm, [8.07e-06; 9.99e-01], 1:0.038 s, tol=3.91e-02 68.00 % -> 4.92 s 274 : 13.700 µm, [8.09e-06; 9.99e-01], 1:0.037 s, tol=4.15e-02 68.25 % -> 4.74 s 275 : 13.750 µm, [8.47e-06; 9.99e-01], 1:0.037 s, tol=4.34e-02 68.50 % -> 4.68 s 276 : 13.800 µm, [8.62e-06; 9.99e-01], 1:0.037 s, tol=3.95e-02 68.75 % -> 4.67 s 277 : 13.850 µm, [8.76e-06; 9.99e-01], 1:0.038 s, tol=3.84e-02 69.00 % -> 4.66 s 278 : 13.900 µm, [8.91e-06; 9.99e-01], 1:0.038 s, tol=4.11e-02 69.25 % -> 4.62 s 279 : 13.950 µm, [9.37e-06; 9.99e-01], 1:0.040 s, tol=4.18e-02 69.50 % -> 4.88 s 280 : 14.000 µm, [9.56e-06; 9.99e-01], 1:0.039 s, tol=3.76e-02 69.75 % -> 4.74 s 281 : 14.050 µm, [9.70e-06; 9.99e-01], 1:0.039 s, tol=3.77e-02 70.00 % -> 4.71 s 282 : 14.100 µm, [1.02e-05; 9.99e-01], 1:0.040 s, tol=4.09e-02 70.25 % -> 4.80 s 283 : 14.150 µm, [1.05e-05; 9.99e-01], 1:0.040 s, tol=3.84e-02 70.50 % -> 4.71 s 284 : 14.200 µm, [1.08e-05; 9.99e-01], 1:0.038 s, tol=3.55e-02 70.75 % -> 4.46 s 285 : 14.250 µm, [1.10e-05; 9.99e-01], 1:0.038 s, tol=3.73e-02 71.00 % -> 4.44 s 286 : 14.300 µm, [1.18e-05; 9.99e-01], 1:0.038 s, tol=3.94e-02 71.25 % -> 4.41 s 287 : 14.350 µm, [1.21e-05; 9.99e-01], 1:0.038 s, tol=3.55e-02 71.50 % -> 4.37 s 288 : 14.400 µm, [1.24e-05; 9.99e-01], 1:0.039 s, tol=3.51e-02 71.75 % -> 4.36 s 289 : 14.450 µm, [1.32e-05; 9.99e-01], 1:0.039 s, tol=3.82e-02 72.00 % -> 4.32 s 290 : 14.500 µm, [1.38e-05; 9.99e-01], 1:0.038 s, tol=3.54e-02 72.25 % -> 4.18 s 291 : 14.550 µm, [1.41e-05; 9.99e-01], 1:0.037 s, tol=3.32e-02 72.50 % -> 4.08 s 292 : 14.600 µm, [1.49e-05; 9.99e-01], 1:0.037 s, tol=3.57e-02 72.75 % -> 4.04 s 293 : 14.650 µm, [1.57e-05; 9.99e-01], 1:0.037 s, tol=3.55e-02 73.00 % -> 4.00 s 294 : 14.700 µm, [1.62e-05; 9.99e-01], 1:0.037 s, tol=3.20e-02 73.25 % -> 4.01 s 295 : 14.750 µm, [1.68e-05; 9.99e-01], 1:0.038 s, tol=3.34e-02 73.50 % -> 4.02 s 296 : 14.800 µm, [1.78e-05; 9.99e-01], 1:0.039 s, tol=3.49e-02 73.75 % -> 4.12 s 297 : 14.850 µm, [1.83e-05; 9.99e-01], 1:0.038 s, tol=3.16e-02 74.00 % -> 3.96 s 298 : 14.900 µm, [1.89e-05; 9.99e-01], 1:0.038 s, tol=3.20e-02 74.25 % -> 3.89 s 299 : 14.950 µm, [1.95e-05; 9.99e-01], 1:0.038 s, tol=3.36e-02 74.50 % -> 3.91 s 300 : 15.000 µm, [2.00e-05; 9.99e-01], 1:0.038 s, tol=3.07e-02 74.75 % -> 3.80 s 301 : 15.050 µm, [2.04e-05; 9.99e-01], 1:0.038 s, tol=3.08e-02 75.00 % -> 3.85 s 302 : 15.100 µm, [2.08e-05; 9.99e-01], 1:0.038 s, tol=3.21e-02 75.25 % -> 3.78 s 303 : 15.150 µm, [2.11e-05; 9.99e-01], 1:0.037 s, tol=2.95e-02 75.50 % -> 3.64 s 304 : 15.200 µm, [2.14e-05; 9.99e-01], 1:0.037 s, tol=2.96e-02 75.75 % -> 3.61 s 305 : 15.250 µm, [2.18e-05; 9.99e-01], 1:0.037 s, tol=3.05e-02 76.00 % -> 3.60 s 306 : 15.300 µm, [2.20e-05; 9.99e-01], 1:0.037 s, tol=2.84e-02 76.25 % -> 3.52 s 307 : 15.350 µm, [2.23e-05; 9.99e-01], 1:0.038 s, tol=2.87e-02 76.50 % -> 3.55 s 308 : 15.400 µm, [2.26e-05; 9.99e-01], 1:0.037 s, tol=2.87e-02 76.75 % -> 3.45 s 309 : 15.450 µm, [2.28e-05; 9.99e-01], 1:0.037 s, tol=2.73e-02 77.00 % -> 3.43 s 310 : 15.500 µm, [2.30e-05; 9.99e-01], 1:0.037 s, tol=2.78e-02 77.25 % -> 3.36 s 311 : 15.550 µm, [2.32e-05; 9.99e-01], 1:0.037 s, tol=2.69e-02 77.50 % -> 3.34 s 312 : 15.600 µm, [2.33e-05; 9.99e-01], 1:0.038 s, tol=2.61e-02 77.75 % -> 3.36 s 313 : 15.650 µm, [2.34e-05; 9.99e-01], 1:0.037 s, tol=2.67e-02 78.00 % -> 3.29 s 314 : 15.700 µm, [2.35e-05; 9.99e-01], 1:0.038 s, tol=2.52e-02 78.25 % -> 3.27 s 315 : 15.750 µm, [2.36e-05; 9.99e-01], 1:0.038 s, tol=2.48e-02 78.50 % -> 3.31 s 316 : 15.800 µm, [2.37e-05; 9.99e-01], 1:0.037 s, tol=2.54e-02 78.75 % -> 3.18 s 317 : 15.850 µm, [2.37e-05; 9.99e-01], 1:0.037 s, tol=2.39e-02 79.00 % -> 3.15 s 318 : 15.900 µm, [2.38e-05; 9.99e-01], 1:0.037 s, tol=2.32e-02 79.25 % -> 3.11 s 319 : 15.950 µm, [2.39e-05; 9.99e-01], 1:0.039 s, tol=2.41e-02 79.50 % -> 3.22 s 320 : 16.000 µm, [2.39e-05; 9.99e-01], 1:0.040 s, tol=2.23e-02 79.75 % -> 3.22 s 321 : 16.050 µm, [2.39e-05; 9.99e-01], 1:0.039 s, tol=2.20e-02 80.00 % -> 3.09 s 322 : 16.100 µm, [2.40e-05; 9.99e-01], 1:0.037 s, tol=2.28e-02 80.25 % -> 2.95 s 323 : 16.150 µm, [2.40e-05; 9.99e-01], 1:0.037 s, tol=2.10e-02 80.50 % -> 2.91 s 324 : 16.200 µm, [2.40e-05; 9.99e-01], 1:0.038 s, tol=2.09e-02 80.75 % -> 2.90 s 325 : 16.250 µm, [2.40e-05; 9.99e-01], 1:0.039 s, tol=2.12e-02 81.00 % -> 2.94 s 326 : 16.300 µm, [2.40e-05; 9.99e-01], 1:0.038 s, tol=1.93e-02 81.25 % -> 2.89 s 327 : 16.350 µm, [2.40e-05; 9.99e-01], 1:0.039 s, tol=1.93e-02 81.50 % -> 2.89 s 328 : 16.400 µm, [2.41e-05; 9.99e-01], 1:0.040 s, tol=1.98e-02 81.75 % -> 2.94 s 329 : 16.450 µm, [2.41e-05; 9.99e-01], 1:0.038 s, tol=1.81e-02 82.00 % -> 2.72 s 330 : 16.500 µm, [2.41e-05; 9.99e-01], 1:0.038 s, tol=1.66e-02 82.25 % -> 2.70 s 331 : 16.550 µm, [2.41e-05; 9.99e-01], 1:0.038 s, tol=1.73e-02 82.50 % -> 2.69 s 332 : 16.600 µm, [2.41e-05; 9.99e-01], 1:0.038 s, tol=1.72e-02 82.75 % -> 2.63 s 333 : 16.650 µm, [2.41e-05; 9.99e-01], 1:0.038 s, tol=1.66e-02 83.00 % -> 2.58 s 334 : 16.700 µm, [2.41e-05; 9.99e-01], 1:0.037 s, tol=1.63e-02 83.25 % -> 2.51 s 335 : 16.750 µm, [2.41e-05; 9.99e-01], 1:0.037 s, tol=1.71e-02 83.50 % -> 2.43 s 336 : 16.800 µm, [2.41e-05; 9.99e-01], 1:0.037 s, tol=1.88e-02 83.75 % -> 2.40 s 337 : 16.850 µm, [2.41e-05; 9.99e-01], 1:0.037 s, tol=2.21e-02 84.00 % -> 2.36 s 338 : 16.900 µm, [2.42e-05; 9.99e-01], 1:0.038 s, tol=2.61e-02 84.25 % -> 2.37 s 339 : 16.950 µm, [2.42e-05; 9.99e-01], 1:0.039 s, tol=3.05e-02 84.50 % -> 2.40 s 340 : 17.000 µm, [2.42e-05; 1.00e+00], 1:0.039 s, tol=3.18e-02 84.75 % -> 2.40 s 341 : 17.050 µm, [2.42e-05; 1.00e+00], 1:0.040 s, tol=3.25e-02 85.00 % -> 2.39 s 342 : 17.100 µm, [2.42e-05; 1.00e+00], 1:0.038 s, tol=3.29e-02 85.25 % -> 2.27 s 343 : 17.150 µm, [2.43e-05; 1.00e+00], 1:0.039 s, tol=2.48e-02 85.50 % -> 2.25 s 344 : 17.200 µm, [2.43e-05; 1.01e+00], 1:0.039 s, tol=1.33e-02 85.75 % -> 2.25 s 345 : 17.250 µm, [2.43e-05; 1.01e+00], 1:0.039 s, tol=8.31e-03 86.00 % -> 2.20 s 346 : 17.300 µm, [2.43e-05; 1.01e+00], 1:0.040 s, tol=6.99e-03 86.25 % -> 2.20 s 347 : 17.350 µm, [2.43e-05; 1.01e+00], 1:0.038 s, tol=6.52e-03 86.50 % -> 2.06 s 348 : 17.400 µm, [2.43e-05; 1.01e+00], 1:0.039 s, tol=6.25e-03 86.75 % -> 2.06 s 349 : 17.450 µm, [2.43e-05; 1.01e+00], 1:0.038 s, tol=6.08e-03 87.00 % -> 2.00 s 350 : 17.500 µm, [2.43e-05; 1.01e+00], 1:0.038 s, tol=5.96e-03 87.25 % -> 1.92 s 351 : 17.550 µm, [2.43e-05; 1.01e+00], 1:0.038 s, tol=5.88e-03 87.50 % -> 1.88 s 352 : 17.600 µm, [2.43e-05; 1.01e+00], 1:0.037 s, tol=5.83e-03 87.75 % -> 1.80 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.71 s 356 : 17.800 µm, [2.43e-05; 1.01e+00], 1:0.037 s, tol=5.71e-03 88.75 % -> 1.67 s 357 : 17.850 µm, [2.43e-05; 1.01e+00], 1:0.037 s, tol=5.69e-03 89.00 % -> 1.64 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.53 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.43 s 363 : 18.150 µm, [2.43e-05; 1.01e+00], 1:0.037 s, tol=5.59e-03 90.50 % -> 1.40 s 364 : 18.200 µm, [2.43e-05; 1.01e+00], 1:0.037 s, tol=5.58e-03 90.75 % -> 1.36 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.038 s, tol=5.55e-03 91.25 % -> 1.31 s 367 : 18.350 µm, [2.43e-05; 1.01e+00], 1:0.038 s, tol=5.54e-03 91.50 % -> 1.29 s 368 : 18.400 µm, [2.43e-05; 1.01e+00], 1:0.037 s, tol=5.52e-03 91.75 % -> 1.23 s Saved: /home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/v1.7.3/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.3/examples/PhaseField/results/Shear2D/Test/ElasIsot_Miehe_AT1_DP_optimMesh/simulation.pickle Saved: /home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/v1.7.3/examples/PhaseField/results/Shear2D/Test/ElasIsot_Miehe_AT1_DP_optimMesh/summary.txt Generate movie 01/21 (4.76 %) 2.24 s Generate movie 02/21 (9.52 %) 1.74 s Generate movie 03/21 (14.29 %) 1.66 s Generate movie 04/21 (19.05 %) 1.66 s Generate movie 05/21 (23.81 %) 1.48 s Generate movie 06/21 (28.57 %) 1.49 s Generate movie 07/21 (33.33 %) 1.37 s Generate movie 08/21 (38.10 %) 1.20 s Generate movie 09/21 (42.86 %) 1.11 s Generate movie 10/21 (47.62 %) 967.87 ms Generate movie 11/21 (52.38 %) 907.62 ms Generate movie 12/21 (57.14 %) 803.68 ms Generate movie 13/21 (61.90 %) 721.30 ms Generate movie 14/21 (66.67 %) 646.41 ms Generate movie 15/21 (71.43 %) 541.47 ms Generate movie 16/21 (76.19 %) 439.06 ms Generate movie 17/21 (80.95 %) 352.58 ms Generate movie 18/21 (85.71 %) 281.08 ms Generate movie 19/21 (90.48 %) 182.70 ms Generate movie 20/21 (95.24 %) 93.85 ms Generate movie 21/21 (100.00 %) 0.00 µs Mesh : 45.696 ms Boundary Conditions : 4.744 ms Matrix : 5.346 s Split : 1.467 s Solver : 9.036 s PostProcessing : 100.527 ms Resolution phase field : 13.708 s Display : 435.456 ms PyVista_Interface : 3.067 s | .. code-block:: Python :lineno-start: 15 import matplotlib.pyplot as plt import numpy as np from EasyFEA import ( Display, Folder, Models, 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 folder = Folder.Results_Dir() + f"{dim}D" 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 = Simulations.PhaseField.Folder( folder, "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.Elastic.Isotropic( 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.PhaseField(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, _, Ku, 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(Ku[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.PhaseField = 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() # Display.plt.close("all") Tic.Clear() return folder_save if __name__ == "__main__": # generates configs Splits = [] Regus = [] for split in splits.copy(): for regu in regus.copy(): Splits.append(split) Regus.append(regu) list_folder: list[str] = [] 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): list_folder.append(result) else: for split, regu in zip(Splits, Regus): list_folder.append(DoSimu(split, regu)) .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 18.470 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 `_