User login

Navigation

You are here

Smoothed Particle Hydrodynamics

Hi all,

         I am new to Imehcanica. I am a second year PhD student in the department of mechanical engineering, IIT Kanpur. I have just started my PhD thesis. The area of work is high velocity impact simulations. I am writing a fortran code for a meshless method smoothed particle hydrodynamics. To validate the code, I am doing a compression test with the code I have written so far. I created a model of 2 cm cube with 1000 particles in it, 10 particles in x,y,z directions respectively. I am compressing the box at a speed of 50m/s in z-direction. The time step I used is 1e-8 sec and Young's modulus, poisson ratio and density are 2e9Pa, 0.25 and 7800 kg/m^3 respectively. The sound speed, c for these material properties comes out to be 506.4 m/s.

        The sph equations that I am using are as given in Libersky's paper - Journal of computational physics 109, 67-75(1993). First I did not use the artificial viscosity term in the momentum equation. I was able to compress the cube upto 30% strain. But after that the particles next to the velocity specified boundary particles(second from topmost and second from bottom most z-plane particles) penetrated through the boundary particles and after 50% strain most of the particles lost neighbor particles and they did not undergo any stress update. I am using Jaumann rate of constitutive equation. I did not model plasticity. I am not using reflecting planes either to simulate boundaries, as given in the paper.

        When I added the artificial viscosity term in the momentum equation, the results were even worse. The particles started flying in the very beginning of the simulation.

        Please help me in this regard. Please let me know if more information is needed from my side.

              Thanking in advance,

Siva Prasad A V S.

ME PhD Student,

IIT Kanpur.

Subscribe to Comments for "Smoothed Particle Hydrodynamics"

Recent comments

More comments

Syndicate

Subscribe to Syndicate
Error | iMechanica

Error

The website encountered an unexpected error. Please try again later.