跪求|隧道及地下工程 - 领先的岩土技术社区，岩土领域的媒体、社区与应用平台！

离线yuzhenzhai

实习生

发帖: 3

土币: 1580

威望: 13

原创币: 0

只看楼主倒序阅读使用道具楼主发表于: 2011-06-01

盾构隧道施工flac3d模拟程序

分享到 淘江湖新浪 QQ微博 QQ空间开心人人豆瓣网易微博百度鲜果白社会飞信

离线wdd-888

助理工程师

发帖: 24

土币: 28

威望: 84

原创币: 0

只看该作者 1楼发表于: 2011-09-30

回楼主(yuzhenzhai) 的帖子

ansys中建模导入flac3d

/prep7
NUMSTR,KP,1
NUMSTR,LINE,1
NUMSTR,AREA,1
NUMSTR,VOLU,1
CSYS,0
k, 1  ,9.495,2.148,.000,
k, 2  ,1.508,2.148,.000,
l, 1  ,2
k, 3  ,9.405,2.398,.000,
k, 4  ,1.598,2.398,.000,
l, 3  ,4
k, 5  ,5.501,.848,.000,
k, 6  ,9.474,2.211,.000,
LARC,1  ,3  ,6
k, 7  ,7.922,4.280,.000,
k, 8  ,9.141,2.944,.000,
LARC,3  ,7  ,8
k, 9  ,2.491,3.777,.000,
k, 10  ,2.218,3.467,.000,
LARC,9  ,4  ,10
k, 11  ,1.574,2.336,.000,
LARC,4  ,2  ,11
k, 12  ,1.603,-.421,.000,
k, 13  ,9.892,.000,.000,
k, 14  ,9.540,2.004,.000,
k, 15  ,9.850,.508,.000,
LARC,13  ,14  ,15
k, 16  ,5.501,5.548,.000,
k, 17  ,1.964,-3.806,.000,
k, 18  ,2.453,-3.977,.000,
k, 19  ,2.085,-3.851,.000,
LARC,17  ,18  ,19
k, 20  ,2.602,-4.023,.000,
k, 21  ,6.612,-4.391,.000,
k, 22  ,3.586,-4.267,.000,
LARC,20  ,21  ,22
k, 23  ,1.751,-4.367,.000,
k, 24  ,2.453,-4.605,.000,
k, 25  ,1.925,-4.431,.000,
LARC,23  ,24  ,25
k, 26  ,2.602,-4.648,.000,
k, 27  ,6.733,-4.981,.000,
k, 28  ,3.618,-4.884,.000,
LARC,26  ,27  ,28
k, 29  ,2.064,4.198,.000,
k, 30  ,.973,2.439,.000,
k, 31  ,1.721,3.806,.000,
LARC,29  ,30  ,31
k, 32  ,12.823,3.231,.000,
k, 33  ,9.939,7.140,.000,
k, 34  ,12.348,4.370,.000,
LARC,32  ,33  ,34
k, 35  ,1.751,-2.802,.000,
k, 36  ,1.751,-3.723,.000,
l, 35  ,36
l, 23  ,17
k, 37  ,.824,-2.802,.000,
k, 38  ,1.449,-2.802,.000,
l, 37  ,38
l, 38  ,35
k, 39  ,9.251,-2.952,.000,
k, 40  ,9.251,-3.723,.000,
l, 39  ,40
k, 41  ,9.251,-4.367,.000,
k, 42  ,10.493,-4.367,.000,
l, 41  ,42
k, 43  ,9.039,-3.806,.000,
l, 41  ,43
k, 44  ,9.507,-2.952,.000,
l, 44  ,39
l, 35  ,36
l, 23  ,17
k, 45  ,.000,11.948,.000,
k, 46  ,13.330,11.948,.000,
l, 45  ,46
k, 47  ,17.000,11.948,.000,
k, 48  ,33.000,11.948,.000,
l, 47  ,48
k, 49  ,33.000,-10.052,.000,
k, 50  ,33.000,-27.052,.000,
l, 49  ,50
k, 51  ,.000,-10.052,.000,
k, 52  ,.000,-27.052,.000,
l, 51  ,52
k, 53  ,17.000,-27.052,.000,
l, 50  ,53
k, 54  ,17.000,-10.052,.000,
l, 49  ,54
k, 55  ,13.393,-10.052,.000,
l, 54  ,55
k, 56  ,9.251,-10.052,.000,
l, 55  ,56
k, 57  ,1.751,-10.052,.000,
l, 56  ,57
l, 57  ,51
k, 58  ,33.000,.000,.000,
l, 48  ,58
k, 59  ,17.000,.000,.000,
l, 47  ,59
l, 54  ,53
k, 60  ,13.393,-27.052,.000,
l, 53  ,60
k, 61  ,9.251,-27.052,.000,
l, 60  ,61
k, 62  ,1.751,-27.052,.000,
l, 61  ,62
l, 62  ,52
k, 63  ,.000,6.235,.000,
l, 45  ,63
k, 64  ,.000,2.603,.000,
l, 63  ,64
k, 65  ,.000,-2.802,.000,
k, 66  ,.000,-3.723,.000,
l, 65  ,66
k, 67  ,2.602,.000,.000,
k, 68  ,6.993,.000,.000,
l, 67  ,68
k, 69  ,10.493,.000,.000,
l, 13  ,69
k, 70  ,10.743,.000,.000,
l, 69  ,70
k, 71  ,13.393,.000,.000,
l, 71  ,59
k, 72  ,.193,.371,.000,
k, 73  ,1.382,2.234,.000,
LARC,2  ,30  ,73
k, 74  ,.121,.088,.000,
k, 75  ,2.545,.625,.000,
LARC,67  ,2  ,75
k, 76  ,10.493,-2.952,.000,
l, 69  ,76
k, 77  ,10.493,-3.723,.000,
l, 76  ,77
l, 77  ,42
k, 78  ,13.393,-4.367,.000,
l, 71  ,78
l, 47  ,32
l, 78  ,55
l, 55  ,60
l, 41  ,56
l, 56  ,61
l, 23  ,57
l, 57  ,62
k, 79  ,2.602,-2.819,.000,
l, 67  ,79
l, 79  ,20
l, 20  ,26
l, 59  ,58
k, 80  ,2.453,-2.816,.000,
l, 35  ,80
k, 81  ,6.488,-2.897,.000,
l, 79  ,81
l, 37  ,65
l, 36  ,66
l, 44  ,76
l, 40  ,77
k, 82  ,1.719,-.392,.000,
k, 83  ,13.336,.821,.000,
LARC,71  ,32  ,83
l, 9  ,29
l, 46  ,47
l, 46  ,33
k, 84  ,8.268,4.770,.000,
l, 84  ,7
k, 85  ,55.000,11.948,.000,
l, 48  ,85
k, 86  ,.000,-49.052,.000,
l, 52  ,86
k, 87  ,55.000,-10.052,.000,
l, 49  ,87
k, 88  ,55.000,-27.052,.000,
l, 50  ,88
k, 89  ,33.000,-49.052,.000,
l, 50  ,89
k, 90  ,17.000,-49.052,.000,
l, 53  ,90
l, 87  ,88
k, 91  ,55.000,-49.052,.000,
l, 88  ,91
l, 86  ,90
l, 90  ,89
l, 89  ,91
l, 59  ,54
l, 58  ,49
k, 92  ,55.000,.000,.000,
l, 58  ,92
l, 85  ,92
l, 92  ,87
k, 93  ,9.400,-.421,.000,
k, 94  ,.910,2.249,.000,
k, 95  ,.615,1.004,.000,
LARC,94  ,37  ,95
k, 96  ,2.453,.000,.000,
l, 67  ,96
k, 97  ,1.446,1.950,.000,
k, 98  ,1.187,.774,.000,
LARC,97  ,38  ,98
k, 99  ,.146,.251,.000,
k, 100  ,1.320,2.039,.000,
LARC,97  ,94  ,100
k, 101  ,.171,.056,.000,
k, 102  ,2.396,.568,.000,
LARC,96  ,97  ,102
l, 96  ,80
l, 18  ,24
l, 80  ,18
k, 103  ,.129,.405,.000,
k, 104  ,.736,2.522,.000,
k, 105  ,.915,2.462,.000,
LARC,30  ,104  ,105
l, 30  ,94
k, 106  ,.000,2.387,.000,
l, 64  ,106
l, 80  ,79
k, 107  ,2.491,-3.988,.000,
LARC,18  ,20  ,107
l, 106  ,65
k, 108  ,7.494,8.285,.000,
LARC,33  ,63  ,108
k, 109  ,6.741,5.485,.000,
LARC,84  ,29  ,109
k, 110  ,6.585,4.905,.000,
LARC,7  ,9  ,110
k, 111  ,.000,-4.367,.000,
l, 66  ,111
l, 23  ,111
l, 111  ,51
l, 2  ,97
l, 17  ,36
l, 24  ,26
l, 40  ,43
k, 112  ,10.743,-4.367,.000,
l, 70  ,112
k, 113  ,8.412,4.975,.000,
k, 114  ,1.887,4.374,.000,
k, 115  ,6.806,5.726,.000,
LARC,113  ,114  ,115
l, 104  ,64
k, 116  ,.675,2.339,.000,
l, 116  ,106
l, 78  ,112
l, 112  ,42
l, 70  ,71
k, 117  ,10.303,2.411,.000,
l, 32  ,117
l, 33  ,113
l, 113  ,84
l, 63  ,114
l, 114  ,29
k, 118  ,5.505,.844,.000,
k, 119  ,1.525,3.961,.000,
LARC,114  ,104  ,119
l, 104  ,116
l, 94  ,116
k, 120  ,9.703,-2.230,.000,
LARC,44  ,13  ,120
k, 121  ,4.496,.106,.000,
k, 122  ,10.357,2.272,.000,
k, 123  ,10.726,.579,.000,
LARC,70  ,122  ,123
k, 124  ,11.542,-2.768,.000,
k, 125  ,10.060,2.347,.000,
k, 126  ,10.242,2.396,.000,
LARC,117  ,125  ,126
k, 127  ,15.216,-2.927,.000,
k, 128  ,6.775,-.706,.000,
LARC,68  ,81  ,128
k, 129  ,7.170,.000,.000,
l, 68  ,129
l, 129  ,13
k, 130  ,6.638,-2.900,.000,
l, 81  ,130
l, 130  ,39
k, 131  ,6.761,-4.373,.000,
k, 132  ,7.342,-4.282,.000,
LARC,131  ,43  ,132
k, 133  ,6.649,-4.386,.000,
LARC,21  ,131  ,133
k, 134  ,6.881,-4.962,.000,
k, 135  ,7.486,-4.865,.000,
LARC,134  ,41  ,135
k, 136  ,10.344,2.307,.000,
LARC,122  ,117  ,136
k, 137  ,9.959,3.138,.000,
LARC,117  ,113  ,137
k, 138  ,4.500,.109,.000,
k, 139  ,9.732,3.033,.000,
LARC,125  ,84  ,139
k, 140  ,15.053,-2.948,.000,
k, 141  ,6.637,-2.900,.000,
k, 142  ,6.941,-.705,.000,
LARC,129  ,141  ,142
k, 143  ,6.758,-4.373,.000,
k, 144  ,6.785,-4.521,.000,
LARC,143  ,134  ,144
k, 145  ,6.638,-4.539,.000,
LARC,21  ,27  ,145
k, 146  ,6.495,-3.272,.000,
LARC,81  ,21  ,146
k, 147  ,6.643,-3.270,.000,
LARC,141  ,143  ,147
k, 148  ,8.750,1.766,.000,
LARC,1  ,68  ,148
k, 149  ,11.567,-2.798,.000,
k, 150  ,8.840,1.644,.000,
LARC,14  ,129  ,150
k, 151  ,9.917,2.303,.000,
LARC,125  ,1  ,151
l, 134  ,27
l, 1  ,14
k, 152  ,10.114,2.208,.000,
k, 153  ,10.476,.563,.000,
LARC,69  ,152  ,153
l, 152  ,122
k, 154  ,9.969,2.163,.000,
LARC,152  ,14  ,154
l, 125  ,152

