商翔宇　劉崇梁　郝飛　趙光思

摘要：為了彌補(bǔ)當(dāng)前土力學(xué)課程實(shí)踐教學(xué)環(huán)節(jié)的不足，在前期開(kāi)發(fā)的土力學(xué)教學(xué)有限元仿真軟件基礎(chǔ)上，新增了三個(gè)計(jì)算模塊，使得該軟件涵蓋了平板載荷試驗(yàn)、三軸試驗(yàn)、恒定滲流、簡(jiǎn)單邊坡穩(wěn)定性和擋墻土壓力分析等土力學(xué)課程中所有重要問(wèn)題的分析。另外，在每個(gè)仿真模塊中添加了關(guān)于該模塊的土力學(xué)原理、有限元基本理論和實(shí)現(xiàn)方法，以及模型計(jì)算參數(shù)的含義與經(jīng)驗(yàn)取值范圍的介紹，提高了軟件的易用性。優(yōu)化后的軟件涵蓋內(nèi)容全面，具有更好的容錯(cuò)性和易用性，能夠服務(wù)于土力學(xué)課程的本科教學(xué)，達(dá)到以“虛”補(bǔ)“實(shí)”，激發(fā)學(xué)生學(xué)習(xí)興趣和提高學(xué)生學(xué)習(xí)能力的目的。

關(guān)鍵詞：土力學(xué)；有限元仿真；課程教學(xué)

中圖分類號(hào)：G642.0；TU43 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1005-2909（2016）05-0164-05

土力學(xué)是高校土木工程專業(yè)的基礎(chǔ)核心課程之一，也是地質(zhì)工程、工程力學(xué)等本科專業(yè)的必修課程，其具有理論與實(shí)踐并重的顯著特點(diǎn)[1]。然而，一方面由于實(shí)驗(yàn)室條件和教學(xué)資源的制約，當(dāng)下的高校教學(xué)無(wú)法開(kāi)展土力學(xué)課程教學(xué)中諸如現(xiàn)場(chǎng)淺層平板載荷實(shí)驗(yàn)和三軸剪切實(shí)驗(yàn)等重要土力學(xué)實(shí)驗(yàn)，難以達(dá)到預(yù)期理想的教學(xué)效果。另一方面，由于恒定滲流、邊坡穩(wěn)定以及土壓力問(wèn)題等經(jīng)典土力學(xué)問(wèn)題往往是授課重點(diǎn)，但由于這些問(wèn)題涉及到的模型假設(shè)眾多，公式推導(dǎo)和計(jì)算方法復(fù)雜，嚴(yán)重降低了學(xué)生的學(xué)習(xí)興趣。近年來(lái)，將計(jì)算機(jī)數(shù)值仿真技術(shù)引入到教學(xué)中[2-3]，改變了傳統(tǒng)的教學(xué)模式，豐富了傳統(tǒng)的教學(xué)方式，成為當(dāng)前的教學(xué)改革方向。

迄今為止，如Ansys、Abaqus等可用于有限元仿真模擬的大型商用軟件較多。這些軟件，除了需要購(gòu)買知識(shí)產(chǎn)權(quán)費(fèi)用外，還因功能強(qiáng)大，要求使用者具有系統(tǒng)的力學(xué)基礎(chǔ)和有限元知識(shí)才能掌握；另一方面這些軟件對(duì)計(jì)算機(jī)的硬件要求較高，因此并不適用于本科教學(xué)。鑒于此，本文第一作者在美國(guó)Smith和Griffiths教授提供的開(kāi)源簡(jiǎn)易有限元代碼基礎(chǔ)上[4]，基于Visual Fortran的集成編譯環(huán)境[5]，開(kāi)發(fā)了能夠初步滿足人機(jī)互動(dòng)及可視化教學(xué)要求的土力學(xué)課程有限元軟件1.0版本，并取得了一定的令人滿意的教學(xué)效果[6]。

