Numerical models of the wind-blown bubble of massive stars usually only account for the wind of a single star. However, since massive stars are usually formed in clusters, it would be more realistic to follow the evolution of a bubble created by several stars. We develope a two-dimensional (2D) model of the circumstellar bubble created by two massive stars, a 40 solar mass star and a 25 solar mass star, and follow its evolution. The stars are separated by approximately 16 pc and surrounded by a cold medium with a density of 20 particles per cubic cm. We use the MPI-AMRVAC hydrodynamics code to solve the conservation equations of hydrodynamics on a 2D cylindrical grid using time-dependent models for the wind parameters of the two stars. At the end of the stellar evolution (4.5 and 7.0 million years for the 40 and 25 solar mass stars, respectively), we simulate the supernova explosion of each star. Each star initially creates its own bubble. However, as the bubbles expand they merge, creating a combined, aspherical bubble. The combined bubble evolves over time, influenced by the stellar winds and supernova explosions. The evolution of a wind-blown bubble created by two stars deviates from that of the bubbles around single stars. In particular, once one of the stars has exploded, the bubble is too large for the wind of the remaining star to maintain and the outer shell starts to disintegrate. The lack of thermal pressure inside the bubble also changes the behavior of circumstellar features close to the remaining star. The supernovae are contained inside the bubble, which reflects part of the energy back into the circumstellar medium.