nummrg,all,,,,low  ! 合并
numcap,all         !压缩

lsel,s,line,,80,84,1     ！生成面
lsel,a,line,,88,89,1
lsel,a,line,,22,23,1
lsel,a,line,,74,75,1
lsel,a,line,,39,41,1
lsel,a,line,,100,103,3
lsel,a,line,,107,109,2
lsel,a,line,,25,71,46
al,all
allsel,all
asbl,1a,all,,delete,keep
numcmp,all

!********定义单元类型******四级围岩
ET,1,PLANE42
ET,2,solid45
!******定义材料模型*******

Type,1
MSHKEY,2
MSHAPE,0,2D
lsel,s,line,,22,104,82,0         !划分网格(第一行)
lesize,all,,,10
lsel,s,line,,39,72,33,0
lesize,all,,,5,1.2
amesh,13
allsel

lsel,s,line,,115
lesize,all,,,10
lsel,s,line,,122,124,2,0
lesize,all,,,3,1.2
amesh,41
allsel

lsel,s,line,,105
lesize,all,,,10
lsel,s,line,,123,125,2,0
lesize,all,,,2
amesh,42
allsel

lsel,s,line,,106
lesize,all,,,10
lsel,s,line,,70,73,3,0
lesize,all,,,2
amesh,36
allsel

