In [1]:

    

from PD1D import PD_Problem
%matplotlib inline
import matplotlib.pyplot as plt


    

In [2]:

    

fixed_length = 40
delta_x = 0.5
fixed_horizon = 3 * delta_x

problem1 = PD_Problem(bar_length=fixed_length, number_of_elements=(fixed_length/delta_x), 
                      horizon=fixed_horizon, constitutive_model_flag='LPS', randomization_factor=0.3)
problem1.solve(prescribed_displacement=0.1)


    

In [3]:

    

nodes1 = problem1.get_nodes()
def_grad1 = problem1.get_deformation_gradient()
plt.plot(nodes1, def_grad1, 'k-');
plt.grid(b=True);
plt.axis();


    

In [4]:

    

fixed_horizon = 4 * delta_x

problem2 = PD_Problem(bar_length=fixed_length, number_of_elements=(fixed_length/delta_x), 
                      horizon=fixed_horizon, constitutive_model_flag='LPS', randomization_factor=0.3)
problem2.solve(prescribed_displacement=0.1)
nodes2 = problem2.get_nodes()
def_grad2 = problem2.get_deformation_gradient()
plt.plot(nodes2, def_grad2, 'k-');
plt.grid(b=True);
plt.axis();


    

In [5]:

    

fixed_horizon = 5 * delta_x

problem3 = PD_Problem(bar_length=fixed_length, number_of_elements=(fixed_length/delta_x), 
                      horizon=fixed_horizon, constitutive_model_flag='LPS', randomization_factor=0.3)
problem3.solve(prescribed_displacement=0.1)
nodes3 = problem3.get_nodes()
def_grad3 = problem3.get_deformation_gradient()
plt.plot(nodes3, def_grad3, 'k-');
plt.grid(b=True);
plt.axis();


    

In [6]:

    

fixed_horizon = 6 * delta_x

problem4 = PD_Problem(bar_length=fixed_length, number_of_elements=(fixed_length/delta_x), 
                      horizon=fixed_horizon, constitutive_model_flag='LPS', randomization_factor=0.3)
problem4.solve(prescribed_displacement=0.1)
nodes4 = problem4.get_nodes()
def_grad4 = problem4.get_deformation_gradient()
plt.plot(nodes4, def_grad4, 'k-');
plt.grid(b=True);
plt.axis();


    

In [7]:

    

plt.plot(nodes1, def_grad1, 'k-', label='$3 \\Delta x$')
plt.plot(nodes2, def_grad2, 'r-', label='$4 \\Delta x$')
plt.plot(nodes3, def_grad3, 'b-', label='$5 \\Delta x$')
plt.plot(nodes4, def_grad4, 'g-', label='$6 \\Delta x$')
plt.legend(loc='upper left')
plt.grid(b=True);
plt.axis('tight');
plt.show()