We present long-time simulations of expanding ultracold neutral plasmas, including a full treatment of the strongly coupled ion dynamics. Thereby, the relaxation of the expanding laser-cooled plasma is studied, taking into account elastic as well as inelastic collisions. It is demonstrated that, depending on the initial conditions, the ionic component of the plasma may exhibit short-range order or even a superimposed long-range order resulting in concentric ion shells. In contrast to ionic plasmas confined in traps, the shell structures build up from the center of the plasma cloud rather than from the periphery.