這個(gè)算例可以看出離散元在模擬散體顆粒力學(xué)性質(zhì)方面的優(yōu)越性。因?yàn)槲易约捍T士大論文做的是有限元循環(huán)塑性本構(gòu)，用的是下負(fù)荷面模型實(shí)現(xiàn)的循環(huán)特性，主要的模型參數(shù)多達(dá)12個(gè)。而在離散元中，我實(shí)際上的材料參數(shù)只有emod、kratio、fric這三個(gè)，當(dāng)然如果孔隙率、級(jí)配也算的話(huà)就比較多了。下面講一下實(shí)現(xiàn)不排水循環(huán)雙軸的實(shí)現(xiàn)思路。

    首先，不排水采用的方法為體積不變法。

    dx=(-A/y^2)*dy

用這個(gè)可以根據(jù)豎向速度計(jì)算橫向速度，速度的方向如下圖。

如下圖為計(jì)算過(guò)程中面積的變化，有少許變化，但是這個(gè)量級(jí)可以看一下，很小。

如下圖為計(jì)算過(guò)程中的應(yīng)力路徑變化。

這里已經(jīng)算了一天了，基本上可以看出比較經(jīng)典的滯回圈和液化現(xiàn)象。當(dāng)p為0基本上可以認(rèn)為是液化。

應(yīng)力應(yīng)變曲線(xiàn)也是比較經(jīng)典的。這個(gè)可以自行和砂土實(shí)驗(yàn)對(duì)比一下。

這里記錄了墻上力的變化，可以根據(jù)正應(yīng)力的損失去計(jì)算靜水壓力。

注意：此處循環(huán)實(shí)現(xiàn)用的是勻速加載，監(jiān)測(cè)應(yīng)力到達(dá)一定值后反向。