Orthotropy in MFront solid material models #2774
Conversation
| OGSToMFront(Q.transpose() | ||
| .template topLeftCorner< | ||
| KelvinVectorDimensions<DisplacementDim>::value, | ||
| KelvinVectorDimensions<DisplacementDim>::value>() * |
There was a problem hiding this comment.
is the topLeftCorner stuff needed now that we have a 2D (4x4) variant of Q?
| sigma = Q.template topLeftCorner< | ||
| KelvinVectorDimensions<DisplacementDim>::value, | ||
| KelvinVectorDimensions<DisplacementDim>::value>() * | ||
| MFrontToOGS(sigma); |
There was a problem hiding this comment.
We rotate sigma back into OGS coordinate system, so epsilon and sigma will be output in xyz. All MFront internal variables (ElasticStrain for example) will be output in the local coordinate system. in the _behaviour.symmetry != mgis::behaviour::Behaviour::Symmetry::ISOTROPIC case we could output a warning that his is the case and users have to be careful. Having the local coordinate system as an output in ParaView for use in transformations would be convenient (but not necessary).
There was a problem hiding this comment.
A todo at the place where the rotation should happen added.
With the assumption of e3 = z-axis being the third basis vector in the 2D case.
The tensors are rotated before being passed to MFront. The results of MFront including the stiffness tensor.
Backed up by OGS orthotropic elasticity tests. The results are copied from the OGS' references.
Somehow the template deduction of MSVC could not handle the additional parameters with default values.
|
OpenGeoSys development has been moved to GitLab. |
Orthotropy in MFront solid material models 2D and 3D cases.
Tested against OGS orthotropic elasticity implementation.