CN105844092A - Prediction and assessment method for two-lead transmission line terminal response - Google Patents

Abstract

The invention discloses a prediction and assessment method for a two-lead transmission line terminal response in a non-analytic distribution source. The method comprises the steps of measuring physical parameters of leads; geometrically discretizing a transmission line; measuring geometric parameters of the transmission line; measuring electric parameters of the transmission line; and finally making a calculation to obtain an induced voltage and an induced current of a parallel or non-parallel two-lead transmission line terminal. In an assessment process, a specific analytic expression of an additional excitation source does not need to be known and only discretized data, distributed along the leads, of an excitation field needs to be known. According to the method, computing errors caused by reasons such as energy loss and the like in an irradiation process of the distribution source do not need to be considered, the transmission line is slightly transformed, the measurement data amount is relatively small, and the accuracy of a measurement result is relatively high; in a measurement process, other additional circuits do not need to be carried, so that the errors caused by additional devices are relatively small; and the method is suitable for various additional electromagnetic field irradiation and is simple in operation and easy to realize.

Description

A kind of prediction and evaluation method of two-conductor line transmission-wire terminal response

Technical field

The present invention relates to the appraisal procedure of Wire termination response, particularly to the prediction and evaluation method of a kind of two-conductor line transmission-wire terminal response.

Background technology

Along with the high speed development of hyundai electronics science and technology, in many electronic equipments, the integrated level of circuit is more and more higher, and the volume of equipment is more and more less, Internal structure also becomes increasingly complex, and the electromagnetic signal strength frequency range used is different, sharply increasing of simultaneous external world radiation interference source category, All make the electromagnetic environment residing for the work of hyundai electronics information equipment the most complicated, often show as spatial density very big, time domain is difficult to estimate, frequently On territory, wavelength band is extremely wide, and type becomes more diverse.

In substantial amounts of high speed integrated circuit and microwave transmission, transmission line is the interconnection line form that a kind of distribution is the most basic.Say accurately, Transmission line is summarized as in electrical engineering and electronic engineering field, makes electromagnetic wave make its medium propagated along certain set direction.For low-frequency transmission For line, the impact of signal is negligible by transmission line itself substantially, but for high frequency transmission line, along with the rising of signal frequency, Transmitting procedure will be affected greatly by the characteristic of transmission line itself.Meanwhile, when the interference of electromagnetic field that transmission line is applied by outside, will produce With induced voltage and the faradic current of exciting field change, the change of these load responses caused by driving source can produce shadow to the transmission usefulness of wire Ring, pull one hair and move the whole body, and then the stability to whole electronic system, functional generation significant impact.First and last, this buterfly effect meeting Have a strong impact on electronic system operating, in some instances it may even be possible to cause whole system paralysis damage, endanger personal safety.

In recent years, numerous researchers have done a lot of work for transmission line coupling problem, but are in most cases directed to the driving source resolved Model, namely need to try to achieve transmission-wire terminal response in the case of known extrinsic motivated source resolution expression formula.But in practical engineering application In the driving source model that relates to the most nonanalytic, often can only obtain the driving source discretization distribution situation along wire, the most urgently propose one The forecast assessment of the two-conductor line transmission-wire terminal response that method is capable of under non-analytic expression distributed source.

Summary of the invention

For overcoming deficiency of the prior art, the present invention proposes the forecast assessment side of the two-conductor line transmission-wire terminal response under a kind of non-analytic expression distributed source Method.The present invention " prediction and evaluation method of a kind of two-conductor line transmission-wire terminal response ", comprises the steps:

Step one, two-conductor line transmission line is carried out the measurement of physical parameter, including: conductor length L, wire radius a, the left end point of two wires Spacing d, wire angle 2 θ, terminator Z₁, Z₂, exciting field polarizing angle α, exciting field angle of incidence φ, ψ；The length phase of described two wires Deng, equal diameters；

Step 2, in OXY coordinate plane, the coordinate of the left and right end points that will be located in the wire 5 above X-axis be set to (0, d/2), (Lcos θ, D/2+Lsin θ), the coordinate of the left and right end points that will be located in the wire 6 below X-axis is set to (0 ,-d/2), (Lcos θ ,-d/2-Lsin θ), two wires About X-axis specular；Along wire direction, wire 5 or wire 6 are split, wire 5 or wire 6 are divided into n line segment, line segment Label is designated as m, m=1,2 ..., n；Measure the length Δ x of each line segment projection in X-axis_m, the projection in X-axis of each line segment Coordinate (the x at midpoint^c _m, y^c _m, z^c _m), and the midpoint of each line segment with it about distance d of the mirror image of X-axis_m；, wherein the value of subscript m is 1,2 ..., n, and assume d₁=d.

