problem with an error

Submitted by dixonarulraj on Thu, 11/28/2013 - 17:51

hi,
i am working on porosity of the soil sample. i have attached the input script with this .
Please help me with this error ERROR: Fix insert/pack (id ins): unknown keyword (fix_insert_pack.cpp:106).

Thanks in advance,
Dixon
atom_style granular
atom_modify map array
boundary p p p
newton off
echo both
communicate single vel yes
units si
region reg block -0.05 0.05 -0.05 0.05 0. 0.15 units box
create_box 1 reg
neighbor 0.02 bin
neigh_modify delay 0

#Material Propoerties
fix m1 all property/global youngsModulus peratomtype 10.e6
fix m2 all property/global poissonsRatio peratomtype 0.25
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.89
fix m4 all property/global coefficientFriction peratomtypepair 1 0.1

#New pair style
pair_style gran/hertz/history
pair_coeff * *
timestep 0.001

fix xwalls1 all wall/gran/hertz/history primitive type 1 xplane -0.05
fix xwalls2 all wall/gran/hertz/history primitive type 1 xplane +0.05
fix ywalls1 all wall/gran/hertz/history primitive type 1 yplane -0.05
fix ywalls2 all wall/gran/hertz/history primitive type 1 yplane +0.05
fix zwalls1 all wall/gran/hertz/history primitive type 1 zplane -0.05
fix zwalls2 all wall/gran/hertz/history primitive type 1 zplane +0.05

fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2500 radius constant 0.005
fix pts2 all particletemplate/sphere 1 atom_type 1 density constant 2500 radius constant 0.008
fix pdd1 all particledistribution/discrete 1. 2 pts1 0.3 pts2 0.7

variable growmax equal 1500
variable growrate equal 1.02

compute dia all property/atom diameter

variable mas equal mass(all)
variable mins equal ${mas}
variable volcube equal 0.1*0.1*0.1
variable nins equal (${volcube}-(${mins}/2650))/${volcube} #porosity of sampls

#particle insertion

group nve_group region reg

fix ins nve_group insert/pack seed 5330 distributiontemplate pdd1 maxattempt 500 &
insert_every 1 overlapcheck yes all_in yes region reg voluraction_regionmef 0.65 ntry_mc 5000
fix integr nve_group nve/sphere

#insert the first particles
run 1
dump dmp all custom 350 post/dump*.growing id type type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
unfix ins

label grow
variable gr loop ${growsmax}
variable grownd atom ${growrate}*c_dia
fix grow all adopt 0 atom diameter v_grownd
print '_____________________________________________________'
print '___${particles} growing...loop: ${gr} with growrate: ${growrate}'
run 1500 post no
print '_____________________________________________________'
if "${g} == ${growsmax}" then "jump error"
if "${nins} > 0.28" then "next gr" "jump SELF grow" else "jump SELF growdone"

label error

print 'Warning: error in growing particles'

label growdone
unfix grow

aaigner's picture

aaigner | Sun, 12/01/2013 - 22:03

Do what the error-message says... check your fix insert/pack.

In your case it is a simple spelling error:
voluraction_regionmef 0.65
should be
volumefraction_region 0.65

Cheers
Andi