/usr/include/gmsh/elasticitySolver.h is in libgmsh-dev 2.10.1+dfsg1-1ubuntu4.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 | // Gmsh - Copyright (C) 1997-2015 C. Geuzaine, J.-F. Remacle
//
// See the LICENSE.txt file for license information. Please report all
// bugs and problems to the public mailing list <gmsh@geuz.org>.
#ifndef _ELASTICITY_SOLVER_H_
#define _ELASTICITY_SOLVER_H_
#include <map>
#include <string>
#include "SVector3.h"
#include "dofManager.h"
#include "functionSpace.h"
template <class scalar> class simpleFunction;
class GModel;
class PView;
class groupOfElements;
struct LagrangeMultiplierField {
int _tag;
groupOfElements *g;
double _tau;
SVector3 _d;
simpleFunction<double> *_f;
LagrangeMultiplierField() : _tag(0), g(0){}
};
struct elasticField {
int _tag; // tag for the dofManager
groupOfElements *g; // support for this field
double _E, _nu; // specific elastic datas (should be somewhere else)
elasticField () : _tag(0), g(0){}
};
struct BoundaryCondition
{
int _tag; // tag for the dofManager
enum location{UNDEF, ON_VERTEX, ON_EDGE, ON_FACE, ON_VOLUME};
location onWhat; // on vertices or elements
groupOfElements *g; // support for this BC
BoundaryCondition() : _tag(0), onWhat(UNDEF), g(0) {}
};
struct dirichletBC : public BoundaryCondition
{
int _comp; // component
simpleFunction<double> *_f;
dirichletBC ():BoundaryCondition(), _comp(0), _f(0){}
};
struct neumannBC : public BoundaryCondition
{
simpleFunction<SVector3> *_f;
neumannBC () : BoundaryCondition(), _f(NULL){}
};
// an elastic solver ...
class elasticitySolver
{
public:
GModel *pModel;
int _dim, _tag;
dofManager<double> *pAssembler;
FunctionSpace<SVector3> *LagSpace;
FunctionSpace<double> *LagrangeMultiplierSpace;
// young modulus and poisson coefficient per physical
std::vector<elasticField> elasticFields;
std::vector<LagrangeMultiplierField> LagrangeMultiplierFields;
// neumann BC
std::vector<neumannBC> allNeumann;
// dirichlet BC
std::vector<dirichletBC> allDirichlet;
public:
elasticitySolver(int tag) : _tag(tag), pAssembler(0), LagSpace(0), LagrangeMultiplierSpace(0) {}
elasticitySolver(GModel *model, int tag);
void addDirichletBC(int dim, int entityId, int component, double value);
void addDirichletBC(int dim, std::string phys, int component, double value);
void addNeumannBC(int dim, int entityId, const std::vector<double> value);
void addNeumannBC(int dim, std::string phys, const std::vector<double> value);
void addElasticDomain(int tag, double e, double nu);
void addElasticDomain(std::string phys, double e, double nu);
virtual ~elasticitySolver()
{
if (LagSpace) delete LagSpace;
if (LagrangeMultiplierSpace) delete LagrangeMultiplierSpace;
if (pAssembler) delete pAssembler;
}
void assemble (linearSystem<double> *lsys);
void readInputFile(const std::string &meshFileName);
void read(const std::string s) {readInputFile(s.c_str());}
virtual void setMesh(const std::string &meshFileName);
void solve();
void postSolve();
void exportKb();
void computeEffectiveStiffness(std::vector<double> stiff);
void computeEffectiveStrain(std::vector<double> strain);
virtual PView *buildDisplacementView(const std::string postFileName);
virtual PView *buildStrainView(const std::string postFileName);
virtual PView *buildStressesView(const std::string postFileName);
virtual PView *buildLagrangeMultiplierView(const std::string posFileName);
virtual PView *buildElasticEnergyView(const std::string postFileName);
virtual PView *buildVonMisesView(const std::string postFileName);
virtual PView *buildVolumeView(const std::string postFileName);
// std::pair<PView *, PView*> buildErrorEstimateView
// (const std::string &errorFileName, double, int);
// std::pair<PView *, PView*> buildErrorEstimateView
// (const std::string &errorFileName, const elasticityData &ref, double, int);
};
#endif
|