public void ResultantForAngle(GlobalDiagram globaldiagram)
{
Boolean flip = false;
Plane planeXY = new Plane(new Point3d(0, 0, 0), new Vector3d(0, 0, 1));
Line resultantFlipped = ResultantForm;
resultantFlipped.Flip();
Vector3d m0 = globaldiagram.MemberLinesForm[0].Direction;
Vector3d m1 = globaldiagram.MemberLinesForm[1].Direction;
if (Vector3d.VectorAngle(m0, m1, planeXY) < Math.PI)
{
if (Vector3d.VectorAngle(m0, resultantFlipped.Direction, planeXY) < Vector3d.VectorAngle(m0, m1, planeXY))
{
flip = true;
}
}
else
{
if (Vector3d.VectorAngle(m1, resultantFlipped.Direction, planeXY) < Vector3d.VectorAngle(m1, m0, planeXY))
{
flip = true;
}
}
ResultantForceForAngle = ResultantForce;
if (flip)
{
ResultantForceForAngle.Flip();
}
return;
}
protected override void SolveInstance(IGH_DataAccess DA)
{
Functions functions = new Functions();
Point3d iStartForceDiagram = new Point3d(0, 0, 0);
double iScalingFactor = 100.0;
List <Point3d> iJoints = new List <Point3d>();
List <int> iMemberStartIndices = new List <int>();
List <int> iMemberEndIndices = new List <int>();
List <int> iForceIndices = new List <int>();
List <double> iForceMagnitudes = new List <double>();
List <double> iForceRotations = new List <double>();
List <int> iSupportIndices = new List <int>();
List <bool> iSupportHorizontals = new List <bool>();
List <bool> iSupportVerticals = new List <bool>();
List <double> iSupportRotations = new List <double>();
bool iDisplayUnifiedDiagram = true;
bool iDisplayForceDiagram = true;
bool iDisplayJointForceDiagram = true;
double iSpacingJoinForceDiagram = 1.0;
bool iDisplayNumericalValues = true;
DA.GetData(0, ref iStartForceDiagram);
DA.GetData(1, ref iScalingFactor);
DA.GetDataList(2, iJoints);
DA.GetDataList(3, iMemberStartIndices);
DA.GetDataList(4, iMemberEndIndices);
DA.GetDataList(5, iForceIndices);
DA.GetDataList(6, iForceMagnitudes);
DA.GetDataList(7, iForceRotations);
DA.GetDataList(8, iSupportIndices);
DA.GetDataList(9, iSupportHorizontals);
DA.GetDataList(10, iSupportVerticals);
DA.GetDataList(11, iSupportRotations);
DA.GetData(12, ref iDisplayUnifiedDiagram);
DA.GetData(13, ref iDisplayForceDiagram);
DA.GetData(14, ref iDisplayJointForceDiagram);
DA.GetData(15, ref iSpacingJoinForceDiagram);
DA.GetData(16, ref iDisplayNumericalValues);
if ((iSupportRotations.Count != iSupportIndices.Count) && (iSupportRotations.Count == 1))
{
iSupportRotations = functions.CreateDefaultList(iSupportIndices.Count, iSupportRotations[0]);
}
if ((iForceRotations.Count != iForceIndices.Count) && (iForceRotations.Count == 1))
{
iForceRotations = functions.CreateDefaultList(iForceIndices.Count, iForceRotations[0]);
}
List <Line> ioSupportLinesForce = new List <Line>();
List <Line> oSupportLinesForm = new List <Line>();
List <Line> oExtForceLinesForm = new List <Line>();
List <Line> oMemberLinesForm = new List <Line>();
List <Line> oSupportLinesForceJoint = new List <Line>();
List <Line> oExtForceLinesForceJoint = new List <Line>();
List <Line> oMemberLinesForceJoint = new List <Line>();
List <Line> oSupportLinesForce = new List <Line>();
List <Line> oExtForceLinesForce = new List <Line>();
List <Line> oMemberLinesForce = new List <Line>();
List <Brep> oDisplayBreps = new List <Brep>();
List <Color> oMemberColors = new List <Color>();
List <Double> oForceMagnitudes = new List <Double>();
List <Plane> oLocationsForceTextTags = new List <Plane>();
List <Point3d> oLocationsJointForceDiagrams = new List <Point3d>();
List <Line> initialSupportLinesForm = new List <Line>();
List <Joint> myJoints = new List <Joint>();
List <ExternalForce> myExternalForces = new List <ExternalForce>();
List <SupportReaction> mySupportReactions = new List <SupportReaction>();
double ratio = 1.0 / 3.0;
Vector3d zPostive = new Vector3d(0, 0, 1);
Vector3d yPositive = new Vector3d(0, 1, 0);
Vector3d yNegative = new Vector3d(0, -1, 0);
Vector3d mySupportVector;
List <int> extForceInputLengths = new List <int>();
extForceInputLengths.Add(iForceIndices.Count);
extForceInputLengths.Add(iForceMagnitudes.Count);
extForceInputLengths.Add(iForceRotations.Count);
List <int> supportInputLengths = new List <int>();
supportInputLengths.Add(iSupportIndices.Count);
supportInputLengths.Add(iSupportHorizontals.Count);
supportInputLengths.Add(iSupportVerticals.Count);
supportInputLengths.Add(iSupportRotations.Count);
List <int> memberInputLengths = new List <int>();
memberInputLengths.Add(iMemberStartIndices.Count);
memberInputLengths.Add(iMemberEndIndices.Count);
List <int> allMaxIndices = new List <int>();
allMaxIndices.Add(iMemberStartIndices.Max());
allMaxIndices.Add(iMemberEndIndices.Max());
allMaxIndices.Add(iForceIndices.Max());
allMaxIndices.Add(iSupportIndices.Max());
bool noErrors = true;
List <Member> allMembers = new List <Member>();
for (int i = 0; i < iMemberStartIndices.Count; i++)
{
allMembers.Add(new Member(iMemberStartIndices[i], iMemberEndIndices[i], iJoints));
oMemberLinesForm.Add(allMembers[i].FormLine);
}
if (functions.MatchingInputs(supportInputLengths) == false)
{
AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "Support related inputs don't have matching list lengths!");
noErrors = false;
}
if (functions.MatchingInputs(extForceInputLengths) == false)
{
AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "External force related inputs don't have matching list lengths!");
noErrors = false;
}
if (functions.MatchingInputs(memberInputLengths) == false)
{
AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "Member related inputs don't have matching list lengths!");
noErrors = false;
}
if (functions.IndicesInBounds(iJoints.Count, allMaxIndices) == false)
{
AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "One or more of the point indices of the members, forces and/or supports inputs are out of bounds!");
noErrors = false;
}
if (noErrors)
{
for (int i = 0; i < iJoints.Count; i++)
{
myJoints.Add(new Joint(iJoints[i]));
}
for (int i = 0; i < iForceIndices.Count; i++)
{
ExternalForce myExternalForce = new ExternalForce(iForceIndices[i], iJoints, yNegative, iForceRotations[i], iForceMagnitudes[i] / iScalingFactor);
myExternalForce.LengthenFormLines(ratio);
myExternalForces.Add(myExternalForce);
oExtForceLinesForm.Add(myExternalForce.FormLine);
}
for (int i = 0; i < iSupportIndices.Count; i++)
{
if (iSupportHorizontals[i])
{
mySupportVector = yPositive;
mySupportVector.Rotate(iSupportRotations[i] / 180.0 * Math.PI, zPostive);
mySupportVector.Rotate(-0.5 * Math.PI, zPostive);
SupportReaction mySupportReaction = new SupportReaction(iSupportIndices[i], iJoints, mySupportVector);
mySupportReactions.Add(mySupportReaction);
initialSupportLinesForm.Add(mySupportReaction.FormLine);
}
if (iSupportVerticals[i])
{
mySupportVector = yPositive;
mySupportVector.Rotate(iSupportRotations[i] / 180.0 * Math.PI, zPostive);
SupportReaction mySupportReaction = new SupportReaction(iSupportIndices[i], iJoints, mySupportVector);
mySupportReactions.Add(mySupportReaction);
initialSupportLinesForm.Add(mySupportReaction.FormLine);
}
}
for (int i = 0; i < mySupportReactions.Count; i++)
{
if ((mySupportReactions[i].JointIndex != mySupportReactions[(i + 1) % mySupportReactions.Count].JointIndex) && (mySupportReactions[i].JointIndex != mySupportReactions[(i + 2) % mySupportReactions.Count].JointIndex))
{
if (mySupportReactions[i].FormLine.DistanceTo(mySupportReactions[(i + 1) % 3].Joint, false) == 0)
{
AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "unstable structure: two support reactions act in the same line!");
noErrors = false;
break;
}
}
}
if (iSupportIndices.Count != 2)
{
AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "invalid structure: only designs with exactly two supports are supported in the current version of this tool!");
noErrors = false;
}
int degreeOfExtDeterminacy = functions.ExternalStaticDeterminacy(mySupportReactions.Count, 3);
if (degreeOfExtDeterminacy < 0)
{
AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "global structure is unstable!");
noErrors = false;
}
else if (degreeOfExtDeterminacy > 0)
{
AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "global structure is statically indeterminate!");
noErrors = false;
}
if (noErrors)
{
Point3d averagePointJoints = functions.CenterPoint(iJoints);
averagePointJoints.Y += 4;
Resultant myResultant = new Resultant(iStartForceDiagram, myExternalForces, ratio, averagePointJoints);
GlobalDiagram myGlobalDiagram = new GlobalDiagram(myResultant, mySupportReactions);
List <Member> myVirtMembers = myGlobalDiagram.Members();
List <HalfMember> myHalfMembers = myGlobalDiagram.GlobalHalfMembers;
for (int i = 0; i < mySupportReactions.Count; i++)
{
mySupportReactions[i].SupportLineForAngle(myHalfMembers);
}
for (int i = 0; i < myGlobalDiagram.GlobalJoints.Count; i++)
{
myGlobalDiagram.SolveForceDiagram(myVirtMembers);
}
for (int i = 0; i < mySupportReactions.Count; i++)
{
ioSupportLinesForce.Add(mySupportReactions[i].ForceLine);
double supportLengthForm = mySupportReactions[i].Force * ratio;
double extendValue = supportLengthForm - 1.0;
mySupportReactions[i].FormLine.Extend(extendValue, 0);
initialSupportLinesForm[i] = mySupportReactions[i].FormLine;
if (mySupportReactions[i].PositiveForce)
{
mySupportReactions[i].FormLine.Flip();
}
oSupportLinesForm.Add(mySupportReactions[i].FormLine);
}
List <Line> iExtForceLinesForce = myResultant.ForceLinesForce;
//Solve Member Forces
int degreeOfIntDeterminacy = functions.InternalStaticDeterminacy(iMemberStartIndices.Count, mySupportReactions.Count, iJoints.Count);
if (degreeOfIntDeterminacy < 0)
{
AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "internal truss structure is unstable!");
oMemberColors.Add(Color.FromName("Black"));
}
else if (degreeOfIntDeterminacy > 0)
{
AddRuntimeMessage(GH_RuntimeMessageLevel.Error, "internal truss structure is statically indeterminate!");
oMemberColors.Add(Color.FromName("Black"));
}
else
{
List <SupportReaction> allSupports = new List <SupportReaction>();
for (int i = 0; i < oSupportLinesForm.Count; i++)
{
allSupports.Add(new SupportReaction(oSupportLinesForm[i], ioSupportLinesForce[i], iJoints));
}
List <ExternalForce> allExternalForces = new List <ExternalForce>();
for (int i = 0; i < oExtForceLinesForm.Count; i++)
{
allExternalForces.Add(new ExternalForce(oExtForceLinesForm[i], iExtForceLinesForce[i], iJoints));
}
GeneralDiagram myGeneralDiagram = new GeneralDiagram(allExternalForces, allSupports, allMembers, iJoints, iSpacingJoinForceDiagram);
for (int i = 0; i < myGeneralDiagram.AllJoints.Count; i++)
{
//myGeneralDiagram.SolveForceDiagram();
myGeneralDiagram.SolveForceDiagramJointBased();
}
allMembers = myGeneralDiagram.AllMembers;
for (int i = 0; i < allMembers.Count; i++)
{
oForceMagnitudes.Add(allMembers[i].Force * iScalingFactor);
if (allMembers[i].PositiveForce == false)
{
oMemberColors.Add(Color.FromArgb(122, Color.FromName("Blue")));
}
else
{
oMemberColors.Add(Color.FromArgb(122, Color.FromName("Red")));
}
}
if (iDisplayForceDiagram)
{
myGeneralDiagram.CreateOverallForceDiagram();
for (int i = 0; i < allMembers.Count; i++)
{
oMemberLinesForce.Add(allMembers[i].ForceLine);
}
for (int i = 0; i < allSupports.Count; i++)
{
oSupportLinesForce.Add(allSupports[i].ForceLine);
}
for (int i = 0; i < allExternalForces.Count; i++)
{
oExtForceLinesForce.Add(allExternalForces[i].ForceLine);
}
}
if (iDisplayJointForceDiagram)
{
oMemberLinesForceJoint = myGeneralDiagram.ForceLinesJoint();
for (int i = 0; i < myGeneralDiagram.AllJoints.Count; i++)
{
oLocationsJointForceDiagrams.Add(myGeneralDiagram.AllJoints[i].InitSolveLocation);
}
for (int i = 0; i < allSupports.Count; i++)
{
oSupportLinesForceJoint.Add(allSupports[i].ForceLineJoint);
}
for (int i = 0; i < allExternalForces.Count; i++)
{
oExtForceLinesForceJoint.Add(allExternalForces[i].ForceLineJoint);
}
}
if (iDisplayUnifiedDiagram)
{
double myScalingFactor = 1.0 / 9.0;
functions.DisplayRectangles(allMembers, myScalingFactor, out oDisplayBreps);
}
for (int i = 0; i < allSupports.Count; i++)
{
oForceMagnitudes.Add(allSupports[i].Force * iScalingFactor);
}
for (int i = 0; i < allExternalForces.Count; i++)
{
oForceMagnitudes.Add(allExternalForces[i].Force * iScalingFactor);
}
if (iDisplayNumericalValues)
{
List <Plane> tempLocations;
functions.DisplayNumericalValues(oMemberLinesForm, out tempLocations);
for (int i = 0; i < tempLocations.Count; i++)
{
oLocationsForceTextTags.Add(tempLocations[i]);
}
functions.DisplayNumericalValues(oSupportLinesForm, out tempLocations);
for (int i = 0; i < tempLocations.Count; i++)
{
oLocationsForceTextTags.Add(tempLocations[i]);
}
functions.DisplayNumericalValues(oExtForceLinesForm, out tempLocations);
for (int i = 0; i < tempLocations.Count; i++)
{
oLocationsForceTextTags.Add(tempLocations[i]);
}
}
}
}
}
DA.SetDataList(0, oSupportLinesForm);
DA.SetDataList(1, oExtForceLinesForm);
DA.SetDataList(2, oMemberLinesForm);
DA.SetDataList(3, oSupportLinesForceJoint);
DA.SetDataList(4, oExtForceLinesForceJoint);
DA.SetDataList(5, oMemberLinesForceJoint);
DA.SetDataList(6, oSupportLinesForce);
DA.SetDataList(7, oExtForceLinesForce);
DA.SetDataList(8, oMemberLinesForce);
DA.SetDataList(9, oDisplayBreps);
DA.SetDataList(10, oMemberColors);
DA.SetDataList(11, oForceMagnitudes);
DA.SetDataList(12, oLocationsForceTextTags);
DA.SetDataList(13, oLocationsJointForceDiagrams);
}