You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
New New New<br />
⎡v ⎤ ⎡ ix v ⎤ ⎡ ix v ⎤ ix<br />
⎢ ⎥ New 1⎢ ⎥ New 1⎢<br />
⎥ New<br />
⎢viy ⎥ = ⎢viy ⎥ + ⎢viy ⎥<br />
⎢ ⎥ 2⎢ ⎥ 2⎢<br />
⎥<br />
⎣v ⎦ ⎣v ⎦ ⎣v ⎦<br />
New New New<br />
iz n iz n+ 1/2 iz n−1/2<br />
ΚΕΦΑΛΑΙΟ 4 o<br />
Εξίσωση 4.100<br />
Παροµοίως, για τον υπολογισµό της επιτάχυνσης πρέπει πρώτα να υπολογισθεί, όπως<br />
αναφέρθηκε και πιο πάνω, η δύναµη. Η ολική δύναµη που ενεργεί πάνω σε κάθε άτοµο<br />
υπολογίζεται αθροίζοντας τα διανύσµατα των δυνάµεων από τα γειτονικά άτοµα.<br />
F<br />
New<br />
ij<br />
r −r r −r<br />
= = ⇒<br />
N New New N New New<br />
i j ∂ New j i New<br />
Pr i( ) r f<br />
New New ij<br />
j= 1 rij ∂r<br />
j=<br />
1 rij<br />
j≠i j≠i ∑ ∑<br />
New<br />
Fij = ∑ − − −<br />
⎡i⎤ f<br />
j<br />
⎢⎣k⎥⎦ N<br />
New<br />
New New New New New New ij<br />
⎡xj xi , yj yi , zj z<br />
⎢ ⎥<br />
⎣<br />
⎤ i ⎦⎢ ⎥ New<br />
j= 1<br />
rij<br />
j≠i Εξίσωση 4.101<br />
New<br />
New New<br />
⎡f⎤ ⎡xj − x ⎤<br />
ix i<br />
New New<br />
⎢ ⎥ N ⎢ ⎥ exp( 2a( rij r0)) exp( a( rijr0)) New New New<br />
fiy 2aD<br />
yjy ⎣<br />
⎡− − − + − − ⎦<br />
⎤<br />
⎢ ⎥ = ∑ ⎢ − i ⎥<br />
Εξίσωση 4.102<br />
New New 2 New New 2 New New 2<br />
⎢ j=<br />
1<br />
( xi xj ) ( yi yj ) ( zi zj<br />
)<br />
New ⎥ ⎢<br />
j i New New ⎥ − + − + −<br />
f ≠<br />
⎣ iz ⎦ ⎢⎣zj−zi⎥⎦ Η νέα επιτάχυνση υπολογίζεται επίσης από το δεύτερο νόµο του Newton:<br />
Άρα:<br />
F r − r<br />
New<br />
∂<br />
a i = =<br />
Εξίσωση 4.103<br />
New N New New<br />
ij 1 i j<br />
∑<br />
Pr () New New New<br />
r = rij<br />
mi mi j=<br />
1 rij ∂r<br />
j≠i ⎡i⎤ N<br />
New<br />
1<br />
f<br />
New<br />
New New New New New New ij<br />
a ⎡xj xi , yj yi , zj z ⎤<br />
⎢ ⎥<br />
i = ∑<br />
i New<br />
m ⎣ − − − ⎦⎢<br />
j<br />
⎥<br />
Εξίσωση 4.104<br />
i j=<br />
1<br />
rij<br />
j≠i ⎢⎣k⎥⎦ New<br />
New New<br />
⎡a⎤ ⎡x ix j − x ⎤ i<br />
N<br />
New<br />
⎢ ⎥ New 1 ⎢ ⎥ f<br />
New New ij<br />
⎢aiy ⎥ = ∑ ⎢yj− yi⎥<br />
New<br />
Εξίσωση 4.105<br />
⎢ mi j=<br />
1<br />
r<br />
New ⎥ ⎢ New New ⎥ ij<br />
a j≠i ⎣ iz ⎦ ⎢⎣zj−zi⎥⎦ 193