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