in.indent.min
a8ae8866676emaster

in.indent.min
View Options

			# 2d indenter simulation with minimization instead of dynamics

			dimension 2
			boundary p s p

			atom_style atomic
			neighbor 0.3 bin
			neigh_modify delay 5

			# create geometry

			lattice hex 0.9
			region box block 0 20 0 10 -0.25 0.25
			create_box 2 box
			create_atoms 1 box

			mass * 1.0

			# LJ potentials

			pair_style lj/cut 2.5
			pair_coeff * * 1.0 1.0 2.5

			# define groups

			region 1 block INF INF INF 1.25 INF INF
			group lower region 1
			group mobile subtract all lower
			set group lower type 2

			# initial velocities

			fix 2 lower setforce 0.0 0.0 0.0

			# indenter

			fix 5 all enforce2d

			# minimize with indenter

			thermo 10

			#dump 1 all atom 10 dump.indent

			#dump 2 all image 25 image.*.jpg type type &
			# zoom 1.6 adiam 1.5
			#dump_modify 2 pad 4

			#dump 3 all movie 25 movie.mpg type type &
			# zoom 1.6 adiam 1.5
			#dump_modify 3 pad 4

			minimize 1.0e-6 1.0e-6 1000 1000

			variable k equal 5000.0/xlat
			variable k1 equal 1000.0/xlat

			fix 4 all indent $k sphere 10 13.0 0 6.0
			fix_modify 4 energy yes
			minimize 1.0e-6 1.0e-6 1000 1000

			fix 4 all indent ${k1} sphere 10 12.5 0 6.0
			fix_modify 4 energy yes
			minimize 1.0e-6 1.0e-6 1000 1000

			fix 4 all indent ${k1} sphere 10 12.0 0 6.0
			fix_modify 4 energy yes
			minimize 1.0e-6 1.0e-6 1000 1000

			fix 4 all indent ${k1} sphere 10 11.4 0 6.0
			fix_modify 4 energy yes
			minimize 1.0e-6 1.0e-6 1000 1000

			fix 4 all indent ${k1} sphere 10 11.2 0 6.0
			fix_modify 4 energy yes
			minimize 1.0e-6 1.0e-6 1000 1000

			fix 4 all indent ${k1} sphere 10 11.0 0 6.0
			fix_modify 4 energy yes
			minimize 1.0e-6 1.0e-6 1000 1000