Elastik qovushqoq plastinkaning kuchlanish-deformatsion holatini matematik modellashtirish

b = Int32.TryParse(Read[1], out m);
// n = 10; m = 8;
// Agar tipga o'tkazishda biror tipga o'tkazishda TryParse da xatolik yuzaga kelsa:
if (!(a && b))
{
// MessageBox.Show(String.Format("{0} satrda son emas !", 2), "Xatolik", MessageBoxButtons.OK, MessageBoxIcon.Exclamation);
return false;
}
{
n1 = n; m1 = m;
DMK.RowCount = n;
DMN.RowCount = m;
MK = new M_K[n];
MN = new M_N[m];
K = new double[2 * n, 2 * n];
A = new double[m];
U = new double[2 * n];
F = new double[2 * n];
/* x = new double[2 * n];
y = new double[n];
px = new double[n];
dt = new double[2 * n];*/
B = new double[4, 6];
Bt = new double[6, 4];
eps = new double[m, 4];
sigma = new double[m, 4];
S = new double[n, n];
F = new double[2 * n];
D = new double[4, 4];
NN = new double[n, 2];
dx = new double[n, 2];
mx = new int[2 * n];
SSigma = new double[n];
SSigmay=new double[n];
tao=new double[n];
}
v = (E*(1-myu)) /( (1 + myu )*(1-2* myu));
D[0, 0] = D[1, 1] = D[2,2] = v;
D[0, 1] = D[1, 0] = v * myu/(1-myu);
D[0, 2] = D[2, 0] = v * myu / (1 - myu);
D[1, 2] = D[2, 1] = v * myu / (1 - myu);
D[3, 3] = v*((1-2*myu)/(2*(1 - myu)));


// for (i = 1; i < 2 * n; i++) dt[i] = 1;
//dt[0] =dt[1]=dt[14]=dt[15]=dt[28]=dt[29]= 0;
//U[0] = U[1] = U[14] = U[15] = U[28] = U[29] = 0;

//---------------MN ni o'qib olish
for (int i = 0; i < m; i++)
{
s = reader.ReadLine();
s = s.Replace(".", ",");