This paper investigates a class of strongly nonlinear elliptic problems with singular terms within the framework of variable exponent Sobolev spaces. The problem is formulated as a differential inclusion \( \zeta (u) + A(u) + H(x,u,\nabla u) \ni \dfrac{f}{u^{\gamma }}, \) where the operator \( A \) is of Leray–Lions type acting between the variable exponent Sobolev space \( W_0^{1,p(\cdot )}(\Omega ) \) and its dual \( W^{-1,p'(\cdot )}(\Omega ) \) . The multivalued mapping \( \zeta \) is maximal monotone with \( 0 \in \zeta (0) \) , and the nonlinearity \( H(x, s, \xi ) \) exhibits natural growth with respect to \( \xi \) of order \( p(x) \) . The singular term involves \( f \in L^{\infty }(\Omega ) \) and \( \gamma > 0 \) . Under appropriate assumptions on the variable exponent \( p(\cdot ) \) , we establish the existence of weak solutions using a comprehensive approximation scheme. The proof employs truncation methods, Yosida approximations of the maximal monotone operator, careful energy estimates, and sophisticated convergence analysis in variable exponent spaces. Our results extend and generalize previous work in fixed exponent Sobolev spaces to the more flexible and physically relevant setting of variable exponents.