lsel,s,line,,116
lesize,all,,,3,1.2
lsel,s,line,,40,126,86
lesize,all,,,4
amesh,26
allsel

lsel,s,line,,12
lesize,all,,,4,
amesh,9
allsel

lsel,s,line,,117
lesize,all,,,3,1.2
lsel,s,line,,100,127,27
lesize,all,,,2
amesh,40
allsel

lsel,s,line,,128
lesize,all,,,2
lsel,s,line,,99
lesize,all,,,2
amesh,39
allsel

lsel,s,line,,93
lesize,all,,,2
amesh,30
allsel

lsel,s,line,,13
lesize,all,,,7
amesh,10
allsel
amesh,11
allsel
amesh,43
allsel

lsel,s,line,,33
lesize,all,,,8
amesh,25
allsel

lsel,s,line,,140
lesize,all,,,2
amesh,44
allsel

lsel,s,,,130,140,10
LCCAT,all
ESIZE,5,0,
ACLEAR,59
AMAP,59,117,32,71,70
allsel
amesh,29
allsel

lsel,s,line,,150
lesize,all,,,2
amesh,48
allsel
amesh,5
allsel

lsel,s,line,,1,2,1
lesize,all,,,17
lsel,s,,,3,6,3
lesize,all,,,2
amesh,2
allsel

lsel,s,,,3,4,1
LCCAT,all
ESIZE,5,0,
ACLEAR,49
AMAP,49,84,7,1,125
allsel

lsel,s,,,5,6,1
LCCAT,all
ESIZE,5,0,
ACLEAR,50
AMAP,50,2,30,29,9
allsel
amesh,4
allsel

lsel,s,line,,148
lesize,all,,,6
amesh,3
allsel

lsel,s,line,,103
lesize,all,,,18
lsel,s,,,117,128,11
LCCAT,all
ESIZE,5,0,
ACLEAR,58
AMAP,58,106,94,37,65
allsel
amesh,12
allsel
amesh,31
allsel
lsel,s,line,,95
lesize,all,,,10
amesh,33
allsel

lsel,s,,,17,63,46
LCCAT,all
ESIZE,5,0,
ACLEAR,53
AMAP,53,97,96,80,38
allsel
amesh,46
allsel

lsel,s,line,,134
lesize,all,,,6
amesh,1
allsel
amesh,27
allsel

LDELE,135
lcomb,162,136
lcomb,147,163
LDELE,147
LDELE,163
LDELE,138
LDELE,144
lcomb,137,161
lcomb,135,164
LDELE,143
amesh,45
allsel

lsel,s,,,21,136,115
LCCAT,all
ESIZE,5,0,
ACLEAR,62
AMAP,62,44,13,129,155
allsel
amesh,28
allsel

