盾构穿越地铁ansys命令流.docx
《盾构穿越地铁ansys命令流.docx》由会员分享，可在线阅读，更多相关《盾构穿越地铁ansys命令流.docx（34页珍藏版）》请在冰点文库上搜索。
盾构穿越地铁ansys命令流
本着三个代表思想，做了此注解。
fini
/cle
*set,r1,2.9！
设参数内径
*set,r2,3.0！
外径
*set,y1,-36.25！
y方向尺寸
*set,y2,-21.25
*set,y3,-17.5
*set,y4,-13.75
*set,y5,-10
*set,y6,-6.25
*set,y7,1.0
*set,y8,3.5
*set,x1,0
*set,x2,15
*set,x3,17.5
*set,x4,32.5
*set,x5,40
*set,x6,42.5
*set,x7,57.5
*set,z1,0
*set,z2,-15
*set,z3,-17.5
*set,z4,-32.5
*set,z5,-40
*set,z6,-42.5
*set,z7,-57.5
/prep7！
设单元类型
et,1,mesh200,6
et,2,solid45
et,3,shell63
r,1,0.35,0.35,0.35,0.35！
常数
mp,ex,1,3.83e6！
材料类型
mp,prxy,1,0.2
mp,dens,1,2000
mp,ex,2,18.83e6
mp,prxy,2,0.28
mp,dens,2,1920
mp,ex,3,20.36e6
mp,prxy,3,0.32
mp,dens,3,1820
mp,ex,4,17.94e6
mp,prxy,4,0.35
mp,dens,4,1840
mp,ex,5,20.36e6
mp,prxy,5,0.32
mp,dens,5,1820.6
mp,ex,6,20.36e6
mp,prxy,6,0.32
mp,dens,6,1820.5
mp,ex,7,35.36e6
mp,prxy,7,0.25
mp,dens,7,1940
mp,ex,8,35.36e6
mp,prxy,8,0.25
mp,dens,8,1940.1
mp,ex,9,20.36e6
mp,prxy,9,0.32
mp,dens,9,1820.2
mp,ex,10,20.36e6
mp,prxy,10,0.32
mp,dens,10,1820.3
mp,ex,11,20.36e6
mp,prxy,11,0.32
mp,dens,11,1820.4
mp,ex,12,3.45e10
mp,prxy,12,0.17
mp,dens,12,1910
mp,ex,13,3.45e10
mp,prxy,13,0.17
mp,dens,13,1910.1
mp,ex,14,3.45e10
mp,prxy,14,0.17
mp,dens,14,1910.2
mp,ex,15,20.36e6
mp,prxy,15,0.28
mp,dens,15,1820.7
mp,ex,16,20.37e6
mp,prxy,16,0.28
mp,dens,16,1820.8
mp,ex,17,20.38e6
mp,prxy,17,0.28
mp,dens,17,1820.9
mp,ex,18,18.38e7
mp,prxy,18,0.24
mp,dens,18,1850
mp,ex,19,18.38e7
mp,prxy,19,0.24
mp,dens,19,1850.1
mp,ex,20,18.38e7
mp,prxy,20,0.24
mp,dens,20,1850.2
mp,ex,21,400e6
mp,prxy,21,0.24
mp,dens,21,1950.3
mp,ex,22,600e6
mp,prxy,22,0.24
mp,dens,22,1950.4
mp,ex,23,400e6
mp,prxy,23,0.24
mp,dens,23,1950.5
mp,ex,24,400e6
mp,prxy,24,0.24
mp,dens,24,1950.6！
材料2~4土体5~6在建隧道周围加固土体未改性质7~8改过性质9~11隧道内土体，9为已建1号线隧道土体10~11为下侧7号线隧道内土体，12~14为衬砌单元材料，12为一号线衬砌单元，13、14对应左右已建隧道，15~17注浆材料未改时，15对应一号线，本命令流中未用，16~17为下侧7号线注浆材料未改。
18~20对应15~17改后材料,所谓未改指性质和土体一样，参数小数点后有变化,21出洞旋喷加固材料22出洞搅拌加固23为7号线下旋喷土体预加固
blc4,x1,y1,57.5,y8-y1!
建立平面
blc4,x1,y1,57.5,y7-y1
blc4,x1,y1,57.5,y6-y1
blc4,x1,y1,57.5,y4-y1
blc4,x1,y1,57.5,y2-y1
blc4,x1,y1,x2,y8-y1
blc4,x1,y1,x6,y8-y1
blc4,x3,y2,7.5,7.5
blc4,x3,y2,3.75,7.5
blc4,x3,y2,7.5,3.75
blc4,x4,y2,7.5,7.5
blc4,x4,y2,3.75,7.5
blc4,x4,y2,7.5,3.75
Cyl4,x3+3.75,y3,r1！
建圆
Cyl4,x3+3.75,y3,r2
Cyl4,x4+3.75,y3,r1
Cyl4,x4+3.75,y3,r2
asel,all！
overlap
aovlap,all
nummrg,all,,,,low！
压缩
numcmp,all
k,59,x1,y1,z2！
建立点为拉伸用
k,60,x1,y1,z6
k,61,x1,y1,z7
l,1,59!
建拉伸线
l,59,60
l,60,61
手工划分///////！
mesh也可用命令流
type,1！
划分面
real,1
mat,10
Asel,s,area,,1,4,1
Amesh,all
type,1
real,1
mat,11
Asel,s,area,,5,8,1
Amesh,all
LPLOT
APLOT
ALLSEL,ALL
APLOT
type,1
real,1
mat,16
Asel,s,area,,9,12,1
Amesh,all
type,1
real,1
mat,17
Asel,s,area,,13,16,1
Amesh,all
type,1
real,1
mat,5
Asel,s,area,,29,32,1
Amesh,all
type,1
real,1
mat,6
Asel,s,area,,33,36,1
Amesh,all
type,1
real,1
mat,2
Asel,s,area,,19,20,1
Asel,a,area,,27
Amesh,all
type,1
real,1
mat,3
Asel,s,area,,21,22,1
Asel,a,area,,28
Amesh,all
type,1
real,1
mat,3
Asel,s,area,,23,24,1
Asel,a,area,,41
Amesh,all
type,1
real,1
mat,3
Asel,s,area,,37,39,1
Asel,a,area,,25,26,1
Amesh,all
type,1
real,1
mat,4
Asel,s,area,,17,18,1
Asel,a,area,,40
Amesh,all
lsel,s,LCCA！
删除连接线
LDELE,all
EXTOPT,ESIZE,5,0！
面拉伸
EXTOPE,ACLEAR,1
type,2
real,1
mat,10
Asel,s,area,,1,4,1
vDRAG,ALL,,,,,,99
EXTOPT,ESIZE,5,0
EXTOPE,ACLEAR,1
type,2
real,1
mat,11
Asel,s,area,,5,8,1
vDRAG,ALL,,,,,,99
EXTOPT,ESIZE,5,0
EXTOPE,ACLEAR,1
type,2
real,1
mat,16
Asel,s,area,,9,12,1
vDRAG,ALL,,,,,,99
EXTOPT,ESIZE,5,0
EXTOPE,ACLEAR,1
type,2
real,1
mat,17
Asel,s,area,,13,16,1
vDRAG,ALL,,,,,,99
EXTOPT,ESIZE,5,0
EXTOPE,ACLEAR,1
type,2
real,1
mat,5
Asel,s,area,,29,32,1
vDRAG,ALL,,,,,,99
EXTOPT,ESIZE,5,0
EXTOPE,ACLEAR,1
type,2
real,1
mat,6
Asel,s,area,,33,36,1
vDRAG,ALL,,,,,,99
EXTOPT,ESIZE,5,0
EXTOPE,ACLEAR,1
type,2
real,1
mat,2
Asel,s,area,,19,20,1
Asel,a,area,,27
vDRAG,ALL,,,,,,99
EXTOPT,ESIZE,5,0
EXTOPE,ACLEAR,1
type,2
real,1
mat,3
Asel,s,area,,21,26,1
Asel,a,area,,28
Asel,a,area,,37,39,1
Asel,a,area,,41
vDRAG,ALL,,,,,,99
EXTOPT,ESIZE,5,0
EXTOPE,ACLEAR,1
type,2
real,1
mat,4
Asel,s,area,,17,18,1
Asel,a,area,,40
vDRAG,ALL,,,,,,99
nummrg,all,,,,low
numcmp,all
EXTOPT,ESIZE,22,0
EXTOPE,ACLEAR,1
type,2
real,1
mat,2
Asel,s,area,,142
Asel,a,area,,147
Asel,a,area,,150
vDRAG,ALL,,,,,,100
type,2
real,1
mat,3
Asel,s,area,,155
Asel,a,area,,160
Asel,a,area,,179
Asel,a,area,,172
Asel,a,area,,176
Asel,a,area,,184
Asel,a,area,,189
Asel,a,area,,196
vDRAG,ALL,,,,,,100
type,2
real,1
mat,4
Asel,s,area,,206
Asel,a,area,,211
Asel,a,area,,220
vDRAG,ALL,,,,,,100
type,2
real,1
mat,10
Asel,s,area,,45
Asel,a,area,,48
Asel,a,area,,51
Asel,a,area,,53
vDRAG,ALL,,,,,,100
type,2
real,1
mat,11
Asel,s,area,,57
Asel,a,area,,60
Asel,a,area,,63
Asel,a,area,,65
vDRAG,ALL,,,,,,100
type,2
real,1
mat,16
Asel,s,area,,70
Asel,a,area,,74
Asel,a,area,,78
Asel,a,area,,81
vDRAG,ALL,,,,,,100
type,2
real,1
mat,17
Asel,s,area,,86
Asel,a,area,,90
Asel,a,area,,94
Asel,a,area,,97
vDRAG,ALL,,,,,,100
type,2
real,1
mat,5
Asel,s,area,,103
Asel,a,area,,108
Asel,a,area,,113
Asel,a,area,,117
vDRAG,ALL,,,,,,100
type,2
real,1
mat,6
Asel,s,area,,123
Asel,a,area,,128
Asel,a,area,,133
Asel,a,area,,137
vDRAG,ALL,,,,,,100
vsel,s,mat,,10,11,1
vsel,r,loc,y,y1,y4
vplot
TYPE,3！
壳单元
REAL,1
MAT,12
ASEL,S,AREA,,457
ASEL,A,AREA,,461
ASEL,A,AREA,,465
ASEL,A,AREA,,468
ASEL,A,AREA,,473
ASEL,A,AREA,,477
ASEL,A,AREA,,481
ASEL,A,AREA,,484
ASEL,A,AREA,,592
ASEL,A,AREA,,595
ASEL,A,AREA,,598
ASEL,A,AREA,,600
ASEL,A,AREA,,604
ASEL,A,AREA,,607
ASEL,A,AREA,,610
ASEL,A,AREA,,612
ASEL,A,AREA,,664
ASEL,A,AREA,,667
ASEL,A,AREA,,670
ASEL,A,AREA,,672
ASEL,A,AREA,,676
ASEL,A,AREA,,679
ASEL,A,AREA,,682
ASEL,A,AREA,,684
APLOT
MSHAPE,0,2D
MSHKEY,1
AMESH,ALL
TYPE,3
REAL,1
MAT,13
ASEL,S,AREA,,44
ASEL,A,AREA,,47
ASEL,A,AREA,,50
ASEL,A,AREA,,52
ASEL,A,AREA,,247
ASEL,A,AREA,,250
ASEL,A,AREA,,253
ASEL,A,AREA,,255
ASEL,A,AREA,,315
ASEL,A,AREA,,318
ASEL,A,AREA,,321
ASEL,A,AREA,,323
ASEL,A,AREA,,753
ASEL,A,AREA,,756
ASEL,A,AREA,,759
ASEL,A,AREA,,761
APLOT
MSHAPE,0,2D
MSHKEY,1
AMESH,ALL
TYPE,3
REAL,1
MAT,14
ASEL,S,AREA,,56
ASEL,A,AREA,,59
ASEL,A,AREA,,62
ASEL,A,AREA,,64
ASEL,A,AREA,,259
ASEL,A,AREA,,262
ASEL,A,AREA,,265
ASEL,A,AREA,,267
ASEL,A,AREA,,327
ASEL,A,AREA,,330
ASEL,A,AREA,,333
ASEL,A,AREA,,335
ASEL,A,AREA,,765
ASEL,A,AREA,,768
ASEL,A,AREA,,771
ASEL,A,AREA,,773
APLOT
MSHAPE,0,2D
MSHKEY,1
AMESH,ALL
nummrg,all,,,,low
numcmp,all
！
！
！
！
（重要一步）手工操作旋转及复制实体，需要先转换为柱坐标，而后旋转并平行移动.方法：
键入csys,5（柱坐标）,然后点击operation下的move/modify，选中要旋转的实体，在xyz中键入旋转角度，例y格填入90，为绕y轴旋转90度。
结合移动、复制、及旋转可建完整个实体。
（注意旋转后移动实体需改回正常坐标，即键入csys,0）
esel,s,mat,,2,4,1！
选定区域土体改性，改为加固
NSLE,S
nplot
nsel,r,loc,z,-6,0
nsel,r,loc,x,x2,x6
nsel,r,loc,y,-28,y8
ESLN,S,1
Eplot
Mpchg,21,all
esel,s,mat,,2,4,1
NSLE,S
nplot
nsel,r,loc,z,-16.25,-15
nsel,r,loc,x,x2,x6
nsel,r,loc,y,-28,y8
ESLN,S,1
Eplot
Mpchg,22,all
esel,s,mat,,2,4,1
NSLE,S
nplot
nsel,r,loc,x,x2,x2+1.25
nsel,a,loc,x,x6-1.25,x6
nsel,a,loc,x,27.5,30
nsel,r,loc,z,-15,-6
nsel,r,loc,y,-28,y8
ESLN,S,1
Eplot
Mpchg,22,all
esel,s,mat,,3
NSLE,S
nplot
nsel,r,loc,z,z2-1.25,z6
nsel,r,loc,x,x2-6,x6+6
nsel,r,loc,y,y2,y3
ESLN,S,1
Eplot
Mpchg,23,all
/solu！
进入求解
Antype,static！
约束边界
nsel,s,loc,x,x1
nsel,a,loc,x,x7
D,all,ux,0
Alls
nsel,s,loc,y,y1
D,all,all,0
Alls
nsel,s,loc,z,z1
nsel,a,loc,z,z7-15
D,all,uz,0
Allsel
Acel,,10
Deltim,0.1,0.05,0.3
Autots,on
Pred,on
Lnsrch,on
Nlgeom,on
Nropt,full
Cnvtol,f,,0.02,2,0.5
esel,s,type,,3！
杀死壳单元计算自重
ekill,all
Allsel,all,all
/PBC,ALL,,1
GPLOT
TIME,1
iswrite,on
lswrite,1
SOLVE
esel,s,mat,,9！
开挖1号线计算
ekill,all
esel,s,mat,,12！
激活1号线壳体
ealive,all
Allsel,all,all
TIME,2
iswrite,on
lswrite,2
SOLVE
Esel,s,mat,,10！
杀死7号线左侧下行线开挖段土体
Eplot
NSLE,S
nplot
nsel,r,loc,z,-6,0
nplot
ESLN,S,1
Eplot
Ekill,all
Esel,s,mat,,13！
激活壳体单元
Eplot
NSLE,S
nplot
nsel,r,loc,z,-6,0
nplot
ESLN,S,1
Eplot
Ealive,all
Esel,s,mat,,16！
改注浆材料
Eplot
NSLE,S
nplot
nsel,r,loc,z,-6,0
nplot
ESLN,S,1
Eplot
Mpchg,19,all
Esel,s,mat,,10！
开挖面加压
Eplot
Esel,r,live
Eplot
Nsle,s
nplot
NSEL,r,loc,z,-6
nplot
sf,all,pres,342000!
加压为node上加压
Allsel,all,all
TIME,3
iswrite,on
lswrite,3
SOLVE
Esel,s,mat,,10
Eplot
NSLE,S
nplot
nsel,r,loc,z,-12,-6
nplot
ESLN,S,1
Eplot
Ekill,all
Esel,s,mat,,13
Eplot
NSLE,S
nplot
nsel,r,loc,z,-12,-6
nplot
ESLN,S,1
Eplot
Ealive,all
Esel,s,mat,,16
Eplot
NSLE,S
nplot
nsel,r,loc,z,-12,-6
nplot
ESLN,S,1
Eplot
Mpchg,19,all
Esel,s,mat,,10
Eplot
Esel,r,live
Eplot
Nsle,s
nplot
NSEL,r,loc,z,-12
nplot
sf,all,pres,342000
Allsel,all,all
TIME,4
iswrite,on
lswrite,4
SOLVE
Esel,s,mat,,10
Eplot
NSLE,S
nplot
nsel,r,loc,z,-21.25,-12
nplot
ESLN,S,1
Eplot
Ekill,all
Esel,s,mat,,13
Eplot
NSLE,S
nplot
nsel,r,loc,z,-21.25,-12
nplot
ESLN,S,1
Eplot
Ealive,all
Esel,s,mat,,16
Eplot
NSLE,S
nplot
nsel,r,loc,z,-21.25,-12
nplot
ESLN,S,1
Eplot
Mpchg,19,all
Esel,s,mat,,10
Eplot
Esel,r,live
Eplot
Nsle,s
nplot
NSEL,r,loc,z,-21.25
nplot
sf,all,pres,342000
Allsel,all,all
TIME,5
iswrite,on
lswrite,5
SOLVE
Esel,s,mat,,10
Eplot
NSLE,S
nplot
nsel,r,loc,z,-36.25,-21.25
nplot
ESLN,S,1
Eplot
Ekill,all
Esel,s,mat,,13
Eplot
NSLE,S
nplot
nsel,r,loc,z,-36.25,-21.25
nplot
ESLN,S,1
Eplot
Ealive,all
Esel,s,mat,,16
Eplot
NSLE,S
nplot
nsel,r,loc,z,-36.25,-21.25
nplot
ESLN,S,1
Eplot
Mpchg,19,all
Esel,s,mat,,10
Eplot
Esel,r,live
Eplot
Nsle,s
nplot
NSEL,r,loc,z,-36.25
nplo
升级会员 