但是隨著該軟件在教學(xué)中的應(yīng)用，發(fā)現(xiàn)其僅能對(duì)平板荷載試驗(yàn)、三軸剪切試驗(yàn)進(jìn)行模擬分析，無(wú)法滿足對(duì)其他諸如滲流、土壓力和邊坡穩(wěn)定等重要經(jīng)典土力學(xué)問(wèn)題的分析。學(xué)生在對(duì)計(jì)算參數(shù)的選擇以及相應(yīng)理論知識(shí)的深入了解方面均存在一定困難。此外，已開(kāi)發(fā)的土力學(xué)課程有限元軟件1.0還存在譬如程序只能向前不能后退、無(wú)法讓用戶在使用過(guò)程中改變參數(shù)等問(wèn)題，以及由于不同計(jì)算模塊合成時(shí)出現(xiàn)的某些數(shù)據(jù)類型錯(cuò)誤和軟件界面亂碼等問(wèn)題，這些都影響了該軟件的應(yīng)用。

針對(duì)上述不足，筆者在前述開(kāi)發(fā)的軟件1.0版本基礎(chǔ)上，增加了土力學(xué)中重要的恒定滲流分析、邊坡穩(wěn)定分析以及土壓力分析三個(gè)計(jì)算模塊，使得新軟件2.0分析模塊涵蓋了土力學(xué)課程本科教學(xué)的所有重要內(nèi)容，并在計(jì)算模塊中增加了對(duì)計(jì)算模型的背景知識(shí)、計(jì)算參數(shù)的物理意義與經(jīng)驗(yàn)取值范圍，以提高程序的易用性。此外，筆者更正了原有程序中的錯(cuò)誤，以提高其容錯(cuò)性，并提供更為完善的人機(jī)交互功能。所開(kāi)發(fā)的優(yōu)化后的土力學(xué)課程有限元仿真軟件2.0能夠更好地進(jìn)行土力學(xué)課程主要實(shí)踐內(nèi)容的仿真模擬，使該課程的教學(xué)內(nèi)容得到了進(jìn)一步豐富，達(dá)到了以“虛”（虛擬仿真）補(bǔ)“實(shí)”（工程實(shí)踐和試驗(yàn)）的預(yù)期教學(xué)效果。

一、軟件的優(yōu)化內(nèi)容

如前所述，本研究所開(kāi)發(fā)的軟件2.0主要在軟件1.0基礎(chǔ)上[6]，新增了三個(gè)計(jì)算模塊，增加了必要的背景知識(shí)和參數(shù)選擇推薦，并修正了部分錯(cuò)誤。限于篇幅，本文主要對(duì)前面兩個(gè)主要優(yōu)化工作進(jìn)行介紹。

優(yōu)化后的程序2.0的啟動(dòng)頁(yè)面如圖1所示，程序主界面中的菜單有平板載荷試驗(yàn)、三軸剪切試驗(yàn)、恒定滲流分析、邊坡穩(wěn)定分析、被動(dòng)土壓力分析、小工具、說(shuō)明等七個(gè)部分。

（一） 提高程序易用性

為了降低本科學(xué)生使用本軟件的知識(shí)門檻，達(dá)到激發(fā)學(xué)生學(xué)習(xí)的興趣，軟件的2.0版本增加各個(gè)計(jì)算模型的分析簡(jiǎn)介、參數(shù)取值參考，并對(duì)有限元軟件使用進(jìn)行必要的介紹，以便于學(xué)生們學(xué)習(xí)和使用。所添加內(nèi)容主要包括如下三個(gè)部分。

1.計(jì)算模型分析簡(jiǎn)介

對(duì)平板載荷試驗(yàn)、三軸剪切試驗(yàn)、恒定滲流、邊坡穩(wěn)定性和被動(dòng)土壓力五個(gè)土力學(xué)中的經(jīng)典試驗(yàn)進(jìn)行了必要總結(jié)。主要介紹了相應(yīng)的實(shí)驗(yàn)?zāi)康暮褪覂?nèi)或原位試驗(yàn)方法，并對(duì)相應(yīng)的土力學(xué)計(jì)算理論和方法進(jìn)行了簡(jiǎn)要的描述。對(duì)于仿真試驗(yàn)，以淺層平板載荷試驗(yàn)為例，本軟件簡(jiǎn)要說(shuō)明了實(shí)驗(yàn)?zāi)康?，然后闡述了其原位試驗(yàn)方法，可使學(xué)生在使用本軟件的同時(shí)再次復(fù)習(xí)書本知識(shí)，達(dá)到溫故而知新的目的。對(duì)于仿真分析，如恒定滲流分析，本模塊首先定義恒定滲流，再介紹常用的計(jì)算理論與經(jīng)驗(yàn)算法，以便學(xué)生提前預(yù)知恒定滲流仿真分析模塊所能實(shí)現(xiàn)的功能。

