Open access peer-reviewed chapter

The Mechanical Behavior of the Cable-in-Conduit Conductor in the ITER Project

Written By

Donghua Yue, Xingyi Zhang and Youhe Zhou

Submitted: 11 September 2018 Reviewed: 31 October 2018 Published: 31 December 2018

DOI: 10.5772/intechopen.82349

From the Edited Volume

Nuclear Fusion - One Noble Goal and a Variety of Scientific and Technological Challenges

Edited by Igor Girka

Chapter metrics overview

1,263 Chapter Downloads

View Full Metrics

Abstract

Cable-in-conduit conductor (CICC) has wide applications, and this structure is often served to undergo heat force-electromagnetic coupled field in practical utilization, especially in the magnetic confinement fusion (e.g., Tokamak). The mechanical behavior in CICC is of relevance to understanding the mechanical response and cannot be ignored for assessing the safety of these superconducting structures. In this chapter, several mechanical models were established to analyze the mechanical behavior of the CICC in Tokamak device, and the key mechanical problems such as the equivalent mechanical parameters of the superconducting cable, the untwisting behavior in the process of insertion, the buckling behavior of the superconducting wire under the action of the thermo-electromagnetic static load, and the Tcs (current sharing temperature) degradation under the thermo-electromagnetic cyclic loads are studied. Finally, we summarize the existing problems and the future research points on the basis of the previous research results, which will help the related researchers to figure out the mechanical behavior of CICC more easily.

Keywords

  • Nb3Sn
  • cable in conduit conductor (CICC)
  • cable stiffness
  • coefficient of thermal expansion
  • untwisting
  • current sharing temperature

1. Introduction

The ITER program is one of the largest and the most influential international energy technology cooperation projects, to verify the engineering feasibility of the magnetic confinement fusion. The core device of the magnetic confinement fusion reactor is the cable-in-conduit conductor (CICC). CICCs were used to build up the superconducting coil for generating strong magnetic fields to confine the high-temperature plasma in a confined space and maintaining the fusion reaction [1]. The ITER superconducting magnet systems mainly consist of four kinds of coils: 6 central solenoids (CS), 18 toroidal field (TF) coils, 6 sets of poloidal field (PF) coils, and 9 pairs of correction fields coils (CC) [2].

As early as the 1960s, the low-temperature superconducting material NbZr was processed into round wire and cables [3]. Subsequently, the superconducting magnets were wound with the structure of an internally cooled conductor (ICS) [4]. The superconducting strand is cooled to the superconducting state by the heat transfer copper tube with liquid helium in it [5]. But, the contact cooling method of ICS is inefficient, and the superconducting material is inclined to have a magnetic flux jump which will make the magnetic system to be quenched. In 1975, Hoenig et al. suggested subdividing the superconductor into strands to suppress the flux jump and twisting them into a cable to reduce the AC losses [6]. In 1980, Lue et al. proposed a cable-in-conduit design, and the innermost part was a perforated copper tube or a high-hardness stainless steel spring to form a liquid helium fast-flowing channel [7]. These two designs are the prototypes of modern CICC conductors. Nowadays, the CS and TF conductors with higher magnetic fields in the ITER project were fabricated by more than 1000 Nb3Sn wires. Disadvantageously, the superconducting properties of Nb3Sn are sensitive to mechanical deformation, which means that the tensile, compressive, and torsional deformations all lead to the reduction of the critical current [8]. Therefore, the strain state of the Nb3Sn strand cannot be ignored. Therefore, during its design, manufacture, and operation stage of the CICC, the mechanical analysis is needed.

Many studies have been published on the equivalent mechanical parameters of the twisted cable with two dimensions, such as the compression modulus of the cross section [9, 10], rather than built a complete three-dimensional model of the twisted cable. Feng et al. have applied the thin rod model to CICC conductor analysis and established the spatial geometry of each superconducting strand in the CICC conductor [11]. Qin et al. have applied the thin rod model to the mechanical analysis of superconducting cables and derived the axial stress–strain curves of primary cables and high-order strands [12]. The influence of pitch on the elastic modulus of the stranded cable and the curvature of the strand has been discussed. The introduction of the copper strand has been found to greatly reduce the axial stiffness of the strand, and the contact deformation between the strands has been found to have a great influence on the stiffness of the strand. The theoretical calculation results have appeared to be in good agreement with the experiments [13, 14]. Yue et al. have conducted a systematic mechanical analysis of the CICC in the design, preparation, and operation stage [15, 16, 17, 18].

In this chapter, first of all, we focus on the equivalent mechanical parameters of the superconducting cable; second, we concentrate on the untwisting behavior in the process of insertion; third, we want to explain the buckling behavior of the superconducting wire under the action of the thermo-electromagnetic static load; and finally, the Tcs degradation under the thermo-electromagnetic cyclic loads is studied. Our goal is to relate the cable stresses and buckling behavior to the thermal and electromagnetic loads so that relations between cable stress and current transport characteristics are built completely.

