![]() |
|
|
|
|
An introduction to the art of Solenoid inductance calculation By David W Knight
|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
The following easy-to-use routine for calculating solenoid inductance by the Rosa-Nagaoka method has been added to the macro library in the spreadsheet Lcalcs.ods . The three functions called are also in the macro library and are discussed in the solenoid incuctance calculation article. Note that the calculation is for external partial inductance only. Internal inductance and lead inductance should be added to the result |
Function
LxRosa(byval h as double, D as double, N as double, dw as double, p as
double) 'External inductance in henrys of a round-wire solenoid, using the Nagaoka-Rosa method. 'D W Knight. v 1.00, 2016-01-11 'Calls functions W82W(), KMGO() and Rosaks90() 'h = coil length / m 'D = coil diameter / m 'N = number of turns 'dw = wire diameter / m 'p = winding pitch / m Dim kNagaoka as double, Lsheet as double, Rosacorr as double kNagaoka = W82W(D/h) Lsheet = 1e-7*pi()*pi()*D*D*N*N*kNagaoka/h Rosacorr = KMGO(N) + Rosaks90( p/dw , p/D , 0) LxRosa =Lsheet -2E-7*pi()*D*N*Rosacorr end function |
Material not yet
transferred. The text below is in the process of being updated and extended. The information will at some point be transferred to the main solenoid inductance article above. |
B1. Solenoid inductance
calculation vs. measurement: Using the techniques outlined in the previous sections, it is possible to calculate the low-frequency external inductance of solenoids to an accuracy of a few parts in 1000 using a micrometer or a good set of engineer's callipers and a personal computer. If a scientific calculator is used, the computation can be simplified by using one of the Wheeler unrestricted formulae with only a small loss of accuracy. In either case, the outcome is usually better than can be obtained using general-purpose measuring equipment. All solenoid inductance calculations are dependent on the effective current-sheet diameter (D) or radius (r=D/2). In section 3 it was noted that at low frequencies D=D0, where D0≈Da and Da is the coil diameter measured from wire-centre to wire-centre, whereas at very high frequencies D=D∞, where D∞ is less than D0 but always greater than the diameter at the inside of the conducting cylinder (i.e., D0-d). Consequently, the measured radio-frequency inductance of a coil (substantially below the SRF at least) will always lie between the low frequency inductance L0, which is the value calculated using D=D0; and the minimum inductance L∞, which is the value which obtains when D=D∞ (presuming that we can determine D∞). The inductance of a solenoid without a magnetic core, presuming that the wire diameter is << the solenoid diameter, is given within a few parts per thousand by:
where D=2r, kL is Nagaoka's coefficient , ks(e) Rosa's round-wire self inductance correction, for external inductance only, and km is Rosa's round-wire mutual inductance correction.. The internal inductance Li , generally < 1% of the total inductance, can be calculated to within 0.016% using the Li-PACAML approximation, see Zint.pdf. (also function Lintern() described on the Zint index page) From the previous discussion, the inductance of a solenoid, neglecting the proximity effect, is:
where kL(0) is Nagaoka's coefficient calculated using ![]() The worst-case minimum inductance, due to reduction in the effective diameter is:
where kL(∞) is Nagaoka's coefficient calculated using ![]() |
Example B1a R G Medhurst, on page 42 of his 1947 paper "H.F Resistance and Self-Capacitance of Single-Layer Solenoids", gives accurate dimensions of a test-coil (#31) used in a study of Q and self-capacitance. The coil is wound on a grooved former using 40 turns of 20SWG bare wire (d=0.9144 mm), with a diameter of 51.9 mm measured at the outside of the winding. The length of the coil is given as 70.1 mm, but the winding pitch is stated to be: p=(70.1-d)/(N-1) from which we may deduce that the measurement was made to the outside of the wire on the side away from the terminations, (as shown in the diagram on the right). Hence the stated length is not the same as the equivalent current sheet length, which is defined as: |
![]() |
![]() Thus it appears that Medhurst uses a definition of coil length that differs from that intended by Rosa and Grover (as discussed in section 2) but is an understandable misinterpretation of the explanation given [Sci169, p119]. Ambiguities of this type can introduce errors of several % in inductance calculations, and there is no use in claiming accuracy in parts per 1000 unless care is taken to avoid them. In this case the equivalent current-sheet length is: ![]() Using this value, with Nagaoka's coefficient kL calculated using Lundin's formula , we have:
The round-wire correction ks(e) is: ks(e) = (3/2) - ln(2p/d)
and km =0.3142. Medhurst measured the inductance of this coil by averaging readings taken over the range 780 to 860 kHz. Taking the average frequency as 820 kHz, the internal inductance is 101.8 nH. The results of the calculation are shown below, where the 'greater than' symbol (>) is used for the L∞ value because the coil inside diameter (D0 - d) was used instead of the slightly larger but undefined D∞. L0 = (μ0 π D0² N² kL(0) /4 ![]() = 43.7331 - 0.5873 + 0.1018 μH
L∞ = (μ0 π D∞² N² kL(∞) /4 ![]() > 42.3694 - 0.5768 + 0.1018 μH
Hence we expect the measured inductance to lie in the range: L = 43.248, +0, -1.354 μH |
Medhurst gives nine inductance measurements for this coil, taken at 10 kHz intervals over the range 780 to 860 kHz and corrected for self capacitance. These measurements are scattered about the mean value, whereas if the inductance were changing significantly with frequency in this region we would expect them to diminish progressively as the frequency increased. Hence we can infer that deviations from the mean are principally due to experimental error, and we can use this information to estimate the standard deviation of the mean value. The measurements are tabulated below: |
|
/ kHz |
Lk / μH |
Lk - Lmean |
(Lk - Lmean)² |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Lmean = 42.3744 |
|
|
The mean (average) of the
measurements is given by:
Where 'Σ' (capital Sigma) means "the sum of" (over the range specified below and above), Lk is the 'kth' measurement of L, and in this case the number of measurements n=9. The square of the standard deviation (known as the variance) can be estimated from:
where we divide the sum of the squares of the deviations by n-1 rather than n because the use of the sample mean rather than the population mean (i.e.the unknown 'true' value of L) removes a degree-of-freedom from the data [see Data Analysis]. Thus, for Medhurst's measurements: σ = √(0.02482222 / 8) = 0.0557 Hence:
Comparing this with our calculations (43.248 > L > 41.894 [μH] ), we can observe that at 820 ±40 kHz, the inductance of the coil is lower than that predicted from physical dimensions, but within the range allowed by variation of effective diameter. It is tempting to think that the high-frequency drop in inductance has been quantified by this comparison (and possibly it has), but such is not necessarily the case. The standard deviation calculated from the scatter in the data is a measure of precision not accuracy, and unknown systematic errors, or simple uncertainties in the coil dimensions, prevent us from drawing firm conclusions. |
Example B1b Once the principle has been established, accurate inductance calculation does not have to be a laborious process.We might note from the earlier discussion however; that the available literature on the subject contains mistakes, with the consequence that existing computer programs and models cannot be trusted without verification of the numbers that they produce. It is a relatively straightforward matter however, to set up a spreadsheet to do the calculation, and once the entered formulae have been verified, new rows can be added and old ones deleted at will. A suitable calculation template is provided by the Open-Document Spreadsheet file: L_calcs.ods (sheet 2), which can be downloaded and amended as required. The spreadsheet as provided is filled with example inductance calculations for seven coils described in the academic literature. Medhurst's coil #31 (M31) from the previous example is there , along with another coil (M32) described in the same paper (these are the only test coils for which Medhurst gave full details). There are also: a coil described in reference [10] (GKMR1), and four coils described in reference [11] (MKG1-4). Interest in the last five coils lies in the fact that the inductance measurements are accompanied by calculations using a method different to the one described here. The results are summarised in the table below: |
|
|
|
|
|
|
|
|
![]() |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L0 /μH |
|
|
[10] |
[11] |
[11] |
[11] |
[11] |
The main point to note is that for
coils MKG1-4, low-frequency (L0)
inductances calculated using the modified National Bureau of Standards
(NBS) method (i.e., the method used here) are in exact agreement with
the method used in ref. [11],
and that for a measurement frequency of 1 kHz, the internal inductance
is the same as the DC value. This provides a vindication, not of the
NBS method which is not in doubt, but of the formula verification
procedures applied here. The calculation of ref. [10] however, is not
in agreement with the NBS method, despite the statement that it was
carried out using the same method as in ref. [11]. Although very
close to the measured value, the calculation of ref. [10] does not
correspond to L0 or L0-Li,
and appears for some reason to be in error. What we find here is that
the measured inductance is sandwiched between our L0
and L∞>
values as in the case of coils M31 and M32, and that the very low p/d
ratio of 1.02 (and probable use of copper tubing rather than wire) has
given rise to an early onset of the high-frequency regime. The results
for coils MKG1-4 show nothing more onerous than a measurement standard
deviation of around 1%, and that the HF regime is not evident for these
coils at a frequency of 1 kHz. If there is anything to be inferred from these comparisons between measurement and calculation, it is that the skin effect and the proximity (effective diameter) effect are linked. All of the coils measured at 1 kHz have an internal inductance component that is the same as the DC value, and show no sign of the diameter effect even though some have low p/d ratios. All of the coils measured at frequencies above 1 kHz have an internal inductance component substantially less than the DC value, and the inductances lie between the L0 and L∞ values. Logically, this is to be expected, because the proximity effect modifies the current distribution within the wire and therefore changes the internal inductance. The proximity effect however, is not simply a perturbation of the internal inductance, because the total variation it causes is greater than the internal inductance. >>>>>> Calculation methods have been refined since 2008. More accurate formulae for effective current sheet diameter have been derived. New spreadsheet: Lcalcs.ods (sheet 1, still a work in progress). >>>>>> |
B2. Additional sources
of deviation: Points to note when comparing inductance measurements with calculations: ![]() ![]() ![]() ![]() ![]() |
[10]
"Stray Capacitances of Single-Layer Solenoid Air-Core
Inductors", G. Grandi, M K Kazimierczuk, A Massarini, U
Reggiani. IEEE Transactions on Industry Applications, Vol 35, No. 5,
Sept/Oct 1999, p1162-1168. Available from: Classic Tesla at time of writing. www.classictesla.com/download/ia99.pdf . [11] "Lumped Parameter Models for Single- and Multiple-Layer Inductors", A Massarini, M K Kazimierczuk, G Grandi. Proc PESC '96, June 1996, p295-301. Available from: Classic Tesla at time of writing. www.classictesla.com/download/pesc96.pdf . |
![]() |
|
|
|
|