发明名称 一种预测场线耦合系统响应的场-路协同模型建立方法
摘要 本发明公开了一种预测场线耦合系统响应的场-路协同模型建立方法,该场路协同模型建立方法是利用ADS平台、HFSS平台与单一干扰源辐射下的修正Taylor模型的结合,建立起从干扰设备辐射的电磁波经由线缆耦合到受扰设备的完整模型,从而获得在该干扰设备对线缆连接的两个终端受扰设备的响应。在明确了等效辐射源的状态下,本发明公开的场路协同模型建立方法可以用于预测任意干扰设备产生的电磁波激励下的线缆上的耦合响应情况,从而为设备级和系统级电磁兼容设计提供信息参考。
申请公布号 CN101697177B 申请公布日期 2012.09.05
申请号 CN200910236514.6 申请日期 2009.10.23
申请人 北京航空航天大学 发明人 陈文青;汪玉梅;吴龙刚;苏东林
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 周长琪
主权项 1.一种预测场线耦合系统响应的场-路协同模型建立方法,所述场线耦合系统是指一种干扰设备辐射的电磁波激励传输线缆终端设备产生的耦合响应的场-路协同模型;其特征在于:所述场-路协同模型是基于数值仿真平台的场-路协同模型建立方法;为了量化在干扰设备辐射的电磁场激励下的线缆连接的终端上的两个受扰设备受到电磁波影响的程度,包括有下列步骤:第一步:获取干扰设备系统级模型C<sub>i</sub>(In<sub>i</sub>,T<sub>i</sub>)干扰设备系统级模型C<sub>i</sub>(In<sub>i</sub>,T<sub>i</sub>)包含有干扰设备的个数i、干扰设备输入信号In<sub>i</sub>和干扰设备传递函数T<sub>i</sub>;第二步:获取干扰设备等效辐射源S<sub>i</sub>(G<sub>i</sub>,L<sub>i</sub>)干扰设备等效辐射源包括有等效辐射源的个数i、等效辐射源的激励G<sub>i</sub>(g<sub>i</sub>,φ<sub>i</sub>)和等效辐射源的尺寸L<sub>i</sub>,则任意一个等效辐射源表达形式为S<sub>i</sub>(G<sub>i</sub>,L<sub>i</sub>);所述等效辐射源的激励G<sub>i</sub>(g<sub>i</sub>,φ<sub>i</sub>)是第一步干扰设备系统级模型C<sub>i</sub>(In<sub>i</sub>,T<sub>i</sub>)的输出信号,且G<sub>i</sub>(G<sub>i</sub>,φ<sub>i</sub>)=In<sub>i</sub>×T<sub>i</sub>;该G<sub>i</sub>包括有激励电压的幅值g<sub>i</sub>和激励电压的相位φ<sub>i</sub>,i表示干扰设备的个数,即表达形式为G<sub>i</sub>(g<sub>i</sub>,φ<sub>i</sub>);所述等效辐射源的尺寸L<sub>i</sub>包括有等效辐射源的长LA<sub>i</sub>、宽LB<sub>i</sub>和高LC<sub>i</sub>或者包括等效辐射源的长LA<sub>i</sub>和半径LR<sub>i</sub>;第三步:获取干扰设备等效辐射源辐射总场信息E<sub>i</sub>(e<sub>i</sub>,θ<sub>i</sub>)将等效辐射源S<sub>i</sub>(G<sub>i</sub>,L<sub>i</sub>)、等效辐射源的激励G<sub>i</sub>(g<sub>i</sub>,φ<sub>i</sub>)和线缆边界KL(l,a,D,σ,ε,μ)引入基于数值仿真平台的场-路协同模型内,解析获得当前的干扰设备等效辐射源辐射总场信息E<sub>i</sub>(e<sub>i</sub>,θ<sub>i</sub>);其中,l表示在场-线耦合系统物理模型中线缆的长度,单位米;a表示在场-线耦合系统物理模型中线缆的半径,单位米;D表示在场-线耦合系统物理模型中线缆与地面的距离,单位米;σ表示线缆的传导率;ε表示线缆的介电常数;μ表示线缆的磁导率;在给定等效辐射源的激励G<sub>i</sub>(g<sub>i</sub>,φ<sub>i</sub>)下,干扰设备等效辐射源为圆柱形单极子天线的场分布为:在r≤λ/2π近场区时,等效辐射源的场分布<maths num="0001"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>H</mi><mo>=</mo><mfrac><mrow><mi>I</mi><mo>&times;</mo><msub><mi>LA</mi><mi>i</mi></msub></mrow><mrow><mn>2</mn><mi>&pi;</mi><msup><mi>r</mi><mn>2</mn></msup></mrow></mfrac></mtd></mtr><mtr><mtd><mi>E</mi><mo>=</mo><mfrac><mrow><mn>30</mn><mo>&times;</mo><mi>I</mi><mo>&times;</mo><msub><mi>LA</mi><mi>i</mi></msub><mo>&times;</mo><mi>&lambda;</mi></mrow><mrow><mi>&pi;</mi><msup><mi>r</mi><mn>3</mn></msup></mrow></mfrac></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>且<img file="FSB00000724504900021.GIF" wi="354" he="185" />其中,r表示观测点到等效辐射源的距离,单位为米;λ表示在仿真频率f下,该频率f所对应的波长;j表示虚数单位;在r>λ/2π远场区时,等效辐射源的场分布<maths num="0002"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>H</mi><mo>=</mo><mfrac><mrow><mi>I</mi><mo>&times;</mo><msub><mi>LA</mi><mi>i</mi></msub></mrow><mi>&lambda;r</mi></mfrac></mtd></mtr><mtr><mtd><mi>E</mi><mo>=</mo><mfrac><mrow><mn>120</mn><mi>&pi;</mi><mo>&times;</mo><mi>I</mi><mo>&times;</mo><msub><mi>LA</mi><mi>i</mi></msub></mrow><mi>&lambda;r</mi></mfrac></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>且<maths num="0003"><![CDATA[<math><mrow><mi>I</mi><mo>=</mo><mi>j</mi><mfrac><mrow><mi>&pi;</mi><msub><mi>G</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>g</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&phi;</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><mn>25</mn><mi>ln</mi><mrow><mo>(</mo><mfrac><msub><mi>LA</mi><mi>i</mi></msub><msub><mi>LR</mi><mi>i</mi></msub></mfrac><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>将线缆边界KL(l,a,D,σ,ε,μ)作为Maxwell方程<maths num="0004"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mo>&dtri;</mo><mo>&times;</mo><msub><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mi>i</mi></msub><mo>=</mo><mfrac><mrow><mo>&PartialD;</mo><mi>&epsiv;</mi><msub><mover><mi>E</mi><mo>&RightArrow;</mo></mover><mi>i</mi></msub></mrow><mrow><mo>&PartialD;</mo><mi>t</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mo>&dtri;</mo><mo>&times;</mo><msub><mover><mi>E</mi><mo>&RightArrow;</mo></mover><mi>i</mi></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mo>&PartialD;</mo><mi>&mu;</mi><msub><mover><mi>H</mi><mo>&RightArrow;</mo></mover><mi>i</mi></msub></mrow><mrow><mo>&PartialD;</mo><mi>t</mi></mrow></mfrac></mtd></mtr></mtable></mfenced></math>]]></maths>的边界条件,求解得到电磁波总场信息E<sub>i</sub>(e<sub>i</sub>,θ<sub>i</sub>);其中,<img file="FSB00000724504900025.GIF" wi="35" he="37" />表示拉普拉斯算子;<img file="FSB00000724504900026.GIF" wi="51" he="63" />表示沿线缆边界KL(l,a,D,σ,ε,μ)的磁场矢量;<img file="FSB00000724504900027.GIF" wi="41" he="64" />表示沿线缆边界KL(l,a,D,σ,ε,μ)的电场矢量;<img file="FSB00000724504900028.GIF" wi="107" he="120" />表示对线缆的介电常数ε和电场矢量<img file="FSB00000724504900029.GIF" wi="41" he="63" />在仿真时间t段内的微分变化量;<img file="FSB000007245049000210.GIF" wi="119" he="120" />表示线缆的磁导率μ和磁场矢量<img file="FSB000007245049000211.GIF" wi="51" he="64" />在仿真时间t段内的微分变化量;e<sub>i</sub>表示等效辐射源辐射总场信息中的电场幅度;θ<sub>i</sub>表示等效辐射源辐射总场信息中的电场相位;第四步:获取线缆终端电压V<sub>i</sub>(l)依据修正Taylor模型<maths num="0005"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mrow><mo>&PartialD;</mo><mi>V</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><mi>X</mi></mrow></mfrac><mo>+</mo><msup><mi>Z</mi><mo>&prime;</mo></msup><mi>I</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>j&omega;</mi><msub><mi>&mu;</mi><mn>0</mn></msub><msubsup><mo>&Integral;</mo><mn>0</mn><mi>D</mi></msubsup><msub><mi>H</mi><mi>Y</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>Z</mi><mo>)</mo></mrow><mi>dz</mi></mtd></mtr><mtr><mtd><mi>V</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>D</mi></msubsup><msub><mi>E</mi><mi>Z</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>Z</mi><mo>)</mo></mrow><mi>dz</mi></mtd></mtr></mtable></mfenced></math>]]></maths>对电磁波总场信息E<sub>i</sub>(e<sub>i</sub>,θ<sub>i</sub>)进行积分运算,从而得到线缆终端电压<img file="FSB000007245049000213.GIF" wi="466" he="88" />其中,<img file="FSB000007245049000214.GIF" wi="159" he="113" />表示线缆沿线电压V(X)沿X轴方向的微分变化量;Z′表示线缆内阻抗;I(X)表示线缆沿线电流;j表示虚数单位;ω表示角频率,且ω=2πf;μ<sub>0</sub>表示空气的磁导率;<img file="FSB000007245049000215.GIF" wi="317" he="94" />表示磁场分量在Z轴方向的积分;H<sub>Y</sub>(X,Z)表示线缆在(X,Z)位置上Y轴方向的总磁场分量;dz表示积分算子;<img file="FSB000007245049000216.GIF" wi="288" he="88" />表示电场分量在Z轴方向的积分;E<sub>Z</sub>(X,Z)表示线缆在(X,Z)位置上Z轴方向的总电场分量;<img file="FSB000007245049000217.GIF" wi="269" he="89" />表示在(l,Z)位置时电场分量在Z轴方向的积分。
地址 100083 北京市海淀区学院路37号