Homec4science
Diffusion lammps (lammm-master-deprecated) Browse

lammps/bench-accel/bench_peri/log.peri.4
60e58c44cf1flammm-master-deprecated

Last Change
Raw File

log.peri.4
View Options

LAMMPS (23 Sep 2010)
# Crack growth in notched 3D Peridynamic block
# Mesh spacing
variable h equal 5.00e-4
# Peridynamic horizon
variable delta equal 3.0*${h}
variable delta equal 3.0*0.0005
# Height of plate (meters)
variable height equal 31.5*${h}
variable height equal 31.5*0.0005
# Width of plate (meters)
variable width equal 39.5*${h}
variable width equal 39.5*0.0005
# Thickness of plate (meters)
variable depth equal 24.5*${h}
variable depth equal 24.5*0.0005
# Height of notch
variable crackheight equal 10*${h}
variable crackheight equal 10*0.0005
# Density of plate
variable mydensity equal 2440.0
# Elastic modulus of material
variable myE equal 72.0e9
# Strain energy release rate at branching
variable myG equal 135.0
# constant, but define it as a variable here
variable pi equal 3.14159265358979323846
units si
boundary s s s
atom_style peri
atom_modify map array
variable myskin equal 2.0*${h}
variable myskin equal 2.0*0.0005
neighbor ${myskin} bin
neighbor 0.001 bin
lattice sc $h
lattice sc 0.0005
Lattice spacing in x,y,z = 0.0005 0.0005 0.0005
variable myxmin equal 0.0
variable myxmax equal ${width}
variable myxmax equal 0.01975
variable myymin equal 0.0
variable myymax equal ${height}
variable myymax equal 0.01575
variable myzmin equal 0.0
variable myzmax equal ${depth}
variable myzmax equal 0.01225
region plate block ${myxmin} ${myxmax} ${myymin} ${myymax} ${myzmin} ${myzmax} units box
region plate block 0 ${myxmax} ${myymin} ${myymax} ${myzmin} ${myzmax} units box
region plate block 0 0.01975 ${myymin} ${myymax} ${myzmin} ${myzmax} units box
region plate block 0 0.01975 0 ${myymax} ${myzmin} ${myzmax} units box
region plate block 0 0.01975 0 0.01575 ${myzmin} ${myzmax} units box
region plate block 0 0.01975 0 0.01575 0 ${myzmax} units box
region plate block 0 0.01975 0 0.01575 0 0.01225 units box
create_box 3 plate
Created orthogonal box = (0 0 0) to (0.01975 0.01575 0.01225)
2 by 2 by 1 processor grid
create_atoms 1 region plate
Created 32000 atoms
pair_style peri/pmb
variable myk equal (2.0/3.0)*${myE}
variable myk equal (2.0/3.0)*7.2e+10
variable myc equal ((18.0*${myk})/(${pi}*(${delta}^4)))
variable myc equal ((18.0*4.8e+10)/(${pi}*(${delta}^4)))
variable myc equal ((18.0*4.8e+10)/(3.141592654*(${delta}^4)))
variable myc equal ((18.0*4.8e+10)/(3.141592654*(0.0015^4)))
variable mydelta equal (${delta}+(${delta}/100.0))
variable mydelta equal (0.0015+(${delta}/100.0))
variable mydelta equal (0.0015+(0.0015/100.0))
variable mys0 equal sqrt((5.0*${myG})/(9.0*${myk}*${delta}))
variable mys0 equal sqrt((5.0*135)/(9.0*${myk}*${delta}))
variable mys0 equal sqrt((5.0*135)/(9.0*4.8e+10*${delta}))
variable mys0 equal sqrt((5.0*135)/(9.0*4.8e+10*0.0015))
variable tmpvar1 equal ${myymax}-${crackheight}
variable tmpvar1 equal 0.01575-${crackheight}
variable tmpvar1 equal 0.01575-0.005
variable tmpvar2 equal 0.5*${width}
variable tmpvar2 equal 0.5*0.01975
region topleft block 0.0 ${tmpvar2} ${tmpvar1} ${myymax} ${myzmin} ${myzmax} units box
region topleft block 0.0 0.009875 ${tmpvar1} ${myymax} ${myzmin} ${myzmax} units box
region topleft block 0.0 0.009875 0.01075 ${myymax} ${myzmin} ${myzmax} units box
region topleft block 0.0 0.009875 0.01075 0.01575 ${myzmin} ${myzmax} units box
region topleft block 0.0 0.009875 0.01075 0.01575 0 ${myzmax} units box
region topleft block 0.0 0.009875 0.01075 0.01575 0 0.01225 units box
region topright block ${tmpvar2} ${myxmax} ${tmpvar1} ${myymax} ${myzmin} ${myzmax} units box
region topright block 0.009875 ${myxmax} ${tmpvar1} ${myymax} ${myzmin} ${myzmax} units box
region topright block 0.009875 0.01975 ${tmpvar1} ${myymax} ${myzmin} ${myzmax} units box
region topright block 0.009875 0.01975 0.01075 ${myymax} ${myzmin} ${myzmax} units box
region topright block 0.009875 0.01975 0.01075 0.01575 ${myzmin} ${myzmax} units box
region topright block 0.009875 0.01975 0.01075 0.01575 0 ${myzmax} units box
region topright block 0.009875 0.01975 0.01075 0.01575 0 0.01225 units box
set region topleft type 2
5000 settings made for type
set region topright type 3
5000 settings made for type
pair_coeff 1 1 ${myc} ${mydelta} ${mys0} 0.0
pair_coeff 1 1 5.432488723e+22 ${mydelta} ${mys0} 0.0
pair_coeff 1 1 5.432488723e+22 0.001515 ${mys0} 0.0
pair_coeff 1 1 5.432488723e+22 0.001515 0.001020620726 0.0
pair_coeff 2 2 ${myc} ${mydelta} ${mys0} 0.0
pair_coeff 2 2 5.432488723e+22 ${mydelta} ${mys0} 0.0
pair_coeff 2 2 5.432488723e+22 0.001515 ${mys0} 0.0
pair_coeff 2 2 5.432488723e+22 0.001515 0.001020620726 0.0
pair_coeff 3 3 ${myc} ${mydelta} ${mys0} 0.0
pair_coeff 3 3 5.432488723e+22 ${mydelta} ${mys0} 0.0
pair_coeff 3 3 5.432488723e+22 0.001515 ${mys0} 0.0
pair_coeff 3 3 5.432488723e+22 0.001515 0.001020620726 0.0
pair_coeff 2 3 ${myc} 0.0 ${mys0} 0.0
pair_coeff 2 3 5.432488723e+22 0.0 ${mys0} 0.0
pair_coeff 2 3 5.432488723e+22 0.0 0.001020620726 0.0
pair_coeff 1 2 ${myc} ${mydelta} ${mys0} 0.0
pair_coeff 1 2 5.432488723e+22 ${mydelta} ${mys0} 0.0
pair_coeff 1 2 5.432488723e+22 0.001515 ${mys0} 0.0
pair_coeff 1 2 5.432488723e+22 0.001515 0.001020620726 0.0
pair_coeff 1 3 ${myc} ${mydelta} ${mys0} 0.0
pair_coeff 1 3 5.432488723e+22 ${mydelta} ${mys0} 0.0
pair_coeff 1 3 5.432488723e+22 0.001515 ${mys0} 0.0
pair_coeff 1 3 5.432488723e+22 0.001515 0.001020620726 0.0
set group all density ${mydensity}
set group all density 2440
32000 settings made for density
variable myvolume equal ($h)^3
variable myvolume equal (0.0005)^3
set group all volume ${myvolume}
set group all volume 1.25e-10
32000 settings made for volume
velocity all set 0.0 0.0 0.0 sum no units box
fix F1 all nve
compute C1 all damage/atom
velocity all ramp vx -10.0 10.0 x ${myxmin} ${myxmax} units box
velocity all ramp vx -10.0 10.0 x 0 ${myxmax} units box
velocity all ramp vx -10.0 10.0 x 0 0.01975 units box
variable mystep equal 0.8*sqrt((2.0*${mydensity})/(512*(${myc}/$h)*${myvolume}))
variable mystep equal 0.8*sqrt((2.0*2440)/(512*(${myc}/$h)*${myvolume}))
variable mystep equal 0.8*sqrt((2.0*2440)/(512*(5.432488723e+22/$h)*${myvolume}))
variable mystep equal 0.8*sqrt((2.0*2440)/(512*(5.432488723e+22/0.0005)*${myvolume}))
variable mystep equal 0.8*sqrt((2.0*2440)/(512*(5.432488723e+22/0.0005)*1.25e-10))
timestep ${mystep}
timestep 2.119314924e-08
thermo 20
run 100
Peridynamic bonds:
total # of bonds = 3457032
bonds/atom = 108.032
Memory usage per processor = 53.3624 Mbytes
Step Temp E_pair E_mol TotEng Press Volume
0 2.0134233e+27 0 0 1.3342785e+09 2.3573796e+14 3.773338e-06
20 1.7695805e+27 1.6163291e+08 0 1.3343188e+09 -1.1615176e+17 3.773338e-06
40 1.3041477e+27 4.6848143e+08 0 1.332729e+09 -1.6325875e+17 3.773338e-06
60 9.8975313e+26 5.7284448e+08 0 1.2287455e+09 -1.3858167e+17 3.773338e-06
80 9.3888573e+26 4.0928092e+08 0 1.0314725e+09 -4.4442267e+16 3.773338e-06
100 8.3930314e+26 3.8522361e+08 0 9.4142265e+08 5.4468637e+16 3.773338e-06
Loop time of 18.7579 on 4 procs for 100 steps with 32000 atoms
Pair time (%) = 18.0539 (96.2472)
Neigh time (%) = 0 (0)
Comm time (%) = 0.625693 (3.33563)
Outpt time (%) = 0.000966132 (0.00515054)
Other time (%) = 0.0772813 (0.411994)
Nlocal: 8000 ave 8000 max 8000 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Nghost: 5125 ave 5125 max 5125 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Neighs: 1.6861e+06 ave 1.77502e+06 max 1.60625e+06 min
Histogram: 2 0 0 0 0 0 0 0 1 1
FullNghs: 3.37221e+06 ave 3.41832e+06 max 3.3261e+06 min
Histogram: 2 0 0 0 0 0 0 0 0 2
Total # of neighbors = 13488836
Ave neighs/atom = 421.526
Neighbor list builds = 0
Dangerous builds = 0
c4science · Help