I have successfully installed the mumax3 coupled with magnetoelasticity module on my computer. I really want to simulate the effect of strain on the hysteresis loop. However, the hysteresis loops are exactly the same under different strains event I run simultaneously various strains at a time or separately.
Here, I would like to attach the code i used. Could you please do me a favor by checking it.
With all the best,
A++
SetGridsize(128, 32, 4)
SetCellsize(4e-9, 4e-9, 4e-9)
setgeom(cuboid(512e-9, 128e-9, 16e-9))
Msat = 1068918 //ok
Aex = 1e-11 //ok
alpha = 0.05
m=Uniform(1,0,0)
u=Uniform(1,0,0)
du=Uniform (1,0,0)
//Elastics parameters
c11 = 229e9 //ok
c12 = 134e9 //ok
c44 = 47.5e9 //ok, c44=(c11_c12)/2
rho = 7.2e3 //ok
eta = 0
//Magnetoelastic parameters
B1 = -9.975e5 // ok, lamda_100=-2/3*((B1)/(c11-c12)) and lamda_111=-1/3*((B2)/(c44)), lamda_100= lamda_111=lamda_s=-7e6
B2 = 9.975e5
FixedLayer = vector(1,0,0)
frozenDispLoc = 1
frozenDispVal = vector(1,0,0)
FreeLayerThickness = 15e-9
Bmax := 100.0e-3
Bstep := 2.0e-3
TableAdd(B_ext)
for B:=Bmax; B>2.0e-3; B-=5*Bstep{
ezz = -0.3
exx = 0.15
eyy = 0.15
B_ext = vector(B, 0, B*0.001)
tablesave()
save(m_full)
save(u)
save(ezz)
save(normstrain)
//SetSolver(5)
fixdt = 5e-13
//running
run(1e-9)
}
for B:=0.0; B>=-Bmax; B-=Bstep{
ezz = -0.3
exx = 0.15
eyy = 0.15
B_ext = vector(B, 0, B*0.001)
tablesave()
save(m_full)
save(u)
save(ezz)
save(normstrain)
//SetSolver(5)
fixdt = 5e-13
//running
run(1e-9)
}
for B:=Bmax; B>2.0e-3; B-=5*Bstep{
ezz = 0
exx = 0
eyy = 0
B_ext = vector(B, 0, B*0.001)
tablesave()
save(m_full)
save(u)
save(ezz)
save(normstrain)
//SetSolver(5)
fixdt = 5e-13
//running
run(1e-9)
}
for B:=0.0; B>=-Bmax; B-=Bstep{
ezz = 0
exx = 0
eyy = 0
B_ext = vector(B, 0, B*0.001)
tablesave()
save(m_full)
save(u)
save(ezz)
save(normstrain)
//SetSolver(5)
fixdt = 5e-13
//running
run(1e-9)
}
for B:=Bmax; B>2.0e-3; B-=5*Bstep{
ezz = 0.3
exx = -0.15
eyy = -0.15
B_ext = vector(B, 0, B*0.001)
tablesave()
save(m_full)
save(u)
save(ezz)
save(normstrain)
//SetSolver(5)
fixdt = 5e-13
//running
run(1e-9)
}
for B:=0.0; B>=-Bmax; B-=Bstep{
ezz = 0.3
exx = -0.15
eyy = -0.15
B_ext = vector(B, 0, B*0.001)
tablesave()
save(m_full)
save(u)
save(ezz)
save(normstrain)
//SetSolver(5)
fixdt = 5e-13
//running
run(1e-9)
}