1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
|
%jose samuel perez huerta
%mayo 2007
%paladium lattice in hydrogen
input featpost3Dplus2D;
% Under the adsorption process,
% the hydrogen molecules are retained in the surface of the
% metal forming a gaseous layer. Nevertheless, the hydrogen molecule is
% too great to get up in the structure of a solid. When the amount of
% hydrogen is increased, the palladium dissociates molecules of the gas
% so that now in its atomic form, they are diffused between the
% interstitial sites of lattice formation of palladium hydrides In this
% process, the lattice undergoes to isotropic expansion, lattice
% constant is increased until a 3\% in each space direction, increasing
% the total volume until in a 10\%
beginfig(1);
numeric sph, sphc, ind, indos, suma, intres, i, j, coos, size, num;
color xdelta, ydelta, zdelta, W[], fro, pro;
path p[], pp[], ppp, pppp;
f := 1.5*(10,5,3);
% ParallelProj := true;
Spread := 50;
sph := 0.15;
sphc := 0.07;
size = 6;
num=200;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
xdelta := (sph,0,0);
ydelta := (0,sph,0);
zdelta := (0,0,sph);
suma:=0;
for ind=1 upto 3:
for indos=1 upto 10:
for intres=1 upto 10:
suma:=suma+1;
V[suma] :=0.5*(indos,intres,ind);
W[suma] :=0.5*((indos-1,intres-1,ind-1)+0.5*(1,1,1));
endfor;
endfor;
endfor;
for ind=1 upto suma:
p[ind] := rigorousfearpath( V[ind], sph );
pp[ind] := rigorousfearpath( W[ind], sphc );
endfor;
for ind=1 upto suma:
fill pp[ind] withcolor red;
draw pp[ind] withcolor black;
fill p[ind] withcolor .75white;
draw p[ind] withcolor blue;
endfor;
%%%%%%%%%%%%%%%%%%%%%%
for i=1 upto num:
for j=1 upto 2:
coos[j] := 1.7*uniformdeviate(size*.5);
endfor;
fro := (coos1,coos2, 4*(.5)+(.3)*normaldeviate);
pro := fro+(0,sph,0);
ppp := rigorousfearpath( fro, sphc);
pppp := rigorousfearpath( pro, sphc);
fill ppp withcolor red;
draw ppp withcolor (.5,.4,0);
fill pppp withcolor red;
draw pppp withcolor black;
endfor;
%%%%%%%%%%%%%%%%%%%%%%%%
endfig;
end;
|