中国地质大学数值模拟上机报告讲解

flaCd数值模拟上机报告计算模型分别如图1、2、3所示，边坡倾角分别为30。、45。、60。，岩土体参数为:容重r=2500kg/m3,弹性模量E=1X108Pa,泊松比=0.3,抗拉强度叫=0.8><106Pa,内聚力C=4.2X104Pa,摩擦角(()=17°试用FLAC3D软件建立单位厚度的计算模型，并进行网格剖分，参数赋值，设定合理的边界条件，利用FLAC3D软件分别计算不同坡角情况下边坡的稳定性，并进行结果分析。附换算公式：1kN/m3=100kg/m3E110°c剪切弹性模量：G=——=1=0.385108(Pa)2(1口)2(10.3)体积弹性模量：K=一E一=—U0——=0.833108(Pa)3(1-2」)3(1-20.3)坡度为30的情况图1倾角为30°的边坡(单位：m)算例分析：命令流：new;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;建立网格模型genzonebrickp0000p110000p2020p30040size50110genzonebrickp040040p1100040p240240p374.64060p4100240&p574.64260p6100060p7100260size30110;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>；定义本构模型modelmohrpropdensity2500.0bulk8.3E7shear3.8E7coh42000.0tens0.8E6friction17;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;设置边界条件fixxyzrangez-0.10.1fixxrangex99.9100.1fixxrangex-0.10.1fixy;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;设置重力加速度setgravity=10.0;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;设定初始条件inisxx0.0syy0.0szz0.0sxy0.0sxz0.0szz0.0inixvel0.0yvel0.0zvel0.0inixdis0.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_loopendcalfos以NUtD3.叔step32345ModelPerspecltve22：障MMtriN圳252013CenterX;5.QQ0e和5V：I.IWfOOOZ：3.C00et{X)1Dist2.775fi+002Ratartion:X:0.000Y：0.000Z：(WTOM蜂Ang.:22.500Magfac-Q.oaDe^OOOItascaCDnsuttngGroup.Inc.Minn就网同MNUSA图1网格剖分图CenterRatation:X;5.QQ0eCenterRatation:X;5.QQ0e和5X:0.000V：I.IWfOOOY：0,000l.m.oooe啾z：o.woDist2.775fi+M2M眇1Ang.:2250a图2速度矢量图rL^C3r>3.uoStep羽※ModelPerspecl^eNw252口心Magfac=0QQOeiOOOVelocityMaximum=3J63e-005Unes^teItascaCDnsuttngGroup.Inc.Minn就网同MNUSACenterRatation:X;5.000CenterRatation:X;5.000斜OSX:0.000Y：1-ITOOefOOOY：0,000z.m.twe必1z：0.(x»Dist2.775a+M2hlsgjAng.:22.500Step32345ModelPerspecltve22:19:59MonNw252013〃月3.叔ContourofVelocityMag.Magfac-QBOQOe*OQOI00000^000132.5000e-0062.500G(h006h5,0000^0065.口OOOfl^m7.5DOD»<K)6圈7.5DOQHKJ6to1.0000&^0«-1QQO(*O05S12500HM5l.25OOMlO5to15000^0051.5DOa&^OSto1750060511.75DDfr<KHtDZDOODMOfi二ZOOOOM蜘口125O*K»12.25DQ^OOStO2.5000^005-25QOONW5S之750M■媾-2.750Qfr^06to3.0DOO&^O5■3.0000frOBb3,1fi27fr00fiInter同2.5e-006ItascaConsuttngGroupJnc.Minneapolis,MNUSA图3速度等值线图Step32345ModelPwspKtive2217:31MmNw252013Rotation:X:D.M0Y0.000Rotation:X:D.M0Y0.000Z:0.000Mag.:1Ang:22.5MContourMDi酒师如tMag.H^scaC^utti^Grcwdm，

USACenterK:S.OOCe+DOlY:1.000&fD002:lDD0e*001CM2775Mo02Magfac=D.OO&e+CO0IO.OXJOe+DOOto2.500De-D012.5(XXte-OO1to5.0000e^0015.0000^001ta7.5(X)(te-GO17.KXXMI01to1.0000什DODII.OOOOfl+DODto1.250te40M-125GG叶口DO忸1.500D€4G00-1.50COe+ODDtot.750te*GO0-1.7500a+DD0toZOOOte4QOO■2.00003+000to2.2500^40002.25C0«+DD0to25001珀号GOO2.5(X)0s+DD0to27500e-000“2.7500s+DDDt3S.OOOte4QOaI3.000Cte+DD0tD3.2阳/000Intertfal=2.5a-001FLACSD3.00St带能"5Mode?Perapet5ive22:13:15MonNov252G13C^larX5000e+00lY:1QOOe+DOOZ:3,000e4-traiDist:2.7756t002Rotation:X：0,0®Y：O.OOGZ:Q.DDDMag.:1Ang:22.500C^larX5000e+00lY:1QOOe+DOOZ:3,000e4-traiDist:2.7756t002Rotation:X：0,0®Y：O.OOGZ:Q.DDDMag.:1Ang:22.500ContourofShearStrainIncteuiemWaghc=0,008*000GradientCdlCulatiDnI4.970(te-QO4ta2,口003022.StX10fi-002to50000Z®馅7.5000+002格1.aOOOe-COIto1.25008-001to1.5tW®1tp5.0000^0027击⑩*猊1.0000e-00tf.25QOfrOOft.50COe-OOt_1.75We-Wlto2.0000e^0f-2QQ00e-W1to2.25006-001-2,2500^001to2.5CKKte-OD1|2.50008-001to2.5而7201lnter;ai=；.5e-0CC图图5剪应变增量云图ItascaCvnsuMhgGjuup,Inc.Minneapolis,mnUSACwierX:5.Q00fi+DD1Yzi.OOOe+OOOCwierX:5.Q00fi+DD1Yzi.OOOe+OOO2:lOMefQOlDisr：W.¥5p+0G2Surface£•,枷=a.iMOei-oOQCoDtourofSXXI.'aghc=O.OOOe+OOOGradientGalculdtionFLAC3n3.00St用期"EMwf^PerspKtive22:51iBTueNov舶2013Rotation;X:D.CXX)Y；O.QM2:G.DDOMag.:1Ang.22.5004日比兄他tn40000e十口的-9.DDD0e+0MtD-S.OOOOe+DOS■9.000064005to^,OOOOs+M5-7OOOOe-KXJSto-600006+005-S.3DDte40D5tD-S.D000fl+DD5-5300064005tc^000Oe+005国DDDDe+OT5to^.D000a+DD5-3.00W6*映匕小叩雎+唯-2.000064005to-1(HJOOe+005O.CKCOe+OOGO.OOQOe+OOOto3E蝌心1曲门球=t.da+005itasc自ConsultingGroup.Inc.Minne叩口限MNUSA图6最小主应力云图CenterRdt加口n:CenterRdt加口n:X;5.00te+l301X:0.000丫:i.crooei-ooQ¥：0.000z：也愉Dist2.775a+M2M呼1Ang.:22.500Step32U&ModelPerspeclive22：51.&5lb?Nov2e2013SurfeceMagf^c-OOQOe<UOOContourofSZZMagfaciO.OOOetOOOGradientC4cuWion4449施梃tD^1.400Gfi+C06d.3<XKte4(Xl€tD-1.2000e4OT6-1.1000e4006CD^1,0000e+(X16-9a»oe-KJO5co-aooooe+cos-7.0000e40Cl51D-6.0000e+C05-5.(X»0a-b005to-I.OCOOs-bOOS21.3003+005■OQOOe+OO5mDODOO^OOOOODOOe+DOOto5,5615卓.003interval-l.0e*005ItascaCdhs眦吗Group,Inc.Minn髓网同MNUSA图7最大主应力云图安全系数：Endingat:15:24:16TueNou262013Fish>printfsFs=1.453515625白白Ge+4白白因此，最终计算边坡稳定性系数为1.453516坡度为45°的情况404061004-4图2倾角为45°的边坡（单位：m）算例分析：命令流：new;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;建立网格模型genzonebrickp0000pl10000p2020p30040size50110genzonebrickp040040p1100040p240240p360060p4100240&p560260p6100060p7100260size30110;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;定义本构模型modelmohrpropdensity2500.0bulk8.3E7shear3.8E7coh42000.0tens0.8E6friction17;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;设置边界条件fixxyzrangez-0.10.1fixxrangex99.9100.1fixxrangex-0.10.1fixy;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;设置重力加速度setgravity=10.0;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;设定初始条件inisxx0.0syy0.0szz0.0sxy0.0sxz0.0szz0.0inixvel0.0yvel0.0zvel0.0inixdis0.0ydis0.0zdis0.0;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;安全系数求解defcalfosminf=0.1

maxf=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_loopendcalfosSisp17M5F.ladelP*rs|>Ktiv&22^5:47NmiIVw252013Cwiter:X:fi.OMe+DDlY：ljMte+0MZ：WDBtMlDist:2.775#+0G2帛利■=D.OOOe+DDDit^caGOfUuttin^GroupInc.MNW4pDli5,MNUSA图1网格剖分图图2图2速度矢量图FLACSD3.00Ste^179^5ModtaPerapeclive22:^1:10MoflNov252Q13CertarRotati口n;TOC\o"1-5"\h\zX.5000e+00lX:0.000Y:1QOOe+DOOY：O.OGOZ:3.00084-001Z:O.DDDDist:2.7758t002Mag.:1Ang:22T500Surfece**=O.OOOei-OOOMagt«-D.oooei-QooVelocityMaximum-1.901e-GG4LinestyIe•tsscaConsuftngGroup.Inc.Minneapolis,MNUSArLA<JSI>3.00St#179犯Mod«Pers0ec5jve22:4*05M画NSW252013RMon:X:0.000RMon:X:0.000Y：0.000Z:O.DOOMag.:fAng:22500X：5000e+00iY:1QOIte+DOQZ:3.0006+001Dist:277sB+002ContourofVelocityMag.UagfK=O.OOOefOOOIQ.QQQOe+OOQ102,0000e-0052.0M(te-CO5to4.0000+0054.00004005to6.0000^00560000fr<X>5tog.OOWZOJ1.0000e^04t.2000^0041.4G0Ge-O041.E00MD41.8000204',9012^X)4S0000e-005toL0000e-QD4to1.2OOQs-OO4t0I.JOOOfi-OMto1GOOOeXmto1aoooe-co4to1.0000e^04t.2000^0041.4G0Ge-O041.E00MD41.8000204',9012^X)4InfterTiM2.0^005图4位移等值线图日图里玄图5剪应变增FLAC3I>3.00St«p17945MwielPerspective23:01:07TueNov2£2013同燃机ion;同燃机ion;X:0.000Y：0.000Z：0.000Mag.:tAng..22.500Magrac=O.OCOe4oMContourofSXXMagfac-O.OWe+OOOGradientGalcii^ion3.0®施种05to“g.OOOOe+OM-9.000064005to-B.ODOOe+DOS百OOO加的力比-7.000Da+005-700GOe+005to-6.0000e+005■€.000064005la^,ODDOfi+OQS-5(XX»e4M5(o4,OOOOe+t)a5400OTe*加5m'S.ODODfl+QOS-3.000064(X)5tD-：.wO0e+00fi-2.00OO6+M5to-I.0000e+O05*1,00008+005toG.OOOOfr+OOOO.OOOOe+OO01&4.537264004Intep/al=I.陪।■哂ItascaConsultingGncup.Ifk.MlnneapKts.MNUSA图6最小主应力云图FLAC3n3.00StBp17945MMWPwrep粗:如日23由忡TueMov舶2013CerterRotation;X:5.D00fl+D01X:D.CXX)Y:1OOOe+DOOY;O.KO2:M*+QQ1Z:0.D00Dot;2.7758+002Mag.:1Ang22.500Mag阻-D.lMCei-OOaContourofSXX护agfe:=O,OOOe+OOOGradientGdlculationI-9.优)9兆他t口400006+。的-9.OD00e+005tD-a.0000e+D05舟.CiOOOfi伸to^.OOOOs+M5-7.00006+005to600006+005--fi.300084005ts-5.D000A+DD5二-5,000034005to^4.00006+005"^.3D!)De4O05to-3.00008+DD5二十0000渊阳to^0000e+005二^.OOOOe+OOStoH0000e+QQ5-toO.OCCOe+OOG■O.OOOOe+OOOlo睢地04lntsr;a(=1.k+D05CoutourofSZZMagfit=£i.QOOe+QOOitascsConsultingGroup.Inc.Uinneapol^MNUSA安全系数:Endingat:15:31:24TueNou262013Fish>printfsFs=1.108984375000e*O0O因此，最终边坡的稳定性系数为1.10894三、坡度为60的情况4048.456100

-图3倾角为60°的边坡（单位：m）算例分析：命令流：new;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;建立网格模型genzonebrickp0000pl10000p2020p30040size50110genzonebrickp040040p1100040p240240p351.55060p4100240&p551.55260p6100060p7100260size30110;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;定义本构模型modelmohrpropdensity2500.0bulk8.3E7shear3.8E7coh42000.0tens0.8E6friction17;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;设置边界条件fixxyzrangez-0.10.1fixxrangex99.9100.1fixxrangex-0.10.1fixy;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;设置重力加速度setgravity=10.0;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;设定初始条件inisxx0.0syy0.0szz0.0sxy0.0sxz0.0szz0.0inixvel0.0yvel0.0zvel0.0inixdis0.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_loopendcalfosStep314^6Model23:00:17WotNov252013RstatiDH：X:0.000Y:D.OOTZ:0.000Mag.:1RstatiDH：X:0.000Y:D.OOTZ:0.000Mag.:1Ang:22.500X:fi.OO&SH-DOlV:I.OMe+DDO2:10Dtef001Dist:2.775^4002SurfaceMaqfac=D.OOOe+COOtta^a口3忡忡gGrwp.Inc.Minneapolis.MNUSA

FL^C3£>3,00Step3i40tModelPerspecirve23：«42MonNw262Q13CenterRotamnCenterRotamnX;5.00De^001X:0.000¥：I.OOTefQOOY：0,000Z：3.000et00lz：&.<x»Dist2.775fi+002M呼Ang.:22.500VelocityMaxkrun=i.306e-oo3皿勒拒ItascaConsutti^Group.rnc.MinneapoJis,MNUSA图2速度矢量图ELACSD3.00Ste^3l40eMS白户的pedive23:Q1:36MonNQV252013CertarX.5000e+00lY1QOQe+OOOZ:3.00084-001Dtst:2.7758+002Rotstion:X：0.000Y：W0Z:Q.DDDMag.:1Ang:22.500CoutourofVelocityMag.t£agtac-O.OOOei-OOO_Q.OOOOe+OOO101D000e-(XMf.tKKKte-QMtD2.0000e-004工DOOOflgta5.0000^0043.0000^)04to4,000X4X)44.0000e-004to&00He@45.aOOOfl-CQ4ta6.0000^0046.00008-004to7.aC(Kte-004T.tMMMte-tXMto6.00Mfl-004SW00e-0O4to9,00We349000*004曲t.QOOOe-OOJI.KXJOfrCOJtat.1(KK}fl-<M3I.IOOOa^OOStai.2MKte-OO31.2000^001tot.3000^00313000e-CC3tot.3060eO«inter闻"tae-C04itS5：sConsultingSroup.Inc.Minneapolis,MNUSAFLAC3L>3.00Center：Rotetton:X:5.0Dte+ai1Center：Rotetton:X:5.0Dte+ai1X:o.owYJ.OOOe-KWOY：0.W0Z：3.0DOe+QOlZ:口.MODist2.775&f0C2Mag.:iAng.:22.500ContourofDisplactmeutMag.blagfac=O.OCOe4OMID.ODDOfi+DOOto5.0000640005.00006+000ID1.0000640011.0DODfl+D01tnF.5000843D11.5000e+00l102.00X)640012.ODD0S+D01to2.5000e40012.5000e+001肉3.(KWe4O0i-3.00QQe+W1103.500064001~3.5DDDfi+tM1to4.CKXX)b-fOD1_J4,0000e+00llo4.5000e*00l~4.50D0fi+Daitn5.0000e+DD1_5.00006+001105.500064001~5.5DD0e+Daito6,Q000e40D1E.0000e+001to6.500064001-C,5D00e+Q0lto7.000064001■7.O0D0a+001io7.44520+301interval-5峥000ItsscaConsultingGroq),Inc.MlnneapOte.MNUSA图4位移等值线图FLAC3D3.00St砰31406Mod9lPerspective23:01:17MonNw252013ConferRation:X:5.0006+011X:0.000Y：t«0e-KK)0Y：MWZ：3.DD0ti+0QlZ：0.000Dist2.775frtW2Mag.:tAng.:22.500ContourofSbearStrainIncrementMagFac=。除谢GradientCHcg口n■-2.0G9Ge-002toO.OQOOeHHM15.000DS-001tDI.DDOOs+OOOI.SOOOe+OOOlo2.000064000"2.50Dte+00Qto3,Q000e400D3.5D006+000忸4.0QOOei<VO4.50006+000IC5.0000640005.5DDtte+OOOto6.00DDe40006.5000e+0001o7.000064000|F.ODOOfi+DOato7.1136B4OD0Interval-5,0e-30iJaItascaConsuttingGroq)iInc.MinneafWfes.MMUSAFLACSD3.00曲学3N时ModtfPerapec5ive23:38:32TueNov262D13RDteion:X：0.0MRDteion:X：0.0MY：0.000Z:Q.DDDMag.:1Ang:22.500CedarX.5.000e+00l七1QOOe+DOOZ:3.00084001Dat:2773*002Surface»/as做=0,008*000ContourofSXXMagtsc-D.oooei-QooGradientCalculdHonI-9.163DS4O05to-9.D000e+D05-9.0000e+005to-3D000e+D£l5-9000064005to-70000e+005-7.300064005to-600006+005-6.3DD0B4OD5tD-5.D000s+DD5■5.000064(X)5tD-4.0000fi+0054MoM心g4晒)06+。的-3000064005to-200006+0053000064005b-1.0000e+005■1.D33!)B4(M5tDO.CCOO^MCO.OOMJfi+DDDto5.3B52a-F(»4Inlen1^-1.0+005ItascaConsyttingGroup.Inc,Minneapolis.MNUSA图6最小主应力云图FLAC3D3.00St.31406Mod?lPerspective23:3S:07TueNov262013RoHton:x：o.ooo忙0.QWRoHton:x：o.ooo忙0.QWZ:0,000Mag.;1Ang.:22.500Y;i.M0e-K>00Z:IDQQe+QOlDist2.775e+QO2Surfaceblagfac=0.(M0e+0MContourofSXXMagfac-O.OOOe+OOOGradientCalcMun目.15W&t85此省如00e+005-3.00006+005to-B.ODDOe+DOS^.OOOOa+005W-7.0DODfi+005-7,Q(XMe+O05E0^,O000e+005■6.000064005ta^,ODODfl+OQS-5000064005to4.W00e+005-4.000064005m^,ODODfl+DOS-3.000064(X)5(O-2.QQO0e+00fi-2.0C(X)640[i5to-1.0000e+Q05■1.0000b*005mG.OOMeH-KO0.0000e+000l&5.静睡+004IntB^al=I.G&4-GQ5ContourofSZZMacrae=o.oOOe+OOO'tascaConsultingGroup.Ifk.WlnneafKts.MNUSA安全系数：Endingat:15:36:29Tuehtou262013Fish>printfsfS=8-95703125O090e-001因此，最终边坡的稳定性系数为0.8957分析：根据三种角度的边坡的位移矢量图，当坡角为30。时，边坡的最大位移发生在边坡后壁某一深度处，这主要是重力梯度的作用结果，在坡脚处的位移很小，可见边坡的角度很小时，边坡基本不会发生沿坡面倾向方向的运动；由位移等值线图可见，边坡的最大位移都发生在坡脚处，且坡角越大的边坡，坡脚处的最终位移越大。通过边坡剪切应变增量和速度矢量图，可以明显看到除30°的边坡外，45°和60°的边坡都存在贯通的塑性区域，即潜在滑动面，速度矢量图有利的佐证了质疑判断，因滑动面外侧区域各网格点的速度明显大于其他区域，说明这一区域出现明显的滑动。由三种角度边坡的安全系数可见：坡脚越大，边坡稳定性越差，当坡脚达到某一角度时，边坡的稳定性将小于1,即可能发生破坏。四、坡度为60°的边坡开挖情况100算例分析：命令流：new;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;建立网格模型genzonebrickp0000pl10000p2020p30040size50110genzonebrickp040040p1100040p240240p350050&p4100240p550250p6100050p7100250size30110genzonebrickp053050p1100050p253250p363060&p4100250p563260p6100060p7100260size15110genzonebrickp045.77050p153050p245.77250p351.55060&p453250p551.55260p663060p763260size15110groupexc1genzonebrickp040040p150050p240240p345.77050&p450250p545.77250&size30110groupexc2groupsection2rangey02groupexc2groupsection1rangey02groupexc1attachface;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;定义本构模型modelasprodensity2500bulk0.83e8she0.38e8;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;设置边界条件fixxyzrangez-.1.1fixxrangex99.9100.1fixxrangex-0.10.1fixy>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;设置重力加速度setgravity00-10.0;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;设定初始条件inixdisp0ydisp0zdisp0inixvel0yvel0zvel0;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>;开挖modmohrprodensity2500bulk0.83e8she0.38e8fric17coh4.2e4ten0.8e6dila20modnullrangegroupsection1modnullrangegroupsection2;安全系数求解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=fsendjfendjoopendcalfos

FLAC3I>3.00Step29634ModelPerspective14：33;20TueN5vZ62013Center.Rotation:X;5OMtetWiX:0.000Y：1.DOOMOOD¥:0.0001.3.0006+00110.000DiSt：2.7756+002Mag/1月叫二Z15ODSurfaceM竭fat=O.OOOh+DQO■JiflxaConsultingGroup.Inc.MinwaptjIis.WNUSA图1网格剖分图FLAC3I>3.00Step296S4ModelPeripactiveii：55：22TueNcv£62013CenterX;5.0005*001Y:1.DOO^*ODDZ3.(X»e+001DiSt：2.775e+CO2Rotation：X:0.000¥:0.000Z:CenterX;5.0005*001Y:1.DOO^*ODDZ3.(X»e+001DiSt：2.775e+CO2CodtourofVelocityHga.M蝴at-O.QOOb+DQOIDOOODe+OOOto5.0000e-0065.000De-005to1.GOODe-OOi1.QOOOb-004to1.5ODDfi-O04_J1$000+004io2.0000e-0042.000DB-G04lo2.5O0Da-OD4'2.SOOOe^)4to3,d0&0e-004二30Koe■CWB3.5M0e-M4150006-004lo4.tJOMfl-OD44.0000e-004to4.5000e-0044.5COOe-OQ4to5.QQWe-W45.DCODB-GO4to5.50DDfi-004I-!5.5000+004to5.8535e-O04interval=5.IWKKSurfaceM/原-0.00064oooVelocityMafflimums5.355fi-OO4KascaConsultingGroupInc.MhmapolSsMNUSA图2速度矢量图FLAC3I>3.00Step296S4ModelPerspscti-ze他沛配MeNov汨加13Center:X;Center:X;50005tW1Y:1.Da)t*0DD1.3.0006+001DiSt：2775frte02Rotation：X:0.000¥:0.00010,000Mag；iArg.:①500ContourofVcIoclI)'Mag.M竭fat=O.OOte+DOOIO.OOOOe+OOOto$.0000+0065.(K)ODe-O05to1.000*541.0000&004toI.SODDe-0041.5000fr<X)4102.00006-004_2.000Dfi-004to2.5O0DH-OD4-2,5000e-004to3.00(幅004_3(HXX)e-CQ4to3.5000e-W4-3.5QODB-QO4to4.0000+0044,00016-004to4.50006-0044SQOO^CWto5.00306-004-S.DOODb-004to5.50000-004I-!5£OOOXX)4to5.85J5S-O04interval=5,O-0MItascaConsultingGroupInc.MhnMp0Ns,MNUSA图3速度等值线图FLAC3I>3.00Step29634ModelPerspective14：53;52TueWov262013Center.Rotation:X;5OOOetWi40.000r1.D«te*OOD¥:0.0001.3.000e+00lZ0,000DiSt：2.775e+0Q2M的.1Arg.:22.5ODContourofDisplacimcntMag.M蝴O.OOOs+DQO■O.OCOOe+OOOt