From 7015cfecb337eff88882902d36d67debe958430b Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Thu, 9 Jan 2014 18:41:09 +0100
Subject: [PATCH] [Algorit] Simplify normal stress computation
---
src/one-body-sample.cc | 21 +++------------------
1 file changed, 3 insertions(+), 18 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 38659cfd..1d2f8b27 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -188,24 +188,9 @@ int main(int argc, char *argv[]) {
};
Vector ell(fineVertexCount);
computeExternalForces(0.0, ell);
-
- double normalStress;
- {
- double volume = 1.0;
- for (size_t i = 0; i < dims; ++i)
- volume *= (myGeometry.C[i] - myGeometry.A[i]);
-
- double area = 1.0;
- for (size_t i = 0; i < dims; ++i)
- if (i != 1)
- area *= (myGeometry.C[i] - myGeometry.A[i]);
-
- // volume * gravity * density / area = normal stress
- // V * g * rho / A = sigma_n
- // m^d * N/kg * kg/m^d / m^(d-1) = N/m^(d-1)
- normalStress = -volume * parset.get<double>("gravity") *
- parset.get<double>("body.density") / area;
- }
+ double normalStress = (myGeometry.A[1] - myGeometry.C[1]) *
+ parset.get<double>("gravity") *
+ parset.get<double>("body.density");
FrictionData const frictionData(parset.sub("boundary.friction"),
normalStress);
--
GitLab