马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
求神解答:本人在做一个车桥耦合系统(移动荷载)的模拟,模拟简支梁上移动质量块-弹簧的作用。求解时出现如下问题,怎么回事呢?下面出现的问题和本人的命令流。
error:Unable to open file file.r001. Unit 29 is already in use. Please send your date to your ansys technical Support Provider
命令流
/config,nres,30000
/PREP7
!model creation
l=20.0 !beam length
h=0.2 !higth
b=0.2 !wigth
EX0=200E9
PRXY0=0.3
DENS0=7850
IM=3*0.2*0.2*0.2*0.2/12
N0=400 !elements number
Nh=10 !number of ele along depth direction
v=2 !velocity is v m/s
F=1000 !load amplitude
l_e=l/N0 !length of each element
time_ele=(l/v)/N0 !time of each element
ET,1,BEAM3$R,1,0.04,IM,0.2 !difine the shell element
!material setting
MPTEMP,1,0
MPDATA,EX,1,,EX0
MPDATA,PRXY,1,,PRXY0
MPDATA,DENS,1,,DENS0
ET,2,mass21,,,4 !difine mass elenent
R,2,F/9.8 !difine mass real constant
KSTIF=1e9$ET,3,14,,,2$R,3,KSTIF !difine the stiffness element
!creat nodes
*DO, i, 1, N0+1
N,i,(i-1)*l_e
*ENDDO
!creat elements
TYPE,1$REAL,1$mat,1
*DO, i, 1, N0
EN,i, i, i+1 !creat beam element
*ENDDO
TYPE,2$REAL,2 !set the MASS element type 2
*DO, i, 1, N0+1
N,1000+i,(i-1)*l_e,0.05 !creat the mass nodes
EN,1000+i,1000+i !creat the mass elements
*ENDDO
TYPE,3$REAL,3 !set the stiffness type 3
*DO, i, 1, N0+1
EN,2000+i,i,1000+i !creat the stiffness elements
*ENDDO
!creat damage
!trans analysis type
!displacement seting
D,1, ,0, , , ,UX,UY, , , ,
D,N0+1, ,0, , , ,UX,UY, , , ,
finish
/SOLU$ANTYPE,4$NLGEOM,ON !open large deformation
NROPT,FULL$ACEL,,9.8 !Newton-Raphson full transient analysis,open gravity
ESEL,S,REAL,,2$EKILL,ALL$ALLSEl,ALL !KILL all the mass element
timint,off$time,1E-6 !close all the transient setting
NSUBST,2$KBC,1$solv
timint,on,all$OUTRES,all,all !open trans
!alphad,0.001 !a damp
!betad,0.0001 !beta damp
EALIVE,10001 !Activation the mass No.10001
time,time_ele
NSUBST, 5
kbc,1$AUTOTS,ON
SOLVE
*do,i,1,N0,1
EKILL,10000+i$EALIVE,10001+i
time,time_ele*(i+1) !time cordinate
NSUBST, 5 !sample frequence
kbc,1$AUTOTS,ON
solve
*enddo
finish
!result extract
/POST26
*DIM,dis,array,5000,2
*DIM,velo,array,5000,2
*DIM,accel,array,5000,2
NUMVAR,200
!store 2Nos results
*do,i,1,2,1
j=i+1
N=(N0+1)*Nh+i*N0*1/4
NSOL,j,N,U,Y !extract displacement signa,j is the name of viration, must j>2
NSOL,j+10,N,ACC,Y !extract acceleration signal
NSOL,j+30,N,VEL,Y
*enddo
!
!
!read the results
*do,j,2,3,1
VGET,dis(1,j-1),j
VGET,accel(1,j-1),j+10
VGET,velo(1,j-1),j+30
*enddo
!
*CFOPEN,1023dis0,TXT
*Vwrite,dis(1,1),dis(1,2)
(2F13.8,' ')
*CFCLOS
finish
*CFOPEN,1023acce0,TXT
*Vwrite,accel(1,1),accel(1,2)
(2F13.6,' ')
*CFCLOS
finish
*CFOPEN,1023velo0,TXT
*Vwrite,velo(1,1),velo(1,2)
(2F16.11,' ')
*CFCLOS
finish
|