diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h index 4799cfc6a3205e67c9bf4bfd38baee0456b27cb5..84cba972d677d52dbb5fa75f53df81865cb0c69d 100644 --- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h +++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h @@ -668,7 +668,8 @@ void RichardsMechanicsLocalAssembler<ShapeFunctionDisplacement, t, x_position, dt, temperature); auto const K_SR = - _ip_data[ip].solid_material.getBulkModulus(t, x_position, &C_el); + _ip_data[ip].solid_material.getBulkModulus(t, x_position, &C_el) / + (1 - alpha); auto const K_LR = liquid_phase.property(MPL::PropertyType::bulk_modulus)