From 104502e2aac803b8cda6842628837067ca40f1fe Mon Sep 17 00:00:00 2001
From: Louis Denis <louis.denis@student.uliege.be>
Date: Wed, 18 May 2022 20:55:53 +0200
Subject: [PATCH] updating first geometries
---
srcs/clampedMicroBeam.geo | 43 +++++++++++++++------------------
srcs/longitudinalCombDevice.geo | 23 ++++++++++--------
2 files changed, 32 insertions(+), 34 deletions(-)
diff --git a/srcs/clampedMicroBeam.geo b/srcs/clampedMicroBeam.geo
index ddbf181..5e94f6d 100644
--- a/srcs/clampedMicroBeam.geo
+++ b/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
diff --git a/srcs/longitudinalCombDevice.geo b/srcs/longitudinalCombDevice.geo
index c28d1f7..2e77d64 100644
--- a/srcs/longitudinalCombDevice.geo
+++ b/srcs/longitudinalCombDevice.geo
@@ -1,7 +1,7 @@
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
--
GitLab