Index: dam engine/trunk/src/Deltares.DamEngine.Calculators.Tests/Deltares.DamEngine.Calculators.Tests.csproj =================================================================== diff -u -r680 -r681 --- dam engine/trunk/src/Deltares.DamEngine.Calculators.Tests/Deltares.DamEngine.Calculators.Tests.csproj (.../Deltares.DamEngine.Calculators.Tests.csproj) (revision 680) +++ dam engine/trunk/src/Deltares.DamEngine.Calculators.Tests/Deltares.DamEngine.Calculators.Tests.csproj (.../Deltares.DamEngine.Calculators.Tests.csproj) (revision 681) @@ -45,6 +45,7 @@ + Index: dam engine/trunk/src/Deltares.DamEngine.Calculators.Tests/KernelWrappers/Common/PlLinesHelperTests.cs =================================================================== diff -u --- dam engine/trunk/src/Deltares.DamEngine.Calculators.Tests/KernelWrappers/Common/PlLinesHelperTests.cs (revision 0) +++ dam engine/trunk/src/Deltares.DamEngine.Calculators.Tests/KernelWrappers/Common/PlLinesHelperTests.cs (revision 681) @@ -0,0 +1,102 @@ +using System; +using System.Collections.Generic; +using System.Linq; +using System.Text; +using System.Threading.Tasks; +using Deltares.DamEngine.Calculators.KernelWrappers.Common; +using Deltares.DamEngine.Data.General; +using Deltares.DamEngine.Data.General.PlLines; +using Deltares.DamEngine.Data.Geometry; +using Deltares.DamEngine.Data.Geotechnics; +using NUnit.Framework; + +namespace Deltares.DamEngine.Calculators.Tests.KernelWrappers.Common +{ + [TestFixture] + public class PlLinesHelperTests + { + [Test] + public void TestCreatePlLines() + { + // expected values retrieved from debugging Dam Classic (rev.663) + // test CanCalculateThePipingFactorUsingSellmeijer4Forces + var location = new Location(); + location.SurfaceLine = CreateSurfaceLineTutorial1(); + var soilProfile1D = CreatePipingSellmeijerProfileWithOneSandlayer(); + var riverLevel = 1.0; + var plLines = PlLinesHelper.CreatePlLines(location, soilProfile1D, riverLevel); + Assert.AreEqual(4, plLines.PLLineCount); + Assert.AreEqual(8, plLines.Lines[PLLineType.PL1].Points.Count); + Assert.IsTrue(new PLLinePoint(0.000, 1.000).LocationEquals(plLines.Lines[PLLineType.PL1].Points[0])); + Assert.IsTrue(new PLLinePoint(14.900, 1.000).LocationEquals(plLines.Lines[PLLineType.PL1].Points[1])); + Assert.IsTrue(new PLLinePoint(34.500, 0.500).LocationEquals(plLines.Lines[PLLineType.PL1].Points[2])); + Assert.IsTrue(new PLLinePoint(40.500, 0.000).LocationEquals(plLines.Lines[PLLineType.PL1].Points[3])); + Assert.IsTrue(new PLLinePoint(50.500, 0.000).LocationEquals(plLines.Lines[PLLineType.PL1].Points[4])); + Assert.IsTrue(new PLLinePoint(58.500, 0.000).LocationEquals(plLines.Lines[PLLineType.PL1].Points[5])); + Assert.IsTrue(new PLLinePoint(61.500, 0.000).LocationEquals(plLines.Lines[PLLineType.PL1].Points[6])); + Assert.IsTrue(new PLLinePoint(75.000, 0.000).LocationEquals(plLines.Lines[PLLineType.PL1].Points[7])); + Assert.AreEqual(0, plLines.Lines[PLLineType.PL2].Points.Count); + Assert.AreEqual(2, plLines.Lines[PLLineType.PL3].Points.Count); + Assert.IsTrue(new PLLinePoint(0.000, 1.000).LocationEquals(plLines.Lines[PLLineType.PL3].Points[0])); + Assert.IsTrue(new PLLinePoint(75.000, 1.000).LocationEquals(plLines.Lines[PLLineType.PL3].Points[1])); + Assert.AreEqual(0, plLines.Lines[PLLineType.PL4].Points.Count); + } + + private static SoilProfile1D CreatePipingSellmeijerProfileWithOneSandlayer() + { + SoilProfile1D soilProfile1D = new SoilProfile1D(); + SoilLayer1D soilLayer1D1 = new SoilLayer1D(); + soilLayer1D1.Name = "L0"; + soilLayer1D1.TopLevel = 10.0; + soilLayer1D1.Soil = new Soil("Topmaterial", 1.0, 1.0); + soilLayer1D1.Soil.PermeabKx = 0.0003; + soilLayer1D1.Soil.DiameterD70 = 0.0003; + soilLayer1D1.Soil.WhitesConstant = 0.5; + soilLayer1D1.Soil.BeddingAngle = 57.0; + soilLayer1D1.IsAquifer = false; + soilProfile1D.Layers.Add(soilLayer1D1); + + SoilLayer1D soilLayer1D2 = new SoilLayer1D(); + soilLayer1D2.Name = "L1"; + soilLayer1D2.TopLevel = -2.0; + soilLayer1D2.Soil = new Soil("Sand", 22.0, 20.0); + soilLayer1D2.Soil.PermeabKx = 0.0001; + soilLayer1D2.Soil.DiameterD70 = 0.0002; + soilLayer1D2.Soil.WhitesConstant = 0.25; + soilLayer1D2.Soil.BeddingAngle = 37.0; + soilLayer1D2.IsAquifer = true; + soilProfile1D.Layers.Add(soilLayer1D2); + soilProfile1D.BottomLevel = -10.0; + return soilProfile1D; + } + + private static SurfaceLine2 CreateSurfaceLineTutorial1(bool includingTraffic = false) + { + SurfaceLine2 surfaceLine2 = new SurfaceLine2(); + surfaceLine2.Name = "Tutorial1"; + surfaceLine2.Geometry = new GeometryPointString(); + surfaceLine2.CharacteristicPoints.GeometryMustContainPoint = true; + SurfaceLine2 line = surfaceLine2; + line.EnsurePointOfType(0.0, 0.0, CharacteristicPointType.SurfaceLevelOutside); + line.EnsurePointOfType(10.0, 0.0, CharacteristicPointType.DikeToeAtRiver); + line.EnsurePointOfType(34.5, 5.0, CharacteristicPointType.DikeTopAtRiver); + if (includingTraffic) + { + line.EnsurePointOfType(35.0, 5.0, CharacteristicPointType.TrafficLoadOutside); + line.EnsurePointOfType(38.5, 5.0, CharacteristicPointType.TrafficLoadInside); + } + line.EnsurePointOfType(40.5, 5.0, CharacteristicPointType.DikeTopAtPolder); + line.EnsurePointOfType(50.5, 0.0, CharacteristicPointType.DikeToeAtPolder); + line.EnsurePointOfType(58.5, 0.0, CharacteristicPointType.DitchDikeSide); + line.EnsurePointOfType(59.5, -2.0, CharacteristicPointType.BottomDitchDikeSide); + line.EnsurePointOfType(61.5, -2.0, CharacteristicPointType.BottomDitchPolderSide); + line.EnsurePointOfType(61.5, 0.0, CharacteristicPointType.DitchPolderSide); + line.EnsurePointOfType(75.0, 0.0, CharacteristicPointType.SurfaceLevelInside); + line.Geometry.SyncCalcPoints(); + return line; + } + + + + } +}