Systems, methods, and software can be used to simulate a fracture treatment. In some aspects, a common solution vector for multiple distinct subsystem models is defined. Each subsystem model represents a distinct subsystem of an injection treatment system. Parameters of the subsystem models are updated based on the solution vector according to predefined relationships between the solution vector and the parameters of the subsystem models. The subsystem models are operated based on the solution vector and the updated parameters.