第5章ANSYS水利工程應(yīng)用實例分析

1、第5章ANSYS水利工程應(yīng)用實例分析 水利工程 / ANSYS重力壩抗震性能分析步驟 ANSYS重力壩抗震性能用實例分析 J本章典型效果圖 ANSYS10.0 土木工程應(yīng)用實例分析 5.1水利工程概述 雖然我國水利資源非常豐富，但河流在地區(qū)和時間分配上很不均衡，許多地區(qū)在枯水季節(jié) 容易出現(xiàn)干早，而在洪水季節(jié)又往往由于水量過多而形成洪澇災(zāi)害。為了解決這一矛盾，人們 修建了許多水利工程來達到防洪、灌溉、發(fā)電、供水、航運等目的，促進國民經(jīng)濟建設(shè)的發(fā)展。 水利工程中各種建筑物按其在水利樞紐中所起的作用，可以分為以下幾類： (1) 擋水建筑物用以攔截河流，形成水庫，如各種壩和水閘以及抵御洪水所用的堤防等

2、。 (2) 泄水建筑物用以宣泄水庫 或渠道)在洪水期間或其它情況下的多余水量，以保證壩(或 渠道)的安全，如各種溢流壩、溢流道、泄洪隧道和泄洪涵管等。 (3) 輸水建筑物為灌溉、發(fā)電或供水，從水庫(或河道)向庫外(或下游)輸水用的建筑物，如 引水隧道、引水涵管、渠道和渡槽等。 (4) 取水建筑物是輸水建筑物的首部建筑，如為灌溉、發(fā)電、供水而建的進水閘、揚水站等。 (5) 整治建筑物用以調(diào)整水流與河床、河岸的相互作用以及防護水庫、湖泊中的波浪和水流 對岸坡的沖刷，如丁壩、順壩、導(dǎo)流堤、護底和護岸等。 由于破壞后果的災(zāi)難性，大型水利工程建設(shè)的首要目標是安全可靠，其次才是經(jīng)濟合理。 所以說研究大壩等

3、水工建筑物的安全分析、評價和監(jiān)控，是工程技術(shù)人員需要解決的課題，正 確分析大壩性態(tài)已經(jīng)成為當務(wù)之急。 當前對各種水利工程評價主要采用有限元分析方法，借助各種有限元軟件對這些水利工程 建筑物進行安全評價，其中應(yīng)用比較廣泛的是ANSY軟件。目前，ANSY軟件在水利工程中主要 應(yīng)用以下幾個方面： (1) 應(yīng)用各種壩體工程的設(shè)計和施工 利用ANSY軟件，模擬各種壩體施工過程以及壩體在使用階段受到各種載荷(如水位變化對 壩體的壓力、地震荷載等)下結(jié)構(gòu)的安全性能進行評價，模擬壩體的溫度場和應(yīng)力場，借助模 擬結(jié)果修改設(shè)計或?qū)误w采取加固措施。 (2) 應(yīng)用于各種引水隧道、引水涵管等設(shè)計和施工 利用ANSY軟

4、件，模擬這些工程開挖、支護、澆注、回填過程，分析結(jié)構(gòu)在載荷作用下的變 形情況、結(jié)構(gòu)的安全可靠度，以及襯砌支護結(jié)構(gòu)在水壓、溫度發(fā)生變化后產(chǎn)生的變形情況和結(jié) 構(gòu)內(nèi)力，依靠ANSY模擬結(jié)果對結(jié)構(gòu)安全性進行評價。 (3) 應(yīng)用于各種水庫閘門的設(shè)計和施工 水庫閘門在上游水作用下將發(fā)生彎曲、扭轉(zhuǎn)、剪切和拉壓等組合變形，利用ANSY中的 SHELL6單元來模擬閘門，禾U用大型結(jié)構(gòu)有限元分析程序 ANSYS對閘門結(jié)構(gòu)進行三維有限元分 析，根據(jù)分析結(jié)果進行強度校核。 5.2 ANSYS重力壩抗震性能分析步驟 重力壩是一種古老而重要的壩型，主要依靠壩體自身重力來維持壩身的穩(wěn)定。巖基上重力 壩的基本剖面呈三角形，上

5、游面通常是垂直的或者稍傾向下游的三角形斷面。 重力壩的具有很多優(yōu)點： (1) 安全可靠。但剖面尺寸較大，抵抗水的滲漏，洪水漫頂，地震或戰(zhàn)爭破壞的能力都比較 強，因而失事率較低。 (2) 對地形、地質(zhì)條件適應(yīng)性強，壩體作用于地基面上的壓應(yīng)力不高，所以對地質(zhì)條件的要 求也較低，低壩甚至可修建在土基上。 (3) 樞紐泄洪容易解決，便于樞紐布置。 (4) 施工方便，便于機械化施工。 (5) 結(jié)構(gòu)作用明確，應(yīng)力計算和穩(wěn)定計算比較簡單。 鑒于重力壩有如此多優(yōu)點，所以它得到了廣泛應(yīng)用。 但是許多大壩都是建在地震多發(fā)和高烈度地區(qū)，并且壩還要承受重力、水壓力等長期載荷 的作用，為確保工程和人民生命財產(chǎn)在偶發(fā)地震

