A numerical model, EC2, is developed to simulate the consolidation of a soil layer subjected to combined electroosmosis and surcharge preloading. EC2 includes the capabilities of a previous model, EC1, in accounting for large-strain-induced nonlinear changes of the physical and geoelectrical properties occurring to the layer, with the additions of two-dimensional consolidation, coupled loading and electric fields, time-dependent loading and current density, and an external hydraulic gradient. A highlight of the EC2 model is the superposition of multiple streams between small elements, which improves the prediction accuracy of the rate of flow and the degree of consolidation. The consolidation algorithm of EC2 is one-dimensional and conducted using finite difference and piecewiselinear methods. This paper presents the theoretical and numerical development of the EC2 model. A companion paper presents the validation of the model and the results of simulations that illustrate the optimization of the consolidating processes for some interesting design scenarios.