Hi Jim,?
The approach I have taken is to gather the 2 port s parameters for a box full of say 3904. Then use the s data with a waveform arithmetic statement in LTSpice to get ft. As ft is an extrapolation of the low frequency gain on a 6dB per octave slope, it is pretty fast to evaluate. Recall, it is just a function of gm and Cbe; unity short circuit current gain crossing point. It does assist in validating the quality of some BIP models. I usually take the s data at the collector current where advertised ft peaks.?
Alan?