//nodes count ox
/// <summary>
/// Разбиение пластинки на трехмерные элементы
/// </summary>
/// <param name="box">Форма пластинки</param>
/// <param name="NcntOX">Количество узлов по иси X</param>
/// <param name="NcntOY">Количество узлов по иси Y</param>
/// <param name="NcntOZ">Количество узлов по иси Z</param>
public void MeshPlate(Box box, int nCntOX, int nCntOY, int nCntOZ)
{
/*формирование сетки для квадратной пластинки*/
//количество узлов на слое
int cntSloi = nCntOX * nCntOY;
//общее количество узлов модели
int obCountUzl = cntSloi * nCntOZ;
//количество КВАДРАТНЫХ элементов на слое
int cntElemSloi = (nCntOX - 1) * (nCntOY - 1);
//общее количество ТЕТРАИДАЛЬНЫХ элементов модели
int countElem = cntElemSloi * (nCntOZ - 1) * 6;
//массив всех узлов
nodes = new Node[obCountUzl];
//массив конечных элементов
elements = new ElementType.Tetra4N[countElem];
bondary = new Bondary();
double stepOX; //шаг по OX
double stepOY; //шаг по OY
double stepOZ; //шаг по ОZ
stepOX = box.lenght / (nCntOX - 1);
stepOY = box.width / (nCntOY - 1);
stepOZ = box.height / (nCntOZ - 1);
double ox = 0,
oy = 0,
oz = 0;
int kk = 0;
// int st = 0;
for (int k = 0; k < nCntOZ; k++) // цикл по OZ
{
oy = 0;
for (int j = 0; j < nCntOY; j++) //цикл по OY
{
ox = 0;
for (int i = 0; i < nCntOX; i++) //цикл по OX
{
//st = 0; // статус узла
if (j == 0 || j == (nCntOY - 1) || i == 0 || i == (nCntOX - 1))
bondary.around.Add(kk);
nodes[kk] = new Node(ox, oy, oz);
kk++;
ox += stepOX;
}
oy += stepOY;
}
oz += stepOZ;
}
// i: cntSloi * iz + iy * nCntOX + ix
// r: cntSloi * iz + iy * nCntOX + (ix+1)
// p: cntSloi * iz + (iy+1) * nCntOX + (ix+1)
// n: cntSloi * iz + (iy+1) * nCntOX + ix
// j: cntSloi * (iz+1) + iy * nCntOX + ix
// s: cntSloi * (iz+1) + iy * nCntOX + (ix+1)
// m: cntSloi * (iz+1) + (iy+1) * nCntOX + (ix+1)
// k: cntSloi * (iz+1) + (iy+1) * nCntOX + ix
kk = 0;
for (int iz = 0; iz < nCntOZ - 1; iz++)
for (int iy = 0; iy < nCntOY - 1; iy++)
{
for (int ix = 0; ix < nCntOX - 1; ix++)
{
int i = cntSloi * iz + iy * nCntOX + ix;
int r = cntSloi * iz + iy * nCntOX + (ix + 1);
int p = cntSloi * iz + (iy + 1) * nCntOX + (ix + 1);
int n = cntSloi * iz + (iy + 1) * nCntOX + ix;
int j = cntSloi * (iz + 1) + iy * nCntOX + ix;
int s = cntSloi * (iz + 1) + iy * nCntOX + (ix + 1);
int m = cntSloi * (iz + 1) + (iy + 1) * nCntOX + (ix + 1);
int k = cntSloi * (iz + 1) + (iy + 1) * nCntOX + ix;
//knip
elements[kk] = new ElementType.Tetra4N(n, k, m, j);
kk++;
//ijkm
elements[kk] = new ElementType.Tetra4N(m, p, n, i);
kk++;
//kmpj
elements[kk] = new ElementType.Tetra4N(n, i, j, m);
kk++;
//rsim
elements[kk] = new ElementType.Tetra4N(s, r, p, i);
kk++;
// mpir
elements[kk] = new ElementType.Tetra4N(i, j, m, s);
kk++;
//rjim
elements[kk] = new ElementType.Tetra4N(i, p, s, m);
kk++;
}
}
}
private void btn_ModelGenerate_Click(object sender, EventArgs e)
{
/* shaft = new Shaft(
readDoubleFromTB(tb_sLenght),
readDoubleFromTB(tb_sWidth),
new RGB(0.8f, 0.8f, 0.8f)*/
box = new Box(
readDoubleFromTB(tb_sLenght),
readDoubleFromTB(tb_sWidth),
readDoubleFromTB(tb_sHeight),
new RGB(0.8f, 0.8f, 0.8f));
Render();
}
