请教j2理论最后更新jacobian 矩阵时候的小问题

magiconly_2010

http://forum.simwe.com/thread-914547-2-1.html Divid 的帖子给出了一个基本例题的注视；程序中这一段有点疑问
C JACOBIAN-------------------[计算雅可比矩阵]+ t+ ?/ m' x, }3 S8 d2 n' l) g" I
C. k# b6 d* R: Q+ }' m
EFFG=EG*SYIELD/SMISES' c; R& T% F% M% |: n; `, H7 t
EFFG2=TWO*EFFG2 G, ^& o; X; \0 n' A5 s+ F
EFFG3=THREE*EFFG2/TWO
EFFLAM=(EBULK3-EFFG2)/THREE, Y# B6 g: x/ n4 ?& f+ S
DO 220 K1=1,NDI, F6 f3 \. \; @. l3 |
DO 210 K2=1,NDI+ |3 H; \7 [" x: A) t2 S
DDSDDE(K2,K1)=EFFLAM
210 CONTINUE
DDSDDE(K1,K1)=EFFG2+EFFLAM( I8 o* r0 @. }5 ]
220 CONTINUE
DO 230 K1=NDI+1,NTENS6 g, k4 `5 l' }; H  j% ~4 j
DDSDDE(K1,K1)=EFFG. f2 @: I) Q+ m0 n0 J8 [
230 CONTINUE9 M9 a" B9 r3 @4 g: g7 X! |# M. L
DO 250 K1=1,NTENS
DO 240 K2=1,NTENS
DDSDDE(K2,K1)=DDSDDE(K2,K1)+FLOW(K2)*FLOW(K1)1 l9 r  a: [: E' u/ Z
1*(EFFHRD-EFFG3)2 o9 H$ C% d$ `' Q# k
240 CONTINUE: d/ h/ ~2 v+ Z& f* K) ?
250 CONTINUE
ENDIF
ENDIF! a5 r- |8 ], Q, S3 e
C这段计算切线刚度矩阵的程序里，这几句是怎么来的，请过来人指点一下
EFFG=EG*SYIELD/SMISES' c; R& T% F% M% |: n; `, H7 t
EFFG2=TWO*EFFG2 G, ^& o; X; \0 n' A5 s+ F
EFFG3=THREE*EFFG2/TWO
EFFLAM=(EBULK3-EFFG2)/THREE
DDSDDE(K2,K1)=DDSDDE(K2,K1)+FLOW(K2)*FLOW(K1)1 l9 r  a: [: E' u/ Z
1*(EFFHRD-EFFG3)2 o
这几句是怎么来的 , Y#

magiconly_2010

显示的有问题，重新发一下
CJACOBIAN-------------------[计算雅可比矩阵]+ t+ ?/ m' x, }3 S8 d2 n' l) g" I
C. k# b6 d* R: Q+ }' m
EFFG=EG*SYIELD/SMISES' c; R& T% F% M% |: n; `, H7 t
EFFG2=TWO*EFFG2 G, ^& o; X; \0 n' A5 s+ F
EFFG3=THREE*EFFG2/TWO
EFFLAM=(EBULK3-EFFG2)/THREE, Y# B6 g: x/ n4 ?& f+ S
DO 220 K1=1,NDI, F6 f3 \. \; @. l3 |
DO 210 K2=1,NDI+ |3 H; \7 [" x: A) t2 S
DDSDDE(K2,K1)=EFFLAM
210 CONTINUE
DDSDDE(K1,K1)=EFFG2+EFFLAM( I8 o* r0 @. }5 ]
220 CONTINUE
DO 230 K1=NDI+1,NTENS6 g, k4 `5 l' }; H  j% ~4 j
DDSDDE(K1,K1)=EFFG. f2 @: I) Q+ m0 n0 J8 [
230 CONTINUE9 M9 a" B9 r3 @4 g: g7 X! |# M. L
DO 250 K1=1,NTENS
DO 240 K2=1,NTENS
DDSDDE(K2,K1)=DDSDDE(K2,K1)+FLOW(K2)*FLOW(K1)1 l9 r  a: [: E' u/ Z
1*(EFFHRD-EFFG3)2 o9 H$ C% d$ `' Q# k
240 CONTINUE: d/ h/ ~2 v+ Z& f* K) ?
250 CONTINUE
ENDIF
ENDIF! a5 r- |8 ], Q, S3 e
C
这这段计算切线刚度矩阵的程序中，下面这几句是怎么来的
EFFG=EG*SYIELD/SMISES'c; R& T% F% M% |: n; `, H7 t
EFFG2=TWO*EFFG2 G, ^& o; X; \0 n' A5 s+ F
EFFG3=THREE*EFFG2/TWO
EFFLAM=(EBULK3-EFFG2)/THREE, Y# B6 g: x/ n4 ?& f+ S
DDSDDE(K2,K1)=DDSDDE(K2,K1)+FLOW(K2)*FLOW(K1)1l9 r  a: [: E' u/ Z

herrliubs

坛子里面有推导
http://forum.simwe.com/viewthrea ... ;highlight=jacobian