*heading
Low Freq Time Domain EM that mimicks Magnetostatics,
with nonlinear B-response using sine wave excitation
Quads and Triangles
***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
*elset,elset=eout
1,4,7, 13,16,19
*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.
*amplitude, name=sine, def=periodic
**N, omega, t0, A0
1, 314.159265, 0.,0.
**c, s
0., 1.0
*step
*electromagnetic,low frequency, transient
**, timeint=b-euler
0.0005,0.02,,0.0005
*d em potential
bcxy0,s,mvp,, 0.d0, 1.,0.,0.
*dsecurrent, amp=sine
bcx1,ck, <tomu>, 0.,1.,0.
bcy1,ck, <tomu>,-1.,0.,0.
*output,field
*element output
EMB,emh,
***output,history
***element output, elset=eout
**emb
*end step