Skip to content
Snippets Groups Projects
Commit 70b972db authored by Denis Louis's avatar Denis Louis
Browse files

generalisation of scale of hybrid geometry

parent 758ed054
No related branches found
No related tags found
No related merge requests found
Lx_poutre = 5;
Ly_poutre = 2;
h_poutre = 2;
h_tot = 10;
width = 10;
nx = 40; // prend beaucoup de temps àpd de 200x40
ny = 16;
scale = 1e-6;
Lx_poutre = 5*scale;
Ly_poutre = 2*scale;
h_poutre = 2*scale;
h_tot = 10*scale;
width = 10*scale;
nx = 10; // prend beaucoup de temps àpd de 200x40
ny = 4;
Point(1) = {0, h_poutre, 0, 2};
Point(2) = {Lx_poutre, h_poutre, 0, 2};
......@@ -75,8 +76,12 @@ SetNumber("Volumic Forces/domain/by",0.); //set to -9.81 for gravity
Physical Curve("BEM_FEM_boundary", 11) = {1,2,3};
Physical Curve("mass", 12) = {10}; //lower electrode
Physical Curve("rest_of_outside", 13) = {5,6,7,8,9,4};
phi_top = 100;
phi_top = 100;
SetNumber("Boundary Conditions/mass/dirichlet", 0);
SetNumber("Boundary Conditions/BEM_FEM_boundary/dirichlet", phi_top);
SetNumber("Boundary Conditions/rest_of_outside/neumann", 0);
Physical Curve("BEM_boundary", 14) = {1,3,5,6,7,8,9,10,2};
\ No newline at end of file
Physical Curve("BEM_boundary", 14) = {1,3,5,6,7,8,9,10,2};
Physical Curve("bottom_edge", 15) = {1};
//SetNumber("Boundary Conditions/bottom_edge/tx", 0);
//SetNumber("Boundary Conditions/bottom_edge/ty", -25000);
\ No newline at end of file
  • Pour Physical Curve("rest_of_outside", 13) = {5,6,7,8,9,4}; supprime le 4, la courbe n'est pas compris dans le domaine BEM :)

0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment