second ordered coupled time delayed differential equations.

here,i varies from 1 to 3. Values of N_i, ma, me,ca,ce, phi(a), phi(b), tau and g are known and constants. How can I solve this system? Is solving with 'dde' feasible?