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