Please check my input script first.:
units metal
boundary f f s
atom_style atomic
pair_style morse 3.4
read_data toola_15724.txt
pair_coeff * * 0.3429 1.35588 2.626
pair_coeff 1 2*4 0.1 1.7 2.2
neighbor 2.0 bin
neigh_modify delay 5
neigh_modify exclude type 1 1
group tool type 1
group newtonian type 2
group thermostat type 3
group boundary type 4
group mobile subtract all boundary tool
fix 1 tool rigid single
compute new mobile temp
velocity mobile create 293 5812775 temp new
velocity tool set 5 0 0 sum yes units box
fix 2 mobile nve
fix 3 boundary setforce 0.0 0.0 0.0
fix 4 tool setforce 0.0 0.0 0.0
velocity thermostat scale 293.0
#fix 5 all enforce2d
thermo 500
thermo_style custom step temp epair pe ke etotal vol press
thermo_modify temp new
thermo_modify lost warn
compute_modify new extra 4
dump 1 all atom 500 dump.testcul5morse
dump_modify 1 scale no
dump 2 all custom 500 dump.forcecul15morsetag xsu ysu zsu ix iy iz
timestep 0.002
run 10000
Okay.The script is working just fine but the problem is that during the cutting( projected through VMD) only the middle part of the work-piece can be cut. The boundary layer remain static. The cutting tool simply by-passing the boundary layer. I need to remove these layers!
Please click the link below which shows the snapshots of the cutting process.
http://i877.photobucket.com/albums/ab335/Dinesh_Balakrishnan/pic4_zps996a641a.png
http://i877.photobucket.com/albums/ab335/Dinesh_Balakrishnan/pic3_zpsd0b8ff6d.png
http://i877.photobucket.com/albums/ab335/Dinesh_Balakrishnan/pic2_zpsc52718d9.png
http://i877.photobucket.com/albums/ab335/Dinesh_Balakrishnan/pic1_zps7e31c600.png