G+Smo  24.08.0
Geometry + Simulation Modules
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
gsMaterialMatrixEval.h
Go to the documentation of this file.
1 
16 #pragma once
17 
18 #include <gsCore/gsFunction.h>
19 #include <gsKLShell/src/gsMaterialMatrixContainer.h>
20 
21 namespace gismo
22 {
23 
24 template <class T, enum MaterialOutput out> class gsMaterialMatrixEvalSingle;
25 
26 template <class T, enum MaterialOutput out>
27 class gsMaterialMatrixEval : public gsFunction<T>
28 {
29 public:
31  gsMaterialMatrixEval( const gsMaterialMatrixContainer<T> & materialMatrices,
32  const gsFunctionSet<T> * deformed,
33  const gsMatrix<T> z)
34  :
35  m_materialMatrices(materialMatrices),
36  m_deformed(deformed),
37  m_z(z)
38  {
39  for (index_t p = 0; p!=deformed->nPieces(); ++p)
40  GISMO_ASSERT(materialMatrices.piece(p)->initialized(),"Material matrix "<<p<<" is incomplete!");
41 
42  this->_makePieces();
43  }
44 
46  gsMaterialMatrixEval( gsMaterialMatrixBase<T> * materialMatrix,
47  const gsFunctionSet<T> * deformed,
48  const gsMatrix<T> z)
49  :
50  m_materialMatrices(deformed->nPieces()),
51  m_deformed(deformed),
52  m_z(z)
53  {
54  GISMO_ASSERT(materialMatrix->initialized(),"Material matrix is incomplete!");
55  for (index_t p = 0; p!=deformed->nPieces(); ++p)
56  m_materialMatrices.set(p,materialMatrix);
57  this->_makePieces();
58  }
59 
61  gsMaterialMatrixEval( const gsMaterialMatrixContainer<T> & materialMatrices,
62  const gsFunctionSet<T> * undeformed,
63  const gsFunctionSet<T> * deformed,
64  const gsMatrix<T> z)
65  :
66  m_materialMatrices(materialMatrices),
67  m_deformed(deformed),
68  m_z(z)
69  {
70  this->_makePieces(undeformed);
71  }
72 
74  gsMaterialMatrixEval( gsMaterialMatrixBase<T> * materialMatrix,
75  const gsFunctionSet<T> * undeformed,
76  const gsFunctionSet<T> * deformed,
77  const gsMatrix<T> z)
78  :
79  m_materialMatrices(deformed->nPieces()),
80  m_deformed(deformed),
81  m_z(z)
82  {
83  for (index_t p = 0; p!=deformed->nPieces(); ++p)
84  m_materialMatrices.set(p,materialMatrix);
85  this->_makePieces(undeformed);
86  }
87 
89  ~gsMaterialMatrixEval()
90  {
91  freeAll(m_pieces);
92  }
93 
95  short_t domainDim() const {return 2;}
96 
104  short_t targetDim() const { return this->piece(0).targetDim(); }
105 
107  const gsFunction<T> & piece(const index_t p) const
108  {
109  return *m_pieces[p];
110  }
111 
113  void eval_into(const gsMatrix<T>& u, gsMatrix<T>& result) const
114  { GISMO_NO_IMPLEMENTATION; }
115 
116 protected:
117  void _makePieces()
118  {
119  m_pieces.resize(m_deformed->nPieces());
120  for (size_t p = 0; p!=m_pieces.size(); ++p)
121  m_pieces.at(p) = new gsMaterialMatrixEvalSingle<T,out>(p,m_materialMatrices.piece(p),m_deformed,m_z);
122  }
123 
124  void _makePieces(const gsFunctionSet<T> * undeformed)
125  {
126  m_pieces.resize(m_deformed->nPieces());
127  for (size_t p = 0; p!=m_pieces.size(); ++p)
128  m_pieces.at(p) = new gsMaterialMatrixEvalSingle<T,out>(p,m_materialMatrices.piece(p),undeformed,m_deformed,m_z);
129  }
130 
131 protected:
132  gsMaterialMatrixContainer<T> m_materialMatrices;
133  const gsFunctionSet<T> * m_deformed;
134  gsMatrix<T> m_z;
135  mutable std::vector<gsMaterialMatrixEvalSingle<T,out> *> m_pieces;
136 };
137 
146 template <class T, enum MaterialOutput out>
147 class gsMaterialMatrixEvalSingle : public gsFunction<T>
148 {
149 public:
150 
152  gsMaterialMatrixEvalSingle( index_t patch,
153  gsMaterialMatrixBase<T> * materialMatrix,
154  const gsFunctionSet<T> * deformed,
155  const gsMatrix<T> z);
156 
158  gsMaterialMatrixEvalSingle( index_t patch,
159  gsMaterialMatrixBase<T> * materialMatrix,
160  const gsFunctionSet<T> * undeformed,
161  const gsFunctionSet<T> * deformed,
162  const gsMatrix<T> z);
163 
164  gsMaterialMatrixEvalSingle(){};
165 
167  short_t domainDim() const {return 2;}
168 
176  short_t targetDim() const { return targetDim_impl<out>(); }
177 
178 
179 private:
181  template<enum MaterialOutput _out>
182  typename std::enable_if<_out==MaterialOutput::Density , short_t>::type targetDim_impl() const { return 1; };
183 
185  template<enum MaterialOutput _out>
186  typename std::enable_if<_out==MaterialOutput::VectorN ||
187  _out==MaterialOutput::CauchyVectorN ||
188  _out==MaterialOutput::VectorM ||
189  _out==MaterialOutput::CauchyVectorM ||
190  _out==MaterialOutput::Generic , short_t>::type targetDim_impl() const { return 3; };
191 
193  template<enum MaterialOutput _out>
194  typename std::enable_if<_out==MaterialOutput::MatrixA ||
195  _out==MaterialOutput::MatrixB ||
196  _out==MaterialOutput::MatrixC ||
197  _out==MaterialOutput::MatrixD , short_t>::type targetDim_impl() const { return 9; };
198 
200  template<enum MaterialOutput _out>
201  typename std::enable_if<_out==MaterialOutput::Stretch , short_t>::type targetDim_impl() const { return 3; };
202 
204  template<enum MaterialOutput _out>
205  typename std::enable_if<_out==MaterialOutput::PStress ||
206  _out==MaterialOutput::PStressN ||
207  _out==MaterialOutput::PCauchyStressN ||
208  _out==MaterialOutput::PStressM ||
209  _out==MaterialOutput::PCauchyStressM ||
210  _out==MaterialOutput::PStrainN ||
211  _out==MaterialOutput::PStrainM , short_t>::type targetDim_impl() const { return 2; };
212 
214  template<enum MaterialOutput _out>
215  typename std::enable_if<_out==MaterialOutput::StretchDir, short_t>::type targetDim_impl() const { return 9; };
216 
218  template<enum MaterialOutput _out>
219  typename std::enable_if<_out==MaterialOutput::PStressDir, short_t>::type targetDim_impl() const { return 9; };
220 
222  template<enum MaterialOutput _out>
223  typename std::enable_if<_out==MaterialOutput::Spec2CovTransform ||
224  _out==MaterialOutput::Spec2ConTransform ||
225  _out==MaterialOutput::Cov2CartTransform ||
226  _out==MaterialOutput::Con2CartTransform, short_t>::type targetDim_impl() const { return 9; };
227 
229  template<enum MaterialOutput _out>
230  typename std::enable_if<_out==MaterialOutput::TensionField, short_t>::type targetDim_impl() const { return 1; };
231 
233  template<enum MaterialOutput _out>
234  typename std::enable_if<_out==MaterialOutput::StretchTransform ||
235  _out==MaterialOutput::PStressTransform, short_t>::type targetDim_impl() const { return 9; };
236 
238  template<enum MaterialOutput _out>
239  typename std::enable_if<_out==MaterialOutput::Theta ||
240  _out==MaterialOutput::Gamma , short_t>::type targetDim_impl() const { return 1; };
241 
243  template<enum MaterialOutput _out>
244  typename std::enable_if<_out==MaterialOutput::Strain ||
245  _out==MaterialOutput::StrainN ||
246  _out==MaterialOutput::StrainM , short_t>::type targetDim_impl() const { return 3; };
247 
249  template<enum MaterialOutput _out>
250  typename std::enable_if<_out==MaterialOutput::Stress ||
251  _out==MaterialOutput::StressN ||
252  _out==MaterialOutput::StressM ||
253  _out==MaterialOutput::CauchyStressN ||
254  _out==MaterialOutput::CauchyStressM , short_t>::type targetDim_impl() const { return 3; };
255 
257  template<enum MaterialOutput _out>
258  typename std::enable_if<_out==MaterialOutput::Thickness, short_t>::type targetDim_impl() const { return 1; };
259 
261  template<enum MaterialOutput _out>
262  typename std::enable_if<_out==MaterialOutput::Parameters, short_t>::type targetDim_impl() const { return m_materialMat->numParameters(); };
263 
265  template<enum MaterialOutput _out>
266  typename std::enable_if<_out==MaterialOutput::Deformation, short_t>::type targetDim_impl() const { return 9; };
267 
268 protected:
270  void setPatch(index_t p) {m_pIndex = p; }
271 
272 public:
273 
275  void eval_into(const gsMatrix<T>& u, gsMatrix<T>& result) const;
276 
277 private:
279  template<enum MaterialOutput _out>
280  typename std::enable_if<_out==MaterialOutput::Density , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
281 
283  template<enum MaterialOutput _out>
284  typename std::enable_if<_out==MaterialOutput::VectorN ||
285  _out==MaterialOutput::VectorM ||
286  _out==MaterialOutput::Generic , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
287 
289  template<enum MaterialOutput _out>
290  typename std::enable_if<_out==MaterialOutput::CauchyVectorN ||
291  _out==MaterialOutput::CauchyVectorM , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
292 
293 
295  template<enum MaterialOutput _out>
296  typename std::enable_if<_out==MaterialOutput::MatrixA ||
297  _out==MaterialOutput::MatrixB ||
298  _out==MaterialOutput::MatrixC ||
299  _out==MaterialOutput::MatrixD , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
300 
302  template<enum MaterialOutput _out>
303  typename std::enable_if<_out==MaterialOutput::PStress ||
304  _out==MaterialOutput::PStressN ||
305  _out==MaterialOutput::PStressM , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
306 
308  template<enum MaterialOutput _out>
309  typename std::enable_if<_out==MaterialOutput::PCauchyStressN ||
310  _out==MaterialOutput::PCauchyStressM , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
311 
313  template<enum MaterialOutput _out>
314  typename std::enable_if<_out==MaterialOutput::PStrainN ||
315  _out==MaterialOutput::PStrainM , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
316 
318  template<enum MaterialOutput _out>
319  typename std::enable_if<_out==MaterialOutput::Stretch , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
320 
322  template<enum MaterialOutput _out>
323  typename std::enable_if<_out==MaterialOutput::StretchDir, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
324 
326  template<enum MaterialOutput _out>
327  typename std::enable_if<_out==MaterialOutput::PStressDir, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
328 
330  template<enum MaterialOutput _out>
331  typename std::enable_if<_out==MaterialOutput::StretchTransform, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
332 
334  template<enum MaterialOutput _out>
335  typename std::enable_if<_out==MaterialOutput::PStressTransform, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
336 
338  template<enum MaterialOutput _out>
339  typename std::enable_if<_out==MaterialOutput::Spec2CovTransform, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
340 
342  template<enum MaterialOutput _out>
343  typename std::enable_if<_out==MaterialOutput::Spec2ConTransform, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
344 
346  template<enum MaterialOutput _out>
347  typename std::enable_if<_out==MaterialOutput::Cov2CartTransform, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
348 
350  template<enum MaterialOutput _out>
351  typename std::enable_if<_out==MaterialOutput::Con2CartTransform, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
352 
354  template<enum MaterialOutput _out>
355  typename std::enable_if<_out==MaterialOutput::TensionField, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
356 
358  template<enum MaterialOutput _out>
359  typename std::enable_if<_out==MaterialOutput::Theta, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
360 
362  template<enum MaterialOutput _out>
363  typename std::enable_if<_out==MaterialOutput::Gamma, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
364 
366  template<enum MaterialOutput _out>
367  typename std::enable_if<_out==MaterialOutput::Strain ||
368  _out==MaterialOutput::StrainN ||
369  _out==MaterialOutput::StrainM , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
370 
372  template<enum MaterialOutput _out>
373  typename std::enable_if<_out==MaterialOutput::Stress ||
374  _out==MaterialOutput::StressN ||
375  _out==MaterialOutput::StressM , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
376 
378  template<enum MaterialOutput _out>
379  typename std::enable_if<_out==MaterialOutput::CauchyStress ||
380  _out==MaterialOutput::CauchyStressN ||
381  _out==MaterialOutput::CauchyStressM , void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
382 
383 
385  template<enum MaterialOutput _out>
386  typename std::enable_if<_out==MaterialOutput::Thickness, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
387 
389  template<enum MaterialOutput _out>
390  typename std::enable_if<_out==MaterialOutput::Parameters, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
391 
393  template<enum MaterialOutput _out>
394  typename std::enable_if<_out==MaterialOutput::Deformation, void>::type eval_into_impl(const gsMatrix<T>& u, gsMatrix<T>& result) const;
395 
396 protected:
397  index_t m_pIndex;
398  gsMaterialMatrixBase<T> * m_materialMat;
399  gsMatrix<T> m_z;
400 };
401 
402 } // namespace
403 
404 
405 #ifndef GISMO_BUILD_LIB
406 #include GISMO_HPP_HEADER(gsMaterialMatrixEval.hpp)
407 #endif
408 
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::MatrixA||_out==MaterialOutput::MatrixB||_out==MaterialOutput::MatrixC||_out==MaterialOutput::MatrixD, short_t >::type targetDim_impl() const
Implementation of targetDim for material tensors.
Definition: gsMaterialMatrixEval.h:197
GISMO_NO_IMPLEMENTATION
#define GISMO_NO_IMPLEMENTATION
Definition: gsDebug.h:129
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::Spec2CovTransform||_out==MaterialOutput::Spec2ConTransform||_out==MaterialOutput::Cov2CartTransform||_out==MaterialOutput::Con2CartTransform, short_t >::type targetDim_impl() const
Implementation of targetDim for transformations.
Definition: gsMaterialMatrixEval.h:226
short_t
#define short_t
Definition: gsConfig.h:35
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::Theta||_out==MaterialOutput::Gamma, short_t >::type targetDim_impl() const
Implementation of targetDim for TFT variables.
Definition: gsMaterialMatrixEval.h:240
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::PStressDir, short_t >::type targetDim_impl() const
Implementation of targetDim for principal stress directions.
Definition: gsMaterialMatrixEval.h:219
gismo::gsMaterialMatrixEvalSingle
This class serves as the evaluator of material matrices, based on gsMaterialMatrixBase.
Definition: gsMaterialMatrixEval.h:24
gismo::gsMaterialMatrixContainer
This class serves as the evaluator of material matrices, based on gsMaterialMatrixBase.
Definition: gsMaterialMatrixContainer.h:33
index_t
#define index_t
Definition: gsConfig.h:32
gismo::gsFunction
A function from a n-dimensional domain to an m-dimensional image.
Definition: gsFunction.h:59
gismo::gsMaterialMatrixEvalSingle::eval_into
void eval_into(const gsMatrix< T > &u, gsMatrix< T > &result) const
Implementation of eval_into, see gsFunction.
Definition: gsMaterialMatrixEval.hpp:72
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::Stress||_out==MaterialOutput::StressN||_out==MaterialOutput::StressM||_out==MaterialOutput::CauchyStressN||_out==MaterialOutput::CauchyStressM, short_t >::type targetDim_impl() const
Specialisation of targetDim for strain.
Definition: gsMaterialMatrixEval.h:254
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::Parameters, short_t >::type targetDim_impl() const
Implementation of targetDim for the parameters.
Definition: gsMaterialMatrixEval.h:262
GISMO_ASSERT
#define GISMO_ASSERT(cond, message)
Definition: gsDebug.h:89
gismo::gsMaterialMatrixBase
This class defines the base class for material matrices.
Definition: gsMaterialMatrixBase.h:32
gismo::gsMaterialMatrixEvalSingle::eval_into_impl
std::enable_if< _out==MaterialOutput::Density, void >::type eval_into_impl(const gsMatrix< T > &u, gsMatrix< T > &result) const
Specialisation of eval_into for densities.
Definition: gsMaterialMatrixEval.hpp:82
gismo::gsMaterialMatrixEvalSingle::gsMaterialMatrixEvalSingle
gsMaterialMatrixEvalSingle(index_t patch, gsMaterialMatrixBase< T > *materialMatrix, const gsFunctionSet< T > *deformed, const gsMatrix< T > z)
Constructor.
Definition: gsMaterialMatrixEval.hpp:24
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::Thickness, short_t >::type targetDim_impl() const
Implementation of targetDim for the thickness.
Definition: gsMaterialMatrixEval.h:258
gismo::gsMaterialMatrixEvalSingle::setPatch
void setPatch(index_t p)
Sets the patch index.
Definition: gsMaterialMatrixEval.h:270
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::PStress||_out==MaterialOutput::PStressN||_out==MaterialOutput::PCauchyStressN||_out==MaterialOutput::PStressM||_out==MaterialOutput::PCauchyStressM||_out==MaterialOutput::PStrainN||_out==MaterialOutput::PStrainM, short_t >::type targetDim_impl() const
Implementation of targetDim for principal stress fields.
Definition: gsMaterialMatrixEval.h:211
gismo::gsMaterialMatrixEvalSingle::domainDim
short_t domainDim() const
Domain dimension, always 2 for shells.
Definition: gsMaterialMatrixEval.h:167
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::StretchTransform||_out==MaterialOutput::PStressTransform, short_t >::type targetDim_impl() const
Implementation of targetDim for principal stress directions.
Definition: gsMaterialMatrixEval.h:235
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::Density, short_t >::type targetDim_impl() const
Implementation of targetDim for densities.
Definition: gsMaterialMatrixEval.h:182
gismo::gsFunctionSet
Interface for the set of functions defined on a domain (the total number of functions in the set equa...
Definition: gsFuncData.h:23
gismo::freeAll
void freeAll(It begin, It end)
Frees all pointers in the range [begin end)
Definition: gsMemory.h:312
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::Strain||_out==MaterialOutput::StrainN||_out==MaterialOutput::StrainM, short_t >::type targetDim_impl() const
Specialisation of targetDim for strain.
Definition: gsMaterialMatrixEval.h:246
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::Stretch, short_t >::type targetDim_impl() const
Implementation of targetDim for principal stretch fields.
Definition: gsMaterialMatrixEval.h:201
gismo::gsFunctionSet::nPieces
virtual index_t nPieces() const
Number of pieces in the domain of definition.
Definition: gsFunctionSet.h:584
gsMaterialMatrixContainer.h
Provides a container for material matrices for thin shells.
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::StretchDir, short_t >::type targetDim_impl() const
Implementation of targetDim for principal stretch directions.
Definition: gsMaterialMatrixEval.h:215
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::Deformation, short_t >::type targetDim_impl() const
Implementation of targetDim for principal stress directions.
Definition: gsMaterialMatrixEval.h:266
gsFunction.h
Provides declaration of Function abstract interface.
gismo::gsMaterialMatrixEvalSingle::targetDim
short_t targetDim() const
Target dimension.
Definition: gsMaterialMatrixEval.h:176
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::TensionField, short_t >::type targetDim_impl() const
Implementation of targetDim for the tension field indicator.
Definition: gsMaterialMatrixEval.h:230
gismo::gsMaterialMatrixEvalSingle::targetDim_impl
std::enable_if< _out==MaterialOutput::VectorN||_out==MaterialOutput::CauchyVectorN||_out==MaterialOutput::VectorM||_out==MaterialOutput::CauchyVectorM||_out==MaterialOutput::Generic, short_t >::type targetDim_impl() const
Implementation of targetDim for stress tensors.
Definition: gsMaterialMatrixEval.h:190