Dear Dr Axel
Thanks alot for your answe
I applied the changes that you mentioned, the test piece I have fixed dimensions and the failure force that was created in the test(ASTM D638), I wrote it in the code, and with the speed, there was no change in the simulation and the fracture wasn’t created. so I exert force , but the piece did not break out of the crack (like the picture)
Is the selected potential type incorrect?
Are the region not correct?
Is the fix section not correct?
thanks for answer
units si
dimension 3
boundary s s p
atom_style peri
atom_modify map array
neighbor 0.0012 bin
create geometry
lattice sc 0.0006
region plate block -0.0125 0.0125 -0.0285 0.0285 -0.004 0.004 units box
create_box 7 plate
create_atoms 1 region plate
region center_circle sphere 0 0 0 0.005 units box
group center_circle region center_circle
delete_atoms group center_circle
region crack block -0.0125 -0.0095 0 0.0002 -0.004 0.004 units box
group crack region crack
delete_atoms group crack
pair_style peri/lps
pair_coeff * * 2.89e9 1.05e9 0.002005 0.03867 0.22
pair_coeff 3 4 2.89e9 1.05e9 0 0.03867 0.22
set group all density 1180
set group all volume 36e-8
velocity all set 0.0 0.0 0.0 sum no units box
region 2 block INF INF INF -0.007 INF INF units box
group lower region 2
region 5 block INF INF 0.007 INF INF INF units box
group upper region 5
region 3 block INF INF INF 0.0002 INF INF units box
group leftlow region 3
region 4 block INF INF 0.0002 INF INF INF units box
group leftup region 4
set group lower type 2
set group leftlow type 3
set group leftup type 4
set group upper type 5
fix 1 all nve
fix 4 upper setforce 0 3.09e8 0
fix 3 lower setforce 0 -3.09e8 0
compute 1 all damage/atom
timestep 3e-8
thermo 20
dump 1 all custom 20 dump.lammpstrj id x y z c_1
dump 2 all xyz 20 dump.xyz
run 5000
