實現(xiàn)步驟：

1. 促進(jìn)學(xué)習(xí)，先放地址：https://github.com/pzp2/d_cmodel.git

2. 首先抄陳育民書第十三章的自定義本構(gòu)的框架

3. 然后看鄧肯張理論

4. 最后在vs2019里面修改生成dll文件

5. 導(dǎo)入flac3d7.0里面查看彈性模量，泊松比或者體積模量和剪切模量的變化情況

6. 這里通過flac自帶的python直接導(dǎo)出所有參數(shù)

7. 比對理論數(shù)值

8. 下面是理論部分以及我主要修改的部分，有些久沒看這些東西了。

? ? ? ? 最重要的修改是陳育民書中的卸載模量表示的中間參數(shù)，原本此參數(shù)在使用的時候是取百分比數(shù)值（判斷卸載與否的時候），但是使用flac3d7.0死活弄不成卸載的樣子，因此我才用達(dá)到最大彈性模量的75%是就判定發(fā)生卸載了，此時在flac3d里面運行得到的卸載模量正確。其次有一個大問題沒有得到解決：陳育民書中單元測試的應(yīng)變得不到書中的樣子，書中是采用位移除以總長度得到的應(yīng)變，這樣得出來的應(yīng)力應(yīng)變曲線斜率是逐漸增大的，這明顯違背理論（書中曲線），采用flac3d自帶的應(yīng)變zone.strain.inc.zz得到的應(yīng)變曲線（應(yīng)力應(yīng)變曲線）是逐漸穩(wěn)定的，既斜率趨于0，這是正常的。

? ? ? ? ? 不想寫了，要打王者榮耀了，還有些注意點是：flac命令流的邊界控制條件，大變形開關(guān)與否，三軸試驗的加載情況。。。

? ? ? ? ? 相關(guān)命令流全部在github里。里面vs2019里要是不想編譯，可以直接用我編譯的dll文件，然后放到.....flac3d_7.0\exe64\plugins\cmodel里面就可以直接運行flac3d7.0里面的p_260_Duncan-Chang.dat文件了，它會自動運行show_prin.py以便直接在flac3d里面顯示需要查看的信息（我原先用matlab查看，發(fā)現(xiàn)太傻了，因為管理浪費時間，能用自帶的就不要用其他的！?。。?/p>

運行好如下：

102670262.18831182——理論加載彈性模量

168265303.51945418——理論卸載彈性模量

102837761.62219529——flac彈性模量

_____________

42286033.691056535——理論體積模量

42317008.51123584——flac體積模量

_____________

46867092.862573996——理論剪切模量

46959137.13429132——flac剪切模量

_____________

-1360040.9212515894——z方向應(yīng)力

-599997.9660351506——x方向應(yīng)力

-599995.9760701974——y方向應(yīng)力

#############——循環(huán)，下面表示的是正確的卸載時的彈性模量

103588186.0152759

168265304.12230998

168525615.53291294

_____________

42286033.76282022

56175205.177637644

_____________

47442813.003809676

84262807.76645648

_____________

-1339939.9078893024

-599997.9687611009

-599995.9817271639

#############

順便把理論貼一貼：