We study numerically the dynamical evolution of models of stellar systems which contain two heavy centers of gravity and extended shells inside and around orbits of these binaries. Initially the binary components surrounded by numerous small-mass particles on shells with small initial velocities. It is shown that at some moment of time there is a collapse of these particles on to the heavier component of the binary. We have the strong interactions of bodies at this time. Frequently, the ‘gravitational slingshot’ effect occurs. Further in time, some part of particles which initially were outside the binary orbit escape from the system. Other particles are captured by binary components forming in the center a ‘dumb-bell’ bar. During evolution of our models different kinds of structures of systems often are very similar to the observed structures of galaxies: spiral and elliptical galaxies, interacting galaxies, different kinds of flows and jets. Totally systems are expanding. Formation of different kinds of structures of galaxies depends on motions of particles with respect to the line of apses of the nuclear binary: 1. in the case of almost orthogonal motion (the gravitational slingshot effect), these particles escape from systems and before an escape form open expanding spirals; 2. in the alternative case (small angles of motion of the particles with respect to this line), these particles are captured by the binary components; 3. in the intermediate cases, particles form different kinds of flows like jets, rings, close spirals etc.