public static void CalculateEffectiveThicknessCoverageLayer_CompleteInput_InputSetOnSubCalculator()
{
// Setup
var validPipingCalculation =
SemiProbabilisticPipingCalculationTestFactory.CreateCalculationWithValidInput <TestSemiProbabilisticPipingCalculation>(
new TestHydraulicBoundaryLocation());
PipingInput input = validPipingCalculation.InputParameters;
var generalInput = new GeneralPipingInput();
using (new PipingSubCalculatorFactoryConfig())
{
// Call
PipingInput inputParameters = validPipingCalculation.InputParameters;
InputParameterCalculationService.CalculateEffectiveThicknessCoverageLayer(
generalInput.WaterVolumetricWeight,
PipingDesignVariableFactory.GetPhreaticLevelExit(inputParameters).GetDesignValue(),
inputParameters.ExitPointL,
inputParameters.SurfaceLine,
inputParameters.StochasticSoilProfile.SoilProfile);
// Assert
var testFactory = (TestPipingSubCalculatorFactory)PipingSubCalculatorFactory.Instance;
EffectiveThicknessCalculatorStub effectiveThicknessCalculator = testFactory.LastCreatedEffectiveThicknessCalculator;
Assert.AreEqual(input.ExitPointL.Value, effectiveThicknessCalculator.ExitPointXCoordinate);
Assert.AreEqual(PipingDesignVariableFactory.GetPhreaticLevelExit(input).GetDesignValue(),
effectiveThicknessCalculator.PhreaticLevel,
input.PhreaticLevelExit.GetAccuracy());
AssertEqualSoilProfiles(input.StochasticSoilProfile.SoilProfile, effectiveThicknessCalculator.SoilProfile);
AssertEqualSurfaceLines(input.SurfaceLine, effectiveThicknessCalculator.SurfaceLine);
Assert.AreEqual(generalInput.WaterVolumetricWeight, effectiveThicknessCalculator.VolumicWeightOfWater);
}
}
public static void CalculateEffectiveThicknessCoverageLayer_WithMultipleCharacteristicTypesOnSamePoint_ReturnsThickness()
{
// Setup
var surfaceLine = new PipingSurfaceLine(string.Empty);
surfaceLine.SetGeometry(new[]
{
new Point3D(0, 0, 10),
new Point3D(20, 0, 10)
});
surfaceLine.SetDikeToeAtRiverAt(surfaceLine.Points.ElementAt(0));
surfaceLine.SetDikeToeAtPolderAt(surfaceLine.Points.ElementAt(0));
surfaceLine.SetDitchDikeSideAt(surfaceLine.Points.ElementAt(0));
surfaceLine.SetBottomDitchPolderSideAt(surfaceLine.Points.ElementAt(1));
surfaceLine.SetBottomDitchDikeSideAt(surfaceLine.Points.ElementAt(1));
surfaceLine.SetDitchPolderSideAt(surfaceLine.Points.ElementAt(1));
var stochasticSoilProfile = new PipingStochasticSoilProfile(
0.0, new PipingSoilProfile(string.Empty, 0, new[]
{
new PipingSoilLayer(5)
{
IsAquifer = true
},
new PipingSoilLayer(20)
{
IsAquifer = false
}
}, SoilProfileType.SoilProfile1D));
var input = new TestPipingInput
{
ExitPointL = (RoundedDouble)10,
SurfaceLine = surfaceLine,
StochasticSoilProfile = stochasticSoilProfile
};
var generalInput = new GeneralPipingInput();
// Call
double thickness = InputParameterCalculationService.CalculateEffectiveThicknessCoverageLayer(
generalInput.WaterVolumetricWeight,
PipingDesignVariableFactory.GetPhreaticLevelExit(input).GetDesignValue(),
input.ExitPointL,
input.SurfaceLine,
input.StochasticSoilProfile.SoilProfile);
// Assert
Assert.AreEqual(5, thickness);
}
private static void UpdateEffectiveThicknessCoverageLayerMean(PipingInput input, GeneralPipingInput generalInput,
LogNormalDistribution effectiveThicknessCoverageLayerDistribution)
{
if (input.SurfaceLine != null && input.StochasticSoilProfile?.SoilProfile != null && !double.IsNaN(input.ExitPointL))
{
var weightedMean = new RoundedDouble(GetNumberOfDecimals(effectiveThicknessCoverageLayerDistribution),
InputParameterCalculationService.CalculateEffectiveThicknessCoverageLayer(
generalInput.WaterVolumetricWeight,
PipingDesignVariableFactory.GetPhreaticLevelExit(input).GetDesignValue(),
input.ExitPointL,
input.SurfaceLine,
input.StochasticSoilProfile.SoilProfile));
if (weightedMean > 0)
{
effectiveThicknessCoverageLayerDistribution.Mean = weightedMean;
}
}
}
public static void CalculateEffectiveThicknessCoverageLayer_InvalidPipingCalculationWithOutput_ReturnsNaN()
{
// Setup
var invalidPipingCalculation =
SemiProbabilisticPipingCalculationTestFactory.CreateCalculationWithValidInput <TestSemiProbabilisticPipingCalculation>(
new TestHydraulicBoundaryLocation());
// Make invalid by having surface line partially above soil profile:
double highestLevelSurfaceLine = invalidPipingCalculation.InputParameters.SurfaceLine.Points.Max(p => p.Z);
double soilProfileTop = highestLevelSurfaceLine - 0.5;
double soilProfileBottom = soilProfileTop - 0.5;
invalidPipingCalculation.InputParameters.StochasticSoilProfile = new PipingStochasticSoilProfile(
0.0,
new PipingSoilProfile("A", soilProfileBottom, new[]
{
new PipingSoilLayer(soilProfileTop)
{
IsAquifer = true
}
}, SoilProfileType.SoilProfile1D));
var generalInput = new GeneralPipingInput();
// Call
PipingInput input = invalidPipingCalculation.InputParameters;
double result = InputParameterCalculationService.CalculateEffectiveThicknessCoverageLayer(
generalInput.WaterVolumetricWeight,
PipingDesignVariableFactory.GetPhreaticLevelExit(input).GetDesignValue(),
input.ExitPointL,
input.SurfaceLine,
input.StochasticSoilProfile.SoilProfile);
// Assert
Assert.IsNaN(result);
}
