我核对了一下，81010中的计算确实错了，你这里给出的公式是正确的。我会在原贴中补充一条，指出这个错误。

81010的问题是，在公式(6)到公式(7)这步直接代入了初速度。因为公式(6)是有“初速度为零”的前提的，把初速代进去就出错了。

关于你更新的结果，我注意到线圈电流的曲线和顶楼不一致，是否是匝数被无意中修改了？和顶楼中给出的“180匝-2*级数”不同

更新后的效率是34.4%，和正确公式计算得到的35.3%的极限似乎过于接近了。直觉上讲，这种考虑了更多不完美因素的模型计算出的结果，不应该和理想条件下的极限那么接近才对。