// Copyright (C) Stichting Deltares 2017. All rights reserved.
//
// This file is part of the Dam Engine.
//
// The Dam Engine is free software: you can redistribute it and/or modify
// it under the terms of the GNU Affero General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU Affero General Public License for more details.
//
// You should have received a copy of the GNU Affero General Public License
// along with this program. If not, see .
//
// All names, logos, and references to "Deltares" are registered trademarks of
// Stichting Deltares and remain full property of Stichting Deltares at all times.
// All rights reserved.
using System;
using System.Collections.Generic;
using System.Linq;
using Deltares.DamEngine.Data.General;
using Deltares.DamEngine.Data.RWScenarios;
namespace Deltares.DamEngine.Calculators.Dikes_Assessment_Regional
{
public class RWScenarioSelector
{
public RWScenarioSelector()
{
}
public PipingModelType PipingModelType { get; set; }
public MStabParameters MStabParameters { get; set; }
private bool WillChoicesResultInScenarios(params Enum[] previousChoices)
{
if (previousChoices.Length >= 2)
{
if (previousChoices[0].Equals(LoadSituation.Dry) && previousChoices[1].Equals(DikeDrySensitivity.None))
{
return false;
}
}
return true;
}
public RWScenarioProfileResult[] GetScenarios(Location location, SoilGeometryProbability soilGeometryProbability)
{
List scenarios = new List();
RWEvaluator[] evaluators = new RWEvaluator[]
{
new DikeMaterialRWEvaluator(),
new HydraulicShortcutRWEvaluator(),
new UpliftRWEvaluator()
};
Enum[] choices = new Enum[evaluators.Length + 1];
foreach (LoadSituation loadSituation in Enum.GetValues(typeof(LoadSituation)))
{
choices[0] = loadSituation;
bool resultsInScenarios = true;
for (int i = 1; i < evaluators.Length + 1; i++)
{
choices[i] = evaluators[i - 1].Evaluate(location, new SoilGeometry(soilGeometryProbability.SoilProfile, null), choices.ToArray());
//Is het een doorlopende tak? Zo niet break af (afbreken als geen veenkade en droge belastingsituatie
resultsInScenarios = WillChoicesResultInScenarios(choices.ToArray());
if (!resultsInScenarios)
{
break;
}
}
if (resultsInScenarios)
{
if (choices[0].Equals(LoadSituation.Wet))
{
if (choices[1].Equals(DikeDrySensitivity.None))
{
if (choices[2].Equals(HydraulicShortcutType.NoHydraulicShortcut))
{
if (choices[3].Equals(UpliftType.NoUplift))
{
scenarios.Add(this.GetScenario(ScenarioType.Scenario08, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.Bishop));
}
else // UpliftType.Uplift
{
scenarios.Add(this.GetScenario(ScenarioType.Scenario04, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.UpliftVan));
scenarios.Add(this.GetScenario(ScenarioType.Scenario08, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.Bishop));
scenarios.Add(this.GetScenario(ScenarioType.Scenario10, choices, location, soilGeometryProbability, FailureMechanismSystemType.Piping, PipingModelType));
}
}
else // HydraulicShortcutType.HydraulicShortcut
{
if (choices[3].Equals(UpliftType.NoUplift))
{
scenarios.Add(this.GetScenario(ScenarioType.Scenario07, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.Bishop));
}
else // UpliftType.Uplift
{
scenarios.Add(this.GetScenario(ScenarioType.Scenario03, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.UpliftVan));
scenarios.Add(this.GetScenario(ScenarioType.Scenario07, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.Bishop));
scenarios.Add(this.GetScenario(ScenarioType.Scenario10, choices, location, soilGeometryProbability, FailureMechanismSystemType.Piping, PipingModelType));
}
}
}
else // DikeDrySensitivity.Dry
{
if (choices[2].Equals(HydraulicShortcutType.NoHydraulicShortcut))
{
if (choices[3].Equals(UpliftType.NoUplift))
{
scenarios.Add(this.GetScenario(ScenarioType.Scenario08, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.Bishop));
}
else // UpliftType.Uplift
{
scenarios.Add(this.GetScenario(ScenarioType.Scenario04, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.UpliftVan));
scenarios.Add(this.GetScenario(ScenarioType.Scenario08, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.Bishop));
scenarios.Add(this.GetScenario(ScenarioType.Scenario10, choices, location, soilGeometryProbability, FailureMechanismSystemType.Piping, PipingModelType));
}
}
else // HydraulicShortcutType.HydraulicShortcut
{
if (choices[3].Equals(UpliftType.NoUplift))
{
scenarios.Add(this.GetScenario(ScenarioType.Scenario07, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.Bishop));
}
else // UpliftType.Uplift
{
scenarios.Add(this.GetScenario(ScenarioType.Scenario03, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.UpliftVan));
scenarios.Add(this.GetScenario(ScenarioType.Scenario07, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.Bishop));
scenarios.Add(this.GetScenario(ScenarioType.Scenario10, choices, location, soilGeometryProbability, FailureMechanismSystemType.Piping, PipingModelType));
}
}
}
}
else // LoadSituation.Dry
{
if (choices[1].Equals(DikeDrySensitivity.Dry))
{
if (choices[2].Equals(HydraulicShortcutType.NoHydraulicShortcut))
{
if (choices[3].Equals(UpliftType.NoUplift))
{
scenarios.Add(this.GetScenario(ScenarioType.Scenario06, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.Bishop));
}
else // UpliftType.Uplift
{
scenarios.Add(this.GetScenario(ScenarioType.Scenario02, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.UpliftVan));
scenarios.Add(this.GetScenario(ScenarioType.Scenario06, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.Bishop));
scenarios.Add(this.GetScenario(ScenarioType.Scenario09, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.HorizontalBalance));
scenarios.Add(this.GetScenario(ScenarioType.Scenario11, choices, location, soilGeometryProbability, FailureMechanismSystemType.Piping, PipingModelType));
}
}
else // HydraulicShortcutType.HydraulicShortcut
{
if (choices[3].Equals(UpliftType.NoUplift))
{
scenarios.Add(this.GetScenario(ScenarioType.Scenario05, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.Bishop));
}
else // UpliftType.Uplift
{
scenarios.Add(this.GetScenario(ScenarioType.Scenario01, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.UpliftVan));
scenarios.Add(this.GetScenario(ScenarioType.Scenario05, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.Bishop));
scenarios.Add(this.GetScenario(ScenarioType.Scenario09, choices, location, soilGeometryProbability, FailureMechanismSystemType.StabilityInside, MStabModelType.HorizontalBalance));
scenarios.Add(this.GetScenario(ScenarioType.Scenario11, choices, location, soilGeometryProbability, FailureMechanismSystemType.Piping, PipingModelType));
}
}
}
}
}
}
return scenarios.ToArray();
}
private RWScenarioProfileResult GetScenario(ScenarioType scenarioType, Enum[] choices, Location location,
SoilGeometryProbability soilGeometryProbability, FailureMechanismSystemType failureMechanismType, MStabModelType modelOption)
{
RWScenarioProfileResult scenario = new RWScenarioProfileResult(location, soilGeometryProbability);
scenario.FailureMechanismType = failureMechanismType;
scenario.MstabModelOption = modelOption;
scenario.PipingModelOption = PipingModelType;
scenario.LoadSituation = (LoadSituation)choices[0];
scenario.DikeDrySensitivity = (DikeDrySensitivity)choices[1];
scenario.HydraulicShortcutType = (HydraulicShortcutType)choices[2];
scenario.UpliftType = (UpliftType)choices[3];
scenario.ScenarioType = scenarioType;
//if (scenario.IsDry)
//{
// // Soils are not transferred to the computational core at the moment, therefore this code has no meaning at the moment
// // You could say the same about soils in general in the DAM C# code
// SoilProfile drySoilProfile = soilProfile.Clone();
// foreach (Layer layer in drySoilProfile.Layers)
// {
// layer.Soil.AbovePhreaticLevel = layer.Soil.DryUnitWeight;
// }
//}
return scenario;
}
private RWScenarioProfileResult GetScenario(ScenarioType scenarioType, Enum[] choices, Location location,
SoilGeometryProbability soilGeometryProbability, FailureMechanismSystemType failureMechanismType, PipingModelType pipingModelType)
{
RWScenarioProfileResult scenario = new RWScenarioProfileResult(location, soilGeometryProbability);
scenario.FailureMechanismType = failureMechanismType;
scenario.MstabModelOption = MStabModelType.Bishop;
scenario.PipingModelOption = pipingModelType;
scenario.LoadSituation = (LoadSituation)choices[0];
scenario.DikeDrySensitivity = (DikeDrySensitivity)choices[1];
scenario.HydraulicShortcutType = (HydraulicShortcutType)choices[2];
scenario.UpliftType = (UpliftType)choices[3];
scenario.ScenarioType = scenarioType;
//if (scenario.IsDry)
//{
// // Soils are not transferred to the computational core at the moment, therefore this code has no meaning at the moment
// // You could say the same about soils in general in the DAM C# code
// SoilProfile drySoilProfile = soilProfile.Clone();
// foreach (Layer layer in drySoilProfile.Layers)
// {
// layer.Soil.AbovePhreaticLevel = layer.Soil.DryUnitWeight;
// }
//}
return scenario;
}
}
}