请教下，想把DDA里面的弹性变形假设修改成弹塑性的，如何改变它的弹性子矩阵 [复制链接]

如题，想模拟下岩体实际变形，所以考虑改变它的本构模型，把 '?1g_C QsS  
            [ 1    u     0 7M<co,"  
Rw^X5ByJE  
E/(1-u^2) *        u    1     0 v.F|8 cG  
>cH}sNHy  
                         0    0     (1-u) /2  ] @PV3G KJ  
{J-Ojw|Y b  
2E0oLl[  
变化为   [e11    e12     e13 \lF-]vz*  
8aTo TA7JA  
                  e21    e22    e23 rC6@ ]  
WNSEc%  
                  e31    e32    e33 ] ;N!opg))d<  
\S>GtlQbn  
请大家指教应如何考虑，如何变化？ |a=7P  
谢谢

一直搞不清楚搞非连续变形和数值流形都是自己在编译程序还是有现成的程序，现有的是什么程序。

是石根华的DDA原程序啊，在他的博士论文里面，他提到： ,r;E[k@  
if non-linear stress-strain relationships are considered,[E] has the following general shape: j9.%(*  
GN+!o($  
[E]=( e11    e12    e13                   |][PbN D  
         e21    e22    e23 dRUmC H  
         e31    e32    e33 )