Advertisement

2. The equivalent mechanical parameters of the CICC

The mechanical behaviors of CICC have two main problems of structure and operation. On the one hand, the equivalent modulus of the cable is dependent on the manufacture parameters such as pitch, porosity, and radius. On the other hand, the electromagnetic load and the extremely low temperature make the internal stress and strain state of the cable difficult to analyze. Therefore, the thin rod model is applied to calculate the equivalent mechanical parameters of CICC conductors.

2.1 The tensile stiffness of the triplet

From the geometry characteristics of the cable, we know that the CICC superconducting cables have a complex structure with five stages of spirals. Therefore, each wire in the cable can be simplified into a thin rod which is elongated in the axial direction under the axial tensile load, and the wires can be contacted tightly or rotate in the lateral direction. The deformation and the force analysis of the wires are shown in Figure 1(a) and (b), respectively.

Figure 1.

The position of a spiral rod (a) and loads acting on the wire (b) [19].

The curvature and torsion of the spiral are as follows:

κ=0;κ=vdtcos2αrds=cos2αr;τ=vdtcosαsinαrds=cosαsinαr.E1

The equilibrium equation of forces and moments can be expressed as

dNds+TκNτ+X=0;dNdsTk++Y=0;dTdsNk+Nk+Z=0;dGdsGτ+HkN+K=0;dG'dsHk++N+K=0;dHdsGk+Gk+Θ=0.E2

Assuming that the thin rod is isotropic and elastic, the moments in any cross section with respect to the axis x,y,z can be written as Ix,Iy,Iz, and the constitutive equations are given by

G=EIxΔκ;G=EIyΔκ;H=GIzΔτ.E3

In the result of solving Eq. (3) with the account for the temperature terms, the expression of the equivalent coefficient of thermal expansion in the axial direction of the triplet can also be given by [15]

αeff=αLαTcos2αsin2α.E4

In Eq. (4), αeff is the equivalent coefficient of thermal expansion of the strand, αL is the coefficient of thermal expansion of the strand in the longitudinal direction, αT is the transverse coefficient of thermal expansion of the strand, and α denotes the helix angle.

2.2 The tensile stiffness of the higher stage strand

Based on the equivalent modulus and thermal expansion of the triplet, the space and the 2D view of the triplet and single wire are shown in Figure 2(a) and (b), respectively.

Figure 2.

Space line of the triplet and single wire (a) and 2D view of the triplet and single wire (b).

The conversion relationship between the local coordinates of the triplet and the higher-level strand can be expressed as

Tk=cosθksinθk0sinθksinαkcosθksinαkcosαksinθkcosαkcosθkcosαksinαk.E5

The curvature and torsion of the secondary stage strand can be given by

κp2κb2κt2=T2T1T1T00cosα2r2+00cosα1r1sinα2.E6

According to the geometric compatibility of the secondary-stage strand, the deformation of the triplet is equal to the tangential strain of the secondary-stage strand, and the torsion of the triplet is equal to the twist angle of the secondary cable. The axial loads and torque of the secondary-stage strand can be obtained. The equilibrium equations can be expressed as

00Ft0=33T1T2Fp2Fb2Ft2E7
Mp1Mb1Mt1=3T2Mp2Mb2Mt2+3r2Ft2sinα2+Fb2cosα2sinθ2r2Ft2sinα2+Fb2cosα2cosθ2r2Ft2cosα2Fb2sinα2.E8

Theoretically, the tensile stiffness of the conductor can be deduced by four times transformations as the CICC conductor has a five-stage twist structure.

Advertisement

3. Rotation analysis of the CICC

In the CICC conductor manufacture process, they twist a superconducting cable and penetrate it into the stainless steel tube. However, due to the friction between the superconducting cable and the stainless steel armor, the drag force of the cable is as high as several tons during the cable penetration. The friction force of the pipe leads to the axial elongation of the superconducting cable, accompanied by the untwisting of the cable, which causes the cable pitch to increase. This makes that the pitch is much larger than the ITER requirement [20]. Therefore, the untwist behavior of the cable must be controlled [17].

In this section, the untwist model is described. The large-scale cable is considered, e.g., ITER TF, CS, and CFETR CSMC. The components of the final cable include petals, central cooling spiral, and wrap, as shown in Figure 3. The model ignores the friction between the jacket and the cable, only modeling de-twists of the cable under the insertion force FInsert.

Figure 3.

The dimensions and parameters of a large-scale cable (e.g., CFETR CSMC).

The cable is divided into three parts in the model: central cooling spiral, six petals, and the wrap. The twist direction of wrap and cooling spiral is left and with the reverse direction for the petal. The torsion constraint is free for the cable when there is undering the uniaxial tension. Therefore, the boundary conditions can be set as F=F0;M=0. The force of the whole cable is from those acts on wrap, sub-cables, and central cooling spiral, which can be described as follows:

F=Fin+Fp+Fst=F0,M=Min+Mp+Mst=0.E9

In Eq. (9), F is the insertion force in the axial direction for the cable. Fst,Fp,Fin are the forces loading on the stainless steel wrap, petals, and inner cooling spiral, respectively; M represents the torque of the whole cable. Mst,Mp,Min denote the torques of stainless steel wrap, petals, and inner cooling spiral.

The force distribution and spatial relationship are shown in Figure 4. Gx,Gx,Hx are sectional moment components of the thin rod. Nx,Nx,Tx are the sectional force components. Xx is the contact force. αx is the spiral angle of each component. rx is the distance between the centroid and the center of the cable, and Lx is the twist pitch.

Figure 4.

Force distribution (A) and spatial relationship (B and C).

We simplified the central cooling spiral into a single helical thin plate. The change of spiral angle, radius, and the axial strain are defined as Δαst,Δrstrst,ξst, respectively. Then, the geometrical equations of the He-inlet spiral, six petals, and stainless steel wrap can be deduced by Costello-Velinsky theory [19]; we can obtain

Heinlet:ε0=ξin+Δαintanαin,β0=ξintanαin+1tanαinΔrinrinΔαin,Petal:ε0=ξP+ΔαPtanαP,β0=ξPtanαP+1tanαPΔrPrPΔαP,Wrapstiffness:ε0=ξst+Δαsttanαst,β0=ξsttanαst+1tanαstΔrstrstΔαst.E10

In the equations, the axial deformation and torsion angle of the cable are expressed as ε0,β0. The only axial tension is considered for central cooling spiral and wrap. Therefore, the equations can be updated as

Heinlet:Gin=0;Gin=0;Hin=0Tin=EinAinξinPetal:GP=0;GP=EPIPxΔκP;HP=GPIPzΔτPTP=EPAPξPWrapstiffness:Gst=0;Gst=0;Hst=0Tst=EstAstξstE11

In the equation systems, E refers to Young’s modulus of each component. A refers to the cross section of each component. EPIPx,GPIP are the bending and torsion stiffness of petals, respectively. ΔκP,ΔτP are the changes in curvature and twist, respectively. ξP is the axial strain of petal.

The balance equations and the compatible equations of central cooling spiral, subcables, and wrap can be obtained as

Nx=Gxτx+Hxκx,Xx=NxτxTxκx,Fx=Txsinαx+Nxcosαx,Mx=Hxsinαx+Gxcosαx+TxrxcosαxNxrxsinαx.E12
rin=Rin;rP=Rin+RP;rst=Rin+2RP,Δrin=ΔRin;ΔrP=ΔRin+ΔRP;Δrst=ΔRin+2ΔRP.E13

where x stands for anyone of in, p, st, which represent inner He-inlet spiral, petal, and stainless steel wrap, respectively; Rp denote the twist radius of petal, and ΔRP=νRPξP, v is Poisson’s ratio of the petal.

The contact force loading on the petal can be written as follows:

6XPsinαP=XstsinαstXinsinαinE14

where XP stand for the line pressure between the petal and inner He-inlet spiral; Xst is the uniform force between the stainless steel wrap and petal; and Xin is the reaction force between inner He-inlet spiral and petal.

From Eqs. (9)(12), one can get

F=F0=6TPsinαP+NPcosαP+Tstsinαst+TinsinαinM=0=6HPsinαP+GPcosαP+TPrPcosαPNPrPsinαPTstrstcosαstTinrincosαinE15

From the above Eqs. (10)(15), we can obtain

(1tanαst+tanαp)ξstυrptanαst2rin+2RPξP+0ξin+1tanαst1rin+2RPΔrintanαstε0β0=00ξst+(1tanαp+tanαpυRPtanαP1rin+RP)ξP+0ξin+1tanαP1rin+RPΔrintanαPε0β0=00ξst+0ξP+(tanαin+1tanαin)ξin+1rintanαinΔrintanαinε0β0=0sinαP6sinαst(EA)stκstξst+[ GBτP+HBκP(EA)PκP ]ξP+sinαP6sinαin(EA)inκinξin+(GDτP+HDκP)Δrin+(GEτP+HEκP)ε0+0β0=0(EA)stsinαstξst+[ 6GBτPcosαP+6HBκPcosαP+6(EA)PsinαP ]ξP+(EA)insinαinξin+6(GDτP+HDκP)cosαPΔrin+6(GEτp+HEκP)cosαPε0+0β0=F0(EA)stcosαstrstξst+ [ 6HBsinαP+6GBcosαP+6(EA)P(rin+rP)cosαP 6(GBτP+HBκP)rPsinαP ]ξP(EA)incosαinrinξin+[ 6HDsinαP+6GDcosαP6(GDτP+HDκP)rPsinαP ]Δrin+[ 6HEsinαP+6GEcosαP6(GEτP+HEκP)rPsinαP ]ε0+0β0=M0E16