2.計(jì)算參數(shù)的取值建議

軟件2.0為用戶嵌入的“參數(shù)含意與取值”這一選項(xiàng)卡，學(xué)生可以隨時(shí)查看各個(gè)模塊中涉及的各計(jì)算模型參數(shù)的含義和取值范圍。其中的土性參數(shù)，參照《工程地質(zhì)手冊(cè)》按照砂土、粉土、粉質(zhì)黏土和一般黏土分類給出了相應(yīng)的經(jīng)驗(yàn)取值范圍，或者提供了可參考的經(jīng)驗(yàn)計(jì)算公式。因?qū)W生缺乏實(shí)踐經(jīng)歷，對(duì)土性參數(shù)理解不到位，其數(shù)據(jù)選取更加不明確，使得仿真試驗(yàn)或分析結(jié)果不具有實(shí)際意義。藉于此部分計(jì)算參數(shù)的取值建議，為普遍缺乏有限元計(jì)算經(jīng)驗(yàn)的本科學(xué)生提供了可靠借鑒。

3.有限元介紹

為使本科學(xué)生了解本程序設(shè)計(jì)思路與算法，加深對(duì)仿真試驗(yàn)和分析的認(rèn)識(shí)，同時(shí)深入淺出地給部分本科學(xué)生提供一個(gè)初步了解有限元思想的入口，程序2.0版本增加了對(duì)有限元相關(guān)知識(shí)的簡(jiǎn)單介紹，主要包括對(duì)有限元基本思路以及相應(yīng)模塊有限元實(shí)現(xiàn)方法的介紹。通過(guò)對(duì)有限元的介紹，使本科學(xué)生初步理解有限元思想，并進(jìn)一步提高深入學(xué)習(xí)的興趣。也使得學(xué)生對(duì)有限元思想及其運(yùn)用有初步認(rèn)識(shí)，破除了其今后學(xué)習(xí)如ABAQUS、 PLAXIS等大型通用軟件時(shí)因不理解有限元理論而形成的心理障礙。

上述部分在程序中的相應(yīng)見(jiàn)面如圖2所示。

（二） 新增計(jì)算模塊

軟件2.0在保留原有軟件平板荷載試驗(yàn)和三軸試驗(yàn)兩個(gè)模塊的基礎(chǔ)上，新增加了恒定滲流分析、邊坡穩(wěn)定分析、被動(dòng)土壓力分析三個(gè)計(jì)算板塊。如下將簡(jiǎn)要介紹新增計(jì)算模塊的具體功能和實(shí)現(xiàn)效果。

在恒定滲流分析、邊坡穩(wěn)定分析、被動(dòng)土壓力分析三個(gè)計(jì)算模塊中，均可以選擇分析簡(jiǎn)介、參數(shù)含義與取值、仿真分析和退出以獲得不同的響應(yīng)。其中，“分析簡(jiǎn)介” 和“參數(shù)含義與取值”主要為了給出必要的介紹，在上文已有介紹?！胺抡娣治觥蹦K為程序的核心模塊，點(diǎn)擊“仿真分析”按鈕后程序?qū)棾觥皡?shù)”、“計(jì)算”和“輸出”三個(gè)子選項(xiàng)卡。其中在“參數(shù)”對(duì)話框中，需要使用者輸入需要的土性參數(shù)和相關(guān)模型幾何參數(shù)（各參數(shù)意義及經(jīng)驗(yàn)取值可從“參數(shù)含義與取值”選項(xiàng)中獲得）；在“計(jì)算”選項(xiàng)卡中，使用者可以選擇不同的控制計(jì)算誤差、計(jì)算網(wǎng)格精度和最大迭代次數(shù)，考慮使用者對(duì)這些參數(shù)的合理選擇需要對(duì)非線性有限元知識(shí)有必要的了解，在程序應(yīng)用界面中給出了相應(yīng)的解釋和推薦值；在“輸出”選項(xiàng)卡中，使用者可根據(jù)需要選擇輸出包含不同數(shù)據(jù)的計(jì)算圖形用于后續(xù)分析。

