Ни чего удивительного. Цитата из вики:
В 1919 г. была издана работа: «Основания технического расчёта пустотных катодных реле малой мощности», в которой излагалась разработанная М. А. Бонч-Бруевичем теория расчёта триода, ставшая основой теории электронных ламп и получившая позже название «теория Бонч-Бруевича — Баркгаузена».
До сих пор все расчётные соотношения получены на основании этой работы. Причём, чтобы, вообще, получить математический аппарат, автор, а так же последующие авторы, включая зарубежных, используют идеальный триод. Т.е. триод, у которого все ветви семейства характеристик представляют собой прямые линии, находящиеся на равном расстоянии друг от друга. У такого триода равное изменение сеточного напряжения приводят к равному изменению анодного тока. Чего, к сожалению, в реальности не существует. Но, иначе мат.аппарат "не вытанцовывается".
До сих пор ни кто не написал работы под названием "теория электровакуумного тетрода" , пентода, а тем более многосеточных ламп.
Вот это и является корнем проблемы даже для графического расчёта, а тем более для симуляции.
Ещё одна проблема: при расчёте гармонический искажений двухтактного УМ, при удачном расположении нагрузочной прямой, гармоники вычитаются, т.е. равно нулю, чего в реальности не может быть и классики предлагают ввести "коэффициент неидентичности". И, какой коэффициент вы вводите в формулу, такой Кг и получаете.
Плюс, многие программы считают Кг по упрощенному "методу трёх ординат". Что тоже не точно, т.к. даёт вычисление только по второй гармонике. Расчёт по "методу пяти ординат" даёт более точное значение, т.к. позволяет вычислить дополнительно амплитуды третьей и четвёртой гармоник. Что весьма актуально для тетродов и пентодов. Особенно, работающих в классе "А".