中国地质大学数值模拟上机报告讲解Word文档格式.docx
- 文档编号:8341445
- 上传时间:2023-05-11
- 格式:DOCX
- 页数:29
- 大小:1.32MB
中国地质大学数值模拟上机报告讲解Word文档格式.docx
《中国地质大学数值模拟上机报告讲解Word文档格式.docx》由会员分享,可在线阅读,更多相关《中国地质大学数值模拟上机报告讲解Word文档格式.docx(29页珍藏版)》请在冰点文库上搜索。
;
设置边界条件
fixxyzrangez-0.10.1
fixxrangex99.9100.1
fixxrangex-0.10.1
fixy
>
设置重力加速度
setgravity=10.0
设定初始条件
inisxx0.0syy0.0szz0.0sxy0.0sxz0.0szz0.0
inixvel0.0yvel0.0zvel0.0
inixdis0.0ydis0.0zdis0.0
安全系数求解defcalfosminf=0.1maxf=2.2loopwhilemaxf-minf>
0.01fs=(maxf+minf)/2.0refric=atan(0.30/fs)*180/3.14recoh=42000/fscommandprofricrefriccohrecohsetmechratio1e-5solvestep5000printfsend_commandaa=mech_ratioifaa<
1e-5thenminf=fselsemaxf=fsend_ifend_loop
end
calfos
rL^C3I>
3.00
Center:
Ratation:
X;
S.OOOe^OOl
X:
0.000
V:
I.IWfOOO
Y:
0,000
Z:
1咖旳01
z:
伽
Dist2.775a+M2
Usgj1
Ang.:
22.&
M
Slep32345ModelPerspective
22.l9.D4MonNpv253013
Migf^c=QQQQe*QOQ
图1网格剖分图
ItascaConsuttr^Group.inc.Minneap^is.MNUSA
l.啲1
o.wo
Dist2.775fi+M2
22.500
图2速度矢量图
rL^C3£
33
Sep32345ModelPerspective22:
i%40MonNw252W
Magfac-0QQOevOOO
Velocity
Maximum^3J63e-005
Unes^te
ItascaCDnsuttngGroup.Inc.NlriiwiHta.MNUSA
ContourofVelocityMag,
Magfac-QBOQOe*OQO
IO.OOOOfr+OOOto2.5000e-0062.500G<
h006h5.0000^0065.DOOH垢曲7.5D0D»
O06
K7.5DOQHM6to1.0000乐蹴
■何0毎005B12500^005[.25OOMO5曲15000^005I.SDOO&
^OStflUSDOeOS
1ITBOMKHtDNDOOdJfi[2.0000e-(KI«
h2250ZK
12.25DQ^OOStO2.5000^005
_2.500(X105102.7500^'
2.7&
Wte^05to3.DD0陆0(再
■a.OOOOfrOBb3.1627fr00fi
Inter间二2.5e-006
ItascaConsuttngGroupJnc.Minneapolis,MNUSA
图3速度等值线图
Step32345ModelPwspKtive2217:
31MmNw252013
Center
K:
S.OOCe+DOl
Y:
1.000&
+D00
2:
10DOe*001
Dist:
2.775e-0G2
Contourof圧朋賊ei玳对Mag.
Magfac=D.OO&
e+CO0
IO.OOCOe+DOOto2.500te-DD1
2.50Xte-001toS.OOOOe^OOl
S.OOQOfrOOlta7.5<
X)(te-GO17.5000*0®
la
II.OOXte+DODto
'
1.2500a+DDDto
I.KXJOa+DDDto
1.75OCte+DO0t3ZOOOte+GOO
=2.00003+000to2.2500^40002.25C(te+DD0to25000^0002.5(X)0s+DD0to27500^-000
"
2.7500s+DDDt33.000te4GOO
I3.0000e+DD0tD3.2449护000
IntBryal=2.5e-001
PLAC3D3,00
Step自加酣Mod^Perapeclive22:
13:
15MonNov252013
Ce^anRDtation:
X.5000e+00lX:
1QOOe+DOOY:
O.OGO
Z:
3.00084-001Z:
Q.DDD
2.7758*002Mag.:
1
Ang:
ContourofSliearStrainIntKdient
Uaghc=0.0如000
4.9700a-OG4ta2.500Cfl-OQ2
2.StX10fi-002ta5.0000^0025.0000^0219tSOOOeg7.5000fr002toi.OOOOea*1.a(KJ0e-C01to1.25CKte-(Mt1.25008-001tai.5000e-0011.5tW®
1tp
175006-C01to2.0000e^»
i
2QQOOe-WIto2.25Wfi-«
122500^001ta站OOO&
ODY
2.50008-001172.5767fr^0t
lntwvai=2.5&
^2
M畀pConsultingGnoup.IncMinneapolis,mnUSA
rLAesn3.oo
StBq32345M«
fe^Perepwtivfi22:
5l:
lSTueNw262013
CwierRotation:
5.D00fl+DD1X:
D.CXX)
1OOOe+OOOY;
O.WO
S.MOefQOl2:
O.DDO
Dot:
2.775BtOG2Mag.:
Ang22.500
Surface
姑妙站sD.OOOei-oOaCoDtcmrofSXX
l.'
agfK^O.OOte+OOO
GradtaniCalculation
-9.0O95e-K»
5tD-500006+005-9.DDDOe+OMtD-S.OOOOe+005■a.DODOB+ODSto^.OOOOfl+M5-7OOOOe-KXJSto-600006+005-S.3DDte40D5tD-S.D000fl+DD5-530006+005tc^000Oe+005J.DDDDe+OT5to^.D000a+DD5-3.{W0e4(Wt?
-2t)000fi+M5-2.000064005to-1OQOOe+M5■I.MDDfi-bOOStoO.CKCOe+OOGO.OOQOe+OOOto3BOGGs-f004
lnlervdi=t.da+005
tascBConsultingGmup.Inc.Hinneapohs.MNUSA
图5剪应变增量云图
图6最小主应力云图
JFLAGD3,00
Surfece
Magf^c-OOQOe<
UOO
ContourofSZZMagfac^O.OOOetOOOGradientC^cuWion
規449曲伽tD^1.400Gfi+M6d.3<
XKte4006tD-1.2000e4OT6-1.1000e4006to^1.0(»
0e+(X16-9ocooe-KJOSco-aooooe+cos-7.0000e40Cl5!
D-6.0000e+C05-5.<
X»
0a-b005to诃.0000§
伽-I.OCOOe-bOOSED-2.3000e4QQ5^tOQOOe^OOStoDODOO^OOOOODOOe+DOOto5,56156#003
interval-l.0e*005
ItascaConsuttngGroup,Inc.MlnnBEfHis.MNUSA
图7最大主应力云图
安全系数:
Endingat:
15:
24:
16TueNou262013
Fish>
printfs
Fs=1,45351562509Oe+090
因此,最终计算边坡稳定性系数为1.453516
坡度为45°
40
100
44
图2倾角为45°
m)算例分析:
命令流:
genzonebrickp0000pl10000p2020p30040size50110
genzonebrickp040040p1100040p240240p360060p4100240&
p560260p6100060p7100260size30110
设置重力加速度
设定初始条件
安全系数求解
defcalfos
minf=0.1
maxf=2.2
loopwhilemaxf-minf>
0.01fs=(maxf+minf)/2.0refric=atan(0.30/fs)*180/3.14recoh=42000/fs
command
profricrefriccohrecohsetmechratio1e-5solvestep5000printfs
end_commandaa=mech_ratioifaa<
1e-5thenminf=fselsemaxf=fs
end」fend_loopendcalfos
Sisp17M5F.ladelP#rs|>
Ktiv&
22:
39:
471.1snNov252013
Cwiter:
fi.OMe+DDl¥
:
1jK»
€+OMZ:
WDBtMl□st:
2.775#+0G2
绛179^6ModmPerapeclive22:
^1:
10MoflNov252Q13
RotatiDH;
X-.5000e+00l
1QOOe+DOO
3.000e+001
Dis*:
27756*002
Mag.:
f
Sutfece
Vagtac=0.000e<
-000
M阿就-D.oooei-Qoo
Maximum=1.9Q1e-Q04
Linastyle
ttascaConsultingGroup,Inc.Uinneapol^MNUSA
FLA<
J3I>
Stei&
179*45M&
ltfPers(?
ec5jve22:
«
:
D5MonNov252Q13
Censn
Rot^iDn:
X:
5000e+00i
¥
1DOOe+DOO
□tst:
2.7758*002
T
CoutourofVelocityMag.
Magtoc=O.OOO^OW
1.0000e^04t.2000^0041.4G0Ge-O041.500MD41.8000e^04■.9Q12^XM
IO.OOOOe+OOOto2.OOOOe^Q52.00000-005to4.0000^0054.0000fi-C05to6.0000e-005&
00Kfeg讪&
000ZW5二SQOOOe-OOSto二I.OOOOe-QCMtO
_1.20003-004toI.JOOOfi-OMto
「1&
QQ0©
伽悟
□Ii.soooe-awto
lntarc^=2.0e^N6
图4位移等值线图
FLAC3H3.00
Step179^5Mod科PerapedMe
22:
46:
19MofiNov252013
£
阳常RDtatlon;
X.5000e+00l¥
、0.000
1QOOe+OOOY:
a.DDD
2.775BtQO2Msg.:
CoGtour&
fDi^laeementMat.
gfac-O.OOOei-lMO
O.OOOOe+DOO102.5000e-0ai2.5QOOS-OO1toS.oaooe-ooii.tXXMfrCOIta7.5000^0017SOOOe-COito1OOOOe+OOO1KJOOe+QQO(a1.25006+000(o
1.SOOOs+WOta
1.TSOOfi+DMlaSDOOOe+OOOU2.25DOe+OQOlo15M0a+DQ0(aitssaConsuMngGroup,Inc.Minneapolis,MNUSA
Step17945Mod4fPersp«
-ive
4147ManNov252013
CoutourofShearStrainIncreiueat
I4.63908-004to2.5aOfta-OQ22.5tX)0frCa2tD5.000052
50000fr<
X^toWOWeas7.5000fr<
x^toi.OOOOfrOOt
■1.aoooe-coita1.2SOOfrOOt
1.25008-001tot.5000e-00l
1.5000^001to1.7500^001
1.7500fr^1to2.00006^01_20000e-001to2.25004-001_2.2500^001to站000XMY|2.50008-001to2.5517s-00t
FLAC3I>
Stop17945MwielPerspective23:
01:
07TueNov2£
2013
Center'
5.0Dte+OHY;
i.OOOe-KWOZ:
3.DQ0e+Q0l
Cist2.775^X2
Ration;
0.OD0
t
Ang..22.500
Magrac=O.OCOe+oM
ContourofSXX
Magfac-O.OWe+OOO
GradientGalciiaSion
-9.00炯<
W5to收QWOe+005-9.000064005to-日.00006*005£
00008+005la-7.000Da+005JOOGOe+OOSto-B.OOOOe+OOi■€.000064005la^.ODDOfi+DaS-5(XX»
e4M5(o4,OMOe+OaS400OTe*005忖^.ODODfl+QOS-3.000064(X)5tD-:
.wO0e+00fi-2.00OO6+M5to-I.0000e+O05*1.00008+005toG.OOOOfr+OOOO.OOOOe+OO01&
4.537264004
lntep/al=I.G^i-QGS
ItascaConsultingGncup.Ifk.MlnneapKts.MNUSA
Step17945M«
f^PerepK*ivfi23:
01.31TueNev262013
Cwfler
Rotation:
X=5.00(te+DD1
x:
d.ooo
1OOOe+DOO
Y;
O.OM
M.畑+QQ1
2:
Q.DOO
Dot;
Mag.^1
Ang:
\!
白gfit=D.QMei-OOO
Graci'
BntG^lculatiDn
HagfK二O.OOte+OOC
I-9.W98e-K»
5tD^OOOOe+OOS-g.oDooe+oos恰-a.ooooe+wjs■a.DOOOfi+ODSto^.OOOOs+MS_-7.00006+005to-6.0000S+005--S.3D0te4O05ts-5.D000fl+DD52-5.000084005to^.OOOOe+OOS-^.DDDDe+OOSto-3.0000a+DD5_-3.MWe4W5t?
-20000e+M52龙DOOOe*(X)5to”1Q000e+QQ5_-1.DKH)fi-f<
M5tDO.OCCOe+OOG
■O.OOOOe+OOOlo3抽404
1血曲=VOe+005
CoutourofSZZ
Magfit=£
i.QOOe+QOO
31:
24TueNou262013
Fs-1・1O8984375000e*60O
因此,最终边坡的稳定性系数为1.10894
三、坡度为60°
48.45
图3倾角为60°
genzonebrickp040040p1100040p240240p351.55060p4100240&
p551.55260p6100060p7100260size30110
inisxx
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 中国地质大学 数值 模拟 上机 报告 讲解
![提示](https://static.bingdoc.com/images/bang_tan.gif)