One approach is to start from $T=\frac12 m[\dot{x}^2+\dot{y}^2+\dot{z}^2]$ and use $x=\rho\cos\theta$, $y=\rho\sin\theta$ so that, e.g. $\dot{x}=\dot{\rho}\cos\theta-\rho\dot{\theta}\sin\theta$. Combine terms and simplify.
Or, for the plane polar part, using $\frac{d\hat{\boldsymbol{\rho}}}{dt} =\dot{\theta}\hat{\boldsymbol{\theta}} $ (consider the change in $\hat{\boldsymbol{\rho}}$ with small $dt$ and compare to the change in $\theta$), we have $\dot{\boldsymbol{\rho}}=\dot{\rho}\hat{\boldsymbol{\rho}}+\rho \frac{d\hat{\boldsymbol{\rho}}}{dt} =\dot{\rho}\hat{\boldsymbol{\rho}}+\rho \dot{\theta}\hat{\boldsymbol{\theta}} $ and use the orthogonality of $\hat{\boldsymbol{\rho}}$ and $\hat{\boldsymbol{\theta}}$ after squaring.
An alternative is to use $T=\frac12 m \left(\frac{ds}{dt}\right)^2 = \frac12 m g_{ij}\dot{q}^i \dot{q}^j$ where the metric $q_{ij}$ is diagonal with $q_{\rho\rho}=q_{zz}=1$ and $q_{\theta\theta}=\rho^2$.