铁球在匀强磁场中的模拟

huangx06

二维情况，铁球在匀强磁场中。采用fmagbc宏施加电磁力。 发现结果与铁球外周一圈的单元密度有关，这个是咋回事？？   就是看后处理的电磁力都不一样，变形也不一样。

附命令流：
/PREP7  
ET,1,PLANE13

KEYOPT,1,1,4
KEYOPT,1,2,0
KEYOPT,1,3,2
KEYOPT,1,4,0
KEYOPT,1,5,0



MP,EX,2,200e9
MP, PRXY, 2, 0.3
MP, MURX, 2, 1e0
MP, MURX, 1, 1
MP, EX,1,1
MP, PRXY,1,0.3

RECTNG,0,0.1,0,0.1,
CYL4,5e-2,5e-2,5e-3,
CYL4,5e-2,5e-2,5.1e-3,
AOVLAP, ALL

ASEL,S,AREA,,2
AATT,2
ASEL,S,AREA,,5
AATT,1
ASEL,S,AREA,,4
AATT,1

LSEL,S,LINE,,1,4,1
LESIZE,ALL,,,24
LSEL,S,LINE,,5,8,1
LESIZE,ALL,,,20       ！这两个是一圈的单元数，如果改变结果就不一样
LSEL,S,LINE,,9,12,1
LESIZE,ALL,,,20       ！这两个是一圈的单元数，如果改变结果就不一样


ASEL,S,AREA,,2
AMESH,ALL
ASEL,S,AREA,,4
AMESH,ALL
ASEL,S,AREA,,5
AMESH,ALL



ALLSEL

!施加磁场
NSEL, S, LOC, X, 0
D, ALL, AZ, 1e-4   
NSEL, S, LOC, X, 0.1
D, ALL, AZ, 0

!!电磁宏加力
esel,s,mat,,2                        ! armature
cm,arm,elem
fmagbc,'arm'




ALLSEL
FINI

/SOLU  
NLGEOM,ON
ANTYPE, STATIC
SOLV

网格份数为20的时候
位移图：
网格份数为100的时候就不变形了。没有位移。

而且两者的电磁力结果也不一样，这是咋回事？

huangx06

自己顶一下。
有没有谁做过这方面的？ 这问题困扰我好久了。
网格很密的时候也会是不同网格有不同的电磁力，这样子怎么敢相信结果

huangx06

怎么没有人呢

huangx06

求高人解答一下