We develop a novel abundance matching method to construct a mock catalog of luminous red galaxies (LRGs) in SDSS, using catalogs of halos and subhalos in N-body simulations for a LCDM model. Motivated by observations suggesting that LRGs are passively-evolving, massive early-type galaxies with a typical age >5Gyr, we assume that simulated halos at z=2 (z2-halo) are progenitors for LRG-host subhalos observed today, and we label the most tightly bound particles in each progenitor z2-halo as LRG ``stars. We then identify the subhalos containing these stars to z=0.3 (SDSS redshift) in descending order of the masses of z2-halos until the comoving number density of the matched subhalos becomes comparable to the measured number density of SDSS LRGs, n=10^{-4} (h/Mpc)^3. Once the above prescription is determined, our only free parameter is the number density of halos identified at z=2 and this parameter is fixed to match the observed number density at z = 0.3. By tracing subsequent merging and assembly histories of each progenitor z2-halo, we can directly compute, from the mock catalog, the distributions of central and satellite LRGs and their internal motions in each host halo at z=0.3. While the SDSS LRGs are galaxies selected by the magnitude and color cuts from the SDSS images and are not necessarily a stellar-mass-selected sample, our mock catalog reproduces a host of SDSS measurements: the halo occupation distribution for central and satellite LRGs, the projected auto-correlation function of LRGs, the cross-correlation of LRGs with shapes of background galaxies (LRG-galaxy weak lensing), and the nonlinear redshift-space distortion effect, the Finger-of-God effect, in the angle-averaged redshift-space power spectrum.