// Copyright (C) Stichting Deltares 2018. 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 Deltares.DamEngine.Calculators.Uplift;
using Deltares.DamEngine.Data.General;
using Deltares.DamEngine.Data.General.PlLines;
using Deltares.DamEngine.Data.Geotechnics;
namespace Deltares.DamEngine.Calculators.KernelWrappers.Common
{
///
/// Helper for Uplift
///
public class UpliftHelper
{
public static UpliftSituation DetermineStabilityUplift(DamKernelInput damKernelInput, bool useRivelLevelLow)
{
var location = damKernelInput.Location;
var upliftCriterion = damKernelInput.DesignScenario.GetUpliftCriterionStability(location.ModelFactors.UpliftCriterionStability);
var riverLevel = damKernelInput.DesignScenario.RiverLevel;
double? riverLevelLow = null;
if (useRivelLevelLow)
{
riverLevelLow = damKernelInput.DesignScenario.RiverLevelLow;
}
var soilGeometry2DName = damKernelInput.SubSoilScenario.FullStiFileName;
UpliftSituation upliftSituation;
var plLines = PlLinesHelper.CreatePlLinesForStability(location, damKernelInput.SubSoilScenario, riverLevel,
soilGeometry2DName, riverLevelLow, out upliftSituation);
if (plLines != null)
{
var surfaceLineWithNewHeight = location.SurfaceLine;
double? upliftFactor = GetLowestUpliftFactor(surfaceLineWithNewHeight,
damKernelInput.SubSoilScenario.SoilProfile1D, soilGeometry2DName, plLines, location);
upliftSituation.IsUplift = (upliftFactor < upliftCriterion);
}
else
{
upliftSituation.IsUplift = false;
}
return upliftSituation;
}
public static double? GetLowestUpliftFactor(SurfaceLine2 surfaceLine, SoilProfile1D soilProfile, string soilGeometry2DName,
PLLines plLines, Location location)
{
var upliftLocationDeterminator = new UpliftLocationDeterminator()
{
SurfaceLine = surfaceLine,
SoilProfile = soilProfile,
SoilGeometry2DName = soilGeometry2DName,
SoilList = location.SoilList,
DikeEmbankmentMaterial = location.GetDikeEmbankmentSoil(),
PLLines = plLines,
XSoilGeometry2DOrigin = location.XSoilGeometry2DOrigin
};
UpliftLocationAndResult upliftLocationAndResult = upliftLocationDeterminator.GetLocationAtWithLowestUpliftFactor();
return upliftLocationAndResult?.UpliftFactor;
}
}
}