simfire/coupled.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
58
59
60
61
62
#ifndef COUPLED_HH
#define COUPLED_HH
// The material parameters for matrix (PU) are taken from
// http://www.iplex.com.au/iplex.php?page=lib&lib=1&sec=2
// The values for carbon fibres are taken from
// http://www.performance-composites.com/carbonfibre/mechanicalproperties_2.asp


#include <deal.II/lac/block_sparse_matrix.h>
#include <deal.II/lac/block_vector.h>


#include <fstream>
#include <iostream>

#include "parameters.hh"


namespace Composite_elasticity_problem
{
  using namespace dealii;


  template<int dim> class ElasticProblem;
  class FiberSubproblem;


  class CoupledProblem
  {
  public:
	  CoupledProblem(const std::string &input_file);
	  ~CoupledProblem();

	  void run ();


  private:
	  void setup_system();
	  void solve();

	  BlockSparsityPattern bsp;
	  BlockSparseMatrix<double> bm;
	  BlockVector<double> brhs;

	  ElasticProblem<3> *elastic;
	  FiberSubproblem *fibers;

	  Parameters::AllParameters parameters;

  };


}

#endif