*heading
Testing basic EM capability - Expected results: A_r + i A_i
Where A_r = A_i = -y e_1 + x e_2
*preprint,model=yes,history=yes
***restart,write,freq=100
*parameter
cond_fct= 1.e-8
cond_cu = 58.0e6
conduct = cond_fct*cond_cu
** mesh related
nelsx = 10
nbl   = 1
nbr   = nelsx+1
**nbr2  = 2*nbr
ntl   = nbl+nelsx*nbr
ntr   = ntl+nbr-1
**
nblpx = nbl+1
nblpy = nbl+nbr
nbrmx = nbr-1
nbrpy = nbr+nbr
ntrmx = ntr-1
ntrmy = ntr-nbr
ntlpx = ntl+1
ntlmy = ntl-nbr
**
n11 = nbl
n12 = n11+1
n13 = n12+nbr
n14 = n13-1
**
x1 = 0.0
y1 = 0.0
x2 = 10.0
y2 = 0.0
x3 = 10.0
y3 = 10.0
x4 = 0.0
y4 = 10.0
**
*node,nset=nbl
<nbl>, <x1>,<y1>
*node,nset=nbr
<nbr>, <x2>,<y2>
*nfill,nset=nbot
nbl,nbr,<nelsx>,1
*node,nset=ntl
<ntl>, <x4>,<y4>
*node,nset=ntr
<ntr>, <x3>,<y3>
*nfill,nset=ntop
ntl,ntr,<nelsx>,1
*nfill,nset=nall
nbot,ntop,<nelsx>,<nbr>
**
*element,type=emc2d3
1, <n11>,<n12>,<n13>
*elgen,elset=right
1, <nelsx>,1,1, <nelsx>,<nbr>,<nelsx>
*element,type=emc2d3
101, <n11>,<n13>,<n14>
*elgen,elset=left
101, <nelsx>,1,1, <nelsx>,<nbr>,<nelsx>
*elset,elset=one
right,left
*elset,elset=bcy0,generate
1,10,1
*elset,elset=bcx1,generate
10,100,10
*elset,elset=bcy1,generate
191,200,1
*elset,elset=bcx0,generate
101,191,10
*surface,name=bcxy0
bcx0,s3
bcy0,s1
*surface,name=bcx1
bcx1,s2
*surface,name=bcy1
bcy1,s2
**
*solid section,elset=one,material=mone
*material,name=mone
*density
1.
*dielectric
8.854187817e-12
*electrical conductivity
<conduct>
*magnetic permeability
1.2566370614e-6
**
*step
*electromagnetic,low frequency, time harmonic
60.d0
*d em potential
bcxy0,s,mvp,, 0.d0, 1.,1.,1.,0.d0, 1.,1.,1.
bcx1, s,mvp,,10.d0, 0.,1.,0.,10.d0, 0.,1.,0.
bcy1, s,mvp,,10.d0,-1.,0.,0.,10.d0,-1.,0.,0.
*decurrent
one,cjnu, 1.0, 1.,1.,1.
*decurrent,imaginary
one,cjnu, 1.0, 1.,1.,1.
*output,field
*element output
EME,EMB
***decurrent
**emc2d4,CJ,1.25,1.0,1.0
**1.,
**2.
*end step