ANSYS基坑开挖模拟命令流

finish
/clear
/filname,ansys2,1
/title,the analysis of gangkou hosipital
/units,si
/prep7
et,1,solid45 !土体
et,2,shell181 !地连墙
et,3,beam4 !支撑
r,1 !土体
r,2,0.62,,,, !地连墙
r,3,0.65,,,,
r,4,0.73,,,,
r,5,0.6,0.6/12,1*0.6*0.6*0.6/12,0.6,1,, !冠梁
r,6,0.4,0.8*0.5*0.5*0.5/12,0.5*0.8*0.8*0.8/12,0.8,0.5,, !腰撑
r,7,0.64,0.8*0.8*0.8*0.8/12,0.8*0.8*0.8*0.8/12,0.8,0.8,, !支撑梁
r,8,1.44,0.8*1.8*1.8*1.8/12,1.8*0.8*0.8*0.8/12,0.8,1.8,, !环梁及对撑梁
r,9,0.2,,,, !加板

mp,dens,1,1830 !第一层土
mp,ex,1,3.5e7
mp,prxy,1,0.25
tb,dp,1
tbdata,,2.34e4,6.8,6.8
mp,dens,2,1870 !第二层土
mp,ex,2,4.62e7
mp,prxy,2,0.3
tb,dp,2
tbdata,,1.52e4,5.3,5.3
mp,dens,3,1930 !第三层土
mp,ex,3,4.62e7
mp,prxy,3,0.3
tb,dp,3
tbdata,,1.35e4,10.8,10.8
mp,dens,4,1710 !第四层土
mp,ex,4,2.86e7
mp,prxy,4,0.2
tb,dp,4
tbdata,,1.1e4,3.3,3.3
mp,dens,5,1780 !第五层土
mp,ex,5,5.4e7
mp,prxy,5,0.2
tb,dp,5
tbdata,,1.35e4,3.5,3.5
mp,dens,6,1830 !第六层土
mp,ex,6,5.4e7
mp,prxy,6,0.2
tb,dp,6
tbdata,,1.9e4,5.7,5.7
mp,dens,7,2040 !第七层土
mp,ex,7,5.23e7
mp,prxy,7,0.3
tb,dp,7
tbdata,,2.18e4,16.3,16.3
mp,dens,8,2500 !地连墙和支撑
mp,ex,8,3e10
mp,prxy,8,0.2

blc4,,,145.2,138.4,-30
/pnum,volu,1
vplot
!*************开始沿Z轴方向切分
wpoff,,,-1.3
vsbw,all
wpoff,,,-0.5
vsbw,all
wpoff,,,-1.1
vsbw,all
wpoff,,,-0.5
vsbw,all
wpoff,,,-1.6
vsbw,all
wpoff,,,-2.5
vsbw,all
wpoff,,,-2.1
vsbw,all
wpoff,,,-6.3
vsbw,all
wpoff,,,-3.4
vsbw,all
wpoff,,,-2.6
vsbw,all
wpoff,,,-1.9
vsbw,all
wpoff,,,-1
vsbw,all
numcmp,all
nummrg,all
allsel
!*************开始沿X轴方向切分
wpcsys,-1
wprota,,,90

wpoff,,,19
vsbw,all
wpoff,,,4
vsbw,all
wpoff,,,2
vsbw,all
*do,i,1,6
wpoff,,,5
vsbw,all
*enddo
*do,i,1,2
wpoff,,,7.4
vsbw,all
*enddo
*do,i,1,9
wpoff,,,5
vsbw,all
*enddo
wpoff,,,5.4
vsbw,all
wpoff,,,9.4
vsbw,all
wpoff,,,8.6
vsbw,all
numcmp,all
nummrg,all
allsel
!*************开始沿Y轴方向切分
wpcsys,-1
wprota,,-90
wpoff,,,7.3
vsbw,all
wpoff,,,13
vsbw,all
wpoff,,,4.7
vsbw,all
*do,i,1,3
wpoff,,,5
vsbw,all
*enddo
wpoff,,,2.8
vsbw,all
wpoff,,,2.2
vsbw,all
*do,i,1,2
wpoff,,,5
vsbw,all
*enddo
*do,i,1,2
wpoff,,,3.7
vsbw,all
*enddo
wpoff,,,3.8
vsbw,all
wpoff,,,6
vsbw,all
wpoff,,,3.8
vsbw,all
wpoff,,,7.4
vsbw,all
*do,i,1,3
wpoff,,,5
vsbw,all
*enddo
wpoff,,,2.4
vsbw,all
wpoff,,,2.6
vsbw,all
*do,i,1,2
wpoff,,,5
vsbw,all
*enddo
wpoff,,,2
vsbw,all
wpoff,,,4
vsbw,all