lsel,s,line,,41
lesize,all,,,2
lsel,s,,,16,17,1
lsel,a,,,65
LCCAT,all
ESIZE,5,0,
ACLEAR,60
AMAP,60,65,66,36,35
allsel

lsel,s,,,8,111,103
LCCAT,all
ESIZE,5,0,
ACLEAR,51
AMAP,51,35,36,18,80
allsel

lsel,s,line,,107
lesize,all,,,2
lsel,s,,,66,111,45
LCCAT,all
ESIZE,5,0,
ACLEAR,52
AMAP,52,66,111,23,17
allsel
amesh,6
allsel
amesh,34
allsel
amesh,35
allsel
amesh,8
allsel
amesh,7
allsel
amesh,64
allsel
amesh,47
allsel

lsel,s,,,113,137,24
LCCAT,all
ESIZE,5,0,
ACLEAR,67
AMAP,67,155,156,40,39
allsel
amesh,56
allsel

lsel,s,,,21,67,46
LCCAT,all
ESIZE,5,0,
ACLEAR,54
AMAP,54,40,39,76,77
allsel

lsel,s,,,68,113,45
LCCAT,all
ESIZE,5,0,
ACLEAR,55
AMAP,55,41,43,77,42
allsel

lsel,s,,,48,50,1
LCCAT,all
ESIZE,5,0,
ACLEAR,63
AMAP,63,112,42,69,70
allsel
amesh,32
allsel

lsel,s,line,,109
lesize,all,,,6
amesh,23
allsel

lsel,s,,,10,11,1
lsel,a,,,112,139,27
lsel,a,,,151
LCCAT,all
ESIZE,5,0,
ACLEAR,66
AMAP,66,57,56,41,23
allsel

lsel,s,,,118,119,1
lsel,a,,,19
LCCAT,all
ESIZE,5,0,
ACLEAR,61
AMAP,61,55,56,41,78
allsel

lsel,s,,,51,53,2
LCCAT,all
ESIZE,5,0,
ACLEAR,57
AMAP,57,54,55,71,59
allsel

lsel,s,line,,23
lesize,all,,,6,1.2
amesh,14
allsel
amesh,19
allsel

lsel,s,line,,74
lesize,all,,,7,1.2
amesh,24
allsel
amesh,37
allsel

lsel,s,line,,25     !第二行
lesize,all,,,6
amesh,17
allsel
amesh,22
allsel
amesh,21
allsel
amesh,20
allsel
amesh,15
allsel
amesh,16
allsel

lsel,s,line,,81    !第三行
lesize,all,,,7
amesh,38
allsel
amesh,18
allsel

lsel,s,,,35,38,1
LCCAT,all
AMAP,65,86,90,53,52
allsel

lsel,r,lcca    !删除前面连接的线元素
ldele,all

arsym,x,all,,0,0     !对称

nummrg,all,,,,low  ! 合并
numcap,all         !压缩

k,1911,0,0,0
k,1912,0,0,1
l,1911,1912

MP,EX,1,1e9
MP,PRXY,1,0.4
MP,DENS,1,1800

MP,EX,2,2e9
MP,PRXY,2,0.4
MP,DENS,2,1800

MP,EX,3,3e9
MP,PRXY,3,0.4
MP,DENS,3,1800

MP,EX,4,4e9
MP,PRXY,4,0.4
MP,DENS,4,1800

MP,EX,5,5e9
MP,PRXY,5,0.4
MP,DENS,5,1800

MP,EX,6,6e9
MP,PRXY,6,0.4
MP,DENS,6,1800

MP,EX,7,7e9
MP,PRXY,7,0.4
MP,DENS,7,1800

MP,EX,8,8e9
MP,PRXY,8,0.4
MP,DENS,8,1800

MP,EX,9,9e9
MP,PRXY,9,0.4
MP,DENS,9,1800

MP,EX,10,10e9
MP,PRXY,10,0.4
MP,DENS,10,1800

MP,EX,11,1.1e10
MP,PRXY,11,0.4
MP,DENS,11,1800

MP,EX,12,1.2e10
MP,PRXY,12,0.4
MP,DENS,12,1800

MP,EX,13,1.3e10
MP,PRXY,13,0.4
MP,DENS,13,1800

MP,EX,14,1.4e10
MP,PRXY,14,0.4
MP,DENS,14,1800

MP,EX,15,1.5e10
MP,PRXY,15,0.4
MP,DENS,15,1800

MP,EX,16,1.6e10
MP,PRXY,16,0.4
MP,DENS,16,1800

MP,EX,17,1.7e10
MP,PRXY,17,0.4
MP,DENS,17,1800

MP,EX,18,1.8e10
MP,PRXY,18,0.4
MP,DENS,18,1800

MP,EX,19,1.9e10
MP,PRXY,19,0.4
MP,DENS,19,1800

MP,EX,20,2e10
MP,PRXY,20,0.4
MP,DENS,20,1800

MP,EX,21,2.1e10
MP,PRXY,21,0.4
MP,DENS,21,1800

MP,EX,22,2.2e10
MP,PRXY,22,0.4
MP,DENS,22,1800

MP,EX,23,2.3e10
MP,PRXY,23,0.4
MP,DENS,23,1800

MP,EX,24,2.4e10
MP,PRXY,24,0.4
MP,DENS,24,1800

MP,EX,25,2.5e10
MP,PRXY,25,0.4
MP,DENS,25,1800

MP,EX,26,2.6e10
MP,PRXY,26,0.4
MP,DENS,26,1800

