Esta es otra forma de verlo, quizás más física.
En primer lugar, con la forma habitual de ver la "fuerza" del péndulo, calculamos la fuerza restauradora en la dirección tangencial, $mg\sin\theta$ . Dejar $m=1$ obtenemos una ecuación diferencial para $\theta$ : $$\ddot{\theta} = -g\sin\theta$$ Utilizando la aproximación $\sin\theta\approx\theta$ válido para las pequeñas $\theta$ Obtenemos la ecuación del movimiento armónico simple (MSS), y un período que es independiente de la amplitud. Ahora $\sin\theta<\theta$ (excepto para $\theta=0$ ), por lo que la fuerza restauradora real es siempre un poco menor. Por lo tanto, esperamos que tome un poco más largo durante un periodo que con el SHM: cuando la varilla del péndulo pasa por el punto más bajo, la fuerza que intenta hacerla volver es un poco más débil que en el caso del péndulo lineal, por lo que tardará más en invertir su dirección.
Ahora comparemos con la derivación que describes. Yo la llamaría "basada en la energía". Observe en primer lugar que $g(1-\cos\theta)$ es la energía potencial del péndulo:
Así que $u$ es la fracción de la energía total $E$ que es cinético: $$Eu = KE$$ $$E(1-u)=PE$$ Como la energía cinética es $\frac{1}{2}(L\dot{\theta})^2$ (ajuste $m=1$ ), tenemos $dt/d\theta = L/\sqrt{2KE}$ , por lo que el periodo es $$T = \int \frac{Ld\theta}{\sqrt{2Eu}} = \int \frac{Ldu}{\sqrt{2Eu}}\frac{d\theta}{du}$$ Dejando de lado las constantes, esta es la fórmula que has dado. Para ver esto, nota que las cosas dentro de las dos últimas raíces cuadradas de tu fórmula son: $$1-u=\frac{g(1-\cos\theta)}{E}$$ y $$1-\frac{E(1-u)}{2g} = \frac{1+\cos\theta}{2}$$ y por eso su producto es $$\frac{g}{2E}(1-\cos\theta)(1+\cos\theta) = \frac{g}{2E}\sin^2\theta$$ Diferenciación de la fórmula $u=1-(g/E)(1-\cos\theta)$ da $du = -(g/E)\sin\theta\,d\theta$ por lo que las fórmulas son equivalentes.
Ahora, punto clave: ¿qué aproximación hacemos al omitir el último factor? Esencialmente, $\cos\theta\approx 1$ para los pequeños $\theta$ . Al igual que la aproximación $\sin\theta\approx\theta$ utilizado en la derivación de la "fuerza", se trata de una aproximación de primer orden. Dado que $\cos\theta<1$ el integrando del péndulo real será mayor que el del péndulo lineal, y el periodo será mayor.
Por cierto, la integral para el caso no lineal es una integral elíptica del primer tipo. Así que no se puede expresar en términos de funciones elementales (es decir, combinaciones finitas de algebraicas, trigonométricas, exponenciales y logarítmicas).