In this paper a collaborative simulation platform for proton exchange membrane fuel cell (PEMFC) power systems is presented, where the stack is simulated by a two-phase distributed parameter model and the auxiliary units by lumped parameter models. By exchanging the dynamic data between the external load/auxiliary units and PEMFC stack, dynamic simulation of PEMFC stack has been carried out during the load changes for various states associated with different characteristic variables. The internal states of the stack can be observed due to variation of external load/auxiliary units. Numerical experiments are provided for a special case with multiple cycles of load changes derived from an acceleration mode of a fuel cell vehicle. The numerical results demonstrate that the "undershoot" of output voltage is due to the response lag of the auxiliary units and liquid water accumulation in the fuel cell stack. (C) 2010 Professor T. Nejat Veziroglu. Published by Elsevier Ltd. All rights reserved.