[发明专利]一种快速预测车辆线束时域辐射敏感度的方法

权利要求书

1.一种快速预测车辆线束时域辐射敏感度的方法，其特征在于：其具体方法如下所述：

过程一：假设线束在无损耗的条件下，用镜像法求取带绝缘层的线束单位长度分布参数，在无绝缘层的导线的单位长度电感矩阵和电容矩阵有关系LC＝CL＝μ₀ε₀，其中，μ₀为真空磁导率，ε₀为真空介电常数，L为单位长度电感矩阵，C为单位长度电容矩阵，其中用镜像法可以求得：

其中，矩阵中对角线上的元素为导线的自电感

Lii=μ02πln(2hiri)---(2)]]>

矩阵中的非对角线上的元素为导线的互电感

Lij=μ02πln(dijdij)---(3)]]>

对于有绝缘层的导线需要重新推导单位长度电容矩阵，由电磁学可知导体电容Ｃ等于导体所带电量与导体电势的比值，即Ｃ＝ｑ/Ｖ，当导体为线电荷时单位长度电容等于线电荷密度ρ与电势V的比值，即

V＝ρ/C＝Sρ  (4)

式中：ρ为线电荷密度；S为电位系数；

电位系数矩阵和电容矩阵的关系为

其中，电位系数矩阵的对角线的元素为

Sii=12πϵ0(1ϵrln1ri+ϵeln1ri+Δri-ln12hi)---(6)]]>

电位系数矩阵的非对角线的元素为

Sij=12πϵ0ln(Xj-Xi)2+(Yj+Yi)2(Xj-Xi)2+(Yj-Yi)2=14πϵ0ln(1+4hihjdij2)---(7)]]>

当导体表面覆盖有绝缘层电介质时，导线上总的电荷密度ρ等于绝缘层外表面束缚电荷密度ρ_r＝a+b与导体与绝缘层上的界面上束缚电荷、自由电荷密度之和ρ_r＝a，即ρ＝ρ_r＝a+b+ρ_r＝a，导体与绝缘层的界面上的束缚电荷、自由电荷密度之和为ρ_r＝a＝ρ/ε_r，将其代入(2)式中得到绝缘层外表面束缚电荷密度ρ_r＝a+b＝ρ-ρ_r＝a＝ρ(ε_r-1)/ε_r，令(ε_r-1)/ε_r＝ε_e，故ρ_r＝a+b＝ρε_e，其中，ε_r是相对介电常数；

过程二：采用等效线束方法对带绝缘层的线束长度分布参数进行处理，得到等效线束的单位长度分布参数，同时对线束的几何参数，终端电阻进行等效；

第一步，将线束进行分组：

线束分组方法如表1所示：

表1线束分组方法

 端口1端口2组1|Z_1i|<Z_cm|Z_2i|<Z_cm组2|Z_1i|<Z_cm|Z_2i|>Z_cm组3|Z_1i|>Z_cm|Z_2i|<Z_cm组4|Z_1i|>Z_cm|Z_2i|>Z_cm

其中，Z_cm是表示线束共模特性阻抗模值，cm表示共模，在分组过程中，如果某线束中导线的终端所接负载非常接近或者等于共模特性阻抗，该导线的终端能量被此类负载完全吸收，从而对线缆的分组可忽略，则可将该线缆分在任意组中；

第二步，等效线束的单位长度电容矩阵及电感矩阵

设线束共有N根导线，以1号线束作为需要仿真的线束，按第一步的规则将其分为5组等效线束，其中，第一组为第1根导线，第2组等效导线由(α-1)根原始导线2-α组成，第2组等效导线由β根原始导线(α+1)-β组成，第3组等效导线由γ根原始导线(β+1)-γ组成，第4组等效导线由N根原始导线(γ+1)-N组成，电容矩阵和电感矩阵的表达式，如下(8)(9)式所示：

Ceq=C11Σq=2αC1qΣq=α+1βC1qΣq=β+1γC1qΣq=γ+1NC1qΣp=2αCp1Σp=2αΣq=2αCpqΣp=2αΣq=α+1βCpqΣp=2αΣq=β+1γCpqΣp=2αΣq=γ+1NCpqΣp=α+1βCp1Σp=α+1βΣq=2αCpqΣp=α+1βΣq=α+1βCpqΣp=α+1βΣq=β+1γCpqΣp=α+1βΣq=γ+1NCpqΣp=β+1γCp1Σp=β+1γΣq=2αCpqΣp=β+1γΣq=α+1βCpqΣp=β+1γΣq=β+1γCpqΣp=β+1γΣq=γ+1NCpqΣp=γ+1NCp1Σp=γ+1NΣq=2αCpqΣp=γ+1NΣq=α+1βCpqΣp=γ+1NΣq=β+1γCpqΣp=γ+1NΣQ=γ+1NCpq---(8)]]>

