diff --git a/src/factories/stackedblocksfactory.cc b/src/factories/stackedblocksfactory.cc
index 2442c3b01475647cb6320e90aa694774720c238d..4ecbf45e03680e23dd488f6d0325589f238f581f 100644
--- a/src/factories/stackedblocksfactory.cc
+++ b/src/factories/stackedblocksfactory.cc
@@ -73,11 +73,13 @@ void StackedBlocksFactory<HostGridType, VectorType>::setBodies() {
 
         for (size_t i=1; i<this->bodyCount_-1; i++) {
             origins[i] = cuboidGeometries_[i-1]->upperLeft();
-            GlobalCoords lowerWeakOrigin_i = lowerWeakOrigin + origins[i];
-            GlobalCoords upperWeakOrigin_i = upperWeakOrigin + origins[i];
+            auto height_i = height/3.0;
 
-            cuboidGeometries_[i] = std::make_shared<CuboidGeometry>(origins[i], length, height);
+            GlobalCoords lowerWeakOrigin_i = lowerWeakOrigin + origins[i];
+            GlobalCoords upperWeakOrigin_i = {upperWeakOrigin[0], height_i};
+            upperWeakOrigin_i += origins[i];
 
+            cuboidGeometries_[i] = std::make_shared<CuboidGeometry>(origins[i], length, height_i);
 
             cuboidGeometries_[i]->addWeakeningPatch(subParset, upperWeakOrigin_i, weakBound(upperWeakOrigin_i));
             cuboidGeometries_[i]->addWeakeningPatch(subParset, lowerWeakOrigin_i, weakBound(lowerWeakOrigin_i));
@@ -85,7 +87,6 @@ void StackedBlocksFactory<HostGridType, VectorType>::setBodies() {
 
         const size_t idx = this->bodyCount_-1;
         origins[idx] = cuboidGeometries_[idx-1]->upperLeft();
-        origins[idx][0] += 0.1;
         lowerWeakOrigin += origins[idx];
 
         cuboidGeometries_[idx] = std::make_shared<CuboidGeometry>(origins[idx], length, height);
diff --git a/src/multi-body-problem-2D.cfg b/src/multi-body-problem-2D.cfg
index 676d1a408726cc0c1a71d3c87010a35182156147..a308c80b63b9cb2b48f0cc096ade7d9cdf464e57 100644
--- a/src/multi-body-problem-2D.cfg
+++ b/src/multi-body-problem-2D.cfg
@@ -1,6 +1,6 @@
 # -*- mode:conf -*-
 [boundary.friction]
-smallestDiameter = 0.5  # 2e-3 [m]
+smallestDiameter = 0.01  # 2e-3 [m]
 
 [timeSteps]
 refinementTolerance = 1e-5 # 1e-5
diff --git a/src/multi-body-problem.cfg b/src/multi-body-problem.cfg
index 7188eed1b5851c795b2903d8000846c4a89c6450..e8a3e97bb00a8c6044364082c4800956d6d2f085 100644
--- a/src/multi-body-problem.cfg
+++ b/src/multi-body-problem.cfg
@@ -11,7 +11,7 @@ vtk.write       = true
 
 [problem]
 finalTime       = 100  # [s] #1000
-bodyCount       = 2
+bodyCount       = 4
 
 [body]
 bulkModulus     = 0.5e5 # [Pa]