I'm modelling a used, linear PSU I have, and in a bias chain, LTS says it would have ~176.6uA of current. The closest I can calculate is 180.3uA, which is still only up by 2%. The bias chain alone, is just 2x 9V1zen, 2x 1N4148, and 160k, all in series between 48V total.
I'm using the zener voltage that LTspice gets, and then a starting guess for Vd of the 1N4148, and solving for the current from all those known voltages and R. Then using that current in the Shockley diode eqn (w/ the param's from LTS, Is=2.52nA , n=1.752, Rs=0.568 ), to get a new Vd. Then iterating that sequence from Vd to Id, until the difference between the new and last Vd, is below ~10E-15, and then the spreadsheet considers it zero.
I'm was using a 9V1 zener BZX84C9V1LY, it has a low Is, and in circuit LTS has Vz=9.0734329V, and my Vd calc gets to +0.143% of the LTS value, but the current is off by +2%
With KDZ9_1B or PTZ9_1B, w/ 1000x and 100,000x the Is of the above, I still only get within 0.3% for Vd and 2-3% for Id
So what other little fudge factor, or parameter, like Early voltage for BJT's , and the (1+VCE/VA) term for BJT's, can I use for a regular diode ?
For DC biasing, I've just worked backwards, calculating "everything" from the output of a 30V PSU, only ignoring a couple of things that I could have added in. And also using the zener voltages from LTS, and the effective beta of 2 driver BJT's.
And going around the voltage control loop, all my calculations are within +/-0.9% of the real model with a 2SC1085, TIP31C, and 2N3055, and LM301A, for ALL currents and voltages, even for nA levels of current.
But the thing that's off the most by about 6% or 6uA (in that circuit), is this bias chain with these zener's and diodes. If I can get all those V's and I's for the BJT's and op-amp within 0.9% or better (most of them are 99.9% or better) using similar methods and equations with the basic BJT model, I'd figure this bias chain would match even better than those.