以恒定滲流分析為例，進(jìn)入“恒定滲流”計(jì)算模塊后，選擇仿真分析選項(xiàng)。程序提醒使用者在“參數(shù)”對(duì)話窗口中輸入必要的物理和模型參數(shù)等：水平滲透系數(shù)、豎直滲透系數(shù)和包括模型的左、右側(cè)水頭的邊界條件，以及包括滲透深度、基坑寬度和止水帷幕深度在內(nèi)的模型幾何參數(shù)。在“計(jì)算”對(duì)話框中，使用者可選擇較粗、中等和較細(xì)三種不同精度的計(jì)算網(wǎng)格精度。在“輸出”對(duì)話框中，使用者則可以選擇輸出計(jì)算域內(nèi)水頭勢(shì)分布圖，還可通過(guò)打開(kāi)程序所在文件夾中的相應(yīng)文件獲取每節(jié)點(diǎn)位勢(shì)和基坑每日涌水量。程序中上述相應(yīng)的參數(shù)輸入和結(jié)果輸出界面見(jiàn)圖3。

根據(jù)輸入的土體參數(shù)，經(jīng)過(guò)計(jì)算，可以選擇較稀疏或密集的水位勢(shì)圖，根據(jù)不同選項(xiàng)得到如圖4所示的兩種仿真分析圖。

圖5和圖6分別給出另外兩個(gè)計(jì)算模塊的參數(shù)輸入界面和計(jì)算結(jié)果輸出界面，由于其界面上信息具有自明性，不再贅述。

二、軟件應(yīng)用

本文所開(kāi)發(fā)的軟件代碼開(kāi)源，免安裝、文件小且易操作，能較好地滿足本科土力學(xué)實(shí)驗(yàn)?zāi)M、課堂教學(xué)和高階學(xué)習(xí)。在軟件1.0基礎(chǔ)上，新增了三個(gè)計(jì)算模塊，并增強(qiáng)了軟件的易用性和容錯(cuò)性。實(shí)踐表明，改進(jìn)后的軟件2.0可以達(dá)到如下三點(diǎn)應(yīng)用目標(biāo)。

（一） 虛擬實(shí)驗(yàn)?zāi)M

淺層平板載荷試驗(yàn)作為地質(zhì)勘查中最常使用的原位試驗(yàn)方法，用于測(cè)量土的變形模量、地基極限承載力等重要的土力學(xué)性質(zhì)，而擋土墻土壓力問(wèn)題和邊坡穩(wěn)定性問(wèn)題均是經(jīng)典的土力學(xué)問(wèn)題。除此之外，滲流問(wèn)題也是土力學(xué)的重點(diǎn)教學(xué)內(nèi)容之一。在目前土力學(xué)課程教學(xué)中，尚難以普遍展開(kāi)上述問(wèn)題的試驗(yàn)；即便少數(shù)學(xué)校能夠有條件進(jìn)行該類實(shí)驗(yàn)，但由于試驗(yàn)影響因素較多，投入較多，難以保證實(shí)驗(yàn)效果。而本文開(kāi)發(fā)的土力學(xué)有限元仿真軟件2.0可實(shí)現(xiàn)對(duì)上述經(jīng)典土力學(xué)問(wèn)題試驗(yàn)的仿真模擬，以數(shù)值模擬的方式彌補(bǔ)目前上述的土力學(xué)實(shí)踐環(huán)節(jié)的不足。

（二） 改善教學(xué)效果

