电磁场模拟 组件的建立 请大家指点
需要电流加载在线圈上 不知道哪里出问题了 总是加不到线圈上 请大家多多指点 已经困扰很久了FINISH
/CLE
/title,MAGNETIC PULSING
/filnam,magnetic field analysis,1
t_loadstep=40
dg=0.1
lg=68
dxi=30
dx0=50
h=5
n=8
a=4
gap=5
murx1=1
p1=1.72e-8
murx2=1
p2=1.72e-8
/PREP7
ET,1,PLANE13 !坯料
ET,2,PLANE13 !qudongian
ET,3,PLANE13 !线圈1
ET,4,PLANE13 !线圈2
ET,5,PLANE13 !waidaokuang
ET,6,PLANE13 !空气
KEYOPT,3,3,1
UIMP,1,MURX,,,MURX1, !pi;iao
UIMP,1,RSVX,,,P1,
UIMP,2,MURX,,,MURX1, !qudongpiian
UIMP,2,RSVX,,,P1,
UIMP,3,MURX,,,MURX2, !xianquan 1
UIMP,3,RSVX,,,P2,
UIMP,4,MURX,,,MURX2, !xianquan 2
UIMP,4,RSVX,,,P2,
UIMP,5,MURX,,,MURX2, !waidaokuang
UIMP,5,RSVX,,,P2,
UIMP,6,MURX,,,1, !kong qi
/PNUM,AREA,1
EMUNIT,MKS
RECTNG,-1.0e-3*Lg/2+0.5,1.0e-3*Lg/2+0.5,-1.0e-3*(Dx0/2+gap+0.17)+0.5,-1.0e-3*(Dx0/2+gap+Dg+0.17)+0.5
RECTNG,-1.0e-3*Lg/2+0.5,1.0e-3*Lg/2+0.5,-1.0e-3*(Dx0/2+gap)+0.5,-1.0e-3*(Dx0/2+gap+0.16)+0.5
*DO,I,1,N/2,1
RECTNG,a*1.0e-3/2+(I-1)*(h+a)*1.0e-3+0.5,a*1.0e-3/2+(I*h+(I-1)*a)*1.0e-3+0.5,-(Dxi*1.0e-3)/2+0.5,-(Dx0*1.0e-3)/2+0.5
RECTNG,a*1.0e-3/2+(I-1)*(h+a)*1.0e-3+0.5,a*1.0e-3/2+(I*h+(I-1)*a)*1.0e-3+0.5,(Dxi*1.0e-3)/2+0.5,(Dx0*1.0e-3)/2+0.5
RECTNG,-(a*1.0e-3/2+(I-1)*(h+a)*1.0e-3)+0.5,-(a*1.0e-3/2+(I*h+(I-1)*a)*1.0e-3)+0.5,-(Dxi*1.0e-3)/2+0.5,-(Dx0*1.0e-3)/2+0.5
RECTNG,-(a*1.0e-3/2+(I-1)*(h+a)*1.0e-3)+0.5,-(a*1.0e-3/2+(I*h+(I-1)*a)*1.0e-3)+0.5,(Dxi*1.0e-3)/2+0.5,(Dx0*1.0e-3)/2+0.5
I=I+1
*ENDDO
RECTNG,-1.0e-3*68/2+0.5,1.0e-3*68/2+0.5,1.0e-3*(Dx0/2+5)+0.5,1.0e-3*(Dx0/2+20)+0.5
RECTNG,-1.0e-3*3*Dx0+0.5,1.0e-3*3*Dx0+0.5,0.5,-1.0e-3*3*Dxi+0.5
RECTNG,-1.0e-3*3*Dx0+0.5,1.0e-3*3*Dx0+0.5,0.5,1.0e-3*3*Dxi+0.5
AOVLAP,ALL
NUMCMP,AREA !压缩编号
APLOT
SAVE,emage-2D-geom.db
SAVE
local,11,0
ASEL,S,AREA,,1,
AATT,1,,1,11
ASEL,S,AREA,,2,
AATT,2,,2,11
ASEL,S,AREA,,4,18,2
AATT,3,,3,11
ASEL,S,AREA,,3,17,2
AATT,4,,4,11
ASEL,S,AREA,,19,
AATT,5,,5,11
ASEL,S,AREA,,20,21 !材料赋给单元
AATT,6,,6,11
allsel,below,area !选择所有实体
/pnum,line,1
lplot
LESIZE,2,,,3 !分割线
LESIZE,4,,,3
LESIZE,1,,,50
LESIZE,3,,,50
asel,s,,,1
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,1
asel,all
LESIZE,6,,,3 !分割线
LESIZE,8,,,3
LESIZE,5,,,50
LESIZE,7,,,50
asel,s,,,2
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,2
asel,all
LESIZE,69,,,5
LESIZE,71,,,5 !分割线
LESIZE,70,,,10
LESIZE,72,,,10
asel,s,,,18
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,18
asel,all
LESIZE,53,,,5 !分割线
LESIZE,54,,,10
LESIZE,55,,,5
LESIZE,56,,,10
asel,s,,,14
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,14
asel,all
LESIZE,37,,,5 !分割线
LESIZE,39,,,5
LESIZE,38,,,10
LESIZE,40,,,10
asel,s,,,10
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,10
asel,all
LESIZE,21,,,5 !分割线
LESIZE,23,,,5
LESIZE,22,,,10
LESIZE,24,,,10
asel,s,,,6
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,6
asel,all
LESIZE,13,,,5 !分割线
LESIZE,15,,,5
LESIZE,14,,,10
LESIZE,16,,,10
asel,s,,,4
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,4
asel,all
LESIZE,29,,,5 !分割线
LESIZE,31,,,5
LESIZE,30,,,10
LESIZE,32,,,10
asel,s,,,8
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,8
asel,all
LESIZE,45,,,5 !分割线
LESIZE,47,,,5
LESIZE,46,,,10
LESIZE,48,,,10
asel,s,,,12
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,12
asel,all
LESIZE,61,,,5 !分割线
LESIZE,63,,,5
LESIZE,62,,,10
LESIZE,64,,,10
asel,s,,,16
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,16
asel,all
LESIZE,65,,,5 !分割线
LESIZE,67,,,5
LESIZE,66,,,10
LESIZE,68,,,10
asel,s,,,17
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,17
asel,all
LESIZE,49,,,5 !分割线
LESIZE,51,,,5
LESIZE,50,,,10
LESIZE,52,,,10
asel,s,,,13
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,13
asel,all
LESIZE,33,,,5 !分割线
LESIZE,35,,,5
LESIZE,34,,,10
LESIZE,36,,,10
asel,s,,,9
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,9
asel,all
LESIZE,17,,,5 !分割线
LESIZE,19,,,5
LESIZE,18,,,10
LESIZE,20,,,10
asel,s,,,5
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,5
asel,all
LESIZE,9,,,5 !分割线
LESIZE,11,,,5
LESIZE,10,,,10
LESIZE,12,,,10
asel,s,,,3
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,3
asel,all
LESIZE,25,,,5 !分割线
LESIZE,27,,,5
LESIZE,26,,,10
LESIZE,28,,,10
asel,s,,,7
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,7
asel,all
LESIZE,41,,,5 !分割线
LESIZE,43,,,5
LESIZE,42,,,10
LESIZE,44,,,10
asel,s,,,11
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,11
asel,all
LESIZE,57,,,5 !分割线
LESIZE,59,,,5
LESIZE,58,,,10
LESIZE,60,,,10
asel,s,,,15
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,15
asel,all
LESIZE,73,,,30 !分割线
LESIZE,75,,,30
LESIZE,74,,,15
LESIZE,76,,,15
asel,s,,,19
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,1
AMESH,19
asel,all
!分割线
LESIZE,83,,,40
LESIZE,85,,,10
LESIZE,79,,,40
LESIZE,86,,,10
asel,s,,,20
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,0
AMESH,20
asel,all
LESIZE,80,,,10
LESIZE,79,,,40
LESIZE,78,,,10
LESIZE,77,,,40
asel,s,,,21
MSHAPE,0,2D !指定单元形状,四边形 2d
MSHKEY,0
AMESH,21
ALLSEL,ALL
esel,s,mat,,3,4,1
cm,coil,elem !define coil component
alls
ALLSEL,ALL
FLST,2,1,4,ORDE,1
FITEM,2,79
DL,P51X, ,ASYM
FINISH
/solu
S=0.005*0.01 !线圈截面积
u0=8000 !定义放电电压
cap=100e-6 !定义电容
R=10e-3 !定义电阻
L=1E-10 !定义电感
bet=R/(2*L)
pi=acos(-1)
omg=1/sqrt(L*cap) !频率1/(L*cap)^(1/2)
total_tim=2*pi/omg !周期(2*pi)*(L*cap)^(1/2)
tim_step=1e-6
*do,i,1,40
time,i*tim_step
j=((u0/(L*omg))*exp(-bet*i*tim_step)*sin(omg*i*tim_step))/(0.005*0.01)
antype,trans
kbc,0
CMSEL,S,coil
bF,coil,js,,,j,
NSUBST,1,1000,0
outres,all,all
allsel,all
solve
*enddo
save
FINISH
页:
[1]