*heading
Testing basic EM capability - Expected results: A_r + i A_i
Where A_r = A_i = -y e_1 + x e_2
*parameter
mu0 = 1.25663706144e-6
omu = 1.0/mu0
tmu = 2.0*omu
*preprint,model=yes,history=yes
***restart,write,freq=100
*parameter
cond_fct= 1.e-8
cond_cu = 58.0e6
conduct = cond_fct*cond_cu
**
*node,nset=nall
7, 0.0, 0.0
8, 5.0, 0.0
9,10.0, 0.0
1, 0.0, 5.0
2, 5.0, 5.0
3,10.0, 5.0
4, 0.0,10.0
5, 5.0,10.0
6,10.0,10.0
**
*element,type=emc2d4,elset=emc2d4
1, 7,8,2,1
2, 8,9,3,2
3, 1,2,5,4
4, 2,3,6,5
*elset,elset=one
emc2d4
**
*orientation,name=oricyl,system=cylindrical
0.,0.,0., 0.,0.,1.
3,
*solid section,elset=one,material=mone
*material,name=mone
*density
1.
*dielectric
8.854187817e-12
*magnetic permeability
1.2566370614e-6
*surface,name=bcxy0
1,s1
2,s1
1,s4
3,s4
*surface,name=bcx1
2,s2
4,s2
*surface,name=bcy1
3,s3
4,s3
**
*step
*Magnetostatic
1.,1.
*d em potential
bcxy0,s,mvp,,0.d0, 1.,1.,1.
*dsecurrent
bcy1,ck,<tmu>, -1.,0.,0.
bcx1,ck,<tmu>,0.,1.,0.
*output,field
*element output
EMB,EMH
***decurrent
**emc2d4,CJ,1.25,1.0,1.0
**1.,
**2.
*end step