土力學(xué)課程教學(xué)中，恒定滲流分析、邊坡穩(wěn)定分析以及土壓力的計(jì)算往往是授課重點(diǎn)，但由于這三個(gè)土力學(xué)問(wèn)題涉及到的模型假設(shè)多，計(jì)算方法復(fù)雜，嚴(yán)重降低了學(xué)生的學(xué)習(xí)興趣。本軟件可以作為課堂教學(xué)的補(bǔ)充，以被動(dòng)土壓力計(jì)算為例，教師課堂上可為學(xué)生講解經(jīng)典的朗肯土壓力理論和摩爾庫(kù)倫理論，介紹其計(jì)算方法，并在課堂上演示本軟件的操作使用，通過(guò)輸出被動(dòng)土壓力—位移曲線，幫助學(xué)生理解土壓力的分析過(guò)程，并通過(guò)計(jì)算后的變形網(wǎng)格讓學(xué)生直觀地感受被動(dòng)土壓力對(duì)土體的影響。

（三）高階學(xué)習(xí)輔助

本文開(kāi)發(fā)的土力學(xué)仿真軟件2.0版本具有較好的易用性，這進(jìn)一步激發(fā)了本科學(xué)生的學(xué)習(xí)興趣，加深學(xué)生對(duì)基本土力學(xué)理論知識(shí)的理解。通過(guò)系統(tǒng)變換土性參數(shù)和試驗(yàn)條件，利用該軟件的模擬結(jié)果，可以與書本上的相應(yīng)經(jīng)典土力學(xué)計(jì)算理論結(jié)果進(jìn)行對(duì)比，通過(guò)這些對(duì)比可以研究相應(yīng)理論中的簡(jiǎn)化假設(shè)的影響。與此同時(shí)，可以引導(dǎo)學(xué)生在土力學(xué)及其相關(guān)知識(shí)上進(jìn)行更高階的學(xué)習(xí)。

三、結(jié)語(yǔ)

在前期開(kāi)發(fā)的土力學(xué)本科教學(xué)有限元軟件1.0的基礎(chǔ)上，筆者開(kāi)發(fā)的軟件2.0版本新增了恒定滲流等三個(gè)計(jì)算模塊，使之能夠分析目前土力學(xué)課程所有重要的、實(shí)踐性很強(qiáng)但又無(wú)法在本科教學(xué)實(shí)現(xiàn)的經(jīng)典問(wèn)題。另外，本研究還補(bǔ)充了模型計(jì)算參數(shù)含義與推薦取值范圍的必要介紹，提高了軟件的易用性。目前的軟件能夠?yàn)楫?dāng)前土力學(xué)課程教學(xué)提供虛擬的實(shí)踐資源，彌補(bǔ)了課堂實(shí)踐教學(xué)環(huán)節(jié)的不足，能更好地激發(fā)了學(xué)生的學(xué)習(xí)興趣。實(shí)踐表明，該程序能夠作為本科土力學(xué)課堂教學(xué)的有益補(bǔ)充，以及高階學(xué)習(xí)的輔助工具，可以提高學(xué)生對(duì)課程知識(shí)的掌握程度，達(dá)到預(yù)期的教學(xué)效果。

參考文獻(xiàn)：

[1]東南大學(xué)，浙江大學(xué)，湖南大學(xué)，等.土力學(xué)[M].北京：中國(guó)建筑工業(yè)出版社，2008.

[2]王常明，王清，范建華，等.計(jì)算機(jī)仿真在土力學(xué)實(shí)驗(yàn)教學(xué)中的應(yīng)用[J].高等建筑教育，2005，14（4）：96-98.

[3]張昭，郭呈周，肖昭然.基于修正劍橋模型的三軸試驗(yàn)仿真模擬[J].巖土力學(xué)，2007，28（9）：1949-1953.

[4]Smith I.M.， Griffiths D.V.有限元方法編程[M].北京：電子工業(yè)出版社，2003.

[5]唐章宏，薛賽男，馮峰，等.Visual Fortran程序設(shè)計(jì)[M].北京：人民郵電出版社，2000.

[6]商翔宇，鄭秀忠，梁恒昌.土力學(xué)本科教學(xué)有限元仿真軟件的開(kāi)發(fā)與應(yīng)用[J].高等建筑育，2014（23）：136-139.