Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- do ii=1,nslip
- !
- k_alpha(ii)=n_k_alpha(ii)-XX(ii)
- phi(ii)=(sum( Mbar*s_x_m(1:9,ii) ) )-&
- ( -para(8)*k_alpha(ii)+kappa_nl(ii) )-para(11)
- if ( phi(ii).gt.0.d0 ) then
- mclaur=phi(ii)
- else
- mclaur=0.d0
- endif
- !Computation of R_\Delta_\lambda_\alpha
- R(ii)= para(12)*XX(ii)-dt*( mclaur )**para(3) ! Eq.(40) Diss Mikkel
- enddo
- .
- .
- .
- do ii=1,nslip
- do jj=1,nslip
- if (phi(ii).gt.0.d0) then
- if (ii.eq.jj) then
- Jakob(ii,jj)=Jakob(ii,jj)-para(8)
- Jakob(ii,jj)=para(12)-dt*para(3)*( mclaur )**(para(3)-1.d0)*Jakob(ii,jj)
- else
- !Jakob(ii,jj)=Jakob(ii,jj)*( XX(ii)/(1.d0-damage)*2.d0/lambda_0**2*(XX(jj)-n_k_alpha(jj)) )
- Jakob(ii,jj)=-dt*para(3)*( mclaur )**(para(3)-1.d0)*Jakob(ii,jj)
- endif
- else
- Jakob(ii,jj)=0.d0
- if (ii.eq.jj) then
- Jakob(ii,jj)=para(12)
- endif
- endif
- enddo
- enddo
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement