1.一種基于裂隙與基質(zhì)網(wǎng)格的非達西滲流應力耦合仿真方法,其特征在于,包括:
2.根據(jù)權(quán)利要求1所述的一種基于裂隙與基質(zhì)網(wǎng)格的非達西滲流應力耦合仿真方法,其特征在于,所述輸入仿真基質(zhì)網(wǎng)格,根據(jù)裂隙的屬性信息建立仿真裂隙網(wǎng)格,獲取所述仿真裂隙網(wǎng)格與所述仿真基質(zhì)網(wǎng)格的網(wǎng)格單元相交參考信息,還包括:
3.根據(jù)權(quán)利要求1所述的一種基于裂隙與基質(zhì)網(wǎng)格的非達西滲流應力耦合仿真方法,其特征在于,所述依據(jù)所述仿真裂隙網(wǎng)格與所述仿真基質(zhì)網(wǎng)格選取網(wǎng)格節(jié)點的初始水壓力和初始位移值,還包括:
4.根據(jù)權(quán)利要求1所述的一種基于裂隙與基質(zhì)網(wǎng)格的非達西滲流應力耦合仿真方法,其特征在于,所述采用裂隙非達西滲流與基質(zhì)滲流時間步迭代求解,得到水壓力容差,判斷時間步迭代計算的收斂性,還包括:
5.根據(jù)權(quán)利要求4所述的一種基于裂隙與基質(zhì)網(wǎng)格的非達西滲流應力耦合仿真方法,其特征在于,所述依據(jù)所述收斂性與時間步信息,得到非達西滲流應力耦合的求解結(jié)果,還包括:
6.根據(jù)權(quán)利要求2所述的一種基于裂隙與基質(zhì)網(wǎng)格的非達西滲流應力耦合仿真方法,其特征在于,所述獲取所述基質(zhì)網(wǎng)格單元與所述裂隙網(wǎng)格單元的相交參考信息時,還包括:
7.根據(jù)權(quán)利要求3所述的一種基于裂隙與基質(zhì)網(wǎng)格的非達西滲流應力耦合仿真方法,其特征在于,依據(jù)所述仿真裂隙網(wǎng)格選取裂隙網(wǎng)格節(jié)點的初始裂隙水壓力以及依據(jù)所述仿真基質(zhì)網(wǎng)格選取基質(zhì)網(wǎng)格節(jié)點的初始基質(zhì)水壓力和初始基質(zhì)位移時,還包括:
8.根據(jù)權(quán)利要求4所述的一種基于裂隙與基質(zhì)網(wǎng)格的非達西滲流應力耦合仿真方法,其特征在于,依次求解基質(zhì)網(wǎng)格體系的滲流方程與求解裂隙網(wǎng)格體系的非達西滲流方程時,還包括:
9.根據(jù)權(quán)利要求8所述的一種基于裂隙與基質(zhì)網(wǎng)格的非達西滲流應力耦合仿真方法,其特征在于,采用位移容差判斷時間步迭代計算是否收斂,如果收斂,轉(zhuǎn)至下一時間步計算,否則轉(zhuǎn)至上一步繼續(xù)迭代計算時,還包括:
10.根據(jù)權(quán)利要求8所述的一種基于裂隙與基質(zhì)網(wǎng)格的非達西滲流應力耦合仿真方法,其特征在于,所述采用自適應時間步長,判斷是否需要重復上一時間步的計算,還包括: