成礦動力學(xué)數(shù)值模擬(第二講)

成礦動力學(xué)數(shù)值模擬

NumericalModellingoftheMetallogenicGeodynamicalProcesses第二講提綱一、數(shù)值模擬的基礎(chǔ)動力學(xué)數(shù)值模擬及其特點怎樣進行數(shù)值模擬二、FLAC軟件的應(yīng)用FLAC的理論基礎(chǔ)和主要特點應(yīng)用FLAC進行數(shù)值模的步驟FLAC的主要命令FLAC自帶的FISH語言FLAC模擬技巧一、動力學(xué)數(shù)值模擬的基礎(chǔ)動力學(xué)數(shù)值模擬及其特點怎樣進行數(shù)值模擬1.動力學(xué)數(shù)值模擬及其優(yōu)缺點理論上說，了解成礦的動力學(xué)過程應(yīng)該有實驗的(experimental)和數(shù)學(xué)的(mathematical)兩種方法，但事實事上，由于極端物理、化學(xué)條件和巨大的時空尺度，我們無法在實驗室以物理和化學(xué)的手段再現(xiàn)實際的成礦過程，而數(shù)學(xué)的方法卻不受這些條件限制；動力學(xué)的核心數(shù)學(xué)問題是動力學(xué)方程組，主要由平衡方程(能量、質(zhì)量和動量)和系統(tǒng)的本構(gòu)方程組成，實際上是一組復(fù)雜的偏微分方程組成，理論上說研究這些方程組應(yīng)該有兩種不同的方法，一是解析法(或稱分析法,analyticmethod)，另一種是數(shù)值法(numericalmethod);

解析方法就是用數(shù)學(xué)上的積分法或積分變換等方法直接求數(shù)學(xué)模型的解，其解稱為解析解，它是數(shù)學(xué)模型的精確解，這種解的最大優(yōu)點在于把表征物質(zhì)運動規(guī)律的變量與激發(fā)條件、時空變化包含在一個表達式中，這樣便于分析物質(zhì)運動規(guī)律。但實際上由于動力學(xué)方程的復(fù)雜性和系統(tǒng)幾何結(jié)構(gòu)和形態(tài)的復(fù)雜性，大多數(shù)情況下，我們是不可能求得分析解的；數(shù)值模擬(numericalmodeling)就是離散方法求解數(shù)學(xué)模型，或者說用數(shù)值方法求解動力學(xué)方程組的數(shù)值解，其解為數(shù)值集合，是數(shù)學(xué)模型的近似解。盡管不如解析法精確，但它對復(fù)雜的動力學(xué)系統(tǒng)而言，是一種切實可行的方法，實際上對地質(zhì)研究而言，其精度已足夠高；數(shù)值模擬又稱數(shù)值實驗(numericalexperiment),實際上它是在計算機的虛擬空間內(nèi)定量地再造某一自然事件或回放真實的過程(quantitativereconstructionofanaturaleventorareplayofa“real-life”process)。數(shù)值模擬的方法分為：有限元(limitedelement)法、有限差分(limiteddifference)法和邊界元法(boundaryelement)等，不管是哪種方法，都得將一個連續(xù)的系統(tǒng)分割成n單元，如果整個系統(tǒng)有m個獨立的變量的話，那么每算一步，就要解m×n個方程。如果系統(tǒng)尺度很大，而我們要保證模擬的準(zhǔn)確度，那么n就必須足夠大，這樣就導(dǎo)致計算工作量相當(dāng)大，在計算機技術(shù)不發(fā)達時，這對數(shù)值模擬來說是一個相當(dāng)大的難題。就現(xiàn)代計算機技術(shù)水平及發(fā)展趨勢來看,這已經(jīng)不是一個大的問題。2.怎樣進行數(shù)值模擬進行數(shù)值模擬的6個關(guān)鍵步驟Sixkeycomponentsformakingamodel問題

A“Story”orkeyquestionThisisaproblemorscenarioyouwanttoexploreingeologicaloranyotherdiscipline.e.g.Canshearbanddevelopinarockblocksubjectedtoshearing?軟件

Anumericalcodetosimulatethisstory

YoucaneitherwriteyourcodesoruseexistentcodesComputercodesmostlywrittenusingFiniteelementorfinitedifferencemethods.3)幾何模型

Constructionofgeometry(mesh)

Geometricalstructuresareapproximatedbyameshthatmayincludeinternalstructuralelements.4)

本構(gòu)關(guān)系

Selectionofappropriaterheology(constitutivelaws)

e.g.elastic;elastic-plastic;viscous;conductiveheattransportetc,eachrepresentsacertainapproximationofrockbehaviours5)模型組成單元的特征參數(shù)

