private static WMesh2D_Mesh Mesh_FreeRim2(ref WEntity2D[] WBs, ref WEntity2D[] Ws, ref GLCommon GLC, string Name)
{
/////找到基础线
WEntity2D Bas = WBs[0];
for (int i = 0; i < WBs.Length; i++)
{
if (((WCurve2D)WBs[i]).Meshed_Check == true)
{
Bas = WBs[i];
break;
}
}
//////找到与Bas相交的夹丝线
WEntity2D[] Wsc = new WEntity2D[0];
for (int i = 0; i < Ws.Length; i++)
{
if (Geos2D_Intersection.Intersection_L2C((WLine2D)Ws[i], Bas).Count > 0)
{
Array.Resize <WEntity2D>(ref Wsc, Wsc.Length + 1);
Wsc[Wsc.Length - 1] = Ws[i];
}
}
/////
return(null);
}
/// 对夹丝边界线进行排序，从左至右输出竖向线
private static WEntity2D[] Sort_WBs(WLine2D AxisX, ref WEntity2D[] WBs, ref GLCommon GLC)
{
int Quan = WBs.Length;
WPoint2D[] Ps = new WPoint2D[0]; /////夹丝边界线与横轴的交点点集
List <WPoint2D> Pst;
WPoint2D P;
int[] Ns = new int[0];
int N = 0;
for (int i = 0; i < Quan; i++)
{
Pst = Geos2D_Intersection.Intersection_L2C(AxisX, WBs[i]);
if (Pst.Count == 0)
{
continue;
}
P = Pst[0];
Array.Resize <WPoint2D>(ref Ps, N + 1);
Array.Resize <int>(ref Ns, N + 1);
Ps[N] = P;
Ns[N] = i;
N++;
}
///////得到竖向线并排序输出
int[] Sort = WMFuncs2D.Ps_Sort(ref Ps, true, true);
WEntity2D[] Out = new WEntity2D[N];
for (int i = 0; i < N; i++)
{
Out[i] = WBs[Ns[Sort[i]]];
}
///////去除竖向线
for (int i = 0; i < N; i++)
{
for (int j = 0; j < WBs.Length; j++)
{
if (Ns[i] == j)
{
for (int k = j; k < WBs.Length - 1; k++)
{
WBs[k] = WBs[k + 1];
}
Array.Resize <WEntity2D>(ref WBs, WBs.Length - 1);
for (int k = i + 1; k < N; k++)
{
Ns[k] = Ns[k] - 1;
}
break;
}
}
}
return(Out);
}
/// 对夹丝线进行排序，从左至右输出
private static WEntity2D[] Sort_Wires(WLine2D AxisX, ref WEntity2D[] Wires, ref GLCommon GLC)
{
List <WPoint2D> Ps_Wires = Geos2D_Intersection.Intersection_L2Cvs(AxisX, ref Wires); ///////得到交点
int Quan = Wires.Length;
int[] PNums = new int[Quan]; /////用于记录排列顺序
for (int i = 0; i < Quan; i++)
{
PNums[i] = i;
}
int Nt;
WPoint2D Pt; /////排列点
for (int i = 0; i <= Quan - 2; i++)
{
for (int j = i + 1; j <= Quan - 1; j++)
{
if (Ps_Wires[i].X > Ps_Wires[j].X)
{
Nt = PNums[i];
PNums[i] = PNums[j];
PNums[j] = Nt;
Pt = Ps_Wires[i];
Ps_Wires[i] = Ps_Wires[j];
Ps_Wires[j] = Pt;
}
}
}
WEntity2D[] ID_Wires_Sorted = new WEntity2D[Quan]; /////按顺序记录排列后的线
int N;
for (int i = 0; i <= Quan - 1; i++)
{
N = PNums[i];
ID_Wires_Sorted[i] = Wires[N];
}
Ps_Wires = null;
return(ID_Wires_Sorted);
}
/// <summary>
/// 生成夹丝区域的节点，单元，输出，绘制等
/// </summary>
/// <param name="AxisX">横向轴线</param>
/// <param name="WBs">夹丝线</param>
/// <param name="Ws">夹丝区域的上下边界线</param>
/// <remarks></remarks>
private static WMesh2D_Mesh Mesh_SingleRim(ref WLine2D AxisX,
ref WEntity2D[] WBs, ref WEntity2D[] Ws,
ref GLCommon GLC)
{
WMesh2D_Mesh Glass_Mesh = new WMesh2D_Mesh("Wires");
Glass_Mesh.Shapes.Add(new WShapeRim2D(Color.Black, 1f));
/////
int N = 0; /////节点编号
int E4 = 0; /////Shell单元编号
int E2 = 0; /////Beam单元编号
WEntity2D[] WBs_V = Sort_WBs(AxisX, ref WBs, ref GLC); /////找到竖向边界排序输出，并在边界线数组中去除竖向边界线
if (WBs_V.Length == 0)
{
return(null);
}
WCurve2D Bd_L = (WCurve2D)WBs_V[0]; /////左侧边界
WCurve2D Bd_R = (WCurve2D)WBs_V[1]; /////右侧边界
WCurve2D[] Bds_O = new WCurve2D[WBs.Length]; /////上下边界线
for (int i = 0; i < WBs.Length; i++)
{
Bds_O[i] = (WCurve2D)WBs[i];
}
///////先求出上下边界
List <WPoint2D> Ps_WiL = new List <WPoint2D>(); /////记录所有线与上边界的交点，从左至右
List <WPoint2D> Ps_WiU = new List <WPoint2D>(); /////记录所有线与下边界的交点，从左至右
List <WPoint2D> Pts; /////辅助点集
for (int i = 0; i <= Ws.Length - 1; i++)
{
Pts = Geos2D_Intersection.Intersection_L2Cvs((WLine2D)Ws[i], ref WBs); /////找与上下边线的交点作为划分的边界
if (Pts.Count == 0)
{
continue; /////夹丝线与上下边界不一定会有交点
}
Add_Ps2_Bnds(ref Pts, ref Bds_O);
Geos2D_Other.Merge_Points(ref Pts); /////去掉多余节点
if (Pts[0].Y > Pts[1].Y)
{
Ps_WiU.Add(Pts[0]);
Ps_WiL.Add(Pts[1]);
}
else
{
Ps_WiU.Add(Pts[1]);
Ps_WiL.Add(Pts[0]);
}
}
///////对上下左右边界进行网格划分//////
Bd_L.Nodes = Geos2D_Modify.DotCurve_Times(Bd_L, GLC.Mesh_WireTimes); /////左边界划分
Bd_R.Nodes = Geos2D_Modify.DotCurve_Times(Bd_R, GLC.Mesh_WireTimes); /////右边界划分
WMFuncs2D.Ps_Sort(Bd_L.Nodes, false, true); /////从下到上排列左边界节点
WMFuncs2D.Ps_Sort(Bd_R.Nodes, false, true); /////从下到上排列右边界节点
//////Shell//////
List <int> Ns_WiL = new List <int>(); //////记录所有线与上边界的交点的节点编号(Shell)，从左至右
List <int> Ns_WiU = new List <int>(); //////记录所有线与下边界的交点的节点编号(Shell)，从左至右
N = WMFuncs2D.NsID_App(ref Bd_L.Nodes, ref Bd_L.Nodes_num, N); //////将左边界的节点进行划分并输入Bnd中(Shell)
N = WMFuncs2D.NsID_App(ref Ps_WiU, ref Ns_WiU, N); //////将上边界的节点进行划分(Shell)
N = WMFuncs2D.NsID_App(ref Bd_R.Nodes, ref Bd_R.Nodes_num, N, true); //////将右边界的节点进行划分并输入Bnd中(Shell)，注意需要反向
N = WMFuncs2D.NsID_App(ref Ps_WiL, ref Ns_WiL, N, true); //////将下边界的节点进行划分，注意需要反向(Shell)
//Add_HT2_Bnds(ref Bd_L.Nodes, ref Bds_O); //////将左边界头尾节点加入上下边界中（为了节点顺序对，右边界在最后输入）
//////输出
Glass_Mesh.Q_FreeNs = Bd_L.Nodes.Count * 2 + Ns_WiL.Count * 2; /////自由节点数量输出(Shell)
Glass_Mesh.Add_Ns(ref Bd_L.Nodes, ref Bd_L.Nodes_num, false, 0); /////输出左边界的节点至Shell节点文件中
Glass_Mesh.Add_Ns(ref Ps_WiU, ref Ns_WiU, false, 0); /////输出上边界的节点至Shell节点文件中
Glass_Mesh.Add_Ns(ref Bd_R.Nodes, ref Bd_R.Nodes_num, true, 0); /////输出右边界的节点至Shell节点文件中，注意反向
Glass_Mesh.Add_Ns(ref Ps_WiL, ref Ns_WiL, true, 0); /////输出下边界的节点至Shell节点文件中，注意反向
/////将边界点放到Shape中
Glass_Mesh.Shapes[0].Add_Ps(ref Bd_L.Nodes, ref Bd_L.Nodes_num, false, 0); /////输出左边界的节点至Shell节点文件中
Glass_Mesh.Shapes[0].Add_Ps(ref Ps_WiU, ref Ns_WiU, false, 0); /////输出上边界的节点至Shell节点文件中
Glass_Mesh.Shapes[0].Add_Ps(ref Bd_R.Nodes, ref Bd_R.Nodes_num, true, 0); /////输出右边界的节点至Shell节点文件中，注意反向
Glass_Mesh.Shapes[0].Add_Ps(ref Ps_WiL, ref Ns_WiL, true, 0); /////输出下边界的节点至Shell节点文件中，注意反向
/////对夹丝线进行网格划分//////
List <WPoint2D> Ps1 = Bd_L.Nodes;
List <WPoint2D> Ps2 = new List <WPoint2D>();
List <int> Ns1 = Bd_L.Nodes_num;
List <int> Ns2 = new List <int>();
WElement2D[] Es = new WElement2D[0];
for (int i = 0; i <= Ws.Length - 1; i++)
{
Ps2 = Geos2D_Modify.DotCurve_Times(Ws[i], Ps_WiU[i], Ps_WiL[i], GLC.Mesh_WireTimes, ref GLC.WGC); /////划分夹丝线
WMFuncs2D.Ps_Sort(Ps2, false, true); /////从下到上排列节点
//Add_HT2_Bnds(ref Ps2, ref Bds_O); /////将头尾节点加入上下边界中
//////Shell//////
N = WMFuncs2D.NsID_App(ref Ps2, ref Ns2, false, N); /////生成Shell单元的节点编号（注意头尾不要）
Ns2[0] = Ns_WiL[i]; /////为生成单元，用下边界节点替换最下点编号
Ns2[Ns2.Count - 1] = Ns_WiU[i]; /////为生成单元，用上边界节点替换最上点编号
E4 = WMFuncs2D.E4sID_App(ref Ns1, ref Ns2, ref Es, E4); /////生成Shell单元
Glass_Mesh.AddCross_NodesPair(ref Ns1, ref Ns2); /////将节点之间架桥，用于显示
/////输出
Glass_Mesh.Add_Ns(ref Ps2, ref Ns2, true, false); /////将Shell的中间节点输出至文件（注意头尾不要）
Glass_Mesh.Add_Es(ref Es); /////输出Shell单元
E2 = WMFuncs2D.E2sID_App(ref Ns2, ref Es, E2, true); /////Beam单元生成
Glass_Mesh.Add_Es(ref Es); /////Beam单元输出
//////循环推进//////
Ps1 = Ps2;
Ns1 = Ns2;
}
///////最右侧单元处理//////
Ps2 = Bd_R.Nodes;
Ns2 = Bd_R.Nodes_num;
//Add_HT2_Bnds(ref Bd_R.Nodes, ref Bds_O); /////将右侧边界头尾节点加入上下边界中
E4 = WMFuncs2D.E4sID_App(ref Ns1, ref Ns2, ref Es, E4); /////输入单元编号
Glass_Mesh.AddCross_NodesPair(ref Ns1, ref Ns2); /////将节点之间架桥，用于显示
Glass_Mesh.Add_Es(ref Es);
return(Glass_Mesh);
}
/// 找到所有的边界线及所有的区域
private static List <int>[] Find_BdsRims(ref WEntity2D[] WBs, ref WEntity2D[] Ws, ref GLCommon GLC, ref WEntity2D[] Bds, WMesh2D_Mesh Mesh)
{
int N = 0; /////记录节点编号
List <WPoint2D> Pts; /////辅助点集
int Q, Q2;
/////找到每条边界线与夹丝线的交点
WCurve2D C;
for (int i = 0; i < WBs.Length; i++)
{
C = (WCurve2D)(WBs[i]);
if (C.Meshed_Check == true)
{
Geos2D_Other.Merge_Points(ref C.Nodes); /////去掉多余节点
Pts = C.Nodes;
if (WMFuncs2D.Length_P2P(C.StartPoint, Pts[0]) > WMFuncs2D.Length_P2P(C.EndPoint, Pts[0]))
{
WMFuncs2D.Points_Reverse(ref Pts, true);
}
if (WMFuncs2D.Near_Check(Pts[0], C.StartPoint) == false)
{
Pts.Insert(0, C.StartPoint);
}
if (WMFuncs2D.Near_Check(Pts[Pts.Count - 1], C.EndPoint) == false)
{
Pts.Add(C.EndPoint);
}
C.Nodes = Pts;
/////
Divide_Ps(Pts, ref N, ref Bds, -1, Mesh, true, GLC.Mesh_Length); /////分割点集为线段集合
}
else
{
Pts = Geos2D_Intersection.Intersection_Cv2Ls(WBs[i], ref Ws); /////找与上下边线的交点作为划分的边界
if (Pts.Count < 2)
{
Array.Resize <WEntity2D>(ref Bds, Bds.Length + 1);
Bds[Bds.Length - 1] = WBs[i];
((WCurve2D)Bds[Bds.Length - 1]).Trace = -1;
WMFuncs2D.Mesh_Curve(Bds[Bds.Length - 1], GLC.Mesh_Length);
((WCurve2D)Bds[Bds.Length - 1]).Nodes_num = Get_NodeNums(((WCurve2D)Bds[Bds.Length - 1]).Nodes, Mesh, ref N);
continue;
}
else
{
Geos2D_Other.Merge_Points(ref Pts); /////去掉多余节点
}
if (WMFuncs2D.Length_P2P(C.StartPoint, Pts[0]) > WMFuncs2D.Length_P2P(C.EndPoint, Pts[0]))
{
WMFuncs2D.Points_Reverse(ref Pts, true);
}
if (WMFuncs2D.Near_Check(Pts[0], C.StartPoint) == false)
{
Pts.Insert(0, C.StartPoint);
}
if (WMFuncs2D.Near_Check(Pts[Pts.Count - 1], C.EndPoint) == false)
{
Pts.Add(C.EndPoint);
}
/////
Divide_Ps(Pts, ref N, ref Bds, i, Mesh, false, GLC.Mesh_Length); /////分割点集为线段集合
}
}
Mesh.Q_FreeNs = N;
/////找到所有夹丝线的交点，形成新的直线加入
for (int i = 0; i < Ws.Length; i++)
{
Pts = Geos2D_Intersection.Intersection_L2Cvs((WLine2D)Ws[i], ref WBs); /////找与上下边线的交点作为划分的边界
if (Pts.Count == 0)
{
continue; /////夹丝线与上下边界不一定会有交点
}
Geos2D_Other.Merge_Points(ref Pts); /////去掉多余节点
if (Pts.Count == 1)
{
continue;
}
/////
Array.Resize <WEntity2D>(ref Bds, Bds.Length + 1);
Q = Bds.Length - 1;
Bds[Q] = new WLine2D(Pts[0], Pts[1]);
Bds[Q].Trace = -2;
WMFuncs2D.Mesh_Curve(Bds[Q], GLC.Mesh_Length); /////生成该线的节点
N = WMFuncs2D.NsID_App(ref ((WCurve2D)Bds[Q]).Nodes, ref ((WCurve2D)Bds[Q]).Nodes_num, false, N); /////Mesh该线,注意头尾不要
Mesh.Add_Ns(ref ((WCurve2D)Bds[Q]).Nodes, ref ((WCurve2D)Bds[Q]).Nodes_num, true, false); /////将节点编号输入至Mesh中,注意头尾不要
((WCurve2D)Bds[Q]).Nodes_num[0] = Get_NodeNum(((WCurve2D)Bds[Q]).Nodes[0], Mesh, ref N);
Q2 = ((WCurve2D)Bds[Q]).Nodes_num.Count - 1;
((WCurve2D)Bds[Q]).Nodes_num[Q2] = Get_NodeNum(((WCurve2D)Bds[Q]).Nodes[Q2], Mesh, ref N);
}
return(Geo2D_CurvesRims.Get_Rims(ref Bds, 1));
}
