马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
算了9千多个子步才出来结果,这是不是就是所谓的收敛慢啊?
在此先谢过了!
题目是一个水平滑块相对于另一个水平滑块滑动(以平面来代替体,也即平面应力问题),不施加压力,那么摩擦力应该其实不起作用。没有重力,滑块两端速度不同,应该类似于弹簧的收缩,本来静态就可以解决的,但我想看下瞬态解决什么样子,当滑块两端速度VL=VR时,收敛速度非常快,26个子步就完成了,但是只是将VR减小了一点,就变成了9千多步,请问有什么办法解决吗?网格划分是否有问题?或者是不是可以增加一个时间非常短的载荷步,速度会快一点?
!!水平滑移,无压力,无重力VL>VR,竟然花去了9千多个子步,剪应力第一子步结果有部分较大(4.0),其余都比较小,处于(e-5)量级,下块体无滑动,上块体y向位移打到了0.2,x向位移呈线性递增,剪应力接近于零,y方向应力接近于0,x向应力随时间成正比。
FINISH
/CLEAR
/CONFIG,NRES,10000
L=100E3
H=10E3
T0=3600*24*365
VL=0.05/T0
VR=0.02/T0
/PREP7
ET,1,PLANE82
ET,2,TARGE169
ET,3,CONTA172
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
KEYOPT,2,2,0
KEYOPT,2,3,0
KEYOPT,2,4,11
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
R,1,,,, ,,
RMORE,,,, ,,
RMORE,,,, ,,
RMORE,,,5, 1E15
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!11
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,2.67E3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.3E10
MPDATA,PRXY,1,,0.28
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MU,1,,0.3
RECTNG,,L,,H
RECTNG,,L/2,H,2*H
!LSEL,S,TAN1,X
!LESIZE,ALL,0.15*H
!LSEL,S,LOC,Y,0
!LSEL,A,LOC,Y,2*H
!LESIZE,ALL,0.1*L
!LSEL,S,LOC,Y,H
!LESIZE,ALL,0.05*L
MSHAPE,0
MSHKEY,1
AMESH,ALL
LSEL,S,LOC,X,0.5*L
LSEL,R,LOC,Y,H
NSLL,S,1
TYPE,2
ESURF
LSEL,S,LOC,X,0.25*L
LSEL,R,LOC,Y,H
NSLL,S,1
TYPE,3
REAL,1
ESURF
ALLSEL
FINISH
/SOLU
ANTYPE,TRANS
NLGEOM,ON
DELTIM,T0,0.01*T0,10*T0
OUTRES,ALL,ALL
TIME,100*T0
AUTOTS,ON
KBC,0
LSEL,S,LOC,Y,0
NSLL,s,1
D,ALL,ALL,0
/PBC,U,,2
LSEL,S,LOC,X,0
LSEL,R,LOC,Y,H+0.5*H
NSLL,S,1
D,ALL,UX,VL*100*T0
LSEL,S,LOC,X,0.5*L
LSEL,R,LOC,Y,H+0.5*H
NSLL,S,1
D,ALL,UX,VR*100*T0
ALLSEL
ALLSEL
SOLVE
!/POST1
!PLNSOL,S,XY
!PLNSOL,U,X
!ANTIME,10,0.2,,,2 |