6、載荷作用下的安全，需對大壩做抗震安全分析。 重力壩抗震性能分析一般分以下五個步驟： 1、創(chuàng)建物理環(huán)境 2、建立模型和劃分網(wǎng)格 3、施加邊界條件和載荷(地震荷載) 4、求解 5、后處理(查看計算結(jié)果) 5.2.1創(chuàng)建物理環(huán)境 在定義壩體抗震性能分析問題的物理環(huán)境時，進入ANSYS前處理器，建立這個壩體抗震 性能分析的數(shù)學仿真模型。按照以下幾個步驟來建立物理環(huán)境： 1、設(shè)置GUT菜單過濾 如果你希望通過 GUI路徑來運行ANSYS，當ANSYS被激活后第一件要做的事情就是選 擇菜單路徑：Main MenuPreferences,執(zhí)行上述命令后，彈出一個如圖5-1所示的對話框出現(xiàn) 后，選擇Struc

7、turalo這樣ANSYS會根據(jù)你所選擇的參數(shù)來對GUI圖形界面進行過濾，選擇 Structural以便在進行壩體抗震性能分析時過濾掉一些不必要的菜單及相應(yīng)圖形界面。 2、定義分析標題(/ TITLE ) 在進行分析前，可以給你所要進行的分析起一個能夠代表所分析內(nèi)容的標題，比如“Dam stability Analysis ”,以便能夠從標題上與其他相似物理幾何模型區(qū)別。用下列方法定義分析標題。 命令：/ TITLE GUI : Utility MenuFileChange Title 3、說明單元類型及其選項（KEYOPT選項） 與ANSYS的其他分析一樣，也要進行相應(yīng)的單元選擇。ANSYS

8、軟件提供了 100種以上的 單元類型，可以用來模擬工程中的各種結(jié)構(gòu)和材料，各種不同的單元組合在一起，成為具體的 物理問題的抽象模型。壩體用PLANE42單元來模擬。 大多數(shù)單元類型都有關(guān)鍵選項 （KEYOPTS），這些選項用以修正單元特性。例如，PLANE42 有如下KEYOPTS : KEYOPT（2）包含或抑制過大位移設(shè)置 KEYOPT（3）平面應(yīng)力、軸對稱、平面應(yīng)變或考慮厚度的平面應(yīng)力設(shè)置 KEYOPT （5）解輸出控制 設(shè)置單元以及其關(guān)鍵選項的方式如下： 命令：ET KEYOPT GUI : Main MenuPreprocessorElement Type Add/Edit/Dele

9、te 圖5-1 GUI圖形界面過濾 4定義單位 結(jié)構(gòu)分析只有時間單位、長度單位和質(zhì)量單位三個基本單位，則所有輸入的數(shù)據(jù)都應(yīng)當是 這三個單位組成的表達方式。如標準國際單位制下，時間是秒（s）,長度是米（m），質(zhì)量是千 克（kg），則導(dǎo)出力的單位是kg? m/s2 （相當于牛頓N），材料的彈性模量單位是 kg/m? s2 （相 當于帕Pa）。 命令：/UNITS 5、定義材料屬性 大多數(shù)單元類型在進行程序分析時都需要指定材料特性，ANSYS程序可方便地定義各種 材料的特性，如結(jié)構(gòu)材料屬性參數(shù)、熱性能參數(shù)、流體性能參數(shù)和電磁性能參數(shù)等。 ANSYS程序可定義的材料特性有以下三種： （1）線性或非線性

10、。 （2）各向同性、正交異性或非彈性。 （3）隨溫度變化或不隨溫度變化。 因為進行壩體抗震性能分析時，ANSYS默認譜分析將忽略材料非線性，因此，壩體抗震 性能分析模型采用彈性模型，只需要定義壩體材料屬性中：容重、彈性模量、泊松比。 命令：MP GUI: Main MenuPreprocessorMaterial PropsMaterial Models 或 Mai n Me nuSolutio nLoad Step OptsOtherCha nge Mat PropsMaterial Models 壩體靜力分析時可考慮材料的非線性，但進行抗震性能分析時，需要將非線性參數(shù)：內(nèi)摩 擦角和內(nèi)聚力刪

11、除。 5.2.2建立模型和劃分網(wǎng)格 創(chuàng)建好物理環(huán)境，就可以建立模型。在進行壩體抗震性能分析時，需要建立模擬壩體的 PLANE82單元。在建立好的模型指定特性（單元類型、選項和材料性質(zhì)等）以后，就可以劃分 有限元網(wǎng)格了。 通過GUI為模型中的各區(qū)賦予特性： 1、選擇 Main Menu PreprocessorMeshing Mesh Attributes Picked Areas 2、點擊模型中要選定的區(qū)域。 3、在對話框中為所選定的區(qū)域說明材料號、實常數(shù)號、單元類型號和單元坐標系號。 通過命令為模型中的各區(qū)賦予特性： ASEL （選擇模型區(qū)域） MAT （說明材料號） TYPE （指定單元類

