In this paper, we discuss univariate and bivariate constrained interpolation problems using the \(\beta \) -rational cubic fractal interpolation function ( \(\beta \) -RCFIF) based solely on function values. The \(\beta \) -RCFIF takes the form \(\tfrac{P_i(\textbf{u})}{Q_i(\textbf{u})}\) , where \(P_i(\textbf{u})\) is a cubic polynomial determined by the interpolation conditions and \(Q_i(\textbf{u})\) is a preassigned quadratic polynomial. We establish suitable conditions on the associated iterated function system (IFS) parameters to ensure that the graph of the corresponding \(\beta \) -RCFIF (the attractor) lies between two specified piecewise functions. For appropriate choices of tension parameters and vertical scaling factors, the \(\beta \) -RCFIF remains confined within a rectangle. Furthermore, we develop a bivariate cubic rational fractal model (BCRFM) by applying the Coons blending technique. By selecting appropriate shape parameters and scaling factors, the BCRFM is shown to lie above a prescribed plane. With properly chosen IFS parameters, the BCRFM is also contained within a cuboid. Several numerical examples are provided to illustrate the effectiveness of the proposed interpolation scheme. Additionally, we incorporate ARIMA-based forecasting to complement the fractal interpolation framework.