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

如题，想模拟下岩体实际变形，所以考虑改变它的本构模型，把 [1OX:O|  
            [ 1    u     0 f{s}[p~  
x
vx5@lx  
E/(1-u^2) *        u    1     0 "eqNd"~  
fk*$}f  
                         0    0     (1-u) /2  ] >_R,^iH"  

;,d^=:S6@  
|=R@nn   
变化为   [e11    e12     e13 woQ UrO(  
j
c%  
                  e21    e22    e23 J.nJ@?O+  
7x6q:4Ep\  
                  e31    e32    e33 ] PVK. %y9  
ccUI\!TD{/  
请大家指教应如何考虑，如何变化？ I_QWdxn  
谢谢

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

是石根华的DDA原程序啊，在他的博士论文里面，他提到： 
~7kIe+V  
if non-linear stress-strain relationships are considered,[E] has the following general shape: e9F\U   
1\hh,s  
[E]=( e11    e12    e13                   q
GAbh  
         e21    e22    e23 Xp_3EQl  
         e31    e32    e33 )