LECSM/output_tools.cpp at master · OptimalDesignLab/LECSM · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
/**
 * \file output_tools.cpp
 * \brief a set of tools to output FEA analysis results
 * \author  Alp Dener <alp.dener@gmail.com>
 * \version 1.0
 */

#include <stdio.h>
#include "./output_tools.hpp"
using namespace std;

// ======================================================================

void output_disp(int nnp, vector<double>& G,
                 vector< vector< vector<int> > >& gm,
                 vector<double>& disp, vector< vector<double> >& nodeDisp)
{
  printf("Directions:\n");
  printf("  0 - x-axis\n");
  printf("  1 - y-axis\n");
  printf("  2 - rotation about z-axis\n");
  vector<double> nd(2);
  for (int A = 0; A < nnp; A++)
  {
    for (int i = 0; i < 3; i++)
    {
      int t = gm[i][A][0];
      double P = gm[i][A][1];
      if (t == 1)             // dof
        {nd[i] = disp[P];}
      else
        {
          if (t == 2)          // dog
            {nd[i] = G[P];}
          else
            {nd[i] = 0;}
        }
      printf("    Node %d displaced %f in direction %d\n", A, nd[i], i);
      nodeDisp[A][i] = nd[i];
    }
  }
}