We present SU(3) gluon propagators calculated on 48*48*48*N_t lattices at beta=6.8 where N_t=64 (corresponding the confinement phase) and N_t=16 (deconfinement) with the bare gauge parameter,alpha, set to be 0.1. In order to avoid Gribov copies, we employ the stochastic gauge fixing algorithm. Gluon propagators show quite different behavior from those of massless gauge fields: (1) In the confinement phase, G(t) shows massless behavior at small and large t, while around 5<t<15 it behaves as massive particle, and (2) effective mass observed in G(z) becomes larger as z increases. (3) In the deconfinement phase, G(z) shows also massive behavior but effective mass is less than in the confinement case. In all cases, slope masses are increasing functions of t or z, which can not be understood as addtional physical poles.