Enhancing LAMMPS Capabilities

LIGGGHTS Models - ReleasedHertz-Mindlin (HM) Granular Contact ModelF = ( k δn− γ vn ) + ( k δ t − γ vt )n ij{n ij{t ij{t ij{normal normaltangential tangentialoverlap relative vel .overlap relative vel .14444244443normalforce14444244443tangentialforcekγknnt43* *= YRδn,5= −2β Snm6* *= 8GRδ5γ t= β m6*−2 S t≥n,*≥ 0,0.* *Sn = 2YR δn, St= 8G1*Rδn22( 1−ν) ( 1−ν)ln( e ) 112β = , = + ,22*ln ( e ) + π Y Y1Y2G*1R2=1=R( 2+ ν)( 1−ν) 2( 2+ ν)( 1−)1 12ν2Y1+R11= ,* *1 2 m+1=mY1+mY...Young's modulus G...Shear modulusν...Poissonνratioe...coeff.ofofrestitution122*Christian Doppler Laboratory on Particulate Flow Modelling www.particulate-flow.at 11