numcmp,all
nummrg,all

/pnum,volu,0
allsel ! 8112个体

!****************保存1
!****************赋予体属性
wpcsys,-1
vsel,s,loc,z,0,-3.4
vatt,1,1,1
vsel,s,loc,z,-3.4,-5
vatt,2,1,1
/color,volu,2
vsel,s,loc,z,-5,-9.6
vatt,3,1,1
vsel,s,loc,z,-9.6,-15.9
vatt,4,1,1
/color,volu,2
vsel,s,loc,z,-15.9,-19.3
vatt,5,1,1
vsel,s,

loc,z,-19.3,-21.9
vatt,6,1,1
/color,volu,2
vsel,s
,loc,z,-21.9,-30
vatt,7,1,1
allsel
!****************赋予面属性
asel,s,loc,x,25 !桩1
asel,r,loc,z,-1.3,-23.8
asel,r,loc,y,25,113.4
aatt,8,2,2
asel,s,loc,y,113.4 !桩1
asel,r,loc,x,25,89.8
asel,r,loc,z,-1.3,-23.8
aatt,8,2,2
asel,s,loc,y,113.4 !桩5
asel,r,loc,x,89.8,109.8
asel,r,loc,z,-1.3,-24.8
aatt,8,3,2
asel,s,loc,y,113.4 !桩1
asel,r,loc,x,109.8,120.2
asel,r,loc,z,-1.3,-23.8
aatt,8,2,2

asel,s,loc,x,120.2 !桩2
asel,r,loc,y,76,113.4
asel,r,loc,z,-1.3,-23.8
aatt,8,3,2
asel,s,loc,y,76 !桩2
asel,r,loc,x,99.8,120.2
asel,r,loc,z,-1.3,-23.8
aatt,8,3,2
asel,s,loc,x,99.8 !桩6
asel,r,loc,y,42.8,76
asel,r,loc,z,-1.3,-24.8
aatt,8,2,2
asel,s,loc,x,99.8 !桩3,4
asel,r,loc,y,25,42.8
asel,r,loc,z,-1.3,-24.8
aatt,8,4,2
asel,s,loc,y,25 !桩3,4
asel,r,loc,x,25,99.8
asel,r,loc,z,-1.3,-24.8
aatt,8,4,2
allsel

numcmp,all
nummrg,all

wpcsys,-1
asel,s,type,,2

!*****************赋予冠梁属性
lplot
lsel,s,loc,z,-1.3
lsel,r,loc,x,25
lsel,r,loc,y,25,113.4
latt,8,5,3

lsel,s,loc,z,-1.3
lsel,r,loc,y,113.4
lsel,r,loc,x,25,120.2
latt,8,5,3

lsel,s,loc,z,-1.3
lsel,r,loc,x,120.2
lsel,r,loc,y,76,113.4
latt,8,5,3

lsel,s,loc,z,-1.3
lsel,r,loc,y,76
lsel,r,loc,x,99.8,120.2
latt,8,5,3

lsel,s,loc,z,-1.3
lsel,r,loc,x,99.8
lsel,r,loc,y,25,76
latt,8,5,3

lsel,s,loc,z,-1.3
lsel,r,loc,y,25
lsel,r,loc,x,25,99.8
latt,8,5,3

lsel,s,real,,5 !验证看一下
allsel
!*****************赋予腰梁属性
wpcsys,-1
asel,s,type,,2
aplot

lplot
lsel,s,loc,z,-2.9
lsel,r,loc,x,25
lsel,r,loc,y,25,113.4
latt,8,6,3

lsel,s,loc,z,-2.9
lsel,r,loc,y,113.4
lsel,r,loc,x,25,120.2
latt,8,6,3

lsel,s,loc,z,-2.9
lsel,r,loc,x,120.2
lsel,r,loc,y,76,113.4
latt,8,6,3

lsel,s,loc,z,-2.9
lsel,r,loc,y,76
lsel,r,loc,x,99.8,120.2
latt,8,6,3

lsel,s,loc,z,-2.9
lsel,r,loc,x,99.8
lsel,r,loc,y,25,76
latt,8,6,3

lsel,s,loc,z,-2.9
lsel,r,loc,y,25
lsel,r,loc,x,25,99.8
latt,8,6,3

lsel,s,real,,6 !验证看一下
allsel

!****************创建内支撑系统
!****************先创建上半圆弧及支撑梁
lsel,s,real,,6
/pnum,kp,1

k,,62.4,76,-2.9 !创建点9451，圆心
circle,9451,37.4,,,180 !此时三个相切点重合

l,8479,9451

numcmp,all !合并重合的关键点
nummrg,all

l,8488,5758
l,8507,5758
l,8516,5758
l,8534,5758
l,8548,5758
l,8576,5758
l,8588,5758
l,8606,5758
l,8618,5758
l,8636,5758
l,8648,5758
l,8662,5758
l,8674,5758
l,8694,5758
l,8706,5758
l,8374,5758
l,8024,5758
l,7324,5758
l,6974,5758
l,6624,5758
l,6274,5758

l,5996,5758
l,6346,5758
l,6696,5758
l,7046,5758
l,7746,5758
l,8096,5758

/pnum,line,1

lovlap,26949,26972,26973,26974,26975,26976,26977
lovlap,26991,26950,26951,26952,26953,26954,26955
lovlap,26948,26956,26957,26958,26959,26960,26961
lovlap,26991,26962,26963,26964,26965
lovlap,26959,26966,26967,

26968,26969,26970,26971

*do,i,26979,26984
ldele,i
*enddo
*do,i,26972,26976
ldele,i
*enddo
ldele,2
6949
*do,i,26950,26958
ldele,i
*enddo
ldele,26948
*do,i,26962,26965
ldele,i
*enddo
ldele,26991
ldele,27027

numcmp,all
nummrg,all

/pnum,line,0

l,5906,6274 !阳角部分斜撑杆
l,5894,6624
l,5874,6974
l,5862,7324

l,8374,8694 !右上角斜撑杆
l,7324,8662
l,6624,8636

l,8488,7746 !左上角斜撑杆
l,8516,6696

/pnum,kp,0
/pnum,line,1

lovlap,27013,26957,26950

lovlap,27014,26955,26956,27019,27020,26963,26964

lovlap,27010,26988

lovlap,27011,26949,26987,26964,26999,26998

lovlap,27012,26981,26948,27039,27040,27041

lovlap,26949,27042,27043,26997,27009,26996,27008

lovlap,27007,26994

lovlap,26981,27008,26995

lovlap,26948,27009,27042

numcmp,all
nummrg,all
!****************再创建下半圆弧及支撑梁

/pnum,kp,1
/pnum,line,0

k,,62.4,62.4,-2.9 !创建点9510，圆心
circle,9510,37.4,,,

l,1479,9510

/pnum,line,1

ldele,27078
ldele,27077

numcmp,all
nummrg,all

/pnum,line,0

l,1488,4708
l,1507,4708
l,1516,4708
l,1534,4708
l,1548,4708
l,1576,4708
l,1588,4708
l,1606,4708
l,1618,4708
l,1636,4708
l,1648,4708
l,2012,4708
l,2362,4708
l,2712,4708
l,3412,4708
l,3762,4708
l,4112,4708

l,1796,4708
l,2146,4708
l,2496,4708
l,3196,4708
l,3546,4708
l,3896,4708

/pnum,line,1

lovlap,27077,27097,27098,27099,27100,27101,27102
lovlap,27103,27079,27080,27081,27082,27083,27084
lovlap,27078,27085,27086,27087,27088,27089,27090
lovlap,27140,27091,27092,27093,27094,27095,27096

*do,i,27111,27116
ldele,i
*enddo

*do,i,27122,27127
ldele,i
*enddo

*do,i,27134,27139
ldele,i
*enddo

*do,i,27146,27150
ldele,i
*enddo

ldele,27090

