all[]
101325
1
1
1
10
2
1e-05
0.5
0
2
xr2 = 1.5^2;
xcen = 5.0;
ycen = 0.;
zcen = 6.0;
xray2 = (x-xcen)^2 + (y-ycen)^2 + (z-zcen)^2;
mesh_vi1 = 1;
if (xray2 < xr2) mesh_vi1 = 1e10;
1
1
0
1
0
1
0.001
1017.24
0.02495
0
0
0
0
0
0
basse
ext
haute
int
0
0
0
0
0
0
0
0
0
0
0
0
m11=1;
m22=1;
m33=1;
m12=0;
m13=0;
m23=0;
m21=0;
m31=0;
m32=0;
c11=1;
c22=1;
c33=1;
c12=0;
c13=0;
c23=0;
c21=0;
c31=0;
c32=0;
k11=0.0001;
k22=0.0001;
k33=0.0001;
k12=0;
k13=0;
k23=0;
k21=0;
k31=0;
k32=0;
fx = 0;
fy = 0;
fz = 0;
0.01
100
0
10
0.7
1