Hello everyone,
I am currently trying to simulate on a cube a coupling of two problems.
For this, I have retrieved the mesh given as an example on the site : https://doc.freefem.org/documentation/mesh-generation.html#the-type-mesh3-in-3-dimension
which is :
load “medit”
load “msh3”
func mesh3 Cube (int[int] &NN, real[int, int] &BB, int[int, int] &L){
real x0 = BB(0,0), x1 = BB(0,1);
real y0 = BB(1,0), y1 = BB(1,1);
real z0 = BB(2,0), z1 = BB(2,1);
int nx = NN[0], ny = NN[1], nz = NN[2];
// 2D mesh
mesh Thx = square(nx, ny, [x0+(x1-x0)*x, y0+(y1-y0)*y]);
// 3D mesh
int[int] rup = [0, L(2,1)], rdown=[0, L(2,0)];
int[int] rmid=[1, L(1,0), 2, L(0,1), 3, L(1,1), 4, L(0,0)];
mesh3 Th = buildlayers(Thx, nz, zbound=[z0,z1],
labelmid=rmid, labelup = rup, labeldown = rdown);
return Th;
}
I’m trying to add boundary conditions that correspond to blocking displacement on certain faces. For this, on the same site (same page) I found this information:
By default the labels are :
1: face y=0,
2: face x=1,
3: face y=1,
4: face x=0,
5: face z=0,
6: face z=1
and the region number is 0.
This information allows me to define the boundary conditions so that my cube is isostatic as follows:
+ on(1,Vy=0)
+ on(4,Vx=0)
+ on(5,Vz=0)
Once I run the calculation I don’t get the desired symmetry.
I would like to know if the information on the face labels corresponds to the defined mesh?
I have tried a new combination :
+ on(3,Vy=0)
+ on(2,Vx=0)
+ on(5,Vz=0)
I get completely inconsistent results when I compare the two possibilities of boundary conditions, when it seems to me that I should obtain a “symmetry” between them.
I hope I have been explicit in my explanations and that someone will know or could help me.
Thank you in advance.
Lisa