Creating a temperature-dependent solar cell in LTSpice (using varistors)
Here is an attempt at a dynamic solution (everything below is some mock-up, tinkered values):



V1 is the reference voltage. E1 serves as a feedback error amplifier that ultimately acts as a source for the load, too. At its negative input, is what would be the variation of the temperature as a function of the load current: Vx is the current sensor, F1 takes this, attenuates it by some value (here 100x), and then feeds it to a simple lowpass RC (can add a ladder afterwards for finer control), whose values make up the thermal constant you're interested in.

Note that I used Rpar=20 and ic=25. These two values represent the temperature at "rest" (I(F1)*Rpar) and the initial temperature (ic). This is what comes up:



If you look at the separate, red V(temp), it starts at a value and decreases, that's because the effect of the imposed initial condition is higher than the "rest" voltage, and, since the load sees little current drawn, the temperature slowly decreases to the equilibrum. The black voltage, V(Vin), starts a bit lower due to the initial warm conditions, but then rises towards the value given by reference-rest.

Then the current ramps up, temperature goes up, and here's the not-so-nice part: the plateau should exponentially rise as a 1-exp(-x), more than it's in the simulation, but not ad infinitum, thats why an RC ladder would have been nicer since you can add series RC, too.

Then the current drops, the temperature slowly decreases, to reach equilibrum. All this time, the voltage drops or rises according to the temperature variation.

Now, I realize this is probably less than a bare-bones show, but it's a start to show how to use the temperature in a feedback loop, seeing as you already have a VCCS in there, which can be put to good use. If you use primitives such as V+F for current sensor, and VCCS or VCVS plus RC elements, your simulation will fly, so the inclusion of such a circuit will barely hinder the rest of the schematic.