I'm trying to put together a set of trajectories for a polar lunar propellant tug; departure LLO @ 90 degrees inclination - TEI - HEEO @ X degrees inclination - circular equatorial - TLI - departure LLO (ideally the same orbit).

I can do all sorts of parts of it:

HHEO@X - Equatorial: easy

Equatorial - LLO@90 by B plane targeting: easy - I think I can hit "clockwise" or anti but haven't tried both,

HEEO - LLO@90: also easy (and indicates there would be an energetically similar trajectory going the other way?)

which is good considering the low base I'm coming off...

But I'm having trouble trying to find a LLO@90 - HEEO - basically getting a low Earth Pe from the TEI, the solver just seems to get stuck jittering between two bad solutions and not really exploring the parameter space (TA for burn, delta-v). Also conceptually how can one find which turn of the corkscrew that is LLO is "the one"?

I've got a bad feeling I'll have to do something like your multi shooting example with two segments? (and repeat for each turn of the screw) a forward TEI and a backward HEEO capture burn and optimize to join them together? Or can I back prop solving for (epoch, deltaV)@(Pe on line of nodes) achieving Luna Pe "at the right place" like reverse B plane targeting?

Thanks - I hope someone can help...