.. 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-287 .. 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.9.0/examples/PhaseField/results/Shear2D/Test/Miehe_AT1_optimMesh 1: 0.050µm, [0.00e+00;0.00e+00], 1:0.068s, tol=1.00e+00 2: 0.100µm, [0.00e+00;0.00e+00], 1:0.059s, tol=7.50e-01 0.25% -> 23.47s 3: 0.150µm, [0.00e+00;0.00e+00], 1:0.058s, tol=5.56e-01 0.50% -> 23.11s 4: 0.200µm, [0.00e+00;0.00e+00], 1:0.058s, tol=4.38e-01 0.75% -> 23.20s 5: 0.250µm, [0.00e+00;0.00e+00], 1:0.059s, tol=3.60e-01 1.00% -> 23.35s 6: 0.300µm, [0.00e+00;0.00e+00], 1:0.059s, tol=3.06e-01 1.25% -> 23.22s 7: 0.350µm, [0.00e+00;0.00e+00], 1:0.059s, tol=2.65e-01 1.50% -> 23.06s 8: 0.400µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.34e-01 1.75% -> 22.90s 9: 0.450µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.10e-01 2.00% -> 22.91s 10: 0.500µm, [0.00e+00;0.00e+00], 1:0.059s, tol=1.90e-01 2.25% -> 22.99s 11: 0.550µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.74e-01 2.50% -> 22.47s 12: 0.600µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.60e-01 2.75% -> 22.47s 13: 0.650µm, [0.00e+00;0.00e+00], 1:0.059s, tol=1.48e-01 3.00% -> 22.84s 14: 0.700µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.38e-01 3.25% -> 22.61s 15: 0.750µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.29e-01 3.50% -> 22.42s 16: 0.800µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.21e-01 3.75% -> 22.50s 17: 0.850µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.14e-01 4.00% -> 22.33s 18: 0.900µm, [0.00e+00;0.00e+00], 1:0.059s, tol=1.08e-01 4.25% -> 22.44s 19: 0.950µm, [0.00e+00;0.00e+00], 1:0.059s, tol=1.02e-01 4.50% -> 22.44s 20: 1.000µm, [0.00e+00;0.00e+00], 1:0.058s, tol=9.75e-02 4.75% -> 22.21s 21: 1.050µm, [0.00e+00;0.00e+00], 1:0.058s, tol=9.30e-02 5.00% -> 22.07s 22: 1.100µm, [0.00e+00;0.00e+00], 1:0.057s, tol=8.88e-02 5.25% -> 21.73s 23: 1.150µm, [0.00e+00;0.00e+00], 1:0.058s, tol=8.51e-02 5.50% -> 21.77s 24: 1.200µm, [0.00e+00;0.00e+00], 1:0.058s, tol=8.16e-02 5.75% -> 21.89s 25: 1.250µm, [0.00e+00;0.00e+00], 1:0.058s, tol=7.84e-02 6.00% -> 21.80s 26: 1.300µm, [0.00e+00;0.00e+00], 1:0.058s, tol=7.54e-02 6.25% -> 21.93s 27: 1.350µm, [0.00e+00;0.00e+00], 1:0.058s, tol=7.27e-02 6.50% -> 21.80s 28: 1.400µm, [0.00e+00;0.00e+00], 1:0.058s, tol=7.02e-02 6.75% -> 21.57s 29: 1.450µm, [0.00e+00;0.00e+00], 1:0.058s, tol=6.78e-02 7.00% -> 21.47s 30: 1.500µm, [0.00e+00;0.00e+00], 1:0.058s, tol=6.56e-02 7.25% -> 21.40s 31: 1.550µm, [0.00e+00;0.00e+00], 1:0.058s, tol=6.35e-02 7.50% -> 21.50s 32: 1.600µm, [0.00e+00;0.00e+00], 1:0.058s, tol=6.15e-02 7.75% -> 21.48s 33: 1.650µm, [0.00e+00;0.00e+00], 1:0.058s, tol=5.97e-02 8.00% -> 21.30s 34: 1.700µm, [0.00e+00;0.00e+00], 1:0.058s, tol=5.80e-02 8.25% -> 21.17s 35: 1.750µm, [0.00e+00;0.00e+00], 1:0.058s, tol=5.63e-02 8.50% -> 21.30s 36: 1.800µm, [0.00e+00;0.00e+00], 1:0.058s, tol=5.48e-02 8.75% -> 21.30s 37: 1.850µm, [0.00e+00;0.00e+00], 1:0.058s, tol=5.33e-02 9.00% -> 21.28s 38: 1.900µm, [0.00e+00;0.00e+00], 1:0.058s, tol=5.19e-02 9.25% -> 21.02s 39: 1.950µm, [0.00e+00;0.00e+00], 1:0.058s, tol=5.06e-02 9.50% -> 21.01s 40: 2.000µm, [0.00e+00;0.00e+00], 1:0.058s, tol=4.94e-02 9.75% -> 20.97s 41: 2.050µm, [0.00e+00;0.00e+00], 1:0.058s, tol=4.82e-02 10.00% -> 20.80s 42: 2.100µm, [0.00e+00;0.00e+00], 1:0.058s, tol=4.71e-02 10.25% -> 20.90s 43: 2.150µm, [0.00e+00;0.00e+00], 1:0.058s, tol=4.60e-02 10.50% -> 20.65s 44: 2.200µm, [0.00e+00;0.00e+00], 1:0.058s, tol=4.49e-02 10.75% -> 20.64s 45: 2.250µm, [0.00e+00;0.00e+00], 1:0.058s, tol=4.40e-02 11.00% -> 20.48s 46: 2.300µm, [0.00e+00;0.00e+00], 1:0.057s, tol=4.30e-02 11.25% -> 20.33s 47: 2.350µm, [0.00e+00;0.00e+00], 1:0.058s, tol=4.21e-02 11.50% -> 20.68s 48: 2.400µm, [0.00e+00;0.00e+00], 1:0.058s, tol=4.12e-02 11.75% -> 20.40s 49: 2.450µm, [0.00e+00;0.00e+00], 1:0.058s, tol=4.04e-02 12.00% -> 20.25s 50: 2.500µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.96e-02 12.25% -> 20.27s 51: 2.550µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.88e-02 12.50% -> 20.17s 52: 2.600µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.81e-02 12.75% -> 20.34s 53: 2.650µm, [0.00e+00;0.00e+00], 1:0.059s, tol=3.74e-02 13.00% -> 20.37s 54: 2.700µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.67e-02 13.25% -> 20.26s 55: 2.750µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.60e-02 13.50% -> 20.05s 56: 2.800µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.54e-02 13.75% -> 19.86s 57: 2.850µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.48e-02 14.00% -> 19.87s 58: 2.900µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.42e-02 14.25% -> 19.91s 59: 2.950µm, [0.00e+00;0.00e+00], 1:0.059s, tol=3.36e-02 14.50% -> 20.02s 60: 3.000µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.31e-02 14.75% -> 19.67s 61: 3.050µm, [0.00e+00;0.00e+00], 1:0.059s, tol=3.25e-02 15.00% -> 19.92s 62: 3.100µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.20e-02 15.25% -> 19.58s 63: 3.150µm, [0.00e+00;0.00e+00], 1:0.057s, tol=3.15e-02 15.50% -> 19.41s 64: 3.200µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.10e-02 15.75% -> 19.71s 65: 3.250µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.05e-02 16.00% -> 19.50s 66: 3.300µm, [0.00e+00;0.00e+00], 1:0.058s, tol=3.01e-02 16.25% -> 19.51s 67: 3.350µm, [0.00e+00;0.00e+00], 1:0.059s, tol=2.96e-02 16.50% -> 19.58s 68: 3.400µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.92e-02 16.75% -> 19.20s 69: 3.450µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.88e-02 17.00% -> 19.23s 70: 3.500µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.84e-02 17.25% -> 19.17s 71: 3.550µm, [0.00e+00;0.00e+00], 1:0.059s, tol=2.80e-02 17.50% -> 19.38s 72: 3.600µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.76e-02 17.75% -> 19.12s 73: 3.650µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.72e-02 18.00% -> 19.02s 74: 3.700µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.68e-02 18.25% -> 18.99s 75: 3.750µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.65e-02 18.50% -> 18.89s 76: 3.800µm, [0.00e+00;0.00e+00], 1:0.059s, tol=2.61e-02 18.75% -> 19.20s 77: 3.850µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.58e-02 19.00% -> 18.76s 78: 3.900µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.55e-02 19.25% -> 18.78s 79: 3.950µm, [0.00e+00;0.00e+00], 1:0.057s, tol=2.52e-02 19.50% -> 18.46s 80: 4.000µm, [0.00e+00;0.00e+00], 1:0.057s, tol=2.48e-02 19.75% -> 18.41s 81: 4.050µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.45e-02 20.00% -> 18.61s 82: 4.100µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.42e-02 20.25% -> 18.38s 83: 4.150µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.40e-02 20.50% -> 18.34s 84: 4.200µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.37e-02 20.75% -> 18.24s 85: 4.250µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.34e-02 21.00% -> 18.36s 86: 4.300µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.31e-02 21.25% -> 18.36s 87: 4.350µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.29e-02 21.50% -> 18.30s 88: 4.400µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.26e-02 21.75% -> 18.15s 89: 4.450µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.23e-02 22.00% -> 18.15s 90: 4.500µm, [0.00e+00;0.00e+00], 1:0.059s, tol=2.21e-02 22.25% -> 18.20s 91: 4.550µm, [0.00e+00;0.00e+00], 1:0.059s, tol=2.19e-02 22.50% -> 18.17s 92: 4.600µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.16e-02 22.75% -> 18.01s 93: 4.650µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.14e-02 23.00% -> 17.96s 94: 4.700µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.12e-02 23.25% -> 17.73s 95: 4.750µm, [0.00e+00;0.00e+00], 1:0.057s, tol=2.09e-02 23.50% -> 17.57s 96: 4.800µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.07e-02 23.75% -> 17.58s 97: 4.850µm, [0.00e+00;0.00e+00], 1:0.057s, tol=2.05e-02 24.00% -> 17.46s 98: 4.900µm, [0.00e+00;0.00e+00], 1:0.058s, tol=2.03e-02 24.25% -> 17.55s 99: 4.950µm, [0.00e+00;0.00e+00], 1:0.057s, tol=2.01e-02 24.50% -> 17.36s 100: 5.000µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.99e-02 24.75% -> 17.45s 101: 5.050µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.97e-02 25.00% -> 17.43s 102: 5.100µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.95e-02 25.25% -> 17.32s 103: 5.150µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.93e-02 25.50% -> 17.17s 104: 5.200µm, [0.00e+00;0.00e+00], 1:0.074s, tol=1.91e-02 25.75% -> 21.91s 105: 5.250µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.90e-02 26.00% -> 17.13s 106: 5.300µm, [0.00e+00;0.00e+00], 1:0.059s, tol=1.88e-02 26.25% -> 17.26s 107: 5.350µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.86e-02 26.50% -> 17.19s 108: 5.400µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.84e-02 26.75% -> 17.14s 109: 5.450µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.83e-02 27.00% -> 17.03s 110: 5.500µm, [0.00e+00;0.00e+00], 1:0.059s, tol=1.81e-02 27.25% -> 17.04s 111: 5.550µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.79e-02 27.50% -> 16.76s 112: 5.600µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.78e-02 27.75% -> 16.82s 113: 5.650µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.76e-02 28.00% -> 16.61s 114: 5.700µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.75e-02 28.25% -> 16.68s 115: 5.750µm, [0.00e+00;0.00e+00], 1:0.057s, tol=1.73e-02 28.50% -> 16.40s 116: 5.800µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.72e-02 28.75% -> 16.49s 117: 5.850µm, [0.00e+00;0.00e+00], 1:0.057s, tol=1.70e-02 29.00% -> 16.29s 118: 5.900µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.69e-02 29.25% -> 16.55s 119: 5.950µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.67e-02 29.50% -> 16.45s 120: 6.000µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.66e-02 29.75% -> 16.33s 121: 6.050µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.65e-02 30.00% -> 16.21s 122: 6.100µm, [0.00e+00;0.00e+00], 1:0.057s, tol=1.63e-02 30.25% -> 16.02s 123: 6.150µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.62e-02 30.50% -> 16.11s 124: 6.200µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.61e-02 30.75% -> 16.14s 125: 6.250µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.59e-02 31.00% -> 15.97s 126: 6.300µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.58e-02 31.25% -> 15.91s 127: 6.350µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.57e-02 31.50% -> 15.83s 128: 6.400µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.56e-02 31.75% -> 15.87s 129: 6.450µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.54e-02 32.00% -> 15.87s 130: 6.500µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.53e-02 32.25% -> 15.81s 131: 6.550µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.52e-02 32.50% -> 15.66s 132: 6.600µm, [0.00e+00;0.00e+00], 1:0.057s, tol=1.51e-02 32.75% -> 15.40s 133: 6.650µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.50e-02 33.00% -> 15.54s 134: 6.700µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.49e-02 33.25% -> 15.36s 135: 6.750µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.48e-02 33.50% -> 15.42s 136: 6.800µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.47e-02 33.75% -> 15.30s 137: 6.850µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.45e-02 34.00% -> 15.28s 138: 6.900µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.44e-02 34.25% -> 15.21s 139: 6.950µm, [0.00e+00;0.00e+00], 1:0.058s, tol=1.43e-02 34.50% -> 15.28s 140: 7.000µm, [6.00e-07;1.59e-03], 1:0.058s, tol=1.42e-02 34.75% -> 15.16s 141: 7.050µm, [2.22e-06;6.20e-03], 1:0.059s, tol=1.39e-02 35.00% -> 15.21s 142: 7.100µm, [4.16e-06;1.23e-02], 1:0.059s, tol=1.38e-02 35.25% -> 15.25s 143: 7.150µm, [6.13e-06;1.92e-02], 1:0.060s, tol=1.37e-02 35.50% -> 15.51s 144: 7.200µm, [7.97e-06;2.65e-02], 1:0.058s, tol=1.35e-02 35.75% -> 15.03s 145: 7.250µm, [9.65e-06;3.41e-02], 1:0.063s, tol=1.35e-02 36.00% -> 16.00s 146: 7.300µm, [1.12e-05;4.18e-02], 1:0.059s, tol=1.34e-02 36.25% -> 14.97s 147: 7.350µm, [1.25e-05;4.96e-02], 1:0.059s, tol=1.33e-02 36.50% -> 14.90s 148: 7.400µm, [1.36e-05;5.75e-02], 1:0.059s, tol=1.32e-02 36.75% -> 14.98s 149: 7.450µm, [1.46e-05;6.55e-02], 1:0.059s, tol=1.31e-02 37.00% -> 14.96s 150: 7.500µm, [1.55e-05;7.35e-02], 1:0.058s, tol=1.30e-02 37.25% -> 14.64s 151: 7.550µm, [1.62e-05;8.16e-02], 1:0.058s, tol=1.29e-02 37.50% -> 14.56s 152: 7.600µm, [1.68e-05;8.98e-02], 1:0.061s, tol=1.28e-02 37.75% -> 15.08s 153: 7.650µm, [1.72e-05;9.80e-02], 1:0.059s, tol=1.28e-02 38.00% -> 14.67s 154: 7.700µm, [1.76e-05;1.06e-01], 1:0.059s, tol=1.27e-02 38.25% -> 14.67s 155: 7.750µm, [1.79e-05;1.15e-01], 1:0.059s, tol=1.26e-02 38.50% -> 14.40s 156: 7.800µm, [1.81e-05;1.23e-01], 1:0.058s, tol=1.25e-02 38.75% -> 14.22s 157: 7.850µm, [1.82e-05;1.32e-01], 1:0.058s, tol=1.24e-02 39.00% -> 14.24s 158: 7.900µm, [1.83e-05;1.40e-01], 1:0.060s, tol=1.24e-02 39.25% -> 14.55s 159: 7.950µm, [1.83e-05;1.49e-01], 1:0.058s, tol=1.23e-02 39.50% -> 14.03s 160: 8.000µm, [1.82e-05;1.58e-01], 1:0.059s, tol=1.22e-02 39.75% -> 14.33s 161: 8.050µm, [1.81e-05;1.67e-01], 1:0.059s, tol=1.21e-02 40.00% -> 14.06s 162: 8.100µm, [1.80e-05;1.76e-01], 1:0.059s, tol=1.21e-02 40.25% -> 14.11s 163: 8.150µm, [1.78e-05;1.85e-01], 1:0.058s, tol=1.20e-02 40.50% -> 13.88s 164: 8.200µm, [1.75e-05;1.94e-01], 1:0.059s, tol=1.19e-02 40.75% -> 14.09s 165: 8.250µm, [1.73e-05;2.03e-01], 1:0.059s, tol=1.18e-02 41.00% -> 13.84s 166: 8.300µm, [1.70e-05;2.12e-01], 1:0.058s, tol=1.18e-02 41.25% -> 13.58s 167: 8.350µm, [1.67e-05;2.22e-01], 1:0.058s, tol=1.17e-02 41.50% -> 13.61s 168: 8.400µm, [1.64e-05;2.31e-01], 1:0.059s, tol=1.16e-02 41.75% -> 13.64s 169: 8.450µm, [1.60e-05;2.41e-01], 1:0.058s, tol=1.16e-02 42.00% -> 13.56s 170: 8.500µm, [1.57e-05;2.51e-01], 1:0.059s, tol=1.15e-02 42.25% -> 13.57s 171: 8.550µm, [1.53e-05;2.60e-01], 1:0.059s, tol=1.14e-02 42.50% -> 13.47s 172: 8.600µm, [1.49e-05;2.70e-01], 1:0.058s, tol=1.14e-02 42.75% -> 13.36s 173: 8.650µm, [1.45e-05;2.80e-01], 1:0.058s, tol=1.13e-02 43.00% -> 13.19s 174: 8.700µm, [1.42e-05;2.90e-01], 1:0.059s, tol=1.13e-02 43.25% -> 13.29s 175: 8.750µm, [1.38e-05;3.01e-01], 1:0.062s, tol=1.12e-02 43.50% -> 13.91s 176: 8.800µm, [1.34e-05;3.11e-01], 1:0.060s, tol=1.11e-02 43.75% -> 13.49s 177: 8.850µm, [1.30e-05;3.21e-01], 1:0.058s, tol=1.11e-02 44.00% -> 13.08s 178: 8.900µm, [1.26e-05;3.32e-01], 1:0.058s, tol=1.10e-02 44.25% -> 13.02s 179: 8.950µm, [1.22e-05;3.43e-01], 1:0.059s, tol=1.10e-02 44.50% -> 13.00s 180: 9.000µm, [1.19e-05;3.56e-01], 1:0.058s, tol=1.09e-02 44.75% -> 12.84s 181: 9.050µm, [1.15e-05;3.71e-01], 1:0.058s, tol=1.08e-02 45.00% -> 12.80s 182: 9.100µm, [1.12e-05;3.86e-01], 1:0.059s, tol=1.07e-02 45.25% -> 12.99s 183: 9.150µm, [1.09e-05;4.01e-01], 1:0.059s, tol=1.07e-02 45.50% -> 12.83s 184: 9.200µm, [1.05e-05;4.17e-01], 1:0.059s, tol=1.06e-02 45.75% -> 12.74s 185: 9.250µm, [1.02e-05;4.32e-01], 1:0.062s, tol=1.06e-02 46.00% -> 13.33s 186: 9.300µm, [9.84e-06;4.47e-01], 1:0.059s, tol=1.06e-02 46.25% -> 12.75s 187: 9.350µm, [9.49e-06;4.62e-01], 1:0.059s, tol=1.05e-02 46.50% -> 12.60s 188: 9.400µm, [9.27e-06;4.77e-01], 1:0.059s, tol=1.04e-02 46.75% -> 12.63s 189: 9.450µm, [9.10e-06;4.94e-01], 1:0.058s, tol=1.03e-02 47.00% -> 12.38s 190: 9.500µm, [8.97e-06;5.14e-01], 1:0.058s, tol=1.03e-02 47.25% -> 12.28s 191: 9.550µm, [8.88e-06;5.39e-01], 1:0.060s, tol=1.02e-02 47.50% -> 12.54s 192: 9.600µm, [8.83e-06;5.71e-01], 1:0.059s, tol=1.03e-02 47.75% -> 12.29s 193: 9.650µm, [8.78e-06;6.08e-01], 1:0.059s, tol=1.04e-02 48.00% -> 12.32s 194: 9.700µm, [8.76e-06;6.49e-01], 1:0.059s, tol=1.07e-02 48.25% -> 12.18s 195: 9.750µm, [8.75e-06;6.95e-01], 1:0.058s, tol=1.13e-02 48.50% -> 11.98s 196: 9.800µm, [9.05e-06;7.43e-01], 1:0.058s, tol=1.23e-02 48.75% -> 11.95s 197: 9.850µm, [9.39e-06;7.94e-01], 1:0.060s, tol=1.47e-02 49.00% -> 12.16s 198: 9.900µm, [9.76e-06;8.46e-01], 1:0.059s, tol=1.92e-02 49.25% -> 11.93s 199: 9.950µm, [9.84e-06;8.94e-01], 1:0.058s, tol=2.57e-02 49.50% -> 11.81s 200: 10.000µm, [9.59e-06;9.28e-01], 1:0.060s, tol=3.11e-02 49.75% -> 12.01s 201: 10.050µm, [9.14e-06;9.49e-01], 1:0.059s, tol=3.15e-02 50.00% -> 11.78s 202: 10.100µm, [9.06e-06;9.58e-01], 1:0.059s, tol=2.98e-02 50.25% -> 11.64s 203: 10.150µm, [9.00e-06;9.60e-01], 1:0.059s, tol=3.01e-02 50.50% -> 11.66s 204: 10.200µm, [9.21e-06;9.60e-01], 1:0.059s, tol=3.37e-02 50.75% -> 11.54s 205: 10.250µm, [9.27e-06;9.60e-01], 1:0.059s, tol=4.01e-02 51.00% -> 11.53s 206: 10.300µm, [9.16e-06;9.71e-01], 1:0.059s, tol=4.63e-02 51.25% -> 11.50s 207: 10.350µm, [8.87e-06;9.82e-01], 1:0.058s, tol=4.81e-02 51.50% -> 11.29s 208: 10.400µm, [8.45e-06;9.89e-01], 1:0.062s, tol=4.32e-02 51.75% -> 11.88s 209: 10.450µm, [8.11e-06;9.92e-01], 1:0.060s, tol=3.51e-02 52.00% -> 11.43s 210: 10.500µm, [8.09e-06;9.94e-01], 1:0.058s, tol=3.08e-02 52.25% -> 11.15s 211: 10.550µm, [8.04e-06;9.95e-01], 1:0.058s, tol=2.91e-02 52.50% -> 11.06s 212: 10.600µm, [7.90e-06;9.95e-01], 1:0.058s, tol=2.80e-02 52.75% -> 10.97s 213: 10.650µm, [7.80e-06;9.95e-01], 1:0.058s, tol=2.75e-02 53.00% -> 10.93s 214: 10.700µm, [7.73e-06;9.94e-01], 1:0.058s, tol=2.85e-02 53.25% -> 10.87s 215: 10.750µm, [7.62e-06;9.94e-01], 1:0.059s, tol=3.07e-02 53.50% -> 11.02s 216: 10.800µm, [7.45e-06;9.93e-01], 1:0.059s, tol=3.36e-02 53.75% -> 10.90s 217: 10.850µm, [7.23e-06;9.93e-01], 1:0.059s, tol=3.59e-02 54.00% -> 10.82s 218: 10.900µm, [6.95e-06;9.92e-01], 1:0.061s, tol=3.61e-02 54.25% -> 11.12s 219: 10.950µm, [6.71e-06;9.92e-01], 1:0.059s, tol=3.49e-02 54.50% -> 10.66s 220: 11.000µm, [6.54e-06;9.92e-01], 1:0.059s, tol=3.27e-02 54.75% -> 10.65s 221: 11.050µm, [6.47e-06;9.92e-01], 1:0.059s, tol=3.08e-02 55.00% -> 10.53s 222: 11.100µm, [6.45e-06;9.92e-01], 1:0.059s, tol=3.00e-02 55.25% -> 10.49s 223: 11.150µm, [6.36e-06;9.93e-01], 1:0.059s, tol=2.96e-02 55.50% -> 10.54s 224: 11.200µm, [6.21e-06;9.93e-01], 1:0.061s, tol=2.93e-02 55.75% -> 10.77s 225: 11.250µm, [6.16e-06;9.93e-01], 1:0.058s, tol=3.12e-02 56.00% -> 10.29s 226: 11.300µm, [6.13e-06;9.93e-01], 1:0.059s, tol=3.48e-02 56.25% -> 10.41s 227: 11.350µm, [6.03e-06;9.94e-01], 1:0.059s, tol=3.87e-02 56.50% -> 10.20s 228: 11.400µm, [5.87e-06;9.94e-01], 1:0.059s, tol=4.11e-02 56.75% -> 10.13s 229: 11.450µm, [5.71e-06;9.94e-01], 1:0.059s, tol=4.15e-02 57.00% -> 10.09s 230: 11.500µm, [5.65e-06;9.95e-01], 1:0.060s, tol=4.01e-02 57.25% -> 10.19s 231: 11.550µm, [5.70e-06;9.95e-01], 1:0.058s, tol=3.84e-02 57.50% -> 9.94s 232: 11.600µm, [5.74e-06;9.95e-01], 1:0.058s, tol=3.71e-02 57.75% -> 9.87s 233: 11.650µm, [5.66e-06;9.95e-01], 1:0.059s, tol=3.53e-02 58.00% -> 9.92s 234: 11.700µm, [5.61e-06;9.95e-01], 1:0.059s, tol=3.48e-02 58.25% -> 9.77s 235: 11.750µm, [5.60e-06;9.95e-01], 1:0.058s, tol=3.73e-02 58.50% -> 9.68s 236: 11.800µm, [5.51e-06;9.96e-01], 1:0.059s, tol=4.00e-02 58.75% -> 9.74s 237: 11.850µm, [5.41e-06;9.96e-01], 1:0.059s, tol=4.25e-02 59.00% -> 9.64s 238: 11.900µm, [5.38e-06;9.96e-01], 1:0.058s, tol=4.26e-02 59.25% -> 9.53s 239: 11.950µm, [5.51e-06;9.96e-01], 1:0.058s, tol=4.20e-02 59.50% -> 9.39s 240: 12.000µm, [5.60e-06;9.96e-01], 1:0.059s, tol=4.03e-02 59.75% -> 9.47s 241: 12.050µm, [5.53e-06;9.96e-01], 1:0.060s, tol=3.75e-02 60.00% -> 9.64s 242: 12.100µm, [5.56e-06;9.96e-01], 1:0.060s, tol=3.75e-02 60.25% -> 9.48s 243: 12.150µm, [5.52e-06;9.97e-01], 1:0.059s, tol=3.99e-02 60.50% -> 9.29s 244: 12.200µm, [5.47e-06;9.97e-01], 1:0.059s, tol=4.29e-02 60.75% -> 9.24s 245: 12.250µm, [5.50e-06;9.97e-01], 1:0.059s, tol=4.42e-02 61.00% -> 9.14s 246: 12.300µm, [5.73e-06;9.97e-01], 1:0.058s, tol=4.43e-02 61.25% -> 9.01s 247: 12.350µm, [5.83e-06;9.97e-01], 1:0.058s, tol=4.23e-02 61.50% -> 8.96s 248: 12.400µm, [5.81e-06;9.97e-01], 1:0.060s, tol=3.95e-02 61.75% -> 9.17s 249: 12.450µm, [5.86e-06;9.97e-01], 1:0.058s, tol=4.02e-02 62.00% -> 8.85s 250: 12.500µm, [5.85e-06;9.97e-01], 1:0.058s, tol=4.32e-02 62.25% -> 8.78s 251: 12.550µm, [5.91e-06;9.97e-01], 1:0.061s, tol=4.52e-02 62.50% -> 9.11s 252: 12.600µm, [6.26e-06;9.97e-01], 1:0.059s, tol=4.57e-02 62.75% -> 8.73s 253: 12.650µm, [6.37e-06;9.97e-01], 1:0.059s, tol=4.31e-02 63.00% -> 8.67s 254: 12.700µm, [6.44e-06;9.98e-01], 1:0.059s, tol=4.07e-02 63.25% -> 8.63s 255: 12.750µm, [6.46e-06;9.98e-01], 1:0.059s, tol=4.15e-02 63.50% -> 8.57s 256: 12.800µm, [6.47e-06;9.98e-01], 1:0.059s, tol=4.43e-02 63.75% -> 8.55s 257: 12.850µm, [6.91e-06;9.98e-01], 1:0.060s, tol=4.65e-02 64.00% -> 8.62s 258: 12.900µm, [7.01e-06;9.98e-01], 1:0.059s, tol=4.45e-02 64.25% -> 8.37s 259: 12.950µm, [6.99e-06;9.98e-01], 1:0.059s, tol=4.13e-02 64.50% -> 8.33s 260: 13.000µm, [6.98e-06;9.98e-01], 1:0.059s, tol=4.14e-02 64.75% -> 8.30s 261: 13.050µm, [6.91e-06;9.98e-01], 1:0.059s, tol=4.41e-02 65.00% -> 8.24s 262: 13.100µm, [7.11e-06;9.98e-01], 1:0.059s, tol=4.64e-02 65.25% -> 8.16s 263: 13.150µm, [7.16e-06;9.98e-01], 1:0.059s, tol=4.35e-02 65.50% -> 8.14s 264: 13.200µm, [7.20e-06;9.98e-01], 1:0.059s, tol=4.08e-02 65.75% -> 8.03s 265: 13.250µm, [7.21e-06;9.98e-01], 1:0.059s, tol=4.22e-02 66.00% -> 8.00s 266: 13.300µm, [7.25e-06;9.98e-01], 1:0.059s, tol=4.48e-02 66.25% -> 7.92s 267: 13.350µm, [7.50e-06;9.98e-01], 1:0.059s, tol=4.48e-02 66.50% -> 7.89s 268: 13.400µm, [7.52e-06;9.98e-01], 1:0.058s, tol=4.09e-02 66.75% -> 7.74s 269: 13.450µm, [7.59e-06;9.98e-01], 1:0.058s, tol=4.04e-02 67.00% -> 7.68s 270: 13.500µm, [7.57e-06;9.98e-01], 1:0.058s, tol=4.28e-02 67.25% -> 7.63s 271: 13.550µm, [7.89e-06;9.98e-01], 1:0.059s, tol=4.45e-02 67.50% -> 7.62s 272: 13.600µm, [7.94e-06;9.98e-01], 1:0.058s, tol=4.04e-02 67.75% -> 7.52s 273: 13.650µm, [8.07e-06;9.99e-01], 1:0.059s, tol=3.91e-02 68.00% -> 7.55s 274: 13.700µm, [8.09e-06;9.99e-01], 1:0.062s, tol=4.15e-02 68.25% -> 7.82s 275: 13.750µm, [8.47e-06;9.99e-01], 1:0.060s, tol=4.34e-02 68.50% -> 7.59s 276: 13.800µm, [8.62e-06;9.99e-01], 1:0.059s, tol=3.95e-02 68.75% -> 7.35s 277: 13.850µm, [8.76e-06;9.99e-01], 1:0.058s, tol=3.84e-02 69.00% -> 7.21s 278: 13.900µm, [8.91e-06;9.99e-01], 1:0.058s, tol=4.11e-02 69.25% -> 7.19s 279: 13.950µm, [9.37e-06;9.99e-01], 1:0.059s, tol=4.18e-02 69.50% -> 7.18s 280: 14.000µm, [9.56e-06;9.99e-01], 1:0.059s, tol=3.76e-02 69.75% -> 7.11s 281: 14.050µm, [9.70e-06;9.99e-01], 1:0.059s, tol=3.77e-02 70.00% -> 7.13s 282: 14.100µm, [1.02e-05;9.99e-01], 1:0.058s, tol=4.09e-02 70.25% -> 6.95s 283: 14.150µm, [1.05e-05;9.99e-01], 1:0.058s, tol=3.84e-02 70.50% -> 6.86s 284: 14.200µm, [1.08e-05;9.99e-01], 1:0.061s, tol=3.55e-02 70.75% -> 7.18s 285: 14.250µm, [1.10e-05;9.99e-01], 1:0.059s, tol=3.73e-02 71.00% -> 6.80s 286: 14.300µm, [1.18e-05;9.99e-01], 1:0.059s, tol=3.94e-02 71.25% -> 6.80s 287: 14.350µm, [1.21e-05;9.99e-01], 1:0.059s, tol=3.55e-02 71.50% -> 6.69s 288: 14.400µm, [1.24e-05;9.99e-01], 1:0.059s, tol=3.51e-02 71.75% -> 6.67s 289: 14.450µm, [1.32e-05;9.99e-01], 1:0.059s, tol=3.82e-02 72.00% -> 6.57s 290: 14.500µm, [1.38e-05;9.99e-01], 1:0.060s, tol=3.54e-02 72.25% -> 6.68s 291: 14.550µm, [1.41e-05;9.99e-01], 1:0.058s, tol=3.32e-02 72.50% -> 6.39s 292: 14.600µm, [1.49e-05;9.99e-01], 1:0.059s, tol=3.57e-02 72.75% -> 6.45s 293: 14.650µm, [1.57e-05;9.99e-01], 1:0.059s, tol=3.55e-02 73.00% -> 6.33s 294: 14.700µm, [1.62e-05;9.99e-01], 1:0.059s, tol=3.20e-02 73.25% -> 6.32s 295: 14.750µm, [1.68e-05;9.99e-01], 1:0.059s, tol=3.34e-02 73.50% -> 6.21s 296: 14.800µm, [1.78e-05;9.99e-01], 1:0.060s, tol=3.49e-02 73.75% -> 6.33s 297: 14.850µm, [1.83e-05;9.99e-01], 1:0.059s, tol=3.16e-02 74.00% -> 6.12s 298: 14.900µm, [1.89e-05;9.99e-01], 1:0.058s, tol=3.20e-02 74.25% -> 6.02s 299: 14.950µm, [1.95e-05;9.99e-01], 1:0.058s, tol=3.36e-02 74.50% -> 5.93s 300: 15.000µm, [2.00e-05;9.99e-01], 1:0.058s, tol=3.07e-02 74.75% -> 5.88s 301: 15.050µm, [2.04e-05;9.99e-01], 1:0.058s, tol=3.08e-02 75.00% -> 5.83s 302: 15.100µm, [2.08e-05;9.99e-01], 1:0.059s, tol=3.21e-02 75.25% -> 5.86s 303: 15.150µm, [2.11e-05;9.99e-01], 1:0.059s, tol=2.95e-02 75.50% -> 5.80s 304: 15.200µm, [2.14e-05;9.99e-01], 1:0.059s, tol=2.96e-02 75.75% -> 5.74s 305: 15.250µm, [2.18e-05;9.99e-01], 1:0.058s, tol=3.05e-02 76.00% -> 5.61s 306: 15.300µm, [2.20e-05;9.99e-01], 1:0.059s, tol=2.84e-02 76.25% -> 5.60s 307: 15.350µm, [2.23e-05;9.99e-01], 1:0.061s, tol=2.87e-02 76.50% -> 5.69s 308: 15.400µm, [2.26e-05;9.99e-01], 1:0.060s, tol=2.87e-02 76.75% -> 5.57s 309: 15.450µm, [2.28e-05;9.99e-01], 1:0.058s, tol=2.73e-02 77.00% -> 5.36s 310: 15.500µm, [2.30e-05;9.99e-01], 1:0.059s, tol=2.78e-02 77.25% -> 5.33s 311: 15.550µm, [2.32e-05;9.99e-01], 1:0.058s, tol=2.69e-02 77.50% -> 5.25s 312: 15.600µm, [2.33e-05;9.99e-01], 1:0.060s, tol=2.61e-02 77.75% -> 5.31s 313: 15.650µm, [2.34e-05;9.99e-01], 1:0.058s, tol=2.67e-02 78.00% -> 5.13s 314: 15.700µm, [2.35e-05;9.99e-01], 1:0.060s, tol=2.52e-02 78.25% -> 5.18s 315: 15.750µm, [2.36e-05;9.99e-01], 1:0.059s, tol=2.48e-02 78.50% -> 5.08s 316: 15.800µm, [2.37e-05;9.99e-01], 1:0.060s, tol=2.54e-02 78.75% -> 5.11s 317: 15.850µm, [2.37e-05;9.99e-01], 1:0.062s, tol=2.39e-02 79.00% -> 5.18s 318: 15.900µm, [2.38e-05;9.99e-01], 1:0.058s, tol=2.32e-02 79.25% -> 4.86s 319: 15.950µm, [2.39e-05;9.99e-01], 1:0.058s, tol=2.41e-02 79.50% -> 4.80s 320: 16.000µm, [2.39e-05;9.99e-01], 1:0.059s, tol=2.23e-02 79.75% -> 4.78s 321: 16.050µm, [2.39e-05;9.99e-01], 1:0.058s, tol=2.20e-02 80.00% -> 4.65s 322: 16.100µm, [2.40e-05;9.99e-01], 1:0.059s, tol=2.28e-02 80.25% -> 4.65s 323: 16.150µm, [2.40e-05;9.99e-01], 1:0.061s, tol=2.10e-02 80.50% -> 4.73s 324: 16.200µm, [2.40e-05;9.99e-01], 1:0.059s, tol=2.09e-02 80.75% -> 4.57s 325: 16.250µm, [2.40e-05;9.99e-01], 1:0.060s, tol=2.12e-02 81.00% -> 4.60s 326: 16.300µm, [2.40e-05;9.99e-01], 1:0.059s, tol=1.93e-02 81.25% -> 4.43s 327: 16.350µm, [2.40e-05;9.99e-01], 1:0.059s, tol=1.93e-02 81.50% -> 4.34s 328: 16.400µm, [2.41e-05;9.99e-01], 1:0.058s, tol=1.98e-02 81.75% -> 4.24s 329: 16.450µm, [2.41e-05;9.99e-01], 1:0.059s, tol=1.81e-02 82.00% -> 4.28s 330: 16.500µm, [2.41e-05;9.99e-01], 1:0.059s, tol=1.66e-02 82.25% -> 4.17s 331: 16.550µm, [2.41e-05;9.99e-01], 1:0.058s, tol=1.73e-02 82.50% -> 4.08s 332: 16.600µm, [2.41e-05;9.99e-01], 1:0.058s, tol=1.72e-02 82.75% -> 4.03s 333: 16.650µm, [2.41e-05;9.99e-01], 1:0.058s, tol=1.66e-02 83.00% -> 3.96s 334: 16.700µm, [2.41e-05;9.99e-01], 1:0.059s, tol=1.63e-02 83.25% -> 3.94s 335: 16.750µm, [2.41e-05;9.99e-01], 1:0.058s, tol=1.71e-02 83.50% -> 3.86s 336: 16.800µm, [2.41e-05;9.99e-01], 1:0.059s, tol=1.88e-02 83.75% -> 3.83s 337: 16.850µm, [2.41e-05;9.99e-01], 1:0.058s, tol=2.21e-02 84.00% -> 3.74s 338: 16.900µm, [2.42e-05;9.99e-01], 1:0.058s, tol=2.61e-02 84.25% -> 3.66s 339: 16.950µm, [2.42e-05;9.99e-01], 1:0.059s, tol=3.05e-02 84.50% -> 3.65s 340: 17.000µm, [2.42e-05;1.00e+00], 1:0.062s, tol=3.18e-02 84.75% -> 3.80s 341: 17.050µm, [2.42e-05;1.00e+00], 1:0.060s, tol=3.25e-02 85.00% -> 3.57s 342: 17.100µm, [2.42e-05;1.00e+00], 1:0.059s, tol=3.29e-02 85.25% -> 3.46s 343: 17.150µm, [2.43e-05;1.00e+00], 1:0.058s, tol=2.48e-02 85.50% -> 3.39s 344: 17.200µm, [2.43e-05;1.01e+00], 1:0.059s, tol=1.33e-02 85.75% -> 3.35s 345: 17.250µm, [2.43e-05;1.01e+00], 1:0.059s, tol=8.31e-03 86.00% -> 3.29s 346: 17.300µm, [2.43e-05;1.01e+00], 1:0.059s, tol=6.99e-03 86.25% -> 3.23s 347: 17.350µm, [2.43e-05;1.01e+00], 1:0.060s, tol=6.52e-03 86.50% -> 3.24s 348: 17.400µm, [2.43e-05;1.01e+00], 1:0.060s, tol=6.25e-03 86.75% -> 3.20s 349: 17.450µm, [2.43e-05;1.01e+00], 1:0.060s, tol=6.08e-03 87.00% -> 3.10s 350: 17.500µm, [2.43e-05;1.01e+00], 1:0.061s, tol=5.96e-03 87.25% -> 3.12s 351: 17.550µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.88e-03 87.50% -> 2.97s 352: 17.600µm, [2.43e-05;1.01e+00], 1:0.076s, tol=5.83e-03 87.75% -> 3.75s 353: 17.650µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.79e-03 88.00% -> 2.82s 354: 17.700µm, [2.43e-05;1.01e+00], 1:0.058s, tol=5.76e-03 88.25% -> 2.74s 355: 17.750µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.74e-03 88.50% -> 2.74s 356: 17.800µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.71e-03 88.75% -> 2.65s 357: 17.850µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.69e-03 89.00% -> 2.61s 358: 17.900µm, [2.43e-05;1.01e+00], 1:0.058s, tol=5.67e-03 89.25% -> 2.51s 359: 17.950µm, [2.43e-05;1.01e+00], 1:0.058s, tol=5.65e-03 89.50% -> 2.44s 360: 18.000µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.64e-03 89.75% -> 2.41s 361: 18.050µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.62e-03 90.00% -> 2.37s 362: 18.100µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.60e-03 90.25% -> 2.29s 363: 18.150µm, [2.43e-05;1.01e+00], 1:0.058s, tol=5.59e-03 90.50% -> 2.20s 364: 18.200µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.58e-03 90.75% -> 2.17s 365: 18.250µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.57e-03 91.00% -> 2.13s 366: 18.300µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.55e-03 91.25% -> 2.07s 367: 18.350µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.54e-03 91.50% -> 2.00s 368: 18.400µm, [2.43e-05;1.01e+00], 1:0.059s, tol=5.52e-03 91.75% -> 1.94s Saved mesh data in: /home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/v1.9.0/examples/PhaseField/results/Shear2D/Test/Miehe_AT1_optimMesh/Meshes Saved simulation and summary in: /home/docs/checkouts/readthedocs.org/user_builds/easyfea/checkouts/v1.9.0/examples/PhaseField/results/Shear2D/Test/Miehe_AT1_optimMesh Generate movie 01/21 (4.76 %) 2.95 s Generate movie 02/21 (9.52 %) 2.55 s Generate movie 03/21 (14.29 %) 2.45 s Generate movie 04/21 (19.05 %) 2.28 s Generate movie 05/21 (23.81 %) 2.15 s Generate movie 06/21 (28.57 %) 2.02 s Generate movie 07/21 (33.33 %) 1.91 s Generate movie 08/21 (38.10 %) 1.75 s Generate movie 09/21 (42.86 %) 1.61 s Generate movie 10/21 (47.62 %) 1.48 s Generate movie 11/21 (52.38 %) 1.35 s Generate movie 12/21 (57.14 %) 1.21 s Generate movie 13/21 (61.90 %) 1.07 s Generate movie 14/21 (66.67 %) 934.34 ms Generate movie 15/21 (71.43 %) 803.96 ms Generate movie 16/21 (76.19 %) 668.08 ms Generate movie 17/21 (80.95 %) 540.04 ms Generate movie 18/21 (85.71 %) 409.39 ms Generate movie 19/21 (90.48 %) 266.78 ms Generate movie 20/21 (95.24 %) 135.12 ms Generate movie 21/21 (100.00 %) 0.00 µs Mesh: 55.255 ms Boundary Conditions: 7.780 ms Matrix: 9.193 s Split: 2.008 s Solver: 13.626 s PostProcessing: 397.647 ms Resolution phase field: 21.612 s Display: 598.827 ms PyVista_Interface: 4.457 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""" 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, folder=folder_save) 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 29.091 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 `_