铁球在匀强磁场中的模拟
二维情况，铁球在匀强磁场中。采用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的时候就不变形了。没有位移。
而且两者的电磁力结果也不一样，这是咋回事？
自己顶一下。
有没有谁做过这方面的？ 这问题困扰我好久了。
网格很密的时候也会是不同网格有不同的电磁力，这样子怎么敢相信结果 怎么没有人呢 求高人解答一下
页:
[1]