/filname,t4,1
/title,the analysis of transient load
/prep7
!define elements and meterial
et,1,117
mp,murx,1,1 !air
mp,murx,2,1 !gap
mp,murx,3,1 !coil
mp,rsvx,3,1.7e-8
mp,rsvx,4,60e-4 !iron
TB,BH,4,,40
TBPT,,355,0.7
,,405,0.8
,,470,0.9
,,555,1.0
,,673,1.1
,,836,1.2
,,1065,1.3
,,1220,1.35
,,1420,1.4
,,1720,1.45
,,2130,1.5
,,2670,1.55
,,3480,1.6
,,4500,1.65
,,5950,1.7
,,7650,1.75
,,10100,1.8
,,13000,1.85
,,15900,1.9
,,21100,1.95
,,26300,2.0
,,32900,2.05
,,42700,2.1
,,61700,2.15
,,84300,2.2
,,110000,2.25
,,135000,2.3
,,200000,2.41
,,400000,2.69
,,800000,3.22
!define parameters
ix=0.045
iy=0.1524
iz=0.127
fy1=0.0254
fy2=0.0508
gz=0.0254
fz=0.0762
gap=0.005
bone=0.0015 !骨架
cw=0.0132
ch=0.024
i=283 !200A的峰值
n=36
jdens1=i*n/(ch*cw)
!create frame
/pnum,volu,1
block,-ix/2,ix/2,-iy/2,iy/2,-iz/2,iz/2-gz
block,-ix/2,ix/2,-iy/2+fy1,-fy2/2,-iz/2+gz,iz/2
vsymm,y,2
vsbv,1,2
vsbv,4,3
numcmp,all
block,-ix/2,ix/2,-fy2/2,fy2/2,fz/2-gap,fz/2
vovlap,all
numcmp,all
block,-ix/2,ix/2,-iy/2,iy/2,iz/2-gz,iz/2
vglue,all
numcmp,all
block,ix/2,ix/2+bone,-fy2/2,fy2/2,-fz/4-ch/2,-fz/4+ch/2
vsymm,x,4
block,-ix/2,ix/2,fy2/2,fy2/2+bone,-fz/4-ch/2,-fz/4+ch/2
vsymm,y,6
vsymm,z,4,7
local,11,1,ix/2,fy2/2,-fz/4+ch/2
wpcsys,-1,11
cyl4,0,0,0,0,bone,90,-ch
local,12,1,-ix/2,fy2/2,-fz/4+ch/2
wpcsys,-1,12
cyl4,0,0,0,90,bone,180,-ch
local,13,1,-ix/2,-fy2/2,-fz/4+ch/2
wpcsys,-1,13
cyl4,0,0,0,180,bone,270,-ch
local,14,1,ix/2,-fy2/2,-fz/4+ch/2
wpcsys,-1,14
cyl4,0,0,0,-90,bone,0,-ch
local,15,1,ix/2,fy2/2,fz/4-ch/2
wpcsys,-1,15
cyl4,0,0,0,0,bone,90,ch
local,16,1,-ix/2,fy2/2,fz/4-ch/2
wpcsys,-1,16
cyl4,0,0,0,90,bone,180,ch
local,17,1,-ix/2,-fy2/2,fz/4-ch/2
wpcsys,-1,17
cyl4,0,0,0,180,bone,270,ch
local,18,1,ix/2,-fy2/2,fz/4-ch/2
wpcsys,-1,18
cyl4,0,0,0,-90,bone,0,ch
vglue,all
numcmp,all
wpcsys,-1,0
csys,0
block,ix/2+bone,ix/2+bone+cw,-fy2/2,fy2/2,-fz/4-ch/2,-fz/4+ch/2
vsymm,x,20
block,-ix/2,ix/2,fy2/2+bone,fy2/2+bone+cw,-fz/4-ch/2,-fz/4+ch/2
vsymm,y,22
vsymm,z,20,23
wpcsys,-1,11
cyl4,0,0,bone,0,bone+cw,90,-ch
wpcsys,-1,12
cyl4,0,0,bone,90,bone+cw,180,-ch
wpcsys,-1,13
cyl4,0,0,bone,180,bone+cw,270,-ch
wpcsys,-1,14
cyl4,0,0,bone,-90,bone+cw,0,-ch
wpcsys,-1,15
cyl4,0,0,bone,0,bone+cw,90,ch
wpcsys,-1,16
cyl4,0,0,bone,90,bone+cw,180,ch
wpcsys,-1,17
cyl4,0,0,bone,180,bone+cw,270,ch
wpcsys,-1,18
cyl4,0,0,bone,-90,bone+cw,0,ch
vglue,all
numcmp,all
wpcsys,-1,0
block,-ix*2,ix*2,-iy*2,iy*2,-iz*2,iz*2
vovlal,all
numcmp,all
!mesh
vsel,s,volu,,1
vatt,2,,1,0
vsel,s,volu,,2,19,17
vatt,4,,1,0
vsel,s,volu,,3,18
vatt,1,,1,0
vsel,s,,,20,27
vatt,3,,1,0
vsel,s,volu,,29
vatt,3,,1,11
vsel,s,volu,,33
vatt,3,,1,15
vsel,s,volu,,31
vatt,3,,1,12
vsel,s,volu,,35
vatt,3,,1,16
vsel,s,volu,,30
vatt,3,,1,13
vsel,s,volu,,34
vatt,3,,1,17
vsel,s,volu,,28
vatt,3,,1,14
vsel,s,volu,,32
vatt,3,,1,18
vsel,s,volu,,36
vatt,1,,1,0
vsel,s,volu,,3,18
vsel,a,volu,,20,35
vmesh,all
vsel,s,volu,,1,2
vsel,a,volu,,19
mshkey,0
mshape,1,3d
esize,0.02
vmesh,all
vsel,s,volu,,36
esize,0.05
vmesh,all
*dim,coils,char,2
*dim,cur,array,2
coils(1)='coil1'
cur(1)=i
coils(2)='coil2'
cur(2)=i
vsel,s,volu,,20,23
vsel,a,volu,,28,31
eslv,s
cm,coils(1),elem
vsel,s,volu,,24,27
vsel,a,volu,,32,35
eslv,s
cm,coils(2),elem
alls
finish
/solu
*cfopen,pt,txt,F:\work\command\chock\test\
antype,trans
nsel,ext
d,all,az,0
t1=0.05
t2=0.005
t=1e-6
c=0
nn=40
pi=2*asin(1)
con=2*pi/nn
neqit,1
*create load
jdens=jdens1*sin(t)
*vwrite,jdens
('jdens=',jdens,)
vsel,s,volu,,20,24,4
eslv,s
bfe,all,js,,0,jdens
vsel,s,volu,,22,26,4
eslv,s
bfe,all,js,,-jdens
vsel,s,volu,,21,25,4
eslv,s
bfe,all,js,,0,-jdens
vsel,s,volu,,23,27,4
eslv,s
bfe,all,js,,jdens
vsel,s,volu,,28,35
eslv,s
bfe,all,js,,0,jdens
t=t+con
c=c+1
outres,all,1
*if,c,eq,(n*0.75+1),then
outpr,,1
*elseif,c,eq,(n+1),then
outpr,,1
*else
outpr,,0
*endif
alls
time,t1
deltim,t2
solve
*end
*do,j,1,41,1
*use,load
*enddo
*cfclos
finish |