找回密码
 注册
Simdroid-非首页
查看: 72|回复: 0

[电磁场] 请教一下大家,关于基于低频电磁的管道漏磁缺陷检测的问题

[复制链接]
发表于 2017-11-6 13:42:06 | 显示全部楼层 |阅读模式 来自 湖北武汉
下面是我的命令流,有几个问题:1、本来铁芯与样板都是q235钢,但因为是谐态的,材料不能为非线性,在网上找的将其材料改为磁导率4000  2、样板的单元选项是AZ+volt,在对成面进行volt耦合,不知道这样对不 3、网格的问题,sweep划分时与free差别很大,为嘛、
最重要的是想知道这样命令流有没什么问题,能力有限,感觉没问题,但是路径出来的结果不对
/title, 3-D  Harmonic Response, 1/2 symmetry
/vie,1,2,-1,3
/pnu,mat,1
/num,1


! *** Model parameters
a1_core=0.02        ! core cross-section width
a2_core=0.04       ! core cross-section depth
w1_core=0.08       ! medium   core width
w2_core=a1_core*2+w1_core       ! overall core width
h_core=0.08       ! overall core height

d_core=0


thk_coil1=0.01               !left right thick
thk_coil2=0.015               !UP thick  
h_coil1=0.04       ! height,left right coils
h_coil2=0.058      ! height,UP coils

d_plate=0.001          !plate
l_plate=0.5
h_plate=0.01
w_plate=0.2


l_flaw=0.0005         !length
h_flaw=0.002      !depth
w_flaw=0.005



d_dmn=0.025         ! depth of surrounding domain

esz1=a1_core/2       ! element size, components
esz2=2*esz1         ! element size, surrounding domain

esz3=thk_coil2/3
esz4=h_flaw/3
esz5=h_plate/2



rs_q235=0.12E-06

rs_coil=1.724e-8

! *** Primary coil
I=10           !  current (A)

N1=100            ! # of turns
R1=0.0668              ! DC resistance (ohms)
care_left=thk_coil1*h_coil1      ! cross-sectional area
volu_left=((thk_coil1*2+a1_core)*(thk_coil1*2+a2_core)-a1_core*a2_core)*h_coil1 ! volume left
js1=n1*I/care_left

N2=200            ! # of turns
R2=0.1336              ! DC resistance (ohms)
care_up=thk_coil2*h_coil2      ! cross-sectional area
volu_up=((thk_coil2*2+a1_core)*(thk_coil2*2+a2_core)-a1_core*a2_core)*h_coil2  ! volume up
js2=n2*I/care_up


frqncy=50               ! operating frequency (Hz)
symm=2                 ! symmetry factor


/nopr
! *** Geometry
/PREP7
wpcs,-1,0
bloc,-(a1_core*2+w1_core)/2,(a1_core*2+w1_core)/2,h_core-a1_core,h_core,-a2_core/2,0  !core
bloc,-(a1_core*2+w1_core)/2,-(a1_core*2+w1_core)/2+a1_core,0,h_core-a1_core,-a2_core/2,0
bloc,(a1_core*2+w1_core)/2-a1_core,(a1_core*2+w1_core)/2,0,h_core-a1_core,-a2_core/2,0



!up coil

wpcs,-1,0                ! up coil ESYS1
wpof,,h_core-a1_core/2
wpro,,,90
cswp,13,1

csys               !up coil
wpcs,-1,13
bloc,0,a2_core/2,-(a1_core/2+thk_coil2),-a1_core/2,-h_coil2/2,h_coil2/2,
bloc,a2_core/2,a2_core/2+thk_coil2,-a1_core/2,a1_core/2,-h_coil2/2,h_coil2/2,
bloc,0,a2_core/2,a1_core/2,(a1_core/2+thk_coil2),-h_coil2/2,h_coil2/2,


wpcs,-1,13               ! up coil ESYS2
wpof,a2_core/2,-a1_core/2
cswp,20,1

cylind,0,thk_coil2,-h_coil2/2,h_coil2/2,-90,0,


wpcs,-1,13               ! up coil ESYS3
wpof,a2_core/2,a1_core/2
cswp,21,1

cylind,0,thk_coil2,-h_coil2/2,h_coil2/2,0,90,
  


!left coil

wpcs,-1,0                ! left coil ESYS1
wpof,-(w1_core+a1_core)/2
wpro,,-90,
cswp,11,1

csys               !left coil
wpcs,-1,11
bloc,-a1_core/2,-(a1_core/2+thk_coil1),0,a2_core/2,0,h_coil1,
bloc,-a1_core/2,a1_core/2,a2_core/2,a2_core/2+thk_coil1,0,h_coil1,
bloc,a1_core/2,(a1_core/2+thk_coil1),0,a2_core/2,0,h_coil1,

wpcs,-1,11              ! left coil ESYS2
wpof,-a1_core/2,a2_core/2
cswp,22,1

cylind,0,thk_coil1,0,h_coil1,90,180,


wpcs,-1,11               ! left coil ESYS3
wpof,a1_core/2,a2_core/2
cswp,23,1

cylind,0,thk_coil1,0,h_coil1,0,90,
  
!right coil

wpcs,-1,0                ! right coil ESYS1
wpof,(w1_core+a1_core)/2
wpro,,90,
cswp,12,1

