之前做過一個真三軸的數(shù)值模擬，這里講一下常規(guī)三軸的基本思路。

1、成樣

    生成圓柱形的試樣，這個比較簡單，直接公開代碼

new
def chicun_par
    sample_rad=0.5
    sample_hight=2
   
    keli_rdmin=0.06
    keli_rdmax=0.09
end
@chicun_par

domain extent [-sample_hight*1.5] [sample_hight*1.5][n=1.4]
wall generate cylinder base 0 0 [-sample_hight*0.5*n] axis 0 0 1 ...
                    height [sample_hight*n] radius [sample_rad] cap false falsewall generate plane position 0 0 [sample_hight*0.5] dip 0 ddir 0
wall generate plane position 0 0 [-sample_hight*0.5] dip 0 ddir 0

ball distribute radius [keli_rdmin] [keli_rdmax] porosity 0.28 ...
    range cylinder end1 0 0 [sample_hight*0.5-keli_rdmin] ...
    end2 0 0 [-sample_hight*0.5+keli_rdmin] radius [sample_rad-keli_rdmin]

cmat default model linear method deform emod 100e6 kratio 1.5 property fric 0.5

ball attribute density 2.7e3 damp 0.7
cycle 2000 calm 50

solve

save sample

2、預壓

    這里是主要難點，圓柱形墻的伺服比較難，如何監(jiān)測側面墻的應力以及如何進行側面墻的加載是個難點。

這里的代碼在付費中公開。下圖為伺服的應力變化。

3、圍壓以及加載，這個和之前三軸一樣的。

這里因為顆粒數(shù)只有五百個，計算效果不是很好，這種單元試驗至少1w個才好看點。

代碼里面的加載速率自己調呀，我給的比較大。