齿轮动态分析是出了问题，求帮助

blues-et

哪个高人给点指点啊，dyna做的网格大了就是惯性张量太大，让该密度，而网格小了就是内存不足
fini
/clear
/prep7
et,1,shell163
et,2,shell163
et,3,solid164
et,4,solid164
keyopt,1,1,7
keyopt,2,1,7
r,1,,,0.1
r,2,,,0.1
mp,dens,3,7.8e3
mp,EX,3,2.1e11  
MP,NUXY,3,0.3
mp,dens,4,7.8e3
MP,EX,4,2.1e11  
MP,NUXY,4,0.3
edmp,rigi,1,7,4
edmp,rigi,2,7,4
MP,DENS,1,7800
mp,EX,1,2.1e11  
MP,NUXY,1,0.3
MP,DENS,2,7800
MP,EX,2,2.1e11  
MP,NUXY,2,0.3

asel,s,,,159,160
aatt,1,1,1
esize,,16
mshape,0,2d
amesh,all
alls
asel,s,,,295,296
aatt,2,2,2
esize,,14
mshape,0,2d
amesh,all
esize,,16
vatt,3,,3
vsweep,1
esize,14
vatt,4,,4
vsweep,2
edpart,create
PARTSEL,'PLOT'
alls
EDCGEN,asts,3,4,0.28,0.18,0,0,0,,,,,0,10000000,0,0
EDCONTACT,0.1,0,2,1,1,1,1,4,0
*dim,time,array,2,1,1
*set,time(1,1,1),0     
*set,time(2,1,1),1
*dim,angular_velocity,array,2,1,1
*set,angular_velocity(1,1,1),300      
*set,angular_velocity(2,1,1),300
*dim,torque,array,2,1,1   
*set,torque(1,1,1),56079   
*set,torque(2,1,1),56079
EDLOAD,ADD,RBOZ,0,1,TIME,angular_velocity, 0, , , , ,
EDLOAD,ADD,RBMZ,0,2,TIME,torque, 0, , , , ,
finish

/sol
time,0.5
edenergy,2,2,2,2
edrst,2000
edhtime,2000
eddump,1  
/status,solu
solve