请问在场路耦合的瞬态分析中如何处理电容器件的参数
如题，我所建立的一个模型中，主要目的是分析次级线圈在z方向上位置不同处的耦合效率，由于模型中需要二极管元器件，故只能采用瞬态分析时。可计算过程中不收敛，我觉得可能是电容参数的设置问题，故想请教各位。/title, analyse of magnet distribute of the ELV type electron accelerator
/FILNAME,ELV_Magnet_Distribute_2D,0
/prep7
et,1,53,,,1
et,2,53,3,,1
mp,murx,1,1
TB,BH,2,1,20, !magnet 30q/130 steel
TBPT,,2,0.05
TBPT,,3,0.06
TBPT,,5,0.1
TBPT,,7,0.17
TBPT,,10,0.34
TBPT,,20,0.88
TBPT,,30,1.4
TBPT,,50,1.55
TBPT,,70,1.6
TBPT,,100,1.67
TBPT,,200,1.72
TBPT,,300,1.75
TBPT,,500,1.805
TBPT,,700,1.82
TBPT,,1000,1.85
TBPT,,2000,1.87
TBPT,,3000,1.88
TBPT,,5000,1.9
TBPT,,7000,1.905
TBPT,,10000,1.93
mp,murx,3,1
mp,rsvx,3,1e-12
mp,murx,4,1
mp,rsvx,4,1e-12
!sencondary coil parameters
sc_l1=329
sc_l2=369
sc_h1=84.5
sc_h2=89.1
sc_h=33
sc_h_delta=43*0
turn_s=3000
fill_s=0.4
res_s=1600
!primary coil parameters
pc_l1=518
pc_l2=659
pc_l=12
pc_h1=70
pc_h2=1800
turn_p=30
fill_p=0.85
res_p=0.01
volt_p=630
!magnet circuit parameter
mp_l1=160
mp_l2=555
mp_l3=696
mp_l4=400
mp_l=20
mp_h1=65
mp_h2=1805
mp_h3=2188
!calculate regions
a_l=800
a_h=2500
!others
circuit_loc=1000
cap=2e-8
res_load=1.2e6
/PNUM ,AREA,1
rectng,0,mp_l1,0,mp_h1
rectng,mp_l1,mp_l2,0,mp_h1
rectng,mp_l2,mp_l2+mp_l,0,mp_h1
K,13,mp_l3,mp_h2,,
K,14,mp_l3+mp_l,mp_h2,,
K,15,mp_l4,mp_h3,,
K,16,mp_l4+mp_l,mp_h3,,
K,17,pc_l1,pc_h1,,
K,18,pc_l1+mp_l,pc_h1,,
K,19,pc_l2,pc_h2,,
K,20,pc_l2+mp_l,pc_h2,,
K,21,sc_l1,sc_h1+sc_h_delta,,
K,22,sc_l1,sc_h1+sc_h+sc_h_delta,,
K,23,sc_l2,sc_h2+sc_h+sc_h_delta,,
K,24,sc_l2,sc_h2+sc_h_delta,,
!cone magnet
A,7,11,14,13
!cone magnet of the top
A,13,14,16,15
!primary coil
A,17,18,20,19
!secondary coils
A,21,22,23,24
rectng,0,a_l,0,a_h
APLOT
aovlap,all
NUMCMP,AREA
APLOT
asel,s,area,,4 ! primary
aatt,3,1,2
asel,s,area,,5 ! secondary
aatt,4,2,2
asel,s,area,,2,3 ! iron
asel,a,area,,6,7
aatt,2,4,1
asel,s,area,,1 ! window
asel,a,area,,8
aatt,1,4,1
asel,s,area,,1,7
esize,15
amesh,all
allsel,all
esize,75
amesh,8
asel,s,area,,4
asum
*get,area_p,area,,area
area_p=area_p/1e6
asel,s,area,,5
asum
*get,area_s,area,,area
area_s=area_s/1e6
ALLSEL,ALL
arscal,all,0,0,0.001,0.001,,,,1
r,1,area_p,turn_p,pc_h2-pc_h1,1,fill_p
r,2,area_s,turn_s,sc_h,1,fill_s
esel,s,mat,,3
eplot
nsle
nplot
*Get,N_pnum,Node,0,Num,Max
esel,s,mat,,4
eplot
nsle
nplot
*Get,N_snum,Node,0,Num,Max
allsel,all
eplot
T=1E-8
PI=2*ASIN(1)
f=400
um=volt_p*1.414
u=um*SIN(2*PI*f*T)
con=2.5e-3/4
c=0
!primary
R,4,u,,
N,,circuit_loc/1e3,0.04,0
*Get,N_cnum,Node,0,Num,Max
N,,circuit_loc/1e3,0.12,0
N,,circuit_loc/1e3+0.03,0.08,0
RMOD,4,15,0.3E-01,1
ET,3,CIRCU124,4,0
TYPE,3
REAL,4
MAT,1
E,n_cnum+1,n_cnum,n_cnum+2
R,5,1,
N,,circuit_loc/1e3,0.08,0
ET,4,CIRCU124,5,0
TYPE,4
REAL,5
MAT,1
E,n_cnum+1,n_cnum+3,n_pnum
R,6,res_p,
ET,5,CIRCU124,0,0
TYPE,5
REAL,6
MAT,1
E,n_cnum+3,n_cnum
!secondary
N,,circuit_loc/1e3+0.15,0.5
N,,circuit_loc/1e3+0.15,0.54
N,,circuit_loc/1e3+0.15,0.58
N,,circuit_loc/1e3+0.075,0.58
N,,circuit_loc/1e3,0.54
ET,6,CIRCU124,0,0 !resistor
R,7,res_load,
RMOD,7,15,-0.3E-01,1
TYPE,6
REAL,7
MAT,1
E,n_cnum+4,n_cnum+6
ET,7,CIRCU124,2,0 !capacitance
R,8,cap
TYPE,7
REAL,8
MAT,1
E,n_cnum+4,n_cnum+5
ET,8,CIRCU124,2,0
R,9,cap
TYPE,8
REAL,9
MAT,1
E,n_cnum+5,n_cnum+6
ET,11,CIRCU125 !diode
R,12
TYPE,11
REAL,12
MAT,1
E,n_cnum+4,n_cnum+7
ET,12,CIRCU125
R,13
TYPE,12
REAL,13
MAT,1
E,n_cnum+7,n_cnum+6
ET,15,CIRCU124,0,0 !resistor
R,16,res_s,
TYPE,15
REAL,16
MAT,1
E,n_cnum+5,n_cnum+8
ET,17,CIRCU124,5,0
R,18,1,
TYPE,17
REAL,18
MAT,1
E,n_cnum+7,n_cnum+8,n_snum
NSEL,S,NODE,,n_cnum
NSEL,A,NODE,,n_cnum+4
d,all,volt
esel,,ename,,53
esel,r,real,,1
nsle
cp,1,curr,all
cp,2,emf,all
esel,,ename,,53
esel,r,real,,2
nsle
cp,3,curr,all
cp,4,emf,all
alls
esel,,ename,,53
nsle
nsel,r,ext
d,all,az
alls
FINISH
/sol
ANTYPE,4
TRNOPT,FULL
*CREATE,ELV
time,t
kbc,0 !ramp
nsubst,5 !substeps
neqit,10
c=c+1
outres,all,last
OUTPR,ALL,ALL
allsel,all
solve
t=t+con
u=um*SIN(2*PI*f*T)
r,4,u
R,8,cap,-0.866510
R,9,cap,-0.866510
*END
*do,c,1,13
*USE,ELV
*enddo
FINISH
次级耦合电路如图：
[[i] 本帖最后由 magnett 于 2008-7-24 21:56 编辑 [/i]]
页:
[1]