flac3d60應(yīng)用案例說明書使用

簡單的坡度在Hoek-Brown材 ExampleApplications在Cam-Clay在Cam-ClayFLAC3D解決的一系列示例應(yīng)用程序問題的文檔。這些示例應(yīng)用程序演示了可以應(yīng)用FLAC3D的各種問題。Index:Examples.指數(shù):FLAC3D潛在應(yīng)用程序的一個(gè)示例。向用戶提供的示例將定期更新。我10°陡比傳統(tǒng)的二1252條件在x=0平面上占優(yōu)勢（見圖3）面下半部分滲流。在穩(wěn)態(tài)條件下，該水位將導(dǎo)致圖5所示的孔隙壓力分布。選取該模型的強(qiáng)度參數(shù)，將FLAC3D結(jié)果與HoekandBray（1981）的圓形破壞圖圖6顯示了使用哪個(gè)圖表作為水流動(dòng)條件的函數(shù)。在我們的例子中，使用的圖表=cγHF0.06。具體重量，γ，25000Nm3和高度，H2537.5kPa的DXF文件。在這種情況下，DXFfish語言“生成水f3dat”，但它也可以很容易地來自任何來源。模型被指定為一個(gè)-庫侖材料模型和幾個(gè)性質(zhì)2002500kg/m32600kg/m3z方向作用于10m/sec2處。幾何計(jì)數(shù)范圍元素用于將飽和密度分配給水位以下的區(qū)域。滾子邊界沿模型兩側(cè)放置，模型底部固定安全系數(shù)采用強(qiáng)度折減法，采用模型安全因子8所示的位移等高線圖來描述由此產(chǎn)生的破壞面;該方案是在修復(fù)了微小的”斜這個(gè)問題也可以在平面應(yīng)變模式（fle”SL-PS.DAT”）和軸對(duì)稱模式（fle”slaxis.dat”）FLACFLAC3D模型匹配垂直部分的幾何（參見圖9。平面應(yīng)變模型的安全系數(shù)計(jì)算與圓形破壞圖的安全系數(shù)計(jì)算吻合，F(xiàn)1.61。破壞時(shí)的型繪制了y=-30垂直面上的位移等值線和矢量圖。DataFile數(shù)據(jù)文件的2.PillarLoadsatIntersectingTunnelsσxx=25MPa（東西）σyy=30MPa（北-南）σzz=17MPa（垂直）1。1:巖體材料性能摩擦角35°（峰值）了使用內(nèi)置的構(gòu)建塊塊創(chuàng)建工具FLAC3D創(chuàng)建的fnal網(wǎng)格。這個(gè)問題的數(shù)據(jù)是”支柱”。本節(jié)末尾列出了f3dat”。Results7顯示了地基平均豎向應(yīng)力（即，垂直的和反應(yīng)的基礎(chǔ)模型的面積除以:13米×7.5米）。注意，平均垂直29mpa89顯PillarStrengthk為強(qiáng)度參數(shù)，V為柱體積，R為柱寬高比，a、b為經(jīng)驗(yàn)參數(shù)。這些表達(dá)式是基其中A是作業(yè)區(qū)C件的1/8方形柱。圖12顯示了該模型底部豎向應(yīng)力的演化。柱內(nèi)的峰值應(yīng)力為域會(huì)導(dǎo)致的軟化，因?yàn)槲覀冊诮o定的軟化曲線的應(yīng)變軸上移動(dòng)得更快。為了修正這種下系數(shù)進(jìn)行縮放。在模擬高質(zhì)量（90）巖體時(shí)，由于使用完全（或接近完全）脆性可能1620之間，因此，假設(shè)=18，軟化速率將縮放0.9倍。關(guān)于軟化模型對(duì)廣義情況的網(wǎng)格依賴性的更詳細(xì)討論見<linktolocationintheUser’sGuidesection>.格林沃爾德，hP?；羧A斯（H.C.Howarth）和哈特曼（I.Hartmann）。礦務(wù)局，技術(shù)文件第605號(hào)（1939）。吉尼亞煤礦學(xué)報(bào)》，第43-66頁（1957年）。Munro55-67（1967）.Data3、ExcavationinaSaturatedSoil以開挖中心點(diǎn)為原點(diǎn)，z軸向下，形成一個(gè)坐標(biāo)系統(tǒng)。在這個(gè)例子中，邊界在|x|=12m和|y|12m1提供了問題條件的示意圖。FLAC3D12米×12米×12米;12×12×121m×1注意，在這個(gè)問題上，ρd+nρw土壤的飽和密度等于墻上密度。因此，初始狀態(tài)的特點(diǎn)是均勻的孔隙壓力分布p=ρwgzσzz（ρdnρw）廣州。水平應(yīng)力發(fā)σxxσyy0.4ρd（0.4（n1）1）ρw]廣州。這些應(yīng)力應(yīng)用于開挖墻體，模擬初流動(dòng)條件。不努力表示固結(jié)效應(yīng)的真實(shí)時(shí)間尺度;我們只對(duì)fnal穩(wěn)態(tài)感。挖掘3米深的切口;Modelingm;隧道底部位于39.5米。以隧道底板原點(diǎn)為原點(diǎn)的坐標(biāo)系統(tǒng);z軸向上，y軸沿隧道軸向上。FLAC3D2所示。該模型約有15000個(gè)區(qū);沿隧道軸線有50個(gè)區(qū)域。質(zhì)量密度為2200kg/m3。==代表預(yù)支護(hù)混凝土和噴射混凝土（30m后）的殼體結(jié)構(gòu)單元?；炷烈r砌如圖4所示為了安裝具有s元素的連續(xù)內(nèi)襯，將為新s元素提供與上一步中現(xiàn)有s元素相同的標(biāo)識(shí)號(hào)。然后，新s將使用連接到現(xiàn)有s的現(xiàn)有節(jié)點(diǎn)。新殼層的初始應(yīng)力為Results隧道開挖30m后，隧道周邊塑性區(qū)范圍如圖6所示。不同的顏分剪切破壞和拉伸正如前面的例子所演示的，F(xiàn)LAC3D非常適合于模擬前進(jìn)的隧道開挖和支護(hù)的行為。這個(gè)例子演示了如何使用GUIFLAC3D構(gòu)建塊窗格構(gòu)建交叉隧道。創(chuàng)建了兩個(gè)示例45°的十字路口的兩個(gè)同等大小的隧道。第二個(gè)例子十字路口的隧道在45°fll，對(duì)隧道和巖體進(jìn)行分組，并沿邊緣劃分 15及材料2中應(yīng)力狀態(tài)沿鹽界面的變化。σxx=-56MPa（east-west）σyy=-28MPa（north-south）σzz=-35MPa（vertical）2質(zhì)見表1。zFLAC3D模型。由于幾何問題網(wǎng)格。由于對(duì)稱的條件，只模擬了1/8的。fgure中的三個(gè)或四個(gè)網(wǎng)格以交互方式顯x0y0z0平面為對(duì)稱平面，具有滾動(dòng)邊界條件。x250m處的遠(yuǎn)場邊界在界（z=250。該模型受到的初始原位應(yīng)力狀態(tài)σxx=-56MPa，σyy=-28MPa，σzz這個(gè)洞將在0.1年內(nèi)開挖。模擬開挖過程中，0.1年后，將沿巖洞施加的正壓力由地應(yīng)力值降至7mpa。然后將壓力保持在7mpa不變，總時(shí)間為1年。地應(yīng)力和鹽的蠕變行為相一致。第二階段采用表格荷載法進(jìn)行洞室開挖，0.1年內(nèi)7mpa。最后，在第三階段，對(duì)模型的蠕變響應(yīng)進(jìn)行為期1年的。即使沒有開挖，鹽在施加的各向異性應(yīng)力場作用下也會(huì)發(fā)生蠕變。這是因?yàn)閮缏傻妮喞?，有一個(gè)應(yīng)力不連續(xù)在鹽和材料2接口。圖5為巖洞3處1年后的位移歷史。圖6和圖7分別為巖洞中心垂直平面和水平x-和y-）方向位移最大值約為0.41m，垂直（-z-）方向位移最大值約為0.05m。7、rtnofBolerenaatmn在600米深的鹽層中開采了兩個(gè)房間和一個(gè)交叉流。柱子上鉆了一個(gè)水平鉆孔，WIPP-reference7.7年，考慮了零時(shí)刻的瞬時(shí)采礦和鉆井。這個(gè)示例問題是Sandia國家進(jìn)行的中等規(guī)模鉆孔測試的簡化模型Ageo，191個(gè)水平面位于|z|53.36m，三個(gè)垂直線位于|y|44.51mx39.62米。房間、平巷、鉆孔的布置如圖1所示，模型前半部分對(duì)稱平面被移除以作說明。試驗(yàn)室寬5.49米，高13.57MPa的壓力，以表示來自覆蓋層的負(fù)載，以及兩側(cè)室、橫移、鉆孔同時(shí)開挖。在整個(gè)7.7年的模擬過程中，溫度保持在300k不變。bulkmodulus，K1.65GPashearmodulus，G1.00GPaWIPP-modelexponent，n4.90WIPP-modelconstant，D5.79×10-36Pa-4.9s-1activationenergy，Q12，000cal/molWIPP-modelconstant，A4.56WIPP-modelconstant，B127.0secondarycreeprate，5.39×10-8s-15），建立系統(tǒng)（在蠕變過程時(shí)間尺度上）的初始”瞬時(shí)”彈性響應(yīng)。上。在分析結(jié)果之前，對(duì)模型進(jìn)行了7.7年的循環(huán)。261MB3.4GHzIni7計(jì)算機(jī)上運(yùn)行完整的模擬大約需要1分鐘。6.6cm）位于柱中心（D點(diǎn)（a點(diǎn)）9cm，鉆孔端（a點(diǎn)）10.5cm沿著漂流點(diǎn)（C點(diǎn)）的中等長度的柱子。567.7年后鉆孔壁的垂x0、x3.5m、x6.1mx9.14m。從柱子中心到房897.7年。雖然由于沿井眼軸線使用的地層數(shù)量相對(duì)較少，因此在分析本例中的實(shí)際數(shù)字時(shí)應(yīng)十分謹(jǐn)慎，但它們似乎反映了在本例開始時(shí)的中等規(guī)模鉆孔試驗(yàn)中觀察到的總體趨勢。但是，應(yīng)該強(qiáng)調(diào)，”三維有限元分析方法之應(yīng)用”，國立大學(xué)土木工程8、Axialand lLoadingofaConcretePile混凝土的軸向和橫向荷1z0處，是一個(gè)自由曲面。模型z=-8mzfxed，模型兩側(cè)|x|=8m處為滾子邊界，y=8m處為滾子邊界。面令中，區(qū)域面是分開的，這樣兩個(gè)通用的網(wǎng)格點(diǎn)也會(huì)分開。一個(gè)水平的水位，粘土的濕密度被分配到水位以下的區(qū)域。在下一步的分析中，0線性增加到期望值。這對(duì)這個(gè)問題特別有幫助，因?yàn)榛炷梁驼惩林g的效應(yīng)將被忽略，使系統(tǒng)穩(wěn)態(tài)響應(yīng)的識(shí)別變得更加。并將其作為歷史記錄。為了方便起見，帽面上的網(wǎng)格點(diǎn)作為映射在符號(hào)帽中。樁頂軸向應(yīng)力與軸向位移的關(guān)系圖如圖3所示，為速度為0.0~5x10-8爬坡，爬坡距離為30000步。注意，在這個(gè)響應(yīng)中最初只觀察到少量的振蕩。100kn時(shí)的響應(yīng)。這種荷載是通過樁頂100kn4cm。這是通過在樁頂施加水帽水平速度為1x10-7。線如圖6所示9、UndrainedCylindricalExpansioninaCam-ClayMedium不排水圓柱的應(yīng)力和孔壓變化。沒有考慮測量裝置fnite長度的影響。aa0Cam-clay材料的性質(zhì)如下（Carteretal.1979）:undrainedcohesion1shearmodulus74×soilconstantslopeofnormalconsolidationlineslopeofelasticswellinglinereferencepressurespecifcvolumeatreference（vλ）σ或=σ'θ=-1.65立方，σ'z=3Cu和初始超孔隙壓力問題=十分之一和高度h（參見圖1中的前視圖。FLAC3Da0/a0=1L/a0=200h/a0=131個(gè)等p0=2.1Cuq0=1.35Cu，預(yù)固結(jié)壓力Camclay屈服函數(shù)計(jì)算（見本構(gòu)模型><section1），（1）2.70Cu。對(duì)于這個(gè)問題，將過固結(jié)比R的值設(shè)為Rpc0/p01.29specifcv0K0默認(rèn)情況下，Biot1。Biot模量（水體積模量除以孔隙率）10010-5100，000步，使空腔半徑在壓力試驗(yàn)結(jié)束時(shí)3q/Cuq/Cu1.777處存點(diǎn)后的屈服路徑為（見Cam-Clay試樣的排水和不排水三軸壓縮變形問題5a=2a0時(shí)有效應(yīng)力和孔隙壓力的徑向分布。在土體處于臨界狀態(tài)的環(huán)形土體方，應(yīng)力和孔隙壓力向它們的原位值演化。這些結(jié)果與Carteretal.（1979）等人結(jié)Wroth19710、SimulationofPull-TestsforFullyBondedRockRein ment全粘結(jié)巖石在現(xiàn)場對(duì)注漿加固的小段進(jìn)行拉拔試驗(yàn)。通常，50厘米或更長的管段被灌入鉆孔。這些米段拉拔試驗(yàn)結(jié)果如圖1所示。這些圖用噸/米表示，變形用毫米表示。方法1:用索單元建模巖石加固的變量所需令在cable命令中進(jìn)行了描述。位長度與變形的關(guān)系表示時(shí)，這很容易做到，如圖1所示。假設(shè)錨索不屈服，則灌漿剪峰值拉拔強(qiáng)度值。這一假設(shè)是合理的，因?yàn)殇摰那?yīng)力通常高于與拉拔力相關(guān)的corr，115mm0.5m1200MPa（以下各節(jié)FLAC3D模型中假設(shè)）的典型電纜，最大屈服力為212kN50%大于抗拔=15.2mm鋼絲的灌漿內(nèi)聚力值為噸/米的峰值抗剪承載力。在這種情況下，grout-數(shù)據(jù)略”01。f3datFLAC3D進(jìn)行拉拔測試的簡單表示。電纜終端節(jié)點(diǎn)以315.2mm電纜的拉拔力與電纜位移（分別為歷史力和disp）4610mm、16.5mm、17.5mm時(shí)電纜軸向力分布情況。注意，在開始，并迅速沿整個(gè)電纜長度（圖。在那個(gè)階段，力電纜一端只是之和3%左右，電纜應(yīng)該斷裂。在拉伸斷裂拉拔試驗(yàn)中給出了元件拉伸斷裂的實(shí)例。一般來說，隨著作用在電纜上的有效壓力（p'）的增加，電纜的剪切粘結(jié)強(qiáng)度也會(huì)增7830°confning壓力，即p'=24MN/m2fgure3進(jìn)行比較。這些fgures表明，隨著初始澆注壓力的增加，破壞程度增加，說明了索-巖界面的摩擦特性。采用數(shù)據(jù)膝關(guān)節(jié)”Pull02”15.2mm斷線電纜11、Approach2:ModelingRockRein mentUsingModifedPileElements方法2:利用修正樁單元對(duì)巖石加固進(jìn)行建模在接下來的章節(jié)中，使用FLAC3D中修改后的樁單元邏輯對(duì)巖石加固進(jìn)行建在Pile命令中描述了創(chuàng)建和查看與modifed元素相關(guān)的變量所需令。注意，修改后的樁單元是通過在命令下給出結(jié)構(gòu)樁的性能rockbolt-flag來激活的。數(shù)據(jù)表”Pulltest04”。f3dat”創(chuàng)建一個(gè)沒有的拉測試。該問題相當(dāng)于在不承認(rèn)的拉flePulltest01中的模型進(jìn)行對(duì)比。f3dat”，目前的模型考慮最大拉拔荷載與接近1時(shí)的無壓拔強(qiáng)度相當(dāng)（這一微小差異可能是由于模型中包含的螺數(shù)據(jù)稍為”Pulltest05”。f3dat”defnes的拉拔試驗(yàn)與前一節(jié)描述的拉拔試驗(yàn)類似，拉拔強(qiáng)度無粘結(jié)，但峰值后剪切粘結(jié)強(qiáng)度減弱。用coupl-on-table對(duì)灌漿的粘結(jié)強(qiáng)度軟化進(jìn)行了控制，如錨桿行為中所述。通過表cct給出了剪切位移與內(nèi)聚力弱化之間的關(guān)系該模型相當(dāng)于糖果的拉拔強(qiáng)度討論，這一次的糖果應(yīng)力，4MN/m2。根據(jù)剪力耦合彈簧行。在8數(shù)據(jù)略”Pulltest07。f3dat”創(chuàng)建了一個(gè)模型，該模型演示了couplinglimitingtable選項(xiàng)151010n/m/m，粘結(jié)強(qiáng)度Hyetta·JWFBawden和R.DReichert。”巖體配制對(duì)全灌漿錨索螺栓粘結(jié)強(qiáng)度的影響”，巖體配制。分鐘?？茖W(xué)。&Geomech。Abstr29（5、503-524（1992）12、WheelLoadoveraBuriedPipey01.25m處施加車輪載荷。車輪載荷作用于表面的四個(gè)網(wǎng)格FLAC3D模型假設(shè)通過管道中心的垂直平面和輪載荷之間的垂直平面是對(duì)稱的。模型1所示。x0平面對(duì)應(yīng)于通過管道中心的垂直面，y0平面對(duì)應(yīng)于車輪載荷之間的平面。該模型包含1344個(gè)區(qū)域，其中fner區(qū)域位于應(yīng)用的車輪載荷附近。mz方向由結(jié)構(gòu)（殼）2s448個(gè)結(jié)構(gòu)單元和255個(gè)結(jié)構(gòu)節(jié)點(diǎn)組成。鋼的楊氏模量是227GPa泊松比是0.25。先在重力載荷作用下達(dá)到平衡狀態(tài);1.11mm0.92mm。將位移重置為零，然后應(yīng)用z速度來模擬車輪負(fù)載。力持續(xù)增加:6000步后壓力達(dá)到750kPa。地表沉降由位移等值線和矢量圖如圖4所示。這顯示了直接在車輪荷載下的最大管內(nèi)表面應(yīng)力分布用最小（即圖9中的主應(yīng)力。FISH函數(shù)mises計(jì)算管道中von13、EmbankmentLoadingonaCam-Clay直有效應(yīng)力之比為6/13。路基的重量是模擬的應(yīng)用附加費(fèi)和排水發(fā)生在土壤表面。k10-12（m/s）/（Pa/m。土模量是平均有效壓力和土樣體積的函Pa的順序，或兩個(gè)數(shù)量級(jí)低于水的體積彈性模量（2×108Pa千瓦。因此，在本例中，擴(kuò)散系數(shù)c由土壤材料控制。其大小可由c=k（k+4/3G）估計(jì)，約為10-6m2/s。tc=L2/cL為模型高度。用L10mtc是圖1所示的模型利用了半對(duì)稱性。它的尺寸是20米寬，10米深。請注意，模型的寬于分析平面兩側(cè)的滾子邊界（y方向、沿對(duì)稱線的滾子邊界和模型遠(yuǎn)邊界（x方向以及模型基底處x、y、z方向的fxed位移。粘土的最大體積彈性模量（bulk- 置為5×106Pa，大約是兩倍的值初始值的實(shí)際體積彈性模量（體積）的底部粘土層。450kPa的壓力，模擬路堤的可能需要調(diào)整這些參數(shù)。有關(guān)這些的討論，請參閱耦合流和力學(xué)計(jì)算。在計(jì)算過程中監(jiān)測應(yīng)力、孔隙壓力和垂直位移。本頁末列出了這個(gè)問題的數(shù)據(jù)。不透水混凝土沉箱墻支撐。[*]開挖長約84米，寬36.5米，開挖深度23.5米。墻體厚和鐵釘支撐。在這個(gè)例子中，坐標(biāo)軸系統(tǒng)被破壞，使得原點(diǎn)在開挖的中心以下，z軸向上。開挖深度z=26.5m，地表z=50m。26.5m的深度，建立初始平衡。初始狀態(tài)包括高度z=39.5m的水位和不透水墻的孔隙壓力分布。這個(gè)問題的K0值（靜土壓力系數(shù)）0.6。2FLAC3D中的流體流動(dòng)計(jì)算，得到與脫水階段應(yīng)用的邊界條件密切對(duì)應(yīng)的分布。這初始水位位于地表以下10.5米（z=39.5米3?；炷翂w、樁、后幫的性能分別見表3、表4、表5。士兵樁沿墻間距為2.25張安裝后在加載表5中列出2中的土壤性質(zhì)作為”干”性質(zhì)輸化，壁面周圍的孔隙壓力保持水壓不變。模型初始分布如圖3所示。士兵樁和電纜樁的位置。這個(gè)描述可以像CAD數(shù)據(jù)一樣容易地導(dǎo)入。樁位置如圖4所示。接?；亟訋缀误w如圖4所示。意，最靠近士兵樁的后接段未粘結(jié)（結(jié)構(gòu)電纜性能grout-00。3、4、50，在開挖深度以下1m。模型對(duì)脫水的響應(yīng)分兩步計(jì)算。首先，僅對(duì)流體流動(dòng)進(jìn)行計(jì)算，使模型在變化外產(chǎn)生孔隙壓力。圖6、7、8為切割3、4、5段脫水后孔隙壓力分布。58.1mm17所示。這大致與經(jīng)歷最高拉伸載荷的第三級(jí)回縮階段表現(xiàn)為水平線段。開挖遵循降水階段，墻體向開挖方向移動(dòng)，直至電纜預(yù)發(fā) Janes合作編寫的演示模型。15、DewateredConstructionofaBracedExcavation析了三個(gè)施工階段:1）20m深度;2）2米深;3）安裝水平支撐，開挖至10米深度。根據(jù)需要，可以在FLAC3D分析中加入額外的挖掘階段。層。初始水位位于地表。支撐開挖的施工墻附近的地表沉降。采用兩種不同的材料模型，-庫侖模型和CYSoil模型，說明了CYSoil模型的制定與PH模型有共同的組成部分。如表5所示，硬化CYSoil特性與PH模型）因此，模型響應(yīng)不應(yīng)該是相同的。5給出了硬化后的CYSoil6所示。FISH函數(shù)設(shè)置從硬化土的性質(zhì)推導(dǎo)出CYSoil的性質(zhì)。注意，cap-yield表面參數(shù)，被設(shè)置為1。[3]則=FISHini_cy中分配帽和初始?jí)毫Γ琍H值的性質(zhì)，和應(yīng)力有關(guān)。這些有效的主應(yīng)力在FISHini_ph中初始化，10倍。通過這樣做，界面的變形性將對(duì)整個(gè)模型的遵從性和計(jì)算速度產(chǎn)生最小的影響。于正截面和剪力筋，選用550MPa/m的代表性值。確定了水的性質(zhì)、孔隙度和滲透率。請注意，F(xiàn)LAC3D所要求的”滲透性”實(shí)際上10.0m/21000kgm3。滾子邊界位于網(wǎng)格的左fxed。在保存初始狀態(tài)之前，給出了地下水位初始預(yù)開挖應(yīng)力狀態(tài)的靜力平衡。在初始狀態(tài)下，孔隙壓力和總（有效）應(yīng)隙壓力分布后，fx沿頂部邊界和側(cè)邊邊界的孔隙壓力，以及沿頂部邊界的飽和度，滿足流動(dòng)條件。初始分布如圖3所示。上述步驟對(duì)于-庫侖、CYSoil和PH示例來說是常見的。CYSoil和PH模型需要額外的設(shè)置。在CYSoilPH模型中，如果需要，有必要指定應(yīng)力相關(guān)的特性和硬化表函數(shù)。如”BracedExcavation-cy.f3dat”中所述，這些函數(shù)用作FISH函數(shù)。首先執(zhí)行FISHcysoil的相關(guān)性質(zhì)。CYSoilmdoel的應(yīng)力相關(guān)特FISHini_cy中。類似地，在FISHini_ph中輸入PH模型的應(yīng)力相關(guān)圖4、圖5、圖6為現(xiàn)階段-庫侖模型、CYSoil模型和PH模型的總垂直應(yīng)力分布。與庫侖材料相比，CYSoil材料在墻體周圍的應(yīng)力分布略有不同。這可以歸因于CYSoil材料的應(yīng)力依賴剛度特性和與墻體相鄰的界面常數(shù)一致的剛度特性之間的差異。fxed孔隙壓力條件釋放，使孔隙壓力在降水過程中pore-pressure0命令完成的。當(dāng)此命令用于修改現(xiàn)有的孔隙壓力時(shí)，更改將被并應(yīng)用區(qū)域gridpointinitializepore-pressure命令，總應(yīng)力將保持不變，但有效應(yīng)力將進(jìn)行調(diào)整15厘米。模型狀態(tài)保存為摩爾庫侖材料的”mc-part-4”。對(duì)CYSoil材料重復(fù)這個(gè)脫水過程。fnal分布如圖9所示，與摩爾庫侖材料的孔對(duì)于PH材料，同樣重復(fù)這個(gè)脫水過程。fnal11所示，脫水引起的位移如圖12所示。含PH材料的模型在此階段保存為”PH-part-4”。2米現(xiàn)在可以開始了。從”mc-part-4”、”cy-part-4”和”ph-part-4”開始，設(shè)置流域。圖13或圖14顯示了挖掘2m時(shí)網(wǎng)格部分的近景。型和PH模型。在-庫侖模型中，開挖底部最大隆起約為3.4cm。在CYSoil模型中，開挖底部的最大隆起約為1.8cm。在PH模型中，開挖底部最大隆起約為2.7cm。墻的響應(yīng)也可以計(jì)算出來。例如，第一次開挖后墻體彎矩分布如圖16所示。請注意，可以使用FLAC3D中的線性繪圖項(xiàng)繪制墻體響應(yīng)的各種結(jié)果（例如，墻移、剪對(duì)于fnal開挖步驟，在墻體頂部安裝水平支柱，然后開挖至10米深。在這個(gè)練習(xí)（）=（10，1，40這就限制了光束只在垂直方向上運(yùn)動(dòng);旋轉(zhuǎn)xy位移受到抑制。此時(shí)，從表面（z40m）到粘土-1層底部（z10m）的所有區(qū)域都被刪除，我們準(zhǔn)備，fle”cypart-，fle”PH-料。此外，在CYSoil和PH材料中挖掘引起的位移程度比摩爾庫侖材料更容易被理CYSoil材料在墻后的表面沉降與PH21所示的地中，基坑底部發(fā)生的隆起比CYSoilPH材料中的隆起要大得多。這主要?dú)w因于CYSoil18CYSoil19PH材料位移等值線圖的對(duì)比中可以明顯看出。[1]改編自xis手冊（2002）的表4.1 xisBV。 （2002xisBV。 XIS版本8，手冊。:《中國社會(huì)科學(xué)》（2002。16、InstallationofaTriple-AnchoredExcavationWall 進(jìn)行了比較。實(shí)例驗(yàn)證了CYSoil模型和塑性硬化（PH）模型在深基坑問題中的適用性。20m，底層（2）100m6016.8橫膈膜墻采用區(qū)域建模（楊氏模量30GPa，泊松比0.15，厚度0.8第2層為0.38。第1步-初始化應(yīng)力狀態(tài)，包括水位，3米。第2步-激活防滲墻和低水位在坑?17.90米。第5步——開挖步驟2（?9.30米。第7步-開挖步驟3（?14.35米。第9步-開挖步驟4（?16.80米。100，z∈（-和y∈（0，2.7，其中x=0表示對(duì)稱平面的開挖。x0處，滿足對(duì)稱條件;fxedboundaryz100和x150處表示遠(yuǎn)端條件。在y軸上不允許變形來模擬平面應(yīng)變條件提是重力、水密度和表信息由zonewater命令設(shè)置。水位以上的土壤為干密度，地始化-應(yīng)力進(jìn)行初始化，摩擦力的計(jì)算K0是內(nèi)聯(lián)進(jìn)行的。庫侖模型賦于土體，僅使用求解彈性令模型，以避免可能出現(xiàn)的不均勻沉降。這個(gè)階段的數(shù)據(jù)是”AnchoredExcavation-1”。f3dat”和結(jié)果在fle”part-1.f3sav”中。特性如表2所示。界面元素被分配在墻體的側(cè)面和底面，以解釋墻體和土壤之間的滑移和分離。請注意，zoneattachdelete命令用于釋放早先分離的墻體和土壤之間的邊界，并zone接口節(jié)點(diǎn)initialize-stress初始化界面應(yīng)力，以加快收斂速度。界面摩擦角是假周圍的土壤中進(jìn)行調(diào)整。模型狀態(tài)在fle”part-2.f3sav”中。3.f3dat的FISH函數(shù)ini_p來完成這個(gè)FISH函數(shù)必須在分配了 sticHardening模型后的第一次循環(huán)之前執(zhí)行。注意，（-4.f3dat”設(shè)置為零。開挖下方左側(cè)邊界和開挖右側(cè)模型頂部的fxed孔隙壓力條件釋放，使孔隙壓力在降水過程中發(fā)生變化。模型中的位移初始化為零，這樣我們就可以監(jiān)water命令中用于修改孔隙壓力的effective關(guān)鍵字完成的。脫水誘導(dǎo)位移如圖4所示。這表明脫水引起的沉降量約為0.01m。模型狀態(tài)保存刪除區(qū)域，模擬開挖第1步（4.8m深度）開挖過程，并將模型循環(huán)至力學(xué)平衡5.f3dat”f3dat”100mpa901.6kN，已從原來的調(diào)整預(yù)應(yīng)力錨力（768kN）2.7/2.3倍，specifcation7.f3dat”8.f3dat9.f3dat”10.f3dat”-1.f3dat”錨固力如圖10所示，開挖至14.35m深度引起的位移如圖11所示。以向下為主，與g值一致。會(huì)議，第1卷，第305-314頁。e /LCPC，巴黎（2002） dTunnel 加強(qiáng)了隧道開能/效果仍然是一項(xiàng)的任務(wù)。工作面本身作為空腔的臨時(shí)支護(hù)，配筋單元提供的配筋Janinetal.（2012，2015）使用fnite單元平臺(tái)的 。Cheng和Lucarelli（2016）再次使用FLAC3D[*]中開發(fā)的塑性硬化模型進(jìn)行了這個(gè)例子。鉆孔表明，基巖具有較強(qiáng)的蝕變作用，具有較好的水平地層學(xué)特征。圖1顯示了是break默認(rèn)值。器和收斂目標(biāo)。這條隧道是用所謂的”ADECO-RS”方法挖掘的（Lunardi2008。該方法的噴射混凝土進(jìn)行了數(shù)值模擬。隧道壁面、工作面和仰拱特性見表2。（半寬度）×70（深度100FLAC3D內(nèi)置擠壓工在橫向上的擴(kuò)展為100m。隧道頂部的深度約為25米。網(wǎng)格已經(jīng)在感的區(qū)域周圍進(jìn)將土體暫設(shè)為庫侖模型得到原位豎向應(yīng)力，命令模型僅求解彈性問題。然后使用FISH函數(shù)和系數(shù)重置原位水平應(yīng)力。在實(shí)際地應(yīng)力的作用下，將土體模型轉(zhuǎn)化為塑性增也是空的。采用預(yù)埋樁單元模擬面板加固（fberglass單元）18m。每隔9米翻新一次屋面加固和屋面前撐。在幾何邏輯中，可以使用CAD程序創(chuàng)建一次面強(qiáng)化和前拱的模式，然后在隧道開挖時(shí)創(chuàng)建一個(gè)偏移量。然而，在本例中，我們選擇使用FISH函數(shù)GenFrontBolt和30m處開始沉降明顯，然后加0.20~0.25m5為垂程，Z。，A.Lucarelli?！彼苄杂不Ｐ虸I:標(biāo)定與驗(yàn)證”，《地質(zhì)力學(xué)應(yīng)用數(shù)值模擬-2016（第四屆Itasca應(yīng)用數(shù)值模擬研討會(huì) Gomez，C.Detournay，R.Hart，M.Nelson，ed。明尼阿波利斯:Itasca咨詢亞寧J.P.DDiasFEmeriaultR.KastnerHLeBissonnaisand195期，42-53。彈簧（2015。亞寧J.P.DDiasRKastnerFEmeriaultAGuilloux，和HLebissonnais?！蹦希?012（200818、SubsidenceaboveHorizontalCut這是一個(gè)平面應(yīng)變問題模型的尺寸（長度）×60100（深度。30米寬、1米高的水平0.75，初始化水平應(yīng)力。在切割的地板和屋頂之間設(shè)置一個(gè)接口，這樣切割可以在沒有屋頂和地板的情況下閉合。材料密度為1600kg/m3，假設(shè)重力為10m/s2。模型網(wǎng)格如圖1所示。Mohr-CoulombMohrT模型的拉伸破壞區(qū)域分別如圖23所示。Mohr-6）預(yù)測的沉降遠(yuǎn)大于本例中的Mohr-Coulomb模型（圖7）預(yù)測的沉降?？赡胬焖苄詰?yīng)變會(huì)裂紋的閉合。（見MohrT模型單層加卸載試驗(yàn)）的背壓（由壓應(yīng)力引起）將減少豎向變形和沉降。MohrT模型允許水平裂縫的完全閉合，并防止了MohrCoulomb模型的人為操作，該模型導(dǎo)致了對(duì)沉降的低估。除沉降規(guī)律19、PunchIndentationofaBondedMaterial種性質(zhì)的模型可以使用PFC來研究:1）PFC球組件，以確保材料在遠(yuǎn)離壓痕處的真實(shí)響應(yīng);2）PFC壁代表沖頭的主體。這種方法將是相當(dāng)昂貴的計(jì)算，因?yàn)樾枰蛞酝獾臎_壓壓痕區(qū)域。此外，沖頭應(yīng)力/變形在PFC中不能被評(píng)估，因?yàn)閴w只起到很小，而且球也很大，但是這個(gè)示例有效地演示了耦合PFCFLAC3D模型的機(jī)制。正如在PFC插件部分中所提到的，PFC墻壁是為表面區(qū)域所的。圖2顯示了所有PFC模型組件都存在的PFC墻壁和模型域。PFC球與壁面存在接觸;每個(gè)球面接觸處的接觸力耦合的PFC-FLAC3D在創(chuàng)建耦合的PFC-flac3d模型之前，必須加載PFC模塊。這些dll包含負(fù)責(zé)實(shí)例化所有與PFC相關(guān)的模型組件的代碼。此外，存在一組模塊將PFC墻壁與FLAC3D區(qū)域和sFLAC3D應(yīng)用程序打開時(shí)，這些模塊不會(huì)被加載。一個(gè)可以通過工具加載這些模塊?負(fù)載PFC菜單項(xiàng)或使用以下命令:載。請注意，在關(guān)閉FLAC3D應(yīng)用程序之前，其他模塊一直處于加載狀態(tài)。球云，使用模型隨機(jī)命令將隨機(jī)初始化為指定值。改變隨機(jī)將導(dǎo)致不同的球的實(shí)現(xiàn)，雖然統(tǒng)計(jì)上，實(shí)現(xiàn)將是相同的。與區(qū)域和結(jié)構(gòu)元素不同，所有PFC模型組件（即空間搜索。由于PFC總是在大應(yīng)變模式下運(yùn)行，因此必須顯式地為區(qū)域和結(jié)構(gòu)元素激活大應(yīng)變模式（modellargestrain命令。在默認(rèn)情況下，PFC使用球和團(tuán)塊的真實(shí)慣性特性來求解動(dòng)態(tài)運(yùn)動(dòng)方程。在這個(gè)模型中，我們希望在PFC有任意的，所以當(dāng)球重新排列時(shí)，所有的格點(diǎn)速度都是fxed。wall-zonecreate命令用于將墻壁facet包裝在surfacezone表面。如圖2所示，墻壁不會(huì)超出模型域。使用wallgenerate命令創(chuàng)建一個(gè)附加的墻來包圍空腔。PFC球和團(tuán)塊既可以在沒有初始的情況下創(chuàng)建，也可以具有任意初始。在這種情況下，balldistribution命令用于創(chuàng)建具有任意到平行六面體中的特定體積分?jǐn)?shù)的modelcmat命令指定接觸模型的行為。采取這些步驟使用contactmethod命令安裝并行鍵。此外，使用相同令設(shè)置并行鍵屬性。contactlin_mode=1意味著，對(duì)于這個(gè)BPM，所有的力量都是增量累積的;在此命令之前，由于球和球面接觸的絕對(duì)而形成的力量。這使得我們可0。因此，BPM處于平衡狀態(tài)，沒有內(nèi)應(yīng)力。關(guān)的區(qū)域網(wǎng)格點(diǎn)，從而允許BPM和continuum模型自然地交互。截面位移。請注意BPM和連續(xù)區(qū)域之間位移的連續(xù)性。此外，位移場是對(duì)稱的，因?yàn)?顯示了通過額外的壓痕破壞的鍵的位置。在這一點(diǎn)上，批量BPM由于局部的鍵斷裂而失敗。結(jié)示了體塊破壞前和體塊破壞后的測量應(yīng)力。BPM在這一點(diǎn)上能夠抵抗縮進(jìn)的斷點(diǎn)明顯更20、SleevedTriaxialTestofaBondedMaterial 數(shù)量的魚來模仿彈性膜的行為。在這里，從s 創(chuàng)建一個(gè)套，在套中合成BPM，如圖1所示。正如在PFC插件一節(jié)中提到的，PFC墻壁是基于外殼的結(jié)構(gòu)元素。PFC球與壁面的情況下創(chuàng)建的，并且平板以恒定的速度移動(dòng)。使用此示例作為基礎(chǔ)，可以輕松地確定BPM范圍。耦合的PFC-FLAC3D在創(chuàng)建耦合的PFC-flac3d模型之前，必須加載PFC模塊。這些dll包含負(fù)責(zé)實(shí)例化所有與PFC相關(guān)的模型組件的代碼。此外，存在一組模塊將PFC墻壁與FLAC3D區(qū)域和加載。一個(gè)可以通過工具加載這些Tools?LoadPFC菜單項(xiàng)或使用以下命令:載。請注意，在關(guān)閉FLAC3D應(yīng)用程序之前，其他模塊一直處于加載狀態(tài)。PFCFLAC3D外殼組成的圓柱形套筒內(nèi)創(chuàng)建的。為了在每次運(yùn)行數(shù)據(jù)觸發(fā)器時(shí)致不同的球的實(shí)現(xiàn)，雖然統(tǒng)計(jì)上，實(shí)現(xiàn)將是相同的。與區(qū)域和結(jié)構(gòu)元素不同，所有PFC化接觸檢測和空間搜索。由于PFC總是在大應(yīng)變模式下運(yùn)行，因此必須顯式地為區(qū)域和塊的真實(shí)慣性特性來求解動(dòng)態(tài)運(yùn)動(dòng)方程。在這個(gè)模型中，我們希望在PFC中使用時(shí)間步長縮放操作（modelmechanical命令FLAC3D行為的方法。垂直創(chuàng)建圓柱形套管。不使用”@”運(yùn)算符，可以使用括號(hào)符號(hào)（例如，'[]'）來防御FISH變量，并在命令中使用這些值。z方向直接指向圓筒的中心，yfxed，使球可以平衡內(nèi)套。wallresolutionwall分辨率模式設(shè)置為full。一般情況下，PFCnone，則會(huì)勾畫出面與球之間的比將壁面理想化為無小面時(shí)產(chǎn)生的接觸。整個(gè)方案被用來緩解這種情況，并且切割角被用來控制墻的”平滑度”。在pfc中的面墻部分詳細(xì)討論了墻的分辨率。PFC球和團(tuán)塊既可以在沒有初始的情況下創(chuàng)建，也可以具有任意初始。在這種情況下，balldistribution命令用于創(chuàng)建具有任意到平行六面體中的特定體積分?jǐn)?shù)的modelcmat命令指定接觸模型的行為。采取這些步驟提供一些摩擦，以消除額外的能量，從系統(tǒng)的球是循環(huán)到一個(gè)小的平均比例。contactcontact方法是一組操作，用于在聯(lián)系人上執(zhí)行任務(wù)，可能修改多個(gè)聯(lián)系人屬性。在鍵安fxing，幾周后所有接觸力均為零。因此，BPM處于平衡狀態(tài)，沒有內(nèi)應(yīng)力。最后，釋放套筒中部的殼體fgurestessBPM的初始聯(lián)系人。創(chuàng)建了無壓力BPM之后，需要進(jìn)行一些額外的工作來構(gòu)建執(zhí)行三軸測試的環(huán)境。首先是一個(gè)FISH函數(shù)，用于計(jì)算基于壓板的應(yīng)力和應(yīng)變。在這里，力在z方向上的壓板是用來計(jì)算平均應(yīng)力的壓板。應(yīng)變也根據(jù)試件高度計(jì)使用以下FISH函數(shù)遞增試樣的配制應(yīng)力:該功能的是逐步增加壓板和套筒上的澆注應(yīng)力。由于我們已經(jīng)為s指定了本地系統(tǒng)，因此使用結(jié)構(gòu)sapply命令來提供conning應(yīng)力非常簡單。同時(shí)，采用內(nèi)置壁面首先，通過拆卸套筒并對(duì)壓板的速度進(jìn)行固定，來進(jìn)行無側(cè)限抗壓強(qiáng)度（UCS）測figure所示。破裂的紐帶也顯示了BPM在失敗時(shí)的普遍損Pa675e4Pa1e5Pa的應(yīng)力應(yīng)變響應(yīng)和斷裂接觸鍵。下表顯示了BPM

?在Cam-Clay21、CylindricalHoleinanInfniteMohr-CoulombMaterial 該問題利用FLAC3D中的平面應(yīng)變條件對(duì)-庫侖塑性模型進(jìn)行了測試。）與圓柱體的長度相比，孔的半徑a較小，因此適用平面應(yīng)變條件。In式中，v為泊松比，為膨脹角G為剪切模量出于建模的目的，圖1所示的領(lǐng)域解決了這個(gè)問題，其中利用了四分之一對(duì)稱幾何的xz邊界位于距離孔軸線的fve孔直徑的距離。定義域的厚度選擇為孔徑的十分之被標(biāo)繪的半徑和歸一化的半徑，如圖4和圖5所示，歸一化位移分 表示 如圖6和圖7所示。注意，兩個(gè)流動(dòng)規(guī)則的數(shù)值應(yīng)力值不能在4.6%可以通過更適當(dāng)?shù)靥幚磉h(yuǎn)場條件來減少錯(cuò)誤，例如，使用Lame解決方案處理厚環(huán)。）通過”有關(guān)。用于生成關(guān)聯(lián)案例的f3dat”如下所示”是非關(guān)聯(lián)的。f3dat”只在分配的位移值上有所不同，可以在項(xiàng)目中找到。（30°在相關(guān)的情況下。通過”associated-nmd。f3dat”（如下所示）和”非關(guān)聯(lián)nmd”。f3dat”FLAC3D（”NMD”--12顯示了NMD48直接對(duì)應(yīng)（使用全十六進(jìn)制網(wǎng)格的結(jié)果?？傻木幠晔?，4，231-236（1969。22、CylindricalHoleinanInfniteHoek-Brown本文對(duì)英菲特Hoek-Brown介質(zhì)中的圓柱孔在地應(yīng)力場作用下，在內(nèi)部壓力作用下的應(yīng)230兆帕的infnite體中形成的。支撐由=5MPa的內(nèi)部壓力提供。 和按比例的內(nèi)部壓力由以下兩個(gè)方程確定FLAC3D模型為平面應(yīng)變模型，分析平面垂直于孔軸線。由于對(duì)稱20半徑處。本構(gòu)模型為Hoek-Brown-PAC在”check.f3dat”FISH函數(shù)中，對(duì)上述方程提供的徑向應(yīng)力、切向應(yīng)力和徑向位移的解析解進(jìn)行了編程。然后將分析結(jié)果和FLAC3D結(jié)果到表中進(jìn)行比較。表4顯示FISHFLAC3D結(jié)果進(jìn)行比較所需的輸入。將應(yīng)力和位移沿從孔中心測量的徑向路徑進(jìn)行比較。FISHradius_tunnel是隧道的半徑，rb_maxdefnes規(guī)范化路徑的長度并53.2m（1）的Fairhurst下巖體開挖的彈塑性響應(yīng)”36，777，809（1999。助理行政，106，1013-1035（1980。（1982（第15屆巖石力學(xué)研討會(huì) 集，1988年10月第31-38頁。多倫多:多倫多大學(xué)土木工程系（1988年。23、RoughStripFootingonaCohesiveFrictionlessMaterial粗糙的條形基礎(chǔ)上對(duì)于精確模擬的數(shù)值模型來說，預(yù)測穩(wěn)定塑性流動(dòng)條件下的倒塌荷載是的（Sloan材料（Tresca模型）的條形基礎(chǔ)承載力的確定。承載力的值是在地基下形成穩(wěn)定的塑性作為”Prandtl’swedge”問題解決方案的一部分而獲得的承載力由Terzaghi和m，遠(yuǎn)端z邊界位于基礎(chǔ)下方10m。域的厚度選擇為1m。y方向上，基礎(chǔ)最右格點(diǎn)在x方向自由。這個(gè)條件是合理的，因?yàn)樵谶吔缟系奈锢砑s束是模糊的，下一個(gè)網(wǎng)格點(diǎn)的速度為零，得到平均基腳壓力。半寬度a是原因。如果假設(shè)線性變化，那么對(duì)于這個(gè)問題，A0.5，a=3.5m。在結(jié)果和討論中討論數(shù)值模擬對(duì)應(yīng)的荷載-5所示，荷載為歸一化平均地基壓力，p/c，dispq，與514.2kPa的解析值相比，相對(duì)誤差為1.72%。制邊界段結(jié)束后半個(gè)區(qū)域?qū)挾忍帲?20.5）;注意，如果假設(shè)變異系數(shù)為0.63，F(xiàn)LAC3D的節(jié)點(diǎn)混合離散化（NMD）特性（文件”NMD.f3dat”）再次運(yùn)行相同的”nmd”f3grid。tet網(wǎng)格的網(wǎng)格點(diǎn)與十六進(jìn)制網(wǎng)格相同。在NMD下，承載力數(shù)值為524.1kPa1.9%514.2kPa。如前所述，承載力的誤差與基礎(chǔ)表觀寬度，sW。和先生?！庇糜邢拊椒ㄓ?jì)算倒塌荷載”，國立大學(xué)土木工程研究所。在Geomech方法。，6，47-76（1982。子（1967）24、SmoothCircularFootingonanAssociatedMohrCoulombMaterial光滑的失敗基礎(chǔ)摩擦20°角被發(fā)現(xiàn)2所示的四分之一圓柱。選取圓柱上底面的xyz軸為坐標(biāo)系。平板由半徑為a的圓盤段表示，域的半徑為15m，高度為10m。應(yīng)用于該領(lǐng)域的邊界條件如圖3所示。對(duì)稱邊界在x=0y=0處的位移分別限制在到網(wǎng)格點(diǎn)，代表在負(fù)z方向上的基腳的范圍。a3.165m。2.0級(jí)的速度×純米/步驟是應(yīng)用在基礎(chǔ)節(jié)點(diǎn)總共9600步伐。為歸一化承載力解析值，q/cdisp為歸一化豎向位移，u_z/aqkPa2010kPa0.4%展，第7卷，第295-340頁，紐約:愛思維爾科學(xué)（1975。Hopkins變形。反式。Soc。倫敦，A組，254（1036，1-45（1961）25、SmoothSquareFootingonaCohesiveFrictionlessMaterial 1975數(shù)值例子的基礎(chǔ)是廣場，由一個(gè)區(qū)域與半寬度（b。優(yōu)點(diǎn)是對(duì)稱的季度，和一個(gè)平的x軸和y軸，垂直方向上向上的-軸為坐標(biāo)系x、y、z邊界的位移在各個(gè)方向上受到限制，x0、y0平面對(duì)應(yīng)的對(duì)稱邊界的位移分別在x、y方向上受到限制。板光滑:xyz方向內(nèi)現(xiàn)3m×3m區(qū)域模擬加載的基礎(chǔ)。型中，a=3.5m，b=3.5m。（1數(shù)值模擬對(duì)應(yīng)的荷載-位移曲線如圖4所示，荷載為歸一化后的平均地基壓力，p/c;fle”createtteme.f3dat”tet67顯示NMD45直接對(duì)應(yīng)（allhex網(wǎng)格的結(jié)果?？梢钥闯?，結(jié)果是非常展，第7卷，第295-340頁，紐約:愛思維爾科學(xué)（1975。Drucker-26、UniaxialCompressiveStrengthofaJointedMaterialSample節(jié)理試樣1平面上，當(dāng)為k>0時(shí)，Cook（19791m4m的圓柱。選擇一種以位于柱體基部的x軸和z軸以及沿柱體軸的y48個(gè)徑向分布的區(qū)域（2所示y方向上施加均勻速度，誘導(dǎo)試樣壓縮。β090°。輸入文件使用一條魚函數(shù)（solveAll）βFISHsigmav計(jì)算的最終垂直應(yīng)力在每次運(yùn)行結(jié)束時(shí)添加到表中。這種方法允許整個(gè)參數(shù)分析保存在一個(gè)文件中。在區(qū)域面板應(yīng)用命令中使用的servo關(guān)鍵字控制應(yīng)用速度的大小。循環(huán)進(jìn)行，直到達(dá)到halt意，相對(duì)誤差小于2%的β值。Cook（197927、DrainedandUndrainedTriaxialCompressionTestonaCam-ClaySamplecam粘土試樣的排水和不排水三軸壓縮試驗(yàn)clay材料制成:初始試樣處于各向同性壓縮狀 Wood（1990）Cam-clay塑性理論的詳細(xì)討clay材料制成，常規(guī)三軸試驗(yàn)的平均壓力和偏應(yīng)力p可以表示為:式 為軸向應(yīng)力，為單元壓力 ， 隨著測試的進(jìn)行，路徑會(huì)向平面中臨界點(diǎn)的交點(diǎn)所表示的臨界狀體積應(yīng)變不再發(fā)生。因此，屈服面不會(huì)發(fā)生軟化或硬化;最終屈服面對(duì)應(yīng)于的值，對(duì)超孔隙壓力u由5×（）xy初始比規(guī)范體積在內(nèi)部計(jì)算（默認(rèn)值P’c0和平均壓力P0下膨脹100000.54×m/秒用于排水的例子。不排水試驗(yàn)，壓縮速0.5×106米/204020步。總共使’1p’1，q，vu23所示表示相對(duì)誤差小于2%。度過固結(jié)試樣孔隙壓力穩(wěn)定減小的對(duì)比（見圖15、圖16。 重度合并的情況。利用FISH函數(shù)對(duì)速度邊界條件進(jìn)行了應(yīng)用，并對(duì)仿真結(jié)果進(jìn)行了誤差:.:

（199030m600kpa300kpa125毫米厚噴射混凝土襯砌接觸應(yīng)力。這些計(jì)算值與和（1979）的解析解進(jìn)行了比較.應(yīng)力、剪切應(yīng)力和剪切應(yīng)力組成FLAC3D模型模擬了在平面應(yīng)變條件下存在各向異性雙軸應(yīng)力場的無限彈性地中圓形隧道的薄片。FLAC3D模型的幾何結(jié)構(gòu)如圖2所示。遠(yuǎn)場邊界放置在距離為隧道半徑的20倍的地方，以近似無限邊界。地應(yīng)力被安裝在所有區(qū)域，并作為作用于遠(yuǎn)場邊界的通過重新組合和融合適當(dāng)?shù)墓?jié)點(diǎn)本地系統(tǒng)，并指定適當(dāng)?shù)乃俣萬條件來實(shí)現(xiàn)的——為沿隧道邊界的24個(gè)區(qū)域。襯結(jié)構(gòu)元素與陰影網(wǎng)狀模式（見Liner-Rein 45tu7內(nèi)聚-剪切的性質(zhì)設(shè)置為零，襯墊與地面之間發(fā)生相對(duì)剪切運(yùn)動(dòng)（見圖。T和全滑移M912所示。當(dāng)襯墊x-軸沿著隧道軸（在全局y方向上z-軸垂直于殼體中部表面（指向內(nèi)。這對(duì)應(yīng)于z正極方向破壞的外光纖處產(chǎn)生正應(yīng)力。）與解析解相反。標(biāo)繪的剪應(yīng)力只描述大小;可以通過FISH函數(shù)打印或 解的定量比較。FISHFLAC3D模型中提取必要的值，并與解析值進(jìn)行比較。通過no_slipFISH10，可以使用無滑移或全滑移解析解。在節(jié)點(diǎn)處使用（2025中的異常值。分配給兩端兩個(gè)節(jié)點(diǎn)的面積是不相同的，因?yàn)樗?，H.H。，以及c.w.（c.w.Schwartz?！彼淼乐ёo(hù)的簡化分析”，J.29、Developmentof sticHingesinaStaticallyLoadedBeam靜載梁塑性定，這將產(chǎn)生一個(gè)運(yùn)動(dòng)學(xué)機(jī)制，從而導(dǎo)致系統(tǒng)的。問題如圖1所示。.MA=M1maxDMD=M2max對(duì)于給定的對(duì)稱幾何問題，在點(diǎn)C和點(diǎn)A之間還有另一個(gè)點(diǎn)D——圖中沒有。許點(diǎn)B上的滾子自由地向右移動(dòng)。作為FLAC3D模型的輸入。點(diǎn)A的彎矩是恒定的，等于Mpl。求出導(dǎo)致第二塑性鉸發(fā)生破壞時(shí)的XD最大彎矩坐標(biāo)。點(diǎn)處A反作用力由點(diǎn)B這個(gè)問題是通過使用兩種不同的方法來確定FLAC3D中的塑料鉸鏈來解決的（參見Properties單元（3所示21個(gè)節(jié)點(diǎn)，每個(gè)節(jié)點(diǎn)的極限塑性力矩為50mn-m。有19個(gè)連結(jié)（即，圖3中除節(jié)點(diǎn)外的所有節(jié)點(diǎn)均已;每個(gè)節(jié)點(diǎn)對(duì)解除一個(gè)鏈在兩種方法中應(yīng)用統(tǒng)一附加費(fèi)-統(tǒng)一附加費(fèi)的效果（2）由施加于節(jié)點(diǎn)的力表示。通過發(fā)布模型求解命令，找到了系統(tǒng)的平衡點(diǎn).繼續(xù)循環(huán)，直到局部不平衡力比低于態(tài)（wmax，現(xiàn)時(shí)附加費(fèi)水平是否達(dá)到均衡。這個(gè)孔是用來確定產(chǎn)生梁倒塌的載荷的。情況是5.80MN/m的負(fù)載。這些結(jié)果綁定了用解析解計(jì)算的值。flie”塑性矩580”可以得到坍塌前彎矩和剪力的分布。f3sav”—45。敗的位于距離6米的位置。這與預(yù)測值為5.86m的閉合形式解的解析解是合理一致Inc.（198930、SimplySupportedIsotropicRectangular teunderCombinedLa landDirectLoads簡支各向同性矩形板在橫向和直接聯(lián)合荷載作用下?lián)隙让娲_定為簡支矩形各向同性板的情況下，a和b，受均布橫向荷載P0和均力作用。問題條件如圖1所示。解析解由Ugural（1981pp.155-156）給出。該解決方案擴(kuò)展了板的小變形理論（即FLAC3D模型具有4×8的網(wǎng)格，包含128個(gè)s 元素，如圖2所示