{"id":126,"date":"2022-08-23T11:41:57","date_gmt":"2022-08-23T02:41:57","guid":{"rendered":"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en2\/?post_type=case&#038;p=126"},"modified":"2022-08-23T14:02:15","modified_gmt":"2022-08-23T05:02:15","slug":"coil_force","status":"publish","type":"case","link":"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en\/case\/coil_force\/","title":{"rendered":"COIL inductance and electromagnetic force calculations"},"content":{"rendered":"<h3>Summary<\/h3>\n<p>This section introduces how to calculate the inductance and electromagnetic force of an external current magnetic field source (COIL).  <\/p>\n<h3>Explanation<\/h3>\n<p>COIL is a mesh-independent feature that is easy to use. An example of an axisymmetric model is presented in <a href=\"\/product\/EMSolution\/en\/case\/coil_inductance\/\" target=\"_blank\" rel=\"noopener noreferrer\" style=\"display:inline\"><font color=\"Red\">&quot;Inductance calculation for COIL (external current magnetic field source)&quot;<\/font><\/a>, but here we will use a model that is 3-dimensional and includes magnetic materials as an example. The accuracy is verified by comparing the two-potential method (A-Ar method), in which the iron core portion is the total potential region and the coil and air portions are the reduced potential regions, with the case in which the entire region is solved with the total potential and the coil is defined with the surface-defined current source SDEFCOIL (A method).  <\/p>\n<p>Fig. 1 shows the excitation coil. The coil consists of four rectangular cross sectional straight sections (GCE) and arc sections (ARC). When determining the inductance and electromagnetic force, the integral regions in Fig. 1 are defined by the integral elements GCE- and ARC-. The analysis is performed in the 1\/8 domain for symmetry. The integral element is defined as the part of COIL within the analysis domain. The magnetomotive force is set to 3000 AT and a linear static field analysis is performed.  <\/p>\n<div class=\"img col1\">\n<div>\n        <a href=\"\/product\/EMSolution\/en\/wp-content\/uploads\/coil_force01.png\" class=\"modal\"><br \/>\n        <img decoding=\"async\" src=\"\/product\/EMSolution\/en\/wp-content\/uploads\/coil_force01.png\" alt=\"\" \/><\/a><\/p>\n<p class=\"text01\">Fig.1\u3000Excitation coil and integral domain<\/p>\n<\/p><\/div>\n<\/div>\n<p>Fig. 2 shows the magnetic flux density intensity distribution obtained with the iron core (specific permeability 1000) and by using the A method. The result is almost the same as obtained using COIL. Table 1 shows the inductance calculation results using the A method. We analyzed the same problem with the boundary condition that Bn=0 and Ht=0 at far boundary. The results are almost the same, and the boundary condition is sufficient. Coil inductance is shown as a one turn equivalent value. The inductance can be calculated either by integrating A\u30fbJ in the coil region and assuming it equals $LI$, or by taking the total magnetic energy (the total spatial integral of $B^2\/2m_0$) as $LI^2\/2$. Both methods are in good agreement. In EMSolution, the flux amount of $LI$ is output as the flux of the coil series. We also provide an option to calculate the magnetic energy (MAGNETIC_ENERGY option).  <\/p>\n<p>On the other hand, Table 2 shows the results obtained by the A-Ar method. The inductance is calculated and added internally, so the coil inductance includes the iron core effect. In this case, too, the analysis was performed with two types of boundary conditions, and the results are consistent. The results are also consistent with the results obtained by the A method. Without regularization, the ICCG method does not converge to $10^{-8}$, but the results are almost identical.  <\/p>\n<div class=\"img col1\">\n<div>\n        <a href=\"\/product\/EMSolution\/en\/wp-content\/uploads\/coil_force02.png\" class=\"modal\"><br \/>\n        <img decoding=\"async\" src=\"\/product\/EMSolution\/en\/wp-content\/uploads\/coil_force02.png\" alt=\"\" \/><\/a><\/p>\n<p class=\"text01\">Fig.2\u3000Iron core and magnetic flux density intensity distribution (the A method)<\/p>\n<\/p><\/div>\n<\/div>\n<p>Table 1.\u3000Inductance calculation results using the A method<br \/>\n\n<table id=\"tablepress-5\" class=\"tablepress tablepress-id-5\">\n<thead>\n<tr class=\"row-1 odd\">\n\t<th class=\"column-1\">boundary condition<\/th><th class=\"column-2\">unknown number<\/th><th class=\"column-3\">Number of non-zero elements<\/th><th class=\"column-4\">ICCG<\/th><th class=\"column-5\">CPU Time<\/th><th class=\"column-6\">L (H) A\u30fbJ<br>Integral<\/th><th class=\"column-7\">L (H) B\u00d7B<br>Integral<\/th>\n<\/tr>\n<\/thead>\n<tbody class=\"row-hover\">\n<tr class=\"row-2 even\">\n\t<td class=\"column-1\">Bn=0<\/td><td class=\"column-2\">883,920<\/td><td class=\"column-3\">14,699,871<\/td><td class=\"column-4\">599<\/td><td class=\"column-5\">427.6<\/td><td class=\"column-6\">4.29417E-07<\/td><td class=\"column-7\">4.28800E-07<\/td>\n<\/tr>\n<tr class=\"row-3 odd\">\n\t<td class=\"column-1\">Ht=0<\/td><td class=\"column-2\">900,260<\/td><td class=\"column-3\">14,973,253<\/td><td class=\"column-4\">431<\/td><td class=\"column-5\">281.0<\/td><td class=\"column-6\">4.30673E-07<\/td><td class=\"column-7\">4.30062E-07<\/td>\n<\/tr>\n<tr class=\"row-4 even\">\n\t<td colspan=\"4\" class=\"column-1\"><\/td><td class=\"column-5\">average<\/td><td class=\"column-6\">4.30045E-07<\/td><td class=\"column-7\">4.29431E-07<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<!-- #tablepress-5 from cache --><\/p>\n<p>Table 2.\u3000Inductance calculation results by A-Ar method<br \/>\n\n<table id=\"tablepress-6\" class=\"tablepress tablepress-id-6\">\n<thead>\n<tr class=\"row-1 odd\">\n\t<th class=\"column-1\">boundary condition<\/th><th class=\"column-2\">normalization<\/th><th class=\"column-3\">unknown number<\/th><th class=\"column-4\">Number of non-zero elements<\/th><th class=\"column-5\">ICCG<\/th><th class=\"column-6\">CPU Time<\/th><th class=\"column-7\">L (H) A\u30fbJ<br>Integral<\/th><th class=\"column-8\"><\/th>\n<\/tr>\n<\/thead>\n<tbody class=\"row-hover\">\n<tr class=\"row-2 even\">\n\t<td class=\"column-1\">Brn=0<\/td><td class=\"column-2\">\u00d7<\/td><td class=\"column-3\">883,920<\/td><td class=\"column-4\">14,699,871<\/td><td class=\"column-5\">545<\/td><td class=\"column-6\">333.1<\/td><td class=\"column-7\">4.30087E-07<\/td><td class=\"column-8\">Divergence at 4.2E-07<\/td>\n<\/tr>\n<tr class=\"row-3 odd\">\n\t<td class=\"column-1\">Brn=0<\/td><td class=\"column-2\">\u3007<\/td><td class=\"column-3\">883,920<\/td><td class=\"column-4\">14,699,871<\/td><td class=\"column-5\">534<\/td><td class=\"column-6\">342.3<\/td><td class=\"column-7\">4.30087E-07<\/td><td class=\"column-8\"><\/td>\n<\/tr>\n<tr class=\"row-4 even\">\n\t<td class=\"column-1\">Htrt=0<\/td><td class=\"column-2\">\u3007<\/td><td class=\"column-3\">900,260<\/td><td class=\"column-4\">14,973,253<\/td><td class=\"column-5\">388<\/td><td class=\"column-6\">295.3<\/td><td class=\"column-7\">4.30590E-07<\/td><td class=\"column-8\"><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<!-- #tablepress-6 from cache --><\/p>\n<p>Table 3 and 4 show the results of electromagnetic force calculations for both methods, showing the total force in each direction for each section \u2460, \u2461, and \u2462 in Fig. 1. It can be seen that the results of both methods are almost identical.  <\/p>\n<p>Table 3.\u3000Electromagnetic force calculated by the A method<br \/>\n\n<table id=\"tablepress-7\" class=\"tablepress tablepress-id-7\">\n<thead>\n<tr class=\"row-1 odd\">\n\t<th class=\"column-1\"><\/th><th class=\"column-2\">Part<\/th><th class=\"column-3\">Fx ( N )<\/th><th class=\"column-4\">Fy ( N )<\/th><th class=\"column-5\">Fz ( N )<\/th>\n<\/tr>\n<\/thead>\n<tbody class=\"row-hover\">\n<tr class=\"row-2 even\">\n\t<td rowspan=\"3\" class=\"column-1\">Bn=0<\/td><td class=\"column-2\">\u2460<\/td><td class=\"column-3\">-1.28170E-01<\/td><td class=\"column-4\">2.17279E-07<\/td><td class=\"column-5\">-7.67471E-01<\/td>\n<\/tr>\n<tr class=\"row-3 odd\">\n\t<td class=\"column-2\">\u2461<\/td><td class=\"column-3\">-7.38606E-03<\/td><td class=\"column-4\">-7.38606E-03<\/td><td class=\"column-5\">-8.16667E-01<\/td>\n<\/tr>\n<tr class=\"row-4 even\">\n\t<td class=\"column-2\">\u2462<\/td><td class=\"column-3\">2.17479E-07<\/td><td class=\"column-4\">-1.28170E-01<\/td><td class=\"column-5\">-7.67471E-01<\/td>\n<\/tr>\n<tr class=\"row-5 odd\">\n\t<td rowspan=\"3\" class=\"column-1\">Ht=0<\/td><td class=\"column-2\">\u2460<\/td><td class=\"column-3\">-1.26507E-01<\/td><td class=\"column-4\">4.01467E-07<\/td><td class=\"column-5\">-7.68262E-01<\/td>\n<\/tr>\n<tr class=\"row-6 even\">\n\t<td class=\"column-2\">\u2461<\/td><td class=\"column-3\">-5.87456E-03<\/td><td class=\"column-4\">-5.87456E-03<\/td><td class=\"column-5\">-8.17471E-01<\/td>\n<\/tr>\n<tr class=\"row-7 odd\">\n\t<td class=\"column-2\">\u2462<\/td><td class=\"column-3\">4.01467E-07<\/td><td class=\"column-4\">-1.26507E-01<\/td><td class=\"column-5\">-7.68262E-01<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<!-- #tablepress-7 from cache --><\/p>\n<p>Table 4.\u3000Electromagnetic force calculated by the A-Ar method<br \/>\n\n<table id=\"tablepress-8\" class=\"tablepress tablepress-id-8\">\n<thead>\n<tr class=\"row-1 odd\">\n\t<th class=\"column-1\"><\/th><th class=\"column-2\">Part<\/th><th class=\"column-3\">Fx ( N )<\/th><th class=\"column-4\">Fy ( N )<\/th><th class=\"column-5\">Fz ( N )<\/th>\n<\/tr>\n<\/thead>\n<tbody class=\"row-hover\">\n<tr class=\"row-2 even\">\n\t<td rowspan=\"3\" class=\"column-1\">Brn=0<\/td><td class=\"column-2\">\u2460<\/td><td class=\"column-3\">-1.26026E-01<\/td><td class=\"column-4\">0.00000E+00<\/td><td class=\"column-5\">-7.67754E-01<\/td>\n<\/tr>\n<tr class=\"row-3 odd\">\n\t<td class=\"column-2\">\u2461<\/td><td class=\"column-3\">-6.03243E-03<\/td><td class=\"column-4\">-6.03243E-03<\/td><td class=\"column-5\">-8.17493E-01<\/td>\n<\/tr>\n<tr class=\"row-4 even\">\n\t<td class=\"column-2\">\u2462<\/td><td class=\"column-3\">0.00000E+00<\/td><td class=\"column-4\">-1.26026E-01<\/td><td class=\"column-5\">-7.67754E-01<\/td>\n<\/tr>\n<tr class=\"row-5 odd\">\n\t<td rowspan=\"3\" class=\"column-1\">Hrt=0<\/td><td class=\"column-2\">\u2460<\/td><td class=\"column-3\">-1.27201E-01<\/td><td class=\"column-4\">0.00000E+00<\/td><td class=\"column-5\">-7.68253E-01<\/td>\n<\/tr>\n<tr class=\"row-6 even\">\n\t<td class=\"column-2\">\u2461<\/td><td class=\"column-3\">-6.74528E-03<\/td><td class=\"column-4\">-6.74528E-03<\/td><td class=\"column-5\">-8.18000E-01<\/td>\n<\/tr>\n<tr class=\"row-7 odd\">\n\t<td class=\"column-2\">\u2462<\/td><td class=\"column-3\">0.00000E+00<\/td><td class=\"column-4\">-1.27201E-01<\/td><td class=\"column-5\">-7.68253E-01<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<!-- #tablepress-8 from cache --><\/p>\n<p>As shown above, the results of the A method with SDEFCOIL and the A-Ar method with COIL are almost in good agreement and both analyses are equivalent. In the above example, the mesh is sufficiently fine and the far boundary is sufficiently far away. However, please note that differences may appear if such conditions are not used.  <\/p>\n<p><!--more--><\/p>\n<h3>How to use<\/h3>\n<p>To calculate the inductance of COIL internally in EMSolution, set the CALC_IND option to 1. The inductance calculation is performed by the MAKE_SYSTEM_MATRICES process.  <\/p>\n<p class=\"slideText\">* NODE_ORDER * EDGE_ORDER * METRIC_MAD * QUAD_TRI <font color=\"Red\">* CALC_IND *<\/font><br \/>\n               1                         1                          0                      0                    <font color=\"Red\">1<\/font>\n<\/p>\n<p>To calculate the electromagnetic force of COIL, set COIL_FORCE in &quot;10. Input\/Output File&quot; to 1. Please note that previously &quot;11.1. Output Options&quot; FORCE_J_B=2 was used, but this has been changed. The INTEG_OPT option was added as a new input data. With INTEG_OPT=1, the element integration point in a finite element is set to the element center point; with INTEG_OPT=0, the element Gaussian integration point is used. When INTEG_OPT =1 is used, computation time is significantly reduced. This option is also applicable for the calculation of conventional evaluation point magnetic fields. For calculations at a sufficient distance from magnetic materials or conductors, INTEG_OPT =1 is sufficient.  <\/p>\n<p class=\"slideText\">* MESH * CURRENT * MAGNETIC * FORCE_J_B * FORCE_NODAL * DISP * ELEM * HEAT * MAGNETIZATION * IRON_LOSS *  E  <font color=\"Red\">* COIL_FORCE *<\/font><br \/>\n       1               0                    0                      0                         0                   0          0           0                      0                          0              0               <font color=\"Red\">1<\/font><br \/>\n* MESH *  A  *  V  *  B  * B_INTEG * J_Q *<br \/>\n       0         0      0      0            0            0<br \/>\n* Q_AVERAGE * HEAT * MAG_FLUX * CUR_FLUX *<br \/>\n              0               0                0                    0<br \/>\n* FORCE_J_B * FORCE_NODAL * MAGNETIC_ENERGY * IRON_LOSS *<br \/>\n             0                          0<br \/>\n* XSIMMETRY * YSIMMETRY * ZSIMMETRY * ROTATION *<br \/>\n             1                      1                    -1                     1<br \/>\n* NO_POINTS * NO_MAT_IDS * NO_COIL_SERIES <font color=\"Red\">* INTEG_OPT *<\/font><br \/>\n             0                        1                            1                           <font color=\"Red\">1<\/font><br \/>\n* MAT_IDS *<br \/>\n         15<br \/>\n* COIL_SERIES *<br \/>\n              1\n<\/p>\n<p>When calculating inductance with CALC_IND=1 or electromagnetic force with COIL_FORCE=1, the definition of integral elements is required. These definitions are done by LOOP-, GCE-, and ARC-. These have the same input format as LOOP, GCE, and ARC, respectively. However, the number of divisions and the number of Gaussian integral points must be entered at the end. These data are included in the coil definition. For example, the integral region in Fig. 1 is represented as follows. Since the upper half is defined here, CURRENT is also half. In this example, the number of coil turns is set to 100T.<br \/>\nNDIV (default value 1) is the number of divisions that divide the integral element equally in the current direction. Three-dimensional Gaussian integration is performed in the divided region. The electromagnetic force is output for each division; INT_X, INT_Y, and INT_Z are the number of integration points in each direction. The default values are 5, 5, and 3, respectively. INT_X and INT_Y are the number of integration points in the cross-sectional direction and INT_Z is the number of integration points in the current direction, with a maximum value of 12. The number of divisions and integration points required depends on the distribution of the magnetic field within the element, so increase the number of divisions and integration points as needed to make sure the results do not change. Since the electromagnetic forces are processed in POST_PROCESSING, you can change the number of integration points, etc., by performing post-processing in the restart calculation.  <\/p>\n<p class=\"slideText\">* GCE * CURRENT(A) *<br \/>\n   GCE          50<br \/>\n* XS(m) *  YS(m)  *  ZS(m)  *  XE(m)  *  YE(m)  *  ZE(m)  *<br \/>\n  0.0875      0.0          0.025      0.0875      0.05        0.025<br \/>\n*  W1X(m)  *  W1Y(m)  *  W1Z(m)  *  W2X(m)  *  W2Y(m)  *  W2Z(m)  *<br \/>\n      0.0                0.0           0.025          0.0125          0.0              0.0<br \/>\n* ARC * CURRENT(A)  *<br \/>\n   ARC            50<br \/>\n*  X(m)  *  Y(m)  *  Z(m)  *  RADIUS(m)  *  AXIAL_W(m)  *  RADIAL_W(m)  *<br \/>\n    0.05       0.05       0.025       0.0375                 0.05                 0.025<br \/>\n* ALPHA(deg) * BETA(deg) * PHI1(deg) * PHI2(deg) *<br \/>\n          0.0                  0.0                 0.0                 90<br \/>\n*  GCE  *  CURRENT(A)  *<br \/>\n   GCE            50<br \/>\n*  W1X(m)  *  W1Y(m)  *  W1Z(m)  *  W2X(m)  *  W2Y(m)  *  W2Z(m)  *<br \/>\n      0.05           0.0875        0.025          0.00           0.0875         0.025<br \/>\n*  W1X(m)  *  W1Y(m)  *  W1Z(m)  *  W2X(m)  *  W2Y(m)  *  W2Z(m)  *<br \/>\n      0.0                0.0            0.025            0.0           0.0125           0.0<br \/>\n<font color=\"Red\">*  GCE  *<\/font>  CURRENT(A)  *<br \/>\n<font color=\"Red\">   GCE-<\/font>         50<br \/>\n* XS(m) *  YS(m)  *  ZS(m)  *  XE(m)  *  YE(m)  *  ZE(m)  *<br \/>\n  0.0875      0.0          0.025      0.0875      0.05        0.025<br \/>\n*  W1X(m)  *  W1Y(m)  *  W1Z(m)  *  W2X(m)  *  W2Y(m)  *  W2Z(m)  *<br \/>\n      0.0                0.0           0.025          0.0125          0.0              0.0<br \/>\n<font color=\"Red\">*  NDIV * INT_X * INT_Y * INT_Z *<\/font><br \/>\n<font color=\"Red\">     1              5            5             3<\/font><br \/>\n<font color=\"Red\">* ARC *<\/font> CURRENT(A)  *<br \/>\n<font color=\"Red\">   ARC-<\/font>           50<br \/>\n*  X(m)  *  Y(m)  *  Z(m)  *  RADIUS(m)  *  AXIAL_W(m)  *  RADIAL_W(m)  *<br \/>\n    0.05       0.05       0.025       0.0375                 0.05                 0.025<br \/>\n* ALPHA(deg) * BETA(deg) * PHI1(deg) * PHI2(deg) *<br \/>\n          0.0                  0.0                 0.0                 90<br \/>\n<font color=\"Red\">*  NDIV * INT_X * INT_Y * INT_Z *<\/font><br \/>\n<font color=\"Red\">     1              5            5             3<\/font><br \/>\n<font color=\"Red\">*  GCE  *<\/font>  CURRENT(A)  *<br \/>\n<font color=\"Red\">   GCE-<\/font>           50<br \/>\n*  W1X(m)  *  W1Y(m)  *  W1Z(m)  *  W2X(m)  *  W2Y(m)  *  W2Z(m)  *<br \/>\n      0.05           0.0875        0.025          0.00           0.0875         0.025<br \/>\n*  W1X(m)  *  W1Y(m)  *  W1Z(m)  *  W2X(m)  *  W2Y(m)  *  W2Z(m)  *<br \/>\n      0.0                0.0            0.025            0.0           0.0125           0.0<br \/>\n<font color=\"Red\">*  NDIV * INT_X * INT_Y * INT_Z *<\/font><br \/>\n<font color=\"Red\">     1              5            5             3<\/font>\n<\/p>\n<p>As described above, when calculating the inductance and electromagnetic force of a COIL, the COIL information must be entered twice. On the one hand, the entire domain of the coil must be input as a conventional one, the coil that generates the magnetic field. On the other, the integral element does not create a magnetic field, but basically only the part within the finite element mesh region should be entered. For CIRCUIT or NETWORK, REGION<em>FACTOR is entered, and the inductance is integrated in the integral element, and using symmetry, REGION<\/em> FACTOR is multiplied. The electromagnetic force is integrated only within the divided integral elements and output to the output file without REGION_FACTOR multiplication. The number of integrating elements is also included in the NO_ELEMENTS of the COIL data. This function can be applied to LOOP, GCE, and ARC, but not to line current elements such as FGCE and FARC.  <\/p>\n<h3>Download<\/h3>\n<p><button type=\"button\" class=\"btn btn-danger btn-lg\"><a href=\"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en\/wp-content\/uploads\/coil-force.zip\">Sample data DL<\/a><\/button>  <\/p>\n<p>\u30fb inputSDEFCOIL<br \/>\n\u30fb inputCOIL<br \/>\n\u30fb pre_geom2D.neu<br \/>\n\u30fb 2D_to_3D  <\/p>\n<p>*This data cannot be run on the trial version due to the large number of nodes.  <\/p>\n","protected":false},"excerpt":{"rendered":"<p>Summary This section introduces how to calculate the inductance and electromagnetic force of an external current magnetic field source (COIL). Explanation COIL is a mesh-independent feature that is easy to use. An example of an axisymmetric model is presented in &quot;Inductance calculation for COIL (external current magnetic field source)&quot;, but here we will use a [&hellip;]<\/p>\n","protected":false},"featured_media":0,"template":"","tags":[],"case_cat":[6],"class_list":["post-126","case","type-case","status-publish","hentry","case_cat-gaibudenryu"],"acf":[],"aioseo_notices":[],"_links":{"self":[{"href":"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en\/wp-json\/wp\/v2\/case\/126"}],"collection":[{"href":"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en\/wp-json\/wp\/v2\/case"}],"about":[{"href":"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en\/wp-json\/wp\/v2\/types\/case"}],"version-history":[{"count":11,"href":"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en\/wp-json\/wp\/v2\/case\/126\/revisions"}],"predecessor-version":[{"id":2358,"href":"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en\/wp-json\/wp\/v2\/case\/126\/revisions\/2358"}],"wp:attachment":[{"href":"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en\/wp-json\/wp\/v2\/media?parent=126"}],"wp:term":[{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en\/wp-json\/wp\/v2\/tags?post=126"},{"taxonomy":"case_cat","embeddable":true,"href":"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en\/wp-json\/wp\/v2\/case_cat?post=126"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}