Because of the inherent nonlinearities involving the behavior of CNTs when excited by electrostatic forces, modeling and simulating their behavior is challenging. The complicated form of the electrostatic force describing the interaction of their cylindrical shape, forming upper electrodes, to lower electrodes poises serious computational challenges. This presents an obstacle against applying and using several nonlinear dynamics tools typically used to analyze the behavior of complicated nonlinear systems undergoing large motion, such as shooting, continuation, and integrity analysis techniques. This works presents an attempt to resolve this issue. We present an investigation of the nonlinear dynamics of carbon nanotubes when actuated by large electrostatic forces. We study expanding the complicated form of the electrostatic force into enough number of terms of the Taylor series. Then, we utilize this form along with an Euler-Bernoulli beam model to study for the first time the dynamic behavior of CNTs when excited by large electrostatic force. The geometric nonlinearity and the nonlinear electrostatic force are considered. An efficient reduced-order model (ROM) based on the Galerkin method is developed and utilized to simulate the static and dynamic responses of the CNTs. Several results are generated demonstrating softening and hardening behavior of the CNTs near their primary and secondary resonances. The effects of the DC and AC voltage loads on the behavior have been studied. The impacts of the initial slack level and CNT diameter are also demonstrated.