Hi everybody
I’m working on a test problem related to our dissertation. In this problem we want to apply distributed load to limited surface at the middle of beam in order to bend it. I wrote several script for this purpose but unfortunately it is not work properly!
Here is my questions:
-Which atom style we must use in this simulation?
-How can we apply distributed force at the middle of the beam?
-How can we set clamped boundary condition in both ends of the beam?
Here is one of my scripts for simulating the problem:
#atom defination
units metal
dimension 3
boundary m m m
atom_style atomic
variable latconst equal 3.615
lattice fcc ${latconst}
region box block -16 16 -16 16 -20 110
region nanobeam block 0 4 0 4 0 80
create_box 1 box
lattice fcc ${latconst} orient x 1 0 0 orient y 0 1 0 orient z 0 0 1
create_atoms 1 region nanobeam
#force field
pair_style eam
pair_coeff * * Cu_u3.eam
#setting
variable z1 equal 0
variable z2 equal 10.845
region boundary1 block INF INF INF INF {z1} {z2} units box
group boundary1 region boundary1
variable z1 equal 278.355
variable z2 equal 289.2
region boundary2 block INF INF INF INF {z1} {z2} units box
group boundary2 region boundary2
variable z1 equal 108.45
variable z2 equal 180.75
region disp block INF INF INF INF {z1} {z2} units box
group disp region disp
velocity all create 298 511124 rot yes mom yes
variable harmonic equal 0.0139
compute 1 disp displace/atom
thermo 100
#run the simulation
fix relax all nvt temp 298 298 0.01
variable t equal 0.005
timestep ${t}
dump 1 all xyz 10 visual2.xyz
run 1000
unfix relax
fix zeroing_force1 boundary1 setforce 0.0 0.0 0.0
fix zeroing_force2 boundary2 setforce 0.0 0.0 0.0
fix force disp addforce 0 harmonic 0
fix fix1 all nvt temp 298 298 0.01
run 1000000
print “All done”
Thanks,
Feisal Asbaghian