PROGRAM Resol

uti
on_
Str
uct
ure ;

USEScr
t ;

CONST

nx =5 ;

ny =5 ;

nz =6 ;

nt=20 ;

nd =nx*
ny*
nz;

epsi
=10E-
11;

i
nf=300 ;

TYPE

i
nt =0.
.i
nf ;

mat
1d =ar
ray
[int] ofi
nteger ;

mat
2d =ar
ray
[int,
int] ofi
nteger ;

mat
3d =ar
ray
[int,
int,
int
] ofi
nteger ;

mat
1dr
e=ar
ray
[int] ofr
eal ;

mat
2dr
e=ar
ray
[int,
int] ofr
eal ;

VAR

dat
a :
text ;

i
,j,
k,n,
l :
int
eger ;

num_
nd,
type_
nd :
mat
3d ;

i
nd_
i,i
nd_
j,i
nd_
k :
mat
1d ;

ndv :
mat
2d ;
 Unm1,
Vnm1,
Wnm1,

Un ,
Vn ,
Wn ,

Unp1,
Vnp1,
Wnp1,

F1 ,
F2 ,
F3 :
mat
1dr
e ;

I
d,

U,
V,W,

A1u,
A1v,
A1w,

A2u,
A2v,
A2w,

A3u,
A3v,
A3w :
mat
2dr
e ;

E,
nu,
rho ,

Fs1,
Fs2,
Fs3,
Fs4,
fv,

Lx,
Ly,
Lz,
Lt,

dx,
dy,
dz,
dt,

A ,
B,C,
D :
real ;

{
***
***
***
***
***
***
*PROCEDURE*
***
***
***
***
***
***
***
***
*}

{
***
***
***
***
***
*INI
TIALI
SATI
ON*
***
***
***
***
***
***
***
*}

Pr
ocedur
eIni
ti
ali
sat
ion;

Begi
n
/
///
///
///
///
///
/Fi
chi
er/
///
///
///
///
///

assi
gn(
dat
a,'
j
oel
_resul
t.
dat
'
) ;

r
ewr
it
e(dat
a) ;

cl
ose(
dat
a) ;

/
///
///
///
///
///Mat
ri
ce/
///
///
///
///
///
/

Fori
:=0t
oinfdo

begi
n

i
nd_
i[
i]:
=0;
ind_
j[
i]:
=0;
ind_
k[i
]:=0;

Unm1[
i]:
=0;
Vnm1[
i]:
=0;
Wnm1[
i]:
=0;

Un[
i]:
=0;
Vn[
i]:
=0;
Wn[
i]:
=0;

Unp1[
i]:
=0;
Vnp1[
i]:
=0;
Wnp1[
i]:
=0;

Forj
:=0t
oinfdo

begi
n

I
d[i
,
j]:
=0;

ndv
[i
,j
]:=0;

Fork:
=0t
oinfdo

begi
n

num_
nd[
i,
j
,k]:
=0 ;

t
ype_
nd[
i,
j
,k]:
=0 ;

end;

end;

end;

end;
Pr
ocedur
eIni
ti
al_
mat
;

Begi
n

Fori
:=0t
oinfdo

begi
n

F1[
i]:
=0;
F2[
i]:
=0;
F3[
i]:
=0;

Forj
:=0t
oinfdo

begi
n

A1u[
i,
j
]:=0;
A1v
[i
,j
]:=0;
A1w[
i,
j
]:=0;

A2u[
i,
j
]:=0;
A2v
[i
,j
]:=0;
A2w[
i,
j
]:=0;

A3u[
i,
j
]:=0;
A3v
[i
,j
]:=0;
A3w[
i,
j
]:=0;

end;

end;

end;

Pr
ocedur
eaf
fi
ch_
resul
t(
X,Y,
Z:mat
1dr
e);

begi
n

assi
gn(
dat
a,'
j
oel
_resul
t.
dat
'
) ;

Append(
dat
a) ;

wr
it
eln(
dat
a,'
nd'
:
3,'
U'
:17,
'
_'
,n,
'
V':
25,
'
_'
,n,
'
W':
25,
'
_'
,n)
;

f
orl
:=1t
onddo

begi
n

wr
it
eln(
dat
a,l
:
3,X[
l]
:27,
Y[l
]:
27,
Z[l
]:
27);

end;

wr
it
eln(
dat
a) ;

cl
ose(
dat
a) ;
 end;

/
///
///
///
///
///
///
///Mai
l
lageetNoeud/
///
///
///
///
///
///
///
///
///
///

pr
ocedur
enumer
otat
ion_
nd;

begi
n

assi
gn(
dat
a,'
num_
type_
nd.
dat
'
) ;

r
ewr
it
e(dat
a) ;

wr
it
eln(
dat
a,'
nx='
,
nx,
''
,
'ny
=',
ny,
''
,
'nz='
,
nz)
;

f
ork:
=1t
onzdo

begi
n

f
ori
:=1t
onxdo

begi
n

f
orj
:=1t
onydo

begi
n

/
/Numer
otat
iondesnoeuds

num_
nd[
i,
j
,k] :
=nx*
ny*
(k-
1)+nx*
(j-
1)+i
;

l :
=num_
nd[
i,
j
,k] ;

i
nd_
i[
l] :
=i ;

i
nd_
j[
l] :
=j ;

i
nd_
k[l
] :
=k ;

/
/Ty
pesdesnoeuds

casekof

1:
casej
of

1:casei
of

1:t
ype_
nd[
i,
j
,k]:
=19;
 nx:t
ype_
nd[
i,
j
,k]:
=20;

el
set
ype_
nd[
i,
j
,k]:
=11;

end;
//casei

ny
:casei
of

1:t
ype_
nd[
i,
j
,k]:
=22;

nx:t
ype_
nd[
i,
j
,k]:
=21;

el
set
ype_
nd[
i,
j
,k]:
=13;

end;
//casei

el
secasei
of

1:t
ype_
nd[
i,
j
,k]:
=14;

nx:t
ype_
nd[
i,
j
,k]:
=12;

el
set
ype_
nd[
i,
j
,k]:
=5;

end;
//casei

end;
//casej

nz
:casej
of

1:casei
of

1:t
ype_
nd[
i,
j
,k]:
=23;

nx:t
ype_
nd[
i,
j
,k]:
=24;

el
set
ype_
nd[
i,
j
,k]:
=15;

end;
//casei

ny
:casei
of

1:t
ype_
nd[
i,
j
,k]:
=26;

nx:t
ype_
nd[
i,
j
,k]:
=25;

el
set
ype_
nd[
i,
j
,k]:
=17;

end;
//casei
 el
secasei
of

1:t
ype_
nd[
i,
j
,k]:
=18;

nx:t
ype_
nd[
i,
j
,k]:
=16;

el
set
ype_
nd[
i,
j
,k]:
=6;

end;
//casei

end;
//casej

el
secasej
of

1:casei
of

1:t
ype_
nd[
i,
j
,k]:
=7;

nx:t
ype_
nd[
i,
j
,k]:
=8;

el
set
ype_
nd[
i,
j
,k]:
=1;

end;
//casei

ny
:casei
of

1:t
ype_
nd[
i,
j
,k]:
=10;

nx:t
ype_
nd[
i,
j
,k]:
=9;

el
set
ype_
nd[
i,
j
,k]:
=3;

end;
//casei

el
secasei
of

1:t
ype_
nd[
i,
j
,k]:
=4;

nx:t
ype_
nd[
i,
j
,k]:
=2;

el
set
ype_
nd[
i,
j
,k]:
=27;

end;
//casei

end;
//casej

end;
//casek

wr
it
e(dat
a,l
:
5,'
('
,
ind_
i[
l]
,
',
',
i
nd_
j[
l]
,
',
',
i
nd_
k[l
],
'
,'
,
't
p:'
,
type_
nd[
i,
j
,k]
:3,
'
)'
);
 end;

wr
it
eln(
dat
a);