12、型號） 建立大壩模型時，壩體和地基是賦予不同材料屬性，本文只進行壩體抗震性能分析。 進行大壩3-D模擬分析時，用SOLID65模擬混凝土單元和SOLID45模擬巖石單元的，2-D 模擬分析時，只需用一個 PLANE42單元就可以。 5.2.3施加約束和荷載 在施加邊界條件和荷載時，既可以給實體模型（關(guān)鍵點、線、面）也可以給有限元模型（節(jié) 點和單元）施加邊界條件和荷載。在求解時，ANSYS程序會自動將加到實體模型上的邊界條件 和載荷轉(zhuǎn)遞到有限元模型上。 重力壩抗震性能分析中，主要是給壩體底部施加自由度約束。 命令：D 作用在重力壩上的荷載包含水壓力、冰壓力、泥沙壓力、地震力及壩體自重荷載等。 (

13、1) 自重荷載。由壩體體積和材料的容重算出。 (2) 靜水壓力。作用在壩面上的靜水壓力可根據(jù)靜水力學原理計算，分為水平力及垂直力。 1 2 水平力：Psp Hi(5-1) 1 2 _ nHr(5-2) 2 n代表上游壩面坡度系數(shù)。 2 垂直力：Pcz 式中：為水容重，H1代表上游水深， 同理可求得下游壩面的總靜水壓力的水平向分力及垂直分力。 (3) 揚壓力。重力壩壩體混凝土或漿砌石砌體不是絕對不透水的，它們的表面及內(nèi)部存在 著無數(shù)微小的孔隙，壩基巖石本身孔隙雖然很少，但是也存在著節(jié)理、裂隙。重力壩建成擋水 后，在上下游水位差的長期作用下，上游的水將通過這些孔隙及壩體和壩基的接觸面、壩基的 節(jié)

14、理裂隙等向下游滲透，從而使得壩體內(nèi)和壩底面產(chǎn)生滲透水壓力。 (4) 動水壓力。在溢流面上作用有動水壓力，壩頂曲線和下游面直線段上的動水壓力很小， 可忽略不計。只計算反弧段上的動水壓力。 (5) 冰壓力。在寒冷地區(qū)水庫表面冬季結(jié)成冰蓋，當氣溫回升時，冰蓋發(fā)生膨脹，因而對 擋水建筑物上游面產(chǎn)生冰壓力。 (6) 泥沙壓力。水庫蓄水后，入庫水流挾帶泥沙，逐年淤積在壩前，對壩面產(chǎn)生泥沙壓力。 (7) 浪壓力。浪壓力與風速和水庫吹程有關(guān)，但在荷載中所占比重較小，通常忽略。 (8 )地震荷載。主要是由建筑物質(zhì)量引起的地震慣性力、地震動水壓力和動土壓力。 5.2.4求解 1.靜力求解 首先對重力壩進行靜力求解

15、：在ANSYS程序根據(jù)現(xiàn)有選項的設(shè)置，從數(shù)據(jù)庫獲取模型和 載荷信息并進行計算求解，將結(jié)果數(shù)據(jù)寫入到結(jié)果文件和數(shù)據(jù)庫中。得到壩體在靜力荷載作用 下的位移場與應(yīng)力，了解壩體在設(shè)計條件下的工作形態(tài)，對混凝土重力壩方案的可靠性進行評 價評價大壩 命令：SOLVE GUI: Main MenuSolutionSolveCurrent LS 2動力分析求解 由于地震時的地面運動以水平方向為主，在地震力作用下結(jié)構(gòu)的振動也以水平振動為主， 故本次分析只考慮了水平方向的地震載荷的作用。對重力壩抗震性能計算分析可以采用以下幾 種方法： 1) 擬靜力法 擬靜力法是一種把地震的影響用一種折算的靜載荷來表示，求出這種地

16、震荷載后，按照常 規(guī)的靜力法進行壩體的各項應(yīng)力、位移的抗震分析方法。它是假定地震時與地面加速度相同的 加速度作用在壩體各部位，求出地震時的慣性力，然后根據(jù)慣性力來評價大壩的安全性。 根據(jù)擬靜力分析方法，大壩的水平地震慣性力可簡化為： Qh KhCzFW 式中： Kh 水平向地震系數(shù)，為地面最大水平加速度代表值與重力加速度的比值； Cz 綜合影響系數(shù)，重力壩取為1/4； F 地震慣性力系數(shù)； W產(chǎn)生慣性力的建筑物的總重量。 采用擬靜力法計算重力壩的地震作用效應(yīng)時，水深h處的地震動水壓力的代表值的計算: Pw(h) ah (h) wH o 式中： FW (h)作用在直立迎水壩面水深 h處的動水壓力

17、代表值； ah水平向地震加速度的代表值，地震烈度為8時對應(yīng)的值是0.2g； 地震作用效應(yīng)折減系數(shù)，除另有規(guī)定外，取0.25; (h)水深h處的地震動水壓力分布系數(shù)； w水體質(zhì)量密度的標準值； H 0 水總深度。 與水平面夾角為的傾斜迎水壩面，按公式(5-4)計算的動水壓力代表值乘以折減系數(shù)： c /90 2) 反應(yīng)譜分析法 反應(yīng)譜分析法是以單質(zhì)點彈性體系在實際地震過程中的反應(yīng)為基礎(chǔ)，來進行結(jié)構(gòu)反應(yīng)的分 析，它通過反應(yīng)譜巧妙地將動力問題靜力化，使得復(fù)雜的結(jié)構(gòu)地震反應(yīng)計算變得簡單易行。按 照這一理論，應(yīng)用地震譜曲線，就可以按照實際地面運動來計算建筑物的反應(yīng)。反應(yīng)譜是單點 彈性體系對于實際地面運動的

18、最大反應(yīng)和體系自振周期的函數(shù)關(guān)系。對于復(fù)雜的結(jié)構(gòu)可以簡化 為若干振型的疊加，每個振型又可轉(zhuǎn)化為一個單質(zhì)點來考慮。使用已經(jīng)確定的設(shè)計反應(yīng)譜計算 重力壩在地震作用下的反應(yīng)，就歸結(jié)為尋求壩體的自振特性。 (5-3 ) (5-4 ) (5-5 ) 地震產(chǎn)生的破壞，與受力大小和受頻譜的最大振動的持續(xù)時間的都有關(guān)系。在進行譜分析 計算前，首先要計算大壩的自振特性。模態(tài)分析用于確定結(jié)構(gòu)的振動特性，即結(jié)構(gòu)的固有頻率 和振型，它們是結(jié)構(gòu)承受動態(tài)荷載設(shè)計中的重要參數(shù)，也是更詳細的動力分析的基礎(chǔ)。 模態(tài)分析計算中采用了子空間迭代法提取模態(tài)。水深h處的地震動水壓力的作用按公式 (5-6) 轉(zhuǎn)化為相應(yīng)的壩面附加質(zhì)量。

19、7； Pw(h) 8 ah w.H0h( 5-6 ) 根據(jù)如圖5-2所示的大壩設(shè)計的反應(yīng)譜曲線圖，可得大壩反應(yīng)譜曲線方程： 110T( max 1)，0T 0.1 0.1 T Tg (5-7) max， 圖5-2大壩設(shè)計反應(yīng)譜 本次重力壩抗震性能分析中， max取值為2，Tg取值為0.2,特征周期T。取值為0.2S。 3) 時程分析法 時程分析方法是將地震動記錄或人工波作用在結(jié)構(gòu)上.直接對結(jié)構(gòu)運動方程進行積分，求得 結(jié)構(gòu)任意時刻地震反應(yīng)的分析方法，所以動態(tài)時程分析方法也稱為直接積分法。 本次大壩抗震性能實例分析采用反應(yīng)譜分析方法。 5.2.5后處理 后處理的目的是以圖和表的形式描述計算結(jié)果。對

20、于大壩抗震性能分析中，進入后處理器 后，查看大壩變形圖和節(jié)點的位移和應(yīng)力。通過研究大壩的變形、位移和應(yīng)力情況，來綜合判 斷大壩的抗震性能及安全性能。 命令：/ P0ST1 GUI: Main MenuGeneral Postproc 首先查看大壩靜力分析求解結(jié)果，再查看大壩動力分析求解結(jié)果。 5.3 ANSYS重力壩抗震性能實例分析 5.3.1實例介紹 圖5-3重力壩斷面結(jié)構(gòu) 實例選取應(yīng)用非常廣泛的重力壩，斷面結(jié)構(gòu)如圖5-3所示。壩高120米，壩底寬為76米，壩 頂為10米，上游壩面坡度和下游壩面坡度如圖中所示。 因為重力壩結(jié)構(gòu)比較簡單，垂直于長度方向的斷面結(jié)構(gòu)受力分布情況也基本相同，并且大

21、壩的縱向長度遠大于其橫斷面，因此大壩抗震性能分析選用單位斷面進行平面應(yīng)變分析是可行 的。 大壩抗震性能分析的計算條件如下： 1）假設(shè)大壩的基礎(chǔ)是嵌入到基巖中，地基是剛性的。 3 2） 大壩采用的材料參數(shù)為：彈性模量E=35GPa,泊松比 =0.2，容重 =25KN/ m 。 3）計算分析大壩水位為120米。 3 4）水的質(zhì)量密度1000kg/m。 5）大壩設(shè)防地震烈度為8，水平方向地震加速度值為0.2g。 5.3.2 GUI操作方法 5.3.2.1創(chuàng)建物理環(huán)境 1）在【開始】菜單中依次選取【所有程序】/【ANSYS10.0】/【ANSYS Product Launcher】, 得到10.0AN