numcmp,all
nummrg,all

/pnum,line,0

l,1488,2146
l,1516,3196

l,1636,2362
l,1606,3412

/pnum,line,1

lovlap,27127,27097,27077

lovlap,27128,27099,27098,27135,27136,27090,27091

lovlap,27129,27096,27078

lovlap,27130,27083,27084,27128,27135,27085,27086

numcmp,all
nummrg,all

!****************再创建中间支撑梁

/pnum,line,0

k,,30,72.2,-2.9 !9551-9557
k,,40.8,72.2,-2.9
k,,51.6,72.2,-2.9
k,,62.4,72.2,-2.9
k,,73.2,72.2,-2.9
k,,84,72.2,-2.9
k,,94.8,72.2,-2.9

k,,30,66.2,-2.9 !9558-9564
k,,40.8,66.2,-2.9
k,,51.6,66.2,-2.9
k,,62.4,66.2,-2.9
k,,73.2,66.2,-2.9
k,,84,66.2,-2.9
k,,94.8,66.2,-2.9

l,5296,9551 !上半对撑27163-27170

*do,i,9551,9556
l,i,i+1
*enddo

l,9557,5512

l,4946,9558 !下半对撑27171-27178

*do,i,9558,9563
l,i,i+1
*enddo

l,9564,5162


l,9551,9558 ! 7个小竖撑
l,9552,9559
l,9553,9560
l,9554,9561
l,9555,9562
l,9556,9563
l,9557,9564

l,5646,9551 !四个小斜杆
l,9558,4596

l,9557,5862
l,9564,4812 !线27189

/pnum,line,0

!************开始创建要加设板的面域

a,8548,8558,9462 !上面六个撑板，面25611-25616
a,8534,8548,9462

,9461
a,8516,8534,9461,9460

a,8558,8576,9465
a,8576,8588,9466,9465
a,8588,8606,9467,9466

a,5996,9451,5646 !左面六个撑板，面25617-25622
a,
6346,9454,9451,5996
a,6696,9455,9454,6346

a,4596,9516,3896
a,3896,9516,9515,3546
a,3546,9515,9514,3196

a,5646,9551,5296 !剪力最大处加设面，25623-25624
a,4946,9558,4596

a,9522,1558,1548 !下面六个撑板，面25625-25630
a,9521,9522,1548,1534
a,9520,9521,1534,1516

a,1558,9523,1576
a,9523,9524,1588,1576
a,9524,9525,1606,1588

a,4812,4112,9534 !右边三个撑板，面25631-25633
a,9534,4112,3762,9533
a,9533,3762,3412,9532

a,5862,5512,9557 !剪力最大处加设面，25634-25635
a,9564,5162,4812

asel,s,area,,25611,25635,1 !赋予撑杆属性
aplot
aatt,8,9,2
allsel

!**************赋予内支撑线属性
!**************先赋值环梁及对撑

lsel,s,line,,26950
lsel,a,line,,26999
lsel,a,line,,27079
lsel,a,line,,27091
lsel,a,line,,26957,26961
lsel,a,line,,26967,26974
lsel,a,line,,26981,26985
lsel,a,line,,26988,26992
lsel,a,line,,27000,27004
lsel,a,line,,27099,27122
lsel,a,line,,27163,27178
lplot
latt,8,8,3

!**********************普通撑杆赋予属性

lsel,s,line,,26948,27189,1
lsel,u,real,,8
latt,8,7,3

allsel
aplot

!************到此为止所有体、面、线均赋予属性，保存3
!************挖去多余土体
/pnum,kp,0
/pnum,line,0
/pnum,volu,0
vplot

vsel,s,loc,x,30,138.2 !开始删除土体1
vsel,r,loc,y,7.3,20.3
vsel,r,loc,z,0,-1.3
vdele,all
allsel

vsel,s,loc,x,99.8,145.2 !开始删除土体2
vsel,r,loc,y,42.8,58.7
vsel,r,loc,z,0,-1.8
vdele,all
allsel

vsel,s,loc,x,129.6,145.2 !开始删除土体3
vsel,r,loc,y,66.2,100.8
vsel,r,loc,z,0,-1.3
vdele,all
allsel


!************接下来划分线，划分网格
lplot
!************X方向土体划分设置
allsel
wpcsys,-1

