The conventional method for functional quantile regression (FQR) is to fit the regression model for each quantile of interest separately. Therefore, the slope function of the regression, as a bivariate function of time and quantile, is estimated as a univariate function of time for each fixed quantile. There are several limitations to this conventional strategy. For example, the monotonicity of conditional quantiles can not be guaranteed, and the smoothness of the slope estimator as a bivariate function can not be controlled. In this paper, we develop a new framework for functional quantile regression. We propose to simultaneously fit the functional quantile regression model for multiple quantiles with the help of a bivariate basis under some constraints so that the estimated quantiles satisfy the monotonicity conditions. Meanwhile, the smoothness of the slope estimator is controlled.The proposed estimator for the slope function is shown to be asymptotically consistent. In addition, we also establish its asymptotic normality. Simulation studies are implemented to evaluate the finite sample performance of the proposed method in comparison with the