.. 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-289 .. 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.8.1/examples/PhaseField/results/Shear2D/Test/Miehe_AT1_optimMesh 1 : 0.050 µm, [0.00e+00; 0.00e+00], 1:0.044s, tol=1.00e+00 2 : 0.100 µm, [0.00e+00; 0.00e+00], 1:0.037s, tol=1.00e+00 0.25% -> 14.65s 3 : 0.150 µm, [0.00e+00; 0.00e+00], 1:0.037s, tol=5.56e-01 0.50% -> 14.60s 4 : 0.200 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=4.38e-01 0.75% -> 14.42s 5 : 0.250 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.60e-01 1.00% -> 14.28s 6 : 0.300 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.06e-01 1.25% -> 14.22s 7 : 0.350 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.65e-01 1.50% -> 14.14s 8 : 0.400 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.34e-01 1.75% -> 14.10s 9 : 0.450 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.10e-01 2.00% -> 14.02s 10 : 0.500 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.90e-01 2.25% -> 14.02s 11 : 0.550 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.74e-01 2.50% -> 13.94s 12 : 0.600 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.60e-01 2.75% -> 13.91s 13 : 0.650 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.48e-01 3.00% -> 13.84s 14 : 0.700 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.38e-01 3.25% -> 13.87s 15 : 0.750 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.29e-01 3.50% -> 13.85s 16 : 0.800 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.21e-01 3.75% -> 13.79s 17 : 0.850 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.14e-01 4.00% -> 13.88s 18 : 0.900 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.08e-01 4.25% -> 13.76s 19 : 0.950 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.02e-01 4.50% -> 13.66s 20 : 1.000 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=9.75e-02 4.75% -> 13.68s 21 : 1.050 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=9.30e-02 5.00% -> 13.59s 22 : 1.100 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=8.88e-02 5.25% -> 13.65s 23 : 1.150 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=8.51e-02 5.50% -> 13.55s 24 : 1.200 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=8.16e-02 5.75% -> 13.65s 25 : 1.250 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=7.84e-02 6.00% -> 13.71s 26 : 1.300 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=7.54e-02 6.25% -> 13.60s 27 : 1.350 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=7.27e-02 6.50% -> 13.49s 28 : 1.400 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=7.02e-02 6.75% -> 13.48s 29 : 1.450 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=6.78e-02 7.00% -> 13.47s 30 : 1.500 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=6.56e-02 7.25% -> 13.40s 31 : 1.550 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=6.35e-02 7.50% -> 13.31s 32 : 1.600 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=6.15e-02 7.75% -> 13.29s 33 : 1.650 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=5.97e-02 8.00% -> 13.21s 34 : 1.700 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=5.80e-02 8.25% -> 13.21s 35 : 1.750 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=5.63e-02 8.50% -> 13.11s 36 : 1.800 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=5.48e-02 8.75% -> 13.11s 37 : 1.850 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=5.33e-02 9.00% -> 13.03s 38 : 1.900 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=5.19e-02 9.25% -> 13.00s 39 : 1.950 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=5.06e-02 9.50% -> 12.95s 40 : 2.000 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=4.94e-02 9.75% -> 12.96s 41 : 2.050 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=4.82e-02 10.00% -> 12.86s 42 : 2.100 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=4.71e-02 10.25% -> 12.83s 43 : 2.150 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=4.60e-02 10.50% -> 12.80s 44 : 2.200 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=4.49e-02 10.75% -> 12.77s 45 : 2.250 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=4.40e-02 11.00% -> 12.68s 46 : 2.300 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=4.30e-02 11.25% -> 12.75s 47 : 2.350 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=4.21e-02 11.50% -> 12.75s 48 : 2.400 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=4.12e-02 11.75% -> 12.63s 49 : 2.450 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=4.04e-02 12.00% -> 12.55s 50 : 2.500 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.96e-02 12.25% -> 12.63s 51 : 2.550 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.88e-02 12.50% -> 12.60s 52 : 2.600 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.81e-02 12.75% -> 12.69s 53 : 2.650 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.74e-02 13.00% -> 12.56s 54 : 2.700 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.67e-02 13.25% -> 12.57s 55 : 2.750 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.60e-02 13.50% -> 12.47s 56 : 2.800 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.54e-02 13.75% -> 12.50s 57 : 2.850 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.48e-02 14.00% -> 12.44s 58 : 2.900 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.42e-02 14.25% -> 12.41s 59 : 2.950 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.36e-02 14.50% -> 12.27s 60 : 3.000 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.31e-02 14.75% -> 12.25s 61 : 3.050 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.25e-02 15.00% -> 12.18s 62 : 3.100 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.20e-02 15.25% -> 12.16s 63 : 3.150 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.15e-02 15.50% -> 12.12s 64 : 3.200 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.10e-02 15.75% -> 12.13s 65 : 3.250 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.05e-02 16.00% -> 12.12s 66 : 3.300 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=3.01e-02 16.25% -> 12.06s 67 : 3.350 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.96e-02 16.50% -> 11.95s 68 : 3.400 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.92e-02 16.75% -> 11.95s 69 : 3.450 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.88e-02 17.00% -> 11.86s 70 : 3.500 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.84e-02 17.25% -> 11.87s 71 : 3.550 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.80e-02 17.50% -> 11.83s 72 : 3.600 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.76e-02 17.75% -> 11.84s 73 : 3.650 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.72e-02 18.00% -> 11.76s 74 : 3.700 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.68e-02 18.25% -> 11.78s 75 : 3.750 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.65e-02 18.50% -> 11.71s 76 : 3.800 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.61e-02 18.75% -> 11.70s 77 : 3.850 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.58e-02 19.00% -> 11.62s 78 : 3.900 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.55e-02 19.25% -> 11.63s 79 : 3.950 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.52e-02 19.50% -> 11.65s 80 : 4.000 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.48e-02 19.75% -> 11.68s 81 : 4.050 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.45e-02 20.00% -> 11.55s 82 : 4.100 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.42e-02 20.25% -> 11.58s 83 : 4.150 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.40e-02 20.50% -> 11.47s 84 : 4.200 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.37e-02 20.75% -> 11.53s 85 : 4.250 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.34e-02 21.00% -> 11.38s 86 : 4.300 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.31e-02 21.25% -> 11.34s 87 : 4.350 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.29e-02 21.50% -> 11.43s 88 : 4.400 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.26e-02 21.75% -> 11.31s 89 : 4.450 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.23e-02 22.00% -> 11.18s 90 : 4.500 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.21e-02 22.25% -> 11.18s 91 : 4.550 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.19e-02 22.50% -> 11.11s 92 : 4.600 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.16e-02 22.75% -> 11.10s 93 : 4.650 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.14e-02 23.00% -> 11.03s 94 : 4.700 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.12e-02 23.25% -> 11.04s 95 : 4.750 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.09e-02 23.50% -> 10.99s 96 : 4.800 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.07e-02 23.75% -> 10.98s 97 : 4.850 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.05e-02 24.00% -> 10.91s 98 : 4.900 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.03e-02 24.25% -> 10.89s 99 : 4.950 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=2.01e-02 24.50% -> 10.83s 100 : 5.000 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.99e-02 24.75% -> 10.78s 101 : 5.050 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.97e-02 25.00% -> 10.73s 102 : 5.100 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.95e-02 25.25% -> 10.72s 103 : 5.150 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.93e-02 25.50% -> 10.68s 104 : 5.200 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.91e-02 25.75% -> 10.74s 105 : 5.250 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.90e-02 26.00% -> 10.66s 106 : 5.300 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.88e-02 26.25% -> 10.69s 107 : 5.350 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.86e-02 26.50% -> 10.65s 108 : 5.400 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.84e-02 26.75% -> 10.61s 109 : 5.450 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.83e-02 27.00% -> 10.52s 110 : 5.500 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.81e-02 27.25% -> 10.50s 111 : 5.550 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.79e-02 27.50% -> 10.50s 112 : 5.600 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.78e-02 27.75% -> 10.49s 113 : 5.650 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.76e-02 28.00% -> 10.34s 114 : 5.700 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.75e-02 28.25% -> 10.37s 115 : 5.750 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.73e-02 28.50% -> 10.27s 116 : 5.800 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.72e-02 28.75% -> 10.26s 117 : 5.850 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.70e-02 29.00% -> 10.35s 118 : 5.900 µm, [0.00e+00; 0.00e+00], 1:0.037s, tol=1.69e-02 29.25% -> 10.55s 119 : 5.950 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.67e-02 29.50% -> 10.19s 120 : 6.000 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.66e-02 29.75% -> 10.13s 121 : 6.050 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.65e-02 30.00% -> 10.04s 122 : 6.100 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.63e-02 30.25% -> 10.01s 123 : 6.150 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.62e-02 30.50% -> 9.95s 124 : 6.200 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.61e-02 30.75% -> 9.93s 125 : 6.250 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.59e-02 31.00% -> 9.86s 126 : 6.300 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.58e-02 31.25% -> 9.86s 127 : 6.350 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.57e-02 31.50% -> 9.78s 128 : 6.400 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.56e-02 31.75% -> 9.78s 129 : 6.450 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.54e-02 32.00% -> 9.72s 130 : 6.500 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.53e-02 32.25% -> 9.75s 131 : 6.550 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.52e-02 32.50% -> 9.72s 132 : 6.600 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.51e-02 32.75% -> 9.68s 133 : 6.650 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.50e-02 33.00% -> 9.66s 134 : 6.700 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.49e-02 33.25% -> 9.73s 135 : 6.750 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.48e-02 33.50% -> 9.61s 136 : 6.800 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.47e-02 33.75% -> 9.60s 137 : 6.850 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.45e-02 34.00% -> 9.52s 138 : 6.900 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.44e-02 34.25% -> 9.52s 139 : 6.950 µm, [0.00e+00; 0.00e+00], 1:0.036s, tol=1.43e-02 34.50% -> 9.45s 140 : 7.000 µm, [6.00e-07; 1.59e-03], 1:0.036s, tol=1.42e-02 34.75% -> 9.45s 141 : 7.050 µm, [2.22e-06; 6.20e-03], 1:0.036s, tol=1.39e-02 35.00% -> 9.46s 142 : 7.100 µm, [4.16e-06; 1.23e-02], 1:0.036s, tol=1.38e-02 35.25% -> 9.42s 143 : 7.150 µm, [6.13e-06; 1.92e-02], 1:0.036s, tol=1.37e-02 35.50% -> 9.33s 144 : 7.200 µm, [7.97e-06; 2.65e-02], 1:0.036s, tol=1.35e-02 35.75% -> 9.34s 145 : 7.250 µm, [9.65e-06; 3.41e-02], 1:0.036s, tol=1.35e-02 36.00% -> 9.26s 146 : 7.300 µm, [1.12e-05; 4.18e-02], 1:0.036s, tol=1.34e-02 36.25% -> 9.27s 147 : 7.350 µm, [1.25e-05; 4.96e-02], 1:0.036s, tol=1.33e-02 36.50% -> 9.17s 148 : 7.400 µm, [1.36e-05; 5.75e-02], 1:0.036s, tol=1.32e-02 36.75% -> 9.19s 149 : 7.450 µm, [1.46e-05; 6.55e-02], 1:0.036s, tol=1.31e-02 37.00% -> 9.11s 150 : 7.500 µm, [1.55e-05; 7.35e-02], 1:0.036s, tol=1.30e-02 37.25% -> 9.08s 151 : 7.550 µm, [1.62e-05; 8.16e-02], 1:0.036s, tol=1.29e-02 37.50% -> 9.03s 152 : 7.600 µm, [1.68e-05; 8.98e-02], 1:0.036s, tol=1.28e-02 37.75% -> 9.01s 153 : 7.650 µm, [1.72e-05; 9.80e-02], 1:0.036s, tol=1.28e-02 38.00% -> 8.91s 154 : 7.700 µm, [1.76e-05; 1.06e-01], 1:0.036s, tol=1.27e-02 38.25% -> 8.88s 155 : 7.750 µm, [1.79e-05; 1.15e-01], 1:0.036s, tol=1.26e-02 38.50% -> 8.90s 156 : 7.800 µm, [1.81e-05; 1.23e-01], 1:0.036s, tol=1.25e-02 38.75% -> 8.91s 157 : 7.850 µm, [1.82e-05; 1.32e-01], 1:0.036s, tol=1.24e-02 39.00% -> 8.85s 158 : 7.900 µm, [1.83e-05; 1.40e-01], 1:0.036s, tol=1.24e-02 39.25% -> 8.84s 159 : 7.950 µm, [1.83e-05; 1.49e-01], 1:0.037s, tol=1.23e-02 39.50% -> 8.86s 160 : 8.000 µm, [1.82e-05; 1.58e-01], 1:0.037s, tol=1.22e-02 39.75% -> 8.82s 161 : 8.050 µm, [1.81e-05; 1.67e-01], 1:0.037s, tol=1.21e-02 40.00% -> 8.78s 162 : 8.100 µm, [1.80e-05; 1.76e-01], 1:0.037s, tol=1.21e-02 40.25% -> 8.75s 163 : 8.150 µm, [1.78e-05; 1.85e-01], 1:0.037s, tol=1.20e-02 40.50% -> 8.71s 164 : 8.200 µm, [1.75e-05; 1.94e-01], 1:0.037s, tol=1.19e-02 40.75% -> 8.70s 165 : 8.250 µm, [1.73e-05; 2.03e-01], 1:0.037s, tol=1.18e-02 41.00% -> 8.66s 166 : 8.300 µm, [1.70e-05; 2.12e-01], 1:0.036s, tol=1.18e-02 41.25% -> 8.56s 167 : 8.350 µm, [1.67e-05; 2.22e-01], 1:0.037s, tol=1.17e-02 41.50% -> 8.56s 168 : 8.400 µm, [1.64e-05; 2.31e-01], 1:0.037s, tol=1.16e-02 41.75% -> 8.55s 169 : 8.450 µm, [1.60e-05; 2.41e-01], 1:0.038s, tol=1.16e-02 42.00% -> 8.71s 170 : 8.500 µm, [1.57e-05; 2.51e-01], 1:0.038s, tol=1.15e-02 42.25% -> 8.77s 171 : 8.550 µm, [1.53e-05; 2.60e-01], 1:0.037s, tol=1.14e-02 42.50% -> 8.57s 172 : 8.600 µm, [1.49e-05; 2.70e-01], 1:0.037s, tol=1.14e-02 42.75% -> 8.37s 173 : 8.650 µm, [1.45e-05; 2.80e-01], 1:0.037s, tol=1.13e-02 43.00% -> 8.49s 174 : 8.700 µm, [1.42e-05; 2.90e-01], 1:0.037s, tol=1.13e-02 43.25% -> 8.37s 175 : 8.750 µm, [1.38e-05; 3.01e-01], 1:0.036s, tol=1.12e-02 43.50% -> 8.19s 176 : 8.800 µm, [1.34e-05; 3.11e-01], 1:0.036s, tol=1.11e-02 43.75% -> 8.14s 177 : 8.850 µm, [1.30e-05; 3.21e-01], 1:0.036s, tol=1.11e-02 44.00% -> 8.06s 178 : 8.900 µm, [1.26e-05; 3.32e-01], 1:0.036s, tol=1.10e-02 44.25% -> 8.07s 179 : 8.950 µm, [1.22e-05; 3.43e-01], 1:0.036s, tol=1.10e-02 44.50% -> 8.00s 180 : 9.000 µm, [1.19e-05; 3.56e-01], 1:0.036s, tol=1.09e-02 44.75% -> 8.01s 181 : 9.050 µm, [1.15e-05; 3.71e-01], 1:0.036s, tol=1.08e-02 45.00% -> 7.92s 182 : 9.100 µm, [1.12e-05; 3.86e-01], 1:0.036s, tol=1.07e-02 45.25% -> 7.98s 183 : 9.150 µm, [1.09e-05; 4.01e-01], 1:0.036s, tol=1.07e-02 45.50% -> 7.84s 184 : 9.200 µm, [1.05e-05; 4.17e-01], 1:0.036s, tol=1.06e-02 45.75% -> 7.85s 185 : 9.250 µm, [1.02e-05; 4.32e-01], 1:0.036s, tol=1.06e-02 46.00% -> 7.83s 186 : 9.300 µm, [9.84e-06; 4.47e-01], 1:0.037s, tol=1.06e-02 46.25% -> 7.89s 187 : 9.350 µm, [9.49e-06; 4.62e-01], 1:0.037s, tol=1.05e-02 46.50% -> 7.91s 188 : 9.400 µm, [9.27e-06; 4.77e-01], 1:0.037s, tol=1.04e-02 46.75% -> 7.83s 189 : 9.450 µm, [9.10e-06; 4.94e-01], 1:0.037s, tol=1.03e-02 47.00% -> 7.80s 190 : 9.500 µm, [8.97e-06; 5.14e-01], 1:0.037s, tol=1.03e-02 47.25% -> 7.72s 191 : 9.550 µm, [8.88e-06; 5.39e-01], 1:0.037s, tol=1.02e-02 47.50% -> 7.67s 192 : 9.600 µm, [8.83e-06; 5.71e-01], 1:0.040s, tol=1.03e-02 47.75% -> 8.35s 193 : 9.650 µm, [8.78e-06; 6.08e-01], 1:0.037s, tol=1.04e-02 48.00% -> 7.65s 194 : 9.700 µm, [8.76e-06; 6.49e-01], 1:0.036s, tol=1.07e-02 48.25% -> 7.55s 195 : 9.750 µm, [8.75e-06; 6.95e-01], 1:0.036s, tol=1.13e-02 48.50% -> 7.48s 196 : 9.800 µm, [9.05e-06; 7.43e-01], 1:0.036s, tol=1.23e-02 48.75% -> 7.44s 197 : 9.850 µm, [9.39e-06; 7.94e-01], 1:0.036s, tol=1.47e-02 49.00% -> 7.34s 198 : 9.900 µm, [9.76e-06; 8.46e-01], 1:0.036s, tol=1.92e-02 49.25% -> 7.33s 199 : 9.950 µm, [9.84e-06; 8.94e-01], 1:0.036s, tol=2.57e-02 49.50% -> 7.31s 200 : 10.000 µm, [9.59e-06; 9.28e-01], 1:0.036s, tol=3.11e-02 49.75% -> 7.29s 201 : 10.050 µm, [9.14e-06; 9.49e-01], 1:0.036s, tol=3.15e-02 50.00% -> 7.20s 202 : 10.100 µm, [9.06e-06; 9.58e-01], 1:0.036s, tol=2.98e-02 50.25% -> 7.19s 203 : 10.150 µm, [9.00e-06; 9.60e-01], 1:0.036s, tol=3.01e-02 50.50% -> 7.13s 204 : 10.200 µm, [9.21e-06; 9.60e-01], 1:0.036s, tol=3.37e-02 50.75% -> 7.12s 205 : 10.250 µm, [9.27e-06; 9.60e-01], 1:0.036s, tol=4.01e-02 51.00% -> 7.07s 206 : 10.300 µm, [9.16e-06; 9.71e-01], 1:0.036s, tol=4.63e-02 51.25% -> 7.04s 207 : 10.350 µm, [8.87e-06; 9.82e-01], 1:0.036s, tol=4.81e-02 51.50% -> 7.01s 208 : 10.400 µm, [8.45e-06; 9.89e-01], 1:0.036s, tol=4.32e-02 51.75% -> 6.96s 209 : 10.450 µm, [8.11e-06; 9.92e-01], 1:0.036s, tol=3.51e-02 52.00% -> 6.92s 210 : 10.500 µm, [8.09e-06; 9.94e-01], 1:0.036s, tol=3.08e-02 52.25% -> 6.89s 211 : 10.550 µm, [8.04e-06; 9.95e-01], 1:0.036s, tol=2.91e-02 52.50% -> 6.85s 212 : 10.600 µm, [7.90e-06; 9.95e-01], 1:0.036s, tol=2.80e-02 52.75% -> 6.86s 213 : 10.650 µm, [7.80e-06; 9.95e-01], 1:0.036s, tol=2.75e-02 53.00% -> 6.80s 214 : 10.700 µm, [7.73e-06; 9.94e-01], 1:0.036s, tol=2.85e-02 53.25% -> 6.82s 215 : 10.750 µm, [7.62e-06; 9.94e-01], 1:0.037s, tol=3.07e-02 53.50% -> 6.81s 216 : 10.800 µm, [7.45e-06; 9.93e-01], 1:0.037s, tol=3.36e-02 53.75% -> 6.75s 217 : 10.850 µm, [7.23e-06; 9.93e-01], 1:0.036s, tol=3.59e-02 54.00% -> 6.70s 218 : 10.900 µm, [6.95e-06; 9.92e-01], 1:0.036s, tol=3.61e-02 54.25% -> 6.67s 219 : 10.950 µm, [6.71e-06; 9.92e-01], 1:0.036s, tol=3.49e-02 54.50% -> 6.61s 220 : 11.000 µm, [6.54e-06; 9.92e-01], 1:0.037s, tol=3.27e-02 54.75% -> 6.64s 221 : 11.050 µm, [6.47e-06; 9.92e-01], 1:0.036s, tol=3.08e-02 55.00% -> 6.54s 222 : 11.100 µm, [6.45e-06; 9.92e-01], 1:0.037s, tol=3.00e-02 55.25% -> 6.59s 223 : 11.150 µm, [6.36e-06; 9.93e-01], 1:0.041s, tol=2.96e-02 55.50% -> 7.23s 224 : 11.200 µm, [6.21e-06; 9.93e-01], 1:0.036s, tol=2.93e-02 55.75% -> 6.44s 225 : 11.250 µm, [6.16e-06; 9.93e-01], 1:0.036s, tol=3.12e-02 56.00% -> 6.36s 226 : 11.300 µm, [6.13e-06; 9.93e-01], 1:0.036s, tol=3.48e-02 56.25% -> 6.34s 227 : 11.350 µm, [6.03e-06; 9.94e-01], 1:0.036s, tol=3.87e-02 56.50% -> 6.29s 228 : 11.400 µm, [5.87e-06; 9.94e-01], 1:0.036s, tol=4.11e-02 56.75% -> 6.27s 229 : 11.450 µm, [5.71e-06; 9.94e-01], 1:0.036s, tol=4.15e-02 57.00% -> 6.21s 230 : 11.500 µm, [5.65e-06; 9.95e-01], 1:0.036s, tol=4.01e-02 57.25% -> 6.20s 231 : 11.550 µm, [5.70e-06; 9.95e-01], 1:0.036s, tol=3.84e-02 57.50% -> 6.12s 232 : 11.600 µm, [5.74e-06; 9.95e-01], 1:0.036s, tol=3.71e-02 57.75% -> 6.11s 233 : 11.650 µm, [5.66e-06; 9.95e-01], 1:0.036s, tol=3.53e-02 58.00% -> 6.06s 234 : 11.700 µm, [5.61e-06; 9.95e-01], 1:0.036s, tol=3.48e-02 58.25% -> 6.04s 235 : 11.750 µm, [5.60e-06; 9.95e-01], 1:0.036s, tol=3.73e-02 58.50% -> 6.05s 236 : 11.800 µm, [5.51e-06; 9.96e-01], 1:0.036s, tol=4.00e-02 58.75% -> 6.02s 237 : 11.850 µm, [5.41e-06; 9.96e-01], 1:0.036s, tol=4.25e-02 59.00% -> 5.96s 238 : 11.900 µm, [5.38e-06; 9.96e-01], 1:0.037s, tol=4.26e-02 59.25% -> 5.96s 239 : 11.950 µm, [5.51e-06; 9.96e-01], 1:0.036s, tol=4.20e-02 59.50% -> 5.90s 240 : 12.000 µm, [5.60e-06; 9.96e-01], 1:0.037s, tol=4.03e-02 59.75% -> 5.99s 241 : 12.050 µm, [5.53e-06; 9.96e-01], 1:0.037s, tol=3.75e-02 60.00% -> 5.96s 242 : 12.100 µm, [5.56e-06; 9.96e-01], 1:0.037s, tol=3.75e-02 60.25% -> 5.87s 243 : 12.150 µm, [5.52e-06; 9.97e-01], 1:0.036s, tol=3.99e-02 60.50% -> 5.75s 244 : 12.200 µm, [5.47e-06; 9.97e-01], 1:0.036s, tol=4.29e-02 60.75% -> 5.72s 245 : 12.250 µm, [5.50e-06; 9.97e-01], 1:0.037s, tol=4.42e-02 61.00% -> 5.70s 246 : 12.300 µm, [5.73e-06; 9.97e-01], 1:0.037s, tol=4.43e-02 61.25% -> 5.78s 247 : 12.350 µm, [5.83e-06; 9.97e-01], 1:0.037s, tol=4.23e-02 61.50% -> 5.72s 248 : 12.400 µm, [5.81e-06; 9.97e-01], 1:0.037s, tol=3.95e-02 61.75% -> 5.61s 249 : 12.450 µm, [5.86e-06; 9.97e-01], 1:0.037s, tol=4.02e-02 62.00% -> 5.60s 250 : 12.500 µm, [5.85e-06; 9.97e-01], 1:0.036s, tol=4.32e-02 62.25% -> 5.49s 251 : 12.550 µm, [5.91e-06; 9.97e-01], 1:0.036s, tol=4.52e-02 62.50% -> 5.43s 252 : 12.600 µm, [6.26e-06; 9.97e-01], 1:0.036s, tol=4.57e-02 62.75% -> 5.39s 253 : 12.650 µm, [6.37e-06; 9.97e-01], 1:0.036s, tol=4.31e-02 63.00% -> 5.34s 254 : 12.700 µm, [6.44e-06; 9.98e-01], 1:0.036s, tol=4.07e-02 63.25% -> 5.33s 255 : 12.750 µm, [6.46e-06; 9.98e-01], 1:0.036s, tol=4.15e-02 63.50% -> 5.28s 256 : 12.800 µm, [6.47e-06; 9.98e-01], 1:0.036s, tol=4.43e-02 63.75% -> 5.24s 257 : 12.850 µm, [6.91e-06; 9.98e-01], 1:0.036s, tol=4.65e-02 64.00% -> 5.19s 258 : 12.900 µm, [7.01e-06; 9.98e-01], 1:0.036s, tol=4.45e-02 64.25% -> 5.18s 259 : 12.950 µm, [6.99e-06; 9.98e-01], 1:0.036s, tol=4.13e-02 64.50% -> 5.12s 260 : 13.000 µm, [6.98e-06; 9.98e-01], 1:0.036s, tol=4.14e-02 64.75% -> 5.10s 261 : 13.050 µm, [6.91e-06; 9.98e-01], 1:0.036s, tol=4.41e-02 65.00% -> 5.05s 262 : 13.100 µm, [7.11e-06; 9.98e-01], 1:0.036s, tol=4.64e-02 65.25% -> 5.02s 263 : 13.150 µm, [7.16e-06; 9.98e-01], 1:0.038s, tol=4.35e-02 65.50% -> 5.19s 264 : 13.200 µm, [7.20e-06; 9.98e-01], 1:0.040s, tol=4.08e-02 65.75% -> 5.47s 265 : 13.250 µm, [7.21e-06; 9.98e-01], 1:0.036s, tol=4.22e-02 66.00% -> 4.94s 266 : 13.300 µm, [7.25e-06; 9.98e-01], 1:0.036s, tol=4.48e-02 66.25% -> 4.91s 267 : 13.350 µm, [7.50e-06; 9.98e-01], 1:0.037s, tol=4.48e-02 66.50% -> 4.94s 268 : 13.400 µm, [7.52e-06; 9.98e-01], 1:0.037s, tol=4.09e-02 66.75% -> 4.89s 269 : 13.450 µm, [7.59e-06; 9.98e-01], 1:0.036s, tol=4.04e-02 67.00% -> 4.82s 270 : 13.500 µm, [7.57e-06; 9.98e-01], 1:0.037s, tol=4.28e-02 67.25% -> 4.79s 271 : 13.550 µm, [7.89e-06; 9.98e-01], 1:0.036s, tol=4.45e-02 67.50% -> 4.72s 272 : 13.600 µm, [7.94e-06; 9.98e-01], 1:0.036s, tol=4.04e-02 67.75% -> 4.69s 273 : 13.650 µm, [8.07e-06; 9.99e-01], 1:0.036s, tol=3.91e-02 68.00% -> 4.65s 274 : 13.700 µm, [8.09e-06; 9.99e-01], 1:0.037s, tol=4.15e-02 68.25% -> 4.66s 275 : 13.750 µm, [8.47e-06; 9.99e-01], 1:0.036s, tol=4.34e-02 68.50% -> 4.57s 276 : 13.800 µm, [8.62e-06; 9.99e-01], 1:0.036s, tol=3.95e-02 68.75% -> 4.53s 277 : 13.850 µm, [8.76e-06; 9.99e-01], 1:0.036s, tol=3.84e-02 69.00% -> 4.49s 278 : 13.900 µm, [8.91e-06; 9.99e-01], 1:0.036s, tol=4.11e-02 69.25% -> 4.47s 279 : 13.950 µm, [9.37e-06; 9.99e-01], 1:0.036s, tol=4.18e-02 69.50% -> 4.42s 280 : 14.000 µm, [9.56e-06; 9.99e-01], 1:0.036s, tol=3.76e-02 69.75% -> 4.38s 281 : 14.050 µm, [9.70e-06; 9.99e-01], 1:0.036s, tol=3.77e-02 70.00% -> 4.33s 282 : 14.100 µm, [1.02e-05; 9.99e-01], 1:0.036s, tol=4.09e-02 70.25% -> 4.31s 283 : 14.150 µm, [1.05e-05; 9.99e-01], 1:0.036s, tol=3.84e-02 70.50% -> 4.26s 284 : 14.200 µm, [1.08e-05; 9.99e-01], 1:0.036s, tol=3.55e-02 70.75% -> 4.23s 285 : 14.250 µm, [1.10e-05; 9.99e-01], 1:0.036s, tol=3.73e-02 71.00% -> 4.17s 286 : 14.300 µm, [1.18e-05; 9.99e-01], 1:0.036s, tol=3.94e-02 71.25% -> 4.15s 287 : 14.350 µm, [1.21e-05; 9.99e-01], 1:0.036s, tol=3.55e-02 71.50% -> 4.10s 288 : 14.400 µm, [1.24e-05; 9.99e-01], 1:0.036s, tol=3.51e-02 71.75% -> 4.07s 289 : 14.450 µm, [1.32e-05; 9.99e-01], 1:0.036s, tol=3.82e-02 72.00% -> 4.03s 290 : 14.500 µm, [1.38e-05; 9.99e-01], 1:0.036s, tol=3.54e-02 72.25% -> 4.02s 291 : 14.550 µm, [1.41e-05; 9.99e-01], 1:0.036s, tol=3.32e-02 72.50% -> 3.98s 292 : 14.600 µm, [1.49e-05; 9.99e-01], 1:0.036s, tol=3.57e-02 72.75% -> 3.96s 293 : 14.650 µm, [1.57e-05; 9.99e-01], 1:0.036s, tol=3.55e-02 73.00% -> 3.92s 294 : 14.700 µm, [1.62e-05; 9.99e-01], 1:0.036s, tol=3.20e-02 73.25% -> 3.89s 295 : 14.750 µm, [1.68e-05; 9.99e-01], 1:0.036s, tol=3.34e-02 73.50% -> 3.86s 296 : 14.800 µm, [1.78e-05; 9.99e-01], 1:0.037s, tol=3.49e-02 73.75% -> 3.84s 297 : 14.850 µm, [1.83e-05; 9.99e-01], 1:0.036s, tol=3.16e-02 74.00% -> 3.79s 298 : 14.900 µm, [1.89e-05; 9.99e-01], 1:0.036s, tol=3.20e-02 74.25% -> 3.76s 299 : 14.950 µm, [1.95e-05; 9.99e-01], 1:0.036s, tol=3.36e-02 74.50% -> 3.71s 300 : 15.000 µm, [2.00e-05; 9.99e-01], 1:0.036s, tol=3.07e-02 74.75% -> 3.68s 301 : 15.050 µm, [2.04e-05; 9.99e-01], 1:0.036s, tol=3.08e-02 75.00% -> 3.63s 302 : 15.100 µm, [2.08e-05; 9.99e-01], 1:0.036s, tol=3.21e-02 75.25% -> 3.60s 303 : 15.150 µm, [2.11e-05; 9.99e-01], 1:0.037s, tol=2.95e-02 75.50% -> 3.60s 304 : 15.200 µm, [2.14e-05; 9.99e-01], 1:0.036s, tol=2.96e-02 75.75% -> 3.52s 305 : 15.250 µm, [2.18e-05; 9.99e-01], 1:0.036s, tol=3.05e-02 76.00% -> 3.47s 306 : 15.300 µm, [2.20e-05; 9.99e-01], 1:0.037s, tol=2.84e-02 76.25% -> 3.47s 307 : 15.350 µm, [2.23e-05; 9.99e-01], 1:0.036s, tol=2.87e-02 76.50% -> 3.40s 308 : 15.400 µm, [2.26e-05; 9.99e-01], 1:0.036s, tol=2.87e-02 76.75% -> 3.38s 309 : 15.450 µm, [2.28e-05; 9.99e-01], 1:0.036s, tol=2.73e-02 77.00% -> 3.32s 310 : 15.500 µm, [2.30e-05; 9.99e-01], 1:0.036s, tol=2.78e-02 77.25% -> 3.29s 311 : 15.550 µm, [2.32e-05; 9.99e-01], 1:0.036s, tol=2.69e-02 77.50% -> 3.25s 312 : 15.600 µm, [2.33e-05; 9.99e-01], 1:0.036s, tol=2.61e-02 77.75% -> 3.22s 313 : 15.650 µm, [2.34e-05; 9.99e-01], 1:0.036s, tol=2.67e-02 78.00% -> 3.18s 314 : 15.700 µm, [2.35e-05; 9.99e-01], 1:0.036s, tol=2.52e-02 78.25% -> 3.15s 315 : 15.750 µm, [2.36e-05; 9.99e-01], 1:0.036s, tol=2.48e-02 78.50% -> 3.10s 316 : 15.800 µm, [2.37e-05; 9.99e-01], 1:0.036s, tol=2.54e-02 78.75% -> 3.07s 317 : 15.850 µm, [2.37e-05; 9.99e-01], 1:0.036s, tol=2.39e-02 79.00% -> 3.06s 318 : 15.900 µm, [2.38e-05; 9.99e-01], 1:0.036s, tol=2.32e-02 79.25% -> 3.00s 319 : 15.950 µm, [2.39e-05; 9.99e-01], 1:0.036s, tol=2.41e-02 79.50% -> 2.98s 320 : 16.000 µm, [2.39e-05; 9.99e-01], 1:0.036s, tol=2.23e-02 79.75% -> 2.95s 321 : 16.050 µm, [2.39e-05; 9.99e-01], 1:0.037s, tol=2.20e-02 80.00% -> 2.93s 322 : 16.100 µm, [2.40e-05; 9.99e-01], 1:0.037s, tol=2.28e-02 80.25% -> 2.89s 323 : 16.150 µm, [2.40e-05; 9.99e-01], 1:0.037s, tol=2.10e-02 80.50% -> 2.85s 324 : 16.200 µm, [2.40e-05; 9.99e-01], 1:0.036s, tol=2.09e-02 80.75% -> 2.80s 325 : 16.250 µm, [2.40e-05; 9.99e-01], 1:0.036s, tol=2.12e-02 81.00% -> 2.76s 326 : 16.300 µm, [2.40e-05; 9.99e-01], 1:0.036s, tol=1.93e-02 81.25% -> 2.73s 327 : 16.350 µm, [2.40e-05; 9.99e-01], 1:0.036s, tol=1.93e-02 81.50% -> 2.70s 328 : 16.400 µm, [2.41e-05; 9.99e-01], 1:0.037s, tol=1.98e-02 81.75% -> 2.67s 329 : 16.450 µm, [2.41e-05; 9.99e-01], 1:0.036s, tol=1.81e-02 82.00% -> 2.60s 330 : 16.500 µm, [2.41e-05; 9.99e-01], 1:0.036s, tol=1.66e-02 82.25% -> 2.58s 331 : 16.550 µm, [2.41e-05; 9.99e-01], 1:0.036s, tol=1.73e-02 82.50% -> 2.54s 332 : 16.600 µm, [2.41e-05; 9.99e-01], 1:0.036s, tol=1.72e-02 82.75% -> 2.50s 333 : 16.650 µm, [2.41e-05; 9.99e-01], 1:0.036s, tol=1.66e-02 83.00% -> 2.45s 334 : 16.700 µm, [2.41e-05; 9.99e-01], 1:0.036s, tol=1.63e-02 83.25% -> 2.44s 335 : 16.750 µm, [2.41e-05; 9.99e-01], 1:0.036s, tol=1.71e-02 83.50% -> 2.40s 336 : 16.800 µm, [2.41e-05; 9.99e-01], 1:0.037s, tol=1.88e-02 83.75% -> 2.38s 337 : 16.850 µm, [2.41e-05; 9.99e-01], 1:0.036s, tol=2.21e-02 84.00% -> 2.31s 338 : 16.900 µm, [2.42e-05; 9.99e-01], 1:0.036s, tol=2.61e-02 84.25% -> 2.28s 339 : 16.950 µm, [2.42e-05; 9.99e-01], 1:0.036s, tol=3.05e-02 84.50% -> 2.23s 340 : 17.000 µm, [2.42e-05; 1.00e+00], 1:0.036s, tol=3.18e-02 84.75% -> 2.21s 341 : 17.050 µm, [2.42e-05; 1.00e+00], 1:0.036s, tol=3.25e-02 85.00% -> 2.17s 342 : 17.100 µm, [2.42e-05; 1.00e+00], 1:0.036s, tol=3.29e-02 85.25% -> 2.13s 343 : 17.150 µm, [2.43e-05; 1.00e+00], 1:0.036s, tol=2.48e-02 85.50% -> 2.09s 344 : 17.200 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=1.33e-02 85.75% -> 2.06s 345 : 17.250 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=8.31e-03 86.00% -> 2.02s 346 : 17.300 µm, [2.43e-05; 1.01e+00], 1:0.037s, tol=6.99e-03 86.25% -> 2.01s 347 : 17.350 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=6.52e-03 86.50% -> 1.96s 348 : 17.400 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=6.25e-03 86.75% -> 1.93s 349 : 17.450 µm, [2.43e-05; 1.01e+00], 1:0.037s, tol=6.08e-03 87.00% -> 1.90s 350 : 17.500 µm, [2.43e-05; 1.01e+00], 1:0.037s, tol=5.96e-03 87.25% -> 1.88s 351 : 17.550 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.88e-03 87.50% -> 1.82s 352 : 17.600 µm, [2.43e-05; 1.01e+00], 1:0.037s, tol=5.83e-03 87.75% -> 1.79s 353 : 17.650 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.79e-03 88.00% -> 1.75s 354 : 17.700 µm, [2.43e-05; 1.01e+00], 1:0.037s, tol=5.76e-03 88.25% -> 1.72s 355 : 17.750 µm, [2.43e-05; 1.01e+00], 1:0.037s, tol=5.74e-03 88.50% -> 1.68s 356 : 17.800 µm, [2.43e-05; 1.01e+00], 1:0.037s, tol=5.71e-03 88.75% -> 1.64s 357 : 17.850 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.69e-03 89.00% -> 1.60s 358 : 17.900 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.67e-03 89.25% -> 1.56s 359 : 17.950 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.65e-03 89.50% -> 1.52s 360 : 18.000 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.64e-03 89.75% -> 1.49s 361 : 18.050 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.62e-03 90.00% -> 1.45s 362 : 18.100 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.60e-03 90.25% -> 1.42s 363 : 18.150 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.59e-03 90.50% -> 1.37s 364 : 18.200 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.58e-03 90.75% -> 1.34s 365 : 18.250 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.57e-03 91.00% -> 1.30s 366 : 18.300 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.55e-03 91.25% -> 1.27s 367 : 18.350 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.54e-03 91.50% -> 1.23s 368 : 18.400 µm, [2.43e-05; 1.01e+00], 1:0.036s, tol=5.52e-03 91.75% -> 1.20s Saved simulation and summary in: /home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/v1.8.1/examples/PhaseField/results/Shear2D/Test/Miehe_AT1_optimMesh Generate movie 01/21 (4.76 %) 2.09 s Generate movie 02/21 (9.52 %) 1.67 s Generate movie 03/21 (14.29 %) 1.55 s Generate movie 04/21 (19.05 %) 1.47 s Generate movie 05/21 (23.81 %) 1.41 s Generate movie 06/21 (28.57 %) 1.29 s Generate movie 07/21 (33.33 %) 1.21 s Generate movie 08/21 (38.10 %) 1.14 s Generate movie 09/21 (42.86 %) 1.05 s Generate movie 10/21 (47.62 %) 987.96 ms Generate movie 11/21 (52.38 %) 872.30 ms Generate movie 12/21 (57.14 %) 800.42 ms Generate movie 13/21 (61.90 %) 705.62 ms Generate movie 14/21 (66.67 %) 607.35 ms Generate movie 15/21 (71.43 %) 516.65 ms Generate movie 16/21 (76.19 %) 434.23 ms Generate movie 17/21 (80.95 %) 345.66 ms Generate movie 18/21 (85.71 %) 260.99 ms Generate movie 19/21 (90.48 %) 177.74 ms Generate movie 20/21 (95.24 %) 91.87 ms Generate movie 21/21 (100.00 %) 0.00 µs Mesh : 45.504 ms Boundary Conditions : 4.525 ms Matrix : 5.019 s Split : 1.410 s Solver : 8.830 s PostProcessing : 247.411 ms Resolution phase field : 13.350 s Display : 432.959 ms PyVista_Interface : 2.913 s | .. code-block:: Python :lineno-start: 14 import matplotlib.pyplot as plt import numpy as np from EasyFEA import ( Display, Folder, Models, Tic, ElemType, Simulations, Paraview, PyVista, ) from EasyFEA.Geoms import Point, Points, Domain, Line, Contour if __name__ == "__main__": Display.Clear() # ---------------------------------------------- # Configuration # ---------------------------------------------- dim = 2 # simu options doSimu = True meshTest = True openCrack = True optimMesh = True # outputs folder = Folder.Results_Dir() + f"{dim}D" plotMesh = False plotEnergy = False makeParaview = False makeMovie = True # phasefield maxIter = 1000 tolConv = 1e-0 # 1e-1, 1e-2, 1e-3 pfmSolver = Models.PhaseField.SolverType.History # Available splits: Bourdin, Amor, Miehe, Stress (isotropic) # He, AnisotStrain, AnisotStress, Zhang (anisotropic) split = Models.PhaseField.SplitType.Miehe # Available regus: AT1, AT2 regu = Models.PhaseField.ReguType.AT1 # ---------------------------------------------- # Geometry # ---------------------------------------------- L = 1e-3 # m l0 = 1e-5 # m thickness = 1 if dim == 2 else 0.1 / 1000 # ---------------------------------------------- # Material # ---------------------------------------------- E = 210e9 # Pa v = 0.3 Gc = 2.7e3 # J/m2 folder_save = Simulations.PhaseField.Folder( folder, "", split, regu, "", tolConv, pfmSolver, meshTest, optimMesh, not openCrack, ) Display.MyPrint(folder_save, "green", end="\n") if doSimu: # ---------------------------------------------- # Mesh # ---------------------------------------------- clC = l0 * 2 if meshTest else l0 / 2 if optimMesh: clD = clC * 4 gap = L * 0.05 h = L if split == Models.PhaseField.SplitType.Bourdin else L / 2 + gap refineDomain = Domain(Point(L / 2 - gap, 0), Point(L, h, thickness), clC) else: clD = clC refineDomain = None 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)] elif 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)] 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], ) # Nodes 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) ) nodes_edges = np.concatenate([nodes_lower, nodes_right, nodes_upper]) # ---------------------------------------------- # Material # ---------------------------------------------- material = Models.Elastic.Isotropic( dim, E=E, v=v, planeStress=False, thickness=thickness ) 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""" E = {E:.2e} Pa; v = {v}; Gc = {Gc:.2e} J/m2; l0 = {l0:.2e} m u_inc = {u_inc:.1e}; N = {N} loadings = np.linspace(u_inc, u_inc*N, N, endpoint=True) for iter, dep in enumerate(loadings): 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()) """ def Loading(dep): 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"]) N = len(loadings) nDetect = 0 list_dep = [] list_f = [] for iter, dep in enumerate(loadings): Loading(dep) u, _, converg = simu.Solve(tolConv, maxIter, convOption=2) simu.Save_Iter() simu.Results_Set_Iteration_Summary(iter, dep * 1e6, "µm", iter / N, True) if not converg: break f = np.sum(simu.Calc_Reaction(dofsX_upper, problemType="elastic")) list_dep.append(dep) list_f.append(f) if simu.Detect_Damage(nodes_edges, 1): nDetect += 1 if nDetect == 10: break # ---------------------------------------------- # Saving # ---------------------------------------------- print() Simulations.Save_pickle((list_f, list_dep), folder_save, "force-displacement") simu.Save(folder_save) else: simu: Simulations.PhaseField = Simulations.Load_Simu(folder_save) list_f, list_dep = Simulations.Load_pickle(folder_save, "force-displacement") # ---------------------------------------------- # Results # ---------------------------------------------- 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) ax = Display.Init_Axes() ax.plot(np.abs(list_dep) * 1e6, np.abs(list_f) * 1e-6, c="blue") ax.set_xlabel("ud [µm]") ax.set_ylabel("f [kN/mm]") ax.grid() Display.Save_fig(folder_save, "force-displacement") if plotMesh: 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 makeParaview: Paraview.Save_simu(simu, folder_save, 400) if makeMovie: simu.Set_Iter(-1) deformFactor = L * 0.05 / simu.Result("displacement_norm").max() iterations = np.arange(0, simu.Niter, simu.Niter // 20) def Func(plotter, iter): simu.Set_Iter(iterations[iter]) thresh = PyVista._pvMesh(simu, "damage", deformFactor).threshold((0, 0.8)) PyVista.Plot(thresh, "damage", 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) plt.show() .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 17.818 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 `_