97单元磁铁边界的问题

lijiliang

以前有贴子说
97单元不适合处理磁铁边界的问题
然后发了段小程序
我计算了三种情况
117单元
97单元未处理边界
97单元处理边界

我发现处理边界后效果反而没有原先的好
例子比较简单
有什么问题希望大家多多交流
/prep7
et,1,97,1
et,2,97
et,3,97,1

mp,murx,1,1
mp,rsvx,1,1e-4
mp,murx,3,1000
mp,rsvx,3,1e-5
mp,murx,2,1

blc5,,,1,1,1
blc5,,,2,2,1
blc5,,,4,4,1

vsel,all
vovlap,all
numcmp,all

vsel,s,,,1
vatt,1,1,1
vsel,s,,,2
vatt,2,2,2
vsel,s,,,3
vatt,3,3,3

vsel,all
vsel,s,,,1,2
aslv
lsla
lesize,all,,,10
vmesh,1
vsel,all
vsel,s,,,2,3
mshape,1,3d
mshkey,0
vmesh,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vsel,s,,,3
asel,s,ext
lsla,s,ext
nsll,s,1
cm,mubdry2n,node
nsla,s,all
norl,all,,1                                                 !这一条是我加上的，因为相对于下面nora的那一指令
cmse,u,mubdry2n
cm,mubdry1n,node
nora,all,1
alls
*get,nmax,node,,num,max
vsel,s,,,3
eslv,s,1
nsle
cm,scrap_e,elem
cm,scrap_n,node
ngen,2,nmax,all
egen,2,nmax,all
cmse,s,scrap_n
cmse,s,scrap_e
VCLEAR,     all
  
vdele,all,,,1   
alls
numc,node
alls                ! UNSELECT "FLAT SURFACE" NODES, LEAVING CORNER
cmse,u,mubdry2n        ! AND EDGE NODES. COUPLE AX/Y/Z   ？？？？？这一条为什么不是cmse,s,mubdry2n
cpintf,ax,1e-8
cpintf,ay,1e-8
cpintf,az,1e-8
alls                ! UNSELECT EDGE/CORNER NODES, LEAVING ONLY FLAT
cmse,u,mubdry1n        ! SURFACE NODES. COUPLE AY/Z ONLY？？？？？这一条为什么不是cmse,s,mubdry1n
cpintf,ay,1e-8
cpintf,az,1e-8

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!1
allsel,all
esel,s,mat,,1
nsle
nsel,r,loc,z,1
d,all,volt,0

nsle
nsel,r,loc,z,0
cp,1,volt,all
*get,nmin,node,,num,min
f,nmin,amps,10000


asel,all
asel,s,ext
da,all,asym
finish
/solu
antype,harmic
harfrq,5000
allsel,all
solve
!!!!!!!!!!!!!!!!!!!!!!!!
/prep7
et,1,117,1
et,2,117
et,3,117,1

mp,murx,1,1
mp,rsvx,1,1e-4
mp,murx,3,1000
mp,rsvx,3,1e-5
mp,murx,2,1

blc5,,,1,1,1
blc5,,,2,2,1
blc5,,,4,4,1

vsel,all
vovlap,all
numcmp,all

vsel,s,,,1
vatt,1,1,1
vsel,s,,,2
vatt,2,2,2
vsel,s,,,3
vatt,3,3,3

vsel,all
vsel,s,,,1,2
aslv
lsla
lesize,all,,,10
vmesh,1
vsel,all
vsel,s,,,2,3
mshape,1,3d
mshkey,0
vmesh,all

allsel,all
esel,s,mat,,1
nsle
nsel,r,loc,z,1
d,all,volt,0

nsle
nsel,r,loc,z,0
cp,1,volt,all
*get,nmin,node,,num,min
f,nmin,amps,10000


asel,all
asel,s,ext
da,all,az,0
finish
/solu
antype,harmic
harfrq,5000
allsel,all
solve

[ 本帖最后由 infrom 于 2008-8-27 10:20 编辑 ]

infrom

在http://forum.simwe.com/thread-774283-1-1.html中

5.关于solid97铁和空气（两种磁导率相差很大的物质）交界面上的处理,solid97算得不准，怎么让他算得准呢？
原因：3D节点 mvp 不准。原因：从电磁学中知道，两种物质交接面处的B法向连续，H切向连续。但是由于A的单值，B=curl（A）的法向和切向都是连续的。

给出的solid97在铁磁材料和空气交界面上计算不准的理由，从原理上是不成立的。
原理上两种不同媒质分界面上的边界条件如果用场量表示的话：B法向连续，H切向连续
而若用计算辅助量矢量磁位A表示的话，与场量相对应的分别是：A矢量相等，A矢量的旋度的切向分量乘上对应媒质的磁阻率是相等的。

lijiliang

那就是97单元本身没有问题可以对这种电磁环境模拟
那个帖子是错误的