Difference between revisions of "Homework Assignment"

From Maxwell
Jump to: navigation, search
(Calculating the force on the solenoid valve by finite element method)
(Part II of the Assignment)
 
(86 intermediate revisions by the same user not shown)
Line 2: Line 2:
 
|-
 
|-
 
| colspan=2 align=center |
 
| colspan=2 align=center |
<font color='blue' size='+2'>Analysis of Solenoid Valve</font>
+
<font color='blue' size='+2'>Analysis of Shunt Resistor</font>
 
{| align="left"
 
{| align="left"
 
  | __TOC__
 
  | __TOC__
Line 9: Line 9:
 
| width=50% |
 
| width=50% |
 
'''Instructor'''
 
'''Instructor'''
* Dániel Marcsa (lecturer)
+
* [http://wiki.maxwell.sze.hu/index.php/Marcsa Dániel Marcsa] (lecturer)
* Lectures: Monday, 14:50 - 16:25 (D201), 16:30 - 17:15 (D105)
+
* Lectures: -
 
* Office hours: by request
 
* Office hours: by request
 
| width=50% |
 
| width=50% |
Line 23: Line 23:
 
=== Knowledge needed to solve the problem ===
 
=== Knowledge needed to solve the problem ===
 
* The main steps of the finite element method;
 
* The main steps of the finite element method;
* Theoretical knowledge of the static magnetic field (for defining materials, for excitation);
+
* Theoretical knowledge of the electromagnetic field simulation (for defining materials, for excitation);
 
* Knowledge of CAD system to create geometry;
 
* Knowledge of CAD system to create geometry;
* Download and install [https://www.ansys.com/academic/free-student-products ANSYS AIM Student].
+
* Download and install [https://www.ansys.com/academic/students/ansys-electronics-desktop-student Ansys Electronics Desktop Student].
  
 
=== The Semester Assignment ===
 
=== The Semester Assignment ===
  
The task consists of two parts, a basic task, with a faultless solution of up to 70%, and an extra task, with an additional maximum of 30%.
+
The task consists of two parts, a basic task, with a faultless solution of up to 80%, and an extra task, with an additional maximum of 20%.
  
 
{| class = "wikitable"
 
{| class = "wikitable"
 
| '''Deadline'''
 
| '''Deadline'''
| May 25, 2019, 12:00
+
| -
 
|-
 
|-
 
| '''Output Form:'''
 
| '''Output Form:'''
Line 42: Line 42:
 
|-
 
|-
 
| '''Place of submission:'''
 
| '''Place of submission:'''
| marcsa@maxwell.sze.hu email address
+
| <span style="color:red">In [https://szelearning.sze.hu/sze-login/index.html '''Moodle system'''].</span>
 
|-
 
|-
 
| '''Late submission:'''
 
| '''Late submission:'''
Line 61: Line 61:
 
|}
 
|}
 
</blockquote>
 
</blockquote>
 +
 
== Part I of the Assignment ==
 
== Part I of the Assignment ==
 
<blockquote>
 
<blockquote>
==== Calculating the force on the solenoid valve by finite element method ====
+
==== Calculating the resistance and the total loss of the shunt resistor by finite element method ====
[[File: ProblemGeometry.png | 360px | thumb | left | alt = The geometry and dimensions of the task. | The geometry and dimensions of the task.]]
+
[[File: CAE_SA_ShuntResistor_2022Spring.png | 800px | thumb | left | alt = ''Figure 1.'' - The geometry and dimensions of the task. | The geometry and dimensions of the task.]]
The task is cylindrical to the vertical (<math>z</math>) axis, and the three-dimensional geometry must be prepared according to the specified dimensions (see figure).
+
 
 +
The geometry dimensions for your task you can find in the following table: [https://docs.google.com/spreadsheets/d/1cieVgmu9ALZBLYRZyxL-eHe27nOpavBuOGftaLXCC9I/edit?usp=sharing '''Semester Assigment'''].
 +
 
 +
This task is a DC current conduction problem. The solved equation is
 +
 
 +
<math>\nabla\cdot\sigma\nabla \varphi=0</math>
 +
 
 +
with following boundary conditions
 +
 
 +
<math>\vec{J}\cdot\vec{n}=-J_{\text{n}}</math> on <math>\Gamma_{\text{J}}</math> (''This is the input.'')
 +
 
 +
and
 +
 
 +
<math>\varphi=U_0 = \text{0 V}</math> on <math>\Gamma_{\text{E}}</math> (''This is the output.''),
 +
 
 +
where <math>J_{\text{n}}</math> is the current density calculated from the specified current excitation.
 +
 
 +
The task: determine the voltage drop, the resistance and the ohmic loss of the problem.<br \>
 +
The voltage drop is the potential difference between the two terminals of the arrangement. You can determine resistance using Ohm's law:
 +
 +
<math>R = \frac{U}{I}</math>,
 +
 
 +
then the ohmic loss
 +
 +
<math>P = I^2\cdot R</math>
 +
 
 +
where <math>U</math> is the voltage drop, <math>I</math> is the current and <math>R</math> is the resistance.
 +
 
 +
{| class="wikitable"  style="text-align: center; width: 1000px; height: 80px;"
 +
|+ Bulk conductivity of materials.
 +
! Material
 +
! Titanium
 +
! Copper
 +
! Aluminum
 +
! Copper-manganin alloy
 +
|-
 +
! <math>\sigma~[\text{MS/m}]</math>
 +
| 1.82 || 58 || 38 || 20.833
 +
|}
  
 
'''Tasks'''
 
'''Tasks'''
 +
* Draw the geometry based on the specified dimensions in [https://www.ansys.com/academic/students/ansys-electronics-desktop-student Ansys Electronics Desktop Student];
 +
* Define the problem based on the given material parameters and boundary conditions;
 +
* Run the FEM simulation;
 +
* Evaluate the results.
 +
 +
The quantities listed in the task can be calculated with the '''Maxwell 3D - DC Conduction''' solver and the '''Q3D Extractor - DC''' solver.
 +
 +
{| width = 100%
 +
|-
 +
| align = center |
 +
[[Image: CAE SA EMSolutions 2022Spring.png | 650px]]
 +
|-
 +
| align = center | <span style="font-size:88%;>''' ''Figure 2.'' - Possible solution (''Left - Maxwell 3D, Right - Q3D Extractor).'''</span>
 +
|}
  
* Define the problem type based on the given parameters;
+
{| class="wikitable"  style="text-align: center; width: 1000px; height: 80px;"
* Creating and specifying the task geometry in [https://www.ansys.com/products/3d-design/ansys-aim ANSYS Discovery AIM];
+
|+ Results of the test example.
* Run the FEM simulation;
+
! Quantities
* Validation of results for a 2mm air gap.
+
! Voltage drop [mV]
{| class="wikitable"  style="text-align: center; width: 600px; height: 100px;"
+
! Resistance [<math>\text{n}\Omega</math>]
|+ The results of the simulations.
+
! Ohmic loss [W]
! Software
+
|-
! Discovery AIM
 
 
! Maxwell 3D
 
! Maxwell 3D
! FEMM
+
| 10.9927 || 18.3211 || 6.5956
! Maxwell 2D
+
|-
 +
! Q3D Extractor
 +
| 10.9759 || 18.2931 || 6.5855
 +
|}
 +
 
 +
{| width=100%
 
|-
 
|-
! Force [N]
+
| align=center |
| 3,586 || 3,582 || 3,542 || 3,587
+
[[File:CAE SA ShuntResistor Maxwell Loss.png|750px]]
 +
| align=center |
 +
[[File:CAE SA ShuntResistor Q3D Loss.png|750px]]
 
|-
 
|-
! Inductance [mH]
+
|align=center | <span style="font-size:88%;">'''''[https://www.ansys.com/products/electronics/ansys-maxwell Ansys Maxwell 3D]'' - Ohmic loss on the surface of shunt resistor.'''</span>
| 39,88
+
|align=center | <span style="font-size:88%;">'''''[https://www.ansys.com/products/electronics/ansys-q3d-extractor Ansys Q3D Extractor]'' - Ohmic loss on the surface of shunt resistor.'''</span>
| 39,84
 
| 39,71
 
| 39,84
 
 
|}
 
|}
  
* Test of Solution Performance (''Solution performance tuning'', ''Curved surface meshing'') as a function of force;
+
'''Step by step tutorials'''
* Number of finite elements (''Tetrahedra''), Energy (''Total Energy''), ''Energy Error'' and ''Delta Energy'' in the function of adaptive steps in two different cases;
+
* [https://drive.google.com/file/d/1M1zkJJA5_rUskylsV1hZ6FXvjUOcesgv/view?usp=sharing Geometry creation]
 +
* [https://drive.google.com/file/d/1-opGHP6VoefJtuioH8osPnr_zXM0aWeN/view?usp=sharing Ansys Maxwell3D]
 +
* [https://drive.google.com/file/d/1ln4kZIZzvCVbRTD-QmCvwI6-trYigTGm/view?usp=sharing Ansys Q3D Extractor]
  
* Postprocessing [at least two of the below]:
+
</blockquote>
:: 1) determination of force and inductance;
 
:: 2) displaying equipotential lines;
 
:: 3) display the minimum and maximum of the magnetic flux density;
 
:: 4) display the magnetic field strength;
 
:: 5) displaying magnetic flux density vectors.
 
  
* Preparation of a Technical Report based on the use of the above results.
+
== Part II of the Assignment ==
 +
<blockquote>
 +
In the student version of Ansys EM, Icepak provides an opportunity to study the thermal phenomena of the task. '''[https://www.ansys.com/products/electronics/ansys-icepak Ansys Icepak]''' is a general CFD solver with specific capabilities for testing the heating and cooling of electronic circuits (PCB / power module).
  
Note: The “Part I of the Assignment” elements will be reviewed in detail during the exercise to prepare and run a simulation model for an electromagnetic task. On this basis, students can easily complete Part I of the Assignment by attending the exercises.
+
The task has a natural convection cooling. The excitation is the ohmic loss from the electromagnetic simulation.
  
'''Note:''' <span style="color:red">''A feladatban a vasmag ('''Core''') és a tekercs ('''Coil''') érintkezik. Ez kétdimenziós (hengerszimmetrikus) feladatnál nem okoz problémát, de a háromdimenziós példánál emiatt az áram nem csak a tekercsben hanem a vasmagban is folyik. Ennek kiküszöbölésére az ANSYS Discovery AIM szoftverben a '''Physics Definition - Electromagnetic Conditions - Insulating''' peremfeltételt kell alkalmazni a tekercsre.''</span> A valóságban is szigetelő veszi körül a tekercset, amit jellemzően nem modellezünk a feladatban, helyette peremfeltételt alkalmazunk.
+
{| class="wikitable"  style="text-align: center; width: 1000px; height: 80px;"
 +
|+ Thermal properties of materials.
 +
! Material
 +
! Titanium
 +
! Copper
 +
! Aluminum
 +
! Copper-manganin alloy
 +
|-
 +
! <math>\rho~[\text{kg}/\text{m}^3]</math>
 +
| 4500 || 8933 || 2689 || 8400
 +
|-
 +
! <math>c_{\text{P}}~[\text{J}/(\text{kg}\cdot\text{°C})]</math>
 +
| 522 || 385 || 951 || 410
 +
|-
 +
! <math>\lambda~[\text{W}/(\text{m}\cdot\text{°C})]</math>
 +
| 21 || 400 || 237.5 || 22
 +
|}
  
; '''Parameters for this Assignment:'''
+
The table below shows the result of the thermal simulation.
: Relative permeability of air and coil <math>\mu_r = 1</math>;
+
 
: The excitation of the coil ('''To validate the results'''): <math>I=0.76~\text{A}</math>, <math>N=789~\text{turns}</math> (''DC'');
+
{| class="wikitable" style="text-align: center; width: 1000px; height: 80px;"
: Magnetic curve of the core and plunge:
+
|+ Results of the test example.
: {| class = "wikitable" style = "text-align: center;"
+
! Quantities
{| width = 100%
+
! Max. temperature [°C]
 +
! Min. temperature [°C]
 +
! Max. velocity [m/s]
 +
|-
 +
! Maxwell 3D + Icepak
 +
| 75.9284 || 71.3658 || 0.2881
 +
|-
 +
! Q3D Extractor + Icepak
 +
| 75.8669 || 71.3002 || 0.2880
 +
|}
 +
 
 +
{| width=100%
 
|-
 
|-
| align = center |
+
| align=center |
[[Image: BHCurve Steel1008.PNG | 650px]]
+
[[File:CAE SA ShuntResistor Maxwell Temp.png|750px]]
 +
| align=center |
 +
[[File:CAE SA ShuntResistor Q3D Temp.png|750px]]
 
|-
 
|-
| align = center | <span style="font-size:88%;>''' ''Figure 1.'' - Steel 1008 steel magnetization curve.'''</span>
+
|align=center | <span style="font-size:88%;">''' ''[https://www.ansys.com/products/electronics/ansys-maxwell Ansys Maxwell 3D] + [https://www.ansys.com/products/electronics/ansys-icepak Ansys Icepak]'' - Temperature distribution on the surface of shunt resistor.'''</span>
 +
|align=center | <span style="font-size:88%;">''' ''[https://www.ansys.com/products/electronics/ansys-q3d-extractor Ansys Q3D Extractor] + [https://www.ansys.com/products/electronics/ansys-icepak Ansys Icepak]'' - Temperature distribution on the surface of shunt resistor.'''</span>
 
|}
 
|}
</blockquote>
 
  
== Part II of the Assignment ==
+
'''Step by step tutorials'''
<blockquote>
+
* [https://drive.google.com/file/d/11b-sW497e8b5vkleeBU70HGv4AGSmg9u/view?usp=sharing Ansys Icepak]
Determining the force acting on the moving part of the problem ('''plunger''') and the inductance of the coil as a function of displacement.
+
 
 +
Here you can find the archive file of test example: [https://drive.google.com/file/d/1-lZzCZwFssZxERUNK1P4iHuyRgSSROd3/view?usp=sharing Shunt Resistor] (''Ansys EM Student 2021 R2'').
  
Note: This task corresponds to the solenoid valve in internal combustion engines that controls the injection.
 
; Specific tasks:
 
* Considering the plunger movement to the initial state (2mm air gap) from <math>-1.8\text{mm} \text{ to } 10\text{mm}</math> ('''at least''' in 7 positions). Detailed instructions on the solution will not be available.
 
* This part is intended to measure the degree of autonomy, initiative and diligence of the student, i.e.:
 
*# is the student able to do independent work;
 
*# to design, assemble and run the simulation alone.
 
 
</blockquote>
 
</blockquote>
  
 
==References==
 
==References==
 
{{reflist}}
 
{{reflist}}

Latest revision as of 17:09, 20 March 2022

Analysis of Shunt Resistor

Instructor

Teaching Assistants:

  • -
  • Office hours: -.

Aim of Assignment

Students will learn the basics of electromagnetic field calculations, their main steps, and gain practice in evaluating results and writing a Technical Report that meets international expectations.

Knowledge needed to solve the problem

  • The main steps of the finite element method;
  • Theoretical knowledge of the electromagnetic field simulation (for defining materials, for excitation);
  • Knowledge of CAD system to create geometry;
  • Download and install Ansys Electronics Desktop Student.

The Semester Assignment

The task consists of two parts, a basic task, with a faultless solution of up to 80%, and an extra task, with an additional maximum of 20%.

Deadline -
Output Form: PDF format. Color drawings should be made so that their contents are clear to the reader in black and white.
Language English
Place of submission: In Moodle system.
Late submission: After every day started, a 5% deduction from the achieved result (i.e. after 5 days delay 100% can only be obtained up to 100% - 5x5% = 75%).
Evaluation: 0 - 48% - Insufficient [F] (1)
50 - 59% - Sufficient [D] (2)
60 - 70% Satisfactory [C] (3)
71 - 84% Good [B] (4)
85 - 100% - Very good [A] (5)
For the formal requirements, the requirements of CFD and mechanics are also valid here.

Part I of the Assignment

Calculating the resistance and the total loss of the shunt resistor by finite element method

The geometry and dimensions of the task.

The geometry dimensions for your task you can find in the following table: Semester Assigment.

This task is a DC current conduction problem. The solved equation is

[math]\nabla\cdot\sigma\nabla \varphi=0[/math]

with following boundary conditions

[math]\vec{J}\cdot\vec{n}=-J_{\text{n}}[/math] on [math]\Gamma_{\text{J}}[/math] (This is the input.)

and

[math]\varphi=U_0 = \text{0 V}[/math] on [math]\Gamma_{\text{E}}[/math] (This is the output.),

where [math]J_{\text{n}}[/math] is the current density calculated from the specified current excitation.

The task: determine the voltage drop, the resistance and the ohmic loss of the problem.
The voltage drop is the potential difference between the two terminals of the arrangement. You can determine resistance using Ohm's law:

[math]R = \frac{U}{I}[/math],

then the ohmic loss

[math]P = I^2\cdot R[/math]

where [math]U[/math] is the voltage drop, [math]I[/math] is the current and [math]R[/math] is the resistance.

Bulk conductivity of materials.
Material Titanium Copper Aluminum Copper-manganin alloy
[math]\sigma~[\text{MS/m}][/math] 1.82 58 38 20.833

Tasks

  • Draw the geometry based on the specified dimensions in Ansys Electronics Desktop Student;
  • Define the problem based on the given material parameters and boundary conditions;
  • Run the FEM simulation;
  • Evaluate the results.

The quantities listed in the task can be calculated with the Maxwell 3D - DC Conduction solver and the Q3D Extractor - DC solver.

CAE SA EMSolutions 2022Spring.png

Figure 2. - Possible solution (Left - Maxwell 3D, Right - Q3D Extractor).
Results of the test example.
Quantities Voltage drop [mV] Resistance [[math]\text{n}\Omega[/math]] Ohmic loss [W]
Maxwell 3D 10.9927 18.3211 6.5956
Q3D Extractor 10.9759 18.2931 6.5855

CAE SA ShuntResistor Maxwell Loss.png

CAE SA ShuntResistor Q3D Loss.png

Ansys Maxwell 3D - Ohmic loss on the surface of shunt resistor. Ansys Q3D Extractor - Ohmic loss on the surface of shunt resistor.

Step by step tutorials

Part II of the Assignment

In the student version of Ansys EM, Icepak provides an opportunity to study the thermal phenomena of the task. Ansys Icepak is a general CFD solver with specific capabilities for testing the heating and cooling of electronic circuits (PCB / power module).

The task has a natural convection cooling. The excitation is the ohmic loss from the electromagnetic simulation.

Thermal properties of materials.
Material Titanium Copper Aluminum Copper-manganin alloy
[math]\rho~[\text{kg}/\text{m}^3][/math] 4500 8933 2689 8400
[math]c_{\text{P}}~[\text{J}/(\text{kg}\cdot\text{°C})][/math] 522 385 951 410
[math]\lambda~[\text{W}/(\text{m}\cdot\text{°C})][/math] 21 400 237.5 22

The table below shows the result of the thermal simulation.

Results of the test example.
Quantities Max. temperature [°C] Min. temperature [°C] Max. velocity [m/s]
Maxwell 3D + Icepak 75.9284 71.3658 0.2881
Q3D Extractor + Icepak 75.8669 71.3002 0.2880

CAE SA ShuntResistor Maxwell Temp.png

CAE SA ShuntResistor Q3D Temp.png

Ansys Maxwell 3D + Ansys Icepak - Temperature distribution on the surface of shunt resistor. Ansys Q3D Extractor + Ansys Icepak - Temperature distribution on the surface of shunt resistor.

Step by step tutorials

Here you can find the archive file of test example: Shunt Resistor (Ansys EM Student 2021 R2).

References