lsel,s,loc,x,1,15
lesize,all,4 !每4m一段

lsel,s,loc,x,20,22
lesize,all,2 !每2m一段

lsel,s,loc,x,23.5,24.5
lesize,all,2 !每2m一段

lsel,s,loc,x,26,29
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,31,34
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,36,39
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,41,44
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,46,49
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,51,54
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,56,62
lsel,u,real,,7
lsel,u,real,,8
lesize,all,3.7 !每3.7m一段

lsel,s,loc,x,63,69
lsel,u,real,,7
lsel,u,real,,8
lesize,all,3.7 !每3.7m一段

lsel,s,loc,x,70,74
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,75,79
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,80,84
lsel,u,real,,7

lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,85,89
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,90,94
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一
段

lsel,s,loc,x,95,99
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,100,104
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,105,109
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,110,114
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,x,115,120
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.7 !每2.7m一段

lsel,s,loc,x,121,129
lesize,all,4.7

lsel,s,loc,x,130,138
lesize,all,4.3

lsel,s,loc,x,139,145
lesize,all,3.5


!************Y方向土体划分设置
lplot
allsel

lsel,s,loc,y,1,7
lesize,all,3.65 !每3.65m一段

lsel,s,loc,y,8,20
lesize,all,,,4 !划分成4段

lsel,s,loc,y,21,24
lesize,all,,,2 !划分成2段

lsel,s,loc,y,26,29
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,y,31,34
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,y,36,39
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,y,41,42.5
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.8 !每2.8m一段

lsel,s,loc,y,43,44
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.2 !每2.2m一段

lsel,s,loc,y,46,49
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,y,51,54
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,y,56,58
lsel,u,real,,7
lsel,u,real,,8
lesize,all,3.7 !每3.7m一段

lsel,s,loc,y,59,62
lsel,u,real,,7
lsel,u,real,,8
lesize,all,3.7 !每3.7m一段

lsel,s,loc,y,63,66
lsel,u,real,,7
lsel,u,real,,8
lesize,all,1.9 !每1.9m一段

lsel,s,loc,y,67,72
lsel,u,real,,7
lsel,u,real,,8
lesize,all,3 !每3m一段

lsel,s,loc,y,73,75
lsel,u,real,,7
lsel,u,real,,8
lesize,all,1.9 !每1.9m一段

lsel,s,loc,y,77,83
lsel,u,real,,7
lsel,u,real,,8
lesize,all,3.7 !每3.7m一段


lsel,s,loc,y,84,88
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,y,89,93
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,y,94,98
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,y,99,100
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.4 !每2.4m一段


lsel,s,loc,y,101,103
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.6 !每2.6m一段

lsel,s,loc,y,104,108
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,y,109,113
lsel,u,real,,7
lsel,u,real,,8
lesize,all,2.5 !每2.5m一段

lsel,s,loc,y,114,115
lesize,all,2 !每2m一段

lsel,s,loc,y,116,119
lesize,all,2 !每2m一段

lsel,s,loc,y,120,138
lesize,all,4 !每4m一段
allsel

!************Z方向土体划分设置
lsel,s,loc,z,-0.

1,-1.2
lesize,all,,,1
lsel,s,loc,z,-1.4,-1.7
lesize,all,,,1

lsel,s,loc,z,-1.9,-2.8
lesize,all,,,1

lsel,s,loc,z,-3,-3.3
lesize,all,,,1

lsel,s,loc,z,-3.5,-4.9
lesize,all,,,2

lsel,s,loc,z,-5.1,-7.4
lesize,all,,,2

lsel,s,loc,z,-7.6,-9.5
lesize,all,,
,1

lsel,s,loc,z,-9.7,-15.8
lesize,all,,,3

lsel,s,loc,z,-16,-19.2
lesize,all,,,2

lsel,s,loc,z,-19.4,-21.8
lesize,all,,,1

lsel,s,loc,z,-22,-23.7
lesize,all,,,1

lsel,s,loc,z,-23.9,-24.7
lesize,all,,,1

lsel,s,loc,z,-24.9,-29
lesize,all,,,2
allsel
!************环梁及支撑梁划分设置

lsel,s,real,,7
lsel,a,real,,8