Step 3, measure the direction of propagation of two-conductor line transmission line extrinsic motivated field, calculate all kinds of transmission line electrical quantity: the long inductance of unitThe long electric capacity of unitThe long impedance Z of unit_m=R_m+jωL_m, unit long admittance Y_m=G_m+jωC_m, TransmissionCharacteristic impedanceWherein, R_mFor the long resistance of unit, G_mFor the long conductance of unit, ε is dielectric Constant, μ is pcrmeability, and j is imaginary unit, and ω is wave frequency variable, subscript m=1,2 ..., n；

Step 4, the induced voltage V (0), V (L) of the two-conductor line transmission-wire terminal utilized under formula (1), (2) calculating extrinsic motivated field radiation parameter With faradic current I (0), I (L):

$\left(\begin{array}{c}V\left(0\right)\\ V\left(L\right)\end{array}\right)=\left(\begin{array}{cc}1+{\mathrm{\ρ}}_1& 0\\ 0& 1+{\mathrm{\ρ}}_2\end{array}\right){\left(\begin{array}{cc}-{\mathrm{\ρ}}_1& e^{\mathrm{\γ}L}\\ e^{\mathrm{\γ}L}& -{\mathrm{\ρ}}_2\end{array}\right)}^{-1}\left(\begin{array}{c}{S}_1\\ S_2\end{array}\right)---\left(1\right)$

$\left(\begin{array}{c}I\left(0\right)\\ I\left(L\right)\end{array}\right)=\frac{1}{Z_c}\left(\begin{array}{cc}1+{\mathrm{\ρ}}_1& 0\\ 0& 1+{\mathrm{\ρ}}_2\end{array}\right){\left(\begin{array}{cc}-{\mathrm{\ρ}}_1& e^{\mathrm{\γ}L}\\ e^{\mathrm{\γ}L}& -{\mathrm{\ρ}}_2\end{array}\right)}^{-1}\left(\begin{array}{c}{S}_1\\ S_2\end{array}\right)---\left(2\right)$

Wherein

$\left(\begin{array}{c}{S}_1\\ S_2\end{array}\right)=\left(\begin{array}{c}\frac{1}{2}{\∫}_0^L{e}^{\mathrm{\γ}x}{V}_{s2}\left(x\right)dx-\frac{V_1}{2}+\frac{V_2}{2}{e}^{\mathrm{\γ}L}\\ -\frac{1}{2}{\∫}_0^L{e}^{\mathrm{\γ}(L-x)}{V}_{s2}\left(x\right)dx-\frac{V_2}{2}+\frac{V_1}{2}{e}^{\mathrm{\γ}L}\end{array}\right)$

$V_1=-E_z^{inc}(0,0,0)d_1$

$V_2=-E_z^{inc}(L_x,0,0)d_n$

WhereinFor coordinate (x, y, z) projection in X-axis of place's incident electric fields,For coordinate, (x, y, z) place's incident electric fields exists Projection on Z axis, L_xFor wire in the length of the projection of X-axis,For m-th line segment at the x coordinate of the central point of the projection of X-axis,Respectively load Z₁And Z₂The reflection coefficient at place.

According to calculated induced voltage and faradic current, secure threshold to transmission line assessment can be predicted further.

The present invention chooses the induced voltage of target two-conductor line any terminal with electric current to the terminal field circuit method effect analyzing under the irradiation of extrinsic motivated field, carries Going out a kind of measuring method calculating two-conductor line transmission-wire terminal inductive voltage and current, the two-conductor line that can effectively measuring be caused by additional electromagnetic field is eventually The inductive voltage and current of end.The operation of these computational methods is relatively easy easy, can reuse in a large number, can be that certain required precision is interior by additional The induced voltage of the Wire termination that Electromagnetic Field Irradiation causes provides a kind of effective means with faradic measurement, for analyzing additional electromagnetism for transmission The coupling effect of line and terminal unit thereof provides important parameter.

Compared with the prior art, present invention have the advantage that

(1) present invention requires no knowledge about extrinsic motivated field concrete analytical expressions, it is only necessary to acquisition driving source along the discretization data of arrangement of conductors is Can calculate, more closing to reality engineer applied；And it is generalized to the situation of non-parallel two-conductor line transmission line；

(2) present invention is based on transmission line rationale, directly obtains Wire termination sensing by extrinsic motivated field along the discretization data of arrangement of conductors Voltage and faradic current, thus the calculating error caused by reasons such as energy losses in distributed source irradiation process need not be considered, transmission line is changed Making less, measurement data is less but measurement result accuracy is higher；

