simfire/fiber_linear.hh at master · jstebel/simfire · 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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
#ifndef FIBER_LINEAR_HH
#define FIBER_LINEAR_HH


#include "fiber.hh"


namespace Composite_elasticity_problem
{
  using namespace dealii;


  /**
   * Linear response in 3 directions.
   */
  class FiberLinear : public FiberSubproblem
  {
  public:
	  FiberLinear(const std::string &mesh_file);
	  ~FiberLinear() override;

	  void attach_matrix_handler(const DoFHandler<3> &dh_) override;

	  void assemble_constraint_mat(SparseMatrix<double> &cm0t, SparseMatrix<double> &cm1t, SparseMatrix<double> &cm0, SparseMatrix<double> &cm1) override;

	  void assemble_fiber_matrix(SparseMatrix<double> &fiber_matrix, Vector<double> &fiber_rhs, double E_fiber, double nu_fiber, double fiber_volume) override;

	  void output_results(const std::string &base_path, const Vector<double> &solution_3d, const Vector<double> &solution_1d) const override;

	  DoFHandler<1,3> &dof_handler()  override { return dh; }

	  unsigned int get_constraint_matrix_n_rows() override;

	  void set_sparsity_pattern(SparsityPattern &sp) override;

	  void set_constraint_matrix_sparsity_pattern(SparsityPattern&, SparsityPattern&, SparsityPattern&, SparsityPattern&) override;

  private:

	  Triangulation<1,3> tri;
	  FESystem<1,3> fe;
	  DoFHandler<1,3> dh;
	  const DoFHandler<3> *dh_3d;
	  Quadrature<1> q_constraint;

	  std::vector<std::pair<DoFHandler<3>::active_cell_iterator,Point<3> > > node_to_cell;

  };


}

#endif