/// <summary>
/// Create FdCoordinateSystem from Rhino plane on surface mid u/v-point.
/// </summary>
internal static FdCoordinateSystem FromRhinoSurface(this Rhino.Geometry.Surface obj)
{
// reparameterize if necessary
if (obj.Domain(0).T0 == 0 && obj.Domain(1).T0 == 0 && obj.Domain(0).T1 == 1 && obj.Domain(1).T1 == 1)
{
// pass
}
else
{
obj.SetDomain(0, new Rhino.Geometry.Interval(0, 1));
obj.SetDomain(1, new Rhino.Geometry.Interval(0, 1));
}
Rhino.Geometry.Plane plane;
obj.FrameAt(0.5, 0.5, out plane);
return(plane.FromRhinoPlane());
}
public Rhino.Commands.Result AddPoints(Rhino.RhinoDoc doc)
{
// Input interval
double _interval = 0.1;
double _noise = 1; // 1% noise
var input = new Rhino.Input.Custom.GetNumber();
input.SetCommandPrompt("Input interval(0.0 ~ 1.0)<0.1>");
Rhino.Input.GetResult res = input.Get();
if (res == Rhino.Input.GetResult.Number)
{
_interval = input.Number();
}
if (_interval == 0.0)
{
_interval = 0.1;
}
input.SetCommandPrompt("Noise factor(0.0 ~ 100.0)");
res = input.Get();
if (res == Rhino.Input.GetResult.Number)
{
_noise = input.Number();
}
// Select surface
var go = new Rhino.Input.Custom.GetObject();
go.SetCommandPrompt("Select surface");
go.GeometryFilter = Rhino.DocObjects.ObjectType.Surface;
res = go.GetMultiple(1, 1024);
if (res == Rhino.Input.GetResult.Nothing)
{
return(Rhino.Commands.Result.Failure);
}
Utility.SetOutputCount(10);
for (int i = 0; i < go.ObjectCount; i++)
{
Rhino.Geometry.Surface surfaceA = go.Object(i).Surface();
if (surfaceA == null)
{
return(Rhino.Commands.Result.Failure);
}
Rhino.Geometry.Interval domU = surfaceA.Domain(0);
Rhino.Geometry.Interval domV = surfaceA.Domain(1);
double u, v;
u = v = 0.0;
for (u = domU.Min; u <= domU.Max; u += _interval)
{
for (v = domV.Min; v <= domV.Max; v += _interval)
{
Rhino.Geometry.Point3d pt = surfaceA.PointAt(u, v);
Rhino.Geometry.Vector3d n = surfaceA.NormalAt(u, v);
Rhino.Geometry.Point3d pt2 = randomPoint(surfaceA, u, v, _interval, _noise);
doc.Objects.AddPoint(pt2);
}
}
}
return(Rhino.Commands.Result.Success);
}
protected override void SolveInstance(Grasshopper.Kernel.IGH_DataAccess DA)
{
if (!DA.GetData(3, ref v))
{
return;
}
if (!FriedChiken.isInitialized)
{
Rhino.Geometry.Surface s = null;
if (!DA.GetData(0, ref s))
{
return;
}
Rhino.Geometry.Interval uDomain = s.Domain(0);
Rhino.Geometry.Interval vDomain = s.Domain(1);
int[] nEdgeNodes = new int[_dim];
DA.GetData(1, ref nEdgeNodes[0]);
DA.GetData(2, ref nEdgeNodes[1]);
for (int i = 0; i < _dim; i++)
{
if (nEdgeNodes[i] < 2)
{
AddRuntimeMessage(Grasshopper.Kernel.GH_RuntimeMessageLevel.Error, "Integers must be greater than or equal to 2");
return;
}
}
newNodes.Clear();
for (int i = 0; i < nEdgeNodes[1]; i++)
{
for (int j = 0; j < nEdgeNodes[0]; j++)
{
newNodes.Add(s.PointAt(uDomain.T0 + (uDomain.T1 - uDomain.T0) / (nEdgeNodes[0] - 1) * j, vDomain.T0 + (vDomain.T1 - vDomain.T0) / (nEdgeNodes[1] - 1) * i));
}
}
int nNewNodes = newNodes.Count;
mikity.NumericalMethodHelper.particle[] particles = new mikity.NumericalMethodHelper.particle[nNewNodes];
for (int i = 0; i < nNewNodes; i++)
{
particles[i] = new mikity.NumericalMethodHelper.particle(newNodes[i][0], newNodes[i][1], newNodes[i][2]);
}
pS = new GH_particleSystem(particles);
node[] lNodes = new node[nNewNodes];
for (int i = 0; i < nNewNodes; i++)
{
lNodes[i] = new node(i);
lNodes[i].copyFrom(pS.Value.particles);
}
nF = new nodalForce(v.X, v.Y, v.Z);
for (int i = 0; i < nNewNodes; i++)
{
nF.addNode(lNodes[i]);
}
pS.Value.addObject(nF);
lGeometry = new Rhino.Geometry.Point3d[nNewNodes];
for (int i = 0; i < nNewNodes; i++)
{
lGeometry[i] = new Rhino.Geometry.Point3d(particles[i][0], particles[i][1], particles[i][2]);
}
}
else
{
nF.forceX = v.X;
nF.forceY = v.Y;
nF.forceZ = v.Z;
}
DA.SetData(0, pS);
}
protected override void SolveInstance(Grasshopper.Kernel.IGH_DataAccess DA)
{
if (!FriedChiken.isInitialized)
{
Rhino.Geometry.Surface s = null;
if (!DA.GetData(0, ref s))
{
return;
}
Rhino.Geometry.Interval uDomain = s.Domain(0);
Rhino.Geometry.Interval vDomain = s.Domain(1);
int[] nEdgeNodes = new int[_dim];
DA.GetData(1, ref nEdgeNodes[0]);
DA.GetData(2, ref nEdgeNodes[1]);
for (int i = 0; i < _dim; i++)
{
if (nEdgeNodes[i] < 2)
{
AddRuntimeMessage(Grasshopper.Kernel.GH_RuntimeMessageLevel.Error, "Integers must be greater than or equal to 2");
return;
}
}
Rhino.Geometry.Mesh m = new Rhino.Geometry.Mesh();
//メッシュノード構築
newNodes.Clear();
for (int i = 0; i < nEdgeNodes[1]; i++)
{
for (int j = 0; j < nEdgeNodes[0]; j++)
{
newNodes.Add(s.PointAt(uDomain.T0 + (uDomain.T1 - uDomain.T0) / (nEdgeNodes[0] - 1) * j, vDomain.T0 + (vDomain.T1 - vDomain.T0) / (nEdgeNodes[1] - 1) * i));
}
}
m.Vertices.AddVertices(newNodes);
int nNewNodes = newNodes.Count;
GH_material mat = null;
GH_gravity gvt = null;
if (!DA.GetData(3, ref mat))
{
return;
}
if (!DA.GetData(4, ref gvt))
{
return;
}
el = MathUtil.isoparametricElements(nEdgeNodes);
int nElements = el.Length;
//メッシュ構築
for (int i = 0; i < nElements; i++)
{
m.Faces.AddFace(el[i][0], el[i][1], el[i][3], el[i][2]);
}
mikity.NumericalMethodHelper.particle[] particles = new mikity.NumericalMethodHelper.particle[nNewNodes];
for (int i = 0; i < nNewNodes; i++)
{
particles[i] = new mikity.NumericalMethodHelper.particle(newNodes[i][0], newNodes[i][1], newNodes[i][2]);
}
eM = new generalSpring();
pS = new GH_particleSystem(particles);
for (int i = 0; i < nElements; i++)
{
if (_subdv == subdivide.quad)
{
eM.addElement(new mikity.NumericalMethodHelper.elements.isoparametricElement(el[i]));
}
if (_subdv == subdivide.triA)
{
eM.addElement(new mikity.NumericalMethodHelper.elements.simplexElement(new int[3] {
el[i][0], el[i][1], el[i][3]
}));
eM.addElement(new mikity.NumericalMethodHelper.elements.simplexElement(new int[3] {
el[i][0], el[i][3], el[i][2]
}));
}
if (_subdv == subdivide.triB)
{
int S = i % (nEdgeNodes[0] - 1);
int T = (i - S) / (nEdgeNodes[0] - 1);
if (T % 2 == 1)
{
S++;
}
if (S % 2 == 0)
{
eM.addElement(new mikity.NumericalMethodHelper.elements.simplexElement(new int[3] {
el[i][0], el[i][1], el[i][3]
}));
eM.addElement(new mikity.NumericalMethodHelper.elements.simplexElement(new int[3] {
el[i][0], el[i][3], el[i][2]
}));
}
else
{
eM.addElement(new mikity.NumericalMethodHelper.elements.simplexElement(new int[3] {
el[i][0], el[i][1], el[i][2]
}));
eM.addElement(new mikity.NumericalMethodHelper.elements.simplexElement(new int[3] {
el[i][2], el[i][1], el[i][3]
}));
}
}
if (_subdv == subdivide.triC)
{
int S = i % (nEdgeNodes[0] - 1);
if (S % 2 == 0)
{
eM.addElement(new mikity.NumericalMethodHelper.elements.simplexElement(new int[3] {
el[i][0], el[i][1], el[i][3]
}));
eM.addElement(new mikity.NumericalMethodHelper.elements.simplexElement(new int[3] {
el[i][0], el[i][3], el[i][2]
}));
}
else
{
eM.addElement(new mikity.NumericalMethodHelper.elements.simplexElement(new int[3] {
el[i][0], el[i][1], el[i][2]
}));
eM.addElement(new mikity.NumericalMethodHelper.elements.simplexElement(new int[3] {
el[i][2], el[i][1], el[i][3]
}));
}
}
}
if (_subdv == subdivide.quad)
{
nElements *= 1;
}
else
{
nElements *= 2;
}
lGeometry = new Rhino.Geometry.Mesh();
lGeometry2 = new Rhino.Geometry.Mesh();
lGeometry.Vertices.Clear();
lGeometry.Faces.Clear();
lGeometry2.Faces.Clear();
for (int i = 0; i < pS.Value.__N; i++)
{
lGeometry.Vertices.Add(particles[i][0], particles[i][1], particles[i][2]);
}
for (int i = 0; i < nElements; i++)
{
if (_subdv == subdivide.quad)
{
lGeometry.Faces.AddFace(eM.elemList[i].el[0], eM.elemList[i].el[1], eM.elemList[i].el[3], eM.elemList[i].el[2]);
}
else
{
lGeometry.Faces.AddFace(eM.elemList[i].el[0], eM.elemList[i].el[1], eM.elemList[i].el[2]);
}
}
for (int i = 0; i < nElements; i++)
{
if (_subdv == subdivide.quad)
{
lGeometry2.Faces.AddFace(eM.elemList[i].el[0], eM.elemList[i].el[1], eM.elemList[i].el[3], eM.elemList[i].el[2]);
}
else
{
lGeometry2.Faces.AddFace(eM.elemList[i].el[0], eM.elemList[i].el[1], eM.elemList[i].el[2]);
}
}
eM.setMaterial(mat.Value, gvt.Value);
pS.Value.addObject(eM);
this.DVPW = GetDVPW(lGeometry);
this.BKGT = GetBKGT(lGeometry);
pS.DVPW = GetDVPW(lGeometry2);
pS.UPGR = GetUPGR(lGeometry2);
pS.BKGT = GetBKGT(lGeometry2);
DA.SetData(0, pS);
DA.SetDataList(1, newNodes);
}
}