子步的案例
workbench分析如何設置載荷步與子步
下面通過一個例子進行說明以上三個參數(shù)的關系和含義，例如總載荷為10000N，初始子步為10，最小子步為5，最大子步為50，則程序求解時，首先在10000/10=1000N求解，該子步求解完成后，第二個子步，程序根據(jù)內(nèi)部評估，增量范圍是10000/5=2000到10000/50=200，之間變化，則第二個子步位置的載荷為1000+（2000~200）之間。
總之，子步是為了提高收斂性而設置的參數(shù)，一般情況，子步數(shù)量越大，則越有利于收斂。但是過多的子步會增加求解時間和占用更多的硬盤空間，可以適當調(diào)試。
展開 多分析步下的裂紋擴展
一個使用所有子步的瞬態(tài)事件載荷計劃將不會受到影響。然而，如果選擇了一個特定的子步，圖10.42，而這個子步在每個裂紋擴展步中都不存在，那么裂紋擴展和/或疲勞周期計數(shù)就會受到影響；下一節(jié)節(jié)描述了一個這方面的例子。
圖10.41 使用所有分析步和子步的瞬態(tài)載荷計劃事件。
圖10.42 為一個分析步選擇特定的子步。
在上面的例子中，在定義載荷計劃之前，對初始裂紋進行了靜態(tài)裂紋分析。因此，F(xiàn)RANC3D已經(jīng)從.dtp文件中讀取了結(jié)果，并知道了子步。
如果用戶不做靜態(tài)分析，F(xiàn)RANC3D就不會對子步有任何了解。如果試圖在插入裂紋后立即進行自動裂紋擴展分析，F(xiàn)RANC3D將呈現(xiàn)疲勞載荷計劃對話框，只識別出分析步，圖10.43。子步標記為"---"，表示子步是未定義的。
圖10.43 選擇分析步對話框-子步ID未定義。
通常建議在插入裂紋后做靜態(tài)裂紋分析，因為這將使用戶能夠確保裂紋模型的結(jié)果與無裂紋模型的結(jié)果一致，并提供分析步的子步SIFs。
無子步的分析步
使用上一節(jié)中的模型，完成了六步裂紋擴展。定義了一個簡單的載荷計劃，使用所有三個分析步的最后一個子步的總和，圖10.44。
展開 ANSYS分析類型與求解器控制選項（1）
子步數(shù)和時間步長
命令：NSUBST, NSBSTP, NSBMX, NSBMN, Carry
NSBSTP - 當前荷載步的子步數(shù)。如果使用了自動時間步（即 AUTOTS,ON）則該子步數(shù)僅用于第一子步，也即第一子步的荷載增量用 NSBSTP 求得，其余子步的荷載增量由程序自動確定。
NSBMX - 當 AUTOTS 打開時，NSBMX 為最大子步數(shù)。
NSBMN - 當 AUTOTS 打開時，NSBMN 為最小子步數(shù)。
Carry - 時間步長繼承控制參數(shù)，其值可取：
=OFF：使用 NSBSTP 確定每個荷載步開始時的時間步長；
=ON：如果 AUTOTS 打開，使用前一荷載步的最后時間步長作為該荷載步開始的時間步長。
該命令中的 NSBSTP 參數(shù)用于確定在當前荷載步內(nèi)，每個子步（或時間步）荷載增量的大?。ㄐ逼潞奢d，如為階躍荷載則一個子步到全值）。
最小和最大子步數(shù)在采用自動時間步時，影響結(jié)果點的多少和收斂控制。例如問題容易收斂，程序會采用較小的子步數(shù)（時間步長大，荷載增量大）得到的結(jié)果點就少。如果問題收斂困難，程序會采用較大的子步數(shù)（時間步長小，荷載增量?。?，可得到較多的結(jié)果點；但是如果問題特別難以收斂，程序會采用最大子步數(shù)（最小時間步長）求解以獲得收斂結(jié)果，
通過平衡迭代一定次數(shù)后（NEQIT 命令設置）仍然不能收斂，則程序判定為不收斂并結(jié)束求解。
建議對該命令的各個參數(shù)都要設置，但是對于一類問題設置多大的數(shù)目是合適的呢？這點只能靠求解控制經(jīng)驗或試算確定。一般可采用缺省的設置選項，不能收斂時可不斷調(diào)整參數(shù)并逐步逼近收斂。
展開 ANSYS Workbench結(jié)構(gòu)瞬態(tài)分析
圖2
載荷步設置、求解器設置、重啟設置、非線性設置、輸出設置、阻尼設置、分析數(shù)據(jù)管理、可視化。
（1）載荷步設置
圖3
Number Of Steps：載荷步，表示分幾步施加載荷。
Current Step Number：當前載荷步。
Step End Time：當前載荷步結(jié)束時間。
Auto Time Stepping：自動時間步長是否打開，這個大多數(shù)情況下打開，讓程序自動決定迭代時間。
Define By：定義載荷子步的方式?？梢酝ㄟ^時間和載荷子步數(shù)來定義，用時間定義的意思是每一載荷子步經(jīng)歷的時間是多少，用載荷子步定義的意思是一個載荷步有多少個載荷子步。
Initial Time Step：初始載荷子步的時間，也就是計算時第一個載荷子步的時間。
Minimum Time Step：最小的載荷子步時間
Maximum Time Step：最大的載荷子步時間
Time integration：時間積分是否打開，如果打開，表示考慮時間對計算結(jié)果的影響，也就是考慮動力響應，如果不打開，那么計算過程中相當于是準靜態(tài)過程。
本次例子定義三個載荷子步，分別施加壓力載荷，如圖所示。
圖4
每個載荷步都可以通過以上選項設置不同的載荷子步。
（2）求解器設置和重啟動設置
圖5
Solver Type：求解器類型，有直接法和迭代法兩種，這個具體的意義以前文章有介紹，這里不說。
Weak Springs：弱彈簧效應，相當于施加一個地面與模型的一根弱彈簧，目的是方便數(shù)值計算。
Large Deflection：大變形是否打開，如果有非線性的大變形，可以打開。
重啟動部分不做介紹。
（3）非線性設置
圖6
Newton-Raphson：牛頓法迭代設置，有完全法等，具體的含義建議自己查資料。
展開 
ANSYS Workbench非線性分析收斂曲線解讀
Solution Output選項
力收斂曲線如下圖所示:
力收斂曲線圖
判斷收斂的方法很簡單，只要“計算的力收斂曲線”落在“力收斂準則”曲線之下，就表示該載荷步或子步收斂了。
該模型中有兩個載荷步，分析設置中時間步長設置為“Program Contrlled”.
除了看上述的力收斂曲線圖，我們可以設置“Solution Output= Solve Output”查看計算輸出信息，從其中可以更詳細地看到收斂情況。
可以將計算輸出的信息與力收斂曲線圖對比起來看，就更容易理解力收斂圖了。
第1個載荷步中，第1個分析子步經(jīng)過了15次迭代收斂（圖中每個圓點代表一次迭代）。
經(jīng)過4個分析子步，第1個載荷步完成加載并收斂。第2個載荷步程序自動設置的信息如下：
初始子步數(shù)量為5，載荷步的分析時間為1s，因此初始的時間步長為0.2s。
第2個載荷步的第1個分析子步，經(jīng)過25次計算迭代后，還不收斂。程序進行自動二分，將時間步長除以2，變?yōu)?.1s。
自動二分是一種用于解決非線性分析過程中收斂困難的策略。當收斂失敗發(fā)生在某個子步中，程序會自動減小時間步長，通常是前一個步長的一半左右。然后，程序會從前一個成功收斂的時間子步繼續(xù)求解。如果再次遇到收斂失敗，程序會繼續(xù)減小時間步長并繼續(xù)求解，直到達到收斂或達到指定的最小時間步長值。這種方法有助于逐步逼近正確解，并確保分析的穩(wěn)定性和準確性。
第2載荷步的第4個子步中，進行了18次迭代未收斂，預測需要50次迭代，超過了程序允許的25次，再次進行二分，將時間步長改為0.05s。
在第5個子步采用0.05s的時間步長，經(jīng)過3次迭代收斂了，程序認為，可以加大一點時間步長，自動改成了0.75s，增長比例為1.5。
展開 從不收斂的結(jié)果中識別正確塑性極限載荷
第1載荷步
AUTOTS,ON !自動載荷步
NSUBST,20,200,5 !載荷子步
OUTRES,ALL,1 !輸出每個子步的結(jié)果
FK,2,MZ,39.270e3 !施加彈性極限扭矩
SOLVE ! 第1載荷步求解
TIME,2 !第2載荷步
NSUBST,20,200,5 !載荷子步
OUTRES,ALL,1 !輸出每個子步的結(jié)果
FK,2,MZ,1.2*52.360e3 !1.2倍塑性極限扭矩
SOLVE ! 第2載荷步求解
FINISH
/POST1
SET,LIST !查看計算的載荷步
SET,LAST !讀入最后載荷子步的結(jié)果
SET,PREVIOUS !提取上一個載荷子步的結(jié)果
PLNSOL,S,EQV,0,1.0 !第四強度SEQV
PLNSOL,S,INT,0,1.0 !第三強度SINT
PLNSOL,S,XY,0,1.0 !切應力SXY
PLNSOL,S,XZ,0,1.0 !切應力SXY
轉(zhuǎn)載自好學ANSYS，詳細內(nèi)容請移步公眾號哦！鏈接：https://mp.weixin.qq.com/s/WJqoXKwTuolAfXuC1a7-mg
展開 『分享』ANSYS結(jié)構(gòu)分析基礎知識
· 修正的（NROPT，MODI）：程序使用修正的牛頓－拉普森方法，在這種方法中正切剛度矩陣在每一子步中都被修正。在一個子步的平衡迭 代期間矩陣不被改變。這個選項不適用于大變形分析。自適應下降是不可用的 · 初始剛度（NROPT，INIT）：程序在每一次平衡迭代中都使用初始剛度矩陣這一選項比完全選項似乎較不易發(fā)散，但它經(jīng)常要求更多次的迭代來得到收斂。它不適用于大變形分析。自適應下降是不可用的。
選項：方程求解器 J對于非線性分析，使用前面的求解器（缺省選項）。 3、在模型上加載，記住在大變型分析中慣性力和點載荷將保持恒定的方向，但表面力將“跟隨”結(jié)構(gòu)而變化。 4、指定載荷步選項。這些選項可以在任何載荷步中改變。下列選項對非線性靜態(tài)分析是可用的：
普通選項
普通選項包括下列：
· Time(TIME)
ANSYS程序借助在每一個載荷步末端給定的TIME參數(shù)識別出載荷步和子步。使用TIME命令來定義受某些實際物理量（如先后時間，所施加的壓力，等等。）限制的TIME值。程序通過這個選項來指定載荷步的末端時間。
注意──在沒有指定TIME值時，程序?qū)⒁罁?jù)缺省自動地對每一個載荷步按1.0 增加TIME（在第一個載荷步的末端以TIME=1.0開始）。
·時間步的數(shù)目〔NSUBST〕
·時間步長〔DELTIM〕非線性分析要求在每一個載荷步內(nèi)有多個子步（或時間步；這兩個術(shù)語是等效的）從而ANSYS可以逐漸施加所給定的載荷，得到精確的解。NSUBST和DELTIM命令都獲得同樣的效果（給定載荷步的起始，最小，及最大步長）。NSNBST 定義在一個載荷步內(nèi)將被使用的子步的數(shù)目，而DELTIM明確地定義時間步長。如果自動時間步長是關閉的，那么起始子步長用于整個載荷步。缺省時是每個載荷步有一個子步。
展開 案例27-帶有3D網(wǎng)格重劃分的結(jié)構(gòu)鋼熱軋分析
重啟動文件在每個子步中保存，因為重劃分需要重啟動文件，而需要重劃分的子步仍然未知。結(jié)果項也存儲在每個子步中。
出于說明的目的，重啟動文件和結(jié)果項在每個子步都會保存，盡管這樣做需要大量內(nèi)存。在大多數(shù)情況下，只需每隔幾個子步保存一次即可。
以下輸入應用運行的求解設置：
第一步加載在時間time=0.7718開始不收斂，最后一個收斂的子步第40子步如下圖：
因為網(wǎng)格重劃分應該從網(wǎng)格出現(xiàn)變形之前的幾步開始，所以從第30個子步開始啟動網(wǎng)格重劃分是合理的。
下圖顯示了原始變形網(wǎng)格和新的良好網(wǎng)格
以下輸入在第30子步啟動重劃分，并讀取新網(wǎng)格：
在網(wǎng)格重劃分之后，程序?qū)⒈砻孑d荷，力，邊界條件和接觸對（如果存在的話）從原來的變形網(wǎng)格傳遞給新網(wǎng)格，檢查網(wǎng)格重劃分后的模型來驗證所有的數(shù)據(jù)已經(jīng)成功轉(zhuǎn)移是一個好習慣。
從原始網(wǎng)格映射到新網(wǎng)格的求解項
在網(wǎng)格重劃分之后，所有的結(jié)果被映射到新網(wǎng)格，造成的殘余力重新平衡，通過mpsolve來實現(xiàn)映射。下圖顯示了網(wǎng)格重劃分前后的接觸壓力。
映射操作結(jié)束了重劃分過程，標準的多幀重啟動使用新網(wǎng)格恢復求解過程。
使用新網(wǎng)格繼續(xù)分析
在將變量從舊網(wǎng)格映射到新網(wǎng)格并重新平衡殘余力后，多幀重啟動（ANTYPE,,restart,,,CONTINUE）使用新網(wǎng)格恢復非線性求解。
下列命令會讓分析重啟動：
結(jié)果和討論
下圖顯示了最后一個收斂子步驟（TIME=.7718）初始運行的變形圖（USUM）。
成功重劃分后，下圖顯示了建立軋制過程后（第一個加載步結(jié)束時）模型的變形圖（USUM）：
在第二個加載步中，通過旋轉(zhuǎn)輥并使用輥與坯料之間的高摩擦接觸來執(zhí)行熱軋過程。
展開 強度折減系數(shù)法在ａｎｓｙｓ中的命令流
讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖
!邊坡在強度折減系數(shù)F=1.2時結(jié)果分析
Resume,'F1.2','db' !讀入邊坡在強度折減系數(shù)F=1.2時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !繪制邊坡模型水平方向位移云圖
plnsol,eppl,eqv !繪制邊坡模型塑性應變云圖
!邊坡在強度折減系數(shù)F=1.4時結(jié)果分析
Resume,'F1.4','db' !讀入邊坡在強度折減系數(shù)F=1.4時
set,1,last !讀入后一個子步
pldisp,1 !繪制邊坡模型變形圖
plnsol,u,x !
展開 ansys非線性瞬態(tài)結(jié)構(gòu)分析重要命令
3.2 選擇積分時間步的準則
瞬時分析的精度取決于時間步的大小。太大的時間步可能引起不可接受的誤差，太小的時間步浪費計算時間和資源。
4
非線性瞬態(tài)熱應力分析中的重要命令
①輸出控制（結(jié)果輸出到數(shù)據(jù)庫），建立存儲規(guī)格。
間接法計算熱應力時，熱分析的結(jié)果文件要作為結(jié)構(gòu)分析的熱載荷輸入。因此，熱分析的載荷步時間步、結(jié)果存儲設置要適應結(jié)構(gòu)分析。
OUTRES, Item, Freq, Cname
Item: NSOL,節(jié)點結(jié)果；ESOL,單元結(jié)果；ALL，所有。
Freq: n,每第n個子步；-n,均分成n段；NONE,一個也不存；ALL,每一子步；LAST,最后一子步；%array%按數(shù)組提供的時刻來存儲。
例1：
NSUBST,6
OUTRES,ERASE 設置到默認值，對于靜態(tài)和瞬態(tài)分析，默認的是輸出每一載荷步的最后子步的所有結(jié)果；諧態(tài)分析是每一子步。
例2：
NSUBST,6
OUTRES,NSOL,2 每2子步保持節(jié)點結(jié)果，其他不保存
例3：
NSUBST,6
OUTRES,ESOL,4 第4子步和第6子步（last），保存單元結(jié)果
②保存
SAVE, Fname, Ext, --, Slab
Slab：ALL，保存所有數(shù)據(jù)（模型數(shù)據(jù)，結(jié)果數(shù)據(jù)，后處理數(shù)據(jù)），默認值；MODEL，模型數(shù)據(jù)；SOLU，結(jié)果數(shù)據(jù)。
③自動時間步
AUTOTS,ON
DELTIM, DTIME, DTMIN, DTMAX, Carry
如果使用自動時間步，當Carry=OFF，以DTIME為起始時間步長，最小時間步不小于DTMIN，最大不大于DTMAX，當Carry=ON，以上一載荷步的最后子時間步長為起始時間步長。
展開 4_APDL基礎及仿真理論_分析重啟動
ANSYS命令流學習筆記4
重啟動技術(shù)：
和多載荷步一樣能實現(xiàn)不同載荷的先后加載。
第一次分析完成后，想把更多的載荷步加入到分析中，或者在線性分析中加入別的加載條件，或者在瞬態(tài)分析中加入另外的時間里程加載曲線，或者在非線性分析收斂失敗時需要恢復。
多點重啟動，可以保存多個子步的結(jié)果文件，可以從任何一個保存的子步結(jié)果重啟動運算。
WB中也可以方便地完成重啟動。
簡單命令流案例如下：
finish
/clear
/prep7
et,1,42
mp,ex,1,2.1e5
mp,nuxy,1,0.3 !單元及材料屬性
rect,0,100,0,100 !建模
esize,10,
amesh,all !劃分網(wǎng)格
dl,1,all,all
sfl,3,pres,-100
/pbc,all,1 !邊界條件
finish
/solu
nlgeom,on !打開大變形
nsubst,4 !設置4子步
antype,static,new !分析類型結(jié)構(gòu)靜力求解
outres,all,all !記錄每一步結(jié)果
autots,off !關閉自動步長
time,1 !求解時間1s
rescontrol, ,all,1,4 !重啟動設置，每一子步保存文件，共4個重啟動文件?？梢詮倪@4個文件中，任意一個，重新啟動運算。
Solve !可以通過時間里程后處理，查看載荷步對結(jié)果的影響
finish
/solu
antyp, ,rest,1,2 !從第二步開始重新加載
sfl,3,pres,-200 !線載荷壓力改為-200
solve !可以通過時間里程后處理，查看載荷步對結(jié)果的影響，對比不同
finish
展開 
如何在ANSYS WORKBENCH中施加分段函數(shù)激勵
3.劃分網(wǎng)格
4.分析設置
設置兩個時間步，
第一步終止時間為1秒，打開自動時間步長，通過載荷步來定義載荷子步，初始子步10步，最小5步，最多20步。
再定義第二步如下
其含義是
第2步終止時間為2秒，打開自動時間步長，通過載荷步來定義載荷子步，初始子步10步，最小5步，最多20步.
5.固定左端
6.在上面施加分布載荷1
首先定義第一個載荷步內(nèi)的函數(shù)載荷
接著休眠期第二段（1-2秒內(nèi)的部分）
得到結(jié)果如下
7.在上面施加分布載荷2
接著休眠期第1段（0-1秒內(nèi)的部分）
得到結(jié)果如下
這就可以了。
至于后面的求解就不再贅述了。
來源：宋博士的博客，版權(quán)歸作者所有。
展開 橡膠的密封性計算
與常規(guī)的壓力載荷相比，施加到子步上的流體壓力大小取決于壓力是常數(shù)輸入還是表格輸入，同時還取決于是階躍加載還是斜坡加載。
如果瞬時施加所有的流體壓力，則會造成在靠近接觸界面附近產(chǎn)出較大的應力梯度，從而造成計算收斂困難。同理，當流體壓力突然被移除或突然中斷流體滲透，也會導致計算收斂困難。為了穩(wěn)定計算并提高計算的收斂性，ANSYS提供了斜坡施加流體壓力的選項，該選項在一個求解周期中通過多個計算子步線性的增加流體壓力進行加載。為了實現(xiàn)斜坡加載，可以使用接觸單元實常數(shù)FPAT指定流體滲透作用時間。輸入正值表示滲透作用時間等于輸入的時間乘以載荷步的增量時間；輸入負值表示滲透作用時間就是該數(shù)值的絕對值。ANSYS的默認值為0.01，則默認的滲透作用時間等0.01乘以當前載荷步的時間增量。
在每一個滲透點，如果當前子步的增量時間小于流體滲透作用時間，即FPAT>(tn-tn-1)，則流體壓力在滲透作用時間周期中，從前一個子步施加的壓力線性變化增加到總的流體壓力，如圖2-3所示。如果FPAT<(tn-tn-1)，則施加當前全部的流體壓力如圖2-4所示。在壓力滲透關閉點，如果當前子步的時間增量小于流體滲透作用時間，則流體壓力在滲透作用時間周期中從前一個子步施加的壓力值線性減小到0，否則流體壓力將被瞬時刪除。
6. 工程實例-O型橡膠密封裝置的壓力滲透有限元分析
1-材料參數(shù)
活塞和缸套的材料為結(jié)構(gòu)鋼， 0型密封圈的材料為橡膠，材料模型參數(shù)如下：Neo-
Hookean模型，Mu=20E6Pa，
D1=0.015.
2-邊界條件
完全固定約束缸套的底邊和側(cè)邊，在活塞頂面上施加向下的位移0.007m，然后第二個載荷步施加流體穿透壓力20MPa。
展開 網(wǎng)格自適應技術(shù)！
6.分析設置
注意：子步的設置之前的幾篇文章就已經(jīng)講過，子步可以設置大，但不能無限大，本例最大設置為10000，讀者可以通過使用不同的Max Substep對比一下對收斂情況的影響。對于超彈性本構(gòu)必須打開大變形開關，否則無法計算。
7.求解
注意：求解至0.57484s時候，出現(xiàn)不收斂情況，也可以考慮打開重啟動，打開方式如下圖所示，然后當不收斂情況出現(xiàn)時，在分析設置中選擇要重啟的子步，然后打開穩(wěn)定性開關，設置如圖所示，點擊Solve繼續(xù)求解，一般情況下，重啟動對收斂幫助很大，但是遇到較強的非線性，重啟動有時也作用不大。
注意：重啟點需要往前移幾個子步，比如在第83個子步處不收斂了，重啟點不能設置在從子步83處重啟，必須設置為83的前移數(shù)字，比如80、78等等，否則也無法收斂。
注意：當使用重啟動功能時，再次點擊Solve，軟件會根據(jù)設定的重啟點開始接著上回的計算繼續(xù)，這不僅是節(jié)約時間的辦法，也是幫助收斂的功能之一。
8.查看不收斂前的變形情況
此上案例是設定橡膠參數(shù)，通過壓頭下壓來使得正方體變形，但是由于大變形問題，因此存在不收斂問題，計算至0.57484s，計算自動終止。
案例二：通過調(diào)用U-P（雜交單元）解決橡膠的大應變問題
由于橡膠表現(xiàn)為不可壓縮性，泊松比近似為0.5，易發(fā)生體積自鎖，因此需要使用U-P單元技術(shù)來解除體積自鎖。
在正方體的Body下插入一條command：Keyopt，matid，6,1代表的是打開U-P雜交單元選項。
展開 ABAQUS中的瞬態(tài)滲流和穩(wěn)態(tài)滲流 ￥10
（附件slpoe_allfix）
瞬態(tài)滲流不考慮固結(jié)沉降
（2）瞬態(tài)滲流考慮固結(jié)沉降時（采用Soil，Transient分析步，只約束邊界節(jié)點位移），邊坡水平面采用*Sflow邊界和只設置水平面零孔壓邊界均只需43子步完成計算，中間只報1U，收斂效果完全相同，孔壓隨時間動態(tài)演變，直至平衡。（附件slope_sflow2、slpoe_pore）
邊坡孔壓 /Pa（每個Frame0.5小時）
（3）瞬態(tài)滲流考慮固結(jié)沉降時，邊坡水平面和斜坡面均采用*Sflow邊界和設置水平面零孔壓邊界、斜坡面*Sflow邊界均只需43子步完成計算，中間只報1U，收斂效果完全相同，孔壓隨時間動態(tài)演變，直至平衡。（附件slope_sflow12、slpoe_pore_sflow1）
邊坡孔壓 /Pa（每個Frame0.5小時）
4. 隧洞算例（小三維C3D8P）
隧洞尺寸
（1）瞬態(tài)滲流不考慮固結(jié)沉降時，洞壁采用*Sflow邊界和采用零孔壓邊界收斂效果完全相同，但孔壓在第一子步就達到穩(wěn)狀態(tài)定，沒有隨時間的變化過程。（附件tunnel_allfix）
第一子步孔壓計算結(jié)果 /Pa
（2）瞬態(tài)滲流考慮固結(jié)沉降時，洞壁采用*Sflow邊界和采用零孔壓邊界收斂效果完全相同，而且孔壓均隨時間動態(tài)演變，逐漸穩(wěn)定。（附件tunnel_pore0、tunnel_sflow）
隧洞孔壓 /Pa（每個Frame0.1小時）
隧洞位移 /m（每個Frame0.1小時）
展開 