Singular anisotropic elliptic equations with gradient-dependent lower order terms

We prove the existence of a solution to a singular anisotropic elliptic equation in a bounded open subset $\Omega$ of $\mathbb R^N$ with $N\ge 2$, subject to a homogeneous boundary condition: \begin{equation} \label{eq0} \left\{ \begin{array}{ll} \mathcal A u+ \Phi(u,\nabla u)=\Psi(u,\nabla u)+ \mathfrak{B} u \quad& \mbox{in } \Omega,\\ u=0 & \mbox{on } \partial\Omega. \end{array} \right. \end{equation} Here $ \mathcal A u=-\sum_{j=1}^N |\partial_j u|^{p_j-2}\partial_j u$ is the anisotropic $\overrightarrow{p}$-Laplace operator, while $\mathfrak B$ is an operator from $W_0^{1,\overrightarrow{p}}(\Omega)$ into $W^{-1,\overrightarrow{p}'}(\Omega)$ satisfying suitable, but general, structural assumptions. $\Phi$ and $\Psi$ are gradient-dependent nonlinearities whose models are the following: \begin{equation*} \label{phi}\Phi(u,\nabla u):=\left(\sum_{j=1}^N \mathfrak{a}_j |\partial_j u|^{p_j}+1\right)|u|^{m-2}u, \quad \Psi(u,\nabla u):=\frac{1}{u}\sum_{j=1}^N |u|^{\theta_j} |\partial_j u|^{q_j}. \end{equation*} We suppose throughout that, for every $1\leq j\leq N$, \begin{equation*}\label{ass} \mathfrak{a}_j\geq 0, \quad \theta_j>0, \quad 0\leq q_j更多