MP,EX,27,2.7e10
MP,PRXY,27,0.4
MP,DENS,27,1800

MP,EX,28,2.8e10
MP,PRXY,28,0.4
MP,DENS,28,1800

MP,EX,29,2.9e10
MP,PRXY,29,0.4
MP,DENS,29,1800

MP,EX,30,3.0e10
MP,PRXY,30,0.4
MP,DENS,30,1800

MP,EX,31,3.1e10
MP,PRXY,31,0.4
MP,DENS,31,1800

MP,EX,32,3.2e10
MP,PRXY,32,0.4
MP,DENS,32,1800

MP,EX,33,3.3e10
MP,PRXY,33,0.4
MP,DENS,33,1800

MP,EX,34,3.4e10
MP,PRXY,34,0.4
MP,DENS,34,1800

MP,EX,35,3.5e10
MP,PRXY,35,0.4
MP,DENS,35,1800

MP,EX,36,3.6e10
MP,PRXY,36,0.4
MP,DENS,36,1800

MP,EX,37,3.7e10
MP,PRXY,37,0.4
MP,DENS,37,1800

MP,EX,38,3.8e10
MP,PRXY,38,0.4
MP,DENS,38,1800

MP,EX,39,3.9e10
MP,PRXY,39,0.4
MP,DENS,39,1800

MP,EX,40,4e10
MP,PRXY,40,0.4
MP,DENS,40,1800

MP,EX,41,4.1e10
MP,PRXY,41,0.4
MP,DENS,41,1800

MP,EX,42,4.2e10
MP,PRXY,42,0.4
MP,DENS,42,1800

type,2   !拉伸
EXTOPT,ESIZE,1,0,   !中到洞顶，面生成体单元的控制命令
EXTOPT,ACLEAR,1
mat,1
asel,s,area,,106,107,1
asel,a,area,,39,40,1
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,    !中到洞右上-jiao
EXTOPT,ACLEAR,1
mat,2
asel,s,area,,30
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,  !中到洞右上
EXTOPT,ACLEAR,1
mat,3
asel,s,area,,31
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,    !中到洞右下
EXTOPT,ACLEAR,1
mat,4
asel,s,area,,33,34,1
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,    !中到洞右下-jiao
EXTOPT,ACLEAR,1
mat,5
asel,s,area,,35
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !中到洞左上-jiao
EXTOPT,ACLEAR,1
mat,6
asel,s,area,,97
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,  !中到洞左上
EXTOPT,ACLEAR,1
mat,7
asel,s,area,,98
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !中到洞左下
EXTOPT,ACLEAR,1
mat,8
asel,s,area,,100,101,1
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,  !中到洞左下-jiao
EXTOPT,ACLEAR,1
mat,9
asel,s,area,,102
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !中隔墙
EXTOPT,ACLEAR,1
mat,10
asel,s,area,,52,58,6
asel,a,area,,60,119,59
asel,a,area,,125,127,2
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !中洞左二衬
EXTOPT,ACLEAR,1
mat,11
asel,s,area,,73,79,6
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !中洞右二衬
EXTOPT,ACLEAR,1
mat,12
asel,s,area,,6,12,6
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !中洞内zuo围眼
EXTOPT,ACLEAR,1
mat,13
asel,s,area,,118,120,2
VDRAG,all,,,,,,304
allsel

!右到洞

type,2
EXTOPT,ESIZE,1,0,   !边一衬
EXTOPT,ACLEAR,1
mat,14
asel,s,area,,29
asel,a,area,,63
asel,a,area,,44
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !上一衬-交
EXTOPT,ACLEAR,1
mat,15
asel,s,area,,48
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !上一衬
EXTOPT,ACLEAR,1
mat,16
asel,s,area,,46
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !下一衬
EXTOPT,ACLEAR,1
mat,17
asel,s,area,,45,64,19
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !下一衬-交
EXTOPT,ACLEAR,1
mat,18
asel,s,area,,47
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !二衬
EXTOPT,ACLEAR,1
mat,19
asel,s,area,,54,56,1
asel,a,area,,28
asel,a,area,,5
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !围眼
EXTOPT,ACLEAR,1
mat,20
asel,s,area,,1
asel,a,area,,62,67,5
VDRAG,all,,,,,,304
allsel

!左导洞

type,2
EXTOPT,ESIZE,1,0,   !边一衬
EXTOPT,ACLEAR,1
mat,21
asel,s,area,,96
asel,a,area,,130
asel,a,area,,111
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !上一衬-交
EXTOPT,ACLEAR,1
mat,22
asel,s,area,,115
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !上一衬
EXTOPT,ACLEAR,1
mat,23
asel,s,area,,113
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !下一衬
EXTOPT,ACLEAR,1
mat,24
asel,s,area,,112,131,19
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !下一衬-交
EXTOPT,ACLEAR,1
mat,25
asel,s,area,,114
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !二衬
EXTOPT,ACLEAR,1
mat,26
asel,s,area,,121,123,1
asel,a,area,,72,95,23
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !围岩
EXTOPT,ACLEAR,1
mat,27
asel,s,area,,68
asel,a,area,,129,134,5
VDRAG,all,,,,,,304
allsel

!右洞

type,2
EXTOPT,ESIZE,1,0,   !右上一衬
EXTOPT,ACLEAR,1
mat,28
asel,s,area,,9
asel,a,area,,42,43,1
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !右上二衬
EXTOPT,ACLEAR,1
mat,29
asel,s,area,,36
asel,a,area,,49,50,1
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !道班
EXTOPT,ACLEAR,1
mat,30
asel,s,area,,2
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !右上围岩
EXTOPT,ACLEAR,1
mat,31
asel,s,area,,3,4,1
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !右下围岩
EXTOPT,ACLEAR,1
mat,32
asel,s,area,,8,27,19
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !右下二衬
EXTOPT,ACLEAR,1
mat,33
asel,s,area,,7
VDRAG,all,,,,,,304
allsel

!左洞

type,2
EXTOPT,ESIZE,1,0,   !左上一衬
EXTOPT,ACLEAR,1
mat,34
asel,s,area,,76
asel,a,area,,109,110,1
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !左上二衬
EXTOPT,ACLEAR,1
mat,35
asel,s,area,,103
asel,a,area,,116,117,1
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !道班
EXTOPT,ACLEAR,1
mat,36
asel,s,area,,69
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !左上围岩
EXTOPT,ACLEAR,1
mat,37
asel,s,area,,70,71,1
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !左下围岩
EXTOPT,ACLEAR,1
mat,38
asel,s,area,,75,94,19
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !左下二衬
EXTOPT,ACLEAR,1
mat,39
asel,s,area,,74
VDRAG,all,,,,,,304
allsel

!外围

type,2
EXTOPT,ESIZE,1,0,   !加固区
EXTOPT,ACLEAR,1
mat,40
asel,s,area,,11,26,15
asel,a,area,,41,59,18
asel,a,area,,78,93,15
asel,a,area,,108,126,18
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !外围围岩
EXTOPT,ACLEAR,1
mat,41
asel,s,area,,13,25,1
asel,a,area,,65,66,1
asel,a,area,,37,38,1
asel,a,area,,10,32,22
asel,a,area,,57,61,4
asel,a,area,,80,92,1
asel,a,area,,104,105,1
asel,a,area,,132,133,1
asel,a,area,,77,99,22
asel,a,area,,124,128,4
VDRAG,all,,,,,,304
allsel

type,2
EXTOPT,ESIZE,1,0,   !中洞内you围眼
EXTOPT,ACLEAR,1
mat,42
asel,s,area,,51,53,2
VDRAG,all,,,,,,304

k,1000,5,0,1
CSKP,20,0,1912,1000,1911,1,1,  !定义一个局部坐标，三点式
DSYS,20,                       !激活局部坐标

nummrg,all,,,,low  ! 合并
numcap,all         !压缩

flac 3d模拟隧道开挖命令流

def moduli_data
s_mod=y_mod/(2.0*(1.0+p_ratio))
b_mod=y_mod/(3.0*(1.0-2.0*p_ratio))
end
set y_mod 5e14
set p_ratio 0.25
moduli_data

mo el ;定义本构模型
prop bulk b_mod shear s_mod

fix x range x  -55.1 -54.9
fix x range x  54.9  55.1
fix z range z  -49.152 -48.952
fix y range y  -0.1 0.1
fix y range y  0.9 1.1

ini density 2479
set gravity 0 0 -10
set large
apply szz -0.115e6 range z 11.848 12.048
;释放0
solve
save 地应力.sav

;围岩
def moduli_data11
s_mod=y_mod/(2.0*(1.0+p_ratio))
b_mod=y_mod/(3.0*(1.0-2.0*p_ratio))
end
set y_mod 6e9
set p_ratio 0.25
moduli_data11

mo mohr   ;所有的
prop bulk b_mod shear s_mod fric 39 coh 7e5 ten 1e5
ini xdis=0 ydis=0 zdis=0
set large
hist unbal

;开挖中导洞
model null range group 1
model null range group 2
model null range group 3
model null range group 4
model null range group 5
model null range group 6
model null range group 7
model null range group 8
model null range group 9
model null range group 10
model null range group 11
model null range group 12
model null range group 13
model null range group 42
step 100

;zhong导洞初期支护
mo el range group 1  ;必须说明范围
mo el range group 2
mo el range group 3
mo el range group 4
mo el range group 5
mo el range group 6
mo el range group 7
mo el range group 8
mo el range group 9
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 1
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 2
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 3
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 4
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 5
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 6
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 7
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 8
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 9
step 1000
save 中导洞初期支护.sav

;施作中隔墙
mo el range group 10
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 10
step 500
save 中隔墙.sav

;施作中导洞二衬
mo el range group 11
mo el range group 12
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 11
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 12
step 500
save 中导洞二衬.sav

;开挖you导洞
model null range group 14
model null range group 15
model null range group 16
model null range group 17
model null range group 18
model null range group 19
model null range group 20
step 100

;you导洞初期支护
mo el range group 14  ;必须说明范围
mo el range group 15
mo el range group 16
mo el range group 17
mo el range group 18
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 14
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 15
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 16
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 17
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 18
step 1000
save 右导洞初期支护.sav

;施作you-bian墙
mo el range group 19
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 19
step 1000
save 右导洞边墙.sav

;开挖zuo导洞
model null range group 21
model null range group 22
model null range group 23
model null range group 24
model null range group 25
model null range group 26
model null range group 27
step 100

;zuo导洞初期支护
mo el range group 21 ;必须说明范围
mo el range group 22
mo el range group 23
mo el range group 24
mo el range group 25
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 21
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 22
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 23
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 24
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 25
step 1000
save 左导洞初期支护.sav

;施作zuo-bian墙
mo el range group 26
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 26
step 1000
save 左导洞边墙.sav

;回填
def moduli_data112
s_mod=y_mod/(2.0*(1.0+p_ratio))
b_mod=y_mod/(3.0*(1.0-2.0*p_ratio))
end
set y_mod 6e9
set p_ratio 0.25
moduli_data112

mo mohr range group 42  ;所有的
prop bulk b_mod shear s_mod fric 39 coh 7e5 ten 1e5 range group 42
step 100

;zuo上开挖
model null range group 34
model null range group 35
model null range group 36
model null range group 37
model null range group 6
model null range group 7
model null range group 22
model null range group 23
step 100

;zuo上初期支护
mo el range group 34  ;必须说明范围
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 34
step 1000

;施作zuo-shang二次衬砌
mo el range group 6
mo el range group 22
mo el range group 35
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 6
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 22
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 35
step 1000
save 左上二衬.sav

;右上开挖
model null range group 28
model null range group 29
model null range group 30
model null range group 31
model null range group 2
model null range group 3
model null range group 15
model null range group 16
model null range group 42
step 100

;右上初期支护
mo el range group 28  ;必须说明范围
pro bulk 1.6e10 shear 1.2e10 dens 2.2e3 range group 28
step 1000

;施作you-上二次衬砌
mo el range group 2
mo el range group 15
mo el range group 29
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 2
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 15
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 29
step 1000
save 右上二衬.sav

;zuo下开挖
model null range group 38
model null range group 39
model null range group 8
model null range group 9
model null range group 24
model null range group 25
step 100

;施作zuo-xia二次衬砌
mo el range group 9
mo el range group 25
mo el range group 36
mo el range group 39
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 9
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 25
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 36
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 39
step 1000
save 左下二衬.sav

;右下开挖
model null range group 32
model null range group 33
model null range group 4
model null range group 5
model null range group 17
model null range group 18
step 100

;施作you-下二次衬砌
mo el range group 5
mo el range group 18
mo el range group 30
mo el range group 33
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 5
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 18
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 30
pro bulk 2.24e10 shear 1.68e10 dens 2.3e3 range group 33
solve
save 右二衬.sav

apply szz -0.005e6 range group 30 x 3.5013 7.5013 z 2.2976 2.4976
apply szz -0.005e6 range group 36 x -7.5013 -3.5013 z 2.2976 2.4976
solve
save 面力.sav

wdd

离线wdd-888

助理工程师

发帖: 24

土币: 28

威望: 84

原创币: 0

只看该作者 2楼发表于: 2011-10-12

new
config fluid ;设置流体场
;set fluid implict on  ;设置隐式算法，只适合饱和的渗流情况，即饱和度始终为1.,详见p222
;定义模型参数，左右边界取dx_s=35m,开挖方向取dy_s=60m,上取至地表dz_ss=ms+dz/2,下取隧道底部向下30m.
def para_set
  sw=2.00         ;河水水位，位于地面以上多少m
  ms=23.00        ;埋深，地表至隧道顶部距离m
  dd=6.00         ;盾构直径
  rr=3.00         ;盾构半径
  dx_bj=35.00         ;x边界距离
  dy_bj=60.00         ;y边界距离
  dz_sbj=ms+dd/2      ;z上边界，从原点算起(埋深+盾构半径)
  dz_xbj=-(30.00+dd/2)  ;z下边界，从原点算起(30m+盾构半径)
  hf_1=9.00              ;为了便于建模，将模型从9m处划分开
  hf_11=-9.00            ;为了便于建模，将模型从-9m处划分开
  y_size=24            ;开挖方向分割为24段
  r_size=3            ;隧道半径分割为3段,n1
  rh_size=8           ;隧道四分之一弧上划分为8段,n3
  rxy_size=4          ;1、2块radcyl中的n4.
;设置流体参数
  fl_mod1=2e8         ;设置流体模量
  fl_dens1=1000       ;设置流体密度
  fl_perm1=5e-7       ;设置渗透系数
;设置摩尔准则参数
  soi_dens1=1.55e3    ;设置土体干密度
  soi_bulk1=25e6      ;设置土体体积模量
  soi_shear1=12e6     ;设置土体剪切模量
  soi_fric1=35
  soi_coh1=2e3
  soi_tens1=1e4
end
para_set