Specificationofrelevantrockproperties(parameters)foryourmodel6)邊界條件和初始條件

Definingnecessaryboundaryconditionsforthemodel–reflectionofnaturaldeformationorthermalenvironment

e.g.

Initialstatus;boundarydisplacementrates;boundarytemperatureetc.7)模型運算及結(jié)果輸出

Runthemodel–computationaliteration

generally,a“time-marching”forwardmodelling二、FLAC軟件的應(yīng)用FLAC的理論基礎(chǔ)及主要特點應(yīng)用FLAC進行數(shù)值模的步驟FLAC的主要命令FLAC自帶的FISH語言FLAC模擬技巧1.FLAC的理論基礎(chǔ)及主要特點FLAC是FastLagrangianAnalysisofContinua的縮寫，美國Itasca公司開發(fā)的,意為連續(xù)介質(zhì)的快速拉格朗日法分析,它是一種以顯式有限差分法(explicitfinitedifference)

進行連續(xù)介質(zhì)物理作用模擬的商業(yè)軟件(commercialcode)。它源自數(shù)學(xué)力學(xué)的拉格朗日元法。最初由Minnesota大學(xué)的PeterCundall開發(fā)的。目前，F(xiàn)LAC已由二維發(fā)展到三維，二維計算程序V3.0以前的為DOS版本，V2.5版本僅僅能夠使用計算機的基本內(nèi)存（64K），所以，程序求解的最大結(jié)點數(shù)僅限于2000個以內(nèi)。1995年，F(xiàn)LAC2D已升級為V3.3的版本，其程序能夠使用護展內(nèi)存，大大發(fā)護展了計算規(guī)模。最新發(fā)布的FLAC已到了5.0版本FLAC的基本原理和算法與離散元相似，但它卻像有限元那樣適用于多種材料的模擬和邊界條件非規(guī)則區(qū)域的連續(xù)問題求解。在求解過程中，F(xiàn)LAC采用了離散元的動力態(tài)松弛法，不需要求解大型聯(lián)立方程組（無需形成剛度矩陣）；另一方面，同以往的差分分析法相比，F(xiàn)LAC不但可以對連續(xù)介質(zhì)進行大變形分析，而且還可以模擬巖土體沿某一軟弱面產(chǎn)生的滑動變形；FLAC還能針對不同的材料，實現(xiàn)相應(yīng)的本構(gòu)方程來比較真實地反映實際材料的動態(tài)行為。程序?qū)⒂嬎銋^(qū)域內(nèi)的介質(zhì)劃分為若干個二維單元,單元之間用節(jié)點相互連接。上述過程中,網(wǎng)格的劃分與有限元不同之處在于其網(wǎng)格分物理網(wǎng)格和數(shù)學(xué)網(wǎng)格且互為影射,所劃分的網(wǎng)格只要有序也可具有不規(guī)則的形狀。計算循環(huán)圖所示.拉格朗日差分法計算循環(huán)FLAC的主要特點1)

對硬件配置較低由于FLAC采用的是顯示有限差分法,在內(nèi)存較小的低檔機上亦可進行較大規(guī)模的計算,16MB時,可計算60000的單元數(shù);(2)強大的前后處理功能

FLAC具有很強的前后處理功能。只要設(shè)置某些控制點的坐標(biāo),軟件就以自動生成計算網(wǎng)格,。用戶可以根據(jù)實際情況通過某些命令修改網(wǎng)格。各階段的計算結(jié)果均可以數(shù)據(jù)文件的形式存盤,一旦需要,可用Restart命令恢復(fù)全部現(xiàn)場,使用起來非常方便。

(3)實現(xiàn)對多種材料和多種過程的模擬

可以模擬彈性模型,莫爾-庫侖模型,橫觀各向同性、遍有節(jié)理模型,應(yīng)變硬化和軟化等多種材料模型。此外該軟件還提供了交界面模型,用戶可以用滑動面來模擬斷層和節(jié)理。還可以根據(jù)實際情況采用某一種模型,也在計算范圍內(nèi)定義若干子區(qū)域,賦予不同的材料不同參數(shù)值,以模擬復(fù)雜的地質(zhì)條件；2.應(yīng)用FLAC進行數(shù)值模的步驟編寫一個完整的FLAC模擬運算程序至少包括如下六個最基本的部分：1）構(gòu)建幾何模型：2）確立本構(gòu)關(guān)系：3）模型參數(shù)賦值：4）確立邊界條件和初始條件：5）模擬過程控制：6）模擬結(jié)果輸出：3.FLAC的主要命令FLAC有兩種基本的動行模式，command-driven和manu-driven，進行科學(xué)計算，一般都是用