22、SYS Product Lau ncher ”對話框。 2）選中【File Management】，在 “Working Directory ” 欄輸入工作目錄 “D:ansysexample5-1”, 在Job Name欄輸入文件名Dam ”。 3）單擊“ RUN ”按鈕，進入 ANSYS10.0的GUI操作界面。 4） 過濾圖形界面： Main MenuPreferences,彈出“ Preferences for GUI Filtering”對話框, 選中“ Structural來對后面的分析進行菜單及相應(yīng)的圖形界面過濾。 5）定義工作標題：Utility Menu File Chang

23、e Title，在彈出的對話框中輸入“ Dam seismic Analysis”，單擊“ OK”，如圖 5-4。 田 Ch曰nge Till巳 2d /TITLE Enter new tide am seismiiq Analysis l-lelp Cancel 圖5-4定義工作標題 6）定義單元類型： a. 定義 PLANE42 單元：Ma in Me nu Preprocessor Eleme nt Type Add/Edit/Delete，彈出一 個單元類型對話框，單擊“ Add”按鈕。彈出如圖5-5所示對話框。在該對話框左面滾動欄中選 擇“ Solid ”，在右邊的滾動欄中選擇“ Q

24、uad 4node 42”單擊“ Apply ”，就定義了“ PLANE42 ” 單元。 圖5-5定義PLANE42單元對話框 b. 設(shè)定 PLANE42 單元選項：Main Menu PreprocessorElement TypeAdd/Edit/Delete，彈 出一個單元類型對話框，選中“Type 2 PLANE42 ”單擊“Options”按鈕，彈出一個“PLANE42 element Type options”對話框，如圖 5-6所示。在“ Element behavior K3 欄后面的下拉菜單中 選取“ Plane strain”其它欄后面的下拉菜單采用ANSYS默認設(shè)置就可以

25、，單擊“ OK”按鈕。 xJ | No eitra output Options for PLAJJE42, Elsnsnt Type Ref No. 1 Element ooorl system defined KI Exj a disp lace rrent sTipe$煜 Extra stress oucutK5 OK Cancel Ejtra surface outputK6 圖5-6 PLANE42單元庫類型選項對話框 通過設(shè)置PLANE42單元選項“ K3 ”為“ Plane strain”來設(shè)定本實例分析采取平面應(yīng)變模 型進行分析。因為大壩是縱向很長的實體，故計算模型可以簡化為平

26、面應(yīng)變問題。 7）定義材料屬性 執(zhí)行 Main MenuPreprocessor Material Props Material Models，彈出 “ Define Material Model Behavior ”對話框，如圖 5-7所示。 圖5-7定義材料本構(gòu)模型對話框 在圖5-7中右邊欄中連續(xù)雙擊StructuralLinear Elasticlsotropic”后，又彈出如圖 5-8 所示Lin ear Isotropic Properties for Material Number 1 ” 對話框，在該對話框中EX ”后面的輸 入欄輸入“ 3.5E10”，在“PRXY”后面的輸入欄

27、輸入“ 0.2”，單擊“ OK”。再在選中“ Density” 并雙擊，彈出如圖5-9所示“ Density for Material Number 1 對話框，在DENS”后面的欄中輸 按鈕。 入邊坡土體材料的密度“ 2500”，單擊“ OK ” !3l ineaknlr叩杞 Properties fohlaterla L izar soTp c loafer tai Proper D AcjaJy *13 圖5-36壩體變形圖 ANSA J1.W 4 tdCrt 21： 5-37節(jié)點解云圖繪制對話框 c. 顯示壩體 x 方向位移云圖：Ma in Men uGe neral Postproc

28、Plot ResultsCo ntour PlotNodal Solu，彈出一個“ Con tour Nodal Solution Data ”對話框，如圖5-37所示，用鼠標依次點擊“ Nodal Solution/DOF Solution/ X-Compoment of displacement ”，再單擊OK” 按鈕，就得到壩體 X 方向 位移云圖，如圖5-38所示。此時，壩體水平方向最大位移為 12.111mm，位置發(fā)生在壩頂。 圖5-38 壩體x方向位移云圖 d. 顯示壩體 y 方向位移云圖： Main MenuGeneral PostprocPlot ResultsContour

29、PlotNodal Solu，彈出一個“ Con tour Nodal Solution Data ”對話框，如圖5-37所示，用鼠標依次點擊“ Nodal Solution/Stress/ Y-Compoment of displacement”，再單擊OK” 按鈕，就得到壩體 Y 方向位移云 圖，如圖5-39所示。此時，壩體垂直方向最大位移為 0.44mm，發(fā)生在壩腹處。 5a- -s- TlWii nr 4： -xu心 i:,- j-31- .:tlf I 圖5-39壩體y方向位移云圖 e. 顯示壩體 x 方向應(yīng)力云圖： Ma in Men uGe neral PostprocPlot