The twist angle β0 can be computed from the equation system (16). So, the rotation of cable can be evaluated [17].

The experimental and numerical results are shown in Figure 5. First, it is easily found that the cabling tension has less impact on cable rotation. Taking the wrapping tension with 200 N, for example, there is no deviation between two different cables with cabling tension 200 N and 800 N, respectively. This result is in good agreement with the numerical model results. Second, the untwisting of the cable was mainly controlled by wrapping tension. Therefore, we can reduce the rotation significantly by increasing the wrapping tension. When insertion force is about 40 kN, the untwisting angle of cables with 600 N wrapping tension is about half of those cases with 200 N.

Figure 5.

Rotated angle per meter as a function of the force: Numerical and experimental results.

Advertisement

4. Buckling behavior in the CICC

4.1 Analytical model without the electromagnetic force

It is known that all the ITER CS and TF coils are wounded by CICCs, which made up of five-stage sub-cables formed around a central cooling tube. The petals and the sub-cables are wrapped with stainless steel tapes. Then, let the wrapped cable inserted into a stainless steel tube, which act as an amor. We assume that the total length of free segment of the superconducting strand on the surface of the cable is L (equal to the twist pitch), and set the fixed constraints on both sides, the wrap band as a uniform spring constraint. The schematic of this analytical model is illustrated in Figure 6. The lengths of spring constraint on both sides are equal to L1 and L3, respectively. The length of the free fraction of the strand is of L2; we can get L1+L2+L3=L.

Figure 6.

Schematic of the strand model ignores the EM force. L,L1,L2 and L3 denote the length of the twist pitch, left spring constraint, free segment, and right spring constraint, respectively. k denotes the rigidity of the bandaging. EI is the bending stiffness of the strand. εT is the total compression strain, which is equal to the sum of εThermal and εSlid.

Since the coefficient of thermal expansion of stainless steel between 923 and 4.2 K is approximately twice that of the Nb3Sn strand, then the superconducting cable is in compression at the end of the cooling. The thermal shrinkage of the cable is denoted by εThermal. Other than this, the strands of the cable can be squeezed into another side of the conduit by the large electromagnetic force; then, there will generate a large void on the other side of the conduit. Due to the gap, the friction force between the strands and the stainless steel armor decrease significantly. As there is no lateral restraint by wrap or friction, the surface strands around the void will show bending deformation by the thermal mismatch. In addition, the strand will slide into the high magnetic field region that will aggravate its bending behavior. εSlid is the stand for this slid strain. Therefore, the total compression strain of the strand εT is the sum of εSlid and εThermal. In this case, εT=εThermal+εSlid.

The mechanical analysis of the strand microelement is shown in Figure 7. The equilibrium equations for the moments are as follows:

dMdx+NdydxQv=0E17

Figure 7.

Mechanical analysis for the strand microelement.

In Eq. (17), Qv,M,N denote the vertical shear force, the bending moment, and the compression force along the axial direction, respectively. After submitting M=EId2y/dx2 into Eq. (17), and making a substitution dQv/dx=ky, we can get the differential equation of the rod with the spring constraints:

EId4ydx4Nd2ydx2+ky=0E18

If the strand has no wrapping, one can get

EId4ydx4Nd2ydx2=0E19

where k denotes the stiffness of the wrapping.

The continuity of displacements, twist angles, bending moments, and shear forces across the coverage/free strand require

wi=wi+1dwidxi=dwi+1dxi+1d2widxi2=d2wi+1dxi+12Ndwidxi+d3widxi+13=Ndwi+1dxi+1+d3wi+1dxi+13E20

In Eq. (20) the subscript i represents the number of wrapping/free segments. The fixed boundary conditions at the coverage fraction have the form:

wx=0=wx=L=0dwdxx=0=dwdx|x=L=0E21

Thus, the general solution of Eq. (18) is obtained as

wi=ai1cosλxi+ai2sinλxi+ai3xi+ai4E22

This general solution can be divided into three situations with a variation of the external force P [21, 22, 23]:

  1. If P2EIk, the general solution has the form as

    wi=ai1coshλ1xi+ai2coshλ2xi+ai3sinhλ1xi+ai4sinhλ2xiE23

  2. If P=2EIk, it becomes

    wi=ai1cosλxi+ai2cosλxi+ai3xisinλxi+ai4xisinλxiE24

  3. If P2EIk, it becomes

    wi=ai1cosλ1xi+ai2cosλ2xi+ai3sinλ1xi+ai4sinλ2xiE25

Here, λ,λ1,λ2 and λ3 are the eigenvalues, respectively; and aij is the constant coefficients.

The axial compression strain can be neglected when the strand gets into the buckling state, as the compression force N is small. Therefore, the total compression of the strand can be expressed as

120Lw2dx=LεTNLEAE26

Here, w denotes the curvature of the strand transverse deflection. The radius of curvature ρ has the form

