Rearranging transfer function for Bode diagram
One way to look at this circuit is to apply the fast analytical circuits techniques or FACTs formalized by Dr. Middlebrook some years ago. This 2nd-order filter is actually the cascading of two 1st-order filters that can separately be solved:



Considering an infinite impedance at the NINV input (the op-amp is considered perfect), the first circuit is a classical low-pass filter governed by the following transfer function: \$H_1(s)=\frac{1}{1+\frac{s}{\omega_{p1}}}\$ with \$\omega_{p1}=\frac{1}{R_3C_a}\$.

For the second filter, we can show that this structure obeys the following expression: \$H_2(s)=\frac{H_0+H^1s\tau_1}{1+s\tau_1}\$ In this expression, \$H_0\$ designates the gain determined for \$s=0\$ and obtained by open-circuiting capacitor \$C_b\$. If you compute the gain in this non-inverting configuration, you should find: \$H_0=\frac{R_2}{R_1}+1\$.

The time constant \$\tau_1\$ is obtained by reducing the excitation to 0 V (shorting the \$V_{in}\$ source) and finding the resistance "seen" from the connections of \$C_b\$. This is the \$R?\$ exercise in the below picture:



By inspection, the only resistance you see if the open-loop gain is infinite is \$R_2\$ leading to \$\tau_1=R_2C_b\$. It gives a second pole located at \$\omega_{p2}=\frac{1}{R_2C_b}\$.

The zero can be obtained by calculating a gain \$H^1\$ in which capacitor \$C_b\$ is set in its high-frequency state (a short circuit). In this case, we have \$H^1=1\$ and the numerator of \$H_2\$ can be written as: \$N(s)=H_0(1+\frac{H^1}{H_0}s\tau_1)\$ which if you develop and rearrange this gives you the second intermediate transfer function: \$H_2(s)=H_0\frac{1+\frac{s}{\omega_{z1}}}{1+\frac{s}{\omega_{p2}}}\$ with \$\omega_{z1}=\frac{1}{(R_2||R_1)C_b}\$ and \$\omega_{p2}=\frac{1}{R_2C_b}\$.

Finally, cascading the two transfer functions leads to the final expression you want, written in a so-called low-entropy format:

\$H(s)=H_0\frac{1+\frac{s}{\omega_{z1}}}{(1+\frac{s}{\omega_{p1}})(1+\frac{s}{\omega_{p2}})}\$ with 

\$\omega_{z1}=\frac{1}{(R_2||R_1)C_b}\$
\$\omega_{p1}=\frac{1}{R_3C_a}\$
\$\omega_{p2}=\frac{1}{R_2C_b}\$
\$H_0=\frac{R_2}{R_1}+1\$

We can check these calculations with a simple Mathcad sheet comparing the brute-force approach and the FACTs: results are identical.




The adopted approach can appear a bit obscure to some readers but it's not and the technique is well formalized nowadays. The book given in the introduction is a possibility to explore these FACTs as well as the book from Vorpérian and all the papers found on-line. Learning the FACTs in 2019 is certainly a good project!