吳子牛 王莎 盧欣奇

摘要：ABAQUS軟件界面操作復(fù)雜，且前后數(shù)據(jù)處理需花費(fèi)大量時(shí)間，針對以上問題，采用C#語言對ABAQUS進(jìn)行二次開發(fā)，突破傳統(tǒng)的單純使用Python腳本配合ABAQUS內(nèi)置GUI插件工具箱開發(fā)的局限性，利用C#語言編寫程序，直接調(diào)用ABAQUS求解器，實(shí)現(xiàn)在不打開軟件界面的情況下進(jìn)行計(jì)算并讀取結(jié)果，并將強(qiáng)度折減系數(shù)顯示在WinForm窗體程序上。開發(fā)后的系統(tǒng)大大簡化了ABAQUS中求解特定巖土邊坡工程問題的復(fù)雜操作，降低了軟件操作難度。

關(guān)鍵詞：ABAQUS;二次開發(fā);C#語言;邊坡穩(wěn)定性;強(qiáng)度折減法

中圖分類號(hào)：TD76

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1001-1277（2022）02-0061-04

doi：10.11792/hj20220211

引 言

邊坡是巖土工程的重要組成部分，對邊坡穩(wěn)定性分析的研究歷來是巖土工程界的重點(diǎn)和難點(diǎn)課題。隨著數(shù)值模擬技術(shù)的發(fā)展，基于有限單元法的強(qiáng)度折減技術(shù)在邊坡穩(wěn)定性分析中得到了廣泛的應(yīng)用[1]。ABAQUS作為大型通用有限元分析軟件的一種，在求解非線性問題方面十分優(yōu)秀，在模擬巖土工程高度非線性問題中，具有優(yōu)秀的分析能力和模擬復(fù)雜系統(tǒng)的可靠性。因此，工程界常將ABAQUS用于巖土工程的邊坡問題分析[2]。

雖然ABAQUS擁有強(qiáng)大的計(jì)算分析功能，但其在實(shí)際應(yīng)用過程中卻對操作人員的經(jīng)驗(yàn)技巧有一定的要求。首先，ABAQUS操作建模界面非常復(fù)雜，操作十分繁瑣，前處理需要花費(fèi)大量時(shí)間;其次，由于其計(jì)算結(jié)果種類眾多，顯示類型豐富，提取所需結(jié)果需要大量后處理操作?；谝陨蟽牲c(diǎn)，導(dǎo)致大部分初學(xué)者無法在短時(shí)間內(nèi)熟練操作軟件，同時(shí)在實(shí)際工程應(yīng)用中也需要花費(fèi)大量的時(shí)間用于各類參數(shù)調(diào)整。

對于上述在計(jì)算前后數(shù)據(jù)處理時(shí)遇到的問題，本文探討利用C#語言和WinForm窗體程序，配合ABAQUS優(yōu)秀的二次開發(fā)特性，編寫開發(fā)友好的圖形用戶界面，簡化利用ABAQUS進(jìn)行邊坡穩(wěn)定性分析計(jì)算強(qiáng)度折減系數(shù)的操作，提升軟件易用性和便利性。

1 傳統(tǒng)ABAQUS二次開發(fā)簡介

ABAQUS僅提供了基于Python腳本的二次開發(fā)接口，其在繼承Python自身函數(shù)庫的基礎(chǔ)上，又對函數(shù)庫進(jìn)行了進(jìn)一步的擴(kuò)展，使用戶通過操作調(diào)用這些函數(shù)，可以跳過ABAQUS/CEA界面直接操作其運(yùn)算內(nèi)核，實(shí)現(xiàn)快速計(jì)算處理。由于選擇了Python語言作為操作腳本，ABAQUS軟件可以通過以下多種方式與自身內(nèi)核進(jìn)行交互：圖形用戶界面（GUI）、界面窗體下部命令行接口（CLI）、Python腳本（Script）[3-5]，命令流和ABAQUS/CEA內(nèi)核關(guān)系如圖1所示。

