material property for each element

[sskalati]

Dear MOOSE Community,
I am trying to run a simulation using ComputeLinearElasticPhaseFieldFractureStress with a custom degradation function defined as:
expression = '((1.0-c)^2)/((1.0-c)^2+a1c(1-0.5*c))'
where a1 is an element-wise property that varies for each element. I provide these values as input to MOOSE, expecting the software to compute the degradation function dynamically based on the elemental data. However, it does not seem to work as intended.

Could you please guide me on how to properly implement this functionality? Are there any specific settings or methods required to ensure that MOOSE correctly utilizes the element-wise values of a1 in the degradation function?

I have attached the a1 values and the MOOSE script related to a1 values as well.
Any help or suggestions would be greatly appreciated!

[UserObjects]
    [a1_read]
        type = PropertyReadFile
        prop_file_name = 'a1_values.csv'
        nprop = 1
        read_type = element
    []
[]


[Functions]
    [a1_fun]
        type = PiecewiseConstantFromCSV
        read_prop_user_object = 'a1_read'
        read_type = 'element'
        column_number = '0'
    []
[]

[Materials]
    [./a1]
        type = GenericFunctionMaterial
        prop_names = a1
        prop_values = a1_fun
    [../]    

    [./degradation]
        type = DerivativeParsedMaterial
        property_name = degradation
        coupled_variables = 'c'
        material_property_names = 'a1'
        expression = '((1.0-c)^2)/((1.0-c)^2+a1*c*(1-0.5*c))'
        derivative_order = 2
    [../]
[]

a1_values.csv

[GiudGiud]

Hello

can you paste your Mesh block? To save me a minute.
This setup looks fine, we just need to check that a1 loads as expected

[sskalati]

[Mesh]
[gmsh]
type = GeneratedMeshGenerator
dim = 2
xmax = 4.1e-3
ymax = 4.1e-3
nx = 126
ny = 126
elem_type = QUAD4
[]

[xaxis]
    type = BoundingBoxNodeSetGenerator
    new_boundary = xaxis
    bottom_left = '0 2.05e-3 0'
    top_right = '4.1e-3 2.05e-3 0'
    input = gmsh
[]

[yaxis]
    type = BoundingBoxNodeSetGenerator
    new_boundary = yaxis
    bottom_left = '2.05e-3 0 0'
    top_right = '2.05e-3 4.1e-3 0'
    input = xaxis
[]

[]

[GiudGiud]

does a1 seem properly loaded?

with c = x + y
I get this for degradation. it does seem to work?