The fractional Zener constitutive law is frequently used as a model of solid-like viscoelastic behavior. In this work, a class of linear viscoelastic models of Zener type, which generalize the fractional Zener model, is studied by the use of Bernstein functions technique. We prove that the corresponding relaxation moduli are completely monotone functions under appropriate thermodynamic restrictions on the parameters. Based on this property, we study the propagation function and establish the subordination principle for the corresponding Zener-type wave equation, which provides an integral representation of the solution in terms of the propagation function and the solution of a related classical wave equation. The analytical findings are supported by numerical examples.