使用ABAQUS自帶的RSG對話框構(gòu)造器（RSG dialog Builder）創(chuàng)建ABAQUS GUI插件程序有其自身的局限性。首先，插件程序只是一個(gè)依附于ABAQUS/CEA界面的懸掛式工具，其使用仍然需要打開主軟件界面，同時(shí)由于ABAQUS軟件占用內(nèi)存較大，打開界面需要占用大量時(shí)間和電腦內(nèi)存，會(huì)降低提交計(jì)算分析的速度;其次，雖然ABAQUS可以編寫GUI插件程序簡化建模操作，但是對于初學(xué)者或是缺少相關(guān)領(lǐng)域?qū)I(yè)知識(shí)的非專業(yè)從業(yè)者來說，計(jì)算完成后的后處理及找到并提取需要的計(jì)算結(jié)果（如安全系數(shù)等）又是一大難點(diǎn)。

2 基于C#語言的ABAQUS二次開發(fā)

2.1 開發(fā)技術(shù)路線

C#是微軟公司開發(fā)的一種面向?qū)ο蟮木幊陶Z言，用其編寫的WinForm窗體程序運(yùn)行在Windows操作系統(tǒng)下，具有速度快、效率高等特點(diǎn)。結(jié)合傳統(tǒng)開發(fā)思路存在的問題和局限性，運(yùn)用C#語言編寫WinForm窗體程序，跳過ABAQUS/CEA界面直接將建模參數(shù)傳送給ABAQUS求解器，從而實(shí)現(xiàn)只需要打開WinForm窗體程序，輸入用于建立邊坡模型的幾何和物理力學(xué)參數(shù)，通過計(jì)算就可以直接在窗體程序上顯示計(jì)算好的強(qiáng)度折減系數(shù)，達(dá)到不打開ABAQUS軟件也可操作其進(jìn)行模擬計(jì)算的目的，使本行業(yè)外的非專業(yè)人士也能操作專業(yè)模擬軟件。

由于ABAQUS只提供了Python語言的二次開發(fā)接口，沒有提供能夠利用C#語言在.NET平臺(tái)下開發(fā)設(shè)計(jì)的通用程序接口和函數(shù)。因此，核心問題就是解決使用C#語言編寫WinForm窗體程序和實(shí)現(xiàn)ABAQUS軟件數(shù)據(jù)通信。

在Windows操作系統(tǒng)下，打開ABAQUS軟件可以使用系統(tǒng)自帶的命令提示符cmd.exe，在Windows命令提示符的控制臺(tái)中輸入命令：Abaqus cae noGUI =*.py，即可直接命令A(yù)BAQUS軟件運(yùn)行指定的Python腳本，同時(shí)不需要打開ABAQUS/CEA界面。因此，可以將Python腳本作為數(shù)據(jù)溝通的中間橋梁，使用C#語言編寫的WinForm窗體程序修改預(yù)設(shè)好的Python腳本模板，再通過調(diào)用計(jì)算機(jī)進(jìn)程，啟動(dòng)命令控制符，控制ABAQUS軟件運(yùn)行修改好的指定腳本代碼，計(jì)算結(jié)束后再將計(jì)算結(jié)果打印輸出到WinForm窗體程序上，軟件運(yùn)行流程如圖2所示。

2.2 Python腳本編寫思路

將Python腳本文件作為數(shù)據(jù)交流的中間媒介，利用WinForm窗體程序把該腳本文件當(dāng)作文本文件，以字符串形式進(jìn)行讀寫修改。為了方便修改，可以將Python腳本文件中的代碼寫成一個(gè)函數(shù)，預(yù)留需要輸入的參數(shù)，后續(xù)可以直接在腳本文件最后添加運(yùn)行特定參數(shù)的函數(shù)代碼，令A(yù)BAQUS讀取運(yùn)行，便可計(jì)算得出結(jié)果，解決C#語言無法直接操作ABAQUS的問題。

所有Python代碼都寫在名為CalSafety_factor的函數(shù)中，函數(shù)參數(shù)包括：地基寬、地基高、邊坡下底寬、邊坡上底寬、邊坡高度、邊坡土體容重、內(nèi)聚力、內(nèi)摩擦角、彈性模量、泊松比。函數(shù)代碼主要包括7個(gè)部分：創(chuàng)建劃分模型，設(shè)置材料屬性，裝配，設(shè)置分析步，設(shè)置載荷邊界條件，網(wǎng)格劃分（內(nèi)置默認(rèn)：網(wǎng)格單元形狀（Element sharp）為四邊形（Quad），劃分技術(shù)（Technique）為掃略（sweep），網(wǎng)格寬度為1單位），任務(wù)和提交（內(nèi)置默認(rèn)：強(qiáng)度折減系數(shù)K取0.5～2）。

2.3 軟件編程思路和代碼

通過C#語言編寫的WinForm窗體程序，提供給使用者界面交互。通過預(yù)留文本框讓用戶輸入建模所需的幾何和物理力學(xué)參數(shù)。設(shè)置按鈕控件，通過按鈕的Button_Click點(diǎn)擊事件，觸發(fā)計(jì)算，最終將計(jì)算結(jié)果（強(qiáng)度折減系數(shù)）用label標(biāo)簽控件顯示在此界面上。

編程思路如下：

1）創(chuàng)建以Result命名的文件夾，在其中創(chuàng)建不同的以時(shí)間命名的文件夾，用于存放修改好的Python腳本和計(jì)算結(jié)果。通過用戶輸入數(shù)據(jù)給Python腳本模板中添加包含指定參數(shù)的代碼觸發(fā)函數(shù)運(yùn)行，同時(shí)添加輸出結(jié)果的代碼，將計(jì)算結(jié)果（強(qiáng)度折減系數(shù)）打印在Python腳本文件最后的注釋中，最后將腳本存放在特定文件夾中。

2）通過創(chuàng)建Process類型對象啟動(dòng)Windows命令控制臺(tái)程序。通過cmd.exe命令控制符找到指定文件夾中改寫好的文件，并通過Abaqus cae noGUI = *.py命令控制ABAQUS軟件運(yùn)行指定腳本，進(jìn)行計(jì)算。

3）計(jì)算完成后，讀取打印在Python腳本文件注釋中的結(jié)果（強(qiáng)度折減系數(shù)），將其顯示在WinForm窗體程序的界面上。

3 邊坡穩(wěn)定性算例分析

3.1 強(qiáng)度折減有限元法基本原理

強(qiáng)度折減有限元法是由ZIENKIEWICZ等[6]提出的一種邊坡穩(wěn)定性分析方法，其基本原理是將邊坡強(qiáng)度參數(shù)中的初始內(nèi)聚力（C）和內(nèi)摩擦角（φ）同時(shí)除以一個(gè)強(qiáng)度折減系數(shù)（K），得到一組新的強(qiáng)度參數(shù)值;然后以這組新參數(shù)值作為輸入?yún)?shù)進(jìn)行有限元計(jì)算。通過不斷改變強(qiáng)度折減系數(shù)，反復(fù)代入，循環(huán)計(jì)算[7]，當(dāng)計(jì)算至邊坡達(dá)到極限平衡狀態(tài)時(shí)，對應(yīng)的強(qiáng)度折減系數(shù)值就可以認(rèn)為是邊坡的最小安全系數(shù)，發(fā)生破壞的滑動(dòng)面就是最危險(xiǎn)的滑動(dòng)面。鄭宏等[8]討論了邊坡穩(wěn)定有限元分析中安全系數(shù)的定義問題，指出了強(qiáng)度折減有限元法的合理性。

折減后的邊坡強(qiáng)度參數(shù)中的初始內(nèi)聚力（C′）和內(nèi)摩擦角（φ′）可分別用式（1）、式（2）求得[6]。

3.2 可靠性驗(yàn)證

此處選擇的算例是DAWSON等[9]（1999）分析的一個(gè)均質(zhì)邊坡。確定以下參數(shù)：地基寬（Foundation_Length），地基高（Foundation_Height），邊坡上底寬（slope_topLen），邊坡下底寬（slope_bottomLen），邊坡高度（slope_height），邊坡土體容重（Volumetric_weight），內(nèi)聚力（Cohesion），內(nèi)摩擦角，彈性模量（Elastic），泊松比（Poisson）。具體參數(shù)值如表1所示，邊坡幾何要素如圖3所示，其兩邊和底部有位移邊界。

