/TITLE,Mechanical analysis on railway tunnel 2nd lining ! 确定分析标题
/NOPR !菜单过滤设置
/PMETH,OFF,0
KEYW,PR_SET,1
KEYW,PR_STRUC,1 !保留结构分析部分菜单
/COM,
/COM,Preferences for GUI filtering have been set to display:
/COM, Structural
/PREP7 !进入前处理器
ET,1,BEAM3 !设置梁单元类型
ET,2, COMBIN14 !设置弹簧单元类型
R,1,0.4,0.0053333,0.4,,,, !设置梁单元几何常数
R,2,400e6,,, !设置弹簧单元几何常数
MPTEMP,,,,,,,, !设置材料模型
MPTEMP,1,0
MPDATA,EX,1,,30.0e9 !输入弹性模量
MPDATA,PRXY,1,,0.2 !输入泊松比
MPTEMP,,,,,,,, !设置材料模型
MPTEMP,1,0
MPDATA,DENS,1,,2500 !输入密度
SAVE !保存数据库
K, 1, , , , !创建关键点（隧道二次衬砌）
K, 2, 4.71, -1.82, ,
K, 3, -4.71, -1.82, ,
K,4,0,-3.75,,
K, 5, 0, 5.05, ,
K, 20, 5.6, -2.162, , !创建关键点（地层弹簧）
K, 30, -5.6, -2.16, ,
K, 40, 0, -4.75, ,
K, 50, 0, 6, , !创建隧道衬砌线
LARC, 2, 5, 1, 5.05, !创建圆弧线（拱顶部）
LARC, 5, 3, 1, 5.05,
LARC, 2, 3, 4 !创建圆弧线（仰拱部）! 创建地层弹簧线
LARC, 20, 50, 1, 6, !创建圆弧线（拱顶部）
LARC,50, 30, 1, 6,
LARC, 20, 30, 40 !创建圆弧线（仰拱部）
!HPTCREATE, TYPE, ENTITY, NHP, LABEL, VAL1, VAL2, VAL3
!HPTCREATE,LINE,3, 444, COORD,0, -3.75, ,
!HPTCREATE,LINE,2, 445, COORD, -4.71, -1.82, ,
!HPTCREATE,LINE,1, 446, COORD,4.71, -1.82, ,
SAVE !保存数据
LESIZE,ALL, , ,10, ,1, , ,1, !设置单元大小, 每条弧线划分成10个单元
TYPE, 1 !设置将要创建单元的类型
MAT, 1 !设置将要创建单元的材料
REAL, 1 !设置将要创建单元的几何常数
LMESH,1,6,1 !将所有直线划分单元
/PNUM, KP, 0 !以下为显示单元编号和颜色
/PNUM, ELEM, 1
/REPLOT !重新显示
SAVE
TYPE, 2 !设置将要创建单元的类型
MAT, 1 !设置将要创建单元的材料
REAL, 2 !设置将要创建单元的几何常数
E, 1,31 !通过两个节点创建弹簧单元
E, 2,32
E,3,33
E,4,34
E,5,35
E,6,36
E,7,37
E,8,38
E,9,39
E,10,40
E,11,41
E,12,42
E,13,43
E,14,44
E,15,45
E,16,46
E,17,47
E,18,48
E,19,49
E,20,50
E,21,51
E,22,52
E,23,53
E,24,54
E,25,55
E,26,56
E,27,57
E,28,58
E,29,59
E,30,60
Lclear,4,6,1 !清除用于创建地层弹簧单元的外层梁单元
LDELE,4, 6,1 ,1 !删除外层直线
Finish !返回Main Menu 主菜单
/SOL !进入求解器
d,31,
ux,0,0,60,1,uy, !在“Ux”和“Uy”两个方向的施加约束
ACEL, 0, 10, 0, !在Y方向施加重力加速度
F, 2 ,Fx, 0 !在节点上施加X方向集中力
F, 13 ,Fx, 43554.43356
F, 14 ,Fx, 86574.64508
F, 15 ,Fx, 128464.6048
F, 16 ,Fx, 168285.8813
F, 17 ,Fx, 204799.5937
F, 18 ,Fx, 236448.4444
F, 19 ,Fx, 261396.9365
F, 20 ,Fx, 277795.0781
F, 21 ,Fx, 283894.0611
F, 12 ,Fx, 247397.3915
F, 30 ,Fx, 190308.0663
F, 29 ,Fx, 149683.6899
F, 28 ,Fx, 103226.8339
F, 27 ,Fx, 52665.0399
F, 26 ,Fx, 0
F, 25 ,Fx, -52665.0399
F, 24 ,Fx, -103226.8339
F, 23 ,Fx, -149683.6899
F, 22 ,Fx, -190308.0663
F, 1 ,Fx, -247397.3915
F, 3 ,Fx, -283894.0611
F, 4 ,Fx, -277795.0781
F, 5 ,Fx, -261396.9365
F, 6 ,Fx, -236448.4444
F, 7 ,Fx, -204799.5937
F, 8 ,Fx, -168285.8813
F, 9 ,Fx, -128464.6048
F, 10 ,Fx, -86574.64508
F, 11 ,Fx, -43554.43356
F, 2 ,Fy, -449636.88 !在节点上施加Y方向集中力
F, 13 ,Fy, -441210
F, 14 ,Fy, -416229.66
F, 15 ,Fy, -375629.1
F, 16 ,Fy, -320951.4
F, 17 ,Fy, -254238.6
F, 18 ,Fy, -178008.6
F, 19 ,Fy, -95102.7
F, 20 ,Fy, 0
F, 21 ,Fy, 82772.64
F, 12 ,Fy, 259324.92
F, 30 ,Fy, 413398.46
F, 29 ,Fy, 454320.04
F, 28 ,Fy, 484308
F, 27 ,Fy, 502604.08
F, 26 ,Fy, 508719.08
F, 25 ,Fy, 502604.08
F, 24 ,Fy, 484308
F, 23 ,Fy, 454320.04
F, 22 ,Fy, 413398.46
F, 1 ,Fy, 259324.92
F, 3 ,Fy, 82772.64
F, 4 ,Fy, 0
F, 5 ,Fy, -95102.7
F, 6 ,Fy, -178008.6
F, 7 ,Fy, -254238.6
F, 8 ,Fy, -320951.4
F, 9 ,Fy, -375629.1
F, 10 ,Fy, -416229.66
F, 11 ,Fy, -441210
Save
NROPT, FULL,, !采用全牛顿－拉普森法进行求解
Allsel !选择所有内容
Outres, all, all !输出所有内容
Solve !求解计算
Finish !求解结束返回Main Menu 主菜单
SAVE
/POST1 !进入后处理器
PLDISP, 1 !绘制变形和未变形图