lplot
lesize,all,,,2

!******************开始划分体、面、线
allsel
vplot

mshape,0,2d
mshkey,1
vmesh,all !划分8023个体

asel,s,type,,2
aplot

amesh,all !有提示，21个面划分不太好

lsel,s,real,,5
lsel,a,real,,6
lsel,a,real,,7
lsel,a,real,,8
lplot

lmesh,all

allsel
eplot

!esel,s,real,,5
!esel,s,real,,6
!esel,s,real,,7
!esel,s,real,,8
!esel,s,type,,2
allsel

!******************保存2


/psf,pres,norm,2,0,1 !设置荷载显示方式为箭头

asel,s,loc,x,30,138.2
aplot
asel,r,loc,y,7.3,20.3
asel,r,loc,z,-1.3
nsla,s,1
nplot
sf,all,pres,75000 !面1,施加荷载75000
allsel

asel,s,loc,x,99.8,145.2
asel,r,loc,y,42.8,58.7
asel,r,loc,z,-1.8
nsla,s,1
sf,all,pres,40000 !面2,施加荷载40000
allsel

asel,s,loc,x,129.6,145.2
asel,r,loc,y,66.2,100.8
asel,r,loc,z,-1.3
nsla,s,1
sf,all,pres,40000 !面3,施加荷载40000
allsel

asel,s,loc,x,19,23
asel,r,loc,y,25,113.4
asel,r,loc,z,0
nsla,s,1
sf,all,pres,30000 !面4,施加荷载30000
allsel

asel,s,loc,x,25,120.2
asel,r,loc,y,115.4,119.4
asel,r,loc,z,0
nsla,s,1
sf,all,pres,30000 !面5,施加荷载30000


allsel
eplot

!**************被挖去的三个土体施加约束

nsel,s,loc,x,30 !删除体1-1
nsel,r,loc,y,7.3,20.3
nsel,r,loc,z,0,-1.3
d,all,ux

allsel

nsel,s,loc,x,138.2 !删除体1-2
nsel,r,loc,y,7.3,20.3
nsel,r,loc,z,0,-1.3
d,all,ux

allsel

nsel,s,loc,y,7.3 !删除体1-3
nsel,r,loc,x,30,138.2
nsel,r,loc,z,0,-1.3
d,all,uy

allsel

nsel,s,loc,y,20.3 !删除体1-4
nsel,r,loc,x,30,138.2
nsel,r,loc,z,0,-1.3
d,all,uy

allsel

nsel,s,loc,x,99.8 !删除体2-1
nsel,r,loc,y,42.8,58.7
nsel,r,loc,z,0,-1.8
d,all,ux

allsel

nsel,s,loc,x,145.2 !删除体2-2
nsel,r,loc,y,42.8,58.7
nsel,r,loc,z,0,-1.8
d,all,ux

allsel

nsel,s,loc,y,42.8 !删除体2-3
nsel,r,loc,x,99.8,145.2
nsel,r,loc,z,0,-1.8
d,all,uy

allsel

nsel,s,loc,y,58.7 !删除体2-4
nsel,r,loc,x,99.8,145.2
nsel,r,loc,z,0,-1.8
d,all,uy

allsel

nsel,s,loc,x,129.6 !删除体3-1
nsel,r,loc,y,66.2,100.8
nsel,r,loc,z,0,-1.3
d,all,ux

allsel

nsel,s,loc,x,145.2 !删除体3-2
nsel,r,loc,y,66.2,100.8
nsel,r,loc,z,0,-1.3
d,all,ux

allsel

nsel,s,loc,y,66.2 !删除体3-3
nsel,r,

loc,x,129.6,145.2
nsel,r,loc,z,0,-1.3
d,all,uy

allsel

nsel,s,loc,y,100.8 !删除体3-4
nsel,r,loc,x,129.6,145.2
nsel,r,loc,z,0,-1.3
d,all,uy

allsel

nsel,s,loc,z,-30 !施加约束
d,all,all

nsel,s,loc,x,0
nsel,a,loc,x,145.2
d,all,ux

nsel,s,loc,y,0
nsel,a,loc,y,138.4
d,all,uy
save

!*******
********保存3

!******************进入求解器部分

/solu
allsel
acel,,,9.8

