Triaxial test using LIGGGHTS-WITH-BONDS

Hello,
I am trying to do a triaxial test using LIGGGHTS-WITH-BONDS. I have never really used this software before and I am struggling in understanding what the parameters mean. The particles are in a box were one wall is applying pressure through the mesh/command/stress/servo code. I am trying to calculate the stress the particles are causing on the other walls and my input parameters are as follows:
Minimum ball size -2.5e-4 m
Maximum ball size -4.5e-4 m
Density -2650 kg/m3
Ball and bond elastic modulus - 40 e9 Pa
Ratio of normal to shear stiffness - 2.5
Parallel bond normal strength - 44e6 Pa
Parallel bond shear strength - 44e6 Pa
Bond strength standard deviation - 8e6 Pa
Ball-ball and ball-wall friction - 0.5

Thats the code and it is attached too:
#bond

atom_style hybrid granular bond/gran n_bondtypes 1 bonds_per_atom 6
#atom_style granular
atom_modify map array
boundary f f f
newton off

#processors 2 1 1

communicate single vel yes

units si

region reg block -0.0025 0.0025 -0.0025 0.0025 0.0000 0.008 units box
create_box 1 reg #bonds 1 6

neighbor 0.002 bin
neigh_modify delay 0

#New pair style
pair_style gran model hertz tangential history
pair_coeff * *

#bond
bond_style gran

variable stressbreak equal 1

#fix boxwalls_x1 all wall/gran model hertz tangential history primitive type 1 xplane 0
#fix boxwalls_x2 all wall/gran model hertz tangential history primitive type 1 xplane 2
#fix boxwalls_y1 all wall/gran model hertz tangential history primitive type 1 yplane 0
#fix boxwalls_y2 all wall/gran model hertz tangential history primitive type 1 yplane 2
#fix boxwalls_z1 all wall/gran model hertz tangential history primitive type 1 zplane 0
#fix boxwalls_z2 all wall/gran model hertz tangential history primitive type 1 zplane 2

#for stressbreak
bond_coeff 1 0.5 44e6 44e6 ${stressbreak} 2.5e9 1e9

# ==============================================================
# create particles IN A DEFINED REGION to collide against a wall
# ==============================================================
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2650 radius constant 0.00025
fix pdd1 all particledistribution/discrete 1. 1 pts1 1
fix ins all insert/pack seed 5330 distributiontemplate pdd1 &
maxattempt 100000 insert_every once overlapcheck yes all_in yes vel constant 0. 0. 0. &
region reg volumefraction_region 0.82

###TODO blow the particles up (see examples/LIGGGHTS/Tutorial_Public/packing/###

#Material properties required for new pair styles

fix m1 all property/global youngsModulus peratomtype 40e9
fix m2 all property/global poissonsRatio peratomtype 0.3
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.2
fix m4 all property/global coefficientFriction peratomtypepair 1 0.5

mass 1 1.0 #dummy

fix bondcr all bond/create/gran 1 1 1 0.0001 1 6

timestep 0.00001

fix cad all mesh/surface/stress file lwall.stl type 1 stress on
fix cad2 all mesh/surface/stress/servo file uwall.stl type 1 com 0. 0. 0.008 ctrlPV force axis 0. 0. 7000 target_val -7000 vel_max 0.008 kp 1e-2
fix cad3 all mesh/surface/stress file side1.stl type 1 stress on
fix cad4 all mesh/surface/stress file side2.stl type 1 stress on
fix cad5 all mesh/surface/stress file side3.stl type 1 stress on
fix cad6 all mesh/surface/stress file side4.stl type 1 stress on
fix geometry all wall/gran model hertz tangential history mesh n_meshes 6 meshes cad cad2 cad3 cad4 cad5 cad6

#apply nve integration to all particles that are inserted as single particles
fix integr all nve/sphere

#output settings, include total thermal energyg
compute 1 all erotate/sphere
variable etotal equal ke+c_1
thermo_style custom step atoms numbond
thermo 1000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes

#insert the first particles so that dump is not empty
shell mkdir post
dump dmp all custom 100 post/dump*.liggghts id type x y z vx vy vz fx fy fz omegax omegay omegaz radius
dump dumpstress all mesh/gran/VTK 200 post/dump*.vtk stress

compute b1 all property/local batom1x batom1y batom1z batom2x batom2y batom2z batom1 batom2 btype
dump bnd all local 100 post/bonds*.bond c_b1[1] c_b1[2] c_b1[3] c_b1[4] c_b1[5] c_b1[6] c_b1[7] c_b1[8] c_b1[9]

#insert particles
run 1
fix_modify bondcr every 0

# ================================================
# after the bonds are formed,
# exclude pairwise interaction ?
# ================================================
#neigh_modify exclude group bonded bonded

#run

run 300000

#run 3500 upto
#dump_modify dmp every 5
#run 500

I am not really seeing any stress response on the walls and Im not sure if my bond coefficients are correct. I would really appreciate some help!
Thank you!