Elastik qovushqoq plastinkaning kuchlanish-deformatsion holatini matematik modellashtirish

for (int j = 0; j < 4; j++)
{
BtD[i1, j] = 0;
for (int g = 0; g < 4; g++) BtD[i1, j] += Bt[i1, g] * D[g, j];
}
}
if (pech)
{

listBox1.Items.Add("**************************************");
listBox1.Items.Add("\t\tB matritsa");

for (int i1 = 0; i1 < 4; i1++)
{
s = "";

for (int j = 0; j < 6; j++)
s += (String.Format("{0:f6} ", B[i1, j]));
listBox1.Items.Add(s);
}
listBox1.Items.Add("\n**************************************");
listBox1.Items.Add("\t\tB_t matritsa");

for (int i1 = 0; i1 < 6; i1++)
{
s = "";
for (int j = 0; j < 4; j++)
s += (String.Format("{0:f6} ", Bt[i1, j]));
listBox1.Items.Add(s);
}


listBox1.Items.Add("**************************************");
listBox1.Items.Add("\t\tD matritsa");
for (int j = 0; j < 4; j++)
{
s = "";

for (int g = 0; g < 4; g++)
s += String.Format("{0:f6} ", D[j, g]);// / (6.4 * Math.Pow(10, 6.0) / 8.0));
listBox1.Items.Add(s);
}

listBox1.Items.Add("**************************************");
listBox1.Items.Add("\t\tB_t_D matritsa");
for (int j = 0; j < 6; j++)
{
s = "";
for (int g = 0; g < 4; g++)
s += String.Format("{0:f6} ", BtD[j, g]);/// (6 * Math.Pow(10, 6.0) / 7.5));
listBox1.Items.Add(s);
}

}
for (int i1 = 0; i1 < 6; i1++)
{
for (int j = 0; j < 6; j++)
{
k[i1, j] = 0;
for (int g = 0; g < 4; g++) k[i1, j] += BtD[i1, g] * B[g, j];
k[i1, j] = k[i1, j] *V;

}
}
double eps0 = alfa * MK[MN[e].i - 1].T;// *E;
Double[] qs = {eps0,eps0,eps0,0};

for (int i1 = 0; i1 < 6; i1++)
{
fe[i1] = 0;
for (int g = 0; g < 4; g++) fe[i1] += BtD[i1, g] * qs[g];
fe[i1] = fe[i1] * V * (1);// (1 - 2*myu);

}
double zij, zjk, zki;
zij = Math.Sqrt(Math.Pow(MK[MN[e].i - 1].x - MK[MN[e].j - 1].x, 2) + Math.Pow(MK[MN[e].i - 1].y - MK[MN[e].j - 1].y, 2));
zjk = Math.Sqrt(Math.Pow(MK[MN[e].k - 1].x - MK[MN[e].j - 1].x, 2) + Math.Pow(MK[MN[e].k - 1].y - MK[MN[e].j - 1].y, 2));
zki = Math.Sqrt(Math.Pow(MK[MN[e].i - 1].x - MK[MN[e].k - 1].x, 2) + Math.Pow(MK[MN[e].i - 1].y - MK[MN[e].k - 1].y, 2));
Ri = 1; Rj = 3;
//*************** plos. napr *******************

//V = V ;
double pi = Math.PI;
fh[0] = Ri * 2 * pi * dx[MN[e].i - 1, 0] * rt * A[e] / Rj;
fh[1] = Ri * 2 * pi * dx[MN[e].i - 1, 1] * rt * A[e] / Rj;
fh[2] = Ri * 2 * pi * dx[MN[e].j - 1, 0] * rt * A[e] / Rj;
fh[3] = Ri * 2 * pi * dx[MN[e].j - 1, 1] * rt * A[e] / Rj;
fh[4] = Ri * 2 * pi * dx[MN[e].k - 1, 0] * rt * A[e] / Rj;
fh[5] = Ri * 2 * pi * dx[MN[e].k - 1, 1] * rt * A[e] / Rj;
if (pech)
{
listBox1.Items.Add("**************************************");
listBox1.Items.Add(String.Format("\t\tF ^e matritsa", e));
for (int j = 0; j < 6; j++)
{
s = "";//(alfa * E * t * MK[MN[e].i - 1].T) / (2.0 * (1 - myu)) * b[0];
s += String.Format(" = {0} ", fe[j]);
listBox1.Items.Add(s);
}
listBox1.Items.Add("**************************************");
listBox1.Items.Add(String.Format("\t\tF ^h matritsa", e));