time,1
step=1

esel,s,type,,3
ekill,all

esel,s,type,,2
ekill,all

antype,static
nsubst,1
autots,on
lnsrch,on !线性搜索
NLGEOM,on
nropt,full
cnvtol,u,,0.05,0

cutcontrol,plslimit,2 !调整最大当量塑形应变
outres,all,all
output,all,all
pred,on !预测（作用是加速收敛）
eqslv,sparse,1e-4
allsel

solve


time,2
step=2

esel,s,real,,6,9,1
ekill,all

vsel,s,loc,x,25,99.8
vsel,r,loc,y,25,113.4
vsel,r,loc,z,0,-1.8
eslv,s
ekill,all

vsel,s,loc,x,99.8,120.2
vsel,r,loc,y,76,113.4
vsel,r,loc,z,0,-1.8
eslv,s
ekill,all

esel,s,real,,5
ealive,all

esel,s,type,,2
esel,u,real,,9
ealive,all

antype,static
nsubst,1
autots,on
lnsrch,on !线性搜索
NLGEOM,on
nropt,full
cnvtol,u,,0.05,0

cutcontrol,plslimit,2 !调整最大当量塑形应变
outres,all,all
output,all,all
pred,on !预测（作用是加速收敛）
eqslv,sparse,1e-4
allsel

solve

time,3
step=3

vsel,s,loc,x,25,99.8
vsel,r,loc,y,25,113.4
vsel,r,loc,z,0,-3.4
eslv,s
ekill,all

vsel,s,loc,x,99.8,120.2
vsel,r,loc,y,76,113.4
vsel,r,loc,z,0,-3.4
eslv,s
ekill,all

esel,s,type,,2
ealive,all

esel,s,type,,3
ealive,all


antype,static
nsubst,1
autots,on
lnsrch,on !线性搜索
NLGEOM,on
nropt,full
cnvtol,u,,0.05,0

cutcontrol,plslimit,2 !调整最大当量塑形应变
outres,all,all
output,all,all
pred,on !预测（作用是加速收敛）
eqslv,sparse,1e-4
allsel

solve


time,4
step=4

vsel,s,loc,x,25,99.8
vsel,r,loc,y,25,113.4
vsel,r,loc,z,0,-7.5
eslv,s
ekill,all

vsel,s,loc,x,99.8,120.2
vsel,r,loc,y,76,113.4
vsel,r,loc,z,0,-7.5
eslv,s
ekill,all

esel,s,type,,2
ealive,all

esel,s,type,,3
ealive,all


antype,static
nsubst,1
autots,on
lnsrch,on !线性搜索
NLGEOM,on
nropt,full
cnvtol,u,,0.05,0

cutcontrol,plslimit,2 !调整最大当量塑形应变
outres,all,all
output,all,all
pred,on !预测（作用是加速收敛）
eqslv,sparse,1e-4
allsel

solve

save

!********************************通用后处理
!********************************通用后处理

!**********自重，第一工况
!******执行完time,1，计算之后

/post1

allsel !100开挖前土体

esel,s,live

plnsol,u,z,0,1 !101开挖前土体UZ

finish

!**********第一次开挖，第二工况

/solu
antype,,rest

!******执行完time,2，计算之后
/post1

allsel

esel,s,live !200第一次开挖土体

esel,r,type,,1
plnsol,u,z,0,1 !201第一次开挖土体

UZ

esel,s,type,,2 !墙

plnsol,u,x,0,1 !202第一次开挖墙UX
plnsol,u,y,0,1 !203第一次开挖墙UY

finish

!**********第二次开挖，第三工况

/solu
antype,,rest

!******执行完time,3，计算之后
/post1

allsel

esel,s,live
/eshape,1 !300第二次开挖土体
/eshape,0

esel,r,type,,1
plnsol,u,z,0,1 !301第二次开挖土
体UZ

esel,s,type,,2 !墙

plnsol,u,x,0,1 !302第二次开挖墙UX
plnsol,u,y,0,1 !303第二次开挖墙UY

esel,s,real,,6,8,1 !内支撑体系
plnsol,u,x,0,1 !304第二次开挖内撑UX
plnsol,u,y,0,1 !305第二次开挖内撑UY

