模型采用简支梁上移动质量弹簧系统。
命令流如下:!移动质量弹簧系统的动力分析---位移耦合法
!定义结果步数为5000
FINISH
/CLEAR
/CONFIG,NRES,5000
/PREP7
SPANL=32
M1=5000
M2=6000
M3=14000
K1=1E7
K1C1=500
K1C2=70
K2=1E8
K2C1=200
K2C2=400
SPEEDV=120*1000/3600 !设置速度
EM=3.5E10
AREA=4.5
IM=5.5
DENG=2600
GRA=9.8
NE=100
NN=NE+1
NN1=NN+1
DL=SPANL/NE
DT=DL/SPEEDV
F1=ACOS(-1)/2/SPANL/SPANL*SQRT(EM*IM/(AREA*DENG))
ET,1,BEAM3
ET,2,MASS21,,,4
ET,3,COMBIN14,1,2
MP,EX,1,EM
MP,NUXY,1,0.2
MP,DENS,1,DENG
R,1,AREA,IM,3.5
R,2,M1
R,3,M2
R,4,M3
R,5,K1,K1C1,K1C2
R,6,K2,K2C1,K2C2
*DO,I,1,NN
N,I,(I-1)*DL
*ENDDO
TYPE,1
REAL,1
*DO,I,1,NE
E,I,I+1
*ENDDO
N,NN1
TYPE,2
REAL,2
E,NN1
N,NN1+1,,0.5
REAL,3
E,NN1+1
N,NN1+2,,1.0
REAL,4
E,NN1+2
TYPE,3
REAL,5
E,NN1,NN1+1
REAL,6
E,NN1+1,NN1+2
D,1,UX,,,,,,UY
D,NN,ALL
FINISH
/SOLU
ANTYPE,4
TRNOPT,FULL !选择瞬态响应分析类型,完全法
SSTIF,ON !是否开启应力强化---一般情况开启
TIMINT,OFF !是否打开瞬态效应(又称为时间积分效应)不打开则为静态。此处先进行静态一步,以确保瞬态分析初始时间不为0
TIME,1E-5
ACEL,,GRA
NSUBST,2 !设置子步数
KBC,1 !设置子荷载步荷载变化类型0为递增,1为阶跃方式
SOLVE
TIMINT,ON
OUTRES,ALL,ALL !设置输出类型
DELTIM,0.05*DT,,0.1*DT !在本荷载步指定时间步长大小 或同NSUBST,SOLCONTROL(确定是否使用已优化的非线性默认求解设置和某些强化的内部求解算法)可用来提供可靠有效的求解设置
KBC,1
AUTOTS,ON !指定是否使用自动时间步长跟踪或载荷步跟踪
DDELE,NN1,UY
*DO,I,1,NN1
CPDELE,ALL
CDIST=(I-1)*DL
D,NN1+1,UX,CDIST
D,NN1+2,UX,CDIST
D,NN1,UX,CDIST
CP,NEXT,UY,I,NN1
TIME,I*DT
SOLVE
*ENDDO
在计算过程中,提示有两个warning:
Nodes I and J of element 104 ( COMBIN14 ) are not coincident.
*** WARNING *** CP = 11.996 TIME= 23:21:17
Nodes I and J of element 105 ( COMBIN14 ) are not coincident.
哪位大神帮解决一下。。。。 |