updating first geometries

srcs/clampedMicroBeam.geo

@@ -4,7 +4,25 @@ Ly_poutre = 1;
 h_poutre = 1;
 h_tot = 5;
 width = 35;
-n = 6;
+
+n = 6; // FEM MESH DENSITY
+
+// additional parameters given to the solver
+SetNumber("Boundary Conditions/left_edge/ux", 0.); 
+SetNumber("Boundary Conditions/left_edge/uy", 0);
+SetNumber("Materials/domain/Young", 150e9);
+SetNumber("Materials/domain/Poisson", 0.27);
+SetNumber("Materials/domain/rho",2300);
+SetNumber("Volumic Forces/FEM_domain/bx",0.);
+SetNumber("Volumic Forces/FEM_domain/by",0.); //set to -9.81 for gravity
+
+phi_top = 112; // à modifier valou
+SetNumber("Boundary Conditions/mass/BEM_domain_1/dirichlet", 0);
+SetNumber("Boundary Conditions/BEM_FEM_boundary/BEM_domain_1/dirichlet", phi_top);
+SetNumber("Boundary Conditions/rest_of_outside/BEM_domain_1/neumann", 0);
+SetNumber("Materials/BEM_domain_1/Epsilon", 8.8541878128e-12); // dielectric permittivity
+
+SetNumber("Non_linear_solver",1);
 
 Point(1) = {0, h_poutre*scale, 0, 2};
 Point(2) = {Lx_poutre*scale, h_poutre*scale, 0, 2};
@@ -59,33 +77,10 @@ Physical Point("fixed_node", 9) = {1};
 
 Physical Surface("BEM_domain_1", 10) = {2};
 
-// additional parameters given to the solver
-SetNumber("Boundary Conditions/left_edge/ux", 0.); // HERE YOU DO NOT HAVE TO IMPOSE BOTH ux and uy simultaneously ! (permet aussi de simuler appuis à roulettes)
-SetNumber("Boundary Conditions/left_edge/uy", 0);
-SetNumber("Materials/domain/Young", 15e9);
-SetNumber("Materials/domain/Poisson", 0.27);
-SetNumber("Materials/domain/rho",2300); //volumic mass of acier
-SetNumber("Boundary Conditions/top_edge/tx", 0.); // ALWAYS NEED TO IMPOSE BOTH tx AND ty ON A GIVEN EDGE (realiste, OK) !
-SetNumber("Boundary Conditions/top_edge/ty", 0.); //set to some non-zero value to induce vertical deflection
-//SetNumber("Boundary Conditions/right_edge/tx", 21e3); // for simple tension conditions
-SetNumber("Boundary Conditions/right_edge/ty", 0.);
-SetNumber("Volumic Forces/FEM_domain/bx",0.);
-SetNumber("Volumic Forces/FEM_domain/by",0.); //set to -9.81 for gravity
-//SetNumber("Boundary Conditions/fixed_node/uy",2.);
-
-SetNumber("Non_linear_solver",1);
-
 // BEM geometry
 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};
-phi_top = 35.5; 
-SetNumber("Boundary Conditions/mass/BEM_domain_1/dirichlet", 0);
-SetNumber("Boundary Conditions/BEM_FEM_boundary/BEM_domain_1/dirichlet", phi_top);
-SetNumber("Boundary Conditions/rest_of_outside/BEM_domain_1/neumann", 0);
-SetNumber("Materials/BEM_domain_1/Epsilon", 8.8541878128e-12); // dielectric permittivity
-//Physical Curve("BEM_boundary", 14) = {1,3,5,6,7,8,9,10,2};
-
 Physical Curve("Electrode", 15) = {1};
 //SetNumber("Boundary Conditions/bottom_edge/tx", 0);
 //SetNumber("Boundary Conditions/bottom_edge/ty", -25000);
\ No newline at end of file

srcs/longitudinalCombDevice.geo

 scale = 2e-6;
 
 // USE WITH MINIMUM 2 FINS, else use longitudinal_comb.geo
-N_fins = 10; // number of fins on one side of the comb
+N_fins = 7; // number of fins on one side of the comb
 
 // WARNING: when using more fins the pull-in voltage decreases
 
@@ -12,39 +12,42 @@ SetNumber("Boundary Conditions/left/ux", 0.); // encastrement
 SetNumber("Boundary Conditions/left/uy", 0);
 SetNumber("Boundary Conditions/right/ux", 0.); // encastrement
 SetNumber("Boundary Conditions/right/uy", 0);
-SetNumber("Materials/FEM_domain/Young", 210e9); // A DETERMINER PRECISEMENT
-SetNumber("Materials/FEM_domain/Poisson", 0.3);
-SetNumber("Materials/FEM_domain/rho",7800); //volumic mass of acier
+SetNumber("Materials/FEM_domain/Young", 150e9);
+SetNumber("Materials/FEM_domain/Poisson", 0.27);
+SetNumber("Materials/FEM_domain/rho",2300); 
 SetNumber("Volumic Forces/FEM_domain/bx",0);
 SetNumber("Volumic Forces/FEM_domain/by",00); // acceleration of accelerometer
 
 // BEM properties in bottom domain
-phi_1 = 140;
+phi_1 = 50;
 SetNumber("Boundary Conditions/BEM_FEM_boundary_1/BEM_domain_1/dirichlet", 0);
 SetNumber("Boundary Conditions/electrode_1/BEM_domain_1/dirichlet", phi_1);
 SetNumber("Boundary Conditions/outside_1/BEM_domain_1/neumann", 0);
 SetNumber("Materials/BEM_domain_1/Epsilon", 8.8541878128e-12); // dielectric permittivity
 
 // BEM properties in top domain
-phi_2 = 5;
+phi_2 = 0;
 SetNumber("Boundary Conditions/BEM_FEM_boundary_2/BEM_domain_2/dirichlet", 0);
 SetNumber("Boundary Conditions/electrode_2/BEM_domain_2/dirichlet", phi_2);
 SetNumber("Boundary Conditions/outside_2/BEM_domain_2/neumann", 0);
 SetNumber("Materials/BEM_domain_2/Epsilon", 8.8541878128e-12); // dielectric permittivity
 
 h_tot = 20*scale;
-h_base = 1*scale;
+h_base = 0.4*scale;
 h_fin = 2.8*scale;
-h_space = 1*scale; // space between electrode and base of the comb
+h_space = 1*scale; // space between bout de l'electrode and base of the clamped beam of the comb
 h_fin_elec = 2.4*scale; // length of the fins of the electrode, can be longer than the classical fins
 // WARNING, abs(h_fin_elec - h_fin) must be lower than h_space
 
-l_bord = 10*scale;
+//l_bord = 10*scale;
+//l_tot = 34.4*scale;
+l_tot = 64*scale;
 l_fin = 0.8*scale;
 l_space = 0.8*scale;
 t_electrode = 1*scale; // width of one electrode
 l_periodic = l_fin + 2*l_space + t_electrode; // distance between two fins
-l_tot = 2*l_bord + l_fin + (N_fins-1)*l_periodic;
+//l_tot = 2*l_bord + l_fin + (N_fins-1)*l_periodic;
+l_bord = (l_tot - l_fin - (N_fins-1)*l_periodic)/2;
 
 unit_l = 0.2*scale; // reference for the transfinite curves