public Norma(Model model, double eps, double load)
{
InitializeComponent();
AnaliticalMukha analitical = new AnaliticalMukha(model.shape.L, model.shape.H, model.material.E, model.material.V, load * model.shape.H);
Steema.TeeChart.Styles.Line lineAnalit = new Steema.TeeChart.Styles.Line();
lineAnalit.Title = "Analytical";
lineAnalit.Color = Color.Red;
tChart1.Series.Add(lineAnalit);
Steema.TeeChart.Styles.Line line = new Steema.TeeChart.Styles.Line();
line.Title = "Nonlinear FEM";
line.Color = Color.Green;
tChart1.Series.Add(line);
Vector alfa = new Vector(3);
alfa[0] = model.shape.L;
alfa[1] = alfa[2] = 0;
int iter = 0;
dataGridView1.RowCount = 1;
dataGridView1.ColumnCount = 2;
dataGridView1.Columns[0].Name = "N";
dataGridView1.Columns[1].Name = "(Au,u)";
for (int i = 20; i <= 200; i += 20)
{
model.mesh = new Mesh(model.shape, i);
Vector uDiscrete = SolidMechanicSolver.Solve(model, eps, false, out iter);
//Matrix K = SolidMechanicSolver.CreateHardMatrix(model);
//Vector KNonLinear = SolidMechanicSolver.NonLinearPart(model, uDiscrete);
Vector FNonlinear = SolidMechanicSolver.CreateTotalVector(model);
double norm = FNonlinear * uDiscrete;//uDiscrete * K * uDiscrete + 2 * uDiscrete * KNonLinear;
line.Add(i, norm);
List <string> row = new List <string>();
row.Add(i.ToString());
row.Add(norm.ToString("0.000000000000"));
dataGridView1.Rows.Add(row.ToArray());
/*
* ApproximationResult result = new ApproximationResult(model.shape, uDiscrete);
*
* lineAnalit.Add(i, Vector.Norm(analitical.U(alfa)));
* line.Add(i, Vector.Norm(result.U(alfa)));*/
}
/*Steema.TeeChart.Styles.Line lineU1a = new Steema.TeeChart.Styles.Line();
* lineU1a.Title = "U1 аналітичне";
* tChart1.Series.Add(lineU1a);
* Steema.TeeChart.Styles.Line lineU3a = new Steema.TeeChart.Styles.Line();
* lineU3a.Title = "U3 аналітичне";
* tChart1.Series.Add(lineU3a);
*
* Vector alfa = new Vector(Constants.DIMENSION_COUNT);
* alfa[0] = model.mesh[0].Begin;
* alfa[1] = 0;
* alfa[2] = 0;
* lineU1.Add(model.mesh[0].Begin, result.U(alfa)[0]);
* lineU3.Add(model.mesh[0].Begin, result.U(alfa)[2]);
* lineU1a.Add(model.mesh[0].Begin, analitical.U(alfa)[0]);
* lineU3a.Add(model.mesh[0].Begin, analitical.U(alfa)[2]);
* for (int i = 0; i < N; i++)
* {
* double xCur = model.mesh[i].End;
* alfa = new Vector(Constants.DIMENSION_COUNT);
* alfa[0] = xCur;
* alfa[1] = 0;
* alfa[2] = 0;
* if (xCur > model.shape.L)
* {
* xCur = model.shape.L;
* alfa[0] = xCur;
* }
* Vector resFEM = result.U(alfa);
* Vector resAnalit = analitical.U(alfa);
* lineU1.Add(xCur, resFEM[0]);
* lineU3.Add(xCur, resFEM[2]);
* lineU1a.Add(xCur, resAnalit[0]);
* lineU3a.Add(xCur, resAnalit[2]);
* }
*/
}
private void button1_Click(object sender, EventArgs e)
{
tChart1.Series.Clear();
double E = double.Parse(Etxt.Text); //модуль Юнга
double v = double.Parse(Vtxt.Text); //коефіцієнт Пуасона
double h = double.Parse(Htxt.Text); //товщина
int N = int.Parse(Ntxt.Text); //к-сть розбиттів
double l = double.Parse(Ltxt.Text); //довжина
double K = double.Parse(Ktxt.Text); //кривизна
double eps = double.Parse(EPStxt.Text); //точність
sigma = double.Parse(Ftxt.Text);
int count = 6 * (N + 1);
model = new Model(E, v, h, l, K, N);
load = sigma;//(-1)*CountLoad(sigma);
Force F = new Force(load, h);
model.load.F1[0] = F.F10;
model.load.F1[1] = F.F11;
model.load.F1[2] = F.F12;
model.load.F3[0] = F.F30;
model.load.F3[1] = F.F31;
model.load.F3[2] = F.F32;
bool linear = false;
if (checkBox1.Checked)
{
linear = true;
}
int iter = 0;
Vector uDiscrete = SolidMechanicSolver.Solve(model, eps, linear, out iter);
txtNorm.Text += "Ітерацій: " + iter + "\r\nРозбиття: " + N + "\r\n";
ApproximationResult result = new ApproximationResult(model.shape, uDiscrete);
AnaliticalMarchuk analitical = new AnaliticalMarchuk(l, h, E, v, load);
//AnaliticalMukha analitical = new AnaliticalMukha(l, h, E, v, load);
Steema.TeeChart.Styles.Line lineU1 = new Steema.TeeChart.Styles.Line();
lineU1.Title = "U1";
lineU1.Color = Color.Red;
tChart1.Series.Add(lineU1);
Steema.TeeChart.Styles.Line lineU3 = new Steema.TeeChart.Styles.Line();
lineU3.Title = "U3";
lineU3.Color = Color.Green;
tChart1.Series.Add(lineU3);
/*Steema.TeeChart.Styles.Line lineU1a = new Steema.TeeChart.Styles.Line();
* lineU1a.Title = "U1 аналітичне";
* lineU1a.Color = Color.Aqua;
* tChart1.Series.Add(lineU1a);*/
Steema.TeeChart.Styles.Line lineU3a = new Steema.TeeChart.Styles.Line();
lineU3a.Title = "U3 аналітичне";
lineU3a.Color = Color.Blue;
tChart1.Series.Add(lineU3a);
Vector alfa = new Vector(Constants.DIMENSION_COUNT);
alfa[0] = model.mesh[0].Begin;
alfa[1] = 0;
alfa[2] = 0;
lineU1.Add(model.mesh[0].Begin, result.U(alfa)[0]);
lineU3.Add(model.mesh[0].Begin, result.U(alfa)[2]);
//lineU1a.Add(model.mesh[0].Begin, analitical.U(alfa)[0]);
lineU3a.Add(model.mesh[0].Begin, -0.01);//W2(alfa[0])); //analitical.U(alfa)[2]); //W2(alfa[0]));
for (int i = 0; i < N; i++)
{
double xCur = model.mesh[i].End;
alfa = new Vector(Constants.DIMENSION_COUNT);
alfa[0] = xCur;
alfa[1] = 0;
alfa[2] = 0;
if (xCur > model.shape.L)
{
xCur = model.shape.L;
alfa[0] = xCur;
}
Vector resFEM = result.U(alfa);
//Vector resAnalit = W2(xCur);
lineU1.Add(xCur, resFEM[0]);
lineU3.Add(xCur, resFEM[2]);
//lineU1a.Add(xCur, analitical.U(alfa)[0]);
//lineU3a.Add(xCur, analitical.U(alfa)[2]);
//lineU1a.Add(xCur, resAnalit[0]);
//lineU3a.Add(xCur, W2(xCur));
}
dataGridView1.RowCount = 1;
dataGridView1.ColumnCount = 5;
dataGridView1.Columns[0].Name = "alfa1";
dataGridView1.Columns[1].Name = "U1";
dataGridView1.Columns[2].Name = "U1_analitical";
dataGridView1.Columns[3].Name = "U3";
dataGridView1.Columns[4].Name = "U3_analitical";
double alfaCur = model.mesh[0].Begin;
alfa[0] = alfaCur;
alfa[1] = 0;
alfa[2] = 0;
List <string> row = new List <string>();
row.Add(alfaCur.ToString());
row.Add(result.U(alfa)[0].ToString("0.0000000000"));
row.Add(analitical.U(alfa)[0].ToString("0.0000000000"));
row.Add(result.U(alfa)[2].ToString("0.0000000000"));
row.Add(analitical.U(alfa)[2].ToString("0.0000000000"));
//row.Add(W2(alfa[0]).ToString("0.0000000000"));
dataGridView1.Rows.Add(row.ToArray());
for (int i = 0; i < N; i++)
{
alfaCur = model.mesh[i].End;
alfa = new Vector(Constants.DIMENSION_COUNT);
alfa[0] = alfaCur;
alfa[1] = 0;
alfa[2] = 0;
if (alfaCur > model.shape.L)
{
alfaCur = model.shape.L;
alfa[0] = alfaCur;
}
row = new List <string>();
row.Add(alfaCur.ToString());
row.Add(result.U(alfa)[0].ToString("0.0000000000"));
row.Add(analitical.U(alfa)[0].ToString("0.0000000000"));
row.Add(result.U(alfa)[2].ToString("0.0000000000"));
//row.Add(analitical.U(alfa)[2].ToString("0.0000000000"));
row.Add(W2(alfa[0]).ToString("0.0000000000"));
dataGridView1.Rows.Add(row.ToArray());
}
}
public EquilibriumPath(Model model, double eps)
{
InitializeComponent();
Steema.TeeChart.Styles.Line lineNL = new Steema.TeeChart.Styles.Line();
lineNL.Title = "Нелінійна задача";
lineNL.Color = Color.Red;
tChart1.Series.Add(lineNL);
Steema.TeeChart.Styles.Line lineL = new Steema.TeeChart.Styles.Line();
lineL.Title = "Лінійна задача";
lineL.Color = Color.Blue;
tChart1.Series.Add(lineL);
Vector alfa = new Vector(3);
alfa[0] = model.shape.L;
alfa[1] = alfa[2] = 0;
int iter = 0;
for (double load = 0; load < 1000; load += 20)
{
Force F = new Force((-1) * load, model.shape.H);
model.load.F1[0] = F.F10;
model.load.F1[1] = F.F11;
model.load.F1[2] = F.F12;
model.load.F3[0] = F.F30;
model.load.F3[1] = F.F31;
model.load.F3[2] = F.F32;
Vector uDiscreteL = SolidMechanicSolver.Solve(model, eps, true, out iter);
Vector uDiscreteNL = SolidMechanicSolver.Solve(model, eps, false, out iter);
if (iter < 1000)
{
ApproximationResult resultNL = new ApproximationResult(model.shape, uDiscreteNL);
lineNL.Add(Vector.Norm(resultNL.U(alfa)), load);
ApproximationResult resultL = new ApproximationResult(model.shape, uDiscreteL);
lineL.Add(Vector.Norm(resultL.U(alfa)), load);
}
else
{
break;
}
}
/*
*/
/*
* for (int i = 20; i <= 200; i += 20)
* {
* model.mesh = new Mesh(model.shape, i);
*
*
* lineAnalit.Add(i, Vector.Norm(analitical.U(alfa)));
* line.Add(i, Vector.Norm(result.U(alfa)));
* }
*/
}
