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, 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 the static and dynamic behavior of CNTs. 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 effect of DC voltage load and AC voltage load on the nonlinearity has been studied. We also investigated the impact of initial slack level on the natural frequency and the nonlinearity. Small diameter and large initial slacked CNTs has been considered.