Techniques for the simulation of the tanpura have advanced significantly in recent years allowing numerically stable inclusion of bridge contact. In this paper tension modulation is added to a tanpura model containing a stiff lossy string, distributed bridge contact and the thread. The model is proven to be unconditionally stable and the numerical solver used has a unique solution as a result of choices made in the discretisation process. Effects due to the distribution of the bridge contact forces by comparison to a single point bridge and of introducing the tension modulation are studied in simulations. This model is intended for use in furthering the understanding of the physics of the tanpura and for informing the development of algorithms for sound synthesis of the tanpura and similar stringed instruments.