(3) induced voltage of Wire termination is directly calculated by the present invention with faradic current, it is not necessary to carry other adjunct circuits, therefore by attached The error that oil (gas) filling device introduces is less, and test result only relies upon extrinsic motivated field along the discretization data of arrangement of conductors and wire relevant physical parameter, The accuracy ensureing result of calculation can be maximized；

(4) computational methods in the present invention are applicable to various types of additional electromagnetic field irradiation, it is adaptable to each frequency range, each period, have one Fixed universality；

(5) computation complexity of the present invention is little, low cost, repeatable high, simple to operate, it is easy to accomplish.

Accompanying drawing explanation

In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to required in embodiment or description of the prior art The accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only some embodiments of the present invention, common for this area From the point of view of technical staff, on the premise of not paying creative work, it is also possible to obtain other accompanying drawing according to these accompanying drawings.

Fig. 1 is the schematic diagram in the present invention split target two-conductor line transmission line, wherein 1,2 loads for Wire termination, and 3 is virtual dividing line, 4 is the central point of m-th line segment projection in X-axis, and 5 is the wire above X-axis, and 6 is the wire below X-axis；

Fig. 2 is the electric current frequency domain response schematic diagram that in the present embodiment, two-conductor line transmission line model terminal 2 is under different conductor angle；

Fig. 3 is the voltage frequency domain response schematic diagram that in the present embodiment, two-conductor line transmission line model terminal 2 is under different conductor angle；

Specific embodiments

Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that institute The embodiment described is only a part of embodiment of the present invention rather than whole embodiments.Based on the embodiment in the present invention, the common skill in this area The every other embodiment that art personnel are obtained under not making creative work premise, broadly falls into the scope of protection of the invention.

Induced voltage at two-conductor line transmission-wire terminal under additional plane wave exciting field irradiation is calculated by the present embodiment with faradic current numerical value, Obtain the response numerical value of frequency domain lower wire end respectively.

Step one: selected two-conductor line transmission line carries out the measurement of physical parameter, and design parameter is as follows:

Driving source:

$E^{inc}\left(\mathrm{\ω}\right)=\frac{1}{2\mathrm{\π}}(\frac{1}{4\×{10}^6+\mathrm{\ω}i}-\frac{1}{4.76\×{10}^8+\mathrm{\ω}i})$

Conductor length L=30 rice (m), wire left end point distance d=0.2 rice (m), wire radius a=0.0015 rice (m), Wire termination loads Z₁=Z₂=293 Ω, wire angle 2 θ=π/1800, exciting field polarizing angle α=0, exciting field angle of incidence φ=0, ψ=π/3.

Step 2: target two-conductor line transmission line is modeled, as it is shown in figure 1, wire is carried out geometry sliding-model control.By each wire by identical Away from being divided into n=500 section, a length of Δ x of m-th line segment projection in x-axis_m=L_x/ n, (wherein L_xFor wire projection in x-axis Length, m=1,2 ..., n), the center point coordinate that m-th line segment projects in x-axis is The center of m-th line segment and its about the distance between the mirror image of x-axis beAt Δ x_mIn the case of less, can approximate D is thought on ground₁=d.

Step 3: calculate linear electrical parameter, specific as follows: m section unit long resistance R_m=0, m section unit long conductance G_m=0, m The long inductance of Duan DanweiThe long electric capacity of m section unitThe long impedance Z of m section unit_m=R_m+jωL_m, the M section unit long admittance Y_m=G_m+jωC_m, m section transmissionCharacteristic impedanceWherein, dielectric Constant ε=(10^-9)/36 π (F/m), π × 10, magnetic permeability μ=4^-7(H/m), subscript m=1,2 ..., n；.

Step 4: utilize formula (3), (4) calculate distributed source frequency domain:

$\begin{array}{c}{V}_{s2}\left(x_m\right)=E_x^{inc}(x_m^c,0,\frac{d_m}{2})-E_x^{inc}(x_m^c,0,-\frac{d_m}{2})\\ =E^{inc}(\cos \mathrm{\α}\cos \mathrm{\φ}\sin \mathrm{\ψ}+\sin \mathrm{\α}\sin \mathrm{\φ})\\ (e^{jkd\sin \mathrm{\ψ}}-1)e^{-{\mathrm{jkx}}_m^c\cos \mathrm{\ψ}\cos \mathrm{\φ}}\end{array}.---\left(3\right)$

