We report preliminary results on the analysis of the three-body Y(10860)=> B barB pi, Y(10860) => (B barB* + c.c.) pi and Y(10860)=> B* barB* pi decays including an observation of the Y(10860)=> Zb(10610)+- pi-+ => [B barB*+c.c]+- pi-+ and Y(10860)=> Zb(10650)+- pi-+ => [B* barB*]+- pi-+ decays as intermediate channels. We measure branching fractions of the three-body decays to be Br(Y(10860)=> [B barB*+c.c.]+- pi-+)=(28.3+-2.9+-4.6)x10^{-3} and Br(Y(10860)=> [B* barB*]+- pi-+)=(14.1+-1.9+-2.4)x10^{-3} and set 90% C.L. upper limit Br(Y(10860)=> [B barB]+- pi-+)<4.0x10^{-3}. We also report results on the amplitude analysis of the three-body Y(10860)=>Y(nS)pi+pi-, n=1,2,3 decays and the analysis of the internal structure of the three-body Y(10860)=>hb(mP)pi+pi-, m=1,2 decays. The results are based on a 121.4 1/fb data sample collected with the Belle detector at a center-of-mass energy near the Y(10860).