通過啟動(dòng)WinForm窗體程序，輸入指定參數(shù)，點(diǎn)擊計(jì)算，計(jì)算得出的強(qiáng)度折減系數(shù)顯示在軟件界面上，如圖4所示。

通過ABAQUS軟件自帶的云圖查看功能，檢驗(yàn)計(jì)算結(jié)果的可靠性，計(jì)算最大塑性應(yīng)變分布云圖如圖5所示，場變量Fv1強(qiáng)度折減系數(shù)云圖如圖6所示。通過ABAQUS打開計(jì)算結(jié)果文件，查詢得知此算例在運(yùn)行到第二個(gè)分析步的0.371 831 s時(shí)無法計(jì)算收斂，計(jì)算停止。由于強(qiáng)度折減系數(shù)的不斷迭代增大，導(dǎo)致邊坡強(qiáng)度參數(shù)中的內(nèi)聚力和內(nèi)摩擦角逐漸減小，最終導(dǎo)致邊坡不穩(wěn)，發(fā)生滑坡。

結(jié)合云圖計(jì)算結(jié)果可知，滑裂面位置和形狀合理，安全系數(shù)計(jì)算結(jié)果近似等于1.0，與Dawson分析結(jié)果一致，因此，二次開發(fā)系統(tǒng)計(jì)算結(jié)果理想可靠。

4 結(jié) 語

本文針對ABAQUS軟件前、后處理復(fù)雜的問題，采用C#語言編寫WinForm窗體程序，通過計(jì)算機(jī)進(jìn)程啟動(dòng)命令提示符達(dá)到對ABAQUS軟件的控制和調(diào)用，使其讀取指定的Python腳本，并實(shí)現(xiàn)計(jì)算結(jié)果直接在WinForm窗體程序界面顯示的功能。通過自編程序進(jìn)行二次開發(fā)，在計(jì)算巖土邊坡問題時(shí)，省略了打開ABAQUS軟件操作CEA界面的復(fù)雜過程，大大降低了有限元軟件的使用難度。

[參 考 文 獻(xiàn)]

[1] 夏園園.基于ABAQUS有限元強(qiáng)度折減法的邊坡穩(wěn)定性分析[J].價(jià)值工程，2018，37（31）：198-202.

[2] 費(fèi)康，彭劼.ABAQUS巖土工程實(shí)例詳解[M].北京：人民郵電出版社，2017.

[3] 賈利勇，富琛陽子，賀高，等.Abaqus GUI程序開發(fā)指南PYTHON語言[M].北京：人民郵電出版社，2016.

[4] 謝瑞敏，王雪軍.基于Python的ABAQUS二次開發(fā)應(yīng)用[J].中國水運(yùn)（下半月），2017，17（10）：77-79.

[5] 曹金鳳，王旭春，孔亮.Python語言在Abaqus中的應(yīng)用[M].北京：機(jī)械工業(yè)出版社，2011.

[6] ZIENKIEWICZ O C，HUMPHESON C，LEWIS R W.Associated and non-associated visco-plasticity and plasticity in soil mechanics[J].Geotechnique，1975，25（4）：671-689.

[7] 閆俊維，王曙光，陳靜瑜.ABAQUS強(qiáng)度折減法在邊坡穩(wěn)定性分析中的應(yīng)用[J].西部交通科技，2015（8）：21-24.

[8] 鄭宏，田斌，劉德富，等.關(guān)于有限元邊坡穩(wěn)定性分析中安全系數(shù)的定義問題[J].巖石力學(xué)與工程學(xué)報(bào)，2005，24（13）：2 225-2 230.

[9] DAWSON E M，ROTH W H，DRESCHER A.Slope staility analysis by strength reduction[J].Geotechnique，1999，49（6）：835-840.

3392500338263