Elastik qovushqoq plastinkaning kuchlanish-deformatsion holatini matematik modellashtirish

r = Li(e, 1) * MK[MN[e].i - 1].x + Li(e, 2) * MK[MN[e].j - 1].x + Li(e, 3) * MK[MN[e].k - 1].x;
Ri = MK[MN[e].i - 1].x; Rj= MK[MN[e].j - 1].x;
mi[0] = MN[e].i; mi[1] = MN[e].j; mi[2] = MN[e].k;
s = "";
// s += (String.Format("b== {0} {1} {2} c== {3} {4} {5}", b[0], b[1], b[2], c[0], c[1], c[2]));
//listBox1.Items.Add(s);
A[e] =Math.Abs((MK[MN[e].i - 1].x * b[0] + MK[MN[e].j - 1].x * b[1] + MK[MN[e].k - 1].x * b[2]) / 2);
Ni = (1.0 / (2 * A[e])) * (a[0] + rt * b[0] + zt * c[0]);
Nj = (1.0 / (2 * A[e])) * (a[1] + rt * b[1] + zt * c[1]);
Nk = (1.0 / (2 * A[e])) * (a[2] + rt * b[2] + zt * c[2]);
eps[e, 0] = (1 / (2.0 * A[e])) * (c[0] * NN[mi[0] - 1, 1]*1.0 + c[1] *1.0*
NN[mi[1] - 1, 1] + c[2] * NN[mi[2] - 1, 1]*1.0);

eps[e, 1] = (1.0 / (2.0 * A[e])) * (b[0] * NN[mi[0] - 1, 0] * 1.0 + b[1] * 1.0 *
NN[mi[1] - 1, 0] + b[2] * NN[mi[2] - 1, 0] * 1.0);

eps[e, 2] = (1 / (2 * A[e])) * ((Ni * NN[mi[0] - 1, 0] * 1.0 + Nj * NN[mi[1] - 1, 0] * 1.0 +
Nk * NN[mi[2] - 1, 0] * 1.0) / rt);

eps[e, 3] = (1 / (2 * A[e])) * (c[0] * NN[mi[0] - 1, 0] + b[0] *
NN[mi[0] - 1, 1] + c[1] * NN[mi[1] - 1, 0] + b[1] * NN[mi[1] - 1, 1] + c[2] *
NN[mi[2] - 1, 0] + b[2] * NN[mi[2] - 1, 1]);
/* s = "";
s += String.Format("{3} || {0:f5} {1:f5} {2:f5} == {4:f4}", c[0] * NN[mi[0] - 1, 1], c[1] *
NN[mi[1] - 1, 1], c[2] * NN[mi[2] - 1, 1], A[e], (c[0] * NN[mi[0] - 1, 1] + c[1] *
NN[mi[1] - 1, 1] + c[2] * NN[mi[2] - 1, 1]));
s += String.Format("----- {0:f5} {1:f5} {2:f5} ", c[0], c[1], c[2]);
listBox1.Items.Add(s);
s = "";
s += String.Format(" {3} {0:f7} {1:f7} {2:f7} {4:f7} ", eps[e, 0], eps[e, 1], eps[e, 2], e + 1, eps[e, 3]);
listBox1.Items.Add(s);*/
}