I. INTRODUCTION
Shielding the cable is used to disallow the inevitable coupling of unwanted radiated electromagnetic energy in electrical devices. The interference of electromagnetic decrease depends on two main factors, namely the design and the grounding. Each application needs to consider noise frequency parameters, the frequency of the signal, the length of the cable and the termination methodology of the cable, in some cases, does influence the final result. Improper use of the cable shielding may aggravate the problem via increasing the coupling of the noise. The coupling that exists in-and-out of the shield is due the inadequate quality of the cables. Thus, a prediction of the susceptibility of coaxial shielded cables is required for the optimization of system design.
The prediction of these disturbances, which are usually induced by external fields or lumped sources, is a classical issue which can be dealt with in a variety of ways. Hence, it can be treated in the frequency domain, and therefore the induced responses of shielded cables are solved by using multiconductor transmission line (MTL) theory [Reference Aguet, Ianovici and Lin1–Reference Saih, Rouijaa and Ghammaz3].
SPICE equivalent circuit models for the susceptibility analysis of shielded cables have drawn considerable attention due to their ability to directly use the existing SPICE codes of many components such as inductors, transistors, etc. Orlandi and Antonini [Reference Orlandi4–Reference Saih, Rouijaa and Ghammaz6] developed some SPICE models for the time- and frequency-domain analysis of bulk current injection test on lossless shielded cables, where the injection clamp was considered as a longitudinal lumped current source. These models were utilized to concentrate on the impacts of the shield's grounding, geometrical, and electrical asymmetries and the link's length on the terminal voltages presented in [Reference Antonini, Scogna and Orlandi7]. Some SPICE models were introduced in [Reference Caniggia and Maradei8] to investigate the conducted and radiated invulnerability of shielded cables, the inverse Fourier transform, however, is needed to get the transient results. Xie et al. used some SPICE models to analyze the conducted and radiated susceptibilities of lossless coaxial cables [Reference Xie, Wang, Fan and Liu9], However, the proposed models rely on one assumption that the coaxial cables are lossless. After then, some lossless models have been proposed to study the variation effects of the incident plane wave on a MTL without shields [Reference Mejdoub, Rouijaa and Ghammaz10]. These models can be employed directly in the time and frequency domains.
In this paper, an equivalent circuit model for the analyses of the radiated and conducted susceptibilities of lossy coaxial cables. The model is valid in both time and frequency domains with linear and non-linear loads and effectively brought into the circuit test systems, such as SPICE and ESACAP [Reference Saih, Rouijaa and Ghammaz11, Reference Inzoli and Rouijaa12]. The variation effect of incident plane wave on coaxial cable is studied with the proposed model. The method is validated by comparing results with other methods.
II. DESCRIPTION OF SHIELDED COAXIAL CABLE
A) Model of shielded coaxial cable
The Telegrapher's equations for a shielded coaxial cable above the ground, excited by an incident plane-wave electromagnetic field, as shown in Fig. 2, can be described by:
Outer system (Shield)
Inner System (Wire)
where V s is the shield-to-ground voltage, I s is the current flowing between the external shield and the ground, V w is the wire-to-shield voltage and I w is the current of the wire.
L s , Cs, Rs , and G s are the per-unit-length (p.u.l) inductance, capacitance, resistance, and conductance of the outer system, respectively, while L w , Cw, Rw , and G w are the p.u.l inductance, capacitance, resistance, and conductance matrices of the inner system.
Z t is the transfer impedance, In case of braided shield, the transfer impedance is given by the complex expression [Reference Tesche, Ianoz and Karlsson13, Reference Saih, Rouijaa and Ghammaz14]
where Z d is the diffusion term, and L t is the inductance, which accounts for the field penetrating through the braid apertures. The expression of both Z d and L t in terms of the braid weave parameters can be found in [Reference Tesche, Ianoz and Karlsson13, Reference Saih, Rouijaa and Ghammaz14]. In our application, we used a simplified expression [Reference Caniggia and Maradei8] Z t = R t + jωL t , where R t is the constant p.u.l. transfer resistance of the shield. Figure 1 shows the comparison between the theoretical and asymptotic results of the transfer impedance magnitude [Reference Saih, Rouijaa and Ghammaz3]
V f (z,t) and I f (z,t) are distributed sources that represent external excitation of the transmission line. These source terms can be written solely in terms of the incident electric field using Faraday's law. For coaxial shielded cable, shown in Fig. 2, we have [Reference Paul15]:
where h is the height of the line, and $E_{totoal_z} ^{inc} (h,z,t)$ and $E_{totoal_x} ^{inc} (h,z,t)$ are the horizontal and vertical components of the total field (incident and reflected), respectively, as appeared in Fig. 3.
When the line is situated over a ground plane, as shown in Fig. 4, the total incident field is the sum of the incident field and the ground-reflected field,
These field components are as follows:
The total field is defined by:
where e x , ey , and e z are the components of the incident electric field vector along the x, y, and z axes, and are given by:
The angle θ E characterizes the polarization sort. The polarization is horizontal if θ E is equivalent to zero and vertical if it is equivalent to 90°. The angle θ p decides the rise with respect to the ground. This one is generally called the incident angle. The angle ϕ p gives the propagation direction relative to the axis O as shown in Fig. 2.
The components of the phase constant along those coordinate axes are:
The phase constant is related to the frequency and properties of the medium as:
where v 0 is the phase velocity in the space and the medium is characterized by the permeability μ = μ 0 μ r and permittivity ε = ε0ε r .
B) Equivalent circuit model for conducted immunity: outer system
In the case of conducted immunity the terms I f , V f of equation (1) are nulls. In order to solve equations (1) and (2) we use the ‘discrete line’ model. For this reason, the cable is discretized in the form of cell; the length of each cell is Δz = λ/10.
The solution to equations (1) for each cell is as follow [Reference Orlandi4–Reference Mejdoub, Rouijaa and Ghammaz16]:
where φ s11(Δz), φ s12(Δz), φ s12(Δz) and φ s22(Δz) are the elements of chain parameters given by,
Z cs is the characteristic impedance of the outer system, using the first term of the Taylor series expansion, we obtain
where
and
The characteristic impedance in this case, is presented as a characteristic resistance R c and capacity C f , as shown in Fig. 5
With the same estimation, the constant of propagation gets to be:
Substituting (18) into (17) gives
From equations (22) we can deduce the equivalent circuit of Fig. 5
C) Equivalent circuit model for conducted immunity: inner system
The solution to equation (2) for each cell is as follow:
where the frequency-domain chain parameters are
Here Z cw and γ w are characteristic impedance of the inner system and the constant of propagation, respectively. They are given by
where
and
Substituting (24) into (23) give
From equations (29) we can deduce the equivalent circuit of Fig. 6
D) Equivalent circuit model for radiated immunity of coaxial shielded cable
This is the same representation as the conducted immunity by adding generators ‘forced’ of voltage and current V FT and I FT representing the coupling between the shield and the incident wave, as shown in Fig. 7. Equations (23) remain unchanged, equation (17) becomes:
with
III. APPLICATION RESULTS
A) Conducted susceptibility analysis
The analysis of the conducted immunity is carried out on the coaxial cable over a ground plane as shown in Fig. 8. The length L and the height h of the cable are 1 m and 1 cm, respectively. The shield radius is r s = 2.5 mm, the inner wire radius is r w = 0.25 mm, the internal dielectric constant is ε r = 2.375. The values of the transfer resistance and inductance are: R T = 100 mΩ/m and L T = 0.5 nH/m. The terminal loads between the shield and the ground are R S1 = 1G Ω and R S2 = 154.363 Ω, while the inner terminations are matched R w2 = R w1 = 44.012 Ω.
The lumped current source adopted for the transient analysis is a clock wave of unit amplitude characterized by period T clk = 20 ns, rise and fall time, and duty cycle τ/T clk = 0.5. Figure 9 shows the wire-to-shield voltage at the cable ends, which agree well with the results from different methods.
The wire-to-shield voltage at the cable ends acquired by the proposed model is appeared in Fig. 9 together with the outcomes determined by the FDTD [Reference Roden, Paul and Gedney17], where the “FDTD” implies the finite-difference time-domain solution to the transmission-line equations of the cable, and by the compact circuit model proposed in [Reference Caniggia and Maradei8]. They are in very good agreement with each other.
B) Conducted susceptibility analysis with non-linear loads
To protect an electrical signal transport device against external electromagnetic disturbances, non-linear protection devices (voltage li miters here) are always connected to the devices in parallel. As mentioned above, the proposed model can also be used when the loads are non-linear. The second configuration considered is similar to the first one; a voltage limiter is added to the right termination, which is composed of two anti-parallel diodes, is considered and its effect is studied, as shown in Fig. 10. Two cases have been simulated by using the proposed model, no voltage limiters, and a voltage limiter connected in parallel to the loads R w2. Figure 11 shows the induced voltage at the loads R w2 for the two cases.
From the results of Fig. 11, it is noticed that the voltage limiter can effectively limit the magnitudes of the induced voltages at the loads.
C) Conducted radiated susceptibility analysis of coaxial cable
The analysis of the radiated immunity is carried out on the coaxial cable as shown in Fig. 12. The shield radius R sh and the inner wire radius are 0.25 and 0.108 mm, respectively. The cable's characteristic impedance is Z c = 50 Ω, and the relative permittivity ε r of the internal dielectric filling is 1.77. The value of the transfer impedance is set to R T = 1 Ω/m and L T = 0 H/m. The height h and the length L are 5.25 mm and 1 m, respectively. The internal conductor is adapted with Z a = Z b = 50 Ω. The external field oriented along x- and propagating along z-axis (E x –K z ) has amplitude E = 1 V/m.
The shield is short circuited on the right (Z 2 = 0.5 Ω). At the output of the coaxial we recover the current into dBA, which is matched with the canonic results published by Smith [Reference Smith18], as shown in Fig. 13.
When the shield is open on the left (Z 1 = 1 GΩ) the lines resonates at f = n × (3 × 108 4/λ), n = 1,3,5,… (f 1 = 75 MHz, f 2 = 225 MHz, f 3 = 375 MHz, etc.), as shown in Fig. 13. To drastically diminish the coupling to internal wire, a two-side grounded configuration for the shield must be utilized.
D) Radiated susceptibility analysis of coaxial cable: variation effect of the outer loads and incident plane wave
The setup utilized for the radiated susceptibility analysis is appeared in Fig. 12. The height h and the length L are 5.25 mm and 1 m, separately. The shield radius R sh and the inner wire radius are 0.25 and 0.0716 mm, separately, and the relative permittivity ε r of the inner filling dielectric is 2.25. The loads Z a and Z b between the inner wire and the shield at the two terminations are 50 Ω. The value of the transfer impedance is set to R T = 0.01 Ω/m and L t = 1 nH/m. The incident field is modeled by the double exponential pulse E(t) = kE 0 [exp (−At) − exp (−Bt)], where E 0 = 50 kV/m, k = 1.3, A = 6 × 108s−1, and B = 4 × 107s−1.
Figure 14 shows the induced voltages at the internal loads R1 and R2 with the incident wave when the outer terminal loads (Zout = Z 1 = Z 2) have different values. Four different values of the outer loads have been considered. At the point when the line is ended with short circuit at the far-end and the near-end (Zout = Z 1 = Z 2 = 0 Ω), the induced voltages are largest. However, this is the most common case in reality.
The shield configuration with both sides is matched (Z c = Z 1 = Z 2 = 244.5 Ω), the induced voltages are not the smallest but decay most quickly. The other values of the outer terminators result in reflections caused by an impedance mismatch at the ends of the outer system. Even though the inner transmission-line system is matched at both ends, the reflected voltage and current of the outer system can couple into the inner system and generate interferences.
Figure 15 shows the transient voltage responses at the internal loads R1 and R2 obtained by the proposed model due to the external electromagnetic fields when the outer termination is matched and the incident wave, corresponding to ϕ p = 0°, &hbox;ϕ p = −30°, ϕ p = −60° and ϕ p = −90°, with θ E = −90°and the polarization angle θ p = −90°, It is obvious that the induced voltage changes with decreasing the value of ϕ p . The induced voltage reaches its lowest among two cases of ϕ p = 0° and ϕ p = −30°. This is reasonable because ϕ p = 0° means that the incident electric field is vertical to the coaxial cable, which results in the minimum of distributed voltage source along the shield, and then induce the minimal voltage.
IV. CONCLUSION
In this paper, we introduced an efficient coaxial shielded cables model, for conducted and radiated coupling, developed for circuit applications. Its advantage then makes it possible to study the frequency- and time domains with linear and non-linear loads, respectively. A detailed study of coaxial shielded cables has been carried out. The suggested technique is verified by comparing the circuit simulator results with the solutions derived by the other methods have revealed a satisfactory accuracy. The responses of the configuration to the external electromagnetic fields changing with the outer termination have been studied with the proposed model. The results show that the induced voltages are the largest when the shielding enclosure is connected to the ground directly, which is the most common case actually, and adding a resistance between the enclosure and the ground can reduce the induced voltages effectively. The analysis also shows that the induced voltage changes with decreasing the value of ϕ p . The induced voltage reaches its lowest when ϕ p = 0°. It is easy to extend the models to coaxial shield cables excited by a non-uniform incident wave, which require totally a discretization. This inquiry will be examined in a further study.
Mohamed Saih received the Engineer Diploma in Electrical Systems and Telecommunications from Cadi Ayyad University, Marrakech, Morocco, where he is currently working toward the Ph.D. degree at the Department of Applied Physics, Electrical Systems and Telecommunications Laboratory, Cadi Ayyad University of Marrakech, Morocco. His research interests include electromagnetic compatibility, multiconductor transmission lines, numerical electromagnetic methods, and antenna designs.
Hicham Rouijaa is a Professor of Physics, attached to Cadi Ayyad University, Marrakesh, Morocco. He is obtained his Ph.D. Thesis on \” Modeling of Multiconductor Transmission Lines using Pade approximant method: Circuit model \”, in 2004, from Aix Marseille University – France. He is an associate member of Electrical Systems and Telecommunications Laboratory LSET at the Cadi Ayyad University. His current research interests concern electromagnetic compatibility and multiconductor transmission lines.
Abdelilah Ghammaz received the Doctor of Electronic degree from the National polytechnic Institut (ENSEEIHT) of Toulouse, France, in 1993. In 1994, he went back to Cadi Ayyad University of Marrakech – Morroco. Since 2003, he has been a Professor at the Faculty of Sciences and technology, Marrakech, Morroco. He is a member of ELectrical Systems and Telecommunications Laboratory LSET at the Cadi Ayyad University. His research interests in the field of electromagnetic compatibility, multiconductor transmission lines, telecommunications, and antennas.