ρ=d2wd2xE27

Here, w stands for the strand transverse deflection. The corresponding maximum strain is

εmax=D2ρE28

Here, D stands for the diameter of the strand. If the maximum strain is larger than 1%, the strand would be considered as cracking [10, 24, 25, 26]. Based on these equations, the coefficient aij can be calculated, as well as the buckling deflection, the relationship between the radius of curvature, and the thermal compression strain εThermal or slid strain εSlid.

4.2 Analytical model with the electromagnetic force

As we know, the magnetic field is maintained at 13 T in the Tcs test procedure, and the electromagnetic force rises with the increasing the current. Therefore, the electromagnetic force is a cyclic load. The strand on the surface of the cable where the gap is formed is selected in this model. The EM force FEM is perpendicular to the Nb3Sn strand, and its direction points to the inner part of the cable. Thus, the strand at the lower loading side has the least lateral constraint that means the strand most likely to have a buckling behavior in the lateral direction. Therefore, the only thing that can prevent the buckling of the strands is the friction force against the cable. The friction factor can be written as a symbol μ.

In Figure 8, μL and μT are the axial friction factor and lateral friction factor, respectively; Ls and L denote the slipping length and buckling length, respectively; L+2Ls is the twist pitch of the first stage; FEM stands for the EM force, and then in the buckling area μTFEM is the lateral constraint on the strand; and P0 and P are the compression force of the strand in the slipping and buckling area, respectively.

Figure 8.

Schematic of the analytical model with the EM force.

Figure 9.

Schematic view of the mechanical analysis for the microelement.

Accordingly, by the torque balance of the microelement, as shown in Figure 9. The equilibrium equation is as follows:

EIy=Py+μTFEM2L2xL2xμTFEML2L2xE29

In Eq. (29) EI stands for the bending stiffness of the Nb3Sn strand; y stands for the deflection of the buckling. After simplifying formula (29), one can get

EIy+Py+μTFEM2L24x2=0E30

The general solution [27] of formula (30) is given as

y=mn4Acosnx+Bsinnxn2x22+n2L28+1E31

In Eq. (31), m=μTFEMEI,n2=PEI. Two boundary conditions are needed to solve the undetermined coefficients. After applying yx=±L/2=0,yx=0=0 to formula (31), one can get

y=mn4cosnxcosnL/2n2x22+n2L28+1E32

If we substitute the boundary yx=±L/2=0 to Eq. (32), Eq. (32) becomes

y=tannL/2nL/2=0E33

We can deduce that nL=8.9868 from Eq. (33). Therefore, the axial compression force in the buckling area is P=80.76EI/L2. There must be a balance in the axial direction between the buckling and slipping segments; we can get that

P0P=μLFEML2+μLFEMLSE34

The geometric compatibility equation of the strand can be expressed as

P0PLEAL/2L/212y2dx+μLFEMLS2EA=0E35

Submitting Eq. (32) into Eq. (35), one can get the relation between axial force and the length of the buckling area:

P0PLEA1.597×105μTFEMEI2L7+μLFEMLS2EA=0E36

Submitting Eq. (34) into Eq. (36) and eliminating the Ls, one can get P0=P+1.597×105EAμLFEMμTFEMEI2L7L24μLFEM2, in which P=80.76EI/L2.

Consider the special case LS=0 that means the end of the strand is locked and the length of the buckling is equal to the twist pitch. In this case Eq. (36) can be simplified as follows:

P0=80.76EIL2+1.597×105EAμTFEMEI2L6E37

Based on these equations, the critical buckling load P0 can be calculated, as well as the relationship between the buckling length, bending stiffness, and the friction factor.

We can know that the strand buckling behavior is depending on the twist pitch of the first stage; the shorter the lay length, the lower the possibility of the strand buckling. The higher the wrap rigidity, the stronger the strand. Furthermore, if we fixed the coverage rate, and with a narrow wrap, which would lead to almost uniform lateral supporting, it can also prevent the strand slid into buckling. This mechanism of buckling is shown in Figure 10.

Figure 10.

Mechanism of the strand buckling behavior during the cooling and conducting processes.

When there is no thermal load and EM force, the original strand in the CICCs is shown in Figure 10(a). When the temperature is cooled down, the axial strain caused by thermal mismatch and the lateral compression raised by the EM load are applied to the strands. We can see that the initial “straight” strand was keeping its shape by well lateral support. When the strand working in a large magnetic field (including self-field), it bears a huge EM force. A large void is generated, at the same time, some strands bending, which is shown in Figure 10(b). Additionally, this bending strain is not the only factor to make the strand fracture. For the SULTAN measurements [28], the samples is about 3.6 m long, and the high-field region is about 0.4 m. During the EM cyclic loading, the cable can slide into the high-field region, which can accelerate the wire bending and leading the strand to fracture. This process is shown in Figure 10(c).

Advertisement

