2_shear.in
View Options

			read_restart confined.restart

			variable nb_points equal 1000
			variable eps_max equal 0.1
			variable shearrate equal 1e-1
			variable dt equal 0.01

			# Don't touch
			variable eps equal xy/lx
			variable rate equal v_shearratelxv_dt
			variable nb_step equal $(round(v_eps_max/(v_rate*v_dt)))
			variable print_freq equal $(round(v_nb_step/v_nb_points))

			# ===== GEOMETRY ====

			change_box all triclinic

			# ==== DUMP ====

			thermo_style custom step atoms lx xy pxy v_eps v_nb_step
			thermo_modify norm no flush yes
			thermo 5000

			# ==== FIXES ====

			timestep ${dt}

			fix 2 all deform 1 xy erate ${rate}
			fix FIX5 all enforce2d

			fix output_file all print ${print_freq} "${eps} $(pxy)" file StrainStress.txt screen no

			run ${nb_step}

			write_data box_sheared.data

DD_1D-SRA/DATASET/P=120000Pa/2_shear.inf855cb73c1c5main