;划分网格
;生成第1块
gen zone radcyl p0 0 0 0 p1 hf_1 0 0 p2 0 dy_bj 0 p3 0 0 hf_1 dim rr rr rr rr size r_size y_size rh_size rxy_size ratio 1 1 1 1 fill group tunnel
gen zone reflect normal 0 0 -1  ;映射网格，生成第2块
;生成第3块
gen zone brick p0 hf_1 0 hf_11 p1 dx_bj 0 hf_11 p2 hf_1 dy_bj hf_11 p3 hf_1 0 hf_1 size 10 y_size 8 ratio 1.1 1 1
;生成第4块
gen zone brick p0 0 0 hf_1 p1 hf_1 0 hf_1 p2 0 dy_bj hf_1 p3 0 0 dz_sbj size 4 y_size 6 ratio 1 1 1.2
;生成第5块
gen zone brick p0 0 0 dz_xbj p1 hf_1 0 dz_xbj p2 0 dy_bj dz_xbj p3 0 0 hf_11 size 4 y_size 8 ratio 1 1 0.83
;生成第6块
gen zone brick p0 hf_1 0 hf_1 p1 dx_bj 0 hf_1 p2 hf_1 dy_bj hf_1 p3 hf_1 0 dz_sbj size 10 y_size 6 ratio 1.1 1 1.2
;生成第7块
gen zone brick p0 hf_1 0 dz_xbj p1 dx_bj 0 dz_xbj p2 hf_1 dy_bj dz_xbj p3 hf_1 0 hf_11 size 10 y_size 8 ratio 1.1 1 0.83
group tunnel1 range cyl end1 0 0 0 end2 0 30 0 rad 3.0
group soil range group tunnel1 not  ;定义其他土体为soil组
plot set rot 20 0 30
plot set mag 1.2
plot bl gr
plot add axes yellow
;加位移边界条件
fix x range x -0.1 0.1
fix x range x 34.9 35.1
FIX y range y -0.1 0.1
fix y range y 59.9 60.1
fix z range z -32.9 -33.1
;plot gpfix red sk   ;查看模型边界情况
;摩尔准则
model mohr
ini density soi_dens1
prop bulk=soi_bulk1 shear=soi_shear1 fric=soi_fric1 coh=soi_coh1 tension=soi_tens1
;流体参数
model fl_iso
prop perm fl_perm1
ini fdens fl_dens1
ini fmod fl_mod1
ini sat 1.0
model fl_null range group tunnel1
set grav 0 0 -10
;初始定应力和孔压设置,侧压力系数k0=1-sinQ=1-sin35
ini szz -403000 grad 0 0 1.55e4  ;ini szz (-(ms+rr)*10*soi_dens1) grad 0 0 10*soi_dens1
ini szz add -137e3 grad 0 0 0.45e4      ;10000*n*26+2*10000
ini sxx -171848.70 grad 0 0 0.661e4     ;1.55e4*k0,-171848.7=-1.55e4*k0*(26米)
ini sxx add -137e3 grad 0 0 0.45e4
ini syy -171848.70 grad 0 0 0.661e4
ini syy add  -137e3 grad 0 0 0.45e4
ini pp 280000 grad 0 0 -10000
app nstress -20e3 range z 26
fix pp 20000 range z 25.9 26.1
;记录点
hist unbal
hist id=2 gp ydis 0 30 0
hist id=3 zone syy 0 30 0
hist id=4 zone szz 0 30 0
;迭代初始平衡
set fl off mech on
solve
ini xdis=0 ydis=0 zdis=0
save initial_equivalence.sav
hist fltime
save initial_equivalence-fltime.sav
;开挖30m
model null range group tunnel1
;添加管片
sel shell id=1 range cyl end1 0 0 0 end2 0 30 0 rad 3.0 group tunnel1 not
sel node local xdir 1 0 0 ydir 0 -1 0 range x -0.1 0.1
sel node fix lsys range x -0.1 0.1
sel node fix x yr zr range x -0.1 0.1
sel node local xdir 0 0 -1 ydir 0 -1 0 range y -0.1 0.1
sel node fix lsys range y -0.1 0.1
sel node fix y xr zr range y -0.1 0.1
sel shell id=1 prop iso=(3.45e10,0.3) thick=0.35 density 2450
set large   ;设置大变形
;---------------------
def get_gp_maxydispyong  ;收索最大水平位移dxmax,y方向的最大位移，即为最大水平位移
   p_gp=gp_head
   dxmax=0.0        ;设置一个为0的初值，获得的节点y方向水平位移在支护压力《=1的情况，一般<0。
   loop while p_gp # null
      dz3=gp_ydis(p_gp)
         if dxmax>=dz3 then
             dxmax=dz3
             xx=gp_xpos(p_gp)
             yx=gp_ypos(p_gp)
             zx=gp_xpos(p_gp)
             id1=gp_id(p_gp)
          endif
      p_gp=gp_next(p_gp)
   endloop
end
;------------------------
;def sup_stress  ;定义开挖面支护压力参数
;     o_press=-308848.70         ;(-171848.70+-137e3),ini syy -171848.70 grad 0 0 0.661e4
;      o_grad=11110.0           ;(0.661e4+0.45e4),ini syy add  -137e3 grad 0 0 0.45e4
;  loop lumda (10.0,1.0)
;      lumda1=lumda/10.0
;      s_press=o_press*lumda1
;      s_grad=o_grad*lumda1
;  endloop
;end
;循环流固耦合计算，开挖+渗流
def excav    ;一步开挖到30m
       n_step1=10
       o_press=-308848.70         ;(-171848.70+-137e3),ini syy -171848.70 grad 0 0 0.661e4
       o_grad=11110.0           ;(0.661e4+0.45e4),ini syy add  -137e3 grad 0 0 0.45e4
     loop lumda(1,n_step1)
           lumda1=(11.0-lumda)/10.0
           s_press=o_press*lumda1
           s_grad=o_grad*lumda1
        command
           apply nstress s_press grad 0 0 s_grad range cyl end1 0 29.99 0 end2 0 30.01 0 rad 3.0
           set fl off
           set mech force 5e-4
           solve
        end_command
        command
           set fl on
           solve age 14400         ;4*60*60=14400;计算4小时渗流
        end_command
        command
            get_gp_maxydispyong
            print fish
        end_command
        if dxmax<+1.6 then
             exit
        endif
        save lumda.sav
     endloop
end

;设置显示
;create disp1
;plot set back white
;plot set rot 20 0 30
;plot set mag 1
;plot add axes yellow
;plot bl gr
;plot add contour disp out on
;plot show

wdd


	https://bbs.yantuchina.com 访问内容超出本站范围，不能确定是否安全


	关闭选中1篇全选

帖子

[求助]跪求 [复制链接]