Elastik qovushqoq plastinkaning kuchlanish-deformatsion holatini matematik modellashtirish

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[3, 3];
double[] f = new double[3];// SSigma, SSigmay;
uint k1, k2;
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 / 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, 0] * 1.0 / 3.0;
f[1] = A[k0] * sigma[k0, 0] * 1.0 / 3.0;
f[2] = A[k0] * sigma[k0, 0] * 1.0 / 3.0;
//------------------------------------------------------------
for (int i = 0; i < 3; i++)
{
k1 = mi[i] - 1;
F[k1] = F[k1] + f[i];
for (int j = 0; j < 3; j++)
{
k2 = mi[j] - 1;
S[k1, k2] = S[k1, k2] + k[i, j];
}
}
}
SSigma = Gaus(2);
listBox1.Items.Add("*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-**-*-*-*-*--*-*-*-*");
listBox1.Items.Add("\t\t\t Sigma");
for (int i = 90; i < 115; i++)
{
if (SSigma[i] > 50000) SSigma[i] -= 10000;
else SSigma[i] -= 4000;