Discussion

The shear locking has been detected in both cases which include modification of the initial mesh based on the maximum and optimized values of the longitudinal order *p*. One can see that the standard *hp*-adaptation leads to the final error value below the desired admissible value of *γT*. The same refers to the *hp*-adaptation preceded by the modification based on the optimized value of the longitudinal order of approximation *p* = 4. In the case of the *hp*-adaptation performed after the modification of the mesh based on the fixed maximum value of *p* = 8, the admissible error value has not been reached, even though the error has been diminished. The reason is the discrete character of the possible *h*-adaptation. For the estimated error level, this adaption has not been performed.

**Figure 9.** A quarter of a bending-dominated plate—initial mesh.

**Figure 10.** A quarter of a bending-dominated plate—*hp*-adapted (final) mesh.

**Figure 11.** A quarter of a bending-dominated plate—after simple detection, no *hp*-adaptation.

**Figure 12.** A quarter of a bending-dominated plate—mesh after optimized modification.

**Figure 13.** A quarter of a bending-dominated plate—after optimization and *hp*-adaptation.

**Figure 14.** A quarter of a bending-dominated plate—convergence for three adaptation cases.

### 5.2.2. A Bending-Dominated Plate
