帮朋友问的一个算例

zoe122208

朋友要做一个线圈内部嵌套2个不同材质的管件，线圈放电瞬间产生的力使两个管件连接，命令流如下，请大家帮忙看看有何问题：FINISH
/CLEAR
/FILNAM,3 dimention welding-1223fuck
/TITLE,3 dimention welding1-1223fuck
KEYW,MAGEDG,1


/PREP7
ET,1,SOLID97,1     !铜线圈单元
ET,2,SOLID97,1     !管和边界单元


MP,MURX,1,1          !铜线圈的相对磁导率
MP,RSVX,1,1.72E-8    !铜线圈的电阻率
MP,MURX,2,1          !铝管的相对磁导率
MP,RSVX,2,2.83E-8    !铝管的电阻率
MP,MURX,3,1          !铜管的相对磁导率
MP,RSVX,3,1.72E-8    !铜管的电阻率
MP,MURX,4,1          !空气的相对磁导率

KEYOPT,1,1,2
KEYOPT,1,2,0
KEYOPT,1,5,0

N=5
DZ=22
T=6

CSYS,1
K,1,DZ/2,0,-T/2
K,2,DZ/2,180
L,1,2
CSYS,0
LSYMM,Z,1
LSYMM,Y,2,,,,,1
NUMMRG,ALL
CM,L1,LINE

KWPAVE,1
WPROTA,,90
RECTNG,0,7,0,5,

VDRAG,1,,,,,,L1
VGEN,N,ALL,,,,,T
NUMMRG,KP
WPCSYS

VADD,1,2,3,4,5,6,7,8,9,10
NUMCMP,ALL
CYLIND,10,9,0,100,0,360,
CYLIND,8,6,32,-68,0,360,
CYLIND,30,0,-100,150,0,360,

VOVLAP,ALL
ALLSEL,ALL
NUMCMP,ALL



ALLSEL,ALL
VSEL,S,,,1
VATT,1,1,1
VSEL,S,,,2
VATT,2,,2
VSEL,S,,,3
VATT,3,,2
VSEL,S,,,4
VATT,4,,2


VSEL,S,VOLU,,1
SMRT,8               
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL

VSEL,S,VOLU,,2
SMRT,8               
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL

VSEL,S,VOLU,,3
SMRT,8               
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL

VSEL,S,VOLU,,4
SMRT,8               
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL

SAVE

ALLSEL,ALL

ASEL,S,,,55,58,1
CM,a_air,area

cmsel,s,a_air    !选空气外表面
extopt,esize,1    !只生成一层单元
type,2
mat,2  !infin111单元,材料为空气
vext,all,,,,,,2,2,2 !延伸出一层单元

asel,s,,,59,72,1   !选单元的外表面
nsla,s,1    !选外表面的节点
sf,all,inf   !定义无限远

ALLSEL,ALL

ASEL,S,AREA,,1
NSLA,S,1
cp,1,volt,all
ALLSEL,ALL

ASEL,S,AREA,,1
NSLA,S,1
CM,COIL,NODE

ASEL,S,AREA,,42
NSLA,S,1
D,ALL,volt,0
ALLSEL,ALL



LOCAL,20,1,0,0,0,
ESYS,20
EMODIF,ALL
R,1,35,5,12089,0,1,0,

CSYS,2
NSEL,S,LOC,X,100     
SF,ALL,INF

CSYS,0
NSEL,S,LOC,X,0      
D,ALL,AZ,0
ALLSEL,ALL

u0=15000               
cap=100e-6               
R=45.91e-3                 
L=37E-07      
bet=R/(2*L)   
pi=acos(-1)
omg=1/sqrt(L*cap)
height=0.007  
width=0.005  
hc=(N*width-0.001)            
total_tim=2*pi/omg        
tinc=2e-6              
ftime=0.5*total_tim        
startime=1e-12            



/SOLU

ANTYPE,trans

*do,tme,startime,ftime,tinc
KBC,0
TRNOPT,FULL
time,tme
nsubst,2,2,1  
OUTRES,ALL,ALL

j=N*(u0/(L*omg)*exp(-bet*tme)*sin(omg*tme))

CMSEL,S,COIL
F,COIL,AMPS,j
AUTOTS,1

allsel,all
solve
*enddo
save
allsel


FINISH

zoe122208

谁帮帮忙看看哪里有问题

绿色的农民

希望有人帮忙的

admint

额 顶下 希望高手看到