tensile test on copper

LAMMPS LAMMPS Mailing List Mirror
1

hi lammps

I'm new to lammps,and I tried to do tensile test.
this is my code
i want to know this is true or not

# tensile   test
echo both
dimension  3  
boundary   s  p   p
units  metal  
atom_style  atomic

##########create  box#######
region  copperbox   block  -80  80   -40   40   -40  40  units box
create_box  1  copperbox
lattice  fcc  3.61
region  copper  block  -60  60  -20  20  -20  20  units box 
create_atoms  1  region  copper  
mass  1  63.546 
thermo_modify lost ignore  
region		left  block   -60  -50  -20  20   -20   20    units box
group		left region left
region		right  block 50  60   -20   20   -20   20    units box
group		right region right
group   middle   subtract   all   left   right 
timestep   0.002

pair_style   eam  
pair_coeff  * *  cu_u3.eam
velocity  all  create  300  4928459  mom  yes  rot  yes   dist  uniform  
velocity  left  create   0   4928459   mom  yes  rot  yes   dist   uniform   
velocity  right  create   0   4928459   mom  yes  rot  yes   dist   uniform
fix  1  all  nvt  temp  300  300  0.01  
fix  2  left   setforce   0  0  0
fix  3  right  setforce   0   0   0
fix  4   all  nve  
thermo 100
run  1000
2

hi lammps

I'm new to lammps,and I tried to do tensile test.
this is my code
i want to know this is true or not

sorry, but this mailing list is not an input file approval service.
your results will not be the slightest bit more or less correct if any
of us says "yay" or ":nay". you have to do what everybody else does:
find a few relevant publications and reproduce those results. also,
you should consult with your adviser/colleagues for advice on how to
do your work.

if, however, you think and can convincingly demonstrate that LAMMPS
does not work as documented, then please let us know where exactly the
problem is and provide sufficient detail and information to reproduce
it.

axel.