public Curve ToRhinoCrosssectionProfile(sFrame sb, double t)
{
sCrossSection scs = sb.crossSection;
List <sXYZ> vertice = new List <sXYZ>();
sPlane secPlane = new sPlane(sb.axis.PointAt(t), sb.localPlane.Xaxis, sb.localPlane.Yaxis);
if (scs.sectionType == eSectionType.AISC_I_BEAM)
{
vertice = scs.GetWbeamFaceVertices(secPlane).ToList();
}
else if (scs.sectionType == eSectionType.HSS_REC)
{
vertice = scs.GetHSSRecFaceVertices_Simple(secPlane);
}
else if (scs.sectionType == eSectionType.HSS_ROUND)
{
vertice = scs.GetHSSRoundFaceVertices_Simple(secPlane);
}
else if (scs.sectionType == eSectionType.RECTANGLAR)
{
vertice = scs.GetRecFaceVertices_Simple(secPlane);
}
else if (scs.sectionType == eSectionType.SQUARE)
{
vertice = scs.GetSquareFaceVertices_Simple(secPlane);
}
else if (scs.sectionType == eSectionType.ROUND)
{
vertice = scs.GetRoundFaceVertices_Simple(secPlane);
}
Polyline pl = new Polyline();
for (int i = 0; i < vertice.Count; ++i)
{
pl.Add(this.ToRhinoPoint3d(vertice[i]));
}
pl.Add(this.ToRhinoPoint3d(vertice[0]));
Curve profile = pl.ToNurbsCurve();
return(profile);
}
internal Dyn.PolyCurve ToRhinoCrosssectionProfile(sFrame sb, double t)
{
sCrossSection scs = sb.crossSection;
List <sXYZ> vertice = new List <sXYZ>();
sPlane secPlane = new sPlane(sb.axis.PointAt(t), sb.localPlane.Xaxis, sb.localPlane.Yaxis);
if (scs.sectionType == eSectionType.AISC_I_BEAM)
{
vertice = scs.GetWbeamFaceVertices(secPlane).ToList();
}
else if (scs.sectionType == eSectionType.HSS_REC)
{
vertice = scs.GetHSSRecFaceVertices_Simple(secPlane);
}
else if (scs.sectionType == eSectionType.HSS_ROUND)
{
vertice = scs.GetHSSRoundFaceVertices_Simple(secPlane);
}
else if (scs.sectionType == eSectionType.RECTANGLAR)
{
vertice = scs.GetRecFaceVertices_Simple(secPlane);
}
else if (scs.sectionType == eSectionType.SQUARE)
{
vertice = scs.GetSquareFaceVertices_Simple(secPlane);
}
else if (scs.sectionType == eSectionType.ROUND)
{
vertice = scs.GetRoundFaceVertices_Simple(secPlane);
}
List <Dyn.Point> dpts = new List <Autodesk.DesignScript.Geometry.Point>();
for (int i = 0; i < vertice.Count; ++i)
{
dpts.Add(this.ToDynamoPoint(vertice[i]));
}
return(Dyn.PolyCurve.ByPoints(dpts, true));
}
public Plane ToRhinoPlane(sPlane spl)
{
return(new Plane(ToRhinoPoint3d(spl.origin), ToRhinoVector3d(spl.Xaxis), ToRhinoVector3d(spl.Yaxis)));
}
private void AwaresBeamResult(sFrame sb, ref sResultRange bsRe, double dataLenTol)
{
StatBeamResults br = new StatBeamResults();
sStatConverter conv = new sStatConverter();
StatBeam b = sb.extraData as StatBeam;
b.RecoverForces();
sb.beamWeight = b.Weight;
sb.results = new List <sFrameResult>();
double len = sb.axis.length;
int count = (int)(len / dataLenTol);
if (count < 1)
{
count = 1;
}
double step = 1.0 / (double)(count);
for (int i = 0; i < count + 1; ++i)
{
double tNow = step * i;
br.t = tNow;
b.GetInterpolatedResultsAt(new C_vector(0, 0, -1), this.FEsystem.DeadLoadFactor, br);
sFrameResult bre = new sFrameResult();
bre.moment = new sXYZ(br.MomentL.x, br.MomentL.y, br.MomentL.z);
bre.force = new sXYZ(-br.ForceL.x, br.ForceL.y, br.ForceL.z); /// force X negate?????????
bre.deflection_mm = conv.TosXYZ(b.Csys.LocalToGlobalVector(br.DeflL * 1000));
bre.deflectionLocal_mm = conv.TosXYZ(br.DeflL * 1000);
bre.parameterAt = tNow;
sPlane secPlane = new sPlane(sb.axis.PointAt(tNow), sb.localPlane.Xaxis, sb.localPlane.Yaxis);
List <sXYZ> secVertices = new List <sXYZ>();
if (sb.crossSection.sectionType == eSectionType.AISC_I_BEAM)
{
secVertices = sb.crossSection.GetWbeamFaceVertices(secPlane).ToList();
}
else if (sb.crossSection.sectionType == eSectionType.HSS_REC)
{
secVertices = sb.crossSection.GetHSSRecFaceVertices_Simple(secPlane).ToList();
}
else if (sb.crossSection.sectionType == eSectionType.HSS_ROUND)
{
secVertices = sb.crossSection.GetHSSRoundFaceVertices_Simple(secPlane).ToList();
}
else if (sb.crossSection.sectionType == eSectionType.SQUARE)
{
secVertices = sb.crossSection.GetSquareFaceVertices_Simple(secPlane).ToList();
}
else if (sb.crossSection.sectionType == eSectionType.RECTANGLAR)
{
secVertices = sb.crossSection.GetRecFaceVertices_Simple(secPlane).ToList();
}
else if (sb.crossSection.sectionType == eSectionType.ROUND)
{
secVertices = sb.crossSection.GetRoundFaceVertices_Simple(secPlane).ToList();
}
if (secVertices != null)
{
for (int j = 0; j < secVertices.Count; ++j)
{
sXYZ svp = secVertices[j];
sXYZ localDir = svp - sb.axis.PointAt(tNow);
sXYZ ToLocalZ = localDir.ProjectTo(secPlane.Zaxis);
sXYZ ToLocalY = localDir.ProjectTo(secPlane.Yaxis);
double len_ToLocalZ = ToLocalZ.Z; // vertical like
double len_ToLocalY = ToLocalY.Y; // horizontal like
br.zL = len_ToLocalZ; // vertical like
br.yL = len_ToLocalY; // horizontal like
double axialStress_X = (br.ForceL.x / b.CrossSection.Area);
//double axialStress_Y = ??;
//double axialStress_Z = ??;
double MyStress = ((br.MomentL.y * br.zL) / b.CrossSection.Iyy);
double MzStress = ((br.MomentL.z * br.yL) / b.CrossSection.Izz);
//double MxStress = ??
double stressTest = axialStress_X + MyStress - MzStress; //why negate?...(StatBeamResult does this)
//do I need?
//b.GetSectionPointResult(br);
sFrameSectionResult secRe = new sFrameSectionResult();
secRe.location = svp;
secRe.deflection_mm = conv.TosXYZ(b.Csys.LocalToGlobalVector(br.DeflL * 1000));
secRe.stress_Combined = Math.Abs(stressTest);
secRe.stress_Axial_X = axialStress_X;
secRe.stress_Moment_Y = MyStress;
secRe.stress_Moment_Z = MzStress;
bre.sectionResults.Add(secRe);
}
}
bsRe.UpdateMaxValues(bre);
sb.results.Add(bre);
}
}
public screenPlane()
{
this.screenSelect = sPlane.screen;
}
