Elastik qovushqoq plastinkaning kuchlanish-deformatsion holatini matematik modellashtirish

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

f[0] = A[k0] * sigma[k0, 0] * 1.0 / 3;
f[1] = A[k0] * sigma[k0, 0] * 1.0 / 3;
f[2] = A[k0] * sigma[k0, 0] * 1.0 / 3;
f[3] = A[k0] * sigma[k0, 0] * 1.0 / 3;
//------------------------------------------------------------
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 = 0; i < n; i++)
{
s = "";


s += String.Format("{1} {0:f5} ", SSigma[i], i + 1);
listBox1.Items.Add(s);
}

}


public void Sigma()
{
StreamWriter writer = new StreamWriter("D:\\_Doktarantura\\_Misollar\\Temo Elastika 2\\Natijalar\\Sxx.txt");
double ss;E = -E * 10;
uint[] mi = new uint[3];
// double[,] S = new double[n, n];
// double[] F = new double[n];
Console.WriteLine("**************************************");
for (int i = 0; i < m; i++)
{
s = "";
mi[0] = MN[i].i; mi[1] = MN[i].j; mi[2] = MN[i].k;
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));