*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
conduct_0 = 0.5*conduct
conduct_1 = 1.5*conduct
mu      = 1.2566370614e-6
mu2     = 2.0*mu
mu3     = 3.0*mu
omu     = 1.0/mu
tomu    = 2.0*omu
thomu   = 3.0*omu
mu_0    = 0.5*mu
mu_1    = 1.5*mu
*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_0>,0.
<conduct_1>,200.
*magnetic permeability
<mu_0>,0.
<mu_1>,200.
*material,name=mnlbh
*density
1.
*electrical conductivity
<conduct_0>,0.
<conduct_1>,200.
*magnetic permeability, nonlinear
*nonlinear BH, dir=1
0., 0., 0.
500.0, <omu>, 0.
1000.0, <tomu>, 0.
1200.0, <thomu>, 0.
0., 0., 200.
1500.0, <omu>, 200.
2000.0, <tomu>, 200.
2200.0, <thomu>, 200.
*material,name=mortho
*density
1.
*electrical conductivity
<conduct_0>,0.
<conduct_1>,200.
*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., 0.
500.0, <omu>, 0.
1000.0, <tomu>, 0.
1200.0, <thomu>, 0.
0., 0., 200.
1500.0, <omu>, 200.
2000.0, <tomu>, 200.
2200.0, <thomu>, 200.
*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.
*temperature
nall,100.
*output,field,freq=1
*element output
EMB,emh,
**eme,emcd,
**emjh,embf
*end step