Trajectory piecewise-linearization-based reduced-order macromodeling methods have been proposed to characterize the time-domain behaviors of large strongly nonlinear systems. However, all these methods rely on frequency-domain model-order-reduction (MOR) methods for linear systems. Therefore, the ac