Our main purpose is to establish the existence of weak solutions of second order quasilinear elliptic systems ⎧ $-{\Delta }_{p}u+{\left|u\right|}^{p-2}u=f_{1\lambda ₁}{\left(x\right)\left|u\right|}^{q-2}u+2\alpha /(\alpha +\beta )g_{\mu }{\left|u\right|}^{\alpha -2}u{\left|v\right|}^{\beta }$, x ∈ Ω, ⎨ $-{\Delta }_{p}v+{\left|v\right|}^{p-2}v=f_{2\lambda ₂}{\left(x\right)\left|v\right|}^{q-2}v+2\beta /(\alpha +\beta )g_{\mu }{\left|u\right|}^{\alpha }{\left|v\right|}^{\beta -2}v$, x ∈ Ω, ⎩ u = v = 0, x∈ ∂Ω, where 1 < q < p < N and $\Omega \subset {ℝ}^N$ is an open bounded smooth domain. Here λ₁, λ₂, μ ≥ 0 and $f_{i{\lambda }_i}\left(x\right)={\lambda }_i{f}_{i+}\left(x\right)+f_{i-}\left(x\right)$ (i = 1,2) are sign-changing functions, where $f_{i\pm }\left(x\right)=max\pm f_i\left(x\right),0$, $g_{\mu }\left(x\right)=a\left(x\right)+\mu b\left(x\right)$, and ${\Delta }_{p}u={div\left(\right|\nabla u|}^{p-2}\nabla u)$ denotes the p-Laplace operator. We use variational methods.