Displacement functions over a one-dimensional bar corresponding to the predefined phase-field function with different element sizes and a kg = 10^−5

Different predefined phase-field function

\phi = exp(-sqrt(x^2)/l)

representing fully broken state in the process of a developing crack were applied to a one-dimensional bar with L =[−1, 1]. The regularising internal length was chosen to be l = 2. The mechanical part of the phase-field formulation was solved via FE-simulation for the displacement field u(x) across the crack, with the applied boundary conditions u(x =−1) =−1 and u(x = 1) = 1. Linear ansatz functions are used, but each time a different number of finite elements. The resulting displacement functions are depicted.

for each file -> Columns
1: x-coordinates of the 1D bar 
2: displacement result

file fig_03_data_1.csv: l/h =  5, 10 elements
file fig_03_data_2.csv: l/h = 15, 30 elements
file fig_03_data_3.csv: l/h = 25, 50 elements
