*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
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
17, 0.0, 0.0, 2.0
18, 5.0, 0.0, 2.0
19,10.0, 0.0, 2.0
11, 0.0, 5.0, 2.0
12, 5.0, 5.0, 2.0
13,10.0, 5.0, 2.0
14, 0.0,10.0, 2.0
15, 5.0,10.0, 2.0
16,10.0,10.0, 2.0
**
*element,type=emc3d8,elset=emc3d8
1, 7,8,2,1,17,18,12,11
2, 8,9,3,2,18,19,13,12
3, 1,2,5,4,11,12,15,14
4, 2,3,6,5,12,13,16,15
*elset,elset=one
emc3d8
**
*surface,name=bcxy0
1,s3
2,s3
1,s6
3,s6
*surface,name=bcx1
2,s4
4,s4
*surface,name=bcy1
3,s5
4,s5
*solid section,elset=one,material=mone
*material,name=mone
*density
1.
*dielectric
8.854187817e-12
*electrical conductivity
<conduct>
*magnetic permeability
1.2566370614e-6
**
*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.
*decurrent,imaginary
one,cjnu, 1.0, 1.,1.,1.
*output,field
*element output
EME,EMB,
EMCDA
***decurrent
**emc3d8,CJ,1.25,1.0,1.0
**1.,
**2.
*end step