The design of shell-infill structures has been a focal point in the topology optimization community due to their advantages in energy absorption characteristics, strength-to weight ratio and bucking resistance. This paper introduces a phase field-based topology optimization method for designing shell-infill structures. Interface-related issues can be easily addressed through the phase field function. A coupled topology optimization process is proposed to establish the connection between the shell and infill, facilitating the generation of optimized structures. The shell thickness, infill pattern and infill volume percentage, can be naturally controlled by different model parameters. Additionally, multiscale phase field topology optimization integrates the numerical homogenization method to evaluate the effective elasticity matrix of the microstructural infill. The approach is introduced for a uniform, periodical microstructure layout in the infill region, thereby achieving superior mechanical properties. Numerical results indicate the effectiveness of the proposed method in the design of both 2D and 3D shell-infill structures.