From 06fa5bbb94879c4c4b6eae7b1e23b16f34d14bb6 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 22 Oct 2009 14:22:19 +0000
Subject: [PATCH] The rod forces need to be flipped. That used to be in
getResultantForce, but here is a better place
[[Imported from SVN: r5107]]
---
dirneucoupling.cc | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index bfeb0597..15cc2092 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -364,6 +364,10 @@ int main (int argc, char *argv[]) try
FieldVector<double,dim> resultantForce, resultantTorque;
resultantForce = rodAssembler.getResultantForce(couplingBoundary, rodX, resultantTorque);
+ // Flip orientation
+ resultantForce *= -1;
+ resultantTorque *= -1;
+
std::cout << "resultant force: " << resultantForce << std::endl;
std::cout << "resultant torque: " << resultantTorque << std::endl;
--
GitLab