5. Mechanical behavior of the CICC caused by electromagnetic force and thermal mismatch

The CICC qualification test samples show gradual degradation of the current-sharing temperature (Tcs) under several hundreds of EM cycles [29, 30], which leads to the Nb3Sn strand’s bending or compressing deformation. In this section, we focus on the relationship between Tcs and axial strain of the cable.

It is known that the ITER CS and TF coils are wounded by CICCs, which made up of five-stage sub-cables formed around a central cooling tube. The fourth petals and the sub-cables are wrapped with stainless steel tapes. Then, the wrapped cable were inserted into a stainless steel tube, which act as an amor, as illustrated in the cross-sectional view in Figure 11. The inner diameter of the CS armor is equal to 36 and 37 mm for TF, respectively. The side length of CS conductor is 49 mm, and the diameter of the TF conductor is of 40 mm. More than 1000 Nb3Sn strands were in the CICC conductors. In this section, the conductors can be simplified into a rope and the petals with circle cross section, as shown in Figure 11. Thus, the cable could be analyzed by using the thin rod model as shown in Figure 12.

Figure 11.

Schematic illustration of the CS and TF CICC cross sections. The symbols Rin and R denote the radius of the He channel and petal, respectively; r is on behalf of the sum of Rin and R: r=RIn+R.

Figure 12.

(a) The deformation of the cable caused by the axial compression; (b) the cross section of the cable.

During the calculations, the parameters R, Rin, and twist pitch (h) are 6, 6, and 450 mm, respectively. According to the geometric relation h=2πrtanα, as shown in Figure 13(b), we can get α=80.5, where α represents the initial helix angle. Figure 13(a) also displays the loads acting on the petal and the geometric relation of the petal centerline.

Figure 13.

(a) Uniform forces acting on the petal; (b) geometric relation of the petal centerline.

Assume that the petals were contact with each other in the original state. The curvature and the twist per unit length of the petal is k,k,τ; then the changes can be written as Δk,Δk,Δτ. They can be expressed as [19]

Δκ=cos2α¯r¯cos2αr=cos2αrΔrr2sinαcosαrΔαE38
Δτ=cosα¯sinα¯r¯cosαsinαr=sinαcosαrΔrr+12sin2αrΔαE39

In Eqs. (38) and (39), α and α¯ is the initial and final helical angle of the petal, Δα and Δr stand for the change of helical angle and radius of the petal, and r and r¯ denote the original and final radius, respectively. The loads and moments can be deduced as.

G=EπR44Δκ;H=EπR441+νΔτ;T=πER2ξE40

Here, T,G,H stand for the axial load, bending moment, and twist moment of the petal, respectively; E is of the axial stiffness of the petal, v is Poisson’s ratio, and ξ stands for the strain in the axial direction of the petal. If the petal is free in the initial state and there is no contact force under the compression state, then the uniform load between the petals is equal to zero. According to the thin rod model presented by Costello [19], the following equation is satisfied:

X=NsinαcosαrTcos2αr=0E41
N=Hcos2αrGsinαcosαrE42

In Eqs. (41) and (42), X stands for the resultant contact force per unit length of a petal; N is of the shear force acting on the petal. From Eqs. (38)(40) and Eq. (42), the shearing force N acting on the petal can be written as

NER2=πvsinαcos3α4r/R21+υΔrr+πcos2α4r/R212sin2α1+υ+2sin2αΔαE43

Substituting Eqs. (40) and (43) into Eq. (41), the axial strain of the petal ξ which is only related to the Δr and Δα can be written as

ξ=υsin2αcos2α4r/R21+υΔrr+sinαcosα4r/R212sin2α1+υ+2sin2αΔαE44

The deformed configuration of the petal in Figure 13 yields

ε=h¯hh=ξ+ΔαtanαE45

From Eqs. (44) and (45), one can get

ε=υsin2αcos2α4r/R21+υΔrr+cosαsinα+sinαcosα4r/R212sin2α1+υ+2sin2αΔαE46

The angle of twist per unit length τ of the petal can be defined by the expression

rτc=rθ¯θh= r r¯1+ξtanαΔα1tanαE47

That is,

rτc=υsinαcos2α4r/R21+υcosαsinαΔrr+cos2α4r/R212sin2α1+υ+2sin2α1ΔαE48

At the two ends of the petal per twist pitch, the rotation is zero, and the compression strain obeys: ε=ε0; then, Eq. (46) and Eq. (48) become

ε0=υsin2αcos2α4r/R21+υΔrr+cosαsinα+sinαcosα4r/R212sin2α1+υ+2sin2αΔα0=υsinαcos2α4r/R21+υcosαsinαΔrr+cos2α4r/R212sin2α1+υ+2sin2α1ΔαE49

The relation between axial strain ε0 and transverse strain εTrans=Δrr can be expressed as

εTrans=C4C2C3C1C4ε0E50

