摘要： 為解決傳統(tǒng)編程方法在實現(xiàn)非線性本構模型時的可維護性差、二次開發(fā)困難等問題，基于有限狀態(tài)機（FSM）原理，提出狀態(tài)模式驅動和數(shù)表驅動的兩種編程實現(xiàn)方法；通過對兩種編程方法的對比，認為FSM結合數(shù)表驅動的方法更適合編寫非線性本構程序；隨后對Davidenkov本構模型滯回曲線進行修正，提出兩種滯回曲線構造形式：指向失效點型和指向修正點型。改進的Davidenkov本構模型糾正了“n倍法”模型在失效后，剪應變反向達到上限，而剪應變未進入失效狀態(tài)的缺點，使得改進的Davidenkov本構模型在模擬土動應力-應變關系時更加合理，尤其是指向失效點模型在多次加卸載作用下，對土體剪切剛度的模擬情況更為精確;為驗證基于FSM編制的本構程序具有可維護性強、二次開發(fā)簡單等特點，依據(jù)FSM結合數(shù)表驅動的方法編制了Davidenkov及其修正的本構模型程序，以驗證該方法的有效性與正確性。

關鍵詞： 巖土動力學; 非線性本構模型; 有限狀態(tài)機; Davidenkov本構模型; 二次開發(fā)

