When applied to simple models, large N expansion is not very tricky, however I would like to know if it is possible to automatize this process to don't expend much time in such calculations.
The only thing that I am aware of, and I will stress is that I have never used it, is the FeynCalc package. I don't know if you can use those tools to build routines that you need. I suppose that you might you already know of this package.