csys               !right coil
wpcs,-1,12
bloc,-a1_core/2,-(a1_core/2+thk_coil1),-a2_core/2,0,-h_coil1,0,
bloc,-a1_core/2,a1_core/2,-(a2_core/2+thk_coil1),-a2_core/2,-h_coil1,0,
bloc,a1_core/2,(a1_core/2+thk_coil1),-a2_core/2,0,-h_coil1,0,

wpcs,-1,12             ! right coil ESYS2
wpof,-a1_core/2,-a2_core/2
cswp,24,1

cylind,0,thk_coil1,-h_coil1,0,-90,-180,


wpcs,-1,12               ! right coil ESYS3
wpof,a1_core/2,-a2_core/2
cswp,25,1

cylind,0,thk_coil1,-h_coil1,0,-90,0,
vglue,all

vsel,none                 ! plate
wpcs,-1,0
bloc,-l_plate/2,l_plate/2,-h_plate-d_plate,-d_plate,-w_plate/2,0
cm,plate_1,volu

vsel,none                                !flaw
bloc,-l_flaw/2,l_flaw/2,-h_flaw-d_plate,-d_plate,-w_flaw/2,0
cm,flaw_v,volu
vatt,1,1,1

cmse,s,plate_1
cmse,a,flaw_v
vsbv,plate_1,flaw_v,,dele,keep
cmse,u,flaw_v
cm,plate_v,volu
vatt,4,4,4





alls
csys
cm,keep_v,volu
*get,xmin,kp,,mnloc,x
*get,xmax,kp,,mxloc,x
*get,ymax,kp,,mxloc,y
*get,ymin,kp,,mnloc,y
*get,zmin,kp,,mnloc,z


wpcs,-1,0                ! surrounding domain

bloc,xmin-d_dmn,xmax+d_dmn,ymin-d_dmn,ymax+d_dmn,zmin-d_dmn,0
allsel,all
vovlap,all
numcmp,all

vsel,s,volu,,1
vatt,4,1,4


vsel,s,volu,,2
vsel,a,volu,,11,18
vatt,3,1,3

vsel,s,volu,,5
vatt,3,1,3,20

vsel,s,volu,,6
vatt,3,1,3,21

vsel,s,volu,,7
vatt,3,1,3,22

vsel,s,volu,,8
vatt,3,1,3,23

vsel,s,volu,,9
vatt,3,1,3,24

vsel,s,volu,,10
vatt,3,1,3,25



vsel,s,volu,,19,21
vatt,2,1,2

vsel,s,volu,,3,4
vatt,1,1,1

alls
! *** FE model
et,1,236                 ! air
mp,murx,1,1

et,2,236                 ! core (laminated, non-conducting)
mp,murx,2,4000


et,3,236              ! left primary coil
mp,murx,3,1
mp,rsvx,3,rs_coil


et,4,236,1                 ! plate
mp,murx,4,4000
mp,rsvx,4,rs_q235










! *** Mesh
numm,kp,1e-8,1e-8

msha,0
esiz,esz3
vsel,s,mat,,3
vswe,all

esiz,esz3
vsel,s,mat,,2
vswe,all

esiz,esz4
vsel,s,volu,,3
vswe,all


msha,1
esiz,esz5
vsel,s,mat,,4
vmes,all




vsel,s,mat,,1
msha,1
esiz,esz2
vmes,all
alls

! *** Boundary conditions and loads
asel,s,ext             ! flux parallel exterior
csys
da,all,az

vsel,s,mat,,4
aslv,s
asel,r,loc,z,0
nsla,s
d,all,volt,0
nsel,all


vsel,s,volu,,2
alls,belo,volu
bfe,all,js,,0,0,-js2

vsel,s,volu,,5,6
alls,belo,volu
bfe,all,js,,0,js2,0

vsel,s,volu,,11
alls,belo,volu
bfe,all,js,,0,js2,0

vsel,s,volu,,12
alls,belo,volu
bfe,all,js,,0,0,js2

vsel,s,volu,,13
alls,belo,volu
bfe,all,js,,0,0,js1

vsel,s,volu,,7,8
alls,belo,volu
bfe,all,js,,0,js1,0

vsel,s,volu,,14
alls,belo,volu
bfe,all,js,,-js1,0,0

vsel,s,volu,,15
alls,belo,volu
bfe,all,js,,0,0,-js1

vsel,s,volu,,16
alls,belo,volu
bfe,all,js,,0,0,-js1

vsel,s,volu,,9,10
alls,belo,volu
bfe,all,js,,0,js1,0

vsel,s,volu,,17
alls,belo,volu
bfe,all,js,,js1,0,0

vsel,s,volu,,18
alls,belo,volu
bfe,all,js,,0,0,js1

csys
alls
eplo
fini

! *** Solution
/solu
antype,harmonic
harf,frqncy
solve
fini

! *** Post-processing
/post1
set,,,,0               ! Real solution set

/com,
vsel,s,mat,,2,4
alls,belo,volu
plve,jt,,,,vect,,on
plve,b,,,,vect,,on
plnsol,emf
alls

alls

path,path1,2,10,40
ppath,1,,-0.01,-0.001,0
ppath,2,,0.01,-0.001,0
pdef,By,B,y,avg
pdef,Bx,B,x,avg

plpath,BY,BX,
prpath,BY,BX,

fini

希望路过的能帮帮忙,谢谢啦
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 20:36 , Processed in 0.029122 second(s), 12 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表