end;

wr
it
eln(
dat
a);

end;

cl
ose(
dat
a);

end;

pr
ocedur
enoeud_
voi
sin;

v
ar

i
i
,jj
,kk,
im1,
ip1,
jm1,
jp1,
km1,
kp1:
int
eger;

begi
n

f
orl
:=1t
onddo

begi
n

i
m1:
=ind_
i[
l]-1;
ii:
=ind_
i[
l];
ip1:
=ind_
i[
l]+1;

j
m1:
=ind_
j[
l]-1;
jj:
=ind_
j[
l];
jp1:
=ind_
j[
l]+1;

km1:
=ind_
k[l
]-1;
kk:
=ind_
k[l
];kp1:
=ind_
k[l
]+1;

ndv
[1,
l
]:=num_
nd[
ii,
jj,
kk] ;

ndv
[2,
l
]:=num_
nd[
im1,
jj,
kk] ;

ndv
[3,
l
]:=num_
nd[
ip1,
jj,
kk] ;

ndv
[4,
l
]:=num_
nd[
ii,
jm1,
kk] ;

ndv
[5,
l
]:=num_
nd[
ii,
jp1,
kk] ;

ndv
[6,
l
]:=num_
nd[
ii,
jj,
km1] ;

ndv
[7,
l
]:=num_
nd[
ii,
jj,
kp1] ;
 ndv
[8,
l
]:=num_
nd[
im1,
jm1,
kk] ;

ndv
[9,
l
]:=num_
nd[
ip1,
jm1,
kk] ;

ndv
[10,
l
]:=num_
nd[
im1,
jp1,
kk] ;

ndv
[11,
l
]:=num_
nd[
ip1,
jp1,
kk] ;

ndv
[12,
l
]:=num_
nd[
im1,
jj,
km1] ;

ndv
[13,
l
]:=num_
nd[
ip1,
jj,
km1] ;

ndv
[14,
l
]:=num_
nd[
im1,
jj,
kp1] ;

ndv
[15,
l
]:=num_
nd[
ip1,
jj,
kp1] ;

ndv
[16,
l
]:=num_
nd[
ii,
jm1,
km1] ;

ndv
[17,
l
]:=num_
nd[
ii,
jp1,
km1] ;

ndv
[18,
l
]:=num_
nd[
ii,
jm1,
kp1] ;

ndv
[19,
l
]:=num_
nd[
ii,
jp1,
kp1] ;

end;

end;

/
///
///
///
///
///
///
///
/Coef
fi
cient/
///
///
///
///
///
///
///

Pr
ocedur
eCoef
fi
cient
;

begi
n

Lx:
=5 ;
Ly:
=5 ;
Lz:
=50 ;

E:
=210e9;
nu:
=0.
285;
rho:
=7800;

A :
=E*(
1-nu)/(
(1+nu)
*(1-
2*nu)
) ;

B:
=E/(
2*(
1-nu)
) ;

C :
=E/(
2*(
1+nu)*(
1-2*
nu)
) ;
 D :
=E*nu/(
(1+nu)*(
1-2*
nu)
) ;

dx:
=Lx/(
nx-1) ;

dy:
=Ly/(
ny-1) ;

dz:
=Lz/(
nz-1) ;

dt:
=1/
2;

Fs1:
=1e6 ;
Fs2:
=1e6{
1e8};

Fs3:
=1e6 ;
Fs4:
=1E6{
1e8};

end;

Pr
ocedur
eCoef
fi
cient
_mat
ri
ce;

v
ar

k5,

xx,
yy,
zz,
tt,
xy,
xz,
yz:
real ;

Begi
n

Coef
fi
cient;

I
nit
ial
_mat;

k5:
=B-C;

xx:
=dx*
dx;

y
y:=dy
*dy;

zz:
=dz*
dz;

t
t:=dt
*dt;

xy:
=dx*
dy;

xz:
=dx*
dz;

y
z:=dy
*dz;
f
orl
:=1t
onddo

begi
n

F1[
l]:
=(r
ho/
tt)
*(Unm1[
l]-
2*Un[
l]
);

F2[
l]:
=(r
ho/
tt)
*(Vnm1[
l]-
2*Vn[
l]
);

F3[
l]:
=(r
ho/
tt)
*(Wnm1[
l]-
2*Wn[
l]
);

caset
ype_
nd[i
nd_
i[
l],
ind_
j[
l],
ind_
k[l
]]of

1:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*B/
zz-
rho/
tt ;

A1u[
l,
ndv
[2,
l
]]:
=A/
xx ;

A1u[
l,
ndv
[3,
l
]]:
=A/
xx ;

A1u[
l,
ndv
[5,
l
]]:
=2*
k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=B/
zz ;

A1u[
l,
ndv
[7,
l
]]:
=B/
zz ;

A1v
[l
,ndv
[2,
l
]]:
=-2*
k5/
(2*
xy) ;

A1v
[l
,ndv
[3,
l
]]:
=2*
k5/
(2*
xy) ;

A1w[
l,
ndv
[12,
l
]]:
=C/
(4*
xz) ;

A1w[
l,
ndv
[13,
l
]]:
=-C/
(4*
xz) ;

A1w[
l,
ndv
[14,
l
]]:
=-C/
(4*
xz) ;