$V_1=-E_z^{inc}(0,0,0)d_1=-E_0{d}_{1}cos\mathrm{\ψ}cos\mathrm{\α}---\left(4\right)$

$V_2=-E_z^{inc}(L_x,0,0)d_n=-E_0{d}_n{\mathrm{cos\ψcos\αe}}^{-{\mathrm{jkL}}_{x}cos\mathrm{\ψ}cos\mathrm{\φ}}=V_1{e}^{-{\mathrm{jkL}}_{x}cos\mathrm{\ψ}cos\mathrm{\φ}}---\left(5\right)$

Wherein E₀It is electric field amplitude, E₀=1/2 π, wave number

Above-mentioned distributed source expression formula (3)～(5) are substituted in BLT equation source item (6) formula:

$\left(\begin{array}{c}{S}_1\\ S_2\end{array}\right)=\left(\begin{array}{c}\frac{1}{2}{\∫}_0^L{e}^{\mathrm{\γ}x}{V}_{s2}\left(x\right)dx-\frac{V_1}{2}+\frac{V_2}{2}{e}^{\mathrm{\γ}L}\\ -\frac{1}{2}{\∫}_0^L{e}^{\mathrm{\γ}(L-x)}{V}_{s2}\left(x\right)dx-\frac{V_2}{2}+\frac{V_1}{2}{e}^{\mathrm{\γ}L}\end{array}\right)---\left(6\right)$

Again source item (6) formula is substituted into BLT equation expression formula (7), (8), it is possible to measure frequency domain induced voltage V (0) of end, V (L) With faradic current I (0), I (L).

$\left(\begin{array}{c}V\left(0\right)\\ V\left(L\right)\end{array}\right)=\left(\begin{array}{cc}1+{\mathrm{\ρ}}_1& 0\\ 0& 1+{\mathrm{\ρ}}_2\end{array}\right){\left(\begin{array}{cc}-{\mathrm{\ρ}}_1& e^{\mathrm{\γ}L}\\ e^{\mathrm{\γ}L}& -{\mathrm{\ρ}}_2\end{array}\right)}^{-1}\left(\begin{array}{c}{S}_1\\ S_2\end{array}\right)---\left(7\right)$

$\left(\begin{array}{c}I\left(0\right)\\ I\left(L\right)\end{array}\right)=\frac{1}{Z_c}\left(\begin{array}{cc}1+{\mathrm{\ρ}}_1& 0\\ 0& 1+{\mathrm{\ρ}}_2\end{array}\right){\left(\begin{array}{cc}-{\mathrm{\ρ}}_1& e^{\mathrm{\γ}L}\\ e^{\mathrm{\γ}L}& -{\mathrm{\ρ}}_2\end{array}\right)}^{-1}\left(\begin{array}{c}{S}_1\\ S_2\end{array}\right)---\left(8\right)$

The measurement result of the faradic current I (L) at terminal impedance 2 and induced voltage V (L) is the most as shown in Figure 2 and Figure 3.In Fig. 2, long dotted line represents Article two, wire is not parallel, and when its angle is π/1800, the measured value of frequency domain faradic current I (L), short dash line represents when two wires are parallel, frequency domain sense The measured value of induced current I (L), fine line represents secure threshold.In Fig. 3, long dotted line represents that two wires are not parallel, when its angle is π/1800, frequently The measured value in territory induced voltage V (L), short dash line represents that, when two wires are parallel, the measured value of frequency domain induced voltage V (L), fine line represents safety threshold Value.

By being calculated two-conductor line transmission-wire terminal induced voltage and faradic current, can further the secure threshold of transmission line system be predicted Assessment, the secure threshold of transmission line depends primarily on maximum breakdown voltage threshold size on line.Therefore terminal induction voltage and faradic is being obtained On the basis of, it is possible to assessment is predicted with regard to concrete transmission line.Assuming that the breakdown voltage in given embodiment is 6 × 10^-10V, can lead to Cross Fig. 3 and show that the frequency range exceeding transmission line secure threshold is 0～10⁸Hz, namely this can be caused to have the non-parallel biography of consumption at the driving source of this frequency range Induced voltage at defeated line terminal 2 is more than the secure threshold of wire itself, it is possible to cause transmission line system unstable.

The above, the only detailed description of the invention of the present invention, but protection scope of the present invention is not limited thereto, any it is familiar with the art Technical staff in the technical scope that the invention discloses, the change that can readily occur in or replacement, all should contain within protection scope of the present invention.Cause This, protection scope of the present invention should be as the criterion with the protection domain of claims.

Claims (1)

