Accurate Force Computation in Magnetic Circuits

Physics of Magnetic Circuits

A magnetic circuit refers to a physical system comprised of components designed to generate and confine magnetic flux. Typically, this flux is generated by either permanent magnets or coils, with the flux paths being constructed from high-permeability ferromagnetic materials, such as iron. In Figure 1, we observe an illustration of a magnetic circuit [1] that includes two iron components with a specific cross-sectional areaS=4cm. 

One part is stationary and the other is moving; they are separated by two air gaps of length L. N-turn coil wrapped around the stationary part, carries the current I ( N=300; I=1A) and is responsible for the flux in the circuit.

Force on the moving part can be determined by calculating the change in the magnetic energy WW+1 that would be produced by moving the part over a small distance LL+1. Force magnitude is obtained as:

 F=ΔWΔL
   (eq.1)

This simple method is based upon the virtual displacement principle and it is often used to calculate forces in magnetic devices.
Since the permeability of iron is much greater than the permeability of air (μFeμ0
), the magnetic field in the iron (HFe=BμFe
), as well as leakage flux, are neglected. Ampere’s law for the contour in Figure 1. can be written as:

 2LHa=NI
    (eq.2)

where Math input error  represents the magnetic field in the air gap.

With no leakage and no field in the iron (HFe0
), the entire magnetic energy of the system is stored in the volume Math input error of the air gap and can be computed as:

 W=12VBHadV=12Vμ0Ha2dV
   (eq. 3)

Assuming uniform field distribution in the air gap allows to calculate W by simply multiplying energy density 12μ0Ha2
 and volume V. The circuit contains 2 air gaps, therefore, after doubling the energy in eq.3 and substituting eq.2 in eq.3:

 W=14μ0S(NI)2L
   (eq. 4)

For two different air gap lengths  Math input error  and Math input error , the change in magnetic energy can be computed as:

  ΔW=W2W1=14μ0S(NI)2L1L2L1L2
  (eq. 5)

Equation 5, combined with eq. 1, yields the following expression for force magnitude:

  F=14μ0S(NI)2L1L2=3.02N
    (eq. 6)


Magnetic circuit example
Figure 1 - Magnetic circuit example

MODEL

To accurately represent the situation in the analytical example, where force has been computed using energies for  Math input error  and Math input error air gaps, the CAD model of the circuit has been created with two air gaps of different lengths:Math input error and Math input error (Figure 2.). 
 

Model of one half of the magnetic circuit is created in SOLIDWORKS

Figure 2 - Model of one half of the magnetic circuit
 

To exploit the problem's inherent symmetry, it suffices to simulate just one-half of the space, corresponding to half of the magnetic circuit. The simulation is conducted using the EMS Magnetostatic study, with a symmetry factor set at 0.5.
Current current-driven Wound Coil with 300 turns and 1A/turn is added to the coil region. Its Entry and Exit Ports are faces in the plane of symmetry.

Material

The following instructions outline the process for creating a new material library, defining a custom material, and applying that material to an element within your model:

  1. In the EMS Manager tree, right-click the Stationary element icon in the Materials folder and select Apply Material to All Bodies. The Material Task Pane opens.
  2. Select “New Library nl  from the “Material Task Pane”.
  3. Browse to the location where to save the new library file (the file extension will be .emsmtr).
  4. Type MyLib for the name. An empty material library with the MyLib is added to the Material Task Panel.
  5. Click on Add Materiel button addm .
  6. Type Mur1200 for the material's name.
  7. Click OK ok .
  8. Type 1200 for the Relative Permeability and take the default for the rest of the fields.
  9. Click OK ok.

 Copper is assigned to the coil region, and air to the surrounding volume. 

Boundary Conditions

Tangential Flux boundary condition
Figure 3 - Tangential Flux boundary condition
 

The normal component of the flux density at the plane of symmetry ( xy plane in Figure 2.) is zero (all flux lines are parallel to this plane). Therefore, the Tangential Flux boundary condition should be applied on all surfaces that belong to the plane, including air and coil cross-sections. To do so:

  1. In the study, right-click the Load/Restraint l rfolder and select Tangential Flux. The Tangential Flux Property page appears.
  2. Click inside the Faces for Tangential Flux box t f box then select all the symmetry faces.
  3. Click OK ok .

Coils

To add a coil to a Magnetostatic study:

  1. In the EMS manger tree right-click on the Coils 1f icon and select Wound Coilff .
  2. Keep default Coil Type as a Current driven coil.
  3. Click inside the Components or Bodies for Coils box2f.
  4. Click on the (+) sign in the upper left corner of the graphics area to open the components tree.
  5. Click on the Coil-1 icon.  It will appear in the Components and Solid Bodies list.
  6. Click inside the Faces for Entry Port box 4f  then select the entry port face.
  7. Click inside the Faces for Exit Port box4f   then select the exit port face.

General Properties:

  1. Click on General properties tab.
  2. Type 300 in the Turns box5f.
  3. Keep the default value of 1 the Current per Turn  field 6f . The units in Amp-Turns.
  4. Click OK 7f .

 

Entry and exit ports of the coil
Figure 4 - Entry and exit ports of the coil
 

Force

EMS automatically calculates the distribution of nodal forces without requiring any user input. However, when it comes to calculating forces or torques on a rigid body, it is necessary to define the specific part on which the force or torque should be computed before initiating the simulation. 

  1. In the EMS manager tree, right-click the Force/Torques ft folder and select Virtual Work vw .
  2. The Forces/Torques Property Manager appears.
  3. Click inside the Components and Bodies for Forces/Torques box ftb .
  4. Click on the (+) sign in the upper left corner of the graphics area to open the components tree.
  5. Click on Moveable part icon.  It will appear in the Components and Solid Bodies list.
  6. Click OK ok.

Mesh

The quality of the mesh in the air gap region is of critical importance for accurate force calculation. EMS allows user to take full control over the mesh resolution.

  1. In the EMS manager tree, right-click the Mesh mesh icon and select Apply Control  ac . The Mesh Control Property Manager Page appears.   
  2. Click inside the Components and Solid Bodies box ftp .
  3. Click on the (+) sign in the upper left corner of the graphical area to open the components tree.
  4. Click on the Air Gap 1, and Air Gap 2 icons.  They will appear in the Components and Solid Bodies list.
  5. Under Control Parameters click inside the Element Size esbox and type 2.5 mm.
  6. Click OK ok .

To mesh the model:

  1. In the EMS manager tree, right-click on the Mesh mesh icon and select Create Mesh cm.
  2. Type 35 in the box for an average number of mesh elements per diagonal for each solid body sb.
  3. Click OK ok.

Right-click Study icon and select Run run, to execute the simulation. Once the computation is done, the program creates five folders in the EMS manager tree. These folders are Report, Magnetic Flux Density, Magnetic Field Intensity, Current Density, and Force Distribution.

RESULTS

Flux density

It is a good habit to first view the magnetic flux density in the model, including the outer air.  This action indicates whether the outer air boundary is far enough.

  1. In the EMS Manager tree, right-click on the Magnetic Flux Density folder and select 3D Fringe Plot> No Clipping.The Magnetic Flux Density Property Manager Page appears.
  2. In the Display box:

             a. Select Br from the magnetic flux density component br .  Directions are based on the global coordinate system.
             b. Set Units unitsto Tesla.
             c. Select Continuous from Fringe Options fo.

       3. Select OK ok.

Flux density in the magnetic circuit

 

Figure 5 - Flux density in the magnetic circuit
 

By examining Figure 5. it becomes clear that the magnetic flux density is very small on the outer air boundary.  Thus, the air box is large enough.  Had it been otherwise, the air box surrounding the magnetic circuit would have to be larger.

Force

In the EMS Manager tree, Double-click Result Table to display the force result.

 

 

EMS Results Table

 

Figure 6 - EMS Results Table
 

Remember that, because of the symmetry about the xy plane, only half of the problem has been modeled. Thus, Math input error and components must be multiplied by a factor of 2, while the Math input error component cancels out. Since Math input error is very small compared to Math input error , the resultant force is virtually in the X direction with a magnitude of Math input error . The analytical solution compares very well with the EMS result.

 

 Analytical Virtual Work SolutionEMS Result
Force [N]3.023.014

Conclusion

The application note delves into the physics and simulation of magnetic circuits, focusing on force computation within a specific configuration. Through analytical derivations and simulation in Solidworks EMS Magnetostatic study, the note demonstrates the calculation of forces using principles like Ampere's law and virtual displacement. Key components such as iron cores, air gaps, and coils are meticulously modeled, with materials assigned and boundary conditions set accordingly. The simulation results, including magnetic flux density and force distribution, are compared with analytical solutions, showcasing close agreement and validating the simulation approach. By accurately capturing the behavior of magnetic circuits, the note offers valuable insights for engineers working in electromagnetics and device design. Additionally, it provides practical guidance on modeling techniques and simulation setup. 

References

[1] Electromagnetics and calculation of fields, by Nathan Ida and Joao P. A. Bastos, 2nd Edition, page 183-184.  Publisher: Springer-Verlag;
 

Share on

Unlock Your Trial – See EMWORKS in Action

The best way to understand the impact of EMWORKS is to experience it yourself.
Fill out the form to access your trial and let us help you make the most of it—because seeing is believing.

  • United States+1
  • United Kingdom+44
  • Afghanistan (‫افغانستان‬‎)+93
  • Albania (Shqipëri)+355
  • Algeria (‫الجزائر‬‎)+213
  • American Samoa+1
  • Andorra+376
  • Angola+244
  • Anguilla+1
  • Antigua and Barbuda+1
  • Argentina+54
  • Armenia (Հայաստան)+374
  • Aruba+297
  • Australia+61
  • Austria (Österreich)+43
  • Azerbaijan (Azərbaycan)+994
  • Bahamas+1
  • Bahrain (‫البحرين‬‎)+973
  • Bangladesh (বাংলাদেশ)+880
  • Barbados+1
  • Belarus (Беларусь)+375
  • Belgium (België)+32
  • Belize+501
  • Benin (Bénin)+229
  • Bermuda+1
  • Bhutan (འབྲུག)+975
  • Bolivia+591
  • Bosnia and Herzegovina (Босна и Херцеговина)+387
  • Botswana+267
  • Brazil (Brasil)+55
  • British Indian Ocean Territory+246
  • British Virgin Islands+1
  • Brunei+673
  • Bulgaria (България)+359
  • Burkina Faso+226
  • Burundi (Uburundi)+257
  • Cambodia (កម្ពុជា)+855
  • Cameroon (Cameroun)+237
  • Canada+1
  • Cape Verde (Kabu Verdi)+238
  • Caribbean Netherlands+599
  • Cayman Islands+1
  • Central African Republic (République centrafricaine)+236
  • Chad (Tchad)+235
  • Chile+56
  • China (中国)+86
  • Christmas Island+61
  • Cocos (Keeling) Islands+61
  • Colombia+57
  • Comoros (‫جزر القمر‬‎)+269
  • Congo (DRC) (Jamhuri ya Kidemokrasia ya Kongo)+243
  • Congo (Republic) (Congo-Brazzaville)+242
  • Cook Islands+682
  • Costa Rica+506
  • Côte d’Ivoire+225
  • Croatia (Hrvatska)+385
  • Cuba+53
  • Curaçao+599
  • Cyprus (Κύπρος)+357
  • Czech Republic (Česká republika)+420
  • Denmark (Danmark)+45
  • Djibouti+253
  • Dominica+1
  • Dominican Republic (República Dominicana)+1
  • Ecuador+593
  • Egypt (‫مصر‬‎)+20
  • El Salvador+503
  • Equatorial Guinea (Guinea Ecuatorial)+240
  • Eritrea+291
  • Estonia (Eesti)+372
  • Ethiopia+251
  • Falkland Islands (Islas Malvinas)+500
  • Faroe Islands (Føroyar)+298
  • Fiji+679
  • Finland (Suomi)+358
  • France+33
  • French Guiana (Guyane française)+594
  • French Polynesia (Polynésie française)+689
  • Gabon+241
  • Gambia+220
  • Georgia (საქართველო)+995
  • Germany (Deutschland)+49
  • Ghana (Gaana)+233
  • Gibraltar+350
  • Greece (Ελλάδα)+30
  • Greenland (Kalaallit Nunaat)+299
  • Grenada+1
  • Guadeloupe+590
  • Guam+1
  • Guatemala+502
  • Guernsey+44
  • Guinea (Guinée)+224
  • Guinea-Bissau (Guiné Bissau)+245
  • Guyana+592
  • Haiti+509
  • Honduras+504
  • Hong Kong (香港)+852
  • Hungary (Magyarország)+36
  • Iceland (Ísland)+354
  • India (भारत)+91
  • Indonesia+62
  • Iran (‫ایران‬‎)+98
  • Iraq (‫العراق‬‎)+964
  • Ireland+353
  • Isle of Man+44
  • Israel (‫ישראל‬‎)+972
  • Italy (Italia)+39
  • Jamaica+1
  • Japan (日本)+81
  • Jersey+44
  • Jordan (‫الأردن‬‎)+962
  • Kazakhstan (Казахстан)+7
  • Kenya+254
  • Kiribati+686
  • Kosovo+383
  • Kuwait (‫الكويت‬‎)+965
  • Kyrgyzstan (Кыргызстан)+996
  • Laos (ລາວ)+856
  • Latvia (Latvija)+371
  • Lebanon (‫لبنان‬‎)+961
  • Lesotho+266
  • Liberia+231
  • Libya (‫ليبيا‬‎)+218
  • Liechtenstein+423
  • Lithuania (Lietuva)+370
  • Luxembourg+352
  • Macau (澳門)+853
  • Macedonia (FYROM) (Македонија)+389
  • Madagascar (Madagasikara)+261
  • Malawi+265
  • Malaysia+60
  • Maldives+960
  • Mali+223
  • Malta+356
  • Marshall Islands+692
  • Martinique+596
  • Mauritania (‫موريتانيا‬‎)+222
  • Mauritius (Moris)+230
  • Mayotte+262
  • Mexico (México)+52
  • Micronesia+691
  • Moldova (Republica Moldova)+373
  • Monaco+377
  • Mongolia (Монгол)+976
  • Montenegro (Crna Gora)+382
  • Montserrat+1
  • Morocco (‫المغرب‬‎)+212
  • Mozambique (Moçambique)+258
  • Myanmar (Burma) (မြန်မာ)+95
  • Namibia (Namibië)+264
  • Nauru+674
  • Nepal (नेपाल)+977
  • Netherlands (Nederland)+31
  • New Caledonia (Nouvelle-Calédonie)+687
  • New Zealand+64
  • Nicaragua+505
  • Niger (Nijar)+227
  • Nigeria+234
  • Niue+683
  • Norfolk Island+672
  • North Korea (조선 민주주의 인민 공화국)+850
  • Northern Mariana Islands+1
  • Norway (Norge)+47
  • Oman (‫عُمان‬‎)+968
  • Pakistan (‫پاکستان‬‎)+92
  • Palau+680
  • Palestine (‫فلسطين‬‎)+970
  • Panama (Panamá)+507
  • Papua New Guinea+675
  • Paraguay+595
  • Peru (Perú)+51
  • Philippines+63
  • Poland (Polska)+48
  • Portugal+351
  • Puerto Rico+1
  • Qatar (‫قطر‬‎)+974
  • Réunion (La Réunion)+262
  • Romania (România)+40
  • Russia (Россия)+7
  • Rwanda+250
  • Saint Barthélemy+590
  • Saint Helena+290
  • Saint Kitts and Nevis+1
  • Saint Lucia+1
  • Saint Martin (Saint-Martin (partie française))+590
  • Saint Pierre and Miquelon (Saint-Pierre-et-Miquelon)+508
  • Saint Vincent and the Grenadines+1
  • Samoa+685
  • San Marino+378
  • São Tomé and Príncipe (São Tomé e Príncipe)+239
  • Saudi Arabia (‫المملكة العربية السعودية‬‎)+966
  • Senegal (Sénégal)+221
  • Serbia (Србија)+381
  • Seychelles+248
  • Sierra Leone+232
  • Singapore+65
  • Sint Maarten+1
  • Slovakia (Slovensko)+421
  • Slovenia (Slovenija)+386
  • Solomon Islands+677
  • Somalia (Soomaaliya)+252
  • South Africa+27
  • South Korea (대한민국)+82
  • South Sudan (‫جنوب السودان‬‎)+211
  • Spain (España)+34
  • Sri Lanka (ශ්‍රී ලංකාව)+94
  • Sudan (‫السودان‬‎)+249
  • Suriname+597
  • Svalbard and Jan Mayen+47
  • Swaziland+268
  • Sweden (Sverige)+46
  • Switzerland (Schweiz)+41
  • Syria (‫سوريا‬‎)+963
  • Taiwan (台灣)+886
  • Tajikistan+992
  • Tanzania+255
  • Thailand (ไทย)+66
  • Timor-Leste+670
  • Togo+228
  • Tokelau+690
  • Tonga+676
  • Trinidad and Tobago+1
  • Tunisia (‫تونس‬‎)+216
  • Turkey (Türkiye)+90
  • Turkmenistan+993
  • Turks and Caicos Islands+1
  • Tuvalu+688
  • U.S. Virgin Islands+1
  • Uganda+256
  • Ukraine (Україна)+380
  • United Arab Emirates (‫الإمارات العربية المتحدة‬‎)+971
  • United Kingdom+44
  • United States+1
  • Uruguay+598
  • Uzbekistan (Oʻzbekiston)+998
  • Vanuatu+678
  • Vatican City (Città del Vaticano)+39
  • Venezuela+58
  • Vietnam (Việt Nam)+84
  • Wallis and Futuna (Wallis-et-Futuna)+681
  • Western Sahara (‫الصحراء الغربية‬‎)+212
  • Yemen (‫اليمن‬‎)+967
  • Zambia+260
  • Zimbabwe+263
  • Åland Islands+358
By visiting or using this website you agree to the storing of cookies on your device to enhance site navigation, analyze site usage, and assist in our marketing efforts. View cookies details.