Coupling CircuitsAndDynamics - MagnetoDynamics2D (Eddy Currents On)

Numerical methods and mathematical models of Elmer
Post Reply
ajpina
Posts: 7
Joined: 22 Mar 2018, 01:15
Antispam: Yes

Coupling CircuitsAndDynamics - MagnetoDynamics2D (Eddy Currents On)

Post by ajpina »

Hi Forum,
I am trying to do an analysis on skin and proximity effects on stranded conductors by using the coupling between CircuitsAndDynamics and MagnetoDynamics2D Solvers. So I wonder if there is anyway I can keep the eddy currents enabled while using the Component Definition (Coil Type = String Stranded).

Here are part of my settings:

Code: Select all

Body 11
	Name = A1T1
	Equation = 1
	Material = 5
End

Body 12
	Name = A1T2
	Equation = 1
	Material = 5
End

Body 13
	Name = A1T3
	Equation = 1
	Material = 5
End

Body 14
	Name = A1T4
	Equation = 1
	Material = 5
End
Component 1
	Name = String A1
	Master Bodies(4) = Integer 11 12 13 14
	Coil Type = String Stranded
	Number of Turns = Real 4
End
Each body is a bundle of wires as can be seen in the attached picture, and the component defines together 4 bundles. These pattern is repeated several times to account for all coils and turns. The circuit definition is shown below:

Code: Select all

$ Circuits = 3
$ noVariables = 10

!!! Phase A
$ C.1.variables = noVariables					

! ------------------------ Circuit 1 ------------------------
$ C.1.perm = zeros(noVariables)
$ C.1.A = zeros(noVariables,noVariables)	
$ C.1.B = zeros(noVariables,noVariables)
$ C.1.Mre = zeros(noVariables,noVariables)	
$ C.1.Mim = zeros(noVariables,noVariables)	

! Define variables.
$ C.1.name.1 = "i_a"
$ C.1.name.2 = "v_a"
$ C.1.name.3 = "i_component(1)"
$ C.1.name.4 = "v_component(1)"
$ C.1.name.5 = "i_component(2)"
$ C.1.name.6 = "v_component(2)"
$ C.1.name.7 = "i_component(3)"
$ C.1.name.8 = "v_component(3)"
$ C.1.name.9 = "i_component(4)"
$ C.1.name.10 = "v_component(4)"

! Define sources:
!------------------
$ C.1.B(0,0) = 1					
$ C.1.source.1 = "ia"		

! Define network cycles:
!-------------------------
!Voltage relations.
$ C.1.B(1,1) = -1
$ C.1.B(1,3) = 1
$ C.1.B(1,5) = 1
$ C.1.B(1,7) = 1
$ C.1.B(1,9) = 1

!Current relations.
$ C.1.B(2,0) = 1				! ia = i1
$ C.1.B(2,2) = -1
$ C.1.B(3,0) = 1				! ia = i2
$ C.1.B(3,4) = -1
$ C.1.B(4,2) = -1				! i1 = i2
$ C.1.B(4,4) = 1
$ C.1.B(5,0) = 1				! ia = i3
$ C.1.B(5,6) = -1
$ C.1.B(6,4) = -1				! i2 = i3
$ C.1.B(6,6) = 1
$ C.1.B(7,0) = 1				! ia = i4
$ C.1.B(7,8) = -1
$ C.1.B(8,6) = -1				! i3 = i4
$ C.1.B(8,8) = 1
The top figure in the picture attached is the kind of results I am after and the bottom one are the results I am getting with my current settings, so I wonder if anyone may point out the changes I need to do in order to solve this analysis.
StrandedCond.PNG
StrandedCond.PNG (139.85 KiB) Viewed 4164 times
My solver settings are:

Code: Select all

Solver 2
	Exec Solver = Always
	Equation = Circuits
	Variable = X
	Procedure = "CircuitsAndDynamics" "CircuitsAndDynamics"
	No Matrix = Logical True
End

Solver 3
	Equation = MgDyn2D
	Procedure = "MagnetoDynamics2D" "MagnetoDynamics2D"
	Exec Solver = Always
	Variable = A
	Nonlinear System Convergence Tolerance = 1e-06
	Nonlinear System Max Iterations = 100
	Nonlinear System Min Iterations = 1
	Nonlinear System Newton After Iterations = 5
	Nonlinear System Relaxation Factor = 0.9
	Nonlinear System Convergence Without Constraints = Logical True
	Export Lagrange Multiplier = Logical True
	Linear System Abort Not Converged = Logical False
	Linear System Solver = Direct
	Linear System Direct Method = umfpack
	Optimize Bandwidth = Logical True
	Linear System Preconditioning =  ILU2
	Linear System Max Iterations =  5000
	Linear System Residual Output =  20
	Linear System Convergence Tolerance =  1e-07
	Mortar BCs Additive =  Logical True
End
Post Reply