This document refers to the data submitted in conjunction with the below paper:

Title: On the use of domain-based material point methods for problems involving large distortion  [Computer Methods in Applied Mechanics and Engineering (2019)]

Authors: L. Wang, W.M.Coombs, C.E.Augarde


Note: 'x' indicates values with respect to horizontal axis and 
      'y' indicates values with respect to vertical axis

-------------------------------------------------------------------------------------------------
fig3_sigxx.txt and fig3_sigzz.txt contain the data used in producing curves
about sigma_xx and sigma_zz in Figure 3, respectively. In each file,
Column A (x) is stress for analytical solution 
Column B (x) is stress for sMPM 
Column C (x) is stress for CPDI1 
Column D (x) is stress for CPDI2t 
Column E (y) is initial height 

-------------------------------------------------------------------------------------------------
fig5a.txt contains the data used in producing Figure 5(a):
Column A (x) is number of MPs per element 
Column B (y) is relative stress error for sMPM 
Column C (y) is relative stress error for CPDI1	 
Column D (y) is relative stress error for CPDI2t

fig5b.txt contains the data used in producing Figure 5(b):
Column A (x) is number of elements 
Column B (y) is relative stress error for sMPM 
Column C (y) is relative stress error for CPDI1	 
Column D (y) is relative stress error for CPDI2t

-------------------------------------------------------------------------------------------------
fig7a.txt contains the data used in producing Figure 7(a):
Column A (x) is the iteration 
Column B (y) is residual 

fig7b.txt contains the data used in producing Figure 7(b):
Column A (x) is the residual of previous iteration 
Column B (y) is residual 

Notably, Figure 7 uses the logrithm scale for axes. 

-------------------------------------------------------------------------------------------------
fig8.txt contains the data used in producing Figure 8:
Column A (x) is the dispalcement, u/L 
Column B (y) is reaction force for FEM 
Column C (y) is reaction force for sMPM

-------------------------------------------------------------------------------------------------
fig10 contains data file for Figure 10:
sMPM.txt for (a), CPDI1.txt for (b), CPDI2q.txt for (c) and CPDI2t.txt for
(d). Each file includes four columns:
Column A (x) is the rotation angle, alpha 
Column B (y) is relative error with rf=1
Column C (y) is relative error with rf=2
Column D (y) is relative error with rf=4
Column E (y) is relative error with rf=8


-------------------------------------------------------------------------------------------------
fig11.txt contains the data used in producing Figure 11:
Column A (x) is the mesh refinement level, rf 
Column B (y) is mean of relative error for sMPM
Column C (y) is mean of relative error for CPDI1
Column D (y) is mean of relative error for CPDI2q
Column E (y) is mean of relative error for CPDI2t
Notably, Figure 11 uses the logrithm scale for axes. 


-------------------------------------------------------------------------------------------------
fig12 contains data file for Figure 12:
rf1.txt for (a), rf2.txt for (b), rf4.txt for (c) and rf8.txt for
(d). Each file includes four columns:
Column A (x) is the rotation angle, alpha 
Column B (y) is relative error for sMPM
Column C (y) is relative error for CPDI1
Column D (y) is relative error for CPDI2q
Column E (y) is relative error for CPDI2t


-------------------------------------------------------------------------------------------------
fig13.txt contains the data used in producing Figure 13:
Column A (x) is the rotation angle, alpha 
Column B (y) is relative error in u for sMPM
Column C (y) is relative error in u for CPDI1 
Column D (y) is relative error in u for CPDI2q 
Column E (y) is relative error in u for CPDI2t


-------------------------------------------------------------------------------------------------
fig16a.txt contains the data used in producing Figure 16(a):
Column A (x) is the dispalcement 
Column B (y) is nodal reaction force for sMPM
Column C (y) is nodal reaction force for CPDI1 
Column D (y) is nodal reaction force for CPDI2q
Column E (y) is nodal reaction force for CPDI2t

fig16b.txt contains the data used in producing Figure 16(b):
Column A (x) is the dispalcement 
Column B (y) is nodal reaction force for sMPM
Column C (y) is nodal reaction force for CPDI1 
Column D (y) is nodal reaction force for CPDI2q
Column E (y) is nodal reaction force for CPDI2t


-------------------------------------------------------------------------------------------------
fig18b.txt contains the data used in producing Figure 18(b):
Column A (x) is the gravity loads 
Column B (y) is displacement for sMPM
Column C (y) is displacement for CPDI1 
Column D (y) is displacement for CPDI2q
Column E (y) is displacement for CPDI2t


-------------------------------------------------------------------------------------------------
fig20 contains data file for Figure 20:
fem_fine.txt include two columns for (x) and (y) in each subfigure for 'FEM(Fine)'
     fem.txt include two columns for (x) and (y) in each subfigure for 'FEM'
    sMPM.txt include two columns for (x) and (y) for 'sMPMM'
   CPDI1.txt include two columns for (x) and (y) for 'CPDI1'
  CPDI2q.txt include two columns for (x) and (y) for 'CPDI2q'
  CPDI2t.txt include two columns for (x) and (y) for 'CPDI2t'


-------------------------------------------------------------------------------------------------
fig21.txt contains the data used in producing Figure 21:
Column A (x) is the rotation angle, alpha 
Column B (y) is reaction force for FEM
Column C (y) is reaction force for sMPM
Column D (y) is reaction force for CPDI1 
Column E (y) is reaction force for CPDI2q
Column F (y) is reaction force for CPDI2t


-------------------------------------------------------------------------------------------------
fig25.txt contains the data used in producing Figure 25:
Column A (x) is the rotation angle, alpha 
Column B (y) is torque for sMPM
Column C (y) is torque for CPDI1 
Column D (y) is torque for CPDI2q
Column E (y) is torque for CPDI2t


