From 48ab7359e2762112997bd34e724efb2b2df93d30 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 20 Aug 2007 09:08:57 +0000
Subject: [PATCH] write grid and deformation into a single file; fix geometric
moments of the rod, to really model a quadratic cross-section
[[Imported from SVN: r1537]]
---
dirneucoupling.cc | 8 +++++---
1 file changed, 5 insertions(+), 3 deletions(-)
diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index 3eb808e9..f263b96f 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -211,7 +211,7 @@ int main (int argc, char *argv[]) try
// Create a solver for the rod problem
// ///////////////////////////////////////////
RodAssembler<RodGridType> rodAssembler(rodGrid);
- rodAssembler.setShapeAndMaterial(1, 1, 1, 2.5e5, 0.3);
+ rodAssembler.setShapeAndMaterial(1, 1/12, 1/12, 2.5e5, 0.3);
RodSolver<RodGridType> rodSolver;
rodSolver.setup(rodGrid,
@@ -536,8 +536,10 @@ int main (int argc, char *argv[]) try
// //////////////////////////////
// Output result
// //////////////////////////////
- LeafAmiraMeshWriter<GridType>::writeGrid(grid, "grid.result");
- LeafAmiraMeshWriter<GridType>::writeBlockVector(grid, x3d, "grid.sol");
+ LeafAmiraMeshWriter<GridType> amiraMeshWriter(grid);
+ amiraMeshWriter.addVertexData(x3d, grid.leafIndexSet());
+ amiraMeshWriter.write("grid.result");
+
writeRod(rodX, "rod3d.result");
// for (int i=0; i<rodX.size(); i++)
--
GitLab