where C1=υsin2αcos2α4r/R21+υ, C2=cosαsinα+sinαcosα4r/R212sin2α1+υ+2sin2α,

C3=υsinαcos2α4r/R21+υcosαsinα, C4=cos2α4r/R212sin2α1+υ+2sin2α1.

It is found that the coefficient between transverse and axial strains is affected by the helical angle α and Poisson’s ratio υ. When the twist pitch of the fifth stage cable is of 427, 450, and 476 mm, the corresponding helical angle of the petal equals to 80, 80.5, and 81°, respectively. Substituting these values into Eq. (50), one can see that the axial strain of the cable ε0 has a linear relationship with the transverse strain εTrans, which is displayed in Figure 14. We carried out an experiment by compressing a CS cable which was fabricated in CASIPP to validate the analytical model. When the CS cable specimen is compressed along axial direction, it extends transversely with high resolution. The transverse extension can be measured by using a laser sensor. The experimental setup is schematically illustrated in the inset of Figure 14. One can find that the theoretical model shows perfect consistency with the experiment.

Figure 14.

Relationship between the axial strain ε0 and the transverse strain εT. The colored lines are obtained by using the presented theoretical model. The green triangle symbols denote the experimental results.

In Figure 14, we can get that the CS conductors with baseline and LTP (long twist pitch) shows a Tcs degradation after EM cycles, while for the samples with STP (short twist pitch) shows an opposite phenomenon, the Tcs becomes a constant even have an enhanced. These experimental results presented in Ref. [30] have a contrast with analytical prediction for many years. Using the theoretical model proposed in this chapter, the Tcs enhancement and degradation behaviors can be predicted quantitatively. Figure 15 shows a comparison of calculated results and experimental results, where the lines show the theoretical results and the dots are the experimental results. It is found that during the primary stage of the EM cycle there was a quick increase of axial strain that will lead to the Tcs drop dramatically. With the same reason, samples with LTP and baseline show an obviously degradation; this same rule applies to TF conductors. However, several STP CS conductors have an increase of Tcs, which means that the initial axial compression strain has been released by the EM loads, squeezing the cables at the high-field zone into the low-field region.

Figure 15.

A comparison of the experimental results and theoretical results based on the presented model [18].

While it is easy to imagine when the expansions of the cable in the high-field zone get stacked, the strain will not be released and the Tcs will drop (see the black dash line in Figure 15). If the inner part of the conductor is smooth enough, the axial compression strain is released completely, the Tcs reaches its upper boundary which is shown by the orange dash line with circle symbol corresponding to the length of the high-field region is 400 mm. If the high-field region’s length is 1000 mm, then the lower boundary of the enhancing Tcs is given with purple line with triangle symbol. We can find that the theoretical results agree with the experimental results very well.

Advertisement

6. Discussion

In the past half-century, the structure of CICC conductor was under continuous optimization and improvement. The CICC conductors have so many advantages such as good self-support, high operational stability, high current carrying capacity, low AC loss, etc., and they are widely used in the superconducting magnets. The nuclear fusion device CFTER built by China has also chosen the CICC conductors. However, there still have some problems to be solved: (1) there is a necessity for theoretical model to explain the degradation of critical current caused by mechanical deformation for the Nb3Sn strand; (2) accurate description of the mechanical response of complex strand structures in the coupling fields remains a challenging problem; (3) the untwisting behavior during the cable penetration is still not suppressed completely. One needs a more effective model to optimize the manufacturing process; and (4) the long-term stability and real-time monitoring of superconducting magnets are also a challenge for the engineers. Based on this, the equivalent mechanical parameters of CICC conductors and their mechanical behavior under coupled fields will be further studied. The authors hope that these models can provide a valuable reference for the related researchers.

Advertisement

Acknowledgments

This work is supported by the Funds of the National Natural Science Foundation of China (Nos. 11622217, 11802291, and 11872196), the National Key Project of Scientific Instrument and Equipment Development (11327802), and National Program for Special Support of Top-Notch Young Professionals. This work is also supported by the Fundamental Research Funds for the Central Universities (lzujbky-2017-ot18 and lzujbky-2017-k18).

References

  1. 1. Motojima O. The ITER project construction status. Nuclear Fusion. 2015;55(10):104023
  2. 2. Mitchell N, Bessette D, Gallix R, et al. The ITER magnet system. IEEE Transactions on Applied Superconductivity. 2008;18(2):435-440
  3. 3. File J, Martin GD, Mills RG, et al. Stabilized, levitated superconducting rings. Journal of Applied Physics. 1968;39(6):2623-2626
  4. 4. Hoenig MO. Internally cooled cabled superconductors. Cryogenics. 1980;20(7):373-389
  5. 5. Lubell MS, Chandrasekhar BS, Mallick GT. Degradation and flux jumping in solenoids of heat-treated Nb-25% Zr wire. Applied Physics Letters. 1963;3(5):79-80
  6. 6. Hoenig M, Montgomery D. Dense supercritical-helium cooled superconductors for large high field stabilized magnets. IEEE Transactions on Magnetics. 1975;11(2):569-572
  7. 7. Lue J. Review of stability experiments on cable-in-conduit conductors. Cryogenics. 1994;34(10):779-786
  8. 8. Ekin J. Strain scaling law for flux pinning in practical superconductors. Part 1: Basic relationship and application to Nb3Sn conductors. Cryogenics. 1980;20(11):611-624
  9. 9. Nijhuis A. A solution for transverse load degradation in ITER Nb3Sn CICCs: Verification of cabling effect on Lorentz force response. Superconductor Science and Technology. 2008;21(5):054011
  10. 10. Zhu JY, Luo W, Zhou YH, et al. Contact mechanical characteristics of Nb3Sn strands under transverse electromagnetic loads in the CICC cross-section. Superconductor Science and Technology. 2012;25(12):125011
  11. 11. Feng J. A cable twisting model and its application in CSIC multi-stage cabling structure. Fusion Engineering and Design. 2009;84(12):2084-2092
  12. 12. Qin J, Wu Y, Warnet LL, et al. A novel numerical mechanical model for the stress–strain distribution in superconducting cable-in-conduit conductors. Superconductor Science and Technology. 2011;24(6):065012
  13. 13. Qin J-G, Xue T-J, Liu B, et al. Cabling Technology of Nb3Sn conductor for CFETR central solenoid model coil. IEEE Transactions on Applied Superconductivity. 2016;26(3):1-5
  14. 14. Qin J, Dai C, Liu B, et al. Optimization of CFETR CSMC cabling based on numerical modeling and experiments. Superconductor Science & Technology. 2015;28(12):125008
  15. 15. Yue D, Zhang X, Yong H, Zhou J, Zhou YH. Controllable rectification of the axial expansion in the thermally driven artificial muscle. Applied Physics Letters. 2015;107(11):111903
  16. 16. Donghua Y, Xingyi Z, Youhe Z. Buckling behavior of Nb3Sn strand caused by electromagnetic force and thermal mismatch in ITER cable-in-conduit conductor. IEEE Transactions on Applied Superconductivity. 2017;27(7):1
  17. 17. Qin J, Yee D, Zhang X, Wu Y, Liu X, Liu H, et al. Rotation analysis on large complex superconducting cables based on numerical modeling and experiments. Superconductor Science and Technology. 2018;31(2):025001
  18. 18. Yue D, Zhang X, Zhou Y-H. Theoretical analysis for the mechanical behavior caused by an electromagnetic cycle in ITER Nb3Sn cable-in-conduit conductors. Acta Mechanica Sinica. 2018;34(4):614-622
  19. 19. Costello GA. Theory of Wire Rope. New York: Springer Science & Business Media; 1997
  20. 20. Nabara Y, Suwa T, Takahashi Y, et al. Behavior of cable assembled with conduit for ITER central solenoid. IEEE Transactions on Applied Superconductivity. 2015;25(3):1-5
  21. 21. Hetenyi M. Beams on Elastic Foundation: Theory with Applications in the Fields of Civil and Mechanical Engineering. 1971
  22. 22. Niu K. Compressive Behavior of Sandwich Panels and Laminates with Damage, Georgia Institute of Technology, Atlanta. 1998
  23. 23. Niu K, Talreja R. Buckling of a thin face layer on Winkler foundation with debonds. Journal of Sandwich Structures and Materials. 1999;1(4):259-278
  24. 24. Mitchell N, Devred A, Larbalestier DC, et al. Reversible and irreversible mechanical effects in real cable-in-conduit conductors. Superconductor Science and Technology. 2013;26(11):114004
  25. 25. Li Y, Ta W, Gao Y. Prediction for IC degradation of superconducting strand under bending. Cryogenics. 2013;58:20-25
  26. 26. Li Y, Yang T, Zhou Y, et al. Spring model for mechanical–electrical properties of CICC in cryogenic–electromagnetic environments. Cryogenics. 2014;62:14-30
  27. 27. Hobbs RE. Pipeline buckling caused by axial loads. Journal of Constructional Steel Research. 1981;1(2):2-10
  28. 28. March SA, Bruzzone P, Stepanov B, et al. Results of the TFEU6 sample tested in SULTAN. IEEE Transactions on Applied Superconductivity. 2013;23(3):4200204
  29. 29. Martovetsky NN, Reiersen WT, Murdoch GR, Bruzzone P, Stepanov B. Qualification of the US conductors for ITER TF magnet system. ITPS. 2018;46(5):1477-1483
  30. 30. Devred A, Backbier I, Bessette D, et al. Challenges and status of ITER conductor production. Superconductor Science and Technology. 2014;27(4):044001

Written By

Donghua Yue, Xingyi Zhang and Youhe Zhou

Submitted: 11 September 2018 Reviewed: 31 October 2018 Published: 31 December 2018