The spin-1 Ising model with the nearest-neighbour bilinear and biquadratic interactions and single-ion anisotropy is simulated on a cellular automaton which improved from the Creutz cellular automaton (CCA) for a simple cubic lattice. The simulations have been made for several k=K/J and d=D/J in the 0 <= d < 3 and -2 <= k <= 0 parameter regions. We confirm the existence of the re-entrant and the successive re-entrant phase transitions near the phase boundary. The phase diagrams characterizing phase transitions are presented for comparison with those obtained from other calculations. The static critical exponents are estimated within the framework of the finite-size scaling theory at d = 0, 1 and 2 in the interval -2 <= k <= 0. The results are compatible with the universal Ising critical behavior.