The effect of convection in RR Lyrae stars has been investigated using nonlinear models that include the effects of time dependence, turbulent pressure, convective overshooting, and convection-pulsation interaction. The convective structure reduces to that of mixing-length theory in the limit of no time dependence and thick convective zones. Initial tests suggested that convection has an important effect on the stability and pulsation mode of these stars. We now have investigated the nonlinear nature of convection at limiting amplitude. Convection serves as a limiting amplitude mechanism for these stars by causing turbulence near the phase of minimum radius that damps the driving. A comparison with Geneva observations of RR Lyrae shows good agreement in the phase dependence and amplitude of the turbulent motions.