The phase-field method enables simulating the spatiotemporal evolution of the coupled physical-order parameters under externally applied fields in a wide range of materials and devices. Leveraging advanced numerical algorithms for solving the nonlinear partial differential equations and scalable parallelization techniques, the phase-field method is becoming a powerful computational tool to model and design devices operating based on multiple-coupled physical processes. This article will highlight examples of applying phase-field simulations to predict new mesoscale physical phenomena and design new-concept magnetomechanical devices by identifying the desirable combination of the composition, size, and geometry of monolithic materials as well as the device structure. A brief outlook of the opportunities and challenges for modeling and designing magnetomechanical devices with phase-field modeling is also provided.