Dear lammps users
im currently working on buckling of a graphene sheet which the boundaries are simply supported. im trying to find out buckling force of the defected model. i have tried multiple monoatom vacancy tests. the results show that the buckling force are higher for defected models such as the model with 1.8% defects with respect to the pristine one.
is there something that im missing in the code or is it really the behavior of the model at the atom scale that shows these strange results?
thanks in advance