functions.magnetic.circuits: Magnetic circuit calculations
In renpow: Renewable Power Systems and the Environment

Description Usage Arguments Details Value Note Author(s) References Examples

Calculates magnetic elements and circuits

1
2
3

reluctance(x)
inductor(x)
flux(x)

x

argument to reluctance or to inductance. For reluctance it is a list of one or more arrays; for instance, one array c(mu,l,A) where mu is relative permeability, l is path length, and A is cross-area. For inductance it is a list(N,rel) where N is number of turns and rel is reluctance For flux it is a list(N,i,rel)where N is number of turns, i is current, and rel is reluctance

Function reluctance calculates reluctance in MA-turn/Wb from mu relative permeability, l path length, and A cross-area. Function inductor calculates inductance in mH from number of turns N and reluctance Function flux calculates magnetic flux from number of turns, current, and reluctance

`rel`	reluctance value from reluctance, inductor, or flux functions
`prnt`	printout of reluctance,inductance, or flux including value and units
`L`	inductance value from inductor or flux functions
`mmf`	magnetomotive force or mmf
`flux`	magnetic flux value

Functions used in Chapters 5 and 10 of Acevedo (2018)

Miguel F. Acevedo Acevedo@unt.edu

Acevedo, M.F. 2018. Introduction to Renewable Electric Power Systems and the Environment with R. Boca Raton, FL: CRC Press. (ISBN 9781138197343)

Irwin, J.D. and R.M. Nelms. 2011. Basic Engineering Circuit Analysis. 11th edition. 2011: Wiley.

# Chapter 5 simple inductor
mu=1000; l=1*10^-2; A=pi*(1*10^-3)^2
reluc <- reluctance(x=list(c(mu,l,A)))
reluc$prnt
ind <- inductor(x=list(N=159,rel=reluc$rel))
ind$prnt

# Chapter 10 air gaped inductor
mucore=1000; lcore=6*10^-2; Acore=1*10^-4
mugap=1; lgap=1*10^-3; Agap <- Acore
reluc <- reluctance(x=list(c(mucore,lcore,Acore),c(mugap,lgap,Agap)))
reluc$prnt

rel.eq <- reluc$rel[1]+reluc$rel[2] 
ind <- inductor(x=list(N=100,rel=rel.eq))
ind$prnt

magckt <- flux(x=list(N=10,i=1,rel=rel.eq))
magckt$prnt

magckt <- flux(x=list(N=10,i=1,r=reluc$r[2]))
magckt$prnt