CN108345753B - A crosstalk noise prediction method for non-parallel cables - Google Patents

A crosstalk noise prediction method for non-parallel cables Download PDF

Info

Publication number
CN108345753B
CN108345753B CN201810179567.8A CN201810179567A CN108345753B CN 108345753 B CN108345753 B CN 108345753B CN 201810179567 A CN201810179567 A CN 201810179567A CN 108345753 B CN108345753 B CN 108345753B
Authority
CN
China
Prior art keywords
matrix
crosstalk
cable
parallel
electromagnetic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810179567.8A
Other languages
Chinese (zh)
Other versions
CN108345753A (en
Inventor
张振华
颜伟
赵阳
王恩荣
陈昊
黄学军
郭晋安
赵俊平
傅春啸
马润青
朱家明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Normal University
Original Assignee
Nanjing Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Normal University filed Critical Nanjing Normal University
Priority to CN201810179567.8A priority Critical patent/CN108345753B/en
Publication of CN108345753A publication Critical patent/CN108345753A/en
Application granted granted Critical
Publication of CN108345753B publication Critical patent/CN108345753B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Resistance Or Impedance (AREA)
  • Communication Cables (AREA)

Abstract

本发明公开了一种针对非平行线缆的串扰噪声预测方法,根据实际线缆参数建立线绕串扰模型,使用分段级联法分析非平行线缆,分段比拟平直条件下的多导体线缆,以此建立非平行多导体电磁耦合模型;依据线缆串扰的常见解法,建立非平行线缆方程,提出一种等效的电磁参数矩阵,预测非平行线缆的近端和远端电磁串扰。通过预测电磁串扰大小及实际危害程度,对电气设备的高效运作具有重要的理论和工程意义,准确地预测串扰,有助于在设计的早期阶段,快速预测并发现潜在的串扰问题,防患于未然;同时,作为设计串扰抑制措施的依据,也能在电气设备制作完成后的测试阶段,为实施串扰抑制的补救措施提供理论依据,从而提高设备系统的稳定性。

Figure 201810179567

The invention discloses a crosstalk noise prediction method for non-parallel cables. A wire-wound crosstalk model is established according to actual cable parameters, a segmented cascade method is used to analyze the non-parallel cables, and segments are compared to multi-conductors under straight conditions. The non-parallel multi-conductor electromagnetic coupling model is established; according to the common solution method of cable crosstalk, the non-parallel cable equation is established, and an equivalent electromagnetic parameter matrix is proposed to predict the near-end and far-end of the non-parallel cable Electromagnetic Crosstalk. By predicting the magnitude of electromagnetic crosstalk and the actual degree of harm, it has important theoretical and engineering significance for the efficient operation of electrical equipment. Accurately predicting crosstalk is helpful to quickly predict and discover potential crosstalk problems in the early stage of design, preventing problems from At the same time, as a basis for designing crosstalk suppression measures, it can also provide a theoretical basis for implementing crosstalk suppression remedial measures in the testing stage after the electrical equipment is manufactured, thereby improving the stability of the equipment system.

Figure 201810179567

Description

一种针对非平行线缆的串扰噪声预测方法A crosstalk noise prediction method for non-parallel cables

技术领域technical field

本发明涉及电磁兼容技术领域,尤其涉及一种针对非平行线缆的串扰噪声预测方法。The invention relates to the technical field of electromagnetic compatibility, in particular to a crosstalk noise prediction method for non-parallel cables.

背景技术Background technique

在电磁兼容中,一般把电子设备对自身的干扰称为串扰。线缆串扰是指线缆间由于能量通过电磁之间的耦合转化而形成一种对原本信号的干扰。在工程使用中,线缆是大量存在的,近年来工业发展要求电力器件具有高效、高可靠、高性能等特点。那么势必需求电气设备“高集成度”和“高频率”,作为连接件的线缆必然置于一个狭小的空间内并且仍需传递极速的电能或信号,那么器件内线缆串扰往往是不可忽视的。In electromagnetic compatibility, the interference of electronic equipment to itself is generally called crosstalk. Cable crosstalk refers to the interference between cables to the original signal due to the conversion of energy through electromagnetic coupling. In engineering use, cables exist in large numbers. In recent years, industrial development requires power devices to have the characteristics of high efficiency, high reliability, and high performance. Then it is bound to require "high integration" and "high frequency" of electrical equipment. Cables as connectors must be placed in a small space and still need to transmit extremely fast power or signals, so the crosstalk of cables within the device is often not negligible. of.

现有的求解场线耦合产生的串扰包括直接求解麦克斯韦方程或其等价的二阶方程,具体的主流方法包括时域有限差分法、有限元法、矩量法等。时域有限差分法是用场量的中心差商近似逼近其对时间和空间的一阶偏微商,通过递推模拟波的时域传播过程,从而得出场分布。有限元法是将给定的拉普拉斯方程的求解问题,转化为求解泛函的极值问题。矩量法是一种求解线性微分方程的方法,将连续方程离散化为代数方程组,以矩阵变化求解为其数学基础。Existing solutions to crosstalk caused by field-line coupling include directly solving Maxwell's equations or their equivalent second-order equations, and specific mainstream methods include time-domain finite difference method, finite element method, and moment-of-moment method. The time-domain finite difference method is to approximate the first-order partial differential quotient of time and space with the central difference quotient of the field quantity, and to obtain the field distribution by recursively simulating the time-domain propagation process of the wave. The finite element method is to transform the given Laplace equation solution problem into the extreme value problem of solving the functional. The method of moments is a method for solving linear differential equations, which discretizes continuous equations into algebraic equations, and solves them with matrix changes as its mathematical basis.

目前对线缆串扰的研究已成体系,且对线缆串扰的电磁机理和最坏状态都有了较深入的研究。但是研究中大部分以空间布局较为理想(平直排布)的线缆束为研究对象,而现实工程的非常规布局线缆束是大量存在,总结来说缺少针对非平行线缆的研究。此外,中国发明专利(CN104007326A)公开了一种《一种快速预测车用线束串扰频域动态特性的方法》,能够实现车用线缆串扰频域动态特性的快速预测,但该专利技术也无法对非平行线缆的串扰进行预测。At present, the research on cable crosstalk has become a system, and the electromagnetic mechanism and worst state of cable crosstalk have been deeply studied. However, most of the researches take cable bundles with ideal spatial layout (straight arrangement) as the research object, and there are a large number of unconventional layout cable bundles in practical projects. In conclusion, there is a lack of research on non-parallel cables. In addition, the Chinese invention patent (CN104007326A) discloses a "method for rapidly predicting the dynamic characteristics of vehicle wire harness crosstalk in frequency domain", which can realize rapid prediction of the dynamic characteristics of vehicle cable crosstalk in frequency domain, but the patented technology cannot Prediction of crosstalk in non-parallel cables.

发明内容SUMMARY OF THE INVENTION

发明目的:本发明的目的是提供一种可以预测非平行线缆之间串扰噪声的方法。Purpose of the Invention: The purpose of the present invention is to provide a method for predicting crosstalk noise between non-parallel cables.

技术方案:本发明的针对非平行线缆的串扰噪声预测方法,包括以下步骤:Technical solution: The crosstalk noise prediction method for non-parallel cables of the present invention includes the following steps:

第一步:获取非平行线缆之间的空间位置参数;Step 1: Obtain the spatial position parameters between non-parallel cables;

第二步:建立单位长度的多导体模型,并列出线缆相应单位长度的阻抗矩阵Z和导纳矩阵Y关于线缆基本电磁矩阵的解析式,结合多导体传输线理论,代理空间参数,推导线缆基本电磁矩阵元素解析式;The second step: establish a multi-conductor model of unit length, and list the impedance matrix Z and admittance matrix Y of the corresponding unit length of the cable. The analytical formula for the basic electromagnetic matrix of the cable, combined with the theory of multi-conductor transmission lines, proxy space parameters, deduce Analytical formula of basic electromagnetic matrix elements of cables;

第三步:基于分段级联的非平行多导体串扰预测方法,使用类比的方法,对非平行线缆微分化,即分段级联,将非平行电磁参数矩阵用平行电磁参数矩阵等效,求解电磁等效矩阵;The third step: the non-parallel multi-conductor crosstalk prediction method based on the segmented cascade, using the analogy method, differentiate the non-parallel cables, that is, segmented cascade, and use the non-parallel electromagnetic parameter matrix to be equivalent to the parallel electromagnetic parameter matrix. , solve the electromagnetic equivalent matrix;

第四步:结合边界条件矩阵,求解出非平行线缆串扰。Step 4: Combine the boundary condition matrix to solve the non-parallel cable crosstalk.

其中,第一步中所述空间位置参数包括倾斜角度、近端距离、远端距离、线缆长度和离地高度。Wherein, the spatial location parameters in the first step include tilt angle, near-end distance, far-end distance, cable length, and ground clearance.

第二步中,所述阻抗矩阵Z和导纳矩阵Y的解析式为:In the second step, the analytical formulas of the impedance matrix Z and the admittance matrix Y are:

Figure GDA0003162416770000021
Figure GDA0003162416770000021

其中,R为电阻矩阵,L为电感矩阵,G为电导矩阵,C为电容矩阵,ω=2πf,f为频率,R、L、G、C为线缆的基本参数矩阵,由于表示电磁感应的G矩阵很小,一般直接忽略不计,当R、L、C矩阵求出,那么Z、Y矩阵也可以求解。Among them, R is the resistance matrix, L is the inductance matrix, G is the conductance matrix, C is the capacitance matrix, ω=2πf, f is the frequency, and R, L, G, and C are the basic parameter matrices of the cable. The G matrix is very small and is generally ignored directly. When the R, L, and C matrices are calculated, the Z and Y matrices can also be solved.

所述电阻矩阵R的元素计算公式为:The element calculation formula of the resistance matrix R is:

Figure GDA0003162416770000022
Figure GDA0003162416770000022

(2)式依照实芯线的半径与集肤深度的关系分成两部分,其中,σ为线芯导体的电导率,rw为实芯圆导体的半径,δ为集肤深度;对于绞合线,可以认为是多股半径为rws的实芯线的并联(rws为一股线芯的半径)。Formula (2) is divided into two parts according to the relationship between the radius of the solid core wire and the skin depth, where σ is the electrical conductivity of the wire core conductor, rw is the radius of the solid core circular conductor, and δ is the skin depth; , which can be considered as a parallel connection of multiple solid core wires with a radius of rws (rws is the radius of one wire core).

电感矩阵L的元素计算公式为:The element calculation formula of the inductance matrix L is:

Figure GDA0003162416770000023
Figure GDA0003162416770000023

其中,hi和hj为线缆离地高度,sij为两条线缆间实际距离;Among them, h i and h j are the heights of the cables from the ground, and s ij is the actual distance between the two cables;

sij=b+z·tan(a)s ij =b+z·tan(a)

其中,a为线缆相交角度,b为线缆间最小距离,z为横向距离。Among them, a is the intersection angle of cables, b is the minimum distance between cables, and z is the lateral distance.

所述第三步具体为将第二步中求得的电磁矩阵元素代入电阻矩阵R和电感矩阵L的解析式中,其中,电阻矩阵R和电感矩阵L的解析式为:The third step is to substitute the electromagnetic matrix elements obtained in the second step into the analytical expressions of the resistance matrix R and the inductance matrix L, wherein the analytical expressions of the resistance matrix R and the inductance matrix L are:

Figure GDA0003162416770000031
Figure GDA0003162416770000031

Figure GDA0003162416770000032
Figure GDA0003162416770000032

求得电感矩阵L后,代入电容矩阵C的解析式:After obtaining the inductance matrix L, substitute it into the analytical formula of the capacitance matrix C:

C=μεL-1 (6)C=μεL -1 (6)

其中,ε为介电常数,μ为导磁系数;Among them, ε is the dielectric constant, μ is the permeability coefficient;

将求得的电阻矩阵R、电感矩阵L以及电容矩阵C代入第二步中阻抗矩阵Z和导纳矩阵Y的解析式中,求得Z和Y,通过先距离z积分,再除以总长度d,求出电磁等效矩阵Zd和Yd,如下:Substitute the obtained resistance matrix R, inductance matrix L and capacitance matrix C into the analytical formulas of impedance matrix Z and admittance matrix Y in the second step, to obtain Z and Y, by first integrating the distance z, and then dividing by the total length d, find the electromagnetic equivalent matrix Z d and Y d , as follows:

Figure GDA0003162416770000033
Figure GDA0003162416770000033

其中,z为横向距离,d为线缆长度。Among them, z is the lateral distance, and d is the cable length.

所述第四步具体为:The fourth step is specifically:

将(1)式求解一阶微分方程组得:Solving the first-order differential equations from equation (1), we get:

Figure GDA0003162416770000034
Figure GDA0003162416770000034

根据式(8)求解出等效于平行线缆的电磁矩阵;According to formula (8), the electromagnetic matrix equivalent to the parallel cable is solved;

列出边界条件方程:List the boundary condition equations:

Figure GDA0003162416770000035
Figure GDA0003162416770000035

其中,

Figure GDA0003162416770000036
Figure GDA0003162416770000037
为传输线两端z=0和z=d的戴维南等效电压源,同时
Figure GDA0003162416770000038
为非平行线缆远端串扰矩阵;该式显示了端口的戴维南等效,线缆两端的端口连接的各种电路都可以等效为电源和内电阻的形式,从而简化分析过程,in,
Figure GDA0003162416770000036
and
Figure GDA0003162416770000037
is the Thevenin equivalent voltage source with z=0 and z=d at both ends of the transmission line, while
Figure GDA0003162416770000038
is the non-parallel cable far-end crosstalk matrix; this formula shows the Thevenin equivalent of the port, and various circuits connected to the ports at both ends of the cable can be equivalent to the form of power supply and internal resistance, thus simplifying the analysis process,

将等效电磁参数矩阵结合平行线缆串扰的二阶求解矩量法,通过(8)式和(9)式得出:Combining the equivalent electromagnetic parameter matrix with the second-order method of moments to solve the crosstalk of parallel cables, the equations (8) and (9) are used to obtain:

Figure GDA0003162416770000041
Figure GDA0003162416770000041

其中d为线缆长度,e为自然常数,γ为对角阵,其对角线元素为矩阵(YZ)的特征根,P为由矩阵(YZ)的特征根对应的特征向量组成的矩阵;where d is the cable length, e is a natural constant, γ is a diagonal matrix, and its diagonal elements are the eigenvalues of the matrix (YZ), and P is a matrix composed of the eigenvectors corresponding to the eigenvalues of the matrix (YZ);

Figure GDA0003162416770000042
Figure GDA0003162416770000042

Figure GDA0003162416770000043
Figure GDA0003162416770000043

Figure GDA0003162416770000044
Figure GDA0003162416770000044

根据(10)式、(11)式、(12)式和(13)式依次求解得出串扰矩阵

Figure GDA0003162416770000045
According to formula (10), formula (11), formula (12) and formula (13), the crosstalk matrix is obtained by solving in sequence
Figure GDA0003162416770000045

有益效果:与现有技术相比,本发明的优点为:首先,本发明建立基于分段级联非平行线缆模型,提出一种等效的电磁参数矩阵,将非平行的电磁参数矩阵用平行线缆参数等效,从而解决预测非平行线束串扰的问题;其次,通过预测电磁串扰大小及实际危害程度,对电气设备的高效运作具有重要的理论和工程意义,准确地预测串扰,有助于在设计的早期阶段,快速预测并发现潜在的串扰问题,防患于未然;同时,作为设计串扰抑制措施的依据,也能在电气设备制作完成后的测试阶段,为实施串扰抑制的补救措施提供理论依据,从而提高设备系统的稳定性。Beneficial effects: Compared with the prior art, the advantages of the present invention are: firstly, the present invention establishes a non-parallel cable model based on segmented cascade, and proposes an equivalent electromagnetic parameter matrix, which uses the non-parallel electromagnetic parameter matrix as The parameters of parallel cables are equivalent, so as to solve the problem of predicting the crosstalk of non-parallel wire bundles; secondly, by predicting the magnitude of electromagnetic crosstalk and the actual degree of harm, it has important theoretical and engineering significance for the efficient operation of electrical equipment. In the early stage of design, it can quickly predict and discover potential crosstalk problems and prevent them from happening. At the same time, as the basis for designing crosstalk suppression measures, it can also be used to implement crosstalk suppression remedial measures in the testing stage after the electrical equipment is fabricated. Provide theoretical basis to improve the stability of the equipment system.

附图说明Description of drawings

图1为非平行线缆之间的空间位置参数获取示意图;FIG. 1 is a schematic diagram of obtaining spatial position parameters between non-parallel cables;

图2为平行线缆串扰分析模型示意图;Figure 2 is a schematic diagram of a parallel cable crosstalk analysis model;

图3为非平行线缆分段级联模型示意图;3 is a schematic diagram of a non-parallel cable segment cascade model;

图4为线缆串扰边界条件示意图。FIG. 4 is a schematic diagram of boundary conditions of cable crosstalk.

具体实施方式Detailed ways

下面结合附图对发明的技术方案作进一步说明。The technical solution of the invention will be further described below with reference to the accompanying drawings.

一种针对非平行线缆的串扰噪声预测方法,包括以下步骤:A method for predicting crosstalk noise for non-parallel cables, comprising the following steps:

第一步:如图1所示,利用位置传感器获得线缆之间的空间位置信息,包括倾斜角度,近端距离,远端距离,线缆长度,离地高度等。具体传感器型号视预测的线缆类型而定,在本发明中不作出详细介绍,但只需能够测出相关空间参数即可;The first step: As shown in Figure 1, use the position sensor to obtain the spatial position information between the cables, including the inclination angle, the near-end distance, the far-end distance, the cable length, and the height above the ground. The specific sensor model depends on the predicted cable type, which is not described in detail in the present invention, but only needs to be able to measure the relevant spatial parameters;

第二步:建立单位长度的多导体模型,并列出线缆相应单位长度的阻抗矩阵Z和导纳矩阵Y关于线缆基本电磁矩阵的解析式,结合多导体传输线理论,代理空间参数,推导线缆基本电磁矩阵元素解析式;在满足工程精度的基础上,可以利用解析表达式求解每一段的电磁矩阵。如图2所示的平行线缆串扰模型中rii和rjj表示线缆单位长度的电阻,lii和ljj代表线缆单位长度的自感,lij则代表线缆间单位长度的互感;cii、cjj和gii、gjj分别代表线缆单位长度的自容和自导,cij和gij则代表相应的互容和互导。Z代表了线缆横跨总长度,dz代表该段MTLs的微分长度,每根线缆相应的电压、电流相量同样出。The second step: establish a multi-conductor model of unit length, and list the impedance matrix Z and admittance matrix Y of the corresponding unit length of the cable. The analytical formula for the basic electromagnetic matrix of the cable, combined with the theory of multi-conductor transmission lines, proxy space parameters, deduce The basic electromagnetic matrix element analytical formula of the cable; on the basis of satisfying the engineering accuracy, the electromagnetic matrix of each section can be solved by the analytical expression. In the parallel cable crosstalk model shown in Figure 2, r ii and r jj represent the resistance per unit length of the cable, l ii and l jj represent the self-inductance per unit length of the cable, and l ij represent the mutual inductance between the cables per unit length ; c ii , c jj and g ii , g jj represent the self-capacitance and self-conductance per unit length of the cable, respectively, and c ij and g ij represent the corresponding mutual capacitance and mutual conductance. Z represents the total length of the cable spanning, dz represents the differential length of the MTLs in this section, and the corresponding voltage and current phasors of each cable are also output.

所述阻抗矩阵Z和导纳矩阵Y的解析式为:The analytical formulas of the impedance matrix Z and the admittance matrix Y are:

Figure GDA0003162416770000051
Figure GDA0003162416770000051

其中,R为电阻矩阵,L为电感矩阵,G为电导矩阵,C为电容矩阵,ω=2πf,f为频率,R、L、G、C为线缆的基本参数矩阵,由于表示电磁感应的G矩阵很小,一般直接忽略不计,当R、L、C矩阵求出,那么Z、Y矩阵也可以求解。Among them, R is the resistance matrix, L is the inductance matrix, G is the conductance matrix, C is the capacitance matrix, ω=2πf, f is the frequency, and R, L, G, and C are the basic parameter matrices of the cable. The G matrix is very small and is generally ignored directly. When the R, L, and C matrices are calculated, the Z and Y matrices can also be solved.

所述电阻矩阵R的元素计算公式为:The element calculation formula of the resistance matrix R is:

Figure GDA0003162416770000052
Figure GDA0003162416770000052

(2)式依照实芯线的半径与集肤深度的关系分成两部分,其中,σ为线芯导体的电导率,rw为实芯圆导体的半径,δ为集肤深度;对于绞合线,可以认为是多股半径为rws的实芯线的并联(rws为一股线芯的半径)。Formula (2) is divided into two parts according to the relationship between the radius of the solid core wire and the skin depth, where σ is the electrical conductivity of the wire core conductor, rw is the radius of the solid core circular conductor, and δ is the skin depth; , which can be considered as a parallel connection of multiple solid core wires with a radius of rws (rws is the radius of one wire core).

电感矩阵L的元素计算公式为:The element calculation formula of the inductance matrix L is:

Figure GDA0003162416770000053
Figure GDA0003162416770000053

其中,hi和hj为线缆离地高度,sij为两条线缆间实际距离;Among them, h i and h j are the heights of the cables from the ground, and s ij is the actual distance between the two cables;

sij=b+z·tan(a)s ij =b+z·tan(a)

其中,a为线缆相交角度,b为线缆间最小距离,z为横向距离Among them, a is the intersection angle of cables, b is the minimum distance between cables, and z is the lateral distance

第三步:基于分段级联的非平行多导体串扰预测方法,使用类比的方法,对非平行线缆微分化,即分段级联,将非平行电磁参数矩阵用平行电磁参数矩阵等效,求解电磁等效矩阵;如图3所示,将线缆截取成m份。若m值取无限大,那么可以将导体间的每一段导体都是相互平行的。由此可见,每一段的平行多导体传输线方程依然成立。The third step: the non-parallel multi-conductor crosstalk prediction method based on the segmented cascade, using the analogy method, differentiate the non-parallel cables, that is, segmented cascade, and use the non-parallel electromagnetic parameter matrix to be equivalent to the parallel electromagnetic parameter matrix. , solve the electromagnetic equivalent matrix; as shown in Figure 3, cut the cable into m parts. If the value of m is infinite, then each section of the conductors can be parallel to each other. It can be seen that the equation of the parallel multi-conductor transmission line for each segment is still valid.

将第二步中求得的电磁矩阵元素代入电阻矩阵R和电感矩阵L的解析式中,其中,电阻矩阵R和电感矩阵L的解析式为:Substitute the electromagnetic matrix elements obtained in the second step into the analytical expressions of the resistance matrix R and the inductance matrix L, where the analytical expressions of the resistance matrix R and the inductance matrix L are:

Figure GDA0003162416770000061
Figure GDA0003162416770000061

Figure GDA0003162416770000062
Figure GDA0003162416770000062

求得电感矩阵L后,代入电容矩阵C的解析式:After obtaining the inductance matrix L, substitute it into the analytical formula of the capacitance matrix C:

C=μεL-1 (6)C=μεL -1 (6)

其中,ε为介电常数,μ为导磁系数;Among them, ε is the dielectric constant, μ is the permeability coefficient;

将求得的电阻矩阵R、电感矩阵L以及电容矩阵C代入第二步中阻抗矩阵Z和导纳矩阵Y的解析式中,求得Z和Y,通过先距离z积分,再除以总长度d,求出电磁等效矩阵Zd和Yd,如下:Substitute the obtained resistance matrix R, inductance matrix L and capacitance matrix C into the analytical formulas of impedance matrix Z and admittance matrix Y in the second step, to obtain Z and Y, by first integrating the distance z, and then dividing by the total length d, find the electromagnetic equivalent matrix Z d and Y d , as follows:

Figure GDA0003162416770000063
Figure GDA0003162416770000063

其中,z为横向距离,d为线缆长度where z is the lateral distance and d is the cable length

第四步:结合边界条件矩阵,求解出非平行线缆串扰。Step 4: Combine the boundary condition matrix to solve the non-parallel cable crosstalk.

将(1)式求解一阶微分方程组得:Solving the first-order differential equations from equation (1), we get:

Figure GDA0003162416770000064
Figure GDA0003162416770000064

根据式(8)求解出等效于平行线缆的电磁矩阵;According to formula (8), the electromagnetic matrix equivalent to the parallel cable is solved;

根据图4所示的边界条件矩阵,列出边界条件方程:According to the boundary condition matrix shown in Figure 4, the boundary condition equations are listed:

Figure GDA0003162416770000071
Figure GDA0003162416770000071

其中,

Figure GDA0003162416770000072
Figure GDA0003162416770000073
为传输线两端z=0和z=d的戴维南等效电压源,同时
Figure GDA0003162416770000074
为非平行线缆远端串扰矩阵;该式显示了端口的戴维南等效,线缆两端的端口连接的各种电路都可以等效为电源和内电阻的形式,从而简化分析过程,in,
Figure GDA0003162416770000072
and
Figure GDA0003162416770000073
is the Thevenin equivalent voltage source with z=0 and z=d at both ends of the transmission line, while
Figure GDA0003162416770000074
is the non-parallel cable far-end crosstalk matrix; this formula shows the Thevenin equivalent of the port, and various circuits connected to the ports at both ends of the cable can be equivalent to the form of power supply and internal resistance, thus simplifying the analysis process,

将等效电磁参数矩阵结合平行线缆串扰的二阶求解矩量法,通过(8)式和(9)式得出:Combining the equivalent electromagnetic parameter matrix with the second-order method of moments to solve the crosstalk of parallel cables, the equations (8) and (9) are used to obtain:

Figure GDA0003162416770000075
Figure GDA0003162416770000075

其中d为线缆长度,e为自然常数,γ为对角阵,其对角线元素为矩阵(YZ)的特征根,P为由矩阵(YZ)的特征根对应的特征向量组成的矩阵;where d is the cable length, e is a natural constant, γ is a diagonal matrix, and its diagonal elements are the eigenvalues of the matrix (YZ), and P is a matrix composed of the eigenvectors corresponding to the eigenvalues of the matrix (YZ);

Figure GDA0003162416770000076
Figure GDA0003162416770000076

Figure GDA0003162416770000077
Figure GDA0003162416770000077

Figure GDA0003162416770000078
Figure GDA0003162416770000078

根据(10)式、(11)式、(12)式和(13)式依次求解得出串扰矩阵

Figure GDA0003162416770000079
According to formula (10), formula (11), formula (12) and formula (13), the crosstalk matrix is obtained by solving in sequence
Figure GDA0003162416770000079

Claims (1)

1.一种针对非平行线缆的串扰噪声预测方法,其特征在于包括以下步骤:1. A method for predicting crosstalk noise for non-parallel cables, characterized by comprising the following steps: 第一步:获取非平行线缆之间的空间位置参数;Step 1: Obtain the spatial position parameters between non-parallel cables; 第二步:建立单位长度的多导体模型,并列出线缆相应单位长度的阻抗矩阵Z和导纳矩阵Y关于线缆基本电磁矩阵的解析式,推导线缆基本电磁矩阵元素解析式;The second step: establish a multi-conductor model of unit length, and list the analytical formula of the impedance matrix Z and admittance matrix Y of the corresponding unit length of the cable about the basic electromagnetic matrix of the cable, and deduce the analytical formula of the basic electromagnetic matrix element of the cable; 第三步:将非平行电磁参数矩阵用平行电磁参数矩阵等效,求解电磁等效矩阵;The third step: the non-parallel electromagnetic parameter matrix is equivalent to the parallel electromagnetic parameter matrix, and the electromagnetic equivalent matrix is solved; 第四步:结合边界条件矩阵,求解出非平行线缆串扰;Step 4: Combine the boundary condition matrix to solve the non-parallel cable crosstalk; 第一步中,所述空间位置参数包括倾斜角度、近端距离、远端距离、线缆长度和离地高度;In the first step, the spatial position parameters include the inclination angle, the near-end distance, the far-end distance, the cable length and the height above the ground; 第二步中,所述阻抗矩阵Z和导纳矩阵Y的解析式为:In the second step, the analytical formulas of the impedance matrix Z and the admittance matrix Y are:
Figure FDA0003162416760000011
Figure FDA0003162416760000011
其中,R为电阻矩阵,L为电感矩阵,G为电导矩阵,C为电容矩阵,ω=2πf,f为频率;Among them, R is the resistance matrix, L is the inductance matrix, G is the conductance matrix, C is the capacitance matrix, ω=2πf, and f is the frequency; 所述电阻矩阵R的元素计算公式为:The element calculation formula of the resistance matrix R is:
Figure FDA0003162416760000012
Figure FDA0003162416760000012
其中,σ为线芯导体的电导率,rw为实芯圆导体的半径,δ为集肤深度;Among them, σ is the electrical conductivity of the wire core conductor, r w is the radius of the solid core circular conductor, and δ is the skin depth; 电感矩阵L的元素计算公式为:The element calculation formula of the inductance matrix L is:
Figure FDA0003162416760000013
Figure FDA0003162416760000013
其中,hi和hj为线缆离地高度,sij为两条线缆间实际距离;Among them, h i and h j are the heights of the cables from the ground, and s ij is the actual distance between the two cables; sij=b+z·tan(a)s ij =b+z·tan(a) 其中,a为线缆相交角度,b为线缆间最小距离,z为横向距离;Among them, a is the intersection angle of cables, b is the minimum distance between cables, and z is the lateral distance; 所述第三步具体为将第二步中求得的电磁矩阵元素代入电阻矩阵R和电感矩阵L的解析式中,其中,电阻矩阵R和电感矩阵L的解析式为:The third step is to substitute the electromagnetic matrix elements obtained in the second step into the analytical expressions of the resistance matrix R and the inductance matrix L, wherein the analytical expressions of the resistance matrix R and the inductance matrix L are:
Figure FDA0003162416760000021
Figure FDA0003162416760000021
Figure FDA0003162416760000022
Figure FDA0003162416760000022
求得电感矩阵L后,代入电容矩阵C的解析式:After obtaining the inductance matrix L, substitute it into the analytical formula of the capacitance matrix C: C=μεL-1 (6)C=μεL -1 (6) 其中,ε为介电常数,μ为导磁系数;Among them, ε is the dielectric constant, μ is the permeability coefficient; 将求得的电阻矩阵R、电感矩阵L以及电容矩阵C代入第二步中阻抗矩阵Z和导纳矩阵Y的解析式中,求得Z和Y,再求出电磁等效矩阵Zd和Yd,如下:Substitute the obtained resistance matrix R, inductance matrix L and capacitance matrix C into the analytical formulas of impedance matrix Z and admittance matrix Y in the second step, obtain Z and Y, and then obtain electromagnetic equivalent matrix Z d and Y d , as follows:
Figure FDA0003162416760000023
Figure FDA0003162416760000023
其中,z为横向距离,d为线缆长度;Among them, z is the lateral distance, d is the cable length; 所述第四步具体为:The fourth step is specifically: 将(1)式求解一阶微分方程组得:Solving the first-order differential equations from equation (1), we get:
Figure FDA0003162416760000024
Figure FDA0003162416760000024
根据式(8)求解出等效于平行线缆的电磁矩阵;According to formula (8), the electromagnetic matrix equivalent to the parallel cable is solved; 列出边界条件方程:List the boundary condition equations:
Figure FDA0003162416760000025
Figure FDA0003162416760000025
其中,
Figure FDA0003162416760000026
Figure FDA0003162416760000027
为传输线两端z=0和z=d的戴维南等效电压源,同时
Figure FDA0003162416760000028
为非平行线缆远端串扰矩阵;
in,
Figure FDA0003162416760000026
and
Figure FDA0003162416760000027
is the Thevenin equivalent voltage source with z=0 and z=d at both ends of the transmission line, while
Figure FDA0003162416760000028
is the non-parallel cable far-end crosstalk matrix;
将等效电磁参数矩阵结合平行线缆串扰的二阶求解矩量法,通过(8)式和(9)式得出:Combining the equivalent electromagnetic parameter matrix with the second-order method of moments to solve the crosstalk of parallel cables, the equations (8) and (9) are used to obtain:
Figure FDA0003162416760000031
Figure FDA0003162416760000031
其中d为线缆长度,e为自然常数,γ为对角阵,其对角线元素为矩阵(YZ)的特征根,P为由矩阵(YZ)的特征根对应的特征向量组成的矩阵;where d is the cable length, e is a natural constant, γ is a diagonal matrix, and its diagonal elements are the eigenvalues of the matrix (YZ), and P is a matrix composed of the eigenvectors corresponding to the eigenvalues of the matrix (YZ);
Figure FDA0003162416760000032
Figure FDA0003162416760000032
Figure FDA0003162416760000033
Figure FDA0003162416760000033
Figure FDA0003162416760000034
Figure FDA0003162416760000034
根据(10)式、(11)式、(12)式和(13)式依次求解得出串扰矩阵
Figure FDA0003162416760000035
According to formula (10), formula (11), formula (12) and formula (13), the crosstalk matrix is obtained by solving in sequence
Figure FDA0003162416760000035
CN201810179567.8A 2018-03-05 2018-03-05 A crosstalk noise prediction method for non-parallel cables Active CN108345753B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810179567.8A CN108345753B (en) 2018-03-05 2018-03-05 A crosstalk noise prediction method for non-parallel cables

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810179567.8A CN108345753B (en) 2018-03-05 2018-03-05 A crosstalk noise prediction method for non-parallel cables

Publications (2)

Publication Number Publication Date
CN108345753A CN108345753A (en) 2018-07-31
CN108345753B true CN108345753B (en) 2021-09-07

Family

ID=62957840

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810179567.8A Active CN108345753B (en) 2018-03-05 2018-03-05 A crosstalk noise prediction method for non-parallel cables

Country Status (1)

Country Link
CN (1) CN108345753B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109710995A (en) * 2018-12-07 2019-05-03 江苏益邦电力科技有限公司 A kind of crosstalk noise prediction technique for random arrangement cable
CN109992819A (en) * 2019-01-25 2019-07-09 南京师范大学镇江创新发展研究院 A kind of cable cross-talk processing method based on adaptive intersection approximate algorithm
CN110531171B (en) * 2019-08-28 2021-04-06 湖南大学 A calculation method for determining the critical wiring spacing of cable crosstalk
CN110516362A (en) * 2019-08-28 2019-11-29 哈尔滨工程大学 A FDTD Solution Method for Electromagnetic Field Distribution of Multiconductor Transmission Lines
CN110554247A (en) * 2019-08-28 2019-12-10 广州广电计量检测股份有限公司 crosstalk simulation modeling method and simulation model of non-parallel cable
CN110909481B (en) * 2019-11-29 2023-05-23 北京应用物理与计算数学研究所 Complex cable bundle distribution parameter modeling method based on moment method
CN111368436B (en) * 2020-03-06 2023-07-21 重庆邮电大学 Time Domain Modeling and Analysis Method for Electromagnetic Coupling Effect of Bending Lines on Conductive Plate
CN114722552B (en) * 2022-06-09 2022-11-08 深圳荣耀智能机器有限公司 Cable length verification method and electronic equipment

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7733198B1 (en) * 2007-05-15 2010-06-08 Sandia Corporation Microfabricated bulk wave acoustic bandgap device
CN103048563B (en) * 2012-12-19 2016-02-10 南京航空航天大学 The method for quick predicting of a kind of right angle above the ground cable bundle terminal crosstalk
CN104007327B (en) * 2014-06-18 2017-03-15 吉林大学 A kind of method of fast prediction vehicle harness time-domain radiation susceptibility
CN104849572B (en) * 2015-05-24 2017-11-03 浙江大学 A kind of HW High Way cross talk restraining method decomposed based on electromagnetic field mode
CN105956309B (en) * 2016-05-13 2019-05-14 南京航空航天大学 A kind of " encoder " design method for eliminating the stable state crosstalk of transmission line distal end

Also Published As

Publication number Publication date
CN108345753A (en) 2018-07-31

Similar Documents

Publication Publication Date Title
CN108345753B (en) A crosstalk noise prediction method for non-parallel cables
Gustavsen et al. Parameter determination for modeling system transients-Part II: Insulated cables
CN105606899B (en) A kind of extracting method of Variable Frequency Drives motor side common code impedance
CN110531171B (en) A calculation method for determining the critical wiring spacing of cable crosstalk
Xiao et al. Analysis of crosstalk between finite-length microstrip lines: FDTD approach and circuit-concept modeling
CN105243197B (en) Single-core power cables harmonic analysis method based on finite element simulation and equivalent circuit
CN107609251B (en) FDS-based modeling method for series polarization model of oiled paper condenser bushing
CN111382531A (en) Current-carrying capacity and temperature simulation method for high-voltage alternating-current submarine cable
Pommerenke et al. Numerical simulation of partial discharge propagation in cable joints using the finite difference time domain method
Maier et al. A high-frequency simulation model of a MV cable joint
CN207148223U (en) A kind of measuring system of large scale electrical power unit wideband impedance operator
Guo et al. An integrated signal and power integrity analysis for signal traces through the parallel planes using hybrid finite-element and finite-difference time-domain techniques
Lu et al. Calculation of Capacitance Matrix of Non-uniform Multi-conductor Transmission Lines based on S-parameters
Krim et al. Two VHDL-AMS-based models of multi-conductor power cables for EMI simulations
Stevanović et al. Multiconductor cable modeling for EMI simulations in power electronics
Ma et al. Predicting Radiated Emissions from a Complex Transportation System Wiring Harness
CN114325167A (en) Microstrip device, measurement system, determination method, device, and storage medium
McGahay Conformal mapping solution for interdigital comb capacitors between ground planes
Salarieh Electromagnetic transient modelling of power transmission line tower and tower-footing grounding system
Hassen et al. A stranded Unshielded Twisted Pair modeling for online fault location using OMTDR-based diagnosis sensor
CN111999564A (en) Calculation method and calculation device of internal dielectric constant value of cable accessories
CN104007327A (en) Method for quickly predicting time domain radiosensitivity of vehicle harness
Berrabah et al. Estimation of distributed and lumped ohmic losses in electrical cables
Niu et al. Research on the Boundary Element Method Algorithm of Capacitance Matrix of Wire Cable with Open Boundary
Ye et al. Time Domain Modeling Method for the Crosstalk Analysis of Multiple Parallel Microstrip Lines

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant