phase2的案例
邊坡剪切強度折減分析(Shear Strength Reduction Analysis)
(1) 項目設置
與Phase2一樣設置合適的單位，另外需要顯式地設置重力加速度。
(2) 網格劃分
邊坡的幾何形狀可以在ADONIS內直接輸入坐標，也可以輸入從Phase2導出的dxf文件，不過應注意的是，目前這個軟件有缺陷，從Phase2導出的dxf文件必須在AutoCAD中重新保存一下才能正確導入。選擇6節點的三角形單元，最大邊長取1.5。
(3) 邊界條件
邊界使用xyfix命令固定三邊，邊坡面不能固定。
(4) 材料性質
使用與Phase2相同的材料參數，如下圖所示。
(5) 計算結果
ADONIS使用了類似FLAC的Bracket計算方法，這種方法比Phase2的遞減方法效果好。計算出的安全系數是1.14，這個數值與SLIDE的計算結果完全相同。
位移等值線圖不完全相同，下圖所示左圖是Phase2的計算結果，最大位移量5m, 右圖是ADONIS的計算結果，最大位移量是3m左右。這是由于兩個程序塑性變形后的計算方法不一樣導致的。從工程角度來看，這個數值沒有什么實際意義，因為在這種狀態下，邊坡已經發生了破壞。
展開 免費二維有限元分析軟件ADONIS計算步驟
material("create","IsoElastic","matid",1,"matname","Example","density",2500,"shear",2.8e+9,"bulk" 3.9e+9)material("assign","matid",1)
5 網格生成
ADONIS的網格目前只有兩種：三節點的三角形和六節點的三角形單元，Phase2除此之外，還提供了四節點和八節點的四邊形單元。
discretize("maxedge",1.5)triangle("maxedge",1.5)
6 邊界條件
ADONIS提供了多種邊界條件，對于地下開挖，邊界的x,y固定即可。
applybc("xyfix","xlim",-34.384,42.255,"ylim",-32.520,32.520)
7 原巖應力
巖石工程地下開挖模擬，必須提供原巖應力。
initial("sxx",-30e6)initial("syy",-30e6)initial("szz",-30e6)
8 計算設置
計算設置在很大程度上影響著有限元分析的結果。下圖比較了ADONIS與Phase2的計算設置。
solve()
9 結果顯示
ADONIS提供了位移，應力，應變，孔隙壓力的等值線圖。
10 與Phase2比較
ADONIS計算的最大位移量是2.48cm;
Phase2計算的最大位移量是2.09cm.
展開 土工格柵(Geogrid)加固軟土地基的堤壩---兩種模擬方法比較
為了進一步檢驗Geogrid所起的作用，下面分析了一個堤壩的加固，分別使用Plaxis 2D和Phase2進行了模擬，同時也比較了二者在模型設置以及計算結果方面的差異。
2 堤壩模型
這個堤壩模型由兩層土組成，上層為Sand Fill, 下層為Soft Clay，顯然上層的強度比下層高，為了改善土的性能，在兩層的分界面使用Geogrid進行加固。首先分析沒有加固時的穩定性，然后分析加固后的穩定性。
3 Plaxis 2D模擬
當輸入土層材料參數時，為了最大程度地與Phase 2的計算作比較，在"Initial" ko設置時，采用了"Manual" 選項，設置ko為1，即土體處于靜水壓力狀態。下面分三種情形進行分析：
(1) 在沒有支護的情況下，最大位移量是0.088m。下圖所示的是位移云圖和最大剪應變圖。計算的安全系數為1.244。
(2) 安裝Geogrid但不設置界面，這相當于Geogrid與土緊密粘合在一起，不發生滑動或分離。在這種情況下，最大位移量0.012m，計算的安全系數是1.265。可以看出，安裝Geogrid有效地阻止了土層的位移，安全系數得到了提高。剪應變圖顯示出剪切帶被Geogrid分割開，沒有形成貫通的剪切帶。
(3) 安裝Geogrid同時設置界面，即考慮了土-結構的相互作用，在這種情況下，最大位移量0.019m，計算的安全系數是1.268，其結果與不設置界面時差不多。
4 Phase 2
使用相同的模型和參數，在Phase 2下運行。分兩種情形：
(1) 不進行支護。
展開 雙層土邊坡穩定性分析出現的問題
使用Phase2計算出的安全系數是1.43，最大剪應變圖如下所示(右圖)，Phase2計算出來的位移非常合理。接下來將試驗<Phi-C reduction and comparison with Bishop's method>中例子來找出發生位移錯誤的具體原因，感覺問題出現在分階段施工的設置步驟上。
Matlab計算相位差
算例：振動頻率f=10Hz；采樣頻率Fs=1000；相位差phi=π/2；數據長度N=1024；
2.1 頻域法
步驟：FFT→求能量最大點的相位（phase）→計算相位差
2.1.1 代碼
% 測試相位差——頻域法
% 2019-2-18
% By zhipeng feng
% 算例：振動頻率f=10Hz；采樣頻率Fs=1000；相位差phi=π/2；數據長度N=1024；
clc; clear; close all
f=10;
Fs=1000;
phi=pi/2;
N=1024*2;
f1=10;
f2=20;
t=(0:N-1) / Fs; %首先生成時間序列
y1=sin(2*pi*f*t);
y2=sin(2*pi*f*t+phi);
figure(1)
plot(t,y1,'r',t,y2,'b');
legend('信號1','信號2');
grid
%FFT
frequency=Fs*(0:N/2)/N; %單邊功率譜
fft_y1=fft(y1);
fft_y2=fft (y2);
%求能量最大點
[~,idx1]=max(abs(fft_y1));
[~,idx2]=max(abs(fft_y2));
% 求兩信號能量最大點的相位
phase1=phase(fft_y1(idx1))*180/pi
phase2=phase(fft_y2(idx2))*180/pi
fprintf('信號2比信號1相位差=%f度\n',phase2-phase1);
運行結果：
信號2比信號1相位差=88.651114度
其實，對于同頻率的兩個向量，其相位滯后dt
展開 軟土地層開挖和支護模擬(Excavation and Support of Soft Soil)---Part 2
(2) 增加Phase_1。Phase 1激活線性載荷，板單元和界面元,這個階段用來模擬開挖前的準備工作。盡管可以使用"Select multiple objects"工具逐個激活這些分量，但最快捷的方法是在模型瀏覽器中選擇。
(3) 增加Phase_2。Phase 2模擬第一步開挖2m。停用這一層土即表示開挖。
(4) 增加Phase_3。Phase 3模擬安裝支桿。激活支桿即可。
(5) 增加Phase_4。Phase 4模擬第二步開挖8m。停用這一層土即表示開挖。在Plaxis 2D中，當停用一個土層時，孔隙壓力不會自動被停用。因此在這種情況下，水仍然留在開挖區域內，模擬的是一個被水淹沒的開挖。
(6) 增加Phase_5。Phase 5模擬第三步開挖。停用這一層土即表示開挖。
至此完成了所有的計算步驟。
展開 FLUENT氣泡破碎與凝聚模擬
（2）單擊任務欄中 （云圖）按鈕，彈出Insert Contour（創建云圖）對話框，單擊OK按鈕進入云圖設定面板。
（3）在Geometry（幾何）選項卡中Locations選擇symmetry 1，Variable選擇Phase 2.Mean Particle Diameter，單擊Apply按鈕創建腐蝕速率云圖。
FLUENT段塞流模擬
（2）單擊主菜單中Setting Up Physics→Materials→Create/Edit，彈出Create/Edit Materials（材料）對話框。單擊Fluent Database按鈕彈出Fluent Database Materials對話框，選擇water liquid單擊Copy按鈕確認。
5 定義模型
（1）單擊命令結構樹中General按鈕，彈出General（總體模型設定）面板。在SolverTime中選擇Transient。勾選Gravity，在Z中填入-9.81m/s2。
（2）在模型設定面板Models中雙擊Multiphase按鈕，彈出Multiphase Model（多相流模型）對話框，選擇VOF，勾選Implicit Body Force，單擊OK按鈕確認并關閉對話框。
（3）在模型設定面板Models中雙擊Multiphase下的Phases按鈕，彈出Phase（多相流設置）對話框，Phase-1保持默認。
在Phase-2對話框中，Phase Material選擇water-liquid，單擊OK按鈕確認并關閉對話框。
（4）在模型設定面板Models中雙擊Viscous按鈕，彈出Viscous Model（湍流模型）對話框，在Model中選擇SST k-omega(2eqn)，單擊OK按鈕確認并關閉對話框。
6 設置邊界條件
（1）單擊主菜單中Setting Up Physics→Zones→Boundaries按鈕啟動的邊界條件面板。
（2）在邊界條件面板中，雙擊gas_inlet彈出邊界條件設置對話框。Velocity Magnitude輸入0.6，單擊OK按鈕確認退出。
展開 操作教程 | FLUENT直列式油水分離器模擬
（2）單擊主菜單中File→Read→Mesh命令，導入.msh網格文件。
2. 定義模型
（1）單擊命令結構樹中General按鈕，彈出General（總體模型設定）面板，Solver中Time選擇Transient。勾選Gravity，在Z中填入-9.81m/s2。
3. 設置材料
（1）單擊主菜單中Setting Up Physics→Materials→Create/Edit，彈出Create/Edit Materials（材料）對話框。單擊Fluent Database按鈕彈出Fluent Database Materials對話框，選擇water-liquid單擊Copy按鈕確認，修改其密度及粘度。
（2）創建新物質oil-liquid。
4. 設置多相流模型
（1）在模型設定面板Models中雙擊Multiphase按鈕，彈出Multiphase Model（多相流模型）對話框，選擇Eulerian，Number of Eulerian Phases填入2，單擊OK按鈕確認。
（2）在Phase選項卡中，Phase-1的Phase Material選擇water，Phase-2的Phase Material選擇oil-liquid，單擊OK按鈕確認。
（3）在Phase Interaction選項卡中，激活表面張力模型。
5. 設置邊界條件
（1）單擊主菜單中Physics→Zones→Boundaries按鈕啟動的邊界條件面板。
（2）在邊界條件面板中，設置inlet的參數如下圖所示。
展開 Drift-flux 模型（2）
如果要進行改動，需要建立在如下的公式基礎上
Volume Fraction of Phase #2 at Inversion Point
Volume fraction of phase #2 at inversion point 指的是 Phase #2 在網格單元中的最大體積率。
對于空氣，常假設氣泡是不連續的，取 1 較為合適。
而對于液體，一般很難確定具體在網格單元中的最大體積率，常取 0.5 進行初步計算。
GLAD中二元光學元件建模
在二元光學中，每個主階次上的次階次數目通常設置為2、4、8等。GLAD中產生二元光學元件命令如下所示：
binary/lens/surface kbeam xrad yrad level nlevels
int2phase/two kbeam1 kbeam2 [2.*pi*(rindex-1)/lambda]
binary/lens/phasescreen ibeams rindex xrad yrad level nlevels
binary/lens/residual ibeams rindex xrad yrad level nlevels
binary/surface kbeam level nlevels
這些命令執行的是產生二元光學的光柵和透鏡，其二元光學表面可以由binary/surface 命令產生，并直接或者間接依靠 int2phase、int2waves、sfocus起著相位屏的作用。二元光學表面可以圖示為plot的強度項。
binary/surface 命令能夠將任意分布的光場轉化為二元光學器件的面形。
二元光柵表面計算：
這里舉一個二元光柵的例子，它是由二元表面組成，然后對其執行“sfocus”命令。表面通過具有一定高度的強度表現出來，高度的單位為厘米，表面的高度可以通過任意一個常規的強度出圖命令顯示出來。“sfocus”命令將表面的高度作為一個相位差引入，然后乘以α = n – 1，傳輸距離為z，將其設定為1。“sfocus”命令還包括了2π ? λ這一項。次級數量分別為2,4,8和16，在這樣的單位下，16個次級分辨率不好，但是表面卻可以以一個近似光滑的分布顯示出來。
展開 
Longformer---回答問題長文本的Transformer
回答: finite element and the finite difference methods
(2) 問題: What is the most popular codes for rock slope analysis?
回答: Phase2 and RS3
(3) 問題: Where is Chuquicamata mine?
回答: major joints sets, shear planes and fault planes
(4) 問題: What is fracture propagation paths?
回答: ability to consider intact rock deformation and movements, complex hydro-mechanical and dynamic analyses, and additional insight into identification of potential instability mechanisms
(5) 問題: How to simulate step-path brittle fracture?
回答: using a simplified discrete fracture network (DFN) model coupled with a Voronoi polygonal mesh approach
(6) 問題: What is the toppling failure?
回答: rotation of columns or blocks of rocks about some fixed base
展開 HSCAE 3D 7.0 smart 華塑注塑成形流動分析系統\
SDCAD 4.0_無解密不好用.rar
Bentley AutoPIPE V8i v09.03.02.05\
Bentley AutoPIPE V8i 09.01.01.04 CHS\
CADMeister.v5.1-ISO 1DVD\
DIgSILENT PowerFactory v14.0_jusb\
dotNET20\
GMI Caliper v3.1\
GMI Imager v5.6\
GMI ModelBuilder v1.4\
GMI Mohrfracs v2.6\
GMI PressCheck v2.5\
GMI SFIB v5.4\
GMI WellCheck v2.5\
pvcad 3.5\
PVElite.v2011-ISO\
SFTC_DEFORM_V10\
Synplify.Premier.pro v9.61 Linux\
212\
Intercorr.Predict.v4.0.rar
INTERCORR.PREDICTPIPE.V3.0.rar
FF-CAM v6.0-ISO_FFCAM_FF\
FracMan v7.00 1CD(裂縫建模軟件)\
Geomagic.Qualify.v12.12\
Intercorr.Predict.v4.0\
INTERCORR.PREDICTPIPE.V3.0\
Lift.Designer.v5.2.Premium.Suite.MultiLanguage-ISO 1CD(電梯系統設計)\
MASTERCAM.X.V10.2.MR2\
Rocscience.Phase2.v8005.Win32\
展開 GLAD中二元光學元件建模
在二元光學中，每個主階次上的次階次數目通常設置為2、4、8等。GLAD中產生二元光學元件命令如下所示：
binary/lens/surface kbeam xrad yrad level nlevels
int2phase/two kbeam1 kbeam2 [2.*pi*(rindex-1)/lambda]
binary/lens/phasescreen ibeams rindex xrad yrad level nlevels
binary/lens/residual ibeams rindex xrad yrad level nlevels
binary/surface kbeam level nlevels
這些命令執行的是產生二元光學的光柵和透鏡，其二元光學表面可以由binary/surface 命令產生，并直接或者間接依靠 int2phase、int2waves、sfocus起著相位屏的作用。二元光學表面可以圖示為plot的強度項。
binary/surface 命令能夠將任意分布的光場轉化為二元光學器件的面形。
二元光柵表面計算：
這里舉一個二元光柵的例子，它是由二元表面組成，然后對其執行“sfocus”命令。表面通過具有一定高度的強度表現出來，高度的單位為厘米，表面的高度可以通過任意一個常規的強度出圖命令顯示出來。“sfocus”命令將表面的高度作為一個相位差引入，然后乘以α = n – 1，傳輸距離為z，將其設定為1。“sfocus”命令還包括了2π ? λ這一項。次級數量分別為2,4,8和16，在這樣的單位下，16個次級分辨率不好，但是表面卻可以以一個近似光滑的分布顯示出來。
展開 使用Extrusion工具產生非結構化的網格(unstructured Mesh)
2. 導入dxf文件
如前所述，理論上使用FLAC3D內置的Point-edge工具和Control point工具可以構造出任何形狀的二維圖形[extrude set select "tunnel"]，不過，從外部導入dxf文件[geometry import "tunnel.dxf" format dxf set "tunnel"]是一種比較快捷的方法。在本例中，從Phase2 [邊坡剪切強度折減分析(Shear Strength Reduction Analysis)]中導入一個隧道模型的dxf文件，這個模型外部邊界尺寸大約為32m, 隧道尺寸大約為3.5m, 采用3節點三角形單元劃分模型，外部邊界的單元邊長約為1.8m，使用梯度方法劃分單元，即隧道周圍網格劃分密一些，遠離隧道邊界的區域網格劃分稀疏一些。
3. 網格劃分
導入dxf文件后，第一件要做的事情是改變默認的單元長度(Set default edge paramaters )。每條邊的單元數目最小值為1，最大值為10000；如果設置的單元長度太小，導致一條邊的單元數目超過10000，那么使用最大值10000作為該條邊的單元數目。根據dxf文件的外部邊界長度32m，推算出每個三角形單元的長度大約為1.8m, 近似取每個單元的邊長為2m[extrude edge length-default 2], 邊界劃分如下圖所示[extrude edge create auto]。
(1) 由于僅創建非結構化的網格，因此選擇Create unstructured meshes only[extrude mesh type unstructured]。(2) Target zone size 該參數指定在塊體內所需的單元尺寸(目標尺寸)，默認值為0。
展開 