巖土工程講義

巖土工程講義第1頁(yè)，共74頁(yè)，2023年，2月20日，星期四主要內(nèi)容軟件介紹動(dòng)力分析DynamicOption樁-土相互作用分析Interface隧道分析StructureElement流固耦合分析Fluid-MechanicalInteraction學(xué)習(xí)方法及經(jīng)驗(yàn)介紹第2頁(yè)，共74頁(yè)，2023年，2月20日，星期四主要內(nèi)容軟件介紹動(dòng)力分析DynamicOption樁-土相互作用分析Interface隧道分析StructureElement流固耦合分析Fluid-MechanicalInteraction學(xué)習(xí)方法及經(jīng)驗(yàn)介紹第3頁(yè)，共74頁(yè)，2023年，2月20日，星期四軟件介紹FastLagrangianAnalysisofContinua美國(guó)Itasca咨詢公司開發(fā)2D程序(1986)1990年代初引入中國(guó)有限差分法(FDM)DOS版→2.0→2.1→3.0→3.1Itasca其他軟件第4頁(yè)，共74頁(yè)，2023年，2月20日，星期四軟件介紹應(yīng)用巖土力學(xué)分析，例礦體滑坡、煤礦開采沉陷預(yù)測(cè)、水利樞紐巖體穩(wěn)定性分析、采礦巷道穩(wěn)定性研究等巖土工程、交通工程、采礦工程、水利工程、地質(zhì)工程特色大應(yīng)變模擬完全動(dòng)態(tài)運(yùn)動(dòng)方程使得FLAC3D在模擬物理上的不穩(wěn)定過程不存在數(shù)值上的障礙動(dòng)力分析功能——地震工程襯砌功能——地下工程可開發(fā)功能FISHVC++第5頁(yè)，共74頁(yè)，2023年，2月20日，星期四豐富的本構(gòu)模型第6頁(yè)，共74頁(yè)，2023年，2月20日，星期四軟件介紹有限差分法(FDM)古老的方法(上世紀(jì)40年代)用差分格式轉(zhuǎn)化控制方程中的微商格式流體力學(xué)；土工滲流問題；固結(jié)FDM&FEM的混合求解FDM的新進(jìn)展時(shí)間步與時(shí)間時(shí)間步的設(shè)置須滿足數(shù)值求解的穩(wěn)定SpecialOption中有特定的時(shí)間步要求動(dòng)力、滲流、流變、溫度中是真實(shí)時(shí)間第7頁(yè)，共74頁(yè)，2023年，2月20日，星期四Lagrangian網(wǎng)格源自流體力學(xué)中的拉格朗日法跟蹤流體質(zhì)點(diǎn)的運(yùn)動(dòng)狀態(tài)跟蹤固體力學(xué)中結(jié)點(diǎn)，按時(shí)步用Lagrangian法研究網(wǎng)格節(jié)點(diǎn)的運(yùn)動(dòng)跟蹤節(jié)點(diǎn)和單元隨材料移動(dòng)邊界和接觸面與單元的邊緣一致積分點(diǎn)隨材料移動(dòng)本構(gòu)方程在相同的材料點(diǎn)賦值固體力學(xué)大變形理論法國(guó)數(shù)學(xué)家、物理學(xué)家Lagrangian復(fù)雜的邊界條件精確描述材料的發(fā)展第8頁(yè)，共74頁(yè)，2023年，2月20日，星期四大應(yīng)變與小應(yīng)變SETsmall默認(rèn)為小應(yīng)變模式不更新節(jié)點(diǎn)坐標(biāo)不進(jìn)行應(yīng)力旋轉(zhuǎn)修正SETlarge大位移，大位移梯度，大轉(zhuǎn)角巖土工程中的大變形問題軟土的固結(jié)變形土體的開挖軟巖巷道地下硐室第9頁(yè)，共74頁(yè)，2023年，2月20日，星期四空間混合離散技術(shù)結(jié)構(gòu)域離散為可由四面體單元組合形成的五面體或六面體等單元以為基本單元(常應(yīng)力、常應(yīng)變)體應(yīng)變的計(jì)算：偏應(yīng)變的計(jì)算：+)/2=(第10頁(yè)，共74頁(yè)，2023年，2月20日，星期四FLAC3D的求解過程平衡方程(動(dòng)量方程)應(yīng)力—應(yīng)變關(guān)系(本構(gòu)模型)Gauss定律單元積分應(yīng)變率速度節(jié)點(diǎn)力新的應(yīng)力對(duì)所有的網(wǎng)格節(jié)點(diǎn)對(duì)所有單元第11頁(yè)，共74頁(yè)，2023年，2月20日，星期四Lagrangian格式動(dòng)量平衡方程F(t)m牛頓運(yùn)動(dòng)定律對(duì)于連續(xù)體在靜力平衡條件下，加速度項(xiàng)為0，方程變?yōu)槠胶夥匠痰?2頁(yè)，共74頁(yè)，2023年，2月20日，星期四例1：自由落體的模擬G=mgS=1/2gt2=20m命令流：configdyngenzonbrisize111inixmul0.1ym0.1zm0.1modelelaspropbulk3e8shear1e8inidens1000setgrav00-10solveage2第13頁(yè)，共74頁(yè)，2023年，2月20日，星期四例1：自由落體的模擬(movie)第14頁(yè)，共74頁(yè)，2023年，2月20日，星期四不平衡力平衡狀態(tài)F≈0最大不平衡力有所有單元確定SETmechforce???最大不平衡力比與內(nèi)力的比值SETmechratio???工程、計(jì)算中間階段 10-4論文、最終結(jié)果 10-5F≈0v≈0v≠0收斂、平衡不收斂、塑性流動(dòng)第15頁(yè)，共74頁(yè)，2023年，2月20日，星期四BLOCKstate*-n與*-p-n：此時(shí)刻進(jìn)入屈服狀態(tài)-p：曾經(jīng)進(jìn)入屈服狀態(tài)，現(xiàn)已退出程序預(yù)留了多個(gè)狀態(tài)空間供用戶定義液化的判斷溫度破壞的判斷PLOTblockstate*-n活動(dòng)塑性區(qū)貫通結(jié)合速度結(jié)果計(jì)算模型邊界選取的影響第16頁(yè)，共74頁(yè)，2023年，2月20日，星期四分析問題的過程建立網(wǎng)格初始條件邊界條件初始應(yīng)力平衡外荷載求解前處理后處理第17頁(yè)，共74頁(yè)，2023年，2月20日，星期四例2：一個(gè)最簡(jiǎn)單的例子genzonbrisize333;建立網(wǎng)格(前處理)modelelas;材料參數(shù)propbulk3e6shear1e6inidens2000;初始條件fixzranz-.1.1;邊界條件fixxranx-.1.1fixxranx2.93.1fixyrany-.1.1fixyrany2.93.1setgrav00-10solve;求解appnstr-10e4ranz3x12y12solveploconzd;后處理 切片功能RUNFLAC3D第18頁(yè)，共74頁(yè)，2023年，2月20日，星期四FLAC3D3.1的新特征64-bitVersion(64位計(jì)算)ParallelProcessingonMultiprocessorComputers(并行算法)Two-DimensionalGridExtrusionToolNodalMixedDiscretization(節(jié)點(diǎn)混合離散方法)EmbeddedStructuralLiner(內(nèi)置襯砌單元)擋土墻的模擬On-boardHelpFile(內(nèi)置的幫助系統(tǒng))第19頁(yè)，共74頁(yè)，2023年，2月20日，星期四主要內(nèi)容軟件介紹動(dòng)力分析DynamicOption樁-土相互作用分析Interface隧道分析StructureElement流固耦合分析Fluid-MechanicalInteraction學(xué)習(xí)方法及經(jīng)驗(yàn)介紹第20頁(yè)，共74頁(yè)，2023年，2月20日，星期四完全非線性的動(dòng)力分析與等效線性方法的比較動(dòng)力荷載動(dòng)力邊界條件力學(xué)阻尼與滯回阻尼地震波的調(diào)整動(dòng)孔壓的生成第21頁(yè)，共74頁(yè)，2023年，2月20日，星期四FLAC與等效線性方法等效線性方法是巖土地震工程中的常用方法動(dòng)本構(gòu)模型等效線性模型剪切模量的降級(jí)曲線阻尼比隨剪應(yīng)變的變化FLAC常規(guī)模型(MC)多種頻率成分的干涉和混合永久變形計(jì)算彈塑性計(jì)算第22頁(yè)，共74頁(yè)，2023年，2月20日，星期四動(dòng)力荷載動(dòng)力輸入的類型加速度時(shí)程速度時(shí)程應(yīng)力(壓力)時(shí)程力時(shí)程APPLYINTERIOR(內(nèi)部)TABLEFISH第23頁(yè)，共74頁(yè)，2023年，2月20日，星期四Quiet邊界靜態(tài)(quiet,粘性)邊界LysmerandKuhlemeyer(1969)模型邊界法向和切向設(shè)置獨(dú)立的阻尼器性能對(duì)于法向p波和s波能很好的吸收對(duì)于傾斜入射的波和Rayleigh波也有所吸收，但存在反射人工邊界仍應(yīng)當(dāng)足夠遠(yuǎn)第24頁(yè)，共74頁(yè)，2023年，2月20日，星期四Quiet邊界應(yīng)用內(nèi)部振動(dòng)(如隧道中的列車振動(dòng)問題)動(dòng)力荷載直接施加在節(jié)點(diǎn)上使用Quiet邊界減小人工邊界上的反射不需要FF邊界外部荷載的底部邊界軟土地基上的地震荷載不適合用加速度或速度邊界條件使用應(yīng)力條件t=-2Csrvs地震底部輸入的側(cè)向邊界扭曲了入射波quietquietquiet第25頁(yè)，共74頁(yè)，2023年，2月20日，星期四Free-field邊界Cundalletal.(1980)自由場(chǎng)網(wǎng)格與主體網(wǎng)格的耦合粘性阻尼器，自由場(chǎng)網(wǎng)格的不平衡力施加到主體網(wǎng)格邊界上設(shè)置條件底部水平，重力方向?yàn)閦向側(cè)面垂直，法向分別為x,y向其他邊界條件在APPLYff之前自由場(chǎng)邊界示意圖第26頁(yè)，共74頁(yè)，2023年，2月20日，星期四Free-field邊界APPLYff將邊界上單元的屬性、條件和變量全部轉(zhuǎn)移ff單元上；設(shè)置以后主體網(wǎng)格上的改動(dòng)將不會(huì)被FF邊界所響應(yīng)可存在任意的本構(gòu)模型以及流體耦合(僅豎向)FF邊界進(jìn)行小變形計(jì)算，主體網(wǎng)格可大變形，F(xiàn)F邊界上的變形要相對(duì)較小存在attach的邊界將不能設(shè)置FF邊界邊界上的Interface將不能連續(xù)動(dòng)力邊界設(shè)置需在FF邊界設(shè)置之前第27頁(yè)，共74頁(yè)，2023年，2月20日，星期四力學(xué)阻尼瑞利(rayleigh)阻尼假設(shè)阻尼與質(zhì)量、剛度的線性關(guān)系參數(shù)確定簡(jiǎn)單(等價(jià)平均應(yīng)變=60%*emax)中心頻率(共振計(jì)算，地震平均頻率)臨界阻尼比缺點(diǎn)：計(jì)算速度慢局部(local)阻尼FLAC3D的靜力分析阻尼參數(shù)簡(jiǎn)單適合簡(jiǎn)單情況第28頁(yè)，共74頁(yè)，2023年，2月20日，星期四滯回阻尼(HystereticDamping)模擬巖土介質(zhì)的動(dòng)模量衰減曲線initialdamphystereticnamesig3(3參數(shù))sig4(4參數(shù))Hardin(1參數(shù))(哈丁模型)Default（2參數(shù))優(yōu)點(diǎn)直接采用模量降級(jí)曲線阻尼比不會(huì)影響時(shí)間步缺點(diǎn)輸出的曲線會(huì)不一致第29頁(yè)，共74頁(yè)，2023年，2月20日，星期四滯回阻尼(HystereticDamping)低循環(huán)應(yīng)變下得到的阻尼比要小于試驗(yàn)結(jié)果，這會(huì)導(dǎo)致低級(jí)的噪聲，尤其在高頻情況下。可以在中心頻率上增加一個(gè)小量的Rayleigh阻尼(~0.2%剛度比例)，這樣也不會(huì)降低時(shí)步；若初始應(yīng)力不為0，剪應(yīng)力-剪應(yīng)變曲線可能不匹配。因此在生成初始應(yīng)力時(shí)就要調(diào)用Hyst阻尼；Hyst阻尼不僅會(huì)增加能量損失，還會(huì)導(dǎo)致在大循環(huán)應(yīng)變下的平均剪切模量的降低，在輸入波的基頻接近共振頻率的時(shí)候，可能會(huì)導(dǎo)致動(dòng)力反應(yīng)幅值的增大；Hyst阻尼之前要做一次彈性無(wú)阻尼求解，以獲得發(fā)生循環(huán)應(yīng)變的最大水平，若循環(huán)應(yīng)變過大導(dǎo)致剪切模量過多的降低，那么用Hyst阻尼是有問題的；即使應(yīng)變較小，使用屈服模型也會(huì)增大應(yīng)變，因此若有廣泛屈服的現(xiàn)象，則使用屈服模型，不用Hyst阻尼參考了Itasca的中國(guó)培訓(xùn)資料第30頁(yè)，共74頁(yè)，2023年，2月20日，星期四地震波的調(diào)整基線校正對(duì)于地震分析的加速度時(shí)程，其積分得到的速度和位移應(yīng)歸0美國(guó)地質(zhì)調(diào)查研究所BasicStrong-MotionAccelerogramProcessingSoftware(BAP)對(duì)網(wǎng)格施加一個(gè)固定速度從而使殘余的位移變?yōu)?動(dòng)力荷載的頻率與單元尺寸的雙向調(diào)整高頻的輸入要求單元尺寸很小一定的單元尺寸對(duì)應(yīng)輸入的最大頻率一般進(jìn)行濾波處理濾掉低能量的高頻FFT.FISOriginSeismoSignal第31頁(yè)，共74頁(yè)，2023年，2月20日，星期四地震波的調(diào)整El-Centro波FFT修正后的時(shí)程修正后FFT5Hz第32頁(yè)，共74頁(yè)，2023年，2月20日，星期四動(dòng)孔壓的生成——液化干沙剪應(yīng)變循環(huán)加載試驗(yàn)初始加載階段，沙土通常先壓實(shí)再膨脹。卸載時(shí)，沙土遵循與加載相似的路徑，但在零應(yīng)變時(shí)，有些殘余體積應(yīng)變存在。取決于初始孔隙率，這可能代表純粹的壓實(shí)假定孔隙中充滿水對(duì)于常體積測(cè)試，有效應(yīng)力降低，孔隙水壓保持不變對(duì)于常荷載測(cè)試，(例如，盒子上法向荷載固定)，孔隙水壓增加，有效應(yīng)力減小有效應(yīng)力為零時(shí)發(fā)生液化第33頁(yè)，共74頁(yè)，2023年，2月20日，星期四動(dòng)孔壓的生成——液化因此孔隙水壓增加不是液化的基本原因由于顆粒間(重組以后)的低接觸力導(dǎo)致有效應(yīng)力的減小描述液化的模型高級(jí)模型：BSHP(邊界面低塑性本構(gòu)模型,Wangetal.1990)簡(jiǎn)單模型：MC+體積應(yīng)變?cè)隽磕Ｐ虵inn模型：Byrne模型：第34頁(yè)，共74頁(yè)，2023年，2月20日，星期四主要內(nèi)容軟件介紹動(dòng)力分析DynamicOption樁-土相互作用分析Interface隧道分析StructureElement流固耦合分析Fluid-MechanicalInteraction學(xué)習(xí)方法及經(jīng)驗(yàn)介紹第35頁(yè)，共74頁(yè)，2023年，2月20日，星期四Interface單元簡(jiǎn)介用途巖體介質(zhì)中的解理、斷層、巖層面地基與土體的接觸箱、槽及其內(nèi)充填物的接觸空間中無(wú)變形的固定“障礙”原理三角形單元(無(wú)厚度!)8參數(shù)三種工作模式粘結(jié)界面粘接滑移庫(kù)倫滑動(dòng)第36頁(yè)，共74頁(yè)，2023年，2月20日，星期四Interface的建立(1)關(guān)鍵要形成同一位置的兩個(gè)節(jié)點(diǎn)(面)“移來移去”(推薦)建兩個(gè)分開的模型建立接觸單元通過INI*add使模型接觸注意dist的含義NOmerge,NOattach!!!接觸面dist1234第37頁(yè)，共74頁(yè)，2023年，2月20日，星期四Interface的建立(2)“導(dǎo)來導(dǎo)去”利用expgrid,impgrid命令進(jìn)行網(wǎng)格導(dǎo)出與導(dǎo)入配合DELETE命令適于內(nèi)部接觸面的建立，或其他前處理工具建立的網(wǎng)格第38頁(yè)，共74頁(yè)，2023年，2月20日，星期四Interface的建立(2)save1.savdelrangrop2notInterface1facesave2.savrest1.sav

delrangroup2

expgrid1.fac3drest2.sav

impgrid1.flac3d

第39頁(yè)，共74頁(yè)，2023年，2月20日，星期四Interface的建立(3)GENseparateINTERFACEwrap指定正確的group第40頁(yè)，共74頁(yè)，2023年，2月20日，星期四Interface參數(shù)的確定虛構(gòu)的為了合并節(jié)點(diǎn)而設(shè)置的接觸面kn=ks=10*真實(shí)的剛性接觸面如料倉(cāng)下料c,D,Tension重要，kn,ks不重要真實(shí)的柔性接觸面斷層；水力劈裂材料試驗(yàn)得到參數(shù)對(duì)于kn,ks：巖石斷層10~100MPa/m(粘土);100GPa(巖石)反分析方法：通過斷層中巖石的變形與原巖的變形第41頁(yè)，共74頁(yè)，2023年，2月20日，星期四例3：?jiǎn)螛冻休d力分析軟土地基bulk1.6878E6shear3.6167E5coh15E3fric12dens1.73E3樁體bulk5e9shear3.75e9dens2.5e30.5m8m10m20m第42頁(yè)，共74頁(yè)，2023年，2月20日，星期四計(jì)算過程施加樁頂荷載計(jì)算結(jié)果第43頁(yè)，共74頁(yè)，2023年，2月20日，星期四影響因素水平因素knksfriccoh1100e100e203.00E+04210k10g0.7friccoh/0.73kgfric7.50E+03第44頁(yè)，共74頁(yè)，2023年，2月20日，星期四計(jì)算工況設(shè)計(jì)

水平因素工況knksfriccoh(1)1111(2)1222(3)1333(4)2123(5)2231(6)2312(7)3132(8)3213(9)3321第45頁(yè)，共74頁(yè)，2023年，2月20日，星期四計(jì)算結(jié)果Ks取1Ks取2Ks取3水平因素knksfriccohRa差異度(1)111137.8%(2)122242.9%(3)133396.9%(4)212337.8%(5)22317.1%(6)231296.9%(7)313254.1%(8)321322.4%(9)332196.9%第46頁(yè)，共74頁(yè)，2023年，2月20日，星期四最優(yōu)方案水平因素knksfriccohk159.2%43.2%52.4%47.3%k247.3%24.1%59.2%64.6%k357.8%96.9%52.7%52.4%極差11.9%72.8%6.8%17.3%最優(yōu)方案2211102211112222第47頁(yè)，共74頁(yè)，2023年，2月20日，星期四合理步驟單樁分析簡(jiǎn)單網(wǎng)格接觸面參數(shù)多次試算理論、實(shí)測(cè)加密網(wǎng)格接觸參數(shù)理想結(jié)果群樁分析Pile結(jié)構(gòu)單元單元參數(shù)理想結(jié)果Pile結(jié)構(gòu)單元單元參數(shù)理想結(jié)果第48頁(yè)，共74頁(yè)，2023年，2月20日，星期四主要內(nèi)容軟件介紹動(dòng)力分析DynamicOption樁-土相互作用分析Interface隧道分析StructureElement流固耦合分析Fluid-MechanicalInteraction學(xué)習(xí)方法及經(jīng)驗(yàn)介紹第49頁(yè)，共74頁(yè)，2023年，2月20日，星期四FLAC3D中的結(jié)構(gòu)單元有限單元梁(beam)單元錨索(cable)單元樁(pile)單元錨桿:rockbolt殼(shell)單元格柵(geogrid)單元土工織物；土工格柵初襯(liner)單元beamcablepileshellgeogridliner第50頁(yè)，共74頁(yè)，2023年，2月20日，星期四結(jié)構(gòu)單元的應(yīng)用土與結(jié)構(gòu)的相互作用樁基；基坑；邊坡錨固地下硐室的支撐結(jié)構(gòu)；采礦；盾構(gòu)土工織物；土工合成材料結(jié)構(gòu)不宜復(fù)雜巖土工程軟件，不宜單純的結(jié)構(gòu)分析復(fù)雜結(jié)構(gòu)的模擬很困難結(jié)構(gòu)單元仍不完善plot顯示雙向接觸結(jié)構(gòu)(擋土墻)結(jié)構(gòu)單元的厚度第51頁(yè)，共74頁(yè)，2023年，2月20日，星期四Liner結(jié)構(gòu)單元三節(jié)點(diǎn)扁平有限單元每個(gè)節(jié)點(diǎn)有6個(gè)自由度3個(gè)移動(dòng)，3個(gè)旋轉(zhuǎn)能夠抵抗膜及彎矩荷載能夠承受主方向的拉壓應(yīng)力能夠模擬管片與土體之間的分離及隨后的重新接觸能夠模擬管片與土體之間的摩擦相互作用法向切向第52頁(yè)，共74頁(yè)，2023年，2月20日，星期四例4：隧道與土體的相互作用半圓隧道直徑3.25m上覆土層厚度5m計(jì)算范圍3r土體彈性計(jì)算(K=30MPa,G=10MPa)參數(shù)化編程幾何尺寸模型參數(shù)網(wǎng)格形狀hthbrB第53頁(yè)，共74頁(yè)，2023年，2月20日，星期四計(jì)算步驟模型網(wǎng)格初始應(yīng)力生成施加管片計(jì)算結(jié)果第54頁(yè)，共74頁(yè)，2023年，2月20日，星期四管片的連接冷連接彎矩和剪力不能直接在環(huán)與環(huán)間傳遞，只能通過其相鄰的介質(zhì)傳遞全連接相鄰的Liner單元在連接處共用一個(gè)節(jié)點(diǎn)，連接處重疊單元不能發(fā)生移動(dòng)或旋轉(zhuǎn)結(jié)點(diǎn)連接即結(jié)點(diǎn)間的連接在6個(gè)方向的自由度上用彈簧來模擬，每個(gè)自由度都可具有一定的特性通縫拼接錯(cuò)縫拼接第55頁(yè)，共74頁(yè)，2023年，2月20日，星期四連接方式的影響變形結(jié)果zr屬性更改6屬性更改注：變形放大200倍第56頁(yè)，共74頁(yè)，2023年，2月20日，星期四例5：結(jié)構(gòu)的動(dòng)力響應(yīng)configdynselpileid=1beg000end001selpilepropdens2400&Emod1.0e10Nu0.3XCArea0.3&XCJ0.16375XCIy0.00625XCIz0.01575&Per2.8CS_sK1.3e11CS_nK1.3e11&CS_nGapoffselnodefixxyzxryrzrranid=1selsetdampcombineddeff1whilesteppingf0=10000*sin(10*dytime)np=nd_headloopwhilenp#nullifnd_pos(np,1,3)=1nd_apply(np,1)=f0endifnp=nd_next(np)endloopendsolveage1pile10000*sin(10*t)第57頁(yè)，共74頁(yè)，2023年，2月20日，星期四例5：結(jié)構(gòu)的動(dòng)力響應(yīng)第58頁(yè)，共74頁(yè)，2023年，2月20日，星期四主要內(nèi)容軟件介紹動(dòng)力分析DynamicOption樁-土相互作用分析Interface隧道分析StructureElement流固耦合分析Fluid-MechanicalInteraction學(xué)習(xí)方法及經(jīng)驗(yàn)介紹第59頁(yè)，共74頁(yè)，2023年，2月20日，星期四FLAC3D流-固耦合分析(單相流)基本功能理論框架計(jì)算模式滲流邊界條件，初始條件單滲流計(jì)算及滲流耦合計(jì)算第60頁(yè)，共74頁(yè)，2023年，2月20日，星期四基本功能滲流各向同性、各向異性不同的滲流模型和屬性流體壓力，涌入量，滲漏量和不滲水邊界抽水井、點(diǎn)源、體積源飽和滲流可采用顯式差分法、隱式差分法非飽和滲流采用顯式差分法滲流-固體-熱的耦合流體和固體的耦合程度依賴于土體顆粒(骨架)的壓縮程度，用Biot系數(shù)表示顆粒的可壓縮程度循環(huán)荷載引起的動(dòng)水壓力變化和土體液化第61頁(yè)，共74頁(yè)，2023年，2月20日，星期四地下水模擬方法有效應(yīng)力計(jì)算不耦合孔壓為了正確計(jì)算有效應(yīng)力滲流計(jì)算已得到孔壓分布飽和有水面線的部分飽和力學(xué)變形產(chǎn)生孔壓靜力或者動(dòng)力不排水孔壓或液化流固耦合計(jì)算模式時(shí)間比例第62頁(yè)，共74頁(yè)，2023年，2月20日，星期四有效應(yīng)力計(jì)算不設(shè)置CONFIGFluid孔壓不改變?cè)O(shè)置孔壓分布INITIALppWATERtableWATERdensityWATERtablefaceSETgravity手動(dòng)設(shè)置干濕密度設(shè)置CONFIGfluid設(shè)置土體干密度

滲流模型MODELfl_isotropicMODELfl_anisotropicMODELfl_nullSETfluidoffsetWATERbulk=0無(wú)滲流模式滲流模式第63頁(yè)，共74頁(yè)，2023年，2月20日，星期四不耦合計(jì)算CONFIGfluidSETmechoff正確的滲透系數(shù)得到孔壓分布和水面線穩(wěn)態(tài)滲流可以減小KfCONFIGfluidSETfluidoff正確的流體模量Kf不需要滲透系數(shù)單滲流模式單力學(xué)模式第64頁(yè)，共74頁(yè)，2023年，2月20日，星期四流固耦合計(jì)算CONFIGfluidSETfluidon真實(shí)的流體模量Kf和滲透系數(shù)耦合方式DpDevDsDevDp計(jì)算模式手動(dòng)調(diào)整的STEP求解主從進(jìn)程的SOLVE求解自動(dòng)STEP求解第65頁(yè)，共74頁(yè)，2023年，2月20日，星期四流固耦合的計(jì)算方法手動(dòng)調(diào)整的STEP求解

SETfluidonmechoff STEP SETfluidoffmechon STEP 主從進(jìn)程的SOLVE求解SETmechforceSETmechsubstepnauto (從進(jìn)程)SETmechsubstepm (主進(jìn)程)SOLVEage自動(dòng)STEP求解STEP第66頁(yè)，共74頁(yè)，2023年，2月20日，星期四滲流問題(CONFIGfluid)分析步驟時(shí)間比例(ts,t