[Getfem-commits] r4855 - /trunk/getfem/interface/tests/matlab/demo_dynamic_plasticity_with_contac.m

[farshid . dabaghi]

Author: fdabaghi
Date: Mon Feb  9 18:12:27 2015
New Revision: 4855

URL: http://svn.gna.org/viewcvs/getfem?rev=4855&view=rev
Log:
modification

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

Modified: 
trunk/getfem/interface/tests/matlab/demo_dynamic_plasticity_with_contac.m
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/interface/tests/matlab/demo_dynamic_plasticity_with_contac.m?rev=4855&r1=4854&r2=4855&view=diff
==============================================================================
--- trunk/getfem/interface/tests/matlab/demo_dynamic_plasticity_with_contac.m   
(original)
+++ trunk/getfem/interface/tests/matlab/demo_dynamic_plasticity_with_contac.m   
Mon Feb  9 18:12:27 2015
@@ -45,13 +45,13 @@
 % alpha is parametr of the generalized integration algorithms The
 % The choice alpha = 1/2 yields the mid point method and alpha = 1 leads to
 % backward Euler integration
-alpha_method = false;
-alpha = 1.0;
+alpha_method = true;
+alpha = 1;
 
 
 f = [0 15000]';
 
-dirichlet_val = 5;
+dirichlet_val = 7;
 
 
 % transient part.
@@ -292,8 +292,8 @@
 % Iterations
 for t = 0:dt:T
   
-  %coeff = -sin(1-t);
-  coeff=- 1;
+  coeff = -sin(1-(4*t));
+ % coeff=- 1;
   
    disp(sprintf('step %d, coeff = %g', step , coeff)); 
    set(md, 'variable', 'VolumicData', get(mf_data, 
'eval',{f(1,1)*coeff;f(2,1)*coeff}));