Variable-order differential operators can be employed as a powerful tool to modeling nonlinear fractional differential equations and chaotical systems. In this paper, we propose a new generalize numerical schemes for simulating variable-order fractional differential operators with power-law, exponential-law and Mittag-Leffler kernel. The numerical schemes are based on the fundamental theorem of fractional calculus and the Lagrange polynomial interpolation. These schemes were applied to simulate the chaotic financial system and memcapacitor-based circuit chaotic oscillator. Numerical examples are presented to show the applicability and efficiency of this novel method.