24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1377  |  回复: 5

xuuuujiawen

新虫 (初入文坛)

[交流] mumax仿真添加周期性边界条件之后磁滞回线不对称怎么解决?

各位大佬,我最近碰到了一个问题,我在mumax上模拟致密六角钡铁氧体薄膜使用周期性边界条件之后,模拟出的磁滞回线不对称,我之前没有遇到过类似的问题,不清楚是什么原因导致的。我设置了晶粒间和晶粒内的交换作用常数并且考虑了模拟边界带来的影响之后,还是没有消除影响,请问可以帮助我分析原因,希望可以和大家交流得到解决方案,谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuuuujiawen

新虫 (初入文坛)

我自己补充一下模型、结果图

发自小木虫IOS客户端
2楼2025-04-22 15:49:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuuuujiawen

新虫 (初入文坛)

以下是我的结果图,应该怎么分析他?求助!!!!感谢大家!
3楼2025-04-22 15:49:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuuuujiawen

新虫 (初入文坛)

4楼2025-04-22 15:56:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuuuujiawen

新虫 (初入文坛)

5楼2025-04-22 15:56:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuuuujiawen

新虫 (初入文坛)

defregion(1, matrix1)
defregion(2, matrix2)
defregion(3, matrix3)
defregion(4, matrix4)
defregion(5, matrix5)
defregion(6, matrix6)
defregion(7, matrix7)
defregion(8, matrix8)
defregion(9, matrix9)
defregion(10, matrix10)
defregion(11, matrix11)
defregion(12, matrix12)
defregion(13, matrix13)
defregion(14, matrix14)
defregion(15, matrix15)
defregion(16, matrix16)
defregion(17, matrix17)
save(regions)

