*heading
Low Freq Time Domain EM that mimicks
Nonlinear B in Magnetostatics with
Quads and Triangles, and NLBH Iso and Ortho
***preprint,model=yes,history=yes
***restart,write,freq=100
*parameter
cond_fct= 1.e-12
cond_cu = 58.0e6
conduct = cond_fct*cond_cu
mu      = 1.2566370614e-6
mu2     = 2.0*mu
mu3     = 3.0*mu
omu     = 1.0/mu
tomu    = 2.0*omu
thomu   = 3.0*omu
*node,nset=nall
 1, 0., 0.
 2, 5., 0.
 3,10., 0.
 4,15., 0.
 5,20., 0.
 6,25., 0.
 7,30., 0.
 8,35., 0.
 9,40., 0.
11, 0., 5.
12, 5., 5.
13,10., 5.
14,15., 5.
15,20., 5.
16,25., 5.
17,30., 5.
18,35., 5.
19,40., 5.
21, 0.,10.
22, 5.,10.
23,10.,10.
24,15.,10.
25,20.,10.
26,25.,10.
27,30.,10.
28,35.,10.
29,40.,10.
*element,type=emc2d3,elset=enlbh
 1,  1,12,11
 2,  1, 2,12
*element,type=emc2d3,elset=eortho
 7,  7,18,17
 8,  7, 8,18
*element,type=emc2d3,elset=elbh
 5,  5,16,15
 6,  5, 6,16
12, 12,23,22
13, 12,13,23
14, 14,25,24
15, 14,15,25
18, 18,29,28
19, 18,19,29
*element,type=emc2d4,elset=enlbh
 4,  4, 5,15,14
*element,type=emc2d4,elset=elbh
 3,  2, 3,13,12
 9,  8, 9,19,18
11, 11,12,22,21
16, 15,16,26,25
17, 17,18,28,27
*surface,name=bcxy0
***x0
 1,s3
 4,s4
 7,s3
11,s4
14,s3
17,s4
***y0
 2,s1
 3,s1
 4,s1
 6,s1
 8,s1
 9,s1
*surface,name=bcx1
 3,s2
13,s2
 6,s2
16,s2
 9,s2
19,s2
*surface,name=bcy1
11,s3
12,s2
14,s2
16,s3
17,s3
18,s2
**
*solid section,elset=elbh,material=mone
*solid section,elset=enlbh,material=mnlbh
*solid section,elset=eortho,material=mortho,ori=ori1
*material,name=mone
*density
1.
*dielectric
8.854187817e-12
*electrical conductivity
<conduct>
*magnetic permeability
<mu>
*material,name=mnlbh
*density
1.
*electrical conductivity
<conduct>
*magnetic permeability, nonlinear
*nonlinear BH, dir=1
0., 0.
1000.0, <omu>
1500.0, <tomu>
1700.0, <thomu>
*material,name=mortho
*density
1.
*electrical conductivity
<conduct>
*magnetic permeability, nonlinear, type=ortho
*nonlinear BH, dir=1
0.,0.
1.,<omu>
*nonlinear BH, dir=2
0.,0.
1.,<omu>
*nonlinear BH, dir=3
0., 0.
1000.0, <omu>
1500.0, <tomu>
1700.0, <thomu>
*orientation,name=ori1
1.,0.,0., 0.,1.,0.
3,0.
*step,amp=ramp
*electromagnetic,low frequency, transient
**, timeint=b-euler
0.1,1.,,0.1
*d em potential
bcxy0,s,mvp,, 0.d0, 1.,0.,0.
*dsecurrent
bcx1,ck, <tomu>, 0.,1.,0.
bcy1,ck, <tomu>,-1.,0.,0.
*output,field,freq=1
*element output
EMB,emh,
**eme,emcd,
**emjh,embf
*end step