Error in the compilator

Feature Request
1

Dear all,
I am very new to the FreeFem++ program. So I need your help.
I want to continue to do simulations of bending objects.
However, the following error continues to appear.

– mesh: Nb of Triangles = 1440, Nb of Vertices 972
– FESpace: Nb of Nodes 10692 Nb of DoF 10692
– Solve :
min -0.0471218 max 0.0483716
min -1.05565e-30 max 1.64479
min -0.010061 max 0.00684872
Plot:: Sorry no ps version for this type of plot 5
– FESpace: Nb of Nodes 10692 Nb of DoF 10692

Hear is the code,

include “cube.idp”
load “msh3”
load “tetgen”
load “medit”
load “iovtk”

// Parameters
real L1 = 1;
real L2 = 11;
real H1 = 0.035;
real H2 = 0.285;

real Fx=5;
real Fy;

real Mu1 = 1.e11;
real Lambda1 = 1.e11;
real Mu2 = 1.e4;
real Lambda2 = 1.e4;

int move = 1;
int fixed = 2;

// Mesh
mesh Th0;

border a1(t=0., 1.){x=0; y=H1-(2*H1)*t; label=fixed;};
border a2(t=0., 1.){x=L2+(-L2)*t; y=H1; label=move;};
border a3(t=0., 1.){x=L2; y=-H2+(H2+H1)*t; label=move;};
border a4(t=0., 1.){x=L1+(L2-L1)*t; y=-H2; label=move;};
border a5(t=0., 1.){x=L1; y=-H1+(-H2+H1)*t; label=move;};
border a6(t=0., 1.){x=(L1)*t; y=-H1; label=move;};

border b1(t=0., 1.){x=0; y=-H1+(2*H1)t; label=fixed;};
border b2(t=0., 1.){x=-L2t; y=H1; label=move;};
border b3(t=0., 1.){x=-L2; y=H1+(-H2-H1)*t; label=move;};
border b4(t=0., 1.){x=-L2+(L2-L1)*t; y=-H2; label=move;};
border b5(t=0., 1.){x=-L1; y=-H2+(H2-H1)t; label=move;};
border b6(t=0., 1.){x=-L1+L1t; y=-H1; label=move;};

Th0 = buildmesh (a1(10) + a2(100) + a3(30) + a4(90)+ a5(20) + a6(10)
+ b2(100) + b3(30) + b4(90)+ b5(20) + b6(10));

// Construction of a box with one hole and two regions
func zmin = 0.;
func zmax = 1.;
int MaxLayer = 10;

mesh3 Th = buildlayers(Th0, MaxLayer, zbound=[zmin, zmax]);

fespace Vh(Th,P1);
Vh u1,u2,u3;
Vh v1,v2,v3;
Vh lambda = (0.8 < abs(x) && 1.2 > abs(x) && y < 1. && y > -1.) ? Lambda2 : Lambda1;
Vh mu = (0.8 < abs(x) && 1.2 > abs(x) && y < 1. && y > -1.) ? Mu2 : Mu1;
Vh Fxx = (x < 0) ? -Fx : Fx;

real sqrt2=sqrt(2.);
macro epsilon(u1,u2,u3) [dx(u1),dy(u2),dz(u3),(dz(u2)+dy(u3))/sqrt2,(dz(u1)+dx(u3))/sqrt2,(dy(u1)+dx(u2))/sqrt2] // EOM
macro div(u1,u2,u3) ( dx(u1)+dy(u2)+dz(u3) ) // EOM

problem Lame([u1,u2,u3],[v1,v2,v3])=
int3d(Th)(
lambda*div(u1,u2,u3)*div(v1,v2,v3)
+2.mu( epsilon(u1,u2,u3)‘*epsilon(v1,v2,v3) ) //’)
)

  • int3d(Th) (Fxx*v1)
  • on(fixed,u1=0,u2=0,u3=0)
    ;

for (int i; i < 10; ++i){
Lame;
Th=movemesh3(Th,transfo=[x+u1,y+u2,z+u3]);
plot(Th, wait=0, prev=1);
u1=u1; u2=u2; u3=u3;
}

Thank you for reading and for your help.

Best regards,

Gooyoon