Global solutions of quasilinear systems of Klein--Gordon equations in 3D