Msat. setRegion(1, 275e3)
Aex. setRegion(1, 20e-12)
Ku1. setRegion(1, 96e3)
AnisU. setRegion(1, vector(0, 0, 1))
alpha. setRegion(1, 0.2)
Temp. setRegion(1, 300)
m. setRegion(1, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(2, 275e3)
Aex. setRegion(2, 20e-12)
Ku1. setRegion(2, 96e3)
AnisU. setRegion(2, vector(0, 0, 1))
alpha. setRegion(2, 0.2)
Temp. setRegion(2, 300)
m. setRegion(2, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(3, 275e3)
Aex. setRegion(3, 20e-12)
Ku1. setRegion(3, 96e3)
AnisU. setRegion(3, vector(0, 0, 1))
alpha. setRegion(3, 0.2)
Temp. setRegion(3, 300)
m. setRegion(3, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(4, 275e3)
Aex. setRegion(4, 20e-12)
Ku1. setRegion(4, 96e3)
AnisU. setRegion(4, vector(0, 0, 1))
alpha. setRegion(4, 0.2)
Temp. setRegion(4, 300)
m. setRegion(4, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(5, 275e3)
Aex. setRegion(5, 20e-12)
Ku1. setRegion(5, 96e3)
AnisU. setRegion(5, vector(0, 0, 1))
alpha. setRegion(5, 0.2)
Temp. setRegion(5, 300)
m. setRegion(5, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(6, 275e3)
Aex. setRegion(6, 20e-12)
Ku1. setRegion(6, 96e3)
AnisU. setRegion(6, vector(0, 0, 1))
alpha. setRegion(6, 0.2)
Temp. setRegion(6, 300)
m. setRegion(6, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(7, 275e3)
Aex. setRegion(7, 20e-12)
Ku1. setRegion(7, 96e3)
AnisU. setRegion(7, vector(0, 0, 1))
alpha. setRegion(7, 0.2)
Temp. setRegion(7, 300)
m. setRegion(7, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(8, 275e3)
Aex. setRegion(8, 20e-12)
Ku1. setRegion(8, 96e3)
AnisU. setRegion(8, vector(0, 0, 1))
alpha. setRegion(8, 0.2)
Temp. setRegion(8, 300)
m. setRegion(8, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(9, 275e3)
Aex. setRegion(9, 20e-12)
Ku1. setRegion(9, 96e3)
AnisU. setRegion(9, vector(0, 0, 1))
alpha. setRegion(9, 0.2)
Temp. setRegion(9, 300)
m. setRegion(9, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(10, 275e3)
Aex. setRegion(10, 20e-12)
Ku1. setRegion(10, 96e3)
AnisU. setRegion(10, vector(0, 0, 1))
alpha. setRegion(10, 0.2)
Temp. setRegion(10, 300)
m. setRegion(10, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(11, 275e3)
Aex. setRegion(11, 20e-12)
Ku1. setRegion(11, 96e3)
AnisU. setRegion(11, vector(0, 0, 1))
alpha. setRegion(11, 0.2)
Temp. setRegion(11, 300)
m. setRegion(11, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(12, 275e3)
Aex. setRegion(12, 20e-12)
Ku1. setRegion(12, 96e3)
AnisU. setRegion(12, vector(0, 0, 1))
alpha. setRegion(12, 0.2)
Temp. setRegion(12, 300)
m. setRegion(12, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(13, 275e3)
Aex. setRegion(13, 20e-12)
Ku1. setRegion(13, 96e3)
AnisU. setRegion(13, vector(0, 0, 1))
alpha. setRegion(13, 0.2)
Temp. setRegion(13, 300)
m. setRegion(13, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(14, 275e3)
Aex. setRegion(14, 20e-12)
Ku1. setRegion(14, 96e3)
AnisU. setRegion(14, vector(0, 0, 1))
alpha. setRegion(14, 0.2)
Temp. setRegion(14, 300)
m. setRegion(14, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(15, 275e3)
Aex. setRegion(15, 20e-12)
Ku1. setRegion(15, 96e3)
AnisU. setRegion(15, vector(0, 0, 1))
alpha. setRegion(15, 0.2)
Temp. setRegion(15, 300)
m. setRegion(15, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(16, 275e3)
Aex. setRegion(16, 20e-12)
Ku1. setRegion(16, 96e3)
AnisU. setRegion(16, vector(0, 0, 1))
alpha. setRegion(16, 0.2)
Temp. setRegion(16, 300)
m. setRegion(16, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Msat. setRegion(17, 275e3)
Aex. setRegion(17, 20e-12)
Ku1. setRegion(17, 96e3)
AnisU. setRegion(17, vector(0, 0, 1))
alpha. setRegion(17, 0.2)
Temp. setRegion(17, 300)
m. setRegion(17, Uniform(sin(pi/36)*Rand(), sin(pi/36)*sqrt(1-(Rand()*Rand())), cos(pi/36)))

Aintra :=20e-13
ext_InterExchange(1,2,Aintra)
ext_InterExchange(1,3,Aintra)
ext_InterExchange(1,4,Aintra)
ext_InterExchange(1,5,Aintra)
ext_InterExchange(1,6,Aintra)
ext_InterExchange(1,7,Aintra)
ext_InterExchange(2,3,Aintra)
ext_InterExchange(2,7,Aintra)
ext_InterExchange(2,8,Aintra)
ext_InterExchange(2,9,Aintra)
ext_InterExchange(3,9,Aintra)
ext_InterExchange(3,10,Aintra)
ext_InterExchange(3,11,Aintra)
ext_InterExchange(3,4,Aintra)
ext_InterExchange(4,5,Aintra)
ext_InterExchange(4,11,Aintra)
ext_InterExchange(4,12,Aintra)
ext_InterExchange(4,13,Aintra)
ext_InterExchange(5,6,Aintra)
ext_InterExchange(5,13,Aintra)
ext_InterExchange(5,14,Aintra)
ext_InterExchange(6,7,Aintra)
ext_InterExchange(6,14,Aintra)
ext_InterExchange(6,15,Aintra)
ext_InterExchange(6,16,Aintra)
ext_InterExchange(7,8,Aintra)
ext_InterExchange(7,16,Aintra)
ext_InterExchange(7,17,Aintra)
ext_InterExchange(8,17,Aintra)
ext_InterExchange(9,10,Aintra)
ext_InterExchange(10,11,Aintra)
ext_InterExchange(11,12,Aintra)
ext_InterExchange(12,13,Aintra)
ext_InterExchange(14,15,Aintra)
ext_InterExchange(15,16,Aintra)
ext_InterExchange(16,17,Aintra)

ext_InterExchange(8,5,Aintra)
ext_InterExchange(8,6,Aintra)
ext_InterExchange(8,12,Aintra)
ext_InterExchange(8,14,Aintra)
ext_InterExchange(8,15,Aintra)

ext_InterExchange(9,4,Aintra)
ext_InterExchange(9,5,Aintra)
ext_InterExchange(9,12,Aintra)
ext_InterExchange(9,13,Aintra)
ext_InterExchange(9,15,Aintra)

ext_InterExchange(17,10,Aintra)
ext_InterExchange(17,11,Aintra)
ext_InterExchange(17,3,Aintra)
ext_InterExchange(17,9,Aintra)
ext_InterExchange(17,13,Aintra)
ext_InterExchange(17,12,Aintra)
ext_InterExchange(17,15,Aintra)
ext_InterExchange(17,14,Aintra)

ext_InterExchange(2,14,Aintra)
ext_InterExchange(2,5,Aintra)
ext_InterExchange(2,13,Aintra)

ext_InterExchange(10,17,Aintra)
ext_InterExchange(10,16,Aintra)
ext_InterExchange(10,7,Aintra)
ext_InterExchange(10,8,Aintra)
ext_InterExchange(10,14,Aintra)
ext_InterExchange(10,15,Aintra)
ext_InterExchange(10,12,Aintra)
ext_InterExchange(10,13,Aintra)

ext_InterExchange(7,10,Aintra)
ext_InterExchange(7,11,Aintra)
ext_InterExchange(7,14,Aintra)

ext_InterExchange(3,13,Aintra)
ext_InterExchange(3,17,Aintra)
ext_InterExchange(3,16,Aintra)

ext_InterExchange(16,10,Aintra)
ext_InterExchange(16,11,Aintra)
ext_InterExchange(16,12,Aintra)
ext_InterExchange(16,3,Aintra)
ext_InterExchange(16,4,Aintra)

ext_InterExchange(11,15,Aintra)
ext_InterExchange(11,16,Aintra)
ext_InterExchange(11,17,Aintra)
ext_InterExchange(11,6,Aintra)
ext_InterExchange(11,7,Aintra)

ext_InterExchange(6,11,Aintra)
ext_InterExchange(6,12,Aintra)
ext_InterExchange(6,8,Aintra)

ext_InterExchange(4,15,Aintra)
ext_InterExchange(4,16,Aintra)
ext_InterExchange(4,9,Aintra)

ext_InterExchange(15,12,Aintra)
ext_InterExchange(15,11,Aintra)
ext_InterExchange(15,4,Aintra)
ext_InterExchange(15,9,Aintra)
ext_InterExchange(15,13,Aintra)
ext_InterExchange(15,10,Aintra)
ext_InterExchange(15,17,Aintra)
ext_InterExchange(15,8,Aintra)

ext_InterExchange(5,8,Aintra)
ext_InterExchange(5,2,Aintra)
ext_InterExchange(5,9,Aintra)

ext_InterExchange(12,15,Aintra)
ext_InterExchange(12,16,Aintra)
ext_InterExchange(12,6,Aintra)
ext_InterExchange(12,14,Aintra)
ext_InterExchange(12,8,Aintra)
ext_InterExchange(12,17,Aintra)
ext_InterExchange(12,9,Aintra)
ext_InterExchange(12,10,Aintra)

ext_InterExchange(14,8,Aintra)
ext_InterExchange(14,17,Aintra)
ext_InterExchange(14,7,Aintra)
ext_InterExchange(14,2,Aintra)
ext_InterExchange(14,10,Aintra)

ext_InterExchange(13,9,Aintra)
ext_InterExchange(13,2,Aintra)
ext_InterExchange(13,3,Aintra)
ext_InterExchange(13,10,Aintra)
ext_InterExchange(13,17,Aintra)

发自小木虫IOS客户端
6楼2025-04-22 15:57:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xuuuujiawen 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见