*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
**
*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
*electrical conductivity,frequency
0.57,50.
0.59,70.
**<conduct>
*magnetic permeability
1.193805208e-6,0.
1.319468914e-6,100.
**1.2566370614e-6
*initial conditions,type=temperature
nall,50.
*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
First LF TH Analysis
*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., oricyl
*decurrent,imaginary
one,cjnu, 1.0, 1.,1.,1., oricyl
*output,field
*element output
EME,EMB,
TEMP,EVOL,
EMCDA
***decurrent
**emc2d4,CJ,1.25,1.0,1.0
**1.,
**2.
*end step