1. the prediction and evaluation method of a two-conductor line transmission-wire terminal response, it is characterised in that described method comprises the steps:

Step one, two-conductor line transmission line is carried out the measurement of physical parameter, including: conductor length L, wire radius a, the left end point of two wires Spacing d, wire angle 2 θ, terminator Z₁, Z₂, exciting field polarizing angle α, exciting field angle of incidence φ, ψ；The length phase of described two wires Deng, equal diameters；

Step 2, in OXY coordinate plane, the coordinate of the left and right end points that will be located in the wire 5 above X-axis be set to (0, d/2), (Lcos θ, D/2+Lsin θ), the coordinate of the left and right end points that will be located in the wire 6 below X-axis is set to (0 ,-d/2), (Lcos θ ,-d/2-Lsin θ), two wires About X-axis specular；Along wire direction, wire 5 or wire 6 are split, wire 5 or wire 6 are divided into n line segment, line segment Label is designated as m, m=1,2 ..., n；Measure the length Δ x of each line segment projection in X-axis_m, the projection in X-axis of each line segment Coordinate (the x at midpoint^c _m, y^c _m, z^c _m), and the midpoint of each line segment with it about distance d of the mirror image of X-axis_m；, wherein the value of subscript m is 1,2 ..., n, and assume d₁=d.

Step 3, measure the direction of propagation of two-conductor line transmission line extrinsic motivated field, calculate all kinds of transmission line electrical quantity: the long inductance of unitThe long electric capacity of unitThe long impedance Z of unit_m=R_m+jωL_m, unit long admittance Y_m=G_m+jωC_m, TransmissionCharacteristic impedanceWherein, R_mFor the long resistance of unit, G_mFor the long conductance of unit, ε is dielectric Constant, μ is pcrmeability, and j is imaginary unit, and ω is wave frequency variable, subscript m=1,2 ..., n；

Step 4, the induced voltage V (0), V (L) of the two-conductor line transmission-wire terminal utilized under formula (1), (2) calculating extrinsic motivated field radiation parameter With faradic current I (0), I (L):

$\left(\begin{array}{c}V\left(0\right)\\ V\left(L\right)\end{array}\right)=\left(\begin{array}{cc}1+{\mathrm{\ρ}}_1& 0\\ 0& 1+{\mathrm{\ρ}}_2\end{array}\right){\left(\begin{array}{cc}-{\mathrm{\ρ}}_1& e^{\mathrm{\γ}L}\\ e^{\mathrm{\γ}L}& -{\mathrm{\ρ}}_2\end{array}\right)}^{-1}\left(\begin{array}{c}{S}_1\\ S_2\end{array}\right)---\left(1\right)$

$\left(\begin{array}{c}I\left(0\right)\\ I\left(L\right)\end{array}\right)=\frac{1}{Z_c}\left(\begin{array}{cc}1+{\mathrm{\ρ}}_1& 0\\ 0& 1+{\mathrm{\ρ}}_2\end{array}\right){\left(\begin{array}{cc}-{\mathrm{\ρ}}_1& e^{\mathrm{\γ}L}\\ e^{\mathrm{\γ}L}& -{\mathrm{\ρ}}_2\end{array}\right)}^{-1}\left(\begin{array}{c}{S}_1\\ S_2\end{array}\right)---\left(2\right)$

Wherein

$\left(\begin{array}{c}{S}_1\\ S_2\end{array}\right)=\left(\begin{array}{c}\frac{1}{2}{\∫}_0^L{e}^{\mathrm{\γ}x}{V}_{s2}\left(x\right)dx-\frac{V_1}{2}+\frac{V_2}{2}{e}^{\mathrm{\γ}L}\\ -\frac{1}{2}{\∫}_0^L{e}^{\mathrm{\γ}(L-x)}{V}_{s2}\left(x\right)dx-\frac{V_2}{2}+\frac{V_1}{2}{e}^{\mathrm{\γ}L}\end{array}\right)$

$V_{s2}\left(x_m\right)=E_x^{inc}(x_m^c,0,\frac{d_m}{2})-E_x^{inc}(x_m^c,0,-\frac{d_m}{2})$

$V_1=-E_z^{inc}(0,0,0)d_1$

$V_2=-E_z^{inc}(L_x,0,0)d_n$

WhereinFor coordinate (x, y, z) projection in X-axis of place's incident electric fields,For coordinate, (x, y, z) place's incident electric fields exists Projection on Z axis, L_xFor wire in the length of the projection of X-axis,For m-th line segment at the x coordinate of the central point of the projection of X-axis,Respectively load Z₁And Z₂The reflection coefficient at place.