中圖分類號： TU435""""" 文獻標志碼：A"" 文章編號： 1000-0844（2025）02-0351-10

DOI：10.20000/j.1000-0844.20230802002

Implementation of the Davidenkov constitutive model

based on finite state machine

DONG Zhengfang1， LI Haojie1， WANG Yongfeng1， JIN Deyin2

（1.School of Civil Engineering and Architecture， Henan University， Kaifeng 475004， Henan， China;

2.JAPAN-CHINA CONSALTANT Co.， Ltd.， Tokyo 134-0086， Japan）

Abstract：

Two methods based on finite state machine （FSM）， namely state model-driven and table-driven， were investigated in this study to solve the problems of poor maintainability and difficult secondary development in the implementation of nonlinear constitutive models using traditional programming methods. Through comparison， the results reveal that the table-driven method can efficiently describe nonlinear constitutive models. Then， the stress-strain hysteretic curve of the Davidenkov constitutive model was modified， from which two construction forms of the hysteretic curve were proposed： failure point and correction point type. The results reveal that the modified Davidenkov model can correct the shortcoming of the “n-fold method” model， in which the shear strain reversely reaches the upper limit shear strain and is unable to enter the failure state， so it is more reasonable to be used in simulating the dynamic stress-strain relationship of soil， especially using the failure point model to simulate the shear stiffness of soil under cyclic loading and unloading. Based on the programming method of FSMs and the table-driven method， the constitutive programs of the Davidenkov model and its modified models were compiled， thereby verifying the effectiveness and correctness of the proposed method.

Keywords：

geotechnical dynamics; nonlinear constitutive model; finite state machine; Davidenkov constitutive model; secondary development

0 引言

土體的動應力-應變關系具有嚴重的非線性、滯后性和變形累積性等特點，其本構理論大致可分為黏彈性和彈塑性兩種，黏彈性理論因其直觀簡單、參數(shù)少且易獲取等優(yōu)點而成為目前研究和應用的主流理論^［1］。黏彈性模型包含等效線性和曼辛非線性兩種，后者在前者的基礎上由曼辛法則構造出滯回曲線方程，使得這類模型更接近土體動應力-應變曲線^［2-3］。實際上，巖土體本構模型數(shù)量較多且各自具有不同的應用范圍。其中由曼辛法則構造的三參數(shù)Davidenkov本構模型因能較好地模擬各類土體而最具代表性，特別是對軟土^［4-6］、砂土^［7］、泥炭質(zhì)土^［8］等在動荷載作用下的非線性行為模擬。Davidenkov本構模型最早由Martin等^［9］提出，隨后Pyke^［10］以“n倍法”代替“2倍法”限制了模型后繼的滯回圈并控制了滯回曲線的發(fā)展方向;陳國興等^［7］引入上限剪應力對Davidenkov骨架曲線進行修正，使模型曲線與上限剪應力的水平線相遇后可沿水平線方向發(fā)展;趙丁鳳等^［11］在原有骨架曲線的基礎上，提出滯回曲線始終沿當前荷載轉向點指向歷史最值點的前進方式，有效解決了在編程時模型信息存儲量大且具有不確定性的問題，但其并未考慮對骨架曲線進行限制，模型仍然存在剪應力能夠無限增長的不合理的問題。

基于Davidenkov本構模型的出色表現(xiàn)，近些年眾多研究人員在不同有限元平臺對其進行了二次開發(fā)與研究。王沿朝等^［4］、陳國興等^［7］基于ABAQUS平臺進行實現(xiàn)并擬合了相關土體參數(shù);陳斌等^［5］基于ANSYS平臺進行二次開發(fā)得到了深厚軟土的相關參數(shù);張如林等^［6］、王國波等^［12］基于FLAC^3D平臺進行實現(xiàn)并與試驗實測結果進行對比，驗證了Davidenkov本構模型的適用性;趙旭清^［13］基于OpenSees平臺進行實現(xiàn)，分析了黏土在地震作用下的動力響應。針對本構程序的二次開發(fā)，不同研究人員擁有不同的算法流程，并未形成統(tǒng)一的編程實現(xiàn)方法;此外，傳統(tǒng)編程方法在實現(xiàn)復雜非線性本構程序時，通常要對模型的所有情況進行遍歷，并通過大量嵌套if.....else或switch.....case等判斷語句進行實現(xiàn)。因此，模型越復雜，判斷語句越多，最終造成本構程序的邏輯可讀性下降，后續(xù)研究人員進行維護、二次開發(fā)時就越困難。因此，有必要尋求一種便捷、高效且能夠實現(xiàn)本構程序“模板化”開發(fā)的編程方法。

有限狀態(tài)機（Finite State Machine，F(xiàn)SM）作為描述事物邏輯轉換關系的一種數(shù)學模型，其應用十分廣泛。例如：自動駕駛汽車預期功能安全危害識別^［14］、運動機器人編程^［15-16］、系統(tǒng)攻防^［17］、集成電路^［18］以及電力系統(tǒng)^［19］等領域。FSM的特點在于能夠將復雜邏輯關系進行抽象化處理，而非線性本構模型在不同受力狀態(tài)下的應力-應變關系具有靈活多變的特點，因此基于FSM的編程方法在實現(xiàn)非線性本構程序時相較于傳統(tǒng)編程方法具有更簡單、更高效的特點。Jin等^［20］將FSM應用到雙折線本構模型的開發(fā)中，發(fā)現(xiàn)FSM的編程相較于傳統(tǒng)編程，在提升編程效率、規(guī)避程序錯誤、提升代碼復用性等方面具有顯著提升。但對于程序的具體實現(xiàn)其并未進行詳細闡述，相關方面的研究也受制于編程方法的抽象性而較為缺乏。

本文基于有限狀態(tài)機理論提出FSM結合狀態(tài)模式和FSM結合數(shù)表的兩種編程方法，給出兩種編程方法的優(yōu)劣對比及實現(xiàn)步驟，為本構開發(fā)人員提供“模板化”的編程思路。此外，為驗證基于FSM編制的本構程序在二次開發(fā)時的便捷性，本文對“n倍法”Davidenkov本構模型的骨架曲線進行限制，并對其進行修正，提出兩種失效模型：指向失效點型和指向修正點型。通過FSM結合數(shù)表的編程方法對上述模型進行實現(xiàn)，以驗證所提方法的有效性。

1 有限狀態(tài)機理論

有限狀態(tài)機（FSM）可以將事物復雜的邏輯關系抽象為有限個穩(wěn)定狀態(tài)間的轉換關系。圖1為雙折線本構模型的FSM模型，模型狀態(tài)可劃分為3種：S、SP和SN;每種狀態(tài)對應一個函數(shù)行為，分別為Y、Y1、Y2。圖2為雙折線本構模型的狀態(tài)遷移圖，詳細描述了不同狀態(tài)之間的轉換關系。由圖2可以看出，事件作為狀態(tài)遷移的信號負責狀態(tài)更新，新、舊狀態(tài)通過事件建立聯(lián)系，狀態(tài)的更新必然引發(fā)對應動作的執(zhí)行。

1.1 狀態(tài)機結合狀態(tài)模式驅動

狀態(tài)模式本質(zhì)上是一種設計模式，核心思想在于將事物復雜的邏輯關系封裝成不同的對象，進而降低不同對象間的耦合，方便程序的二次開發(fā)。

FSM結合狀態(tài)模式的編程方法不僅可以降低不同狀態(tài)間的耦合，還可以實現(xiàn)不同狀態(tài)間的定向遷移。編程方法如下：

（1） 根據(jù)模型特點枚舉出模型所有狀態(tài)及觸發(fā)狀態(tài)轉換的所有事件。

（2） 建立狀態(tài)控制類，主要負責模型不同狀態(tài)間的切換與函數(shù)執(zhí)行，將不同狀態(tài)子類進行鏈接。

（3） 建立抽象狀態(tài)類，通過抽象狀態(tài)類定義狀態(tài)接口、封裝函數(shù)行為。一個狀態(tài)類可以對應多個函數(shù)行為，針對本構而言，大致可分為狀態(tài)切換與應力計算兩個行為。

（4） 根據(jù)模型狀態(tài)數(shù)目編寫具體狀態(tài)類。具體狀態(tài)類用于實現(xiàn)抽象狀態(tài)類所對應的具體函數(shù)行為，并且負責在有需求的情況下進行狀態(tài)更新。

以雙折線本構模型為例，給出FSM結合狀態(tài)模式所編制程序的UML關系如圖3所示。

FSM結合狀態(tài)模式編制的本構程序將本構模型的各個狀態(tài)分散至不同的狀態(tài)子類，達到了高內(nèi)聚、低耦合的程序執(zhí)行效果。因此，本構程序在進行二次開發(fā)（FSM元素數(shù)量變化）或維護時，研究人員只需將對應的狀態(tài)子類進行增減或修改即可。

1.2 狀態(tài)機結合數(shù)表驅動

基于狀態(tài)模式驅動的本構程序，本構關系分散在各個狀態(tài)子類中，研究人員無法直觀地看出整個模型的邏輯關系，容易引發(fā)邏輯分散的問題。若將程序的邏輯關系以數(shù)組的形式進行表征，其可讀性和可維護性將變得十分可觀。本構程序實現(xiàn)的另一種方式就是建立一個描繪狀態(tài)遷移的數(shù)表，通過FSM與數(shù)表結合的方式快速實現(xiàn)復雜本構程序的編寫。編程方法如下：

（1） 根據(jù)模型特點枚舉出模型所有狀態(tài)及觸發(fā)狀態(tài)轉換的所有事件;

（2） 定義不同狀態(tài)的函數(shù)行為，針對本構模型而言，主要為應力函數(shù)的計算;

（3） 通過結構體數(shù)組實現(xiàn)不同狀態(tài)間遷移關系。

若想進一步使得本構程序的封裝性更強、模塊化程度更高，可以將本構程序拆分為模型子類、狀態(tài)機屬性類（FSMItem）和狀態(tài)機類（FSM），三者互為友元。其中，F(xiàn)SMItem對FSM參數(shù)進行定義，將本構模型中所有的狀態(tài)和事件進行枚舉;FSM則主要負責狀態(tài)切換及動作執(zhí)行，最終形成成員變量為狀態(tài)遷移表（FSMTable）。模型子類通過FSMTable將FSMItem與FSM聯(lián)系起來。

以雙折線本構模型為例，給出基于FSM結合數(shù)表的編程方法所編制程序的流程如圖4所示。

FSM結合數(shù)表編制的本構程序通過二維數(shù)表的形式建立起不同狀態(tài)、事件間的聯(lián)系;通過成員函數(shù)的形式定義不同狀態(tài)的行為函數(shù)。因此，本構程序在進行維護或二次開發(fā)時只需對二維數(shù)表進行修改或增減對應模型元素與行為函數(shù)即可。

1.3 兩種方法對比

基于FSM結合狀態(tài)模式的編程方法和基于FSM結合數(shù)表的編程方法，在實現(xiàn)本構程序時均提供“模板化”的編程思路。兩種編程方法均可通過拆分模型屬性;引入FSM來控制模型遷移關系，從而避免了傳統(tǒng)編程實現(xiàn)的過程中大量嵌套判斷語句的現(xiàn)象，提升了代碼的可讀性與可維護性。兩種編程方法在程序實現(xiàn)方面各具優(yōu)劣，具體如表1所列。

本構方程通常為若干數(shù)學函數(shù)的組合，所以本構曲線大多呈現(xiàn)非線性、走向隨機、形式多樣等特點，用FSM模型可描述為模型狀態(tài)多且隨機，但運算邏輯簡單。因此，采用FSM結合數(shù)表的編程方法更適合本構程序的開發(fā)。為清楚表達本構程序的邏輯轉換關系與二次開發(fā)的便捷性，本文采用FSM結合數(shù)表的編程方法對Davidenkov及其修正本構模型進行實現(xiàn)。

2 Davidenkov及其滯回曲線構造

Hardin等^［21］最早提出動剪切模量比的概念，隨后Martin等^［9］在此基礎上，為更好地模擬各類土體的動剪切模量比，提出具有3個參數(shù)的Davidenkov本構模型，但模型存在剪應力能夠無限增長的不合理現(xiàn)象。為此，陳國興等^［7］在原有本構模型的基礎上，用分段函數(shù)表示不同應變區(qū)間下剪應力的函數(shù)對應關系，使得本構模型的剪應力不隨剪應變的增長而無限增長。此外，Davidenkov本構模型在編程實現(xiàn)時需預留許多的狀態(tài)變量來記憶荷載轉向點。針對該問題，趙丁鳳等^［11］參照“n倍法”提出了修正的Davidenkov本構模型（本文稱“n倍法”Davidenkov本構模型），使得模型的滯回曲線始終由當前荷載轉向點指向歷史最值點，降低信息存儲的同時方便了程序實現(xiàn)。但“n倍法”Davidenkov本構模型并未考慮對模型剪應力進行限制，模型仍然存在剪應力能夠無限增長的不合理現(xiàn)象。因此，本文擬對“n倍法”修正的Davidenkov本構模型的骨架曲線進行限制，使其更加滿足土體動應力-應變關系。

“n倍法”修正的Davidenkov本構模型骨架曲線表達式為：

τ=Gmaxγ［1-H（γ）］ （1）

其中：

H（γ）=γγ0^2B1+γγ0^2BA （2）

式中：A、B、γ0為模型擬合參數(shù);Gmax為初始剪切模量;γ、τ分別為動剪應變、動剪應力。對式（1）中的變量γ求導，得到骨架曲線剪切模量：

G=Gmax1-1+2ABγ^2B0γ^2B0+γ^2BH（γ） （3）

滯回曲線始終指向最值點，表達式為：

τ=τc+Gmax（γ-γc）1-Hγ-γc2n（4）

式中：γc、τc分別為荷載轉向點處剪應變與剪應力。滯回曲線剪切模量：

G=Gmax1-1+2AB（2nγ0）^2B（2nγ0）^2B+γ-γc^2B·Hγ-γc2n（5）

其中（2nγ0）^2B由當前荷載轉向點及歷史最大（最?。c確定：

（2nγ0）^2B=（γex±γc）^2B1-RR （6）

R=1-τex±τcGmax（γex±γc）^1A （7）

式中：γex、τex分別為荷載轉向點處剪應變與剪應力（注：符號“±”在加載時取“-”，卸載時取“+”）。

由式（3）、（5）可知，剪切模量的確定與等效剪應變相關，所以基于增量形式的優(yōu)化等效剪應變算法^［11］如下：

γ^t+Δt=γt+signΔγ^t+Δtincre （8）

Δγ^t+Δtincre=γ^t+Δtgen-γtgen （9）

γtgen=43Jt2（eoij） "（10）

式中：Jt2（eoij）為用張量eoij表示的應變偏量第二不變量。

eoij=eij-eij，c （11）

式中：eij為當前應變張量;eij，c為荷載轉向點應變張量。此時模型可依據(jù)Δγ^t+Δtincrelt;0直接判斷土體發(fā)生加卸載轉變，一旦發(fā)生加卸載變化eij，c將會更新，γtgen將會清零重新計算。

基于“n倍法”改進的Davidenkov本構模型的應力-應變曲線如圖5所示，模型剪應力存在無限增長的情況違背了大多巖土材料遵循的先硬化、后軟化的規(guī)律。實際上，若不考慮土體軟化，土體都存在一個上限剪應力，當土體剪應力大于上限剪應力時，土體就會破壞。不考慮土體軟化，對“n倍法”Davidenkov本構模型骨架曲線進行限制，修正后的骨架曲線為：

τ（γ）=Gmaxγ［1-H（γ）］γc≤γult

Gmaxγult［1-H（γult）］γcgt;γult（12）

τult=Gmaxγult［1-H（γult）］ （13）

加卸載曲線為：

τ（γ）=

τc+Gmax（γ-γc）1-Hγ-γc2n，τ≤τult

±τult，τgt;τult（14）

式中：γult、τult分別為失效點處剪應變與剪應力?；凇皀倍法”構造的失效Davidenkov本構模型，當0lt;γult≤γ時，模型曲線會沿水平線發(fā)展;卸載后模型遵循滯回曲線始終由荷載轉向點指向歷史最值點的方式前進。曲線只有當應變等于最大剪應變時，應力才等于上限剪應力，隨著荷載峰值的不斷增大，模型曲線呈現(xiàn)出圖6中（0→...→10）的走向。觀察圖6可知，模型曲線在多次加卸載以后耗能能力（一次加卸載曲線所圍面積）有降低趨勢，且存在應變反向到達失效剪應變而未進入失效的現(xiàn)象，這顯然與土體特性不符。

針對上述問題，本文采用修正最值點的方式對模型曲線進行修正，最值點確定方式修改為：

方式一：當最大剪應力τmax=τult時，模型卸載時的滯回曲線始終指向失效點（-γult，-τult）;模型卸載再加載時，滯回曲線逐步指向歷史最值點;當最大剪應力τmax=-γult，模型卸載時其滯回曲線指向失效點（γult，τult），卸載再加載其滯回曲線指向歷史最值點。最值點如表2所列。

對比曼辛法則構造的失效模型^［7］可知，曼辛法則構造的失效模型滯回曲線可能出現(xiàn)應變反向未達失效剪應變，而應力先達失效剪應力的現(xiàn)象。因此，本文提出第二種模型最值點確定方式，對模型滯回曲線進行修正。

方式二：當歷史最大剪應力τmax=τult，模型卸載時滯回曲線指向修正最值點;模型卸載再加載時，其滯回曲線指向歷史最值點;當歷史最大剪應力τmax=-γult，模型卸載時其滯回曲線指向修正最值點，卸載再加載時，其滯回曲線指向歷史最值點。最值點確定如表3所列。

3 模型實現(xiàn)與驗證

為驗證模型的適用性及FSM結合數(shù)表的編程方法在實現(xiàn)本構程序時的可行性與便捷性。本文先對“n倍法”改進的Davidenkov本構模型進行實現(xiàn)，隨后在該本構程序的基礎上對兩種修正本構模型進行二次開發(fā)，并對本構走向及剛度變化進行分析。

“n倍法”改進的Davidenkov本構模型的程序實現(xiàn)分為兩個關鍵步驟：（1）根據(jù)模型信息定義模型狀態(tài)（骨架曲線SC、滯回曲線HC）、事件（加載SP、卸載SR、超越歷史最值點SH）和動作函數(shù)（骨架方程和恢復力方程）;（2）通過結構體聲明狀態(tài)遷移表（表4），并通過二維數(shù)組實現(xiàn)狀態(tài)切換。程序對于最值點及荷載轉向點的判斷十分關鍵，基于FSM結合數(shù)表編制的本構子程序可以通過程序上一收斂步狀態(tài)與當前觸發(fā)事件來減少判斷過程，從而實現(xiàn)定向更新最值點及荷載轉向點。

基于FSM結合數(shù)表實現(xiàn)“n倍法”Davidenkov本構程序的流程示意如圖7所示。

為驗證程序的可靠性，本文設置一個1 m×1 m×1 m規(guī)則的一階六面體單元。該單元底部固定約束，頂部一節(jié)點施加如圖8所示的位移荷載^［8］。提取一個積分點處的應力應變，如圖9中實線所示。伴隨著荷載的變化“n倍法”改進型Davidenkov本構模型的加卸載曲線出現(xiàn)明顯的滯回現(xiàn)象，且滯回圈飽滿。隨著荷載的持續(xù)，單元的剪切剛度會發(fā)生衰減，當荷載發(fā)生轉向時單元剪切剛度會發(fā)生突變?？梢钥闯龌贔SM結合數(shù)組的編程方法所編制的改進型Davidenkov本構程序能正確反映土體單元的加卸載情況。

該本構程序的一個優(yōu)勢是：用戶可以在輸出程序不同時刻所處階段，通過與自變量進行對比來檢查程序是否編寫正確。將程序上一收斂步狀態(tài)和當前步觸發(fā)事件輸出，可獲得單元此刻所處階段。隨著荷載的變化，本構模型在原有狀態(tài)的基礎上會觸發(fā)不同事件，該本構程序則會由此進入新的狀態(tài)并執(zhí)行對應動作。

基于FSM編制本構程序的另一個優(yōu)勢在于研究人員在對本構程序進行維護、二次開發(fā)時，具有較高的可操作性與簡便性?；贔SM編制的本構程序在進行二次開發(fā)時，只需添加新增的狀態(tài)及遷移關系即可。下面以修正的Davidenkov本構模型為例，在原有本構程序的基礎上對其進行二次開發(fā)，修正Davidenkov本構模型的狀態(tài)遷移關系如表5所列，程序二次開發(fā)流程示意如圖7所示。

根據(jù)狀態(tài)遷移表編制修正的Davidenkov本構程序（表5），施加如圖8所示的位移和荷載，并假定土體失效剪應變γult=12。圖9中兩種虛線分別對應兩種修正模型的加卸載曲線。

對比“n倍法”所改進的Davidenkov本構模型可知，該本構模型的剪應變在達到上限后，隨著后續(xù)的增長，該模型的剪應力保持不變，曲線沿水平線前進。此外，對比以上兩種修正本構模型發(fā)現(xiàn)，當τmax=τult時，兩種修正的本構模型在卸載再加載階段均指向歷史最值點，兩者的區(qū)別在于模型卸載階段的不同走向，其中指向失效點的模型在卸載階段始終指向模型首次失效點;而“指向修正點”模型構造的滯回曲線會沿當前荷載轉向點指向修正點的前進方式前進（圖10），均符合模型修正效果。

隨著荷載峰值的不斷增大，兩種改進模型在經(jīng)歷多次“上骨架”曲線之后，模型應力-應變關系與剛度會產(chǎn)生顯著差異。施加圖11所示的位移荷載，對比兩種修正本構模型曲線及剛度變化。

指向失效點的模型依據(jù)表2所列完成最值點的信息更新，模型在卸載時始終指向初始失效點，多次加卸載之后應力-應變曲線呈現(xiàn)出圖11所示的變化趨勢。隨著拉平曲線的不斷延長，曲線在轉向以后的初始斜率越低，表明隨著最大剪應變的不斷增大，模型在荷載轉向后的初始剪切剛度越低;此外，隨著最大剪應變的不斷增大，指向失效點模型在反向加載至上限剪應變時模型的剛度越接近于0，極限狀態(tài)下（拉平曲線無限長）剪切剛度為0。模型剛度變化如圖11所示的剛度折損曲線，顯然指向失效點時，模型所描述的剛度變化趨勢與物理現(xiàn)象吻合;

指向修正點模型依據(jù)表3所列的最值點信息進行更新，模型荷載轉向時會根據(jù)歷史最值點對走向進行修正，在多次加卸載后曲線呈現(xiàn)“外擴”現(xiàn)象。隨著水平線的不斷延長，曲線在轉向后的初始斜率并無明顯變化，表明隨著最大剪應變的增大，指向修正點模型在荷載轉向后的初始剛度并無明顯變化;此外，當曲線到達失效點時的斜率也趨于定值，表明隨著最大剪應變的不斷增大，指向修正點模型在進入失效時的剛度趨于定值。顯然指向修正點模型并不具備反向達到失效點時，剪切剛度隨著最大剪應變的增長而趨于0的性質(zhì)，但其曲線走向與曼辛二倍法構造的失效模型曲線走向基本一致^［6］;同時，曲線呈現(xiàn)“等向外擴”的現(xiàn)象，改善了“n倍法”改進型Davidenkov本構模型滯回圈較為扁長的現(xiàn)象^［11］。

4 結論

本文基于有限狀態(tài)機結合數(shù)表的形式編制Davidenkov模型子程序并進行驗證，得出結論如下：

（1） 為復雜非線性本構模型的開發(fā)提供了兩種編程方法：有限狀態(tài)機結合狀態(tài)模式和有限狀態(tài)機結合數(shù)表。認為有限狀態(tài)機結合數(shù)表的編程方法更適合實現(xiàn)本構程序的開發(fā)，為“模板化”添加和維護本構程序提供了新的編程思路。

（2） 基于“n倍法”改進的Davidenkov模型提出了兩種滯回曲線構造形式：指向失效點型和指向修正點型。該模型有效地解決了原有Davidenkov模型在失效后反向達到上限剪應變而未失效的現(xiàn)象。

（3） 編制了Davidenkov模型及其修正模型子程序，驗證了基于有限狀態(tài)機的編程方法在快速實現(xiàn)復雜非線性本構程序的二次開發(fā)、規(guī)避復雜模型的程序錯誤等方面具有重要作用。

實際上，本文所提方法也適用于其他復雜非線性本構模型，例如復合材料、混凝土材料、橡膠材料等。抽象出的本構狀態(tài)越多、遷移越復雜，基于有限狀態(tài)機的編程方法就越方便、越高效。

參考文獻（References）

［1］ 馬曉文，梁慶國，趙濤，等.土動力學研究綜述及思考［J］.世界地震工程，2021，37（4）：217-230.

MA Xiaowen，LIANG Qingguo，ZHAO Tao，et al.Review and further consideration on research of soil dynamics［J］.World Earthquake Engineering，2021，37（4）：217-230.

［2］ 蔡袁強，于玉貞，袁曉銘，等.土動力學與巖土地震工程［J］.土木工程學報，2016，49（5）：9-30.

CAI Yuanqiang，YU Yuzhen，YUAN Xiaoming，et al.Soil dynamics and geotechnical earthquake engineering［J］.China Civil Engineering Journal，2016，49（5）：9-30.

［3］ 孔亮，鄭穎人，王燕昌.土體動本構模型研究評述［J］.寧夏大學學報（自然科學版），2001，22（1）：17-22，40.

KONG Liang，ZHENG Yingren，WANG Yanchang.Development in soil dynamic constitutive model［J］.Journal of Ningxia University （Natural Science Edition），2001，22（1）：17-22，40.

［4］ 王沿朝，陳清軍.基于粒子群優(yōu)化算法的軟土場地Davidenkov模型參數(shù)擬合與應用［J］.振動與沖擊，2019，38（17）：8-16.

WANG Yanchao，CHEN Qingjun.Parametric fitting of soft soil sites' Davidenkov model based on PSO algorithm and its application［J］.Journal of Vibration and Shock，2019，38（17）：8-16.

［5］ 陳斌，張棟梁，張春進，等.寧波深厚軟土動力特性研究及其地震響應分析［J］.地震工程學報，2016，38（3）：338-343.

CHEN Bin，ZHANG Dongliang，ZHANG Chunjin，et al.Analysis of dynamic characteristics and seismic response of deep soft soil in Ningbo area［J］.China Earthquake Engineering Journal，2016，38（3）：338-343.

［6］ 張如林，樓夢麟.基于達維堅科夫骨架曲線的軟土非線性動力本構模型研究［J］.巖土力學，2012，33（9）：2588-2594.

ZHANG Rulin，LOU Menglin.Study of nonlinear dynamic constitutive model of soft soils based on Davidenkov skeleton curve［J］.Rock and Soil Mechanics，2012，33（9）：2588-2594.

［7］ 陳國興，莊海洋.基于Davidenkov骨架曲線的土體動力本構關系及其參數(shù)研究［J］.巖土工程學報，2005，27（8）：860-864.

CHEN Guoxing，ZHUANG Haiyang.Developed nonlinear dynamic constitutive relations of soils based on Davidenkov skeleton curve［J］.Chinese Journal of Geotechnical Engineering，2005，27（8）：860-864.

［8］ 丁祖德，黃娟，袁鐵映，等.昆明泥炭質(zhì)土動剪切模量與阻尼比的試驗研究［J］.巖土力學，2017，38（12）：3627-3634.

DING Zude，HUANG Juan，YUAN Tieying，et al.Experimental study of dynamic shear modulus and damping ratio of peaty soil in Kunming［J］.Rock and Soil Mechanics，2017，38（12）：3627-3634.

［9］ MARTIN P P，SEED H B.One-dimensional dynamic ground response analyses［J］.Journal of the Geotechnical Engineering Division，1982，108（7）：935-952.

［10］ PYKE R M.Nonlinear soil models for irregular cyclicloadings［J］.Journal of the Geotechnical Engineering Division，1979，105（6）：715-726.

［11］ 趙丁鳳，阮濱，陳國興，等.基于Davidenkov骨架曲線模型的修正不規(guī)則加卸載準則與等效剪應變算法及其驗證［J］.巖土工程學報，2017，39（5）：888-895.

ZHAO Dingfeng，RUAN Bin，CHEN Guoxing，et al.Validation of modified irregular loading-unloading rules based on Davidenkov skeleton curve and its equivalent shear strain algorithm implemented in ABAQUS［J］.Chinese Journal of Geotechnical Engineering，2017，39（5）：888-895.

［12］ 王國波，尹驥，楊林德，等.Davidenkov模型在FLAC^3D中的開發(fā)及驗證［J］.武漢理工大學學報，2008，30（8）：143-146.

WANG Guobo，YIN Ji，YANG Linde，et al.Implementation and verification of Davidenkov model in FLAC^3D［J］.Journal of Wuhan University of Technology，2008，30（8）：143-146.

［13］ 趙旭清.Davidenkov土動力非線性模型在OpenSees軟件中的實現(xiàn)與應用［J］.土木工程與管理學報，2015，32（2）：19-22，73.

ZHAO Xuqing.Implementation and application of nonlinear dynamic Davidenkov model of soil in OpenSees［J］.Journal of Civil Engineering and Management，2015，32（2）：19-22，73.

［14］ 熊璐，賈通，陳君毅，等.基于有限狀態(tài)機的預期功能安全危害識別方法［J］.同濟大學學報（自然科學版），2023，51（4）：616-622.

XIONG Lu，JIA Tong，CHEN Junyi，et al.Hazard identification method for safety of the intended functionality based on finite state machine［J］.Journal of Tongji University （Natural Science），2023，51（4）：616-622.

［15］ 許君杰，郭麗峰，劉莉，等.仿人足球機器人視覺系統(tǒng)精確定位的特征提取［J］.機械設計與制造，2023（6）：258-262.

XU Junjie，GUO Lifeng，LIU Li，et al.Feature extraction for accurate localization of humanoid soccer robot vision system［J］.Machinery Design amp; Manufacture，2023（6）：258-262.

［16］ 張秀麗，王琪，黃森威，等.一種多模型融合的仿獵豹四足機器人復雜運動控制方法［J］.機器人，2022，44（6）：682-693，707.

ZHANG Xiuli，WANG Qi，HUANG Senwei，et al.A multi-model fusion based complex motion control approach for a cheetah-mimicking quadruped robot［J］.Robot，2022，44（6）：682-693，707.

［17］ 劉培文，舒輝，呂小少，等.基于有限狀態(tài)機的內(nèi)核漏洞攻擊自動化分析技術［J］.計算機科學，2022，49（11）：326-334.

LIU Peiwen，SHU Hui，L Xiaoshao，et al.Automatic analysis technology of kernel vulnerability attack based on finite state machine［J］.Computer Science，2022，49（11）：326-334.

［18］ SOLOV'EV V V.Synthesis of fast finite state machines on programmable logic integrated circuits by splitting internal states［J］.Journal of Computer and Systems Sciences International，2022，61（3）：360-371.

［19］ CONTRERAS C，QUIRS J C，CASAUCAO I，et al.STATCOM switching technique based on a finite-state machine［J］.Electronics，2023，12（6）：1481.

［20］ JIN D Y，ITO Y.Application offinite state machine theory to the simulation of reversed non-linear hysteretic relationships［J］.Tsinghua Science and Technology，2008，13（S1）：46-52.

［21］ HARDIN B O，DRNEVICH V P.Shear modulus and damping in soils：measurement and parameter effects （terzaghi leture）［J］.Journal of the Soil Mechanics and Foundations Division，1972，98（6）：603-624.

（本文編輯：任 棟）