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

如题，想模拟下岩体实际变形，所以考虑改变它的本构模型，把 :EYu 4Y  
            [ 1    u     0 ,I/2.Q})[  
!-F^VGD(8  
E/(1-u^2) *        u    1     0
?rky6  
Nvi Fq  
                         0    0     (1-u) /2  ] U.HoFf+HN  
D!c1;IHZ  
u,:CJ[3  
变化为   [e11    e12     e13 <'yf|N
!9G  
hEVjeC  
                  e21    e22    e23 B{[f}h.n  
=U=e?AOG2  
                  e31    e32    e33 ] E"G:K`Q  
@X3 gBGY)  
请大家指教应如何考虑，如何变化？ ROPC |  
谢谢

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

是石根华的DDA原程序啊，在他的博士论文里面，他提到： sYAG,r>h  
if non-linear stress-strain relationships are considered,[E] has the following general shape: hW0,5>[7%  
gt2>nTJz.Z  
[E]=( e11    e12    e13                   ]ro1{wm!WU  
         e21    e22    e23 1A|x$j6m  
         e31    e32    e33 )