This paper presents a new code for the two-dimensional discrete element method (DEM) and relevant simulations to quantitatively characterize the contact force behavior of the Nb 3 Sn strands in the ITER CICC cross-section under a transverse electromagnetic load. In order to obtain the essential parameters in the contact force model employed in the DEM, a simulation of the experiments conducted by Nijhuis et al (2004 IEEE Trans. Appl. Supercond. 14 1489 is first performed, where the load-displacement curve predicted by the code is in good agreement with the measurements. After that, the contact force chain between strands and its distribution is quantitatively analyzed by the code. It is found that the contact force distribution among strands is heterogeneous and strongly anisotropic. In other words, the force chain distribution, which determines the behavior of the assembly of strands with discrete media, and the distribution of area average magnitude of the contact force are obviously inhomogeneous. To describe this inhomogeneity, here, the probability density function (PDF) is used in the statistical analysis. The numerical results show that the PDFs of the magnitudes of the resultant contact force, normal contact force, and tangential contact force all decay with an exponential law, and that PDFs of the directions of the contact forces are all anisotropic and exhibit about six periodic changes in which the peak values in the direction parallel to the applied electromagnetic load are appreciably larger than the other peaks.

No abstract