public rm21HorArc(ptsPoint begPt, ptsPoint endPt, Azimuth incomingAzimuth, Double radius) : base(begPt, endPt) { if (0 == utilFunctions.tolerantCompare((endPt - begPt).Length, 0.0, 0.000001)) throw new ArcExceptionZeroLengthNotDefined(); Double tanOffsetToEndPt; ptsRay incomingRay = new ptsRay(); incomingRay.advanceDirection = 1; incomingRay.StartPoint = begPt; incomingRay.HorizontalDirection = incomingAzimuth; tanOffsetToEndPt = incomingRay.getOffset(endPt); if (0 == utilFunctions.tolerantCompare(tanOffsetToEndPt, 0.0, 0.000001)) throw new ArcExceptionZeroDeflectionNotDefined(); this.BeginStation = 0.0; this.Radius = radius; this.BeginAzimuth = incomingAzimuth; deflDirection = Math.Sign(tanOffsetToEndPt); Deflection defl = Deflection.ctorDeflectionFromAngle(90.0, deflDirection); Azimuth azToCenter = incomingAzimuth + defl; ptsVector traverseToCenterVec = new ptsVector(azToCenter, radius); ArcCenterPt = begPt + traverseToCenterVec; this.BeginRadiusVector = this.ArcCenterPt - this.BeginPoint; Azimuth endVecAz = (this.ArcCenterPt - endPt).Azimuth; this.EndRadiusVector = new ptsVector(endVecAz, radius); this.EndPoint = this.ArcCenterPt + this.EndRadiusVector; var deflectionDbl = endVecAz - this.BeginRadiusVector.Azimuth; this.Deflection = new Deflection(deflectionDbl, deflDirection); this.EndAzimuth = this.BeginAzimuth + this.Deflection; // applies to English projects only (for now) this.BeginDegreeOfCurve = HorizontalAlignmentBase.computeDegreeOfCurve(this.Radius); this.EndDegreeOfCurve = this.BeginDegreeOfCurve; this.Length = 100.0 * this.Deflection.getAsRadians() / this.BeginDegreeOfCurve.getAsRadians(); this.Length = Math.Abs(this.Length); this.EndStation = this.BeginStation + this.Length; }
public void appendTangent(ptsPoint TangentEndPoint) { // See "Solving SSA triangles" // http://www.mathsisfun.com/algebra/trig-solving-ssa-triangles.html var lastChainItem = allChildSegments.Last(); if (!(lastChainItem is rm21HorArc)) throw new Exception("Can't add tangent on a tangent segment."); var finalArc = lastChainItem as rm21HorArc; var incomingTanRay = new ptsRay(); incomingTanRay.StartPoint = finalArc.BeginPoint; incomingTanRay.HorizontalDirection = finalArc.BeginAzimuth; int offsetSide = Math.Sign(incomingTanRay.getOffset(TangentEndPoint)); double rad = finalArc.Radius; Azimuth traverseToRevisedCenterPtAzimuth = finalArc.BeginAzimuth + offsetSide * Math.PI/2.0; ptsVector traverseToRevisedCenterPt = new ptsVector(traverseToRevisedCenterPtAzimuth, rad); ptsPoint revCenterPt = finalArc.BeginPoint + traverseToRevisedCenterPt; ptsVector ccToTEPvec = finalArc.ArcCenterPt - TangentEndPoint; ptsAngle rho = Math.Asin(rad / ccToTEPvec.Length); ptsAngle ninetyDegrees = new ptsAngle(); ninetyDegrees.setFromDegreesDouble(90.0); ptsAngle tau = ninetyDegrees - rho; Azimuth ccToPtVecAz = ccToTEPvec.Azimuth; Azimuth arcBegRadAz = finalArc.BeginRadiusVector.Azimuth; ptsCogo.Angle.Deflection outerDefl = ccToPtVecAz.minus(arcBegRadAz); ptsDegree defl = Math.Abs((tau - outerDefl).getAsDegreesDouble()) * offsetSide; Deflection newDeflection = new Deflection(defl.getAsRadians()); finalArc.setDeflection(newDeflection: newDeflection); var appendedLineSegment = new rm21HorLineSegment(finalArc.EndPoint, TangentEndPoint); appendedLineSegment.BeginStation = finalArc.EndStation; appendedLineSegment.Parent = this; allChildSegments.Add(appendedLineSegment); this.EndAzimuth = appendedLineSegment.EndAzimuth; this.EndPoint = appendedLineSegment.EndPoint; restationAlignment(); if (alignmentData.Count > 0) alignmentData.RemoveAt(alignmentData.Count-1); alignmentData.Add(new alignmentDataPacket(alignmentData.Count, finalArc)); alignmentData.Add(new alignmentDataPacket(alignmentData.Count, appendedLineSegment)); }