Abstract: Based on the relativistic mean field (RMF) model with Thomas-Fermi approximation, we investigate the elastic properties of neutron star matter. The elastic constants are estimated by introducing deformations on the nuclear pasta structures in β-equilibrium, where various crystalline configurations are considered in a fully three-dimensional geometry without the Wigner-Seitz approximation. Two scenarios with different symmetry energy slope (L=41.34 and 89.39 MeV) are examined, where the the elastic constants can vary by ten times. By fitting to the numerical results, we improve the analytic formulae for the elastic properties of nuclear pasta by introducing damping factors.