高频电磁场谐波分析，求高手！！

本人最近在做关于高频 电磁场的谐波分析，用于求耦合度的。简单来说，就是对于一个有3个端口的密闭空间。空间表面是非理想导体（电导率已知），我从一个端口放入一带电金属探针，并作为激励（输入功率为Pin=1W），然后在其余的两个口也放入一探针引出功率（分别记为P1和P2）。在输入口可以测到反射回来的功率（p-）。
按理说整个系统应该符合能量守恒，也就是
输入功率（Pin）— 反射功率（P-）= 表面导体的损耗功率(Pc)+引出功率（P1+P2）
但是我算出来竟然不符合能量守恒！！郁闷了很多天，不知道程序究竟哪里错误了。求高手提点！
之前已成功建模并划分网格，都没有问题。已下为边界条件及求解
！*****************定义理想导体面（电壁）*************************
VSEL,S,,,1,,,1 !选择导体空间
aslv !选择导体外表面
asel,U,area,,7
asel,U,area,,28
asel,U,area,,26
asel,U,area,,30
asel,U,area,,32
cm,area1,area
sfa,area1,,shld,0.66e7,0！定义外表面导体的电导率
！******************定义端口***************************************
sfa,26,,port,1 !Pin口
hfport,1,coax,,tem,ext,0.0076,0.0175,0,0,1 ！已同轴波导方式馈入作为激励，输入功率为1W
sfa,28,,port,2 !P1,引出端口
hfport,2,coax,,tem,impd,0.0076,0.0175,0,0
sfa,30,,port,3 !P2，引出端口
hfport,3,coax,,tem,impd,0.0076,0.040,0,0
et,1,hf119
et,2,0
asel,all
allsel,all
!*****************高频计算******************************************
/solu
antype,harmic
harfrq,322801531 ！频率是有高频模态分析仿真得到的，然后代入高频谐波分析
solve
finish
!*********功率计算***************************
set,1,1,
esel,s,sfe,port,2 ！得到P1口的引出功率
etable,pe,nmisc,4
ssum
esel,s,sfe,port,3 ！得到P2口的引出功率
etable,pe,nmisc,4
ssum
esel,s,sfe,port,1 ！得到Pin口的反射功率
etable,pe,nmisc,4
ssum
同时表面损耗也能通过公式计算出来，但是几部分相加竟都不符合能量守恒！请高手帮忙看看可能会存在问题的地方