## 6 Conclusions  
Three types of the transition piezoelectric elements joining the basic piezoelectric elements are possible. The joined basic elements are the three-dimensional piezoelectric elements (or hierarchical symmetric-thickness piezoelectric elements) and the first-order symmetric-thickness piezoelectric elements. The first (classical) transition element guarantees continuity of the displacement and electric potential fields between the basic elements. It also guarantees that the assumption of no elongation of the normals to the mid-surface holds on the boundary between the transition elements and the first-order elements. The second (modified) transition element, additionally, guarantees the gradual change from the plane stress to three-dimensional stress state between the basic elements of the first-order and three-dimensional (or hierarchical) character, respectively. The third (enhanced) transition element, additionally, assures the gradual change of the assumption of no elongation of the normals to the mid-surface between the basic elements. The latter assumption holds on the boundary with the first-order element and is gradually switched off towards the boundary with the three-dimensional (or hierarchical) element. On this boundary, this assumption is completely off. \begin{table}
\begin{tabular}{c c c c c c c} \hline \hline
**Problem** & **Result** & & & **Model Type** & & \\
**Type** & **Quantity * & **MI** & **TR1** & **TR2** & **TR3** & **RM** \\ \hline \multirow{4}{*}{\
\begin{tabular}{c} bending- \\ dominated \\ \end{tabular} } & figure no. & Figures 20–22 & Figure 20 & Figure 21 & Figure 22 & Figures 20–22 \\  & dofs number \(N\)[1] & 13,068 & 11,912 & 11,912 & 11,912 & 8712 \\  & \(\log(E_{r}-E)\) [N\(\cdot\)m] & 2.922 & 2.983 & 2.879 & 2.803 & 0.253 \\  & \((E_{r}-E)/E_{r}\) [\%] & 0.0149 & 0.0162 & 0.0127 & 0.0116 & 0.00003 \\  & \(-\log(\delta_{8}/\delta_{7})/\log(N_{8}/N_{7})\) & 3.330 & 3.301 & 2.867 & 3.789 & 4.504 \\ \hline \multirow{4}{*}{
\begin{tabular}{c} bending- \\ dominated \\ shell \\ \end{tabular} } & figure no. & Figures 23–25 & Figure 23 & Figure 24 & Figure 25 & Figures 23–25 \\  & dofs number \(N\)[1] & 13,068 & 10,824 & 10,824 & 10,824 & 8712 \\  & \(\log(E_{r}-E)\) [N\(\cdot\)m] & 1.701 & 1.749 & 1.626 & 1.478 & 0.511 \\  & \((E_{r}-E)/E_{r}\) [\%] & 0.0046 & 0.0051 & 0.0039 & 0.0018 & 0.0003 \\  & \(-\log(\delta_{8}/\delta_{7})/\log(N_{8}/N_{7})\) & 3.411 & 3.393 & 3.741 & 4.429 & 3.788 \\ \hline \hline \end{tabular}
\end{table} 
Table 2: Error and convergence summary for three model problems. The gradual change of the stress state from the plane to three-dimensional is based on introduction of the blending functions that determine contributions of the basic models to the transition model. Implementation of this idea on the element level needs modifications of the strain vector, and the piezoelectric and dielectric constant matrices as well. The gradual change of the assumption of no elongation of the normals needs utilization of the blending functions, which play the role of gradually switching functions. On the element level, the nodal values of the complement (to unity) of the gradually switching functions appear in the penalty stiffness matrix responsible for the gradual change of the constraints. Comparing the stress distributions within the transition zones for three model problems, one can notice that high stress gradients appear on the boundary between the classical transition elements and the first-order elements. In the case of the enhanced transition elements, the gradual change of the stress state can be observed. In the case of the modified transition element, the stress distribution is intermediate with respect to the previous two distributions. It can be observed that the discrete models of the three analyzed model structures, with three types of the transition elements employed, produce the convergence curves of slightly higher convergence for the case of the modified elements in comparison to the classical elements.