Elastik qovushqoq plastinkaning kuchlanish-deformatsion holatini matematik modellashtirish

}
///Console.WriteLine("**************************************");
for (int k0 = 0; k0 < m; k0++)
{
mi[0] = MN[k0].i - 1; mi[1] = MN[k0].j - 1; mi[2] = MN[k0].k - 1;

k[0, 0] = A[k0] * 2.0 / 12.0; k[0, 1] = A[k0] * 1.0 / 12.0; k[0, 2] = A[k0] * 1.0 / 12.0;
k[1, 0] = A[k0] * 1.0 / 12.0; k[1, 1] = A[k0] * 2.0 / 12.0; k[1, 2] = A[k0] * 1.0 / 12.0;
k[2, 0] = A[k0] * 1.0 / 12.0; k[2, 1] = A[k0] * 1.0 / 12.0; k[2, 2] = A[k0] * 2.0 / 12.0;
f[0] = -A[k0] * sigma[k0, 3] * 1.0 / 3.0;
f[1] = -A[k0] * sigma[k0, 3] * 1.0 / 3.0;
f[2] = -A[k0] * sigma[k0, 3] * 1.0 / 3.0;
//------------------------------------------------------------
for (int i = 0; i < 3; i++)
{
k1 = mi[i];
F[k1] = F[k1] + f[i];
for (int j = 0; j < 3; j++)
{
k2 = mi[j];
S[k1, k2] = S[k1, k2] + k[i, j];
}

}
}
tao = Gaus(2);
//listBox1.Items.Add("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-**-*-*-*-*--*-*-*-*");
//listBox1.Items.Add("\t\t\t Sigma");
//for (int i = 0; i < 115; i++)
//{
// if (SSigmay[i] > 100) SSigmay[i] -= 700;
// else SSigmay[i] -= 300;
//}
double ssigmai,eps_i;
for (int i = 0; i < n; i++)

{
s = "";

ssigmai = Math.Sqrt(1.0 / 2) * Math.Sqrt(Math.Pow(SSigma[i] - SSigmay[i], 2.0) + 4 * tao[i] * tao[i]);
eps_i = Math.Sqrt(Math.Pow(eps[i,0] - eps[i, 1], 2.0) + 4 * eps[i, 3] * eps[i,3]);
s += String.Format("{1} {0:f5} {2:f5} {3:f5}", SSigma[i], i + 1, SSigmay[i], tao[i]);
listBox1.Items.Add(s);
s = String.Format("{0:f5} {1:f5} {2:f5} {3:f5} {4:f5} {5:f5} {6:f5}", NN[i, 0], NN[i, 1], SSigma[i], SSigmay[i], tao[i],ssigmai,eps_i);
writer.WriteLine(s);
}
writer.Close();

}

}
}