diff --git a/web/content/docs/benchmarks/creep-after-excavation-bgra/CreepAfterExcavation.md b/web/content/docs/benchmarks/creep-after-excavation-bgra/CreepAfterExcavation.md
new file mode 100644
index 0000000000000000000000000000000000000000..5b72eacec498bccc4c3811ffccd308a415118ca6
--- /dev/null
+++ b/web/content/docs/benchmarks/creep-after-excavation-bgra/CreepAfterExcavation.md
@@ -0,0 +1,127 @@
++++
+date = "2018-08-10T11:46:45+01:00"
+title = "Creep in a rectangle domain with hole under pressure on its top surface"
+weight = 171
+project = "ThermoMechanics/CreepBGRa/CreepAfterExcavation/CreepAfterExcavation.prj"
+author = "Wenqing Wang"
+
+[menu]
+ [menu.benchmarks]
+ parent = "thermo-mechanics"
+
++++
+
+{{< data-link >}}
+
+-------------------------------------------------------------------------
+Following up the benchmark about the BGRa creep model described on the
+[this page](https://www.opengeosys.org/docs/benchmarks/creepbgra/creepbrga/),
+this example represents the creep in the near field of
+drift in the deep rock salt after excavation. The domain and the
+geometry are shown in the following figure:
+<figure>
+ <img src="../mesh.png" alt="Mesh and Geometry" id="fig_6"
+ style="height:400px;width:490px;">
+ <figcaption>Mesh and Geometry.</figcaption>
+</figure>
+
+The domain has two material groups, which are highlighted by different
+colors. The material group that is in the
+top part of the domain represents a cap rock type, while the other
+material group is for rock salt. The material properties of the rocks are given in the following table:
+
+---------------------- ---------- ----------- ------------
+| |Cap rock | Rock salt| Unit|
+| --- |:---------:| -----:|----:|
+|Density | 2000 | 2170 | kg/m$^{3}$|
+|Young’s Modulus | 7.0 | 7.65 | GPa |
+|Poisson ratio | 0.3 | 2.7 | - |
+|Thermal conductivity | 5 | 5 | W/(mK) |
+---------------------- ---------- ----------- ------------
+
+The parameters of the BGRa creep model are $A=0.18\, \mbox{d}^{-1}$,
+$m=5$, $Q=54 \mbox{ kJ/mol}$.
+
+
+The width
+and the height of of the domain are 300 m and 340 m, respectively. The
+height of the cap rock portion is 40 m. The drift to be excavated has a
+radius of 50 m.
+
+Here we consider the creep of the rock after excavation. Therefore, we
+assume a steady state of stress and temperature of excavation at the
+beginning of the current simulation. For this assumption, the boundary
+conditions are given as:
+
+- top boundary: $\tau_x = \sigma_x=0$, $\tau_y=\sigma_y=70$ MPa,
+ $T=310$ K.
+
+- two lateral boundaries: normal displacement is fixed, and no heat
+ flux.
+
+- bottom boundary: normal displacement is fixed, and $T=320$ K.
+
+- circle of drift surface: traction free
+ (${ \mathbf\sigma}\cdot \mathbf n = 0$), and $T=300$ K as excavation
+ conditions.
+
+The boundary conditions of the mechanical process lead to a distinct
+stress field after excavation. Besides, heat capacity is neglected in
+order to pose a steady state temperature field after excavation. The
+initial stresses are zero and the initial temperature is linearly
+distributed from top to bottom with the boundary values of temperature.
+
+The time step sizes of the simulation are: One step of 0.001 day, 10 steps of 0.1 day, and the
+remaining steps of 1 day.
+
+The following three figures are plotted by using the results of the simulation of 1000 days creep.
+The three figures display the distribution of horizontal and vertical stresses at times of
+108 days, 409 days and 1000 days, respectively.
+
+In theses three figures,
+the left sub-figure show the time variations of horizontal and vertical stresses
+ at a position just close to the top of the drift, and the
+green vertical line in it marks the time of the displayed stress field.
+
+<figure>
+ <img src="../stress_xx_yy_20.png" alt="Stress distribution at the time of 109 days." id="fig_2">
+ <figcaption>Stress distribution at the time of 109 days.</figcaption>
+</figure>
+
+<figure>
+ <img src="../stress_xx_yy_50.png" alt="Stress distribution at the time of 409 days." id="fig_3">
+ <figcaption>Stress distribution at the time of 409 days.</figcaption>
+</figure>
+
+<figure>
+ <img src="../stress_xx_yy_110.png" alt="Stress distribution at the time of 1000 days." id="fig_4">
+ <figcaption>Stress distribution at the time of 1000 days.</figcaption>
+</figure>
+
+
+The above three figures show that the absolute value of the horizontal stress increase
+signification in the areas above and beneath the drift due to the creep,
+and the change of the vertical stress is slow compared to that of the
+horizontal stress.
+
+The following figure shows the strain distribution at the
+end of the simulation time at 1000 days.
+<figure>
+ <img src="../strain.png" alt="Strain distribution at the time of 1000 days." id="fig_5">
+ <figcaption>Strain distribution at the time of 1000 days.</figcaption>
+</figure>
+
+
+The steady-state temperature distribution is displayed in the following figure
+<figure>
+ <img src="../T.png" alt="Temperature distribution at the time of 1000 days." id="fig_6">
+ <figcaption>Temperature distribution at the time of 1000 days.</figcaption>
+</figure>
+
+## Note:
+For the automatic benchmarking, the time duration of creep is reduced to 50 days in order to reduce the run time.
+ If one wants to test this benchmark for 1000 days' creep, please change the end time in the tag of `<t_end> `
+in the project file as
+ ```
+<t_end>1000</t_end>
+```
\ No newline at end of file
diff --git a/web/content/docs/benchmarks/creep-after-excavation-bgra/T.png b/web/content/docs/benchmarks/creep-after-excavation-bgra/T.png
new file mode 100644
index 0000000000000000000000000000000000000000..3ad5f8b981098fb74e5dfda508f009763501d546
--- /dev/null
+++ b/web/content/docs/benchmarks/creep-after-excavation-bgra/T.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:a360404be4d5764b85c0bead5ba92d92faaa9542a12957f72d9d706b4b4e39c6
+size 384898
diff --git a/web/content/docs/benchmarks/creep-after-excavation-bgra/mesh.png b/web/content/docs/benchmarks/creep-after-excavation-bgra/mesh.png
new file mode 100644
index 0000000000000000000000000000000000000000..9b4029da7e61d2fac8e848b5ece92b690ec3a93e
--- /dev/null
+++ b/web/content/docs/benchmarks/creep-after-excavation-bgra/mesh.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c2095fc571a43c62a3eb1816d8aaf0ead9f4c66f1974a544fda5ced0d2beed49
+size 326201
diff --git a/web/content/docs/benchmarks/creep-after-excavation-bgra/strain.png b/web/content/docs/benchmarks/creep-after-excavation-bgra/strain.png
new file mode 100644
index 0000000000000000000000000000000000000000..f0e960169e3b6baa242ca08e5f069fd020fd12b5
--- /dev/null
+++ b/web/content/docs/benchmarks/creep-after-excavation-bgra/strain.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:ff2fb210d6f7d543082f490c55ab265e27095efef30fe60f6dd596df0d87ce67
+size 153593
diff --git a/web/content/docs/benchmarks/creep-after-excavation-bgra/stress_xx_yy_110.png b/web/content/docs/benchmarks/creep-after-excavation-bgra/stress_xx_yy_110.png
new file mode 100644
index 0000000000000000000000000000000000000000..d192fbc4e577658358db571ff8bd1aba005aaff3
--- /dev/null
+++ b/web/content/docs/benchmarks/creep-after-excavation-bgra/stress_xx_yy_110.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:abde5d089d6edc8847dc425c8feec31515832449415bc6bb5f8a414d05dd632a
+size 870467
diff --git a/web/content/docs/benchmarks/creep-after-excavation-bgra/stress_xx_yy_20.png b/web/content/docs/benchmarks/creep-after-excavation-bgra/stress_xx_yy_20.png
new file mode 100644
index 0000000000000000000000000000000000000000..0323c62ccc22eae7c819e42ac8e5a474d32b2900
--- /dev/null
+++ b/web/content/docs/benchmarks/creep-after-excavation-bgra/stress_xx_yy_20.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:3ea312ed37bd6c593ca4b4a7f4a7b1d7d064c83af53bc70c675eadd2517adef2
+size 847877
diff --git a/web/content/docs/benchmarks/creep-after-excavation-bgra/stress_xx_yy_50.png b/web/content/docs/benchmarks/creep-after-excavation-bgra/stress_xx_yy_50.png
new file mode 100644
index 0000000000000000000000000000000000000000..e3cc98ff1c410a7007d97db00dda36ff8b16a6eb
--- /dev/null
+++ b/web/content/docs/benchmarks/creep-after-excavation-bgra/stress_xx_yy_50.png
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:f50ecaaeaaf49a8e6ac371571f55ab0278df5b27a15b43c015ebdf7827b93c32
+size 858276