We present a class of first-order perturbative spacetime solutions within the framework of \(F(R) = R + \alpha R^2\) gravity, coupled with a phenomenological hybrid equation of state \(p = w\rho + \beta \rho ^2\) . Using a covariant (1+1+2) approach and a perturbative expansion in the parameter \(\alpha \) , we obtain metric solutions truncated at \(\mathcal {O}(\alpha )\) for static spherically symmetric configurations. These approximate solutions are designed to isolate and study the specific effects of quadratic gravitational corrections on spacetime geometry, separate from the complexities of realistic nuclear matter. The solutions reveal key theoretical features: (i) a characteristic deformation of the spacetime metric governed by the parameter \(\gamma = (w + 2\beta \rho _c)/(1 + w + \beta \rho _c)\) , (ii) internal curvature transitions when \(1-3w+6\beta \rho _c < 0\) , and (iii) modified stability criteria within this theoretical framework. Our work provides a mathematically tractable perturbative model for studying how strong-field F(R) gravity interacts with dense matter, offering closed-form expressions at first order that illuminate the functional dependence of spacetime geometry on model parameter. These perturbative results provide a new class of benchmark solutions for testing modified gravity theories in the strong-field regime, complementing numerical studies that use microphysical equations of state.