etable,fxi,smisc,1
etable,fyi,smisc,2
etable,Mi,smisc,6
etable,fxj,smisc,7
etable,fyj,smisc,8
etable,Mj,smisc,12

plls,mi,mj,1,0 !306第二次开挖内撑M
plls,fxi,fxj,1,0 !307第二次开挖内撑N
plls,fyi,fyj,1,0 !308第二次开挖内撑V

esel,s,real,,8 !环撑及对撑
plnsol,u,x,0,1 !309第二次开挖环撑及对撑UX
plnsol,u,y,0,1 !310第二次开挖环撑及对撑UY

etable,fxi,smisc,1
etable,fyi,smisc,2
etable,Mi,smisc,6
etable,fxj,smisc,7
etable,fyj,smisc,8
etable,Mj,smisc,12

plls,mi,mj,1,0 !311第二次开挖环撑及对撑M
plls,fxi,fxj,1,0 !312第二次开挖环撑及对撑N
plls,fyi,fyj,1,0 !313第二次开挖环撑及对撑V

finish

!**********第三次开挖，第四工况

/solu
antype,,rest

!******执行完time,4，计算之后
/post1

allsel

esel,s,live
/eshape,1 !400第三次开挖土体
/eshape,0

esel,r,type,,1
plnsol,u,z,0,1 !401第三次开挖土体UZ

esel,s,type,,2 !墙
esel,u,real,,9

plnsol,u,x,0,1 !402第三次开挖墙UX
plnsol,u,y,0,1 !403第三次开挖墙UY

esel,s,real,,6,8,1 !内支撑体系
plnsol,u,x,0,1 !404第三次开挖内撑UX
plnsol,u,y,0,1 !405第三次开挖内撑UY

etable,fxi,smisc,1
etable,fyi,smisc,2
etable,Mi,smisc,6
etable,fxj,smisc,7
etable,fyj,smisc,8
etable,Mj,smisc,12

plls,mi,mj,1,0 !406第三次开挖内撑M
plls,fxi,fxj,1,0 !407第三次开挖内撑N
plls,fyi,fyj,1,0 !408第三次开挖内撑V



esel,s,real,,8 !环撑及对撑
plnsol,u,x,0,1 !409第三次开挖环撑及对撑UX
plnsol,u,y,0,1 !410第三次开挖环撑及对撑UY

etable,fxi,smisc,1
etable,fyi,smisc,2
etable,Mi,smisc,6
etable,fxj,smisc,7
etable,fyj,smisc,8
etable,Mj,smisc,12

plls,mi,mj,1,0 !411第三次开挖环撑及对撑M
plls,fxi,fxj,1,0 !412第三次开挖环撑及对撑N
plls,fyi,fyj,1,0 !413第三次开挖环撑及对撑V

finish


!********接下来相减,2-1,time,1和time,2,相减

/post1
!执行路径1，手动读入第二工况，第一子步

esel,s,type,,1
esel,r,live

LCDEF,1,1,1
LCDEF,2,2,1
LCOPER,SUB,1,,2

plnsol,u,z,0,1 !204第一次相减后土体UZ

esel,s,type,,2 !墙

plnsol,u,x,0,1 !204第一次相减后墙UX
plnsol,u,y,0,1 !205第一次相减后墙UY


!********接

下来相减,3-1,time,1和time,3,相减

/post1
!执行路径1，手动读入第三工况，第一子步

esel,s,type,,1
esel,r,live

LCDEF,1,1,1
LCDEF,3,3,1
LCOPER,SUB,1,,3

plnsol,u,z,0,1 !314第二次相减后土体UZ

esel,s,type,,2 !墙

plnsol,u,x,0,1 !315第二次相减后墙UX
plnsol,u,y,0,1 !316第二次相减后墙UY

!********接下来相减,4-1,time,1和time,4,相减

/post1

!执行路
径1，手动读入第四工况，第一子步

esel,s,type,,1
esel,r,live

LCDEF,1,1,1
LCDEF,4,4,1
LCOPER,SUB,1,,3

plnsol,u,z,0,1 !414第三次相减后土体UZ

esel,s,type,,2 !墙

plnsol,u,x,0,1 !315第三次相减后墙UX
plnsol,u,y,0,1 !316第三次相减后墙UY


finish