30、ResultsCo ntour PlotNodal Solu，彈出一個“ Con tour Nodal Solution Data ”對話框，如圖5-37所示，用鼠標依次點擊“ Nodal Solution/stress/ x-Compoment of stress,再單擊OK”按鈕，就得到壩體 x方向應(yīng)力云圖，女口 圖5-40所示。 f. 顯示壩體 y 方向應(yīng)力云圖：Ma in Men uGe neral PostprocPlot ResultsCo ntour PlotNodal Solu，彈出一個“ Con tour Nodal Solution Data ”對話框，如圖5-37所示，用

31、鼠標依次點擊“ Nodal Solution/stress/ y-Compoment of stress,再單擊OK”按鈕，就得到壩體 x方向應(yīng)力云圖，女口 圖5-41所示。 圖5-40壩體x方向應(yīng)力云圖 圖5-41壩體y方向應(yīng)力云圖 g.顯示壩體第一主應(yīng)力云圖： Mai n Me nuGe neral PostprocPlot ResultsCo ntour PlotNodal 圖5-42第一主應(yīng)力云圖 Solu，彈出一個“ Con tour Nodal Solution Data ”對話框，如圖5-37所示，用鼠標依次點擊“ Nodal Solution/stress/ 1st princ

32、ipal stress,再單擊OK”按鈕，就得到壩體第一主應(yīng)力云圖， 如圖5-42 所示。 2. 抗震性能分析求解結(jié)果 1）繪制壩體振型圖 a. 讀入模態(tài)分析分析求解結(jié)果數(shù)據(jù)：Utility Menu resume，彈出一個Resume Database 對話框，選中剛才保存的文件Dam-modal.db”，單擊OK”按鈕。 b. 繪制壩體第1階震型 讀入第 1 階數(shù)據(jù)：Main Me nuGe neral PostprocRead ResultsFirst set. 繪制壩體第 1 階震型：Main MenuGeneral PostprocPlot ResultsDeformed Shape

33、，彈出一 個Plot Deform Shape對話框，選中Def +undeformed”，單擊OK ”按鈕，得到壩體第 1 階震型，如圖5-43所示。 c. 繪制壩體第2階震型 讀入第 2 階數(shù)據(jù)：Main MenuGeneral PostprocRead ResultsNext set. 繪制壩體第 2 階震型：Main MenuGeneral PostprocPlot ResultsDeformed Shape，彈出一 個Plot Deform Shape對話框，選中Def +undeformed”，單擊OK ”按鈕，得到壩體第 2 階震型，如圖5-44所示。 圖5-43壩體第1階振型

34、依次重復(fù)執(zhí)行c步驟，就可以依次畫出壩體第 圖5-44壩體第2階振型 3到第18階振型，如圖5-45圖5-60所示。 d. 調(diào)出模態(tài)分析結(jié)果的固有頻率和振型：Main Me nu Ge neral Postproc Read Summary,彈 出如圖5-26所示對話框。 ANSYS10.0 土木工程應(yīng)用實例分析 圖5-49壩體第7階振型圖5-50壩體第8階振型 圖5-45壩體第3階振型 圖5-46壩體第4階振型 圖5-47壩體第5階振型 圖5-48壩體第6階振型 第5章ANSYS水利工程應(yīng)用實例分析 圖5-55壩體第13階振型圖5-56壩體第15階振型 AN、 皿4 2 M -b LflX.

35、討律“rtlWT D ul P-Miutaj 圖5-51壩體第9階振型 圖5-52壩體第10階振型 圖5-53壩體第11階振型 圖5-54壩體第12階振型 TK 4 111? inis tJb -13 JTEXL 艱+L4 OSD -E* EE au d 2HDL3IA ANSYS10.0 土木工程應(yīng)用實例分析 圖5-57壩體第15階振型 圖5-58壩體第16階振型 圖5-60壩體第18階振型 Utility Menuresume，彈出一個Resume Dam- comination.db ”，單擊OK 按鈕。 圖5-59壩體第17階振型 2）繪制壩體位移云圖 a讀入抗震性能分析最終求解結(jié)果數(shù)

36、據(jù)： Database對話框，選中剛才保存的文件 b繪制合并模態(tài)求解后壩體位移云圖： 因為響應(yīng)譜分析是在頻域內(nèi)進行的，對于結(jié)構(gòu)動力特性依賴于頻率而變化，因此在模態(tài)分 析后要進行模態(tài)合并求解，才能得到壩體結(jié)構(gòu)真實的總體效應(yīng)。合并模態(tài)求解后，得到壩體在 各階頻率的真實位移云圖。 讀入第 1 階數(shù)據(jù)：Main Me nuGe neral PostprocRead ResultsFirst set. 繪制壩體第 1 階位移云圖：Main MenuGeneral PostprocPlot ResultsContour PlotNodal Solu，彈出一個Con tour Nodal Solution