A1w[
l,
ndv
[15,
l
]]:
=C/
(4*
xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;
 A2v
[l
,ndv
[2,
l
]]:
=k5/
xx ;

A2v
[l
,ndv
[3,
l
]]:
=k5/
xx ;

A2v
[l
,ndv
[5,
l
]]:
=2*
A/y
y ;

A2v
[l
,ndv
[6,
l
]]:
=k5/
zz ;

A2v
[l
,ndv
[7,
l
]]:
=k5/
zz ;

A2u[
l,
ndv
[2,
l
]]:
=-2*
A*(
D/A)
/(2*
xy) ;

A2u[
l,
ndv
[3,
l
]]:
=2*
A*(
D/A)
/(2*
xy) ;

A2w[
l,
ndv
[6,
l
]]:
=-2*
A*(
D/A)
/(2*
yz) ;

A2w[
l,
ndv
[7,
l
]]:
=2*
A*(
D/A)
/(2*
yz) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
B/xx
-2*
k5/
yy-
2*A/
zz-
rho/
tt ;

A3w[
l,
ndv
[2,
l
]]:
=B/
xx ;

A3w[
l,
ndv
[3,
l
]]:
=B/
xx ;

A3w[
l,
ndv
[5,
l
]]:
=2*
k5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=A/
zz ;

A3w[
l,
ndv
[7,
l
]]:
=A/
zz ;

A3u[
l,
ndv
[12,
l
]]:
=C/
(4*
xz) ;

A3u[
l,
ndv
[13,
l
]]:
=-C/
(4*
xz) ;

A3u[
l,
ndv
[14,
l
]]:
=-C/
(4*
xz) ;

A3u[
l,
ndv
[15,
l
]]:
=C/
(4*
xz) ;

A3v
[l
,ndv
[6,
l
]]:
=-2*
k5/
(2*
yz) ;

A3v
[l
,ndv
[7,
l
]]:
=2*
k5/
(2*
yz) ;
 end;

3:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*B/
zz-
rho/
tt ;

A1u[
l,
ndv
[2,
l
]]:
=A/
xx ;

A1u[
l,
ndv
[3,
l
]]:
=A/
xx ;

A1u[
l,
ndv
[4,
l
]]:
=2*
k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=B/
zz ;

A1u[
l,
ndv
[7,
l
]]:
=B/
zz ;

A1v
[l
,ndv
[2,
l
]]:
=2*
k5/
(2*
xy) ;

A1v
[l
,ndv
[3,
l
]]:
=-2*
k5/
(2*
xy) ;

A1w[
l,
ndv
[12,
l
]]:
=C/
(4*
xz) ;

A1w[
l,
ndv
[13,
l
]]:
=-C/
(4*
xz) ;

A1w[
l,
ndv
[14,
l
]]:
=-C/
(4*
xz) ;

A1w[
l,
ndv
[15,
l
]]:
=C/
(4*
xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;

A2v
[l
,ndv
[2,
l
]]:
=k5/
xx ;

A2v
[l
,ndv
[3,
l
]]:
=k5/
xx ;

A2v
[l
,ndv
[4,
l
]]:
=2*
A/y
y ;

A2v
[l
,ndv
[6,
l
]]:
=k5/
zz ;

A2v
[l
,ndv
[7,
l
]]:
=k5/
zz ;
 A2u[
l,
ndv
[2,
l
]]:
=2*
A*(
D/A)
/(2*
xy) ;

A2u[
l,
ndv
[3,
l
]]:
=-2*
A*(
D/A)
/(2*
xy) ;

A2w[
l,
ndv
[6,
l
]]:
=2*
A*(
D/A)
/(2*
yz) ;

A2w[
l,
ndv
[7,
l
]]:
=-2*
A*(
D/A)
/(2*
yz) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
B/xx
-2*
k5/
yy-
2*A/
zz-
rho/
tt ;

A3w[
l,
ndv
[2,
l
]]:
=B/
xx ;

A3w[
l,
ndv
[3,
l
]]:
=B/
xx ;

A3w[
l,
ndv
[4,
l
]]:
=2*
k5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=A/
zz ;

A3w[
l,
ndv
[7,
l
]]:
=A/
zz ;

A3u[
l,
ndv
[12,
l
]]:
=C/
(4*
xz) ;

A3u[
l,
ndv
[13,
l
]]:
=-C/
(4*
xz) ;

A3u[
l,
ndv
[14,
l
]]:
=-C/
(4*
xz) ;

A3u[
l,
ndv
[15,
l
]]:
=C/
(4*
xz) ;

A3v
[l
,ndv
[6,
l
]]:
=2*
k5/
(2*
yz) ;

A3v
[l
,ndv
[7,
l
]]:
=-2*
k5/
(2*
yz) ;

end;

2:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt;
 A1u[
l,
ndv
[2,
l
]]:
=2*
A/xx ;

A1u[
l,
ndv
[4,
l
]]:
=k5/
yy ;

A1u[
l,
ndv
[5,
l
]]:
=k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=k5/
zz ;

A1u[
l,
ndv
[7,
l
]]:
=k5/
zz ;

A1v
[l
,ndv
[4,
l
]]:
=2*
A*(
D/A)
/(2*
xy) ;

A1v
[l
,ndv
[5,
l
]]:
=-2*
A*(
D/A)
/(2*
xy) ;

A1w[
l,
ndv
[6,
l
]]:
=2*
A*(
D/A)
/(2*
xz) ;

A1w[
l,
ndv
[7,
l
]]:
=-2*
A*(
D/A)
/(2*
xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*B/
zz-
rho/
tt;

A2v
[l
,ndv
[2,
l
]]:
=2*
k5/
xx ;

A2v
[l
,ndv
[4,
l
]]:
=A/
yy ;

A2v
[l
,ndv
[5,
l
]]:
=A/
yy ;

A2v
[l
,ndv
[6,
l
]]:
=B/
zz ;

A2v
[l
,ndv
[7,
l
]]:
=B/
zz ;

A2u[
l,
ndv
[4,
l
]]:
=2*
k5/
(2*
xy) ;

A2u[
l,
ndv
[5,
l
]]:
=-2*
k5/
(2*
xy) ;

A2w[
l,
ndv
[16,
l
]]:
=C/
(4*
yz) ;

A2w[
l,
ndv
[17,
l
]]:
=-C/
(4*
yz) ;

A2w[
l,
ndv
[18,
l
]]:
=-C/
(4*
yz) ;

A2w[
l,
ndv
[19,
l
]]:
=C/
(4*
yz) ;
/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*B/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[2,
l
]]:
=2*
k5/
xx ;

A3w[
l,
ndv
[4,
l
]]:
=B/
yy ;

A3w[
l,
ndv
[5,
l
]]:
=B/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=A/
zz ;

A3w[
l,
ndv
[7,
l
]]:
=A/
zz ;

A3u[
l,
ndv
[6,
l
]]:
=2*
k5/
(2*
xz) ;

A3u[
l,
ndv
[7,
l
]]:
=-2*
k5/
(2*
xz) ;

A3v
[l
,ndv
[16,
l
]]:
=C/
(4*
yz) ;

A3v
[l
,ndv
[17,
l
]]:
=-C/
(4*
yz) ;

A3v
[l
,ndv
[18,
l
]]:
=-C/
(4*
yz) ;

A3v
[l
,ndv
[19,
l
]]:
=C/
(4*
yz) ;

end;

4:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt;

A1u[
l,
ndv
[3,
l
]]:
=2*
A/xx ;

A1u[
l,
ndv
[4,
l
]]:
=k5/
yy ;

A1u[
l,
ndv
[5,
l
]]:
=k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=k5/
zz ;

A1u[
l,
ndv
[7,
l
]]:
=k5/
zz ;
 A1v
[l
,ndv
[4,
l
]]:
=-2*
A*(
D/A)
/(2*
xy) ;

A1v
[l
,ndv
[5,
l
]]:
=2*
A*(
D/A)
/(2*
xy) ;

A1w[
l,
ndv
[6,
l
]]:
=-2*
A*(
D/A)
/(2*
xz) ;

A1w[
l,
ndv
[7,
l
]]:
=2*
A*(
D/A)
/(2*
xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*B/
zz-
rho/
tt;

A2v
[l
,ndv
[3,
l
]]:
=2*
k5/
xx ;

A2v
[l
,ndv
[4,
l
]]:
=A/
yy ;

A2v
[l
,ndv
[5,
l
]]:
=A/
yy ;

A2v
[l
,ndv
[6,
l
]]:
=B/
zz ;

A2v
[l
,ndv
[7,
l
]]:
=B/
zz ;

A2u[
l,
ndv
[4,
l
]]:
=-2*
k5/
(2*
xy) ;

A2u[
l,
ndv
[5,
l
]]:
=2*
k5/
(2*
xy) ;

A2w[
l,
ndv
[16,
l
]]:
=C/
(4*
yz) ;

A2w[
l,
ndv
[17,
l
]]:
=-C/
(4*
yz) ;

A2w[
l,
ndv
[18,
l
]]:
=-C/
(4*
yz) ;

A2w[
l,
ndv
[19,
l
]]:
=C/
(4*
yz) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*B/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[3,
l
]]:
=2*
k5/
xx ;

A3w[
l,
ndv
[4,
l
]]:
=B/
yy ;

A3w[
l,
ndv
[5,
l
]]:
=B/
yy ;
 A3w[
l,
ndv
[6,
l
]]:
=A/
zz ;

A3w[
l,
ndv
[7,
l
]]:
=A/
zz ;

A3u[
l,
ndv
[6,
l
]]:
=-2*
k5/
(2*
xz) ;

A3u[
l,
ndv
[7,
l
]]:
=2*
k5/
(2*
xz) ;

A3v
[l
,ndv
[16,
l
]]:
=C/
(4*
yz) ;

A3v
[l
,ndv
[17,
l
]]:
=-C/
(4*
yz) ;

A3v
[l
,ndv
[18,
l
]]:
=-C/
(4*
yz) ;

A3v
[l
,ndv
[19,
l
]]:
=C/
(4*
yz) ;

end;

6:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*B/
yy-
2*k5/
zz-
rho/
tt;

A1u[
l,
ndv
[2,
l
]]:
=A/
xx ;

A1u[
l,
ndv
[3,
l
]]:
=A/
xx ;

A1u[
l,
ndv
[4,
l
]]:
=B/
yy ;

A1u[
l,
ndv
[5,
l
]]:
=B/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=2*
k5/
zz ;

A1v
[l
,ndv
[8,
l
]]:
=C/
(4*
xy) ;

A1v
[l
,ndv
[9,
l
]]:
=-C/
(4*
xy) ;

A1v
[l
,ndv
[10,
l
]]:
=-C/
(4*
xy) ;

A1v
[l
,ndv
[11,
l
]]:
=C/
(4*
xy) ;
 A1w[
l,
ndv
[2,
l
]]:
=2*
k5/
(2*
xz) ;

A1w[
l,
ndv
[3,
l
]]:
=-2*
k5/
(2*
xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
B/xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;

A2v
[l
,ndv
[2,
l
]]:
=B/
xx ;

A2v
[l
,ndv
[3,
l
]]:
=B/
xx ;

A2v
[l
,ndv
[4,
l
]]:
=A/
yy ;

A2v
[l
,ndv
[5,
l
]]:
=A/
yy ;

A2v
[l
,ndv
[6,
l
]]:
=2*
k5/
zz ;

A2u[
l,
ndv
[8,
l
]]:
=C/
(4*
xy) ;

A2u[
l,
ndv
[9,
l
]]:
=-C/
(4*
xy) ;

A2u[
l,
ndv
[10,
l
]]:
=-C/
(4*
xy) ;

A2u[
l,
ndv
[11,
l
]]:
=C/
(4*
xy) ;

A2w[
l,
ndv
[4,
l
]]:
=2*
k5/
(2*
yz) ;

A2w[
l,
ndv
[5,
l
]]:
=-2*
k5/
(2*
yz) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[2,
l
]]:
=k5/
xx ;

A3w[
l,
ndv
[3,
l
]]:
=k5/
xx ;

A3w[
l,
ndv
[4,
l
]]:
=k5/
yy ;

A3w[
l,
ndv
[5,
l
]]:
=k5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=2*
A/zz ;

A3u[
l,
ndv
[2,
l
]]:
=2*
A*(
D/A)
/(2*
xz) ;
 A3u[
l,
ndv
[3,
l
]]:
=-2*
A*(
D/A)
/(2*
xz) ;

A3v
[l
,ndv
[4,
l
]]:
=2*
A*(
D/A)
/(2*
yz) ;

A3v
[l
,ndv
[5,
l
]]:
=-2*
A*(
D/A)
/(2*
yz) ;

end;

7:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt;

A1u[
l,
ndv
[3,
l
]]:
=2*
A/xx ;

A1u[
l,
ndv
[5,
l
]]:
=2*
k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
= k5/
zz ;

A1u[
l,
ndv
[7,
l
]]:
= k5/
zz ;

A1v
[l
,ndv
[1,
l
]]:
=-2*
k5/
(xy
)-2*
A*(
D/A)
/(xy
) ;

A1v
[l
,ndv
[3,
l
]]:
=2*
k5/
(xy
) ;

A1v
[l
,ndv
[5,
l
]]:
=2*
A*(
D/A)
/(xy
) ;

A1w[
l,
ndv
[6,
l
]]:
=-2*
A*(
D/A)
/(2*
xz) ;

A1w[
l,
ndv
[7,
l
]]:
=2*
A*(
D/A)
/(2*
xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;

A2v
[l
,ndv
[3,
l
]]:
=2*
k5/
xx ;

A2v
[l
,ndv
[5,
l
]]:
=2*
A/y
y ;

A2v
[l
,ndv
[6,
l
]]:
= k5/
zz ;

A2v
[l
,ndv
[7,
l
]]:
= k5/
zz ;
 A2u[
l,
ndv
[1,
l
]]:
=-2*
A*(
D/A)
/(xy
)-2*
k5/
(xy
) ;

A2u[
l,
ndv
[3,
l
]]:
=2*
A*(
D/A)
/(xy
) ;

A2u[
l,
ndv
[5,
l
]]:
=2*
k5/
(xy
) ;

A2w[
l,
ndv
[6,
l
]]:
=-2*
A*(
D/A)
/(2*
yz) ;

A2w[
l,
ndv
[7,
l
]]:
=2*
A*(
D/A)
/(2*
yz) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[3,
l
]]:
=2*
k5/
xx ;

A3w[
l,
ndv
[5,
l
]]:
=2*
k5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
= A/
zz ;

A3w[
l,
ndv
[7,
l
]]:
= A/
zz ;

A3u[
l,
ndv
[6,
l
]]:
=-2*
k5/
(2*
xz) ;

A3u[
l,
ndv
[7,
l
]]:
=2*
k5/
(2*
xz) ;

A3v
[l
,ndv
[6,
l
]]:
=-2*
k5/
(2*
yz) ;

A3v
[l
,ndv
[7,
l
]]:
=2*
k5/
(2*
yz) ;

end;

8:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt;

A1u[
l,
ndv
[2,
l
]]:
=2*
A/xx ;

A1u[
l,
ndv
[5,
l
]]:
=2*
k5/
yy ;
 A1u[
l,
ndv
[6,
l
]]:
= k5/
zz ;

A1u[
l,
ndv
[7,
l
]]:
= k5/
zz ;

A1v
[l
,ndv
[1,
l
]]:
=2*
k5/
(xy
)+2*
A*(
D/A)
/(xy
) ;

A1v
[l
,ndv
[2,
l
]]:
=-2*
k5/
(xy
) ;

A1v
[l
,ndv
[5,
l
]]:
=-2*
A*(
D/A)
/(xy
) ;

A1w[
l,
ndv
[6,
l
]]:
=2*
A*(
D/A)
/(2*
xz) ;

A1w[
l,
ndv
[7,
l
]]:
=-2*
A*(
D/A)
/(2*
xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt ;

A2v
[l
,ndv
[2,
l
]]:
=2*
k5/
xx ;

A2v
[l
,ndv
[5,
l
]]:
=2*
A/y
y ;

A2v
[l
,ndv
[6,
l
]]:
= k5/
zz ;

A2v
[l
,ndv
[7,
l
]]:
= k5/
zz ;

A2u[
l,
ndv
[1,
l
]]:
=2*
A*(
D/A)
/(xy
)+2*
k5/
(xy
) ;

A2u[
l,
ndv
[2,
l
]]:
=-2*
A*(
D/A)
/(xy
) ;

A2u[
l,
ndv
[5,
l
]]:
=-2*
k5/
(xy
) ;

A2w[
l,
ndv
[6,
l
]]:
=-2*
A*(
D/A)
/(2*
yz) ;

A2w[
l,
ndv
[7,
l
]]:
=2*
A*(
D/A)
/(2*
yz) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[2,
l
]]:
=2*
k5/
xx ;

A3w[
l,
ndv
[5,
l
]]:
=2*
k5/
yy ;
 A3w[
l,
ndv
[6,
l
]]:
= A/
zz ;

A3w[
l,
ndv
[7,
l
]]:
= A/
zz ;

A3u[
l,
ndv
[6,
l
]]:
=2*
k5/
(2*
xz) ;

A3u[
l,
ndv
[7,
l
]]:
=-2*
k5/
(2*
xz) ;

A3v
[l
,ndv
[6,
l
]]:
=-2*
k5/
(2*
yz) ;

A3v
[l
,ndv
[7,
l
]]:
=2*
k5/
(2*
yz) ;

end;

9:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt;

A1u[
l,
ndv
[2,
l
]]:
=2*
A/xx ;

A1u[
l,
ndv
[4,
l
]]:
=2*
k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
= k5/
zz ;

A1u[
l,
ndv
[7,
l
]]:
= k5/
zz ;

A1v
[l
,ndv
[1,
l
]]:
=-2*
k5/
(xy
)-2*
A*(
D/A)
/(xy
) ;

A1v
[l
,ndv
[2,
l
]]:
=2*
k5/
(xy
) ;

A1v
[l
,ndv
[4,
l
]]:
=2*
A*(
D/A)
/(xy
) ;

A1w[
l,
ndv
[6,
l
]]:
=2*
A*(
D/A)
/(2*
xz) ;

A1w[
l,
ndv
[7,
l
]]:
=-2*
A*(
D/A)
/(2*
xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;
 A2v
[l
,ndv
[2,
l
]]:
=2*
k5/
xx ;

A2v
[l
,ndv
[4,
l
]]:
=2*
A/y
y ;

A2v
[l
,ndv
[6,
l
]]:
= k5/
zz ;

A2v
[l
,ndv
[7,
l
]]:
= k5/
zz ;

A2u[
l,
ndv
[1,
l
]]:
=-2*
A*(
D/A)
/(xy
)-2*
k5/
(xy
) ;

A2u[
l,
ndv
[2,
l
]]:
=2*
A*(
D/A)
/(xy
) ;

A2u[
l,
ndv
[4,
l
]]:
=2*
k5/
(xy
) ;

A2w[
l,
ndv
[6,
l
]]:
=2*
A*(
D/A)
/(2*
yz) ;

A2w[
l,
ndv
[7,
l
]]:
=-2*
A*(
D/A)
/(2*
yz) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[2,
l
]]:
=2*
k5/
xx ;

A3w[
l,
ndv
[4,
l
]]:
=2*
k5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
= A/
zz ;

A3w[
l,
ndv
[7,
l
]]:
= A/
zz ;

A3u[
l,
ndv
[6,
l
]]:
=2*
k5/
(2*
xz) ;

A3u[
l,
ndv
[7,
l
]]:
=-2*
k5/
(2*
xz) ;

A3v
[l
,ndv
[6,
l
]]:
=2*
k5/
(2*
yz) ;

A3v
[l
,ndv
[7,
l
]]:
=-2*
k5/
(2*
yz) ;

end;

10:
begi
n
/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt;

A1u[
l,
ndv
[3,
l
]]:
=2*
A/xx ;

A1u[
l,
ndv
[4,
l
]]:
=2*
k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
= k5/
zz ;

A1u[
l,
ndv
[7,
l
]]:
= k5/
zz ;

A1v
[l
,ndv
[1,
l
]]:
=2*
k5/
(xy
)+2*
A*(
D/A)
/(xy
) ;

A1v
[l
,ndv
[3,
l
]]:
=-2*
k5/
(xy
) ;

A1v
[l
,ndv
[4,
l
]]:
=-2*
A*(
D/A)
/(xy
) ;

A1w[
l,
ndv
[6,
l
]]:
=-2*
A*(
D/A)
/(2*
xz) ;

A1w[
l,
ndv
[7,
l
]]:
=2*
A*(
D/A)
/(2*
xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt ;

A2v
[l
,ndv
[3,
l
]]:
=2*
k5/
xx ;

A2v
[l
,ndv
[4,
l
]]:
=2*
A/y
y ;

A2v
[l
,ndv
[6,
l
]]:
= k5/
zz ;

A2v
[l
,ndv
[7,
l
]]:
= k5/
zz ;

A2u[
l,
ndv
[1,
l
]]:
=2*
A*(
D/A)
/(xy
)+2*
k5/
(xy
) ;

A2u[
l,
ndv
[3,
l
]]:
=-2*
A*(
D/A)
/(xy
) ;

A2u[
l,
ndv
[4,
l
]]:
=-2*
k5/
(xy
) ;

A2w[
l,
ndv
[6,
l
]]:
=2*
A*(
D/A)
/(2*
yz) ;

A2w[
l,
ndv
[7,
l
]]:
=-2*
A*(
D/A)
/(2*
yz) ;
/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[3,
l
]]:
=2*
k5/
xx ;

A3w[
l,
ndv
[4,
l
]]:
=2*
k5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
= A/
zz ;

A3w[
l,
ndv
[7,
l
]]:
= A/
zz ;

A3u[
l,
ndv
[6,
l
]]:
=-2*
k5/
(2*
xz) ;

A3u[
l,
ndv
[7,
l
]]:
=2*
k5/
(2*
xz) ;

A3v
[l
,ndv
[6,
l
]]:
=2*
k5/
(2*
yz) ;

A3v
[l
,ndv
[7,
l
]]:
=-2*
k5/
(2*
yz) ;

end;

15:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt;

A1u[
l,
ndv
[2,
l
]]:
= A/
xx ;

A1u[
l,
ndv
[3,
l
]]:
= A/
xx ;

A1u[
l,
ndv
[5,
l
]]:
=2*
k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=2*
k5/
zz ;

A1v
[l
,ndv
[2,
l
]]:
=-2*
k5/
(2*
xy) ;

A1v
[l
,ndv
[3,
l
]]:
=2*
k5/
(2*
xy) ;

A1w[
l,
ndv
[2,
l
]]:
=2*
k5/
(2*
xz) ;
 A1w[
l,
ndv
[3,
l
]]:
=-2*
k5/
(2*
xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;

A2v
[l
,ndv
[2,
l
]]:
= k5/
xx ;

A2v
[l
,ndv
[3,
l
]]:
= k5/
xx ;

A2v
[l
,ndv
[5,
l
]]:
=2*
A/y
y ;

A2v
[l
,ndv
[6,
l
]]:
=2*
k5/
zz ;

A2u[
l,
ndv
[2,
l
]]:
=-2*
A*(
D/A)
/(2*
xy) ;

A2u[
l,
ndv
[3,
l
]]:
=2*
A*(
D/A)
/(2*
xy) ;

A2w[
l,
ndv
[1,
l
]]:
=2*
k5/
(yz)+2*
A*(
D/A)
/(y
z) ;

A2w[
l,
ndv
[5,
l
]]:
=-2*
k5/
(yz) ;

A2w[
l,
ndv
[6,
l
]]:
=-2*
A*(
D/A)
/(y
z) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[2,
l
]]:
= k5/
xx ;

A3w[
l,
ndv
[3,
l
]]:
= k5/
xx ;

A3w[
l,
ndv
[5,
l
]]:
=2*
k5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=2*
A/zz ;

A3u[
l,
ndv
[2,
l
]]:
=2*
A*(
D/A)
/(2*
xz) ;

A3u[
l,
ndv
[3,
l
]]:
=-2*
A*(
D/A)
/(2*
xz) ;

A3v
[l
,ndv
[1,
l
]]:
=2*
A*(
D/A)
/(y
z)+2*
k5/
(yz) ;

A3v
[l
,ndv
[5,
l
]]:
=-2*
A*(
D/A)
/(y
z) ;
 A3v
[l
,ndv
[6,
l
]]:
=-2*
k5/
(yz) ;

end;

17:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt;

A1u[
l,
ndv
[2,
l
]]:
= A/
xx ;

A1u[
l,
ndv
[3,
l
]]:
= A/
xx ;

A1u[
l,
ndv
[4,
l
]]:
=2*
k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=2*
k5/
zz ;

A1v
[l
,ndv
[2,
l
]]:
=2*
k5/
(2*
xy) ;

A1v
[l
,ndv
[3,
l
]]:
=-2*
k5/
(2*
xy) ;

A1w[
l,
ndv
[2,
l
]]:
=2*
k5/
(2*
xz) ;

A1w[
l,
ndv
[3,
l
]]:
=-2*
k5/
(2*
xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;

A2v
[l
,ndv
[2,
l
]]:
= k5/
xx ;

A2v
[l
,ndv
[3,
l
]]:
= k5/
xx ;

A2v
[l
,ndv
[4,
l
]]:
=2*
A/y
y ;

A2v
[l
,ndv
[6,
l
]]:
=2*
k5/
zz ;

A2u[
l,
ndv
[2,
l
]]:
=2*
A*(
D/A)
/(2*
xy) ;

A2u[
l,
ndv
[3,
l
]]:
=-2*
A*(
D/A)
/(2*
xy) ;
 A2w[
l,
ndv
[1,
l
]]:
=-2*
k5/
(yz)-2*
A*(
D/A)
/(y
z) ;

A2w[
l,
ndv
[4,
l
]]:
=2*
k5/
(yz) ;

A2w[
l,
ndv
[6,
l
]]:
=2*
A*(
D/A)
/(y
z) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[2,
l
]]:
= k5/
xx ;

A3w[
l,
ndv
[3,
l
]]:
= k5/
xx ;

A3w[
l,
ndv
[4,
l
]]:
=2*
k5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=2*
A/zz ;

A3u[
l,
ndv
[2,
l
]]:
=2*
A*(
D/A)
/(2*
xz) ;

A3u[
l,
ndv
[3,
l
]]:
=-2*
A*(
D/A)
/(2*
xz) ;

A3v
[l
,ndv
[1,
l
]]:
=-2*
A*(
D/A)
/(y
z)-2*
k5/
(yz) ;

A3v
[l
,ndv
[4,
l
]]:
=2*
A*(
D/A)
/(y
z) ;

A3v
[l
,ndv
[6,
l
]]:
=2*
k5/
(yz) ;

end;

16:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt ;

A1u[
l,
ndv
[2,
l
]]:
=2*
A/xx ;

A1u[
l,
ndv
[4,
l
]]:
= k5/
yy ;

A1u[
l,
ndv
[5,
l
]]:
= k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=2*
k5/
zz ;
 A1v
[l
,ndv
[4,
l
]]:
=2*
A*(
D/A)
/(2*
xy) ;

A1v
[l
,ndv
[5,
l
]]:
=-2*
A*(
D/A)
/(2*
xy) ;

A1w[
l,
ndv
[1,
l
]]:
=-2*
k5/
(xz)-2*
A*(
D/A)
/(xz) ;

A1w[
l,
ndv
[2,
l
]]:
=2*
k5/
(xz) ;

A1w[
l,
ndv
[6,
l
]]:
=2*
A*(
D/A)
/(xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;

A2v
[l
,ndv
[2,
l
]]:
=2*
k5/
xx ;

A2v
[l
,ndv
[4,
l
]]:
= A/
yy ;

A2v
[l
,ndv
[5,
l
]]:
= A/
yy ;

A2v
[l
,ndv
[6,
l
]]:
=2*
k5/
zz ;

A2u[
l,
ndv
[4,
l
]]:
=2*
k5/
(2*
xy) ;

A2u[
l,
ndv
[5,
l
]]:
=-2*
k5/
(2*
xy) ;

A2w[
l,
ndv
[4,
l
]]:
=2*
k5/
(2*
yz) ;

A2w[
l,
ndv
[5,
l
]]:
=-2*
k5/
(2*
yz) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[2,
l
]]:
=2*
k5/
xx ;

A3w[
l,
ndv
[4,
l
]]:
= k5/
yy ;

A3w[
l,
ndv
[5,
l
]]:
= k5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=2*
A/zz ;

A3u[
l,
ndv
[1,
l
]]:
=-2*
A*(
D/A)
/(xz)-2*
k5/
(xz) ;
 A3u[
l,
ndv
[2,
l
]]:
=2*
A*(
D/A)
/(xz) ;

A3u[
l,
ndv
[6,
l
]]:
=2*
k5/
(xz) ;

A3v
[l
,ndv
[4,
l
]]:
=2*
A*(
D/A)
/(2*
yz) ;

A3v
[l
,ndv
[5,
l
]]:
=-2*
A*(
D/A)
/(2*
yz) ;

end;

18:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt ;

A1u[
l,
ndv
[3,
l
]]:
=2*
A/xx ;

A1u[
l,
ndv
[4,
l
]]:
= k5/
yy ;

A1u[
l,
ndv
[5,
l
]]:
= k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=2*
k5/
zz ;

A1v
[l
,ndv
[4,
l
]]:
=-2*
A*(
D/A)
/(2*
xy) ;

A1v
[l
,ndv
[5,
l
]]:
=2*
A*(
D/A)
/(2*
xy) ;

A1w[
l,
ndv
[1,
l
]]:
=2*
k5/
(xz)+2*
A*(
D/A)
/(xz) ;

A1w[
l,
ndv
[3,
l
]]:
=-2*
k5/
(xz) ;

A1w[
l,
ndv
[6,
l
]]:
=-2*
A*(
D/A)
/(xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;

A2v
[l
,ndv
[3,
l
]]:
=2*
k5/
xx ;

A2v
[l
,ndv
[4,
l
]]:
= A/
yy ;

A2v
[l
,ndv
[5,
l
]]:
= A/
yy ;
 A2v
[l
,ndv
[6,
l
]]:
=2*
k5/
zz ;

A2u[
l,
ndv
[4,
l
]]:
=-2*
k5/
(2*
xy) ;

A2u[
l,
ndv
[5,
l
]]:
=2*
k5/
(2*
xy) ;

A2w[
l,
ndv
[4,
l
]]:
=2*
k5/
(2*
yz) ;

A2w[
l,
ndv
[5,
l
]]:
=-2*
k5/
(2*
yz) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[3,
l
]]:
=2*
k5/
xx ;

A3w[
l,
ndv
[4,
l
]]:
= k5/
yy ;

A3w[
l,
ndv
[5,
l
]]:
= k5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=2*
A/zz ;

A3u[
l,
ndv
[1,
l
]]:
=2*
A*(
D/A)
/(xz)+2*
k5/
(xz) ;

A3u[
l,
ndv
[3,
l
]]:
=-2*
A*(
D/A)
/(xz) ;

A3u[
l,
ndv
[6,
l
]]:
=-2*
k5/
(xz) ;

A3v
[l
,ndv
[4,
l
]]:
=2*
A*(
D/A)
/(2*
yz) ;

A3v
[l
,ndv
[5,
l
]]:
=-2*
A*(
D/A)
/(2*
yz) ;

end;

23:
begi
n

F2[
l]:
=F2[
l]+2*
Fs3/
(dy
);

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt;
 A1u[
l,
ndv
[3,
l
]]:
=2*
A/xx ;

A1u[
l,
ndv
[5,
l
]]:
=2*
k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=2*
k5/
zz ;

A1v
[l
,ndv
[1,
l
]]:
=-2*
k5/
(xy
)-2*
A*(
D/A)
/(xy
) ;

A1v
[l
,ndv
[3,
l
]]:
=2*
k5/
(xy
) ;

A1v
[l
,ndv
[5,
l
]]:
=2*
A*(
D/A)
/(xy
) ;

A1w[
l,
ndv
[1,
l
]]:
=2*
k5/
(xz)+2*
A*(
D/A)
/(xz) ;

A1w[
l,
ndv
[3,
l
]]:
=-2*
k5/
(xz) ;

A1w[
l,
ndv
[6,
l
]]:
=-2*
A*(
D/A)
/(xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;

A2v
[l
,ndv
[3,
l
]]:
=2*
k5/
xx ;

A2v
[l
,ndv
[5,
l
]]:
=2*
A/y
y ;

A2v
[l
,ndv
[6,
l
]]:
=2*
k5/
zz ;

A2u[
l,
ndv
[1,
l
]]:
=-2*
A*(
D/A)
/(xy
)-2*
k5/
(xy
) ;

A2u[
l,
ndv
[3,
l
]]:
=2*
A*(
D/A)
/(xy
) ;

A2u[
l,
ndv
[5,
l
]]:
=2*
k5/
(xy
) ;

A2w[
l,
ndv
[1,
l
]]:
=2*
k5/
(yz)+2*
A*(
D/A)
/(y
z) ;

A2w[
l,
ndv
[5,
l
]]:
=-2*
k5/
(yz) ;

A2w[
l,
ndv
[6,
l
]]:
=-2*
A*(
D/A)
/(y
z) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;
 A3w[
l,
ndv
[3,
l
]]:
=k5/
xx ;

A3w[
l,
ndv
[5,
l
]]:
=K5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=A/
zz ;

A3u[
l,
ndv
[1,
l
]]:
=2*
A*(
D/A)
/(xz)+2*
k5/
(xz) ;

A3u[
l,
ndv
[3,
l
]]:
=-2*
A*(
D/A)
/(xz) ;

A3u[
l,
ndv
[6,
l
]]:
=-2*
k5/
(xz) ;

A3v
[l
,ndv
[1,
l
]]:
=2*
A*(
D/A)
/(y
z)+2*
k5/
(yz) ;

A3v
[l
,ndv
[5,
l
]]:
=-2*
A*(
D/A)
/(y
z) ;

A3v
[l
,ndv
[6,
l
]]:
=-2*
k5/
(yz) ;

end;

24:
begi
n

F2[
l]:
=F2[
l]+2*
Fs4/
(dy
);

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt;

A1u[
l,
ndv
[2,
l
]]:
=2*
A/xx ;

A1u[
l,
ndv
[5,
l
]]:
=2*
k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=2*
k5/
zz ;

A1v
[l
,ndv
[1,
l
]]:
=2*
k5/
(xy
)+2*
A*(
D/A)
/(xy
) ;

A1v
[l
,ndv
[2,
l
]]:
=-2*
k5/
(xy
) ;

A1v
[l
,ndv
[5,
l
]]:
=-2*
A*(
D/A)
/(xy
) ;

A1w[
l,
ndv
[1,
l
]]:
=-2*
k5/
(xz)-2*
A*(
D/A)
/(xz) ;
 A1w[
l,
ndv
[2,
l
]]:
=2*
k5/
(xz) ;

A1w[
l,
ndv
[6,
l
]]:
=2*
A*(
D/A)
/(xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;

A2v
[l
,ndv
[2,
l
]]:
=2*
k5/
xx ;

A2v
[l
,ndv
[5,
l
]]:
=2*
A/y
y ;

A2v
[l
,ndv
[6,
l
]]:
=2*
k5/
zz ;

A2u[
l,
ndv
[1,
l
]]:
=2*
A*(
D/A)
/(xy
)+2*
k5/
(xy
) ;

A2u[
l,
ndv
[2,
l
]]:
=-2*
A*(
D/A)
/(xy
) ;

A2u[
l,
ndv
[5,
l
]]:
=-2*
k5/
(xy
) ;

A2w[
l,
ndv
[1,
l
]]:
=2*
k5/
(yz)+2*
A*(
D/A)
/(y
z) ;

A2w[
l,
ndv
[5,
l
]]:
=-2*
k5/
(yz) ;

A2w[
l,
ndv
[6,
l
]]:
=-2*
A*(
D/A)
/(y
z) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[2,
l
]]:
=k5/
xx ;

A3w[
l,
ndv
[5,
l
]]:
=K5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=A/
zz ;

A3u[
l,
ndv
[1,
l
]]:
=-2*
A*(
D/A)
/(xz)-2*
k5/
(xz) ;

A3u[
l,
ndv
[2,
l
]]:
=2*
A*(
D/A)
/(xz) ;

A3u[
l,
ndv
[6,
l
]]:
=2*
k5/
(xz) ;

A3v
[l
,ndv
[1,
l
]]:
=2*
A*(
D/A)
/(y
z)+2*
k5/
(yz) ;
 A3v
[l
,ndv
[5,
l
]]:
=-2*
A*(
D/A)
/(y
z) ;

A3v
[l
,ndv
[6,
l
]]:
=-2*
k5/
(yz) ;

end;

25:
begi
n

F2[
l]:
=F2[
l]-2*
Fs2/
(dy
);

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt;

A1u[
l,
ndv
[2,
l
]]:
=2*
A/xx ;

A1u[
l,
ndv
[4,
l
]]:
=2*
k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=2*
k5/
zz ;

A1v
[l
,ndv
[1,
l
]]:
=-2*
k5/
(xy
)-2*
A*(
D/A)
/(xy
) ;

A1v
[l
,ndv
[2,
l
]]:
=2*
k5/
(xy
) ;

A1v
[l
,ndv
[4,
l
]]:
=2*
A*(
D/A)
/(xy
) ;

A1w[
l,
ndv
[1,
l
]]:
=-2*
k5/
(xz)-2*
A*(
D/A)
/(xz) ;

A1w[
l,
ndv
[2,
l
]]:
=2*
k5/
(xz) ;

A1w[
l,
ndv
[6,
l
]]:
=2*
A*(
D/A)
/(xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;

A2v
[l
,ndv
[2,
l
]]:
=2*
k5/
xx ;

A2v
[l
,ndv
[4,
l
]]:
=2*
A/y
y ;

A2v
[l
,ndv
[6,
l
]]:
=2*
k5/
zz ;

A2u[
l,
ndv
[1,
l
]]:
=-2*
A*(
D/A)
/(xy
)-2*
k5/
(xy
) ;
 A2u[
l,
ndv
[2,
l
]]:
=2*
A*(
D/A)
/(xy
) ;

A2u[
l,
ndv
[4,
l
]]:
=2*
k5/
(xy
) ;

A2w[
l,
ndv
[1,
l
]]:
=-2*
k5/
(yz)-2*
A*(
D/A)
/(y
z) ;

A2w[
l,
ndv
[4,
l
]]:
=2*
k5/
(yz) ;

A2w[
l,
ndv
[6,
l
]]:
=2*
A*(
D/A)
/(y
z) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[2,
l
]]:
=k5/
xx ;

A3w[
l,
ndv
[4,
l
]]:
=k5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=A/
zz ;

A3u[
l,
ndv
[1,
l
]]:
=-2*
A*(
D/A)
/(xz)-2*
k5/
(xz) ;

A3u[
l,
ndv
[2,
l
]]:
=2*
A*(
D/A)
/(xz) ;

A3u[
l,
ndv
[6,
l
]]:
=2*
k5/
(xz) ;

A3v
[l
,ndv
[1,
l
]]:
=-2*
A*(
D/A)
/(y
z)-2*
k5/
(yz) ;

A3v
[l
,ndv
[4,
l
]]:
=2*
A*(
D/A)
/(y
z) ;

A3v
[l
,ndv
[6,
l
]]:
=2*
k5/
(yz) ;

end;

26:
begi
n

F2[
l]:
=F2[
l]-2*
Fs1/
(dy
);

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*k5/
yy-
2*k5/
zz-
rho/
tt;

A1u[
l,
ndv
[3,
l
]]:
=2*
A/xx ;
 A1u[
l,
ndv
[4,
l
]]:
=2*
k5/
yy ;

A1u[
l,
ndv
[6,
l
]]:
=2*
k5/
zz ;

A1v
[l
,ndv
[1,
l
]]:
=2*
k5/
(xy
)+2*
A*(
D/A)
/(xy
) ;

A1v
[l
,ndv
[3,
l
]]:
=-2*
k5/
(xy
) ;

A1v
[l
,ndv
[4,
l
]]:
=-2*
A*(
D/A)
/(xy
) ;

A1w[
l,
ndv
[1,
l
]]:
=2*
k5/
(xz)+2*
A*(
D/A)
/(xz) ;

A1w[
l,
ndv
[3,
l
]]:
=-2*
k5/
(xz) ;

A1w[
l,
ndv
[6,
l
]]:
=-2*
A*(
D/A)
/(xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
k5/
xx-
2*A/
yy-
2*k5/
zz-
rho/
tt;

A2v
[l
,ndv
[3,
l
]]:
=2*
k5/
xx ;

A2v
[l
,ndv
[4,
l
]]:
=2*
A/y
y ;

A2v
[l
,ndv
[6,
l
]]:
=2*
k5/
zz ;

A2u[
l,
ndv
[1,
l
]]:
=2*
A*(
D/A)
/(xy
)+2*
k5/
(xy
) ;

A2u[
l,
ndv
[3,
l
]]:
=-2*
A*(
D/A)
/(xy
) ;

A2u[
l,
ndv
[4,
l
]]:
=-2*
k5/
(xy
) ;

A2w[
l,
ndv
[1,
l
]]:
=-2*
k5/
(yz)-2*
A*(
D/A)
/(y
z) ;

A2w[
l,
ndv
[4,
l
]]:
=2*
k5/
(yz) ;

A2w[
l,
ndv
[6,
l
]]:
=2*
A*(
D/A)
/(y
z) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
k5/
xx-
2*k5/
yy-
2*A/
zz-
rho/
tt;

A3w[
l,
ndv
[3,
l
]]:
=k5/
xx ;
 A3w[
l,
ndv
[4,
l
]]:
=k5/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=A/
zz ;

A3u[
l,
ndv
[1,
l
]]:
=2*
A*(
D/A)
/(xz)+2*
k5/
(xz) ;

A3u[
l,
ndv
[3,
l
]]:
=-2*
A*(
D/A)
/(xz) ;

A3u[
l,
ndv
[6,
l
]]:
=-2*
k5/
(xz) ;

A3v
[l
,ndv
[1,
l
]]:
=-2*
A*(
D/A)
/(y
z)-2*
k5/
(yz) ;

A3v
[l
,ndv
[4,
l
]]:
=2*
A*(
D/A)
/(y
z) ;

A3v
[l
,ndv
[6,
l
]]:
=2*
k5/
(yz) ;

end;

5,
11,
12,
13,
14,
19,
20,
21,
22 /
/k=1

:
begi
n

A1u[
l,
ndv
[1,
l
]]:
={1;
//}-
rho/
tt;

A2v
[l
,ndv
[1,
l
]]:
={1;
//}-
rho/
tt;

A3w[
l,
ndv
[1,
l
]]:
={1;
//}-
rho/
tt;

end;

27:
begi
n

/
/equat
ion1

A1u[
l,
ndv
[1,
l
]]:
=-2*
A/xx-
2*B/
yy-
2*B/
zz-
rho/
tt;

A1u[
l,
ndv
[2,
l
]]:
=A/
xx ;

A1u[
l,
ndv
[3,
l
]]:
=A/
xx ;

A1u[
l,
ndv
[4,
l
]]:
=B/
yy ;

A1u[
l,
ndv
[5,
l
]]:
=B/
yy ;
 A1u[
l,
ndv
[6,
l
]]:
=B/
zz ;

A1u[
l,
ndv
[7,
l
]]:
=B/
zz ;

A1v
[l
,ndv
[8,
l
]]:
=C/
(4*
xy) ;

A1v
[l
,ndv
[9,
l
]]:
=-C/
(4*
xy) ;

A1v
[l
,ndv
[10,
l
]]:
=-C/
(4*
xy) ;

A1v
[l
,ndv
[11,
l
]]:
=C/
(4*
xy) ;

A1w[
l,
ndv
[12,
l
]]:
=C/
(4*
xz) ;

A1w[
l,
ndv
[13,
l
]]:
=-C/
(4*
xz) ;

A1w[
l,
ndv
[14,
l
]]:
=-C/
(4*
xz) ;

A1w[
l,
ndv
[15,
l
]]:
=C/
(4*
xz) ;

/
/equat
ion2

A2v
[l
,ndv
[1,
l
]]:
=-2*
B/xx-
2*A/
yy-
2*B/
zz-
rho/
tt;

A2v
[l
,ndv
[2,
l
]]:
=B/
xx ;

A2v
[l
,ndv
[3,
l
]]:
=B/
xx ;

A2v
[l
,ndv
[4,
l
]]:
=A/
yy ;

A2v
[l
,ndv
[5,
l
]]:
=A/
yy ;

A2v
[l
,ndv
[6,
l
]]:
=B/
zz ;

A2v
[l
,ndv
[7,
l
]]:
=B/
zz ;

A2u[
l,
ndv
[8,
l
]]:
=C/
(4*
xy) ;

A2u[
l,
ndv
[9,
l
]]:
=-C/
(4*
xy) ;

A2u[
l,
ndv
[10,
l
]]:
=-C/
(4*
xy) ;

A2u[
l,
ndv
[11,
l
]]:
=C/
(4*
xy) ;
 A2w[
l,
ndv
[16,
l
]]:
=C/
(4*
yz) ;

A2w[
l,
ndv
[17,
l
]]:
=-C/
(4*
yz) ;

A2w[
l,
ndv
[18,
l
]]:
=-C/
(4*
yz) ;

A2w[
l,
ndv
[19,
l
]]:
=C/
(4*
yz) ;

/
/equat
ion3

A3w[
l,
ndv
[1,
l
]]:
=-2*
B/xx
-2*
B/y
y-2*
A/zz-
rho/
tt;

A3w[
l,
ndv
[2,
l
]]:
=B/
xx ;

A3w[
l,
ndv
[3,
l
]]:
=B/
xx ;

A3w[
l,
ndv
[4,
l
]]:
=B/
yy ;

A3w[
l,
ndv
[5,
l
]]:
=B/
yy ;

A3w[
l,
ndv
[6,
l
]]:
=A/
zz ;

A3w[
l,
ndv
[7,
l
]]:
=A/
zz ;

A3u[
l,
ndv
[12,
l
]]:
=C/
(4*
xz) ;

A3u[
l,
ndv
[13,
l
]]:
=-C/
(4*
xz) ;

A3u[
l,
ndv
[14,
l
]]:
=-C/
(4*
xz) ;

A3u[
l,
ndv
[15,
l
]]:
=C/
(4*
xz) ;

A3v
[l
,ndv
[16,
l
]]:
=C/
(4*
yz) ;

A3v
[l
,ndv
[17,
l
]]:
=-C/
(4*
yz) ;

A3v
[l
,ndv
[18,
l
]]:
=-C/
(4*
yz) ;

A3v
[l
,ndv
[19,
l
]]:
=C/
(4*
yz) ;

end;

end;

End;

end;
/
///
///
///Pr
ocedur
esetFonct
ionscal
cul
mat
ri
ciel/
///
///
///

pr
ocedur
eident
it
y;

begi
n

f
ori
:=1t
onddo

f
orj
:=1t
onddo

i
fi=j
thenI
d[i
,
j]:
=1

el
seI
d[i
,
j]:
=0;

end;

Funct
ioni
nv(
varY:
mat
2dr
e):
mat
2dr
e;/
/Inv
ersed'
unemat
ri
ce

v
ar

p,
q,
pass:
real ;

l
i :
int
eger;

X :
mat
2dr
e;

begi
n

i
dent
it
y;

X:
=Y;

/
/Di
agonal
i
sat
ion

f
ork:
=1t
onddo

begi
n

p:
=X[
k,
k];

l
i:
=k ;

i
fX[
k,
k]=0t
hen /
/Casdepi
votnul

begi
n
f
ori
:=k+1t
onddo

i
fX[
i,
k]<>0t
hen

begi
n

p:
=X[
i,
k];

l
i:
=i ;

br
eak;

end;

i
fli
<>kt
hen

begi
n

f
orj
:
=1t
onddo

begi
n

pass :
=X[
k,
j];

X[
k,
j]:
=X[
li
,
j];

X[
li
,
j]:
=pass ;

pass :
=Id[
k,
j];

I
d[k,
j
]:=I
d[l
i
,j
];

I
d[l
i
,j
]:=pass ;

end;

end;

end;

i
fX[
k,
k]=0t
hen

begi
n

wr
it
eln(
'Er
reur
:mat
ri
cenon-
inv
ersi
ble'
);

br
eak;
 end;

f
ori
:=1t
onddo

i
fi<>kt
hen

begi
n

q:
=X[
i,
k];

f
orj
:=1t
onddo

i
f(p<>0)and(
q<>0)t
hen

begi
n

X[
i,
j
]:=X[
i,
j
]-X[
k,
j]
*q/
p;

I
d[i
,
j]:
=Id[
i,
j
]-I
d[k,
j
]*q/
p;

end;

end;

end;

/
/Id-
fi
cat
ion

f
ori
:=1t
onddo

begi
n

p:
=X[
i,
i
];

f
orj
:=1t
onddo