  We study the following quasilinear elliptic system for all $i=1,\cdots,m$
\begin{equation*} \label{}
  -div(\Phi'(|\nabla u_i|^2) \nabla u_i) = H_i(u) \quad \text{in} \ \
\mathbb{R}^n
  \end{equation*} where $u=(u_i)_{i=1}^m: \mathbb R^n\to \mathbb R^m$ and the
nonlinearity $ H_i(u) \in C^1(\mathbb R^m)\to \mathbb R$ is a general
nonlinearity. Several celebrated operators such as the prescribed mean
curvature, the Laplacian and the $p$-Laplacian operators fit in the above form,
for appropriate $\Phi$. We establish a Hamiltonian identity of the following
form for all $x_n\in\mathbb R$ \begin{equation*}\label{} \int_{\mathbb R^{n-1}}
\left(\sum_{i=1}^{m} \left[ \frac{1}{2} \Phi\left(|\nabla u_i|^2\right) -
\Phi'\left(|\nabla u_i|^2\right) |\partial_{x_n} u_i|^2 \right] - \tilde H(u)
\right) d x'\equiv C, \end{equation*} where $x=(x',x_n)\in\mathbb R^{n}$ and
$\tilde H$ is the antiderivative of $H=(H_i)_{i=1}^m$. This can be seen as a
counterpart of celebrated pointwise inequalities provided by Caffarelli,
Garofalo and Segala in \cite{cgs} and by Modica in \cite{m}. For the case of
system of equations, that is when $m\ge 2$, we show that as long as $$\alpha
\ge \alpha^*:=\inf_{s>0}\left\{\frac{2 s \Phi'(s)}{\Phi(s)}\right\}$$ the
function $I_\alpha(r):=\frac{1}{r^{n-\alpha}} \int_{B_r} \sum_{i=1}^{m}
\Phi(|\nabla u_i|^2) - 2\tilde H(u)$ is monotone nondecreasing in $r$. We call
this a weak monotonicity formula since for $m=1$ it is shown in \cite{cgs} that
$I_\alpha(r)$ is monotone when $\alpha\ge 1$, under certain conditions on
$\Phi$. We prove De Giorgi type results and Liouville theorems for $H$-monotone
and stable solutions in two and three dimensions when the system is symmetric.