Leq=L11Σq=2αC1qα-1Σq=α+1βC1qβΣq=β+1γC1qγΣq=γ+1NC1qNΣp=2αCp1α-1Σp=2αΣq=2αCpq(α-1)2Σp=2αΣq=α+1βCpq(α-1)βΣp=2αΣq=β+1γCpq(α-1)γΣp=2αΣq=γ+1NCpq(α-1)NΣp=α+1βCp1βΣp=α+1βΣq=2αCpq(α-1)βΣp=α+1βΣq=α+1βCpqβ2Σp=α+1βΣq=β+1γCpqβγΣp=α+1βΣq=γ+1NCpqβNΣp=β+1γCp1γΣp=β+1γΣq=2αCpq(α-1)γΣp=β+1γΣq=α+1βCpqγβΣp=β+1γΣq=β+1γCpqγ2Σp=β+1γΣq=γ+1NCpqγNΣp=γ+1NCp1NΣp=γ+1NΣq=2αCpq(α-1)NΣp=γ+1NΣq=α+1βCpqβNΣp=γ+1NΣq=β+1γCpqγNΣp=γ+1NΣq=γ+1NCpqN2---(9)]]>

第三步，线束等效模型的位置坐标：

求取第i组所有线束的对地高度

hi=h1+h2+···hkK]]>

(10)]]>

等效线束的间距

第四步，对线束的终端负载进行等效：

线束的终端负载可以分为共模负载和差模负载，即单根线束与参考平面之间的负载定义为共模负载，线束之间的负载定义为差模负载；

等效线束的终端的等效阻抗可以被分为以下三类：

第一类，连接在线束终端和参考平面之间的共模负载的等效；

第二类，连接在同一等效线束中的线束之间的差模负载的等效；

第三类，连接在不同等效线束之间的差模负载；

过程三：将过程二所得参数代入传输线辐射敏感度时域有限差分公式进行运算，求取线束的时域辐射敏感度：

θ_E表示入射电场的入射方向与x轴的正向的夹角，θ_p表示入射电场的入射方向在yoz平面上投影与y轴的正向之间的夹角，φ_p表示入射电场的方向与a方向的夹角，分别代表传输线的源端阻抗及负载端阻抗；

入射电场与坐标位置以及频率有关，如(12)式所示入射电场：

E^→incident(x,y,z,ω)=E^o(ω)(exa→x+eya→y+eza→z)×e-jβxxe-jβyye-jβzz---(12)]]>

其中，ω是入射电场的频率，入射电场矢量在直角坐标系下沿x,y和z轴的各个分量为：

e_x＝sinθ_Esinθ_p

e_y＝-sinθ_Ecosθ_pcosφ_p-cosθ_Esinφ_p

e_z＝-sinθ_Ecosθ_psinφ_p+cosθ_Ecosφ_p

相位系数在坐标系各轴方向上的分量为：

β_x＝-βcosθ_p

β_y＝-βsinθ_pcosφ_p

β_z＝-βsinθ_psinφ_p

上式是入射电场关于频率的表达，所以(12)式转变为下式：

ξ→incident(x,y,z,t)=ξo(t-xvx-yvy-zvz)×[exa→x+eya→y+eza→z]---(13)]]>

电场的时间函数由ξ_o(t)表示，其中入射电场沿各坐标轴的传输速度为：

vx=ωβx=-vcosθp]]>

vy=ωβy=-vsinθpcosφp]]>

vz=ωβz=-vsinθpsinφp]]>

式中，v是入射电场的传输速度；

然后，将入射电场的时域表达与多导体传输线方程进行结合得到(14)、(15)式如下：

∂∂zV(z,t)+L∂∂t(z,t)=[1vzAT-AL]∂∂tξo(t-zvz)---(14)]]>

∂∂zI(z,t)+C∂∂tV(z,t)=-CAT∂∂tξo(t-zvz)---(15)]]>

式中，传输线上的电压V(z,t)和电流I(z,t)是与传输线上位置以及时间相关的函数，将(14)和(15)离散化，将整个传输线在空间上分成NDZ段，每段为Δz，即空间步长；将总的求解时间划分为NDT段，每段为Δt，即时间步长；将NDZ+1点电压V₁，V₂，···，V_NDZ,V_NDZ+1与NDT点电流I₁，I₂，···，I_NDZ作交织；每一个电压点和相邻的电流点间隔Δz/2；另外，时间点也必须进行交织，每一个电压时间点和相邻的电流时间点间隔Δt/2，由时域有限差分法方法得到上述方程的迭代式(16)-(19)：

V1n+1=(ΔzΔtRsC+1n)-1{(ΔzΔtRsC-1n)V1n-2RsI1n+1/2+(Vsn+1+Vsn)-ΔzΔtRsCAT[ξo(tn+1)-ξo(tn)]}---(16)]]>

VNDZ+1n+1=(ΔzΔtRLC+1n)-1{(ΔzΔtRLC-1n)VNDZ+1n+2RLINDZn+1/2+(VLn+1+VLn)-ΔzΔtRLCAT[ξo(tn+1-NDZΔtvz)-ξo(tn-NDZΔtvz)}---(17)]]>

Vkn+1=Vkn-ΔtΔzC-1[Ikn+1/2-Ik-1n+1/2]-AT[δo(tn+1-(k-1)Δzvz)-ξo(tn-(k-1)Δzvz)]---(18)]]>

Ikn+3/2=Ikn+1/2-ΔtΔzL-1[Vk+1n+1-Vkn+1]+L-1[1vzAT-AL][ξo(tn+3/2-(k-1/2)Δzvz)-ξo(tn+1/2-(k-1/2)Δzvz)]---(19)]]>

其中

AT=...(exxi+eyyi)...]]>

AL=...(xivx+yivy)ez...]]>

向量A_T，A_L中包含线束的几何位置坐标，其稳定性条件为其中v_iMAX是最大模式速度，将等效后的参数代入到(16)-(19)的迭代格式中就可以计算出线束的辐射敏感度。