G+Smo  25.01.0
Geometry + Simulation Modules
 
Loading...
Searching...
No Matches
gsMaterialMatrixBaseDim.h
Go to the documentation of this file.
1
16#pragma once
17
18#include <gsCore/gsFuncData.h>
19#include <gsUtils/gsThreaded.h>
20#include <gsKLShell/src/gsMaterialMatrixBase.h>
21
22namespace gismo
23{
24
25template <short_t dim, class T>
26class gsMaterialMatrixBaseDimData;
27
35template <short_t dim, class T>
36class gsMaterialMatrixBaseDim : public gsMaterialMatrixBase<T>
37{
38
39public:
40
41 typedef T Scalar_t;
42
43 using Base = gsMaterialMatrixBase<T>;
44
45 typedef typename Base::function_ptr function_ptr;
46
47 // enum {Linear=0}; // If true (1), this property entails that S = C *˙E
48
49 gsMaterialMatrixBaseDim()
50 :
51 Base(nullptr,nullptr,nullptr,nullptr)
52 {
53 this->defaultOptions();
54 membersSetZero();
55 }
56
57 gsMaterialMatrixBaseDim(const gsFunctionSet<T> * mp)
58 :
59 Base(mp,nullptr,nullptr,nullptr)
60 {
61 GISMO_ASSERT(mp->targetDim()==dim,"Geometric dimension and the template dimension are not the same!");
62 this->defaultOptions();
63 membersSetZero();
64 }
65
66 gsMaterialMatrixBaseDim(const gsFunctionSet<T> * mp,
67 const gsFunctionSet<T> * mp_def)
68 :
69 Base(mp,mp_def,nullptr,nullptr)
70 {
71 GISMO_ASSERT(mp->targetDim()==dim,"Geometric dimension and the template dimension are not the same!");
72 GISMO_ASSERT(mp_def->targetDim()==dim,"Geometric dimension and the template dimension are not the same!");
73 this->defaultOptions();
74 membersSetZero();
75 }
76
77 gsMaterialMatrixBaseDim(const gsFunctionSet<T> * mp,
78 const gsFunctionSet<T> * thickness,
79 const gsFunctionSet<T> * Density)
80 :
81 Base(mp,nullptr,thickness,Density)
82 {
83 GISMO_ASSERT(mp->targetDim()==dim,"Geometric dimension and the template dimension are not the same!");
84 this->defaultOptions();
85 membersSetZero();
86 }
87
88 gsMaterialMatrixBaseDim(const gsFunctionSet<T> * mp,
89 const gsFunctionSet<T> * mp_def,
90 const gsFunctionSet<T> * thickness,
91 const gsFunctionSet<T> * Density)
92 :
93 Base(mp,mp_def,thickness,Density)
94 {
95 GISMO_ASSERT(mp->targetDim()==dim,"Geometric dimension and the template dimension are not the same!");
96 GISMO_ASSERT(mp_def->targetDim()==dim,"Geometric dimension and the template dimension are not the same!");
97 this->defaultOptions();
98 membersSetZero();
99 }
100
102 virtual ~gsMaterialMatrixBaseDim() {}
103
104public:
106 virtual void defaultOptions() override;
107
109 virtual void density_into(const index_t patch, const gsMatrix<T>& u, gsMatrix<T>& result) const override;
110
112 virtual void thickness_into(const index_t patch, const gsMatrix<T>& u, gsMatrix<T>& result) const override;
113
115 virtual void parameters_into(const index_t patch, const gsMatrix<T> & u, gsMatrix<T>& result) const override;
116
118 virtual gsMatrix<T> eval3D_spec2cov(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z) const override;
119
121 virtual gsMatrix<T> eval3D_spec2con(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z) const override;
122
124 virtual gsMatrix<T> eval3D_cov2cart(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z) const override;
125
127 virtual gsMatrix<T> eval3D_con2cart(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z) const override;
128
130 virtual gsMatrix<T> eval3D_deformation(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z) const override;
131
133 virtual gsMatrix<T> eval3D_strain(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z) const override;
134 using Base::eval3D_strain;
135
137 virtual gsMatrix<T> eval3D_tensionfield(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z, enum MaterialOutput out) const override;
138 using Base::eval3D_tensionfield;
139
141 virtual gsMatrix<T> eval3D_pstretch(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z) const override;
142
144 virtual gsMatrix<T> eval3D_pstretchDir(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z) const override;
145
147 virtual gsMatrix<T> eval3D_pstrain(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z) const override;
148
149 virtual bool initialized() const override
150 {
151 return ((m_thickness != nullptr)); // && (m_patches!= nullptr) //HV: could be defined without geometry, then call "setUndeformed"
152 }
153
154public:
155
156 void _computePoints(const index_t patch, const gsMatrix<T> & u) const;
157
159 gsMatrix<T> _getAcov_def(index_t k, T z) const;
160
162 gsMatrix<T> _getAcon_def(index_t k, T z) const;
163
165 gsMatrix<T> _getBcov_def(index_t k, T z) const;
166
168 gsMatrix<T> _getncov_def(index_t k, T z) const;
169
171 gsMatrix<T> _getGcov_def(index_t k, T z) const;
172
174 gsMatrix<T> _getGcon_def(index_t k, T z) const;
175
177 gsMatrix<T> _getacov_def(index_t k, T z) const;
178
180 gsMatrix<T> _getacon_def(index_t k, T z) const;
181
183 gsMatrix<T> _getgcov_def(index_t k, T z) const;
184
186 gsMatrix<T> _getgcon_def(index_t k, T z) const;
187
189 gsMatrix<T> _getAcov_ori(index_t k, T z) const;
190
192 gsMatrix<T> _getAcon_ori(index_t k, T z) const;
193
195 gsMatrix<T> _getBcov_ori(index_t k, T z) const;
196
198 gsMatrix<T> _getncov_ori(index_t k, T z) const;
199
201 gsMatrix<T> _getGcov_ori(index_t k, T z) const;
202
204 gsMatrix<T> _getGcon_ori(index_t k, T z) const;
205
207 gsMatrix<T> _getacov_ori(index_t k, T z) const;
208
210 gsMatrix<T> _getacon_ori(index_t k, T z) const;
211
213 gsMatrix<T> _getgcov_ori(index_t k, T z) const;
214
216 gsMatrix<T> _getgcon_ori(index_t k, T z) const;
217
219 void _computeMetricDeformed(const index_t patch, const gsMatrix<T> & u) const;
220
222 void _computeMetricUndeformed(const index_t patch, const gsMatrix<T> & u) const;
223
225 void _getMetric(const index_t k, const T z) const;
226
228 void _getMetric(index_t k, T z, const gsMatrix<T> & C) const;
229
231 void _getMetricDeformed(const index_t k, const T z) const;
232
233 void _getMetricDeformed(const gsMatrix<T> & C) const;
234
236 void _getMetricUndeformed(const index_t k, const T z) const;
237
239 std::pair<gsVector<T>,gsMatrix<T>> _evalStretch(const gsMatrix<T> & C, const gsMatrix<T> & gcon_ori ) const;
240
242 std::pair<gsVector<T>,gsMatrix<T>> _evalPStrain(const gsMatrix<T> & C ) const;
243
245 std::pair<gsVector<T>,gsMatrix<T>> _evalPStress(const gsMatrix<T> & S ) const;
246
248 void _computeStretch(const gsMatrix<T> & C, const gsMatrix<T> & gcon_ori ) const;
249
251 void _computePStrain(const gsMatrix<T> & C ) const;
252
254 void _computePStress(const gsMatrix<T> & C ) const;
255
256
258 gsMatrix<T> _transformation(const gsMatrix<T> & basis1, const gsMatrix<T> & basis2 ) const;
259
260 virtual void setUndeformed(const gsFunctionSet<T> * undeformed) override
261 {
262 GISMO_ASSERT(undeformed->targetDim()==dim,"Geometric dimension and the template dimension are not the same!");
263 function_ptr f_ptr = memory::make_shared_not_owned(undeformed);
264 Base::setUndeformed(f_ptr);
265 }
266
267 virtual void setDeformed(const gsFunctionSet<T> * deformed) override
268 {
269 GISMO_ASSERT(deformed->targetDim()==dim,"Geometric dimension and the template dimension are not the same!");
270 function_ptr f_ptr = memory::make_shared_not_owned(deformed);
271 Base::setDeformed(f_ptr);
272 }
273
274 virtual void setUndeformed(const function_ptr undeformed) override
275 {
276 GISMO_ASSERT(undeformed->targetDim()==dim,"Geometric dimension and the template dimension are not the same!");
277 Base::setUndeformed(undeformed);
278 }
279
280 virtual void setDeformed(const function_ptr deformed) override
281 {
282 GISMO_ASSERT(deformed->targetDim()==dim,"Geometric dimension and the template dimension are not the same!");
283 Base::setDeformed(deformed);
284 }
285
286private:
287
288 template <short_t _type>
289 typename std::enable_if<_type==0, index_t>::type _tensionField(const gsVector<T> &, const gsVector<T> & Ep) const;
290
291 template <short_t _type>
292 typename std::enable_if<_type==1, index_t>::type _tensionField(const gsVector<T> &, const gsVector<T> & Ep) const;
293
294 template <short_t _type>
295 typename std::enable_if<_type==2, index_t>::type _tensionField(const gsVector<T> &, const gsVector<T> & Ep) const;
296
297 template <short_t _type>
298 typename std::enable_if<_type!=0 &&
299 _type!=1 &&
300 _type!=2, index_t>::type _tensionField(const gsVector<T> &, const gsVector<T> & Ep) const
301 {GISMO_NO_IMPLEMENTATION;}
302
304 template<short_t _dim>
305 typename std::enable_if<_dim==2, void>::type _computeMetricDeformed_impl(const index_t patch, const gsMatrix<T> & u) const;
306
308 template<short_t _dim>
309 typename std::enable_if<_dim==3, void>::type _computeMetricDeformed_impl(const index_t patch, const gsMatrix<T> & u) const;
310
312 template<short_t _dim>
313 typename std::enable_if<_dim==2, void>::type _computeMetricUndeformed_impl(const index_t patch, const gsMatrix<T> & u) const;
314
316 template<short_t _dim>
317 typename std::enable_if<_dim==3, void>::type _computeMetricUndeformed_impl(const index_t patch, const gsMatrix<T> & u) const;
318
320 template<short_t _dim>
321 typename std::enable_if<_dim==2, void>::type _getMetric_impl(const index_t k, const T z) const;
322
323 template<short_t _dim>
324 typename std::enable_if<_dim==3, void>::type _getMetric_impl(const index_t k, const T z) const;
325
327 template<short_t _dim>
328 typename std::enable_if<_dim==2, void>::type _getMetricDeformed_impl(const index_t k, const T z) const;
329
331 template<short_t _dim>
332 typename std::enable_if<_dim==3, void>::type _getMetricDeformed_impl(const index_t k, const T z) const;
333
335 template<short_t _dim>
336 typename std::enable_if<_dim==2, void>::type _getMetricUndeformed_impl(const index_t k, const T z) const;
337
339 template<short_t _dim>
340 typename std::enable_if<_dim==3, void>::type _getMetricUndeformed_impl(const index_t k, const T z) const;
341
342 // -------------------------------------------------------------------------------------------------------------
343 // Separate getters
344 // -------------------------------------------------------------------------------------------------------------
345
347 template<short_t _dim>
348 typename std::enable_if<_dim==2, gsMatrix<T>>::type _getBcov_def_impl(index_t k, T z) const;
349
351 template<short_t _dim>
352 typename std::enable_if<_dim==3, gsMatrix<T>>::type _getBcov_def_impl(index_t k, T z) const;
353
355 template<short_t _dim>
356 typename std::enable_if<_dim==2, gsMatrix<T>>::type _getncov_def_impl(index_t k, T z) const;
357
359 template<short_t _dim>
360 typename std::enable_if<_dim==3, gsMatrix<T>>::type _getncov_def_impl(index_t k, T z) const;
361
363 template<short_t _dim>
364 typename std::enable_if<_dim==2, gsMatrix<T>>::type _getGcov_def_impl(index_t k, T z) const;
365
367 template<short_t _dim>
368 typename std::enable_if<_dim==3, gsMatrix<T>>::type _getGcov_def_impl(index_t k, T z) const;
369
371 template<short_t _dim>
372 typename std::enable_if<_dim==2, gsMatrix<T>>::type _getGcon_def_impl(index_t k, T z) const;
373
375 template<short_t _dim>
376 typename std::enable_if<_dim==3, gsMatrix<T>>::type _getGcon_def_impl(index_t k, T z) const;
377
379 template<short_t _dim>
380 typename std::enable_if<_dim==2, gsMatrix<T>>::type _getgcov_def_impl(index_t k, T z) const;
381
383 template<short_t _dim>
384 typename std::enable_if<_dim==3, gsMatrix<T>>::type _getgcov_def_impl(index_t k, T z) const;
385
387 template<short_t _dim>
388 typename std::enable_if<_dim==2, gsMatrix<T>>::type _getBcov_ori_impl(index_t k, T z) const;
389
391 template<short_t _dim>
392 typename std::enable_if<_dim==3, gsMatrix<T>>::type _getBcov_ori_impl(index_t k, T z) const;
393
395 template<short_t _dim>
396 typename std::enable_if<_dim==2, gsMatrix<T>>::type _getncov_ori_impl(index_t k, T z) const;
397
399 template<short_t _dim>
400 typename std::enable_if<_dim==3, gsMatrix<T>>::type _getncov_ori_impl(index_t k, T z) const;
401
403 template<short_t _dim>
404 typename std::enable_if<_dim==2, gsMatrix<T>>::type _getGcov_ori_impl(index_t k, T z) const;
405
407 template<short_t _dim>
408 typename std::enable_if<_dim==3, gsMatrix<T>>::type _getGcov_ori_impl(index_t k, T z) const;
409
411 template<short_t _dim>
412 typename std::enable_if<_dim==2, gsMatrix<T>>::type _getGcon_ori_impl(index_t k, T z) const;
413
415 template<short_t _dim>
416 typename std::enable_if<_dim==3, gsMatrix<T>>::type _getGcon_ori_impl(index_t k, T z) const;
417
419 template<short_t _dim>
420 typename std::enable_if<_dim==2, gsMatrix<T>>::type _getgcov_ori_impl(index_t k, T z) const;
421
423 template<short_t _dim>
424 typename std::enable_if<_dim==3, gsMatrix<T>>::type _getgcov_ori_impl(index_t k, T z) const;
425
426
427protected:
428
429 void membersSetZero()
430 {
431 m_data.mine().membersSetZero();
432 }
433
434 using Base::m_patches;
435 using Base::m_defpatches;
436
437 using Base::m_options;
438
439 using Base::m_pars;
440 using Base::m_thickness;
441 using Base::m_density;
442
443 // Geometric data point
444 mutable util::gsThreaded< gsMaterialMatrixBaseDimData<dim,T> > m_data;
445public:
446
447# define Eigen gsEigen
448 EIGEN_MAKE_ALIGNED_OPERATOR_NEW //must be present whenever the class contains fixed size matrices
449# undef Eigen
450};
451
452template<short_t dim, class T>
453class gsMaterialMatrixBaseDimData
454{
455public:
456
457 void membersSetZero()
458 {
459 m_Acov_ori.setZero(); m_Acon_ori.setZero(); m_Acov_def.setZero(); m_Acon_def.setZero(); m_Bcov_ori.setZero(); m_Bcon_ori.setZero(); m_Bcov_def.setZero(); m_Bcon_def.setZero();
460 m_acov_ori.setZero(); m_acon_ori.setZero(); m_acov_def.setZero(); m_acon_def.setZero();
461 m_ncov_ori.setZero(); m_ncov_def.setZero();
462 m_Gcov_ori.setZero(); m_Gcon_ori.setZero(); m_Gcov_def.setZero(); m_Gcon_def.setZero(); m_Gcov_ori_L.setZero(); m_Gcov_def_L.setZero();
463 m_gcov_ori.setZero(); m_gcov_def.setZero();m_gcon_ori.setZero(); m_gcon_def.setZero();
464 m_Acov_ori_mat.setZero(); m_Acon_ori_mat.setZero(); m_Acov_def_mat.setZero(); m_Acon_def_mat.setZero(); m_Bcov_ori_mat.setZero(); m_Bcov_def_mat.setZero();
465 m_acov_ori_mat.setZero(); m_acon_ori_mat.setZero(); m_acov_def_mat.setZero(); m_acon_def_mat.setZero(); m_ncov_ori_mat.setZero(); m_ncov_def_mat.setZero(); m_normal_ori_mat.setZero(); m_normal_def_mat.setZero();
466
467 }
468 // Material parameters and kinematics
469 mutable gsMatrix<T> m_parmat;
470 mutable gsVector<T> m_parvals;
471 mutable gsMatrix<T> m_Tmat,m_rhomat;
472
473 mutable gsMatrix<T> m_Acov_ori, m_Acon_ori, m_Acov_def, m_Acon_def, m_Bcov_ori, m_Bcon_ori, m_Bcov_def, m_Bcon_def;
474 mutable gsMatrix<T> m_acov_ori, m_acon_ori, m_acov_def, m_acon_def;
475 mutable gsMatrix<T> m_ncov_ori, m_ncov_def;
476 mutable gsMatrix<T> m_Gcov_ori, m_Gcon_ori, m_Gcov_def, m_Gcon_def, m_Gcov_ori_L, m_Gcov_def_L;
477 mutable gsMatrix<T> m_gcov_ori, m_gcov_def,m_gcon_ori, m_gcon_def;
478 mutable gsMatrix<T> m_Acov_ori_mat, m_Acon_ori_mat, m_Acov_def_mat, m_Acon_def_mat, m_Bcov_ori_mat, m_Bcov_def_mat;
479 mutable gsMatrix<T> m_acov_ori_mat, m_acon_ori_mat, m_acov_def_mat, m_acon_def_mat, m_ncov_ori_mat, m_ncov_def_mat, m_normal_ori_mat, m_normal_def_mat;
480
481 mutable gsMatrix<T> m_stretches, m_stretchvec, m_pstress, m_pstressvec, m_pstrain, m_pstrainvec;
482
483 mutable T m_J0_sq, m_J_sq;
484
485 mutable gsVector<T> m_thetas, m_gammas;
486
487# define Eigen gsEigen
488 EIGEN_MAKE_ALIGNED_OPERATOR_NEW
489# undef Eigen //must be present whenever the class contains fixed size matrices
490};
491
492#ifdef GISMO_WITH_PYBIND11
493
497 void pybind11_init_gsMaterialMatrixBaseDim2(pybind11::module &m);
498 void pybind11_init_gsMaterialMatrixBaseDim3(pybind11::module &m);
499
500#endif // GISMO_WITH_PYBIND11
501
502} // namespace
503
504#ifndef GISMO_BUILD_LIB
505#include GISMO_HPP_HEADER(gsMaterialMatrixBaseDim.hpp)
506#endif
gismo::gsFunctionSet
Interface for the set of functions defined on a domain (the total number of functions in the set equa...
Definition gsFunctionSet.h:219
gismo::gsFunctionSet::targetDim
virtual short_t targetDim() const
Dimension of the target space.
Definition gsFunctionSet.h:595
gismo::gsMaterialMatrixBaseDim
This class defines the base class for material matrices.
Definition gsMaterialMatrixBaseDim.h:37
gismo::gsMaterialMatrixBaseDim::_computeMetricUndeformed
void _computeMetricUndeformed(const index_t patch, const gsMatrix< T > &u) const
Computes metric quantities on the undeformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:535
gismo::gsMaterialMatrixBaseDim::eval3D_strain
virtual gsMatrix< T > eval3D_strain(const index_t patch, const gsMatrix< T > &u, const gsMatrix< T > &z) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:271
gismo::gsMaterialMatrixBaseDim::_transformation
gsMatrix< T > _transformation(const gsMatrix< T > &basis1, const gsMatrix< T > &basis2) const
Computes the stretch given deformation tensor C, into class members m_stretches and m_stretchDirs.
Definition gsMaterialMatrixBaseDim.hpp:1502
gismo::gsMaterialMatrixBaseDim::_computePStrain
void _computePStrain(const gsMatrix< T > &C) const
Computes the stretch given deformation tensor C, into class members m_stretches and m_stretchDirs.
Definition gsMaterialMatrixBaseDim.hpp:1492
gismo::gsMaterialMatrixBaseDim::_getncov_def
gsMatrix< T > _getncov_def(index_t k, T z) const
Returns the covariant n tensor on the deformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:671
gismo::gsMaterialMatrixBaseDim::_computeMetricUndeformed_impl
std::enable_if< _dim==2, void >::type _computeMetricUndeformed_impl(const index_t patch, const gsMatrix< T > &u) const
Implementation of _getMetric for planar geometries.
Definition gsMaterialMatrixBaseDim.hpp:543
gismo::gsMaterialMatrixBaseDim::_getGcov_ori
gsMatrix< T > _getGcov_ori(index_t k, T z) const
Returns the covariant metric tensor on the original geometry.
Definition gsMaterialMatrixBaseDim.hpp:751
gismo::gsMaterialMatrixBaseDim::_computeMetricDeformed_impl
std::enable_if< _dim==3, void >::type _computeMetricDeformed_impl(const index_t patch, const gsMatrix< T > &u) const
Implementation of _computeMetricUndeformed for surface geometries.
gismo::gsMaterialMatrixBaseDim::_getMetric_impl
std::enable_if< _dim==2, void >::type _getMetric_impl(const index_t k, const T z) const
Implementation of _getMetricDeformed for planar geometries.
gismo::gsMaterialMatrixBaseDim::_getGcov_def
gsMatrix< T > _getGcov_def(index_t k, T z) const
Returns the covariant metric tensor on the deformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:677
gismo::gsMaterialMatrixBaseDim::eval3D_spec2con
virtual gsMatrix< T > eval3D_spec2con(const index_t patch, const gsMatrix< T > &u, const gsMatrix< T > &z) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:90
gismo::gsMaterialMatrixBaseDim::_getBcov_def
gsMatrix< T > _getBcov_def(index_t k, T z) const
Returns the covariant b tensor on the deformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:665
gismo::gsMaterialMatrixBaseDim::_getMetricDeformed
void _getMetricDeformed(const index_t k, const T z) const
Gets metric quantities on the deformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:1073
gismo::gsMaterialMatrixBaseDim::_computeMetricDeformed_impl
std::enable_if< _dim==2, void >::type _computeMetricDeformed_impl(const index_t patch, const gsMatrix< T > &u) const
Implementation of _computeMetricUndeformed for planar geometries.
Definition gsMaterialMatrixBaseDim.hpp:425
gismo::gsMaterialMatrixBaseDim::eval3D_tensionfield
virtual gsMatrix< T > eval3D_tensionfield(const index_t patch, const gsMatrix< T > &u, const gsMatrix< T > &z, enum MaterialOutput out) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:296
gismo::gsMaterialMatrixBaseDim::density_into
virtual void density_into(const index_t patch, const gsMatrix< T > &u, gsMatrix< T > &result) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:40
gismo::gsMaterialMatrixBaseDim::eval3D_pstretch
virtual gsMatrix< T > eval3D_pstretch(const index_t patch, const gsMatrix< T > &u, const gsMatrix< T > &z) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:200
gismo::gsMaterialMatrixBaseDim::_getBcov_ori_impl
std::enable_if< _dim==2, gsMatrix< T > >::type _getBcov_ori_impl(index_t k, T z) const
Implementation of _getMetricUndeformed for planar geometries.
Definition gsMaterialMatrixBaseDim.hpp:901
gismo::gsMaterialMatrixBaseDim::_getGcon_def
gsMatrix< T > _getGcon_def(index_t k, T z) const
Returns the contravariant metric tensor on the deformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:683
gismo::gsMaterialMatrixBaseDim::_getMetricDeformed_impl
std::enable_if< _dim==2, void >::type _getMetricDeformed_impl(const index_t k, const T z) const
Implementation of _getMetricDeformed for surface geometries.
Definition gsMaterialMatrixBaseDim.hpp:1081
gismo::gsMaterialMatrixBaseDim::eval3D_pstrain
virtual gsMatrix< T > eval3D_pstrain(const index_t patch, const gsMatrix< T > &u, const gsMatrix< T > &z) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:244
gismo::gsMaterialMatrixBaseDim::_getAcon_ori
gsMatrix< T > _getAcon_ori(index_t k, T z) const
Returns the contravariant a tensor on the original geometry.
Definition gsMaterialMatrixBaseDim.hpp:732
gismo::gsMaterialMatrixBaseDim::eval3D_con2cart
virtual gsMatrix< T > eval3D_con2cart(const index_t patch, const gsMatrix< T > &u, const gsMatrix< T > &z) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:135
gismo::gsMaterialMatrixBaseDim::_getgcon_def
gsMatrix< T > _getgcon_def(index_t k, T z) const
Returns the contravariant basis vector g on the deformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:710
gismo::gsMaterialMatrixBaseDim::_getacov_def
gsMatrix< T > _getacov_def(index_t k, T z) const
Returns the covariant basis vector a on the deformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:690
gismo::gsMaterialMatrixBaseDim::_getgcon_ori
gsMatrix< T > _getgcon_ori(index_t k, T z) const
Returns the contravariant basis vector g on the original geometry.
Definition gsMaterialMatrixBaseDim.hpp:784
gismo::gsMaterialMatrixBaseDim::_computeMetricDeformed
void _computeMetricDeformed(const index_t patch, const gsMatrix< T > &u) const
Computes metric quantities on the deformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:417
gismo::gsMaterialMatrixBaseDim::_getBcov_def_impl
std::enable_if< _dim==2, gsMatrix< T > >::type _getBcov_def_impl(index_t k, T z) const
Implementation of _getMetricUndeformed for planar geometries.
Definition gsMaterialMatrixBaseDim.hpp:803
gismo::gsMaterialMatrixBaseDim::_getgcov_def_impl
std::enable_if< _dim==2, gsMatrix< T > >::type _getgcov_def_impl(index_t k, T z) const
Implementation of _getMetricUndeformed for planar geometries.
Definition gsMaterialMatrixBaseDim.hpp:871
gismo::gsMaterialMatrixBaseDim::_getBcov_ori
gsMatrix< T > _getBcov_ori(index_t k, T z) const
Returns the covariant b tensor on the original geometry.
Definition gsMaterialMatrixBaseDim.hpp:739
gismo::gsMaterialMatrixBaseDim::_getacon_ori
gsMatrix< T > _getacon_ori(index_t k, T z) const
Returns the contravariant basis vector a on the original geometry.
Definition gsMaterialMatrixBaseDim.hpp:771
gismo::gsMaterialMatrixBaseDim::defaultOptions
virtual void defaultOptions() override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:34
gismo::gsMaterialMatrixBaseDim::_computeStretch
void _computeStretch(const gsMatrix< T > &C, const gsMatrix< T > &gcon_ori) const
Computes the stretch given deformation tensor C, into class members m_stretches and m_stretchDirs.
Definition gsMaterialMatrixBaseDim.hpp:1476
gismo::gsMaterialMatrixBaseDim::_getacov_ori
gsMatrix< T > _getacov_ori(index_t k, T z) const
Returns the covariant basis vector a on the original geometry.
Definition gsMaterialMatrixBaseDim.hpp:764
gismo::gsMaterialMatrixBaseDim::_evalStretch
std::pair< gsVector< T >, gsMatrix< T > > _evalStretch(const gsMatrix< T > &C, const gsMatrix< T > &gcon_ori) const
Computes the stretch given deformation tensor C, into a pair.
Definition gsMaterialMatrixBaseDim.hpp:1331
gismo::gsMaterialMatrixBaseDim::_computePStress
void _computePStress(const gsMatrix< T > &C) const
Computes the principal stresses of a given stress tensor S, into class members m_pstress and m_pstres...
Definition gsMaterialMatrixBaseDim.hpp:1484
gismo::gsMaterialMatrixBaseDim::eval3D_pstretchDir
virtual gsMatrix< T > eval3D_pstretchDir(const index_t patch, const gsMatrix< T > &u, const gsMatrix< T > &z) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:222
gismo::gsMaterialMatrixBaseDim::_getAcov_def
gsMatrix< T > _getAcov_def(index_t k, T z) const
Returns the covariant a tensor on the deformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:651
gismo::gsMaterialMatrixBaseDim::_getgcov_ori
gsMatrix< T > _getgcov_ori(index_t k, T z) const
Returns the covariant basis vector g on the original geometry.
Definition gsMaterialMatrixBaseDim.hpp:778
gismo::gsMaterialMatrixBaseDim::~gsMaterialMatrixBaseDim
virtual ~gsMaterialMatrixBaseDim()
Destructor.
Definition gsMaterialMatrixBaseDim.h:102
gismo::gsMaterialMatrixBaseDim::_getgcov_def
gsMatrix< T > _getgcov_def(index_t k, T z) const
Returns the covariant basis vector g on the deformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:704
gismo::gsMaterialMatrixBaseDim::_getGcon_ori
gsMatrix< T > _getGcon_ori(index_t k, T z) const
Returns the contravariant metric tensor on the original geometry.
Definition gsMaterialMatrixBaseDim.hpp:757
gismo::gsMaterialMatrixBaseDim::_evalPStrain
std::pair< gsVector< T >, gsMatrix< T > > _evalPStrain(const gsMatrix< T > &C) const
Computes the principal strain given deformation tensor C, into a pair.
Definition gsMaterialMatrixBaseDim.hpp:1428
gismo::gsMaterialMatrixBaseDim::_getGcov_ori_impl
std::enable_if< _dim==2, gsMatrix< T > >::type _getGcov_ori_impl(index_t k, T z) const
Implementation of _getMetricUndeformed for planar geometries.
Definition gsMaterialMatrixBaseDim.hpp:935
gismo::gsMaterialMatrixBaseDim::thickness_into
virtual void thickness_into(const index_t patch, const gsMatrix< T > &u, gsMatrix< T > &result) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:57
gismo::gsMaterialMatrixBaseDim::_getGcon_def_impl
std::enable_if< _dim==2, gsMatrix< T > >::type _getGcon_def_impl(index_t k, T z) const
Implementation of _getMetricUndeformed for planar geometries.
gismo::gsMaterialMatrixBaseDim::_getncov_def_impl
std::enable_if< _dim==2, gsMatrix< T > >::type _getncov_def_impl(index_t k, T z) const
Implementation of _getMetricUndeformed for planar geometries.
Definition gsMaterialMatrixBaseDim.hpp:820
gismo::gsMaterialMatrixBaseDim::_getAcov_ori
gsMatrix< T > _getAcov_ori(index_t k, T z) const
Returns the covariant a tensor on the original geometry.
Definition gsMaterialMatrixBaseDim.hpp:725
gismo::gsMaterialMatrixBaseDim::_getMetricUndeformed_impl
std::enable_if< _dim==3, void >::type _getMetricUndeformed_impl(const index_t k, const T z) const
Implementation of _getMetricUndeformed for surface geometries.
gismo::gsMaterialMatrixBaseDim::eval3D_cov2cart
virtual gsMatrix< T > eval3D_cov2cart(const index_t patch, const gsMatrix< T > &u, const gsMatrix< T > &z) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:112
gismo::gsMaterialMatrixBaseDim::_getGcon_ori_impl
std::enable_if< _dim==2, gsMatrix< T > >::type _getGcon_ori_impl(index_t k, T z) const
Implementation of _getMetricUndeformed for planar geometries.
gismo::gsMaterialMatrixBaseDim::_evalPStress
std::pair< gsVector< T >, gsMatrix< T > > _evalPStress(const gsMatrix< T > &S) const
Computes the principal stress given stress tensor S, into a pair.
Definition gsMaterialMatrixBaseDim.hpp:1384
gismo::gsMaterialMatrixBaseDim::_getGcon_def_impl
std::enable_if< _dim==3, gsMatrix< T > >::type _getGcon_def_impl(index_t k, T z) const
Implementation of _getMetricUndeformed for surface geometries.
gismo::gsMaterialMatrixBaseDim::_getMetricUndeformed_impl
std::enable_if< _dim==2, void >::type _getMetricUndeformed_impl(const index_t k, const T z) const
Implementation of _getMetricUndeformed for planar geometries.
Definition gsMaterialMatrixBaseDim.hpp:1210
gismo::gsMaterialMatrixBaseDim::_getncov_ori
gsMatrix< T > _getncov_ori(index_t k, T z) const
Returns the covariant n tensor on the original geometry.
Definition gsMaterialMatrixBaseDim.hpp:745
gismo::gsMaterialMatrixBaseDim::_computeMetricUndeformed_impl
std::enable_if< _dim==3, void >::type _computeMetricUndeformed_impl(const index_t patch, const gsMatrix< T > &u) const
Implementation of _getMetric for surface geometries.
gismo::gsMaterialMatrixBaseDim::_getgcov_ori_impl
std::enable_if< _dim==2, gsMatrix< T > >::type _getgcov_ori_impl(index_t k, T z) const
Implementation of _getMetricUndeformed for planar geometries.
Definition gsMaterialMatrixBaseDim.hpp:969
gismo::gsMaterialMatrixBaseDim::eval3D_spec2cov
virtual gsMatrix< T > eval3D_spec2cov(const index_t patch, const gsMatrix< T > &u, const gsMatrix< T > &z) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:68
gismo::gsMaterialMatrixBaseDim::_getMetricUndeformed
void _getMetricUndeformed(const index_t k, const T z) const
Gets metric quantities on the undeformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:1202
gismo::gsMaterialMatrixBaseDim::_getMetricDeformed_impl
std::enable_if< _dim==3, void >::type _getMetricDeformed_impl(const index_t k, const T z) const
Implementation of _getMetricDeformed for surface geometries.
gismo::gsMaterialMatrixBaseDim::_getGcov_def_impl
std::enable_if< _dim==2, gsMatrix< T > >::type _getGcov_def_impl(index_t k, T z) const
Implementation of _getMetricUndeformed for planar geometries.
Definition gsMaterialMatrixBaseDim.hpp:837
gismo::gsMaterialMatrixBaseDim::parameters_into
virtual void parameters_into(const index_t patch, const gsMatrix< T > &u, gsMatrix< T > &result) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:157
gismo::gsMaterialMatrixBaseDim::_getAcon_def
gsMatrix< T > _getAcon_def(index_t k, T z) const
Returns the contravariant a tensor on the deformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:658
gismo::gsMaterialMatrixBaseDim::_getMetric
void _getMetric(const index_t k, const T z) const
Gets metric quantities on the deformed and undeformed geometries.
Definition gsMaterialMatrixBaseDim.hpp:1021
gismo::gsMaterialMatrixBaseDim::_getncov_ori_impl
std::enable_if< _dim==2, gsMatrix< T > >::type _getncov_ori_impl(index_t k, T z) const
Implementation of _getMetricUndeformed for planar geometries.
Definition gsMaterialMatrixBaseDim.hpp:918
gismo::gsMaterialMatrixBaseDim::_getGcon_ori_impl
std::enable_if< _dim==3, gsMatrix< T > >::type _getGcon_ori_impl(index_t k, T z) const
Implementation of _getMetricUndeformed for surface geometries.
gismo::gsMaterialMatrixBaseDim::eval3D_deformation
virtual gsMatrix< T > eval3D_deformation(const index_t patch, const gsMatrix< T > &u, const gsMatrix< T > &z) const override
See gsMaterialMatrixBase for details.
Definition gsMaterialMatrixBaseDim.hpp:176
gismo::gsMaterialMatrixBaseDim::_getacon_def
gsMatrix< T > _getacon_def(index_t k, T z) const
Returns the contravariant basis vector a on the deformed geometry.
Definition gsMaterialMatrixBaseDim.hpp:697
gismo::gsMaterialMatrixBase
This class defines the base class for material matrices.
Definition gsMaterialMatrixBase.h:33
gismo::gsMaterialMatrixBase::eval3D_strain
virtual gsMatrix< T > eval3D_strain(const index_t, const gsMatrix< T > &, const gsMatrix< T > &) const
{ function_description }
Definition gsMaterialMatrixBase.h:530
gismo::gsMaterialMatrixBase::eval3D_tensionfield
virtual gsMatrix< T > eval3D_tensionfield(const index_t, const gsMatrix< T > &, const gsMatrix< T > &, enum MaterialOutput) const
to do
Definition gsMaterialMatrixBase.h:601
gismo::gsMatrix
A matrix with arbitrary coefficient type and fixed or dynamic size.
Definition gsMatrix.h:41
gismo::MaterialOutput
MaterialOutput
This class describes the output type.
Definition gsMaterialMatrixUtils.h:99
index_t
#define index_t
Definition gsConfig.h:32
GISMO_NO_IMPLEMENTATION
#define GISMO_NO_IMPLEMENTATION
Definition gsDebug.h:129
GISMO_ASSERT
#define GISMO_ASSERT(cond, message)
Definition gsDebug.h:89
gsFuncData.h
This object is a cache for computed values from an evaluator.
gsMaterialMatrixBase.h
Provides a base class for material matrices.
gsThreaded.h
Wrapper for thread-local data members.
gismo::memory::make_shared_not_owned
shared_ptr< T > make_shared_not_owned(const T *x)
Creates a shared pointer which does not eventually delete the underlying raw pointer....
Definition gsMemory.h:189
gismo
The G+Smo namespace, containing all definitions for the library.