diff --git a/@BearImp/calcLoad.m b/@BearImp/calcLoad.m
index 1560bf58f43c51d169461308548d726517cb6009..4c3ec037e160ede73c889226dc657d17a05d0afd 100644
--- a/@BearImp/calcLoad.m
+++ b/@BearImp/calcLoad.m
@@ -61,7 +61,7 @@ B.delta_s=deltaShaft(psi,F_r,F_a,B,G,L,B.delta_s,AddOn)';
 if L.allBallsSameMaterial
      
     for IndPosBearing=1:B.numOfCagePositions
-        B.Q(IndPosBearing) = totalReactionForcePerBallForXYZ(B,L,G,B.delta_s(:,IndPosBearing),L.IndBall_conductive,psi(IndPosBearing));
+        B.Q(IndPosBearing) = max(totalReactionForcePerBallForXYZ(B,L,G,B.delta_s(:,IndPosBearing),L.IndBall_conductive,psi(IndPosBearing)),eps); % limit to 2.2204e-16 N
     end
     
 else
@@ -75,7 +75,7 @@ else
         B.extendedTotalReactionForce=true;
         
         for IndPosBearing=1:B.numOfCagePositions
-            B.Q(posBall_conductive,IndPosBearing)=totalReactionForcePerBallForXYZ(B,L,G,B.delta_s(:,IndPosBearing),L.IndBall_conductive(posBall_conductive),B.psi_conductive(posBall_conductive,IndPosBearing));
+            B.Q(posBall_conductive,IndPosBearing) = max(totalReactionForcePerBallForXYZ(B,L,G,B.delta_s(:,IndPosBearing),L.IndBall_conductive(posBall_conductive),B.psi_conductive(posBall_conductive,IndPosBearing)),eps); % limit to 2.2204e-16 N
         end
     end
 end