A liquid-fueled molten salt reactor (MSR) can reach a deep burnup based on online reprocessing and continuously refueling, which requires significantly different burnup calculation methods for MSRs compared with those for the traditional reactors. To address the unique burnup features and consider the fidelity of isotopic evolution in an MSR, a fuel depletion code ThorMCB is developed based on the OpenMC coupled with a specific depletion code, MODEC. Furthermore, to lower the computational cost of acquiring the equilibrium state through the time evolution step by step for an MSR, an equilibrium burnup calculation code ThorMCB-eq based on the OpenMC and MODEC is developed, which can obtain the equilibrium burnup efficiently. A single fuel lattice of MSR and an a Molten Salt Fast Reactor (MSFR) benchmark are applied for verifying the correctness of the ThorMCB and ThorMCB-eq codes. Compared with a neutron transport calculation code KENO-VI coupled with MODEC, the maximum deviation of the dominant heavy nuclides (HNs) at equilibrium state by ThorMCB is less than 10%, and that of the total mass of fission products (FPs) is less than 3%. For the MSFR benchmark, the neutronic parameters including temperature reactivity coefficient, the mass evolution of main HNs and FPs and breeding ratio (BR) from ThorMCB agree with the references. The equilibrium behavior can be quickly obtained with ThorMCB-eq, and the relative mass deviations of most nuclides keep around 2% in comparison with the results of step-by-step burnup evolution with ThorMCB. Furthermore, the same fuel contents and micro one-group cross sections at equilibrium are obtained with two different types of start-up fuels and a constant power density and fuel reprocessing scheme. In conclusion, the verified results indicate that ThorMCB and ThorMCB-eq can both provide reliable simulation for depletion evolution and equilibrium burnup for MSR fuel cycle.