[Getfem-commits] r4815 - /trunk/getfem/interface/tests/matlab/demo_plasticity.m

[farshid . dabaghi]

Author: fdabaghi
Date: Wed Nov 12 16:09:10 2014
New Revision: 4815

URL: http://svn.gna.org/viewcvs/getfem?rev=4815&view=rev
Log:
stress-strain fig

Modified:
    trunk/getfem/interface/tests/matlab/demo_plasticity.m

Modified: trunk/getfem/interface/tests/matlab/demo_plasticity.m
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/interface/tests/matlab/demo_plasticity.m?rev=4815&r1=4814&r2=4815&view=diff
==============================================================================
--- trunk/getfem/interface/tests/matlab/demo_plasticity.m       (original)
+++ trunk/getfem/interface/tests/matlab/demo_plasticity.m       Wed Nov 12 
16:09:10 2014
@@ -193,10 +193,18 @@
       % Compute Von Mises or Tresca stress
       VM = get(md, 'compute elastoplasticity Von Mises or Tresca', 'sigma', 
mf_vm, 'Von Mises');
     end
-  
+    
+    
+        expr_epsilon= '((Grad_u+Grad_u'')/2)';
+       Epsilon_u = gf_model_get(md, 'interpolation', expr_epsilon, mim_data); 
+       
+       sigma_fig(1,step)=sigma(90001);
+       Epsilon_u_fig(1,step)=Epsilon_u(90001);
+       
+       
     if (do_plot)
       figure(2)
-      subplot(2,1,1);
+      subplot(3,1,1);
       gf_plot(mf_vm,VM, 'deformation',U,'deformation_mf',mf_u,'refine', 4, 
'deformation_scale',1, 'disp_options', 0); % 'deformed_mesh', 'on')
       colorbar;
       axis([-20 120 -20 40]);
@@ -204,7 +212,7 @@
       n = t(step);
       title(['Von Mises criterion for t = ', num2str(step)]);
 
-      subplot(2,1,2);
+      subplot(3,1,2);
       gf_plot(mf_vm,plast, 'deformation',U,'deformation_mf',mf_u,'refine', 4, 
'deformation_scale',1, 'disp_options', 0);  % 'deformed_mesh', 'on')
       colorbar;
       axis([-20 120 -20 40]);
@@ -212,9 +220,17 @@
       n = t(step);
       title(['Plastification for t = ', num2str(step)]);
     
+      subplot(3,1,3);
+      
+     
+     plot(Epsilon_u_fig, sigma_fig,'r','LineWidth',2)
+        xlabel('Strain');
+        ylabel('Stress')
+        axis([0 0.35 -8000 16000 ]);
+    hold on;
       pause(0.1);
     end
-
+ 
 end;