Issue 34

R. Brighenti et alii, Frattura ed Integrità Strutturale, 34 (2015) 80-89; DOI: 10.3221/IGF-ESIS.34.08

 

 

 

,

( ) n

( ) n W d x x (

( ) (  x x x



)    d

( ) (  x x x

 





)  

f

f

f

) W d

f

f

( ) x

( ) x

( ) x

(6)

i

i

i

i

i

i   x x

where (

) is the Dirac delta function that is numerically approximated by adopting a proper decreasing function i  x x ( W is also called kernel or smoothing function), and h is the smoothing length (or of the distance

 x x

W

(

)

i

support dimension). The numerical evaluation of Eq. (6 1

) is expressed as follows:

f

( ) k x







( W h  x x , )

f

m

( ) x

(7)

i

k

k

i

k 

k

where the index k spans all the particles over the neighbors of that placed at i

x because of the compact support of the

, , k k k m 

are the position, mass and density of the generic particle k , respectively, and

kernel function. The parameters ( ) k f x is the function value at

x

k x . By using a SPH approximation for the material’s stress field and inserting it in the dynamic equilibrium equation,



0      , the following expression can be obtained: i i x b

, ij j

  

  



1, 2, 3 i 

( ) k V V W    x p k ij

( ) , ) h V x V b     p i p

 x x

p i   

,

(8)

(

0

j

k

p

p 

,

k

p m

/ k m V   . By using the potential formulation presented above, k k

p V represents the volume of the particle p and

where we get:

[ ]  k p

V V 







 x x



W

, ) h

( ) k x

(

p

k ij

j

k

p

,

   

   

2

   

   

 

 

s

s

( ) a h A E  

'

x x

x x

1 2

(9)

[ ]  k p

k

p

e kp

k

p

e kp

( )

( )

mk

0,

 







V V

q

p

k

i kp

( )

V V 

s

s

k

p

e kp

e kp

( )

( )







 x x

W

, ) h

( ) k x

(

ij

j

k p

,

where q is the i -th component of the versor connecting the particles k and p , and [ ] k p indicates the generic particle k falling within the neighbor of particle p . ( ) i kp

D YNAMIC EQUILIBRIUM EQUATIONS

he dynamic governing equations of the discretized problem are given by: previous sections), the damping force vector, the external force vector and the vector of actions due to the collision of the particle with the elastic boundaries, respectively. The numerical integration of the motion equations can explicitly be done for each particle, once the force vector acting on it is known. Explicit time integration methods, even if conditionally stable, are usually adopted in particles simulations to avoid the use of large matrices such as the stiffness one. T i b       Mx F F F F 0  d e or T  Mx F  (10) where M is the mass matrix, x  is the vector of the particles center acceleration (the rotation degrees of freedom are neglected since particles with small rotational inertia are assumed). Further, , , i d e F F F , b F are the internal force vector (see

85