37、Data 對話框，用鼠標依次點擊 Nodal Solutio n/DOF Solution/ x-Compoment of displacement得到壩體第1階x方向位移云圖，如圖 5-61所示.同法選 擇y-Compoment of displacement,得到壩體第1階y方向位移云圖，如圖 5-62所示。 讀入下一階數(shù)據(jù)：Main Me nuGe neral PostprocRead ResultsNext set. 繪制壩體下一階位移云圖：Main MenuGeneral PostprocPlot ResultsContour PlotNodal Solu，彈出一個Con tour

38、Nodal Solution Data 對話框，用鼠標依次點擊Nodal Solutio n/DOF 第5章ANSYS水利工程應(yīng)用實例分析 Soluti on/ x-Compome ntof displaceme nt得到壩體下一階 x方向位移云圖，如圖 5-63所示.同法選 擇y-Compoment of displacement,得到壩體下一階階 y方向位移云圖，如圖 5-64所示。重復(fù)讀 入下一個數(shù)據(jù)，可繪制出合并模態(tài)求解后壩體真實位移云圖，如圖5-65圖5-72所示。 圖5-65壩體第3階x方向位移云圖圖5-66壩體第3階y方向位移云圖 TW 5IE .ailMW 圖5-61壩體第1階

39、x方向位移云圖 圖5-62壩體第1階y方向位移云圖 圖5-63壩體第2階x方向位移云圖 b-QA MLlLkC 圖5-64壩體第2階y方向位移云圖 MLtHEOV MUI IM mji-u. I AW 4* * SBC *-.118-1 kHrei*4 *u jwliU c AxwljfjAjr 圖5-69壩體第5階x方向位移云圖 圖5-68壩體第4階y方向位移云圖 m百 EK綸 in: iMiuju rW - fflK - El】 L-IL 圖5-70壩體第5階y方向位移云圖 第5章ANSYS水利工程應(yīng)用實例分析 3）繪制壩體應(yīng)力 讀入第1階數(shù)據(jù): /應(yīng)變云圖 Main Menu Ge ne

40、ral PostprocRead ResultsFirst set. 繪制壩體第 1 階應(yīng)力：Main MenuGeneral PostprocPlot ResultsContour PlotNodal Solu , 彈出一個Contour Nodal Solution Data” 對話框，用鼠標依次點擊Nodal Solution/ stress/ x-Compoment of stress得到壩體第1階x方向應(yīng)力云圖，如圖5-73所示。同法選擇y-Compoment of stress/1st principal stress/,得到壩體第1階y方向應(yīng)力/第 1主應(yīng)力云圖，如圖 5-74和

41、圖5-75 所示。 繪制壩體第 1 階主應(yīng)變云圖：Main Me nuGe neral PostprocPlot ResultsCo ntour PlotNodal Solu，用鼠標依次點擊Nodal Solution/ strain/1st principal strain ”得到壩體第1階第1主應(yīng)變云 圖，如圖5-76所示。 讀入下一階數(shù)據(jù)： Main Me nuGe neral PostprocRead ResultsNext set. 繪制壩體下一階應(yīng)力/應(yīng)變云圖：與繪制第1階一樣，就可以繪制出下一階壩體應(yīng)力/應(yīng)變 云圖。重復(fù)讀入下一個數(shù)據(jù)，可繪制出合并模態(tài)求解后壩體應(yīng)力云圖，如圖5-

42、77圖5-7 圖5-75壩體第1階第1主應(yīng)力云圖圖5-76壩體第1階第1主應(yīng)變云圖 圖5-73壩體第1階x方向應(yīng)力云圖 圖5-74壩體第1階y方向應(yīng)力云圖 ANSYS10.0 土木工程應(yīng)用實例分析 圖5-81壩體第3階第1主應(yīng)力云圖圖5-82壩體第3階第1主應(yīng)變云 圖5-77壩體第2階x方向應(yīng)力云圖 圖5-78壩體第2階y方向應(yīng)力云圖 圖5-79壩體第2階第1主應(yīng)力云圖 圖5-80壩體第2階第1主應(yīng)變云圖 第5章ANSYS水利工程應(yīng)用實例分析 圖5-88壩體第4階第1主應(yīng)力云圖圖5-89壩體第5階第1主應(yīng)變云 U u. jflM；科 糾X Jsm 3 呻 M* d-LCr.: KHL MUM

43、卯 AlllTTIJ 圖5-84壩體第3階x方向應(yīng)力云圖 圖5-85壩體第3階y方向應(yīng)力云圖 Ft 昭-4 到I哋 sto-t :W TSeMS -ZM2I5-LMI-tW ”刪門VI卯I -Z2EW 4UM1 3flT 7UE MW 圖5-86壩體第4階x方向應(yīng)力云圖 圖5-87壩體第4階y方向應(yīng)力云圖 ET1 .-4 IFU$“ ；l ! 傅 nr 士丄醬述 Err?tii|. UMJ4 -J-EalEI .loeinT 5rU NLI* ，:HI：4. ANSYS10.0 土木工程應(yīng)用實例分析 圖5-94壩體第7階x方向應(yīng)力云圖圖5-95壩體第7階y方向應(yīng)力云圖 圖5-90壩體第5階第1

44、主應(yīng)力云圖 圖5-91壩體第5階第1主應(yīng)變云 圖5-92壩體第5階x方向應(yīng)力云圖 圖5-93壩體第5階y方向應(yīng)力云圖 ANSYS10.0 土木工程應(yīng)用實例分析 圖5-96壩體第7階第1主應(yīng)力云圖 - 進行模態(tài)求解時并沒有同時進行模態(tài)擴展, 圖5-97壩體第7階第1主應(yīng)變云圖 而是在獲得譜解后又單獨進行了擴展模態(tài)求解。 !指定結(jié)構(gòu)分析 !定義工作標題 定義工作文件名 !定義壩體單元 !定義平面應(yīng)變 !彈性模量 !泊松比 !密度 !創(chuàng)建關(guān)鍵點 這樣，有明顯意義的模態(tài)為1,2, 3, 4, 5和 7階模態(tài)。 - 大壩水平方向最大位移發(fā)生在壩頂。 5.3.2命令流實現(xiàn) (1) 創(chuàng)建物理環(huán)境 /COM,

45、 Structural /TITLE, Dam seismic An alysis /FILNAM,dam,1 (2) 建立模型 !進入前處理器 /PREP7 !定義單元類型 ET,1,PLANE42 KEYOPT,1,3,2 !定義壩體材料屬性 MP,EX,1,3.5E10 MP,PRXY,1,0.2 MP,DENS,1,2500 !創(chuàng)建壩體線模型 K,1, K,2,76, K,3,15.6,104.1, K,4,15.6,120” K,5,5.6,120, 1,1,2 1.2.3 1.3.4 !連接關(guān)鍵點1、2生成直線 1,4,5 1,5,1 !創(chuàng)建邊坡面模型 /PNUM, kp, on

46、!劃分網(wǎng)格生成有限元模型 !打開關(guān)鍵點號開關(guān) lsel,s,1 !選擇線L1 lesize,all,20 !把所選擇線分為20段 lsel,s,2 !選擇線L2 lesize,all,32 !把所選擇線分為32段 lsel,s,3 !選擇線L1 lesize,all,6 !把所選擇線分為6段 lsel,s,4 !選擇線L2 lesize,all,4 !把所選擇線分為4段 lsel,s,5 !選擇線L2 lesize,all,40 !把所選擇線分為40段 amesh,all !劃分壩體單元 save,dam-grid.db !保存網(wǎng)格模型 a,1,2,3,4,5!由關(guān)鍵點生成面 !劃分壩體單元網(wǎng)

47、格 !設(shè)置網(wǎng)格份數(shù) 施加約束和荷載 !壩體底部施加約束 n sel,s,loc,y,0 !選擇Y=0線上所有節(jié)點 d,all,uy !對所選擇節(jié)點約束X、Y方向位移 d,all,ux acel,9.8 sfl,5,pres,0,1101370 (4)求解 1)靜力分析求解 /solu an type,static !施加重力加速度 !施加水壓力荷載 !設(shè)定為靜力求解 nsubst,100 !設(shè)定最大子步數(shù)為100 In srch, on !打開線性搜索 allsel solve !進行求解 save,dam-static,db !保存求解結(jié)果數(shù)據(jù) fin ish 2）抗震性能分析求解 !模態(tài)分

48、析 /solu an type,modal !設(shè)定為模態(tài)求解 modopt,subsp,18 !米用子空間迭代法提取18階模態(tài) solve !進行模態(tài)求解 save,dam-modal,db !保存求解結(jié)果數(shù)據(jù) set,list !調(diào)出模態(tài)分析各階頻率 fin ish !反應(yīng)譜分析 /solu an type,spectrum !設(shè)定為反應(yīng)譜求解 spopt,sprs,18,yes !設(shè)置單點反應(yīng)譜分析選項 svtyp,2 !設(shè)定為地震加速度譜分析 sed,1,0,0 !定義反應(yīng)譜激磁方向是 x方向 freq,3.5138,8.073,11.118,14.3,21.776,24.98,30.7

49、34,34.119,36.263,39.043,40.852,43.237,47.907,4 9.679,52.521,56.661,58.613 !輸入振動頻率 sv,1.456,2,1.899,1.699,1.459,1.4,1.325,1.293,1.276,1.256,1.245,1.231,1.209,1.201,1.191,1.18,1.176 ,1.171 !輸入對應(yīng)反應(yīng)譜值 solve !進行求解 save,Dam-spectrum,db !保存求解結(jié)果 !模態(tài)擴展求解 /solu an type,modal !設(shè)定為模態(tài)求解 expass, on !定義為模態(tài)擴展分析 mxpa nd,18,yes,0.001 !設(shè)置模態(tài)擴