Elastik qovushqoq plastinkaning kuchlanish-deformatsion holatini matematik modellashtirish

}
public void Sigma_Oss()
{
double ss;
uint[] mi = new uint[3];
// double[,] S = new double[n, n];
// double[] F = new double[n];
listBox1.Items.Add("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-**-*-*-*-*--*-*-*-*");
listBox1.Items.Add("\t\t\t sigma");

for (int i = 0; i < m; i++)
{
s = "";
mi[0] = MN[i].i; mi[1] = MN[i].j; mi[2] = MN[i].k;
E =E* 100;
ss = E / ((1 + myu) * (1 - 2 * myu));
sigma[i, 0] = ss * (eps[i, 0] - alfa * MK[mi[0] - 1].T + myu * (-eps[i, 0] + eps[i, 1] + eps[i, 2] - alfa * MK[mi[0] - 1].T));
sigma[i, 1] = ss * (eps[i, 1] - alfa * MK[mi[0] - 1].T + myu * (eps[i, 0] - eps[i, 1] + eps[i, 2] - alfa * MK[mi[0] - 1].T));
sigma[i, 2] = ss * (eps[i, 2] - alfa * MK[mi[0] - 1].T + myu * (eps[i, 0] + eps[i, 1] - eps[i, 2] - alfa * MK[mi[0] - 1].T));
sigma[i, 3] = (E * eps[i, 3]) / (2*(1+myu));
s += String.Format("{0} {1:f5} {2:f5} {3:f5} {4:f5} ", i + 1, sigma[i, 0],sigma[i, 1],sigma[i, 2],sigma[i, 3]);
listBox1.Items.Add(s);
}
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++) S[i, j] = 0;
F[i] = 0;
}
double[,] k = new double[4, 4];
double[] f = new double[4], SSigma;
uint k1, k2,sq=12;
Console.WriteLine("**************************************");
for (int k0 = 0; k0 < m; k0++)
{
mi[0] = MN[k0].i; mi[1] = MN[k0].j; mi[2] = MN[k0].k;

k[0, 0] = A[k0] * 2.0 / sq; k[0, 1] = A[k0] * 1.0 / sq; k[0, 2] = A[k0] * 1.0 / sq; k[0, 3] = A[k0] * 1.0 / sq;
k[1, 0] = A[k0] * 1.0 / sq; k[1, 1] = A[k0] * 2.0 / sq; k[1, 2] = A[k0] * 1.0 / sq; k[1, 3] = A[k0] * 1.0 / sq;
k[2, 0] = A[k0] * 1.0 / sq; k[2, 1] = A[k0] * 1.0 / sq; k[2, 2] = A[k0] * 2.0 / sq; k[2, 3] = A[k0] * 1.0 / sq;