刚才忘了将求解命令流附上,希望楼主别生气,谢谢
finish !返回上一次主菜单
/solu !进入求解器
csys,0 !改变坐标系
/view,1,-1,0.5,1 !改变视图显示方向
/replot !重新绘制
eplot !显示单元图
!save !保存数据库
nsubst,10 !指定荷载子步
nlgeom,on !指定大变形
nropt,full !指定newton-lapnace迭代
pred,on !启用非线性求解预测器
lnsrch,on !启用线性搜索
outres,all,last !选择将每一子步的最后一个迭代计算步结果保存
!施加约束条件
nsel,s,loc,x,-50 !选择节点
d,all,ux !施加X方向约束
nsel,s,loc,x,25 !选择节点
d,all,ux !施加X方向约束
nsel,s,loc,y,-40 !选择节点
d,all,all !uy !施加Y方向约束
nsel,s,loc,z,15 !选择节点
nsel,a,loc,z,-24.6
d,all,uz !施加Z方向约束
allsel,all
acel,0,10,0 !施加重力
step=0 !用于记录计算步的参数
time,1 !第1个荷载步(时间步)
solve !求解
save !存储数据库
finish !返回上一次主菜单
/solu
antype,,rest !设定分析类型为重启动分析,重启动点默认为上次求解结束时
time,2 !设定时间步
AUTOTS,1 !自动增加计算步
DELTIM,0.1,0.001,0.2,1 !**时间步选项设置
csys,12 !将当前坐标变为前面自定义的极座标
wpcsys,-1
vsel,s,loc,x,0,6 !选择体
vsel,r,loc,z,0,-3.6 !再选择体
eslv,s !选择单元
ekill,all !杀死开挖土体单元,包括核心图、管片层、注浆层
asel,s,loc,x,0,5.99 !选择面
asel,r,loc,z,-3.6
SFA,all,1,PRES,0.3e6 !施加掌子面顶进压力
asel,s,loc,x,6 !选择面
asel,r,loc,z,0,-3.59
SFA,all,1,PRES,0.15e6 !施加注浆压力
allsel,all !选择所有元素
solve !求解
save !存储数据库
finish !返回上一主菜单
!****第1步开挖,第2个计算步******
/solu !进入求解器
antype,,rest !重新启动
time,3 !第三次计算
AUTOTS,1
DELTIM,0.1,0.001,0.2,1
csys,12 !改变坐标系
wpcsys,-1 !改变工作坐标系
vsel,s,loc,x,0,5.2 !选择体
vsel,r,loc,z,0,-3.6
eslv,s !选择单元
ekill,all !杀死核心土单元
vsel,s,loc,x,5.251,5.79 !选择体
vsel,r,loc,z,0,-3.6
eslv,s !选择单元
mpchg,5,all !改变管片单元材料参数
vsel,s,loc,x,5.81,5.99 !选择体
vsel,r,loc,z,0,-3.6
eslv,s !选择单元
mpchg,6,all !改变注浆层材料参数
asel,s,loc,x,0,5.99 !选择面
asel,r,loc,z,-3.6
SFA,all,1,PRES,0.3e6 !在面上施加压力
asel,s,loc,x,6
asel,r,loc,z,0,-3.6
SFA,all,1,PRES,0 !取消注浆压力
allsel,all !选择所有元素
solve !求解
save |