Elastik qovushqoq plastinkaning kuchlanish-deformatsion holatini matematik modellashtirish

K[MN[e].j * 2 - 2, MN[e].j * 2 - 2] += k[2, 2];
K[MN[e].j * 2 - 2, MN[e].j * 2 - 1] += k[2, 3];
K[MN[e].j * 2 - 2, MN[e].k * 2 - 2] += k[2, 4];
K[MN[e].j * 2 - 2, MN[e].k * 2 - 1] += k[2, 5];

K[MN[e].j * 2 - 1, MN[e].i * 2 - 2] += k[3, 0];
K[MN[e].j * 2 - 1, MN[e].i * 2 - 1] += k[3, 1];
K[MN[e].j * 2 - 1, MN[e].j * 2 - 2] += k[3, 2];
K[MN[e].j * 2 - 1, MN[e].j * 2 - 1] += k[3, 3];
K[MN[e].j * 2 - 1, MN[e].k * 2 - 2] += k[3, 4];
K[MN[e].j * 2 - 1, MN[e].k * 2 - 1] += k[3, 5];
//----------------------------------------
K[MN[e].k * 2 - 2, MN[e].i * 2 - 2] += k[4, 0];
K[MN[e].k * 2 - 2, MN[e].i * 2 - 1] += k[4, 1];
K[MN[e].k * 2 - 2, MN[e].j * 2 - 2] += k[4, 2];
K[MN[e].k * 2 - 2, MN[e].j * 2 - 1] += k[4, 3];
K[MN[e].k * 2 - 2, MN[e].k * 2 - 2] += k[4, 4];
K[MN[e].k * 2 - 2, MN[e].k * 2 - 1] += k[4, 5];

K[MN[e].k * 2 - 1, MN[e].i * 2 - 2] += k[5, 0];
K[MN[e].k * 2 - 1, MN[e].i * 2 - 1] += k[5, 1];
K[MN[e].k * 2 - 1, MN[e].j * 2 - 2] += k[5, 2];
K[MN[e].k * 2 - 1, MN[e].j * 2 - 1] += k[5, 3];
K[MN[e].k * 2 - 1, MN[e].k * 2 - 2] += k[5, 4];
K[MN[e].k * 2 - 1, MN[e].k * 2 - 1] += k[5, 5];

F[MN[e].i * 2 - 2] += fe[0];
F[MN[e].i * 2 - 1] += fe[1];
F[MN[e].j * 2 - 2] += fe[2];
F[MN[e].j * 2 - 1] += fe[3];
F[MN[e].k * 2 - 2] += fe[4];
F[MN[e].k * 2 - 1] += fe[5];

}
public void Inversiya_set(int elm, double Pr, double Pz)
{
dx[elm, 0] = Pr;
//dx[elm, 1] = FE;
dx[elm, 1] = Pz;
}

public void Maxkamlash(int elm, int tur)
{
if (tur == 1) mx[2 * elm - 2] = 1;
else if (tur == 2) mx[2 * elm - 1] = 1;