Hi Elmer-community!
I want to simulate the force when a rotor is pulled out of a stator. To do so i want to use its symmetry. My Model is a 2-pole motor. For this reason I cut the full Model into a half model and apply anti-periodic radial Mortar
EDIT: I removed the here provided files, as there were some errors. Please see my problem description in the following post below:
viewtopic.php?p=30311#p30311
The simulation files can be found here:
viewtopic.php?p=30314#p30314
Simulation of Magnetic Force with a half 3D Motormodel
Simulation of Magnetic Force with a half 3D Motormodel
Last edited by FFF on 17 Apr 2024, 09:45, edited 3 times in total.
Re: Simulation of Magnetic Force with a half 3D Motormodel
Last edited by FFF on 17 Apr 2024, 09:40, edited 1 time in total.
Re: Simulation of Magnetic Force with a half 3D Motormodel
Hi,
I highly appreciate, if someone could have a look at the following corrected files and give me hints on my questions below. ( I had to split up my message into separate posts, because i can only upload three files per post...)
I figured out, that the files I provided in my initial post are difficult to reproduce, because I manually had to modify them... (i.e. I zipped the wrong files)! Sorry for the confusion!
Just run the new files (provided in the next post) in the terminal with the below given (example-)commands and you should get exactly my output/results.
If you copy the commands into a shell script you can use: time ./myshellscript.sh in the terminal to measure the total time of the simulation process. With this method I compared the computational time of the symmetric (half) models with the full model.
These are my meshed models where...
...the full_Motormodel consist of 3.6 Million Elements (note, I cut the view (i.e. Clipping in Salome) of the full_Motormodel for a better visualization of its mesh) The model half_motor_mortarBC_on_all_subfaces consist of 1.826 Million Elements The model mesh_half_motor_mortarBC_on_merged_subfaces consist of 1.83 Million Elements
The difference between the two half motor models will be clarified in the next post.
I highly appreciate, if someone could have a look at the following corrected files and give me hints on my questions below. ( I had to split up my message into separate posts, because i can only upload three files per post...)
I figured out, that the files I provided in my initial post are difficult to reproduce, because I manually had to modify them... (i.e. I zipped the wrong files)! Sorry for the confusion!
Just run the new files (provided in the next post) in the terminal with the below given (example-)commands and you should get exactly my output/results.
Code: Select all
time salome -t -w1 full_Motormodel.py
time ElmerGrid 8 2 full_Motormodel.unv >OutputElmerGrid.log
time ElmerSolver full_Motormodel.sif >OutputElmerSolver.log
These are my meshed models where...
...the full_Motormodel consist of 3.6 Million Elements (note, I cut the view (i.e. Clipping in Salome) of the full_Motormodel for a better visualization of its mesh) The model half_motor_mortarBC_on_all_subfaces consist of 1.826 Million Elements The model mesh_half_motor_mortarBC_on_merged_subfaces consist of 1.83 Million Elements
The difference between the two half motor models will be clarified in the next post.
Re: Simulation of Magnetic Force with a half 3D Motormodel
An overview of the full_Motormodel geometry is show below.
The difference between the two half motor models is how I applied the Mortar BC.
In my first approach I merged all subfaces of the cut into two separate faces In my second approach I applied for each subface a individual Mortar BC Please see the results and my questions in the next post.
The difference between the two half motor models is how I applied the Mortar BC.
In my first approach I merged all subfaces of the cut into two separate faces In my second approach I applied for each subface a individual Mortar BC Please see the results and my questions in the next post.
Re: Simulation of Magnetic Force with a half 3D Motormodel
The names of the scalar values are:
The output of the full_Motormodel
The output of the mesh_half_motor_mortarBC_on_merged_subfaces
The output of the half_motor_mortarBC_on_all_subfaces
I know from Peters answer here (https://www.elmerfem.org/forum/viewtopi ... 303#p30303), that the magnitude of the results of the symmetric models are according to their geometrical fractions. Therfeore, comparing the results of the field energy (column 2) and the Axial Force (column 5), it can be seen that the magnitudes of the model half_motor_mortarBC_on_all_subfaces as well as half_motor_mortarBC_on_merged_subfaces do coincide quit well (at least for the former model) with the full_Motormodel. But comparing the two radial forces (column 3, 4) and the torque (column 6), it can be seen that this is not the case.
1. What could be the reason for this? Can the radial forces be neglected because they vectorially add due to the symmetric geometry? What is the case of the torque?
The computation time of the simulation process (total = automatic generation of the geometry and mesh in salome + export of .unv + transformation with ElmerGrid + computation of ElmerSolver) were:
2. What could be the reason that the simulation time of the full_motormodel is almost half of the half_motor_mortarBC_on_all_subfaces?
full_Motormodel.sif
OutputElmerSolver.log of full_Motormodel
half_motor_mortarBC_on_merged_subfaces.sif
Observing the OutputElmerSolver.log of half_motor_mortarBC_on_merged_subfaces.sif, it can be seen that the solver had difficulties to converge, as almost always the Number of max. Iterations were needed.
half_motor_mortarBC_on_all_subfaces.sif
Observing the OutputElmerSolver.log of half_motor_mortarBC_on_all_subfaces.sif, it can be seen that in this simulation the solver had to use the max number of itterations in the first 6 loops.
Comparing the B-Field distribution, only the full_motormodel do match my expectations. See the Figures below.
3. What could be the reason for the different field distribution?
You can download the simulation models in the next post.
Code: Select all
Variables in columns of matrix:
1: res: eddy current power
2: res: electromagnetic field energy
3: component 1: res: magnetic force 1
4: component 1: res: magnetic force 2
5: component 1: res: magnetic force 3
6: component 1: res: magnetic torque
Code: Select all
0.000000000000E+000 6.482809757103E+000 1.310099687527E-002 4.733919862382E-002 -8.146031875554E+001 -1.770470748866E-003
Code: Select all
0.000000000000E+000 3.285102780709E+000 1.797453227471E-001 -1.048604815010E+002 -3.971418110289E+001 1.039687840910E-003
The output of the half_motor_mortarBC_on_all_subfaces
Code: Select all
0.000000000000E+000 3.253456763210E+000 -1.608737810012E-001 -9.685255289219E+001 -4.077288817520E+001 6.147278705764E-003
1. What could be the reason for this? Can the radial forces be neglected because they vectorially add due to the symmetric geometry? What is the case of the torque?
The computation time of the simulation process (total = automatic generation of the geometry and mesh in salome + export of .unv + transformation with ElmerGrid + computation of ElmerSolver) were:
- full_motormodel: 57min in total and and ElmerSolver only was 49min
- half_motor_mortarBC_on_merged_subfaces: 141min in total and ElmerSolver only was 137min
- half_motor_mortarBC_on_all_subfaces: 87min in total and ElmerSolver only was 83min
2. What could be the reason that the simulation time of the full_motormodel is almost half of the half_motor_mortarBC_on_all_subfaces?
full_Motormodel.sif
Code: Select all
Check Keywords "Warn"
Header
Mesh DB "." "full_Motormodel"
Results Directory "result_full_Motormodel"
End
Simulation
Max Output Level = 3
Coordinate System = "Cartesian"
Coordinate Mapping(3) = 1 2 3 !The permutation of the coordinates.
Simulation Type = Steady
Steady State Max Iterations = 1
Output Intervals(1) = 1
Use Mesh Names = Logical True ! Use the mesh.names from Salome/ElmerGrid
End
Constants
Permittivity of Vacuum = 8.8542e-12
Permeability of Vacuum = $4*pi*1.0e-7
End
Body 1
Name = Air
Equation = 1
Material = 3
End
Body 2
Name = Shaft
Equation = 1
Material = 2
End
Body 3
Name = Magnets
Equation = 1
Material = 1
End
Body 4
Name = Airgap
Equation = 1
Material = 3
End
Body 5
Name = Stator
Equation = 1
Material = 2
End
Body 6
Name = CopperAndSlotAir
Equation = 1
Material = 3
End
Material 1
Relative Permeability = 1.0
Relative Permittivity = 1.0
Magnetization 1 = 1.15e6 ! magnetization [A/m]
End
Material 2
H-B Curve = Variable coupled iter
Real Monotone Cubic
0 0
1 663.146
1.1 1067.5
1.2 1705.23
1.3 2463.11
1.4 3841.67
1.5 5425.74
1.6 7957.75
1.7 12298.3
1.8 20462.8
1.9 32169.6
2 61213.4
2.1 111408
2.2 188487.757
2.3 267930.364
2.4 347507.836
End
End
Material 3
Relative Permittivity = 1.0
Relative Permeability = 1.0
Electric Conductivity = 0.0
End
Component 1
Name = string "ShaftAndMagnet"
Master Bodies(2) = integer 2 3
calculate magnetic force = logical true
Calculate Magnetic Torque = True
End
Boundary Condition 1
Name = InfAirFace
AV {e} = real 0 !Vector Potential
AV = real 0 !scalar Potential
End
Equation 1
Active Solvers(3) = 1 2 3 4
End
Solver 1
Equation = MGDynamics
Procedure = "MagnetoDynamics" "WhitneyAVSolver"
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Max Iterations = 200
Nonlinear System Convergence Tolerance = 1.0e-5
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Nonlinear System Consistent Norm = Logical True
Linear System Solver = Iterative
Linear System Iterative Method = idrs
Linear System Robust = True
Linear System Robust Start Iteration = Integer 100
Edge Basis = True
Linear System Preconditioning = ILU1
Linear System Max Iterations = 5000
Linear System Robust = True
Linear System Residual Output = 100
Linear System Convergence Tolerance = 1.0e-7
Linear System Abort Not Converged = False
Nonlinear System Consistent Norm = True
Steady State Convergence Tolerance = 1e-5
Apply Mortar BCs = True
Mortar BCs Additive = True
End
Solver 2
Equation = MGDynamicsCalc
Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields"
Nonlinear System Relaxation Factor = 1
Linear System Solver = Iterative
Linear System Symmetric = True
Linear System Iterative Method = CG
Linear System Max Iterations = 5000
Linear System Convergence Tolerance = 1.0e-8
Linear System Preconditioning = ILU0
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 1
Linear System Precondition Recompute = 1
Discontinuous Bodies = True
Calculate Elemental Fields = True
Calculate Magnetic Field Strength = True
Calculate Nodal Fields = False
Calculate Nodal Forces = False
End
Solver 3
Exec Solver = After Timestep
Equation = "ResultOutput"
Procedure = "ResultOutputSolve" "ResultOutputSolver"
Output File Name = RotorInStator
Vtu format = Logical True
Discontinuous Bodies = Logical True
File Append = Logical True
Save Geometry Ids = Logical True
End
Solver 4
Exec Solver = after all
Equation = SaveScalars
!Filename Numbering = True
Filename = ScalarValues.dat
File Append = Logical True
Procedure = "SaveData" "SaveScalars"
save component results = logical true
Variable 1 = magnetic flux density e
End
Code: Select all
ELMER SOLVER (v 9.0) STARTED AT: 2024/04/16 15:48:11
ParCommInit: Initialize #PEs: 1
MAIN:
MAIN: =============================================================
MAIN: ElmerSolver finite element software, Welcome!
MAIN: This program is free software licensed under (L)GPL
MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd.
MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi
MAIN: Version: 9.0 (Rev: unknown, Compiled: 2024-04-09)
MAIN: Running one task without MPI parallelization.
MAIN: Running with just one thread per task.
MAIN: HYPRE library linked in.
MAIN: MUMPS library linked in.
MAIN: Lua interpreter linked in.
MAIN: =============================================================
MAIN:
MAIN:
MAIN: -------------------------------------
MAIN: Reading Model: full_Motormodel.sif
LoadInputFile: Scanning input file: full_Motormodel.sif
LoadInputFile: Scanning only size info
LoadInputFile: First time visiting
LoadInputFile: Reading base load of sif file
LoadInputFile: Loading input file: full_Motormodel.sif
LoadInputFile: Reading base load of sif file
CheckKeyword: Unlisted keyword: [linear system robust start iteration] in section: [solver 1]
LoadInputFile: Number of BCs: 1
LoadInputFile: Number of Body Forces: 0
LoadInputFile: Number of Initial Conditions: 0
LoadInputFile: Number of Materials: 3
LoadInputFile: Number of Equations: 1
LoadInputFile: Number of Solvers: 4
LoadInputFile: Number of Bodies: 6
ListTagKeywords: Setting weight for keywords!
ListTagKeywords: No parameters width suffix: normalize by area
ListTagKeywords: Setting weight for keywords!
ListTagKeywords: No parameters width suffix: normalize by volume
MAIN:
MAIN: -------------------------------------
MAIN: Steady state iteration: 1
MAIN: -------------------------------------
MAIN:
WARNING:: GetPermittivity: Permittivity not defined in material, defaulting to that of vacuum
100 0.2147E+00
200 0.1152E-01
300 0.4090E-02
400 0.7343E-02
500 0.2736E-03
600 0.9992E-05
700 0.5221E-03
Idrs robust: 700 0.7967E-05 595 0.5221E-03
ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.13740310E-05 2.0000000 ) :: mgdynamics
100 0.2741E-01
200 0.5111E-02
300 0.4403E-03
400 0.3230E-04
500 0.1880E-04
Idrs robust: 555 0.1041E-04 485 0.3380E-03
ComputeChange: NS (ITER=2) (NRM,RELC): ( 0.17094581E-05 0.21756337 ) :: mgdynamics
100 0.1525E+00
200 0.1019E-02
300 0.8677E-03
400 0.1195E-03
500 0.3340E-04
600 0.1599E-04
700 0.1184E-04
800 0.1325E-04
Idrs robust: 825 0.8319E-05 760 0.3918E-03
ComputeChange: NS (ITER=3) (NRM,RELC): ( 0.13499735E-05 0.23500093 ) :: mgdynamics
100 0.1241E-01
200 0.5692E-03
300 0.2752E-03
400 0.3997E-02
500 0.8638E-04
600 0.7877E-05
700 0.1870E-03
800 0.1427E-03
Idrs robust: 840 0.2408E-05 715 0.6164E-03
ComputeChange: NS (ITER=4) (NRM,RELC): ( 0.14186588E-05 0.49616806E-01 ) :: mgdynamics
100 0.1062E-02
200 0.2042E-02
300 0.1400E-03
400 0.1077E-03
500 0.2653E-04
600 0.8078E-05
700 0.6749E-05
800 0.6966E-05
900 0.2719E-04
Idrs robust: 940 0.4828E-05 745 0.5521E-03
ComputeChange: NS (ITER=5) (NRM,RELC): ( 0.14730995E-05 0.37652292E-01 ) :: mgdynamics
100 0.3018E-02
200 0.1091E-03
300 0.9277E-04
400 0.9264E-05
500 0.1668E-04
600 0.6767E-05
700 0.4618E-05
800 0.3568E-05
Idrs robust: 845 0.3201E-05 790 0.1107E-02
ComputeChange: NS (ITER=6) (NRM,RELC): ( 0.15062592E-05 0.22259606E-01 ) :: mgdynamics
100 0.1685E-03
200 0.1825E-03
300 0.1241E-03
400 0.3412E-05
Idrs robust: 445 0.3412E-05 400 0.3311E-03
ComputeChange: NS (ITER=7) (NRM,RELC): ( 0.15271347E-05 0.13763835E-01 ) :: mgdynamics
100 0.2724E-03
200 0.3973E-04
300 0.2599E-04
400 0.3255E-05
Idrs robust: 415 0.2233E-05 390 0.5001E-03
ComputeChange: NS (ITER=8) (NRM,RELC): ( 0.15365431E-05 0.61419090E-02 ) :: mgdynamics
100 0.1612E-03
Idrs robust: 190 0.8954E-05 140 0.6061E-02
ComputeChange: NS (ITER=9) (NRM,RELC): ( 0.15396533E-05 0.20220960E-02 ) :: mgdynamics
100 0.1416E-04
200 0.4693E-04
Idrs robust: 280 0.5552E-05 250 0.7738E-03
ComputeChange: NS (ITER=10) (NRM,RELC): ( 0.15404651E-05 0.52710153E-03 ) :: mgdynamics
Idrs robust: 70 0.4097E-05 55 0.8890E-03
ComputeChange: NS (ITER=11) (NRM,RELC): ( 0.15403883E-05 0.49848167E-04 ) :: mgdynamics
Idrs robust: 50 0.4077E-05 10 0.3752E-03
ComputeChange: NS (ITER=12) (NRM,RELC): ( 0.15403859E-05 0.15304350E-05 ) :: mgdynamics
ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.15403859E-05 2.0000000 ) :: mgdynamics
WARNING:: GetPermittivity: Permittivity not defined in material, defaulting to that of vacuum
ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.0000000 0.0000000 ) :: mgdynamicscalc
MAIN: *** Elmer Solver: ALL DONE ***
MAIN: The end
SOLVER TOTAL TIME(CPU,REAL): 3623.87 2939.49
ELMER SOLVER FINISHED AT: 2024/04/16 16:37:10
Code: Select all
Check Keywords "Warn"
Header
Mesh DB "." "half_motor_mortarBC_on_merged_subfaces"
Results Directory "result_half_motor_mortarBC_on_merged_subfaces"
End
Simulation
Max Output Level = 3
Coordinate System = "Cartesian"
Coordinate Mapping(3) = 1 2 3 !The permutation of the coordinates.
Simulation Type = Steady
Steady State Max Iterations = 1
Output Intervals(1) = 1
Use Mesh Names = Logical True ! Use the mesh.names from Salome/ElmerGrid
End
Constants
Permittivity of Vacuum = 8.8542e-12
Permeability of Vacuum = $4*pi*1.0e-7
End
Body 1
Name = Air
Equation = 1
Material = 3
End
Body 2
Name = Shaft
Equation = 1
Material = 2
End
Body 3
Name = Magnet
Equation = 1
Material = 1
End
Body 4
Name = Airgap
Equation = 1
Material = 3
End
Body 5
Name = Stator
Equation = 1
Material = 2
End
Body 6
Name = CopperAndSlotAir
Equation = 1
Material = 3
End
Material 1
Relative Permeability = 1.0
Relative Permittivity = 1.0
Magnetization 1 = 1.15e6 ! magnetization [A/m]
End
Material 2
H-B Curve = Variable coupled iter
Real Monotone Cubic
0 0
1 663.146
1.1 1067.5
1.2 1705.23
1.3 2463.11
1.4 3841.67
1.5 5425.74
1.6 7957.75
1.7 12298.3
1.8 20462.8
1.9 32169.6
2 61213.4
2.1 111408
2.2 188487.757
2.3 267930.364
2.4 347507.836
End
End
Material 3
Relative Permittivity = 1.0
Relative Permeability = 1.0
Electric Conductivity = 0.0
End
Component 1
Name = string "ShaftAndMagnet"
Master Bodies(2) = integer 2 3
calculate magnetic force = logical true
Calculate Magnetic Torque = True
End
Boundary Condition 1
Name = InfAirFace
AV {e} = real 0 !Vector Potential
AV = real 0 !scalar Potential
End
Boundary Condition 2
Name = Mortarface1
Mortar BC = Integer 3
Mortar BC Static = Logical True
Anti Radial Projector = Logical True
Galerkin Projector = Logical True
End
Boundary Condition 3
Name = Mortarface2
End
Equation 1
Active Solvers(3) = 1 2 3 4
End
Solver 1
Equation = MGDynamics
Procedure = "MagnetoDynamics" "WhitneyAVSolver"
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Max Iterations = 200
Nonlinear System Convergence Tolerance = 1.0e-5
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Nonlinear System Consistent Norm = Logical True
Linear System Solver = Iterative
Linear System Iterative Method = idrs
Linear System Robust = True
Linear System Robust Start Iteration = Integer 100
Edge Basis = True
Linear System Preconditioning = ILU1
Linear System Max Iterations = 5000
Linear System Robust = True
Linear System Residual Output = 100
Linear System Convergence Tolerance = 1.0e-7
Linear System Abort Not Converged = False
Nonlinear System Consistent Norm = True
Steady State Convergence Tolerance = 1e-5
Apply Mortar BCs = True
Mortar BCs Additive = True
End
Solver 2
Equation = MGDynamicsCalc
Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields"
Nonlinear System Relaxation Factor = 1
Linear System Solver = Iterative
Linear System Symmetric = True
Linear System Iterative Method = CG
Linear System Max Iterations = 5000
Linear System Convergence Tolerance = 1.0e-8
Linear System Preconditioning = ILU0
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 1
Linear System Precondition Recompute = 1
Discontinuous Bodies = True
Calculate Elemental Fields = True
Calculate Magnetic Field Strength = True
Calculate Nodal Fields = False
Calculate Nodal Forces = False
End
Solver 3
Exec Solver = After Timestep
Equation = "ResultOutput"
Procedure = "ResultOutputSolve" "ResultOutputSolver"
Output File Name = RotorInStator
Vtu format = Logical True
Discontinuous Bodies = Logical True
File Append = Logical True
Save Geometry Ids = Logical True
End
Solver 4
Exec Solver = after all
Equation = SaveScalars
!Filename Numbering = True
Filename = ScalarValues.dat
File Append = Logical True
Procedure = "SaveData" "SaveScalars"
save component results = logical true
Variable 1 = magnetic flux density e
End
Code: Select all
ELMER SOLVER (v 9.0) STARTED AT: 2024/04/16 15:44:35
ParCommInit: Initialize #PEs: 1
MAIN:
MAIN: =============================================================
MAIN: ElmerSolver finite element software, Welcome!
MAIN: This program is free software licensed under (L)GPL
MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd.
MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi
MAIN: Version: 9.0 (Rev: unknown, Compiled: 2024-04-09)
MAIN: Running one task without MPI parallelization.
MAIN: Running with just one thread per task.
MAIN: HYPRE library linked in.
MAIN: MUMPS library linked in.
MAIN: Lua interpreter linked in.
MAIN: =============================================================
MAIN:
MAIN:
MAIN: -------------------------------------
MAIN: Reading Model: half_motor_mortarBC_on_merged_subfaces.sif
LoadInputFile: Scanning input file: half_motor_mortarBC_on_merged_subfaces.sif
LoadInputFile: Scanning only size info
LoadInputFile: First time visiting
LoadInputFile: Reading base load of sif file
LoadInputFile: Loading input file: half_motor_mortarBC_on_merged_subfaces.sif
LoadInputFile: Reading base load of sif file
CheckKeyword: Unlisted keyword: [linear system robust start iteration] in section: [solver 1]
LoadInputFile: Number of BCs: 3
LoadInputFile: Number of Body Forces: 0
LoadInputFile: Number of Initial Conditions: 0
LoadInputFile: Number of Materials: 3
LoadInputFile: Number of Equations: 1
LoadInputFile: Number of Solvers: 4
LoadInputFile: Number of Bodies: 6
ListTagKeywords: Setting weight for keywords!
ListTagKeywords: No parameters width suffix: normalize by area
ListTagKeywords: Setting weight for keywords!
ListTagKeywords: No parameters width suffix: normalize by volume
MAIN:
MAIN: -------------------------------------
MAIN: Steady state iteration: 1
MAIN: -------------------------------------
MAIN:
WARNING:: GetPermittivity: Permittivity not defined in material, defaulting to that of vacuum
100 0.4365E-01
200 0.2465E+01
300 0.5796E-02
400 0.1682E-01
500 0.1725E-01
600 0.1983E-02
700 0.1903E-02
800 0.9745E-03
900 0.8328E-03
1000 0.5124E-03
1100 0.5005E-03
1200 0.6219E-03
1300 0.1429E-02
1400 0.2821E-03
1500 0.4975E-03
1600 0.2781E-03
1700 0.7177E-02
1800 0.7393E-03
1900 0.9526E-03
2000 0.7619E-03
2100 0.5607E-03
2200 0.3356E-03
2300 0.2998E-03
2400 0.2969E-03
2500 0.3207E-03
2600 0.3832E-03
2700 0.7089E-03
2800 0.2299E-03
2900 0.7487E-03
3000 0.1076E-02
3100 0.1662E-02
3200 0.1387E-02
3300 0.1434E-02
3400 0.1547E-02
3500 0.2139E-02
3600 0.1080E-02
3700 0.4415E-03
3800 0.8697E-03
3900 0.3546E-03
4000 0.9161E-03
4100 0.5834E-01
4200 0.5327E-01
4300 0.6421E-01
4400 0.8989E-03
4500 0.1358E-02
4600 0.5277E-03
4700 0.4220E-03
4800 0.1440E-02
4900 0.3542E-02
5000 0.1555E-01
Idrs robust: 5000 0.1531E-03 2710 0.1555E-01
ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.16528967E-05 2.0000000 ) :: mgdynamics
100 0.5734E-01
200 0.4446E-02
300 0.3908E-02
400 0.3869E-02
500 0.2131E-02
600 0.4540E-03
700 0.8643E-03
800 0.4747E-01
900 0.1219E-02
1000 0.1837E-02
1100 0.2754E-03
1200 0.2166E-03
1300 0.1863E-03
1400 0.3070E-03
1500 0.8037E-03
1600 0.5116E-02
1700 0.2333E-03
1800 0.5998E-03
1900 0.2697E-03
2000 0.5218E-03
2100 0.9310E-03
2200 0.3037E-03
2300 0.1803E-02
2400 0.2327E-02
2500 0.2976E-03
2600 0.3176E-03
2700 0.4282E-03
2800 0.6379E-03
2900 0.1053E-02
3000 0.8196E-03
3100 0.1061E-02
3200 0.2389E-02
3300 0.8147E-03
3400 0.3469E-02
3500 0.4350E-03
3600 0.1462E-02
3700 0.2685E-03
3800 0.2659E-03
3900 0.5485E-03
4000 0.5433E-03
4100 0.3912E-03
4200 0.2004E-02
4300 0.3775E-03
4400 0.5700E-03
4500 0.3699E-03
4600 0.4140E-03
4700 0.2959E-03
4800 0.3398E-03
4900 0.5770E-03
5000 0.1556E-03
Idrs robust: 5000 0.1099E-03 4240 0.1556E-03
ComputeChange: NS (ITER=2) (NRM,RELC): ( 0.20482435E-05 0.21363514 ) :: mgdynamics
100 0.9026E-02
200 0.8163E-02
300 0.7671E-02
400 0.1408E-02
500 0.3596E-03
600 0.1722E-02
700 0.6939E-03
800 0.3250E-03
900 0.1789E-03
1000 0.3605E-03
1100 0.1046E-02
1200 0.1308E-02
1300 0.8795E-02
1400 0.3229E-02
1500 0.6919E-03
1600 0.2980E-03
1700 0.4301E-03
1800 0.4642E-03
1900 0.4293E-03
2000 0.3924E-03
2100 0.4509E-03
2200 0.6151E-03
2300 0.3720E-03
2400 0.9943E-03
2500 0.4824E-03
2600 0.8014E-03
2700 0.7034E-03
2800 0.1611E-02
2900 0.1246E-02
3000 0.3636E-03
3100 0.3066E-03
3200 0.3231E-03
3300 0.1116E-01
3400 0.2796E-03
3500 0.3545E-03
3600 0.4777E-03
3700 0.2875E-03
3800 0.1673E-03
3900 0.1284E-03
4000 0.1430E-03
4100 0.1342E-03
4200 0.2339E-03
4300 0.1374E-03
4400 0.4268E-03
4500 0.1451E-03
4600 0.3799E-02
4700 0.3887E-03
4800 0.1555E-03
4900 0.1929E-03
5000 0.2087E-03
Idrs robust: 5000 0.1130E-03 4830 0.2087E-03
ComputeChange: NS (ITER=3) (NRM,RELC): ( 0.19077118E-05 0.71048180E-01 ) :: mgdynamics
100 0.1161E+00
200 0.5851E-01
300 0.1694E-02
400 0.2746E-02
500 0.1860E-02
600 0.5306E-03
700 0.1849E-03
800 0.1332E-03
900 0.4041E-03
1000 0.2340E-03
1100 0.1134E-02
1200 0.1846E-03
1300 0.1243E-03
1400 0.1092E-02
1500 0.1173E-03
1600 0.1192E-02
1700 0.4368E-03
1800 0.1357E-03
1900 0.1320E-03
2000 0.6304E-04
2100 0.5634E-04
2200 0.1966E-03
2300 0.7902E-04
2400 0.8721E-04
2500 0.1553E-02
2600 0.5325E-04
2700 0.7049E-04
2800 0.5932E-04
2900 0.4700E-04
3000 0.3467E-04
3100 0.6147E-04
3200 0.4901E-04
3300 0.3994E-04
3400 0.4383E-04
3500 0.4748E-04
3600 0.3804E-04
3700 0.1202E-03
3800 0.4043E-04
3900 0.4203E-04
4000 0.6382E-04
4100 0.8822E-04
4200 0.6232E-04
4300 0.5055E-04
4400 0.5857E-04
4500 0.8166E-04
4600 0.6442E-04
4700 0.6366E-04
4800 0.4901E-04
4900 0.1324E+01
5000 0.9108E-04
Idrs robust: 5000 0.2902E-04 3570 0.9108E-04
ComputeChange: NS (ITER=4) (NRM,RELC): ( 0.19245102E-05 0.87669373E-02 ) :: mgdynamics
100 0.3149E-02
200 0.5220E-01
300 0.1087E-02
400 0.1629E-02
500 0.1934E-03
600 0.4853E-03
700 0.1312E-02
800 0.5068E-03
900 0.1603E-03
1000 0.1744E-03
1100 0.9512E-04
1200 0.8161E-04
1300 0.6535E-04
1400 0.9193E-04
1500 0.5678E-04
1600 0.3499E-04
1700 0.5281E-04
1800 0.7005E-04
1900 0.6726E-04
2000 0.3367E-03
2100 0.7675E-04
2200 0.2219E-03
2300 0.2141E-03
2400 0.6548E-03
2500 0.2368E-02
2600 0.5149E-02
2700 0.1281E-02
2800 0.2982E-03
2900 0.2602E-03
3000 0.1420E-02
3100 0.3524E-03
3200 0.3564E-03
3300 0.8199E-04
3400 0.6655E-04
3500 0.6893E-04
3600 0.1131E-03
3700 0.5020E-04
3800 0.1779E-03
3900 0.3890E-03
4000 0.3463E-03
4100 0.1104E-03
4200 0.1504E-03
4300 0.6994E-03
4400 0.1090E-03
4500 0.2590E-03
4600 0.9288E-04
4700 0.1111E-03
4800 0.1060E-03
4900 0.1197E-03
5000 0.2255E-02
Idrs robust: 5000 0.3499E-04 1600 0.2255E-02
ComputeChange: NS (ITER=5) (NRM,RELC): ( 0.19564210E-05 0.16444904E-01 ) :: mgdynamics
100 0.1090E-02
200 0.2231E-03
300 0.2219E-03
400 0.5273E-02
500 0.1396E-03
600 0.1203E-03
700 0.6056E-04
800 0.6092E-03
900 0.1981E-03
1000 0.1677E-03
1100 0.3923E-03
1200 0.1371E-03
1300 0.1324E-03
1400 0.1396E-03
1500 0.5648E-04
1600 0.1010E-03
1700 0.1195E-03
1800 0.3965E-03
1900 0.3453E-03
2000 0.1628E-03
2100 0.1017E-03
2200 0.1022E-03
2300 0.1595E-03
2400 0.5355E-04
2500 0.5135E-04
2600 0.3949E-04
2700 0.4575E-04
2800 0.4639E-04
2900 0.9762E-04
3000 0.4366E-04
3100 0.3693E-04
3200 0.4390E-04
3300 0.4569E-04
3400 0.8356E-04
3500 0.4808E-04
3600 0.4757E-04
3700 0.9120E-04
3800 0.5301E-04
3900 0.5044E-04
4000 0.4954E-04
4100 0.1494E-03
4200 0.8858E-04
4300 0.4058E-04
4400 0.3957E-04
4500 0.4587E-04
4600 0.4886E-04
4700 0.5379E-04
4800 0.4162E-04
4900 0.5408E-04
5000 0.1610E-03
Idrs robust: 5000 0.3567E-04 4910 0.1610E-03
ComputeChange: NS (ITER=6) (NRM,RELC): ( 0.20076728E-05 0.25858024E-01 ) :: mgdynamics
100 0.6953E-03
200 0.1670E-03
300 0.6896E-03
400 0.7258E-04
500 0.6017E-03
600 0.1264E-03
700 0.1723E-03
800 0.1867E-03
900 0.1801E-03
1000 0.3928E-02
1100 0.5021E-03
1200 0.1398E-03
1300 0.8174E-04
1400 0.2035E-03
1500 0.1118E-03
1600 0.1217E-03
1700 0.1594E-02
1800 0.7412E-04
1900 0.6358E-04
2000 0.4964E-04
2100 0.4421E-04
2200 0.1569E-03
2300 0.7016E-04
2400 0.1260E-03
2500 0.6605E-03
2600 0.8444E-04
2700 0.1352E-03
2800 0.1305E-03
2900 0.1016E-03
3000 0.1202E-03
3100 0.1536E-03
3200 0.9191E-04
3300 0.1256E-03
3400 0.8284E-04
3500 0.2425E-03
3600 0.7900E-04
3700 0.1509E-03
3800 0.1327E-03
3900 0.1013E-03
4000 0.2692E-03
4100 0.1654E-03
4200 0.3735E-03
4300 0.2639E-03
4400 0.1558E-03
4500 0.9077E-04
4600 0.2327E-03
4700 0.1102E-03
4800 0.6470E-04
4900 0.6087E-04
5000 0.7569E-04
Idrs robust: 5000 0.4252E-04 2525 0.7569E-04
ComputeChange: NS (ITER=7) (NRM,RELC): ( 0.20979032E-05 0.43955041E-01 ) :: mgdynamics
100 0.2261E-03
200 0.8803E-04
300 0.8937E-04
400 0.1478E-03
500 0.9063E-04
600 0.2219E-03
700 0.4796E-03
800 0.4773E-03
900 0.4985E-04
1000 0.3617E-03
1100 0.1350E-03
1200 0.3352E-04
1300 0.4178E-03
1400 0.2515E-04
1500 0.4275E-04
1600 0.3218E-04
1700 0.4539E-04
1800 0.3907E-04
1900 0.3391E-04
2000 0.7146E-04
2100 0.5988E-04
2200 0.9127E-03
2300 0.6216E-04
2400 0.6740E-04
2500 0.6125E-04
2600 0.5038E-04
2700 0.1243E-03
2800 0.5850E-04
2900 0.4636E-04
3000 0.4027E-03
3100 0.1061E-03
3200 0.1832E-03
3300 0.9466E-03
3400 0.5661E-04
3500 0.1061E-03
3600 0.4466E-04
3700 0.4662E-04
3800 0.5210E-04
3900 0.5407E-04
4000 0.5477E-04
4100 0.6025E-04
4200 0.5434E-04
4300 0.6017E-03
4400 0.4228E-04
4500 0.8160E-04
4600 0.5315E-04
4700 0.4091E-04
4800 0.1136E-03
4900 0.1408E-03
5000 0.9297E-04
Idrs robust: 5000 0.2515E-04 1400 0.9297E-04
ComputeChange: NS (ITER=8) (NRM,RELC): ( 0.21196058E-05 0.10291670E-01 ) :: mgdynamics
100 0.1740E-03
200 0.8981E-04
300 0.1461E-03
400 0.5619E-04
500 0.1334E-03
600 0.6508E-04
700 0.2088E-03
800 0.7560E-03
900 0.4904E-03
1000 0.1542E-03
1100 0.3296E-03
1200 0.2752E-03
1300 0.4075E-03
1400 0.1528E-03
1500 0.5445E-03
1600 0.5775E-02
1700 0.4786E-03
1800 0.6742E-03
1900 0.1583E-03
2000 0.1533E-03
2100 0.9232E-03
2200 0.6532E-03
2300 0.3715E-03
2400 0.3318E-03
2500 0.1797E-02
2600 0.2619E-03
2700 0.2201E-03
2800 0.2370E-03
2900 0.3484E-03
3000 0.4368E-03
3100 0.2436E-03
3200 0.3048E-03
3300 0.4821E-03
3400 0.4444E-03
3500 0.2834E-03
3600 0.4612E-01
3700 0.4678E-01
3800 0.6908E-01
3900 0.8000E-03
4000 0.1619E-03
4100 0.1311E-03
4200 0.1245E-03
4300 0.1376E-03
4400 0.1229E-03
4500 0.3482E-03
4600 0.4219E-03
4700 0.3832E-03
4800 0.3750E-02
4900 0.2498E-01
5000 0.1764E-01
Idrs robust: 5000 0.3170E-04 665 0.1764E-01
ComputeChange: NS (ITER=9) (NRM,RELC): ( 0.21283902E-05 0.41357627E-02 ) :: mgdynamics
100 0.1412E-03
200 0.4192E-04
300 0.3013E-04
400 0.4182E-04
500 0.4341E-04
600 0.6403E-04
700 0.1194E-03
800 0.2392E-03
900 0.8726E-04
1000 0.1836E-03
1100 0.5637E-03
1200 0.2849E-03
1300 0.1471E-02
1400 0.8008E-03
1500 0.1722E-03
1600 0.2360E-03
1700 0.9088E-04
1800 0.9890E-04
1900 0.1949E-03
2000 0.6464E-03
2100 0.8430E-04
2200 0.8259E-04
2300 0.5120E-04
2400 0.7815E-04
2500 0.5250E-04
2600 0.1003E-03
2700 0.5074E-03
2800 0.6467E-04
2900 0.5840E-04
3000 0.3776E-02
3100 0.1437E-03
3200 0.9350E-04
3300 0.4252E-04
3400 0.2693E-03
3500 0.9938E-04
3600 0.8725E-04
3700 0.2901E-03
3800 0.9337E-04
3900 0.3574E-04
4000 0.6219E-04
4100 0.3797E-04
4200 0.3057E-03
4300 0.6535E-04
4400 0.5989E-04
4500 0.4812E-04
4600 0.4745E-04
4700 0.8593E-04
4800 0.5995E-04
4900 0.5802E-04
5000 0.1242E-03
Idrs robust: 5000 0.3417E-04 4575 0.1242E-03
ComputeChange: NS (ITER=10) (NRM,RELC): ( 0.21868417E-05 0.27090796E-01 ) :: mgdynamics
100 0.1744E-03
200 0.1459E-03
300 0.2095E-03
400 0.1006E-02
500 0.5140E-03
600 0.3898E-03
700 0.2323E-03
800 0.2059E-03
900 0.3200E-04
1000 0.5168E-04
1100 0.6973E-04
1200 0.2493E-03
1300 0.4417E-03
1400 0.2829E-03
1500 0.1208E-03
1600 0.7260E-04
1700 0.6018E-04
1800 0.7982E-04
1900 0.1596E-02
2000 0.4483E-04
2100 0.3578E-04
2200 0.4309E-04
2300 0.4472E-04
2400 0.7075E-04
2500 0.1329E-03
2600 0.1340E-02
2700 0.7290E-04
2800 0.6143E-04
2900 0.1513E-03
3000 0.9294E-04
3100 0.8880E-04
3200 0.7480E-04
3300 0.7704E-04
3400 0.7860E-04
3500 0.9248E-04
3600 0.8178E-04
3700 0.7051E-03
3800 0.1134E-03
3900 0.1391E-03
4000 0.4750E-03
4100 0.5589E-03
4200 0.1974E-03
4300 0.1125E-02
4400 0.5846E-03
4500 0.6031E-02
4600 0.4903E-02
4700 0.4657E-01
4800 0.2294E-01
4900 0.1062E+00
5000 0.4345E-01
Idrs robust: 5000 0.4143E-04 3020 0.4345E-01
ComputeChange: NS (ITER=11) (NRM,RELC): ( 0.23505468E-05 0.72158267E-01 ) :: mgdynamics
100 0.5705E-04
200 0.8221E-04
300 0.1581E-03
400 0.9373E-04
500 0.1214E-03
600 0.1067E-03
700 0.6260E-04
800 0.3819E-03
900 0.7946E-04
1000 0.1887E-03
1100 0.1375E-02
1200 0.6685E-04
1300 0.6561E-04
1400 0.3647E-03
1500 0.1714E-02
1600 0.4273E-03
1700 0.6005E-03
1800 0.4780E-03
1900 0.1162E-02
2000 0.8036E-03
2100 0.1226E-02
2200 0.2284E-02
2300 0.2764E-02
2400 0.4127E-03
2500 0.1884E-02
2600 0.7938E-03
2700 0.3305E-03
2800 0.1037E-02
2900 0.5034E-03
3000 0.6103E-03
3100 0.2322E-02
3200 0.1745E-03
3300 0.2865E-03
3400 0.2082E-03
3500 0.1943E-03
3600 0.3372E-03
3700 0.4290E-03
3800 0.4761E-03
3900 0.3324E-03
4000 0.2160E-03
4100 0.3798E-03
4200 0.1285E-03
4300 0.1127E-03
4400 0.1314E-03
4500 0.1377E-03
4600 0.4384E-02
4700 0.4266E-03
4800 0.1059E-02
4900 0.1594E-02
5000 0.2157E-02
Idrs robust: 5000 0.3510E-04 910 0.2157E-02
ComputeChange: NS (ITER=12) (NRM,RELC): ( 0.23569727E-05 0.27300669E-02 ) :: mgdynamics
100 0.8482E-04
200 0.7117E-04
300 0.1292E-03
400 0.3226E-03
500 0.1675E-03
600 0.6269E-03
700 0.1033E-03
800 0.1411E-02
900 0.5997E-03
1000 0.1594E-01
1100 0.3991E-03
1200 0.9130E-03
1300 0.2046E-03
1400 0.1401E-03
1500 0.1815E-03
1600 0.4294E-02
1700 0.9925E-04
1800 0.6496E-04
1900 0.6199E-04
2000 0.5417E-04
2100 0.8082E-04
2200 0.1632E-03
2300 0.9667E-04
2400 0.2035E-03
2500 0.6187E-04
2600 0.1895E-03
2700 0.5746E-03
2800 0.2939E-03
2900 0.2213E-03
3000 0.3280E-03
3100 0.1132E-03
3200 0.3875E-03
3300 0.2985E-03
3400 0.2875E-03
3500 0.2629E-03
3600 0.1210E-03
3700 0.8552E-03
3800 0.1394E-03
3900 0.7870E-04
4000 0.1786E-03
4100 0.6603E-04
4200 0.6633E-04
4300 0.7103E-04
4400 0.5925E-04
4500 0.5526E-04
4600 0.4614E-04
4700 0.3807E-03
4800 0.5208E-02
4900 0.8112E-03
5000 0.3894E-04
Idrs robust: 5000 0.2912E-04 4485 0.3894E-04
ComputeChange: NS (ITER=13) (NRM,RELC): ( 0.25422490E-05 0.75635004E-01 ) :: mgdynamics
100 0.3873E-04
200 0.6708E-04
300 0.4783E-04
400 0.3595E-04
500 0.4847E-04
600 0.3326E-04
700 0.2592E-04
800 0.4575E-04
900 0.5361E-04
1000 0.5779E-04
1100 0.5509E-04
1200 0.9493E-04
1300 0.6002E-04
1400 0.5358E-03
1500 0.2459E-03
1600 0.4190E-03
1700 0.3570E-03
1800 0.3226E-03
1900 0.1296E-03
2000 0.1371E-03
2100 0.4951E-03
2200 0.9278E-04
2300 0.9028E-04
2400 0.7545E-04
2500 0.6475E-04
2600 0.6091E-04
2700 0.1181E-03
2800 0.4974E-03
2900 0.1172E-03
3000 0.1126E-02
3100 0.5904E-04
3200 0.9492E-04
3300 0.9257E-04
3400 0.1116E-03
3500 0.2683E-03
3600 0.4844E-03
3700 0.3717E-02
3800 0.1104E-02
3900 0.1019E-03
4000 0.1016E-03
4100 0.1149E-03
4200 0.1628E-03
4300 0.1131E-03
4400 0.5332E-04
4500 0.1420E-03
4600 0.2113E-03
4700 0.8746E-04
4800 0.5831E-03
4900 0.3601E-04
5000 0.6768E-04
Idrs robust: 5000 0.2447E-04 1165 0.6768E-04
ComputeChange: NS (ITER=14) (NRM,RELC): ( 0.25518963E-05 0.37876125E-02 ) :: mgdynamics
100 0.3516E-03
200 0.2378E-04
Idrs robust: 225 0.1497E-04 220 0.3840E-03
ComputeChange: NS (ITER=15) (NRM,RELC): ( 0.25520763E-05 0.70512554E-04 ) :: mgdynamics
100 0.1438E-04
200 0.2702E-04
Idrs robust: 275 0.1438E-04 100 0.2050E-02
ComputeChange: NS (ITER=16) (NRM,RELC): ( 0.25520805E-05 0.16434050E-05 ) :: mgdynamics
ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.25520805E-05 2.0000000 ) :: mgdynamics
WARNING:: GetPermittivity: Permittivity not defined in material, defaulting to that of vacuum
ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.0000000 0.0000000 ) :: mgdynamicscalc
MAIN: *** Elmer Solver: ALL DONE ***
MAIN: The end
SOLVER TOTAL TIME(CPU,REAL): 20872.46 8278.88
ELMER SOLVER FINISHED AT: 2024/04/16 18:02:34
half_motor_mortarBC_on_all_subfaces.sif
Code: Select all
Check Keywords "Warn"
Header
Mesh DB "." "half_motor_mortarBC_on_all_subfaces"
Results Directory "result_half_motor_mortarBC_on_all_subfaces"
End
Simulation
Max Output Level = 3
Coordinate System = "Cartesian"
Coordinate Mapping(3) = 1 2 3 !The permutation of the coordinates.
Simulation Type = Steady
Steady State Max Iterations = 1
Output Intervals(1) = 1
Use Mesh Names = Logical True ! Use the mesh.names from Salome/ElmerGrid
End
Constants
Permittivity of Vacuum = 8.8542e-12
Permeability of Vacuum = $4*pi*1.0e-7
End
Body 1
Name = Air
Equation = 1
Material = 3
End
Body 2
Name = Shaft
Equation = 1
Material = 2
End
Body 3
Name = Magnet
Equation = 1
Material = 1
End
Body 4
Name = Airgap
Equation = 1
Material = 3
End
Body 5
Name = Stator
Equation = 1
Material = 2
End
Body 6
Name = CopperAndSlotAir
Equation = 1
Material = 3
End
Material 1
Relative Permeability = 1.0
Relative Permittivity = 1.0
Magnetization 1 = 1.15e6 ! magnetization [A/m]
End
Material 2
H-B Curve = Variable coupled iter
Real Monotone Cubic
0 0
1 663.146
1.1 1067.5
1.2 1705.23
1.3 2463.11
1.4 3841.67
1.5 5425.74
1.6 7957.75
1.7 12298.3
1.8 20462.8
1.9 32169.6
2 61213.4
2.1 111408
2.2 188487.757
2.3 267930.364
2.4 347507.836
End
End
Material 3
Relative Permittivity = 1.0
Relative Permeability = 1.0
Electric Conductivity = 0.0
End
Component 1
Name = string "ShaftAndMagnet"
Master Bodies(2) = integer 2 3
calculate magnetic force = logical true
Calculate Magnetic Torque = True
End
Boundary Condition 1
Name = InfAirFace
AV {e} = real 0 !Vector Potential
AV = real 0 !scalar Potential
End
Boundary Condition 2
Name = Mortar1_Air
Mortar BC = Integer 3
Mortar BC Static = Logical True
Anti Radial Projector = Logical True
Galerkin Projector = Logical True
End
Boundary Condition 3
Name = Mortar2_Air
End
Boundary Condition 4
Name = Mortar1_Stator
Mortar BC = Integer 5
Mortar BC Static = Logical True
Anti Radial Projector = Logical True
Galerkin Projector = Logical True
End
Boundary Condition 5
Name = Mortar2_Stator
End
Boundary Condition 6
Name = Mortar1_Copper
Mortar BC = 7
Mortar BC Static = Logical True
Anti Radial Projector = Logical True
Galerkin Projector = Logical True
End
Boundary Condition 7
Name = Mortar2_Copper
End
Boundary Condition 8
Name = Mortar1_SlotAir
Mortar BC = 9
Mortar BC Static = Logical True
Anti Radial Projector = Logical True
Galerkin Projector = Logical True
End
Boundary Condition 9
Name = Mortar2_SlotAir
End
Boundary Condition 10
Name = Mortar1_AirgapStator
Mortar BC = 11
Mortar BC Static = Logical True
Anti Radial Projector = Logical True
Galerkin Projector = Logical True
End
Boundary Condition 11
Name = Mortar2_AirgapStator
End
Boundary Condition 12
Name = Mortar1_AirgapMiddle
Mortar BC = 13
Mortar BC Static = Logical True
Anti Radial Projector = Logical True
Galerkin Projector = Logical True
End
Boundary Condition 13
Name = Mortar2_AirgapMiddle
End
Boundary Condition 14
Name = Mortar1_AirgapRotor
Mortar BC = 15
Mortar BC Static = Logical True
Anti Radial Projector = Logical True
Galerkin Projector = Logical True
End
Boundary Condition 15
Name = Mortar2_AirgapRotor
End
Boundary Condition 16
Name = Mortar1_Magnet
Mortar BC = 17
Mortar BC Static = Logical True
Anti Radial Projector = Logical True
Galerkin Projector = Logical True
End
Boundary Condition 17
Name = Mortar2_Magnet
End
Boundary Condition 18
Name = Mortar1_Shaft
Mortar BC = 19
Mortar BC Static = Logical True
Anti Radial Projector = Logical True
Galerkin Projector = Logical True
End
Boundary Condition 19
Name = Mortar2_Shaft
End
Equation 1
Active Solvers(3) = 1 2 3 4
End
Solver 1
Equation = MGDynamics
Procedure = "MagnetoDynamics" "WhitneyAVSolver"
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Max Iterations = 200
Nonlinear System Convergence Tolerance = 1.0e-5
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Nonlinear System Consistent Norm = Logical True
Linear System Solver = Iterative
Linear System Iterative Method = idrs
Linear System Robust = True
Linear System Robust Start Iteration = Integer 100
Edge Basis = True
Linear System Preconditioning = ILU1
Linear System Max Iterations = 5000
Linear System Robust = True
Linear System Residual Output = 100
Linear System Convergence Tolerance = 1.0e-7
Linear System Abort Not Converged = False
Nonlinear System Consistent Norm = True
Steady State Convergence Tolerance = 1e-5
Apply Mortar BCs = True
Mortar BCs Additive = True
End
Solver 2
Equation = MGDynamicsCalc
Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields"
Nonlinear System Relaxation Factor = 1
Linear System Solver = Iterative
Linear System Symmetric = True
Linear System Iterative Method = CG
Linear System Max Iterations = 5000
Linear System Convergence Tolerance = 1.0e-8
Linear System Preconditioning = ILU0
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 1
Linear System Precondition Recompute = 1
Discontinuous Bodies = True
Calculate Elemental Fields = True
Calculate Magnetic Field Strength = True
Calculate Nodal Fields = False
Calculate Nodal Forces = False
End
Solver 3
Exec Solver = After Timestep
Equation = "ResultOutput"
Procedure = "ResultOutputSolve" "ResultOutputSolver"
Output File Name = RotorInStator
Vtu format = Logical True
Discontinuous Bodies = Logical True
File Append = Logical True
Save Geometry Ids = Logical True
End
Solver 4
Exec Solver = after all
Equation = SaveScalars
!Filename Numbering = True
Filename = ScalarValues.dat
File Append = Logical True
Procedure = "SaveData" "SaveScalars"
save component results = logical true
Variable 1 = magnetic flux density e
End
Code: Select all
ELMER SOLVER (v 9.0) STARTED AT: 2024/04/16 15:43:14
ParCommInit: Initialize #PEs: 1
MAIN:
MAIN: =============================================================
MAIN: ElmerSolver finite element software, Welcome!
MAIN: This program is free software licensed under (L)GPL
MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd.
MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi
MAIN: Version: 9.0 (Rev: unknown, Compiled: 2024-04-09)
MAIN: Running one task without MPI parallelization.
MAIN: Running with just one thread per task.
MAIN: HYPRE library linked in.
MAIN: MUMPS library linked in.
MAIN: Lua interpreter linked in.
MAIN: =============================================================
MAIN:
MAIN:
MAIN: -------------------------------------
MAIN: Reading Model: half_motor_mortarBC_on_all_subfaces.sif
LoadInputFile: Scanning input file: half_motor_mortarBC_on_all_subfaces.sif
LoadInputFile: Scanning only size info
LoadInputFile: First time visiting
LoadInputFile: Reading base load of sif file
LoadInputFile: Loading input file: half_motor_mortarBC_on_all_subfaces.sif
LoadInputFile: Reading base load of sif file
CheckKeyword: Unlisted keyword: [linear system robust start iteration] in section: [solver 1]
LoadInputFile: Number of BCs: 19
LoadInputFile: Number of Body Forces: 0
LoadInputFile: Number of Initial Conditions: 0
LoadInputFile: Number of Materials: 3
LoadInputFile: Number of Equations: 1
LoadInputFile: Number of Solvers: 4
LoadInputFile: Number of Bodies: 6
ListTagKeywords: Setting weight for keywords!
ListTagKeywords: No parameters width suffix: normalize by area
ListTagKeywords: Setting weight for keywords!
ListTagKeywords: No parameters width suffix: normalize by volume
MAIN:
MAIN: -------------------------------------
MAIN: Steady state iteration: 1
MAIN: -------------------------------------
MAIN:
WARNING:: GetPermittivity: Permittivity not defined in material, defaulting to that of vacuum
100 0.3880E-01
200 0.2924E-01
300 0.9412E-02
400 0.6834E-03
500 0.1168E-02
600 0.8027E-01
700 0.2831E-02
800 0.1754E-02
900 0.1246E-02
1000 0.4536E-02
1100 0.2538E-02
1200 0.8762E-03
1300 0.1309E-02
1400 0.1663E-02
1500 0.3396E-02
1600 0.3500E-02
1700 0.1828E-02
1800 0.1454E-02
1900 0.3718E-03
2000 0.2528E-03
2100 0.2672E-03
2200 0.7713E-03
2300 0.4089E-03
2400 0.1160E-03
2500 0.1453E-03
2600 0.1310E-03
2700 0.1527E-03
2800 0.1123E-03
2900 0.2380E-03
3000 0.1708E-03
3100 0.1813E-03
3200 0.2496E-03
3300 0.2815E-03
3400 0.9011E-03
3500 0.1609E-02
3600 0.4387E-03
3700 0.2509E-03
3800 0.2416E-02
3900 0.1861E-03
4000 0.3485E-03
4100 0.2443E-03
4200 0.4948E-03
4300 0.5665E-03
4400 0.5461E-03
4500 0.3993E-03
4600 0.6255E-03
4700 0.2094E-03
4800 0.2185E-03
4900 0.2140E-03
5000 0.1863E-01
Idrs robust: 5000 0.1387E-03 4920 0.1863E-01
ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.16491115E-05 2.0000000 ) :: mgdynamics
100 0.5221E-01
200 0.7850E-02
300 0.3552E-02
400 0.1283E-02
500 0.1173E-02
600 0.3196E-03
700 0.4041E-03
800 0.3319E-03
900 0.3606E-03
1000 0.1074E-02
1100 0.1261E-02
1200 0.3498E-03
1300 0.4866E-03
1400 0.9785E-04
1500 0.7129E-04
1600 0.1058E-03
1700 0.8752E-04
1800 0.8615E-04
1900 0.9749E-04
2000 0.8119E-04
2100 0.8297E-04
2200 0.1105E-03
2300 0.3337E-03
2400 0.1532E-02
2500 0.1479E-03
2600 0.1348E-03
2700 0.1190E-03
2800 0.1008E-01
2900 0.2170E-03
3000 0.2197E-03
3100 0.2770E-03
3200 0.1123E-03
3300 0.8318E-04
3400 0.1131E-03
3500 0.1125E-03
3600 0.2121E-03
3700 0.3657E-03
3800 0.4940E-03
3900 0.3823E-03
4000 0.1540E-03
4100 0.1645E-03
4200 0.1222E-03
4300 0.2954E-03
4400 0.4330E-03
4500 0.1575E-03
4600 0.1283E-03
4700 0.1721E-03
4800 0.5378E-03
4900 0.2385E-02
5000 0.1536E-03
Idrs robust: 5000 0.7444E-04 4865 0.1536E-03
ComputeChange: NS (ITER=2) (NRM,RELC): ( 0.19848751E-05 0.18479079 ) :: mgdynamics
100 0.9706E-02
200 0.3096E-02
300 0.2706E-02
400 0.2288E-02
500 0.3641E-02
600 0.3410E-01
700 0.1246E-02
800 0.1068E-02
900 0.1441E-03
1000 0.1623E-03
1100 0.1005E-03
1200 0.1621E-03
1300 0.1124E-02
1400 0.2210E-03
1500 0.1120E-02
1600 0.8404E-03
1700 0.2061E-03
1800 0.2857E-03
1900 0.2330E-03
2000 0.2200E-02
2100 0.1364E-03
2200 0.1099E-03
2300 0.6600E-04
2400 0.5640E-04
2500 0.1736E-03
2600 0.6518E-04
2700 0.6776E-04
2800 0.2171E-03
2900 0.1219E-03
3000 0.8545E-04
3100 0.8031E-04
3200 0.6520E-04
3300 0.2064E-03
3400 0.6641E-04
3500 0.1579E-03
3600 0.5776E-04
3700 0.8523E-04
3800 0.7015E-04
3900 0.9992E-04
4000 0.1540E-03
4100 0.8005E-04
4200 0.1017E-03
4300 0.1129E-03
4400 0.1020E-03
4500 0.1823E-03
4600 0.2349E-03
4700 0.4124E-03
4800 0.9584E-03
4900 0.1256E-02
5000 0.1064E-03
Idrs robust: 5000 0.6490E-04 4795 0.1064E-03
ComputeChange: NS (ITER=3) (NRM,RELC): ( 0.17651062E-05 0.11721071 ) :: mgdynamics
100 0.4294E-01
200 0.5002E-02
300 0.9649E-02
400 0.2623E-02
500 0.1013E-02
600 0.4991E-03
700 0.1171E-03
800 0.7714E-04
900 0.7476E-04
1000 0.1840E-03
1100 0.7671E-03
1200 0.1859E-03
1300 0.2861E-03
1400 0.8005E-04
1500 0.1170E-03
1600 0.1009E-03
1700 0.8476E-04
1800 0.1060E-03
1900 0.6088E-04
2000 0.9219E-04
2100 0.1648E-03
2200 0.1131E-03
2300 0.3060E-03
2400 0.6591E-03
2500 0.1995E-03
2600 0.8818E-04
2700 0.1643E-03
2800 0.1502E-03
2900 0.1242E-03
3000 0.8042E-04
3100 0.7046E-03
3200 0.2789E-03
3300 0.2133E-03
3400 0.5182E-03
3500 0.1248E-02
3600 0.4604E-03
3700 0.4232E-03
3800 0.1388E-02
3900 0.1796E-02
4000 0.1270E-02
4100 0.1461E-02
4200 0.5225E-03
4300 0.1720E-02
4400 0.1211E-01
4500 0.1505E-03
4600 0.1560E-03
4700 0.6452E-03
4800 0.1071E-03
4900 0.7131E-04
5000 0.2687E-03
Idrs robust: 5000 0.5588E-04 4845 0.2687E-03
ComputeChange: NS (ITER=4) (NRM,RELC): ( 0.18575934E-05 0.51059852E-01 ) :: mgdynamics
100 0.1019E-01
200 0.2745E-03
300 0.2460E-03
400 0.2632E-03
500 0.1951E-03
600 0.5136E-04
700 0.4994E-04
800 0.5374E-04
900 0.2837E-02
1000 0.5507E-04
1100 0.5124E-04
1200 0.5309E-04
1300 0.1204E-02
1400 0.7403E-04
1500 0.2041E-03
1600 0.7761E-03
1700 0.1348E-03
1800 0.1417E-03
1900 0.2398E-03
2000 0.4738E-03
2100 0.1576E-03
2200 0.2213E-02
2300 0.1739E-03
2400 0.1403E-02
2500 0.5035E-03
2600 0.4062E-03
2700 0.3774E-03
2800 0.6790E-03
2900 0.1985E-02
3000 0.2170E-03
3100 0.3747E-03
3200 0.3694E-03
3300 0.1565E-03
3400 0.3292E-03
3500 0.5279E-03
3600 0.4450E-03
3700 0.1620E-02
3800 0.9739E-04
3900 0.1549E-02
4000 0.1134E-03
4100 0.6007E-03
4200 0.6418E-03
4300 0.4444E-03
4400 0.4696E-03
4500 0.1372E-03
4600 0.3901E-03
4700 0.7372E-03
4800 0.1463E-03
4900 0.1225E-03
5000 0.2867E-02
Idrs robust: 5000 0.3223E-04 810 0.2867E-02
ComputeChange: NS (ITER=5) (NRM,RELC): ( 0.19026571E-05 0.23968464E-01 ) :: mgdynamics
100 0.8438E-02
200 0.2140E-03
300 0.1888E-03
400 0.4525E-04
500 0.3663E-04
600 0.2653E-04
Idrs robust: 650 0.2120E-04 630 0.1630E-02
ComputeChange: NS (ITER=6) (NRM,RELC): ( 0.19261776E-05 0.12285997E-01 ) :: mgdynamics
100 0.5418E-03
200 0.2701E-03
300 0.1884E-03
400 0.5524E-04
500 0.4772E-04
600 0.5359E-03
700 0.5677E-04
800 0.5470E-04
900 0.6209E-04
1000 0.5091E-04
1100 0.6698E-04
1200 0.1060E-03
1300 0.1057E-03
1400 0.1152E-03
1500 0.5611E-04
1600 0.1348E-03
1700 0.9440E-04
1800 0.1943E-03
1900 0.8569E-04
2000 0.5772E-04
2100 0.1816E-03
2200 0.4561E-04
2300 0.4205E-04
2400 0.1466E-03
2500 0.4020E-04
2600 0.4485E-04
2700 0.7603E-04
2800 0.1979E-03
2900 0.7292E-04
3000 0.9839E-04
3100 0.8471E-04
3200 0.1220E-03
3300 0.9583E-04
3400 0.9562E-04
3500 0.1277E-03
3600 0.8631E-04
3700 0.8133E-04
3800 0.6820E-04
3900 0.1031E-03
4000 0.1015E-03
4100 0.1384E-03
4200 0.1238E-03
4300 0.4883E-04
4400 0.4781E-04
4500 0.8367E-04
4600 0.1592E-03
4700 0.6709E-04
4800 0.3101E-04
4900 0.3991E-04
5000 0.1409E-03
Idrs robust: 5000 0.3199E-04 4915 0.1409E-03
ComputeChange: NS (ITER=7) (NRM,RELC): ( 0.20906776E-05 0.81904850E-01 ) :: mgdynamics
100 0.4158E-03
200 0.8309E-04
300 0.9866E-04
400 0.1479E-03
500 0.3502E-04
600 0.8174E-04
700 0.9188E-04
800 0.3785E-03
Idrs robust: 880 0.1505E-04 870 0.1005E-02
ComputeChange: NS (ITER=8) (NRM,RELC): ( 0.20975646E-05 0.32887174E-02 ) :: mgdynamics
100 0.7678E-03
200 0.5538E-04
300 0.1616E-04
Idrs robust: 375 0.1394E-04 355 0.4013E-03
ComputeChange: NS (ITER=9) (NRM,RELC): ( 0.20999484E-05 0.11358068E-02 ) :: mgdynamics
100 0.5252E-03
200 0.3394E-04
300 0.1518E-03
Idrs robust: 305 0.1577E-04 270 0.5544E-03
ComputeChange: NS (ITER=10) (NRM,RELC): ( 0.21000401E-05 0.43699473E-04 ) :: mgdynamics
100 0.3583E-04
200 0.5820E-04
300 0.2339E-04
Idrs robust: 380 0.1978E-04 360 0.3581E-03
ComputeChange: NS (ITER=11) (NRM,RELC): ( 0.21008872E-05 0.40326485E-03 ) :: mgdynamics
100 0.6272E-04
Idrs robust: 170 0.1799E-04 165 0.6051E-03
ComputeChange: NS (ITER=12) (NRM,RELC): ( 0.21010539E-05 0.79334162E-04 ) :: mgdynamics
100 0.1233E-04
200 0.1187E-04
300 0.1564E-04
Idrs robust: 320 0.1083E-04 260 0.8151E-03
ComputeChange: NS (ITER=13) (NRM,RELC): ( 0.21008909E-05 0.77556466E-04 ) :: mgdynamics
100 0.1323E-03
200 0.1919E-04
300 0.1695E-04
Idrs robust: 370 0.1434E-04 180 0.5834E-03
ComputeChange: NS (ITER=14) (NRM,RELC): ( 0.21008100E-05 0.38529545E-04 ) :: mgdynamics
100 0.7409E-04
200 0.4570E-04
Idrs robust: 250 0.1273E-04 160 0.4998E-03
ComputeChange: NS (ITER=15) (NRM,RELC): ( 0.21008388E-05 0.13715179E-04 ) :: mgdynamics
100 0.1138E-04
200 0.1827E-04
Idrs robust: 225 0.1161E-04 105 0.4032E-03
ComputeChange: NS (ITER=16) (NRM,RELC): ( 0.21008561E-05 0.82313607E-05 ) :: mgdynamics
ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.21008561E-05 2.0000000 ) :: mgdynamics
WARNING:: GetPermittivity: Permittivity not defined in material, defaulting to that of vacuum
ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.0000000 0.0000000 ) :: mgdynamicscalc
MAIN: *** Elmer Solver: ALL DONE ***
MAIN: The end
SOLVER TOTAL TIME(CPU,REAL): 9973.54 4987.58
ELMER SOLVER FINISHED AT: 2024/04/16 17:06:21
3. What could be the reason for the different field distribution?
You can download the simulation models in the next post.
Re: Simulation of Magnetic Force with a half 3D Motormodel
I would be very thankful if somebody could have a look at my case and answer my three questions in the post above.
- Attachments
-
- half_motor_mortarBC_on_all_subfaces.zip
- (12.9 KiB) Downloaded 8 times
-
- half_motor_mortarBC_on_merged_subfaces.zip
- (12.2 KiB) Downloaded 7 times
-
- full_Motormodel.zip
- (10 KiB) Downloaded 12 times