2. All boundary conditions must be written as `SetNumber("Boundary Conditions/name_of_the_edge/...")` with `name_of_the_edge` the
2. All boundary conditions must be written as `SetNumber("Boundary Conditions/name_of_the_physical_curve/...")`
edge on which the condition is imposed.
with `name_of_the_physical_curve` the physical curve on which the condition is imposed.
3. The Dirichlet boundary conditions, i.e the imposed horizontal and vertical displacement `u_x` and `u_y` of an edge,
3. The Dirichlet boundary conditions, i.e the imposed horizontal and vertical displacement `u_x` and `u_y` of an edge,
must be written as `Setnumber("Boundary Conditions/name_of_the_edge/ux", desired_value);`, same applies for `u_y`.
must be written as `Setnumber("Boundary Conditions/name_of_the_physical_curve/ux", desired_value);`, same applies for `u_y`. Note that the horizontal displacement can be imposed on a physical curve without imposing the vertical
For example, if the user want to imposed an embedding on an the bottom edge:
imposed on an edge, must be written as `SetNumber("Boundary Conditions/name_of_edge/tx", desired_value);`, same applies for
4. The Neumann boundary conditions, i.e surface traction in the horizontal (`t_x`) or vertical (`t_y`) direction imposed on
`t_y`.
an edge, must be written as `SetNumber("Boundary Conditions/name_of_the_physical_curve/tx", desired_value)`,
same applies for `t_y`. Note that on a specific edge, both horizontal and vertical surface tractions must be
prescribed. If the user only wants to prescribe `t_x`, `t_y` must also be set to zero.
5. Dirichlet boundary conditions can also be imposed on a specific point of the domain. In this case, a `Physical Point`
must be created.
-**MATERIAL PROPERTIES**:
-**MATERIAL PROPERTIES**:
All the material properties must be written as `SetNumber("Materials/domain/name_of_property", value_of_property);`.
All the material properties must be written as `SetNumber("Materials/FEM_domain/name_of_property", value_of_property);`.
The different properties that must be specified are:
The different properties that must be specified are:
- Young modulus, named `Young`.
- Young's modulus, named `Young`.
- Poisson ration, named `Poisson`.
- Poisson's ratio, named `Poisson`.
-Volumic density, named `rho`.
-Mass density, named `rho`.
-**VOLUMIC FORCES**:
-**VOLUMIC FORCES**:
Volumic forces, i.e forces applied on all the volume in the horizontal (`b_x`) or vertical (`b_y`) direction
Volumic forces, in the horizontal (`b_x`) or vertical (`b_y`) direction must be written as `SetNumber("Volumic Forces/FEM_domain/b_x",0);`, same applies for `b_y`.
must be written as `SetNumber("Volumic Forces/FEM_domain/b_x",0);`, same applies for `b_y`.
- For the electrostatic BEM part:
- For the BEM domains:
1. The lines corresponding to the boundary of the BEM surfaces must be created in such a way that the area of the BEM surface is located on the left of the Curve Loop.
1. The lines of the BEM surfaces must be created in such a way that the air of the surface is located to the left of the Curve Loop.
2. Multiple BEM domains can be defined. One BEM domain should be defined as a `Physical Surface` with the precise name `BEM_domain_X`, with `X` the identifier of the BEM domain (`X=1` if it is the first domain, `X=2` if it is the second one, ...).
2. BEM domain should be defined as `Physical Curve` with the precise name `BEM_domain_X`, with `X` the number of the BEM domain (`X=1` if it is the first one, `X=2` if it is the seconde one, ...).
3. Boundary conditions and material properties must then be specified for the BEM domain:
For example, if the user want the `Curve Loop(1)` to be a FEM domain:
3. Lines belonging to the BEM and FEM domains must be specified as `Physical Curve("BEM_FEM_boundary_X, x)={Lines_of_the_two_domains}`. For example, if `Line(1)`, `Line(4)` and `Line(5)` belong to the intersection of the BEM domain and the FEM domain:
4. The electrode on line defining the electrode on wich the electric potential is imposed must be defined as `Physical Curve("electrode_X", x)={Lines_corresponding_to_electrode}`.
5. The line which are not defined as an electrode or as a an intersection between the two types of domain must be specified as `Physical Curve("outside_X", x) = {Lines_of_outside};`
6. Boundary conditions and materials properties must then be specified for the Bem domain:
-**BOUNDARY CONDITIONS**:
-**BOUNDARY CONDITIONS**:
- Two Dirichlet boundary conditions need to be imposed:
1. Dirichlet boundary conditions, corresponding to imposing an electric potential, can be imposed on a physical curve as:
1. The electric potential of the electrode, specified as `SetNumber("Boundary Conditions/electrode_X/BEM_domain_X/dirichlet", value_of_potential);`.
2. A condition on the boundary between the FEM and the BEM domain. It is often considered as the place where the electric potential is equal to 0, defined as `SetNumber("Boundary Conditions/BEM_FEM_boundary_X/BEM_domain_X/dirichlet", value_of_potential);`.
2. Neumann boundary conditions, corresponding to imposing the normal component of the electric field, can be imposed on a physical curve as: `SetNumber("Boundary Conditions/name_of_the_physical_curve/BEM_domain_X/neumann", value_of_field);`.
- One Neumann boundary condition on the line which are not defined as an electrode or as a an intersection between the two types of domain, written as `SetNumber("Boundary Conditions/outside_X/BEM_domain_X/neumann", 0);`.
-**MATERIAL PROPERTIES**: The dielectric permittivity of the material inside the BEM domain must be specified. It is done writing `SetNumber("Materials/BEM_domain_X/Epsilon", value_of_Epsilon);`.
-**MATERIAL PROPERTIES**: The dielectric permittivity of the BEM domain must be specified, it is done writing `SetNumber("Materials/BEM_domain_X/Epsilon", value_Epsilon);`.
**Note**: These operations must be done for each BEM domain the user wants to create, replacing `X` by the given identifier of the BEM domain.
**Note**: These operations must be done for each BEM domains the user want to create, replacing `X` by the number of the BEM domain.
- For the coupling between the FEM and BEM domains:
1. All the lines belonging to the interface between the FEM domain and the different BEM domains must be specified as a single `Physical Curve("BEM_FEM_boundary, x)={List_of_the_interface_lines}`. As an example, if `Line(1)`, `Line(4)` and `Line(5)` belong to the intersection of the BEM domains and the FEM domain:
