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

摘要

本发明公开了一种快速预测车辆线束时域辐射敏感度的方法，过程一假设线束在无损耗的条件下，用镜像法求取带绝缘层的线束单位长度分布参数，过程二采用等效线束方法对带绝缘层的线束长度分布参数进行处理，得到等效线束的单位长度分布参数，同时对线束的几何参数，终端电阻进行等效；过程三将过程二所得参数代入传输线辐射敏感度时域有限差分公式进行运算，求取线束的时域辐射敏感度；本发明的有益效果重新推得线束导线间电容矩阵的表达式；相对于完整模型计算时间缩短了50％左右；可以观测到线束上任意点的感应电压，入射电场的入射角度可以取实际过程中的任意值。

主权项

一种快速预测车辆线束时域辐射敏感度的方法，其特征在于：其具体方法如下所述：过程一：假设线束在无损耗的条件下，用镜像法求取带绝缘层的线束单位长度分布参数，在无绝缘层的导线的单位长度电感矩阵和电容矩阵有关系LC＝CL＝μ0ε0，其中，μ0为真空磁导率，ε0为真空介电常数，L为单位长度电感矩阵，C为单位长度电容矩阵，其中用镜像法可以求得：其中，矩阵中对角线上的元素为导线的自电感Lii=μ02πln(2hiri)---(2)]]>矩阵中的非对角线上的元素为导线的互电感Lij=μ02πln(dij`dij)---(3)]]>对于有绝缘层的导线需要重新推导单位长度电容矩阵，由电磁学可知导体电容C等于导体所带电量与导体电势的比值，即C＝q/V，当导体为线电荷时单位长度电容等于线电荷密度ρ与电势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)将线束剩余导线中终端阻抗满足|Z1i|＜Zcm，且|Z2i|＜Zcm的所有导线归为组a；(2)将线束剩余导线中终端阻抗满足|Z1i|＜Zcm，且|Z2i|＞Zcm的所有导线归为组b；(3)将线束剩余导线中终端阻抗满足|Z1i|＞Zcm，且|Z2i|＜Zcm的所有导线归为组c；(4)将线束剩余导线中终端阻抗满足|Z1i|＞Zcm，且|Z2i|＞Zcm的所有导线归为组d，其中，|Z1i|和|Z2i|表示线束中导线i两端的终端阻抗，Zcm是表示线束共模特性阻抗模值，cm表示共模，在分组过程中，如果某线束中导线的终端所接负载非常接近或者等于共模特性阻抗，该导线的终端能量被此类负载完全吸收，从而对线缆的分组可忽略，则可将该线缆分在任意组中；第二步，等效线束的单位长度电容矩阵及电感矩阵设线束共有N根导线，以1号线束作为需要仿真的线束，按第一步的规则将其分为5组等效线束，其中，第一组为1号导线，第二组等效导线由(α‑1)根原始导线2至α组成，第三组等效导线由(β‑α)根原始导线(α+1)至β组成，第四组等效导线由(γ‑β)根原始导线(β+1)至γ组成，第五组等效导线由(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)]]>等效线束的间距dij=4hihjexp(4π·Lij-eqμ0)-1---(11)]]>第四步，对线束的终端负载进行等效：线束的终端负载可以分为共模负载和差模负载，即单根线束与参考平面之间的负载定义为共模负载，线束之间的负载定义为差模负载；等效线束的终端的等效阻抗可以被分为以下三类：第一类，连接在线束终端和参考平面之间的共模负载的等效；第二类，连接在同一等效线束中的线束之间的差模负载的等效；第三类，连接在不同等效线束之间的差模负载；过程三：将过程二所得参数代入传输线辐射敏感度时域有限差分公式进行运算，求取线束的时域辐射敏感度：θE表示入射电场的入射方向与x轴的正向的夹角，θp表示入射电场的入射方向在yoz平面上投影与y轴的正向之间的夹角，φp表示入射电场的方向与方向的夹角，分别代表传输线的源端阻抗及负载端阻抗；入射电场与坐标位置以及频率有关，如(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轴的各个分量为：ex＝sinθE sinθpey＝‑sinθE cosθp cosφp‑cosθE sinφpez＝‑sinθE cosθp sinφp+cosθE cosφp相位系数在坐标系各轴方向上的分量为：βx＝‑βcosθpβy＝‑βsinθp cosφpβz＝‑βsinθp sinφ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点电压V1，V2，···，VNDZ,VNDZ+1与NDT点电流I1，I2，···，INDZ作交织；每一个电压点和相邻的电流点间隔Δ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...]]>向量AT，AL中包含线束的几何位置坐标，其稳定性条件为其中viMAX是最大模式速度，将等效后的参数代入到(16)‑(19)的迭代格式中就可以计算出线束的辐射敏感度。