Source code for felupe._field._planestrain

# -*- coding: utf-8 -*-
"""
 _______  _______  ___      __   __  _______  _______ 
|       ||       ||   |    |  | |  ||       ||       |
|    ___||    ___||   |    |  | |  ||    _  ||    ___|
|   |___ |   |___ |   |    |  |_|  ||   |_| ||   |___ 
|    ___||    ___||   |___ |       ||    ___||    ___|
|   |    |   |___ |       ||       ||   |    |   |___ 
|___|    |_______||_______||_______||___|    |_______|

This file is part of felupe.

Felupe is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.

Felupe is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.

You should have received a copy of the GNU General Public License
along with Felupe.  If not, see <http://www.gnu.org/licenses/>.

"""

import numpy as np
from ._base import Field
from ..math import sym as symmetric


[docs]class FieldPlaneStrain(Field): r"""A plane strain Field on points of a two-dimensional `region` with dimension `dim` (default is 2) and initial point `values` (default is 0). This is a modified :class:`Field` class in which the radial coordinates are evaluated at the numeric integration points. The :meth:`grad`-method is modified in such a way that it does not only contain the in-plane 2d-gradient but also the circumferential stretch as shown in Eq.(1). .. code-block:: | dudX(2d) : 0 | dudX(planestrain) = | ..................| (1) | 0 : 0 | Attributes ---------- region : Region The region on which the field will be created. dim : int (default is 2) The dimension of the field. values : float (default is 0.0) or array A single value for all components of the field or an array of shape (region.mesh.npoints, dim)`. """ def __init__(self, region, dim=2, values=0): """A plane strain Field on points of a two-dimensional `region` with dimension `dim` (default is 2) and initial point `values` (default is 0). Attributes ---------- region : Region The region on which the field will be created. dim : int (default is 2) The dimension of the field. values : float (default is 0.0) or array A single value for all components of the field or an array of shape (region.mesh.npoints, dim)`. """ # init base Field super().__init__(region, dim=dim, values=values) def _interpolate_2d(self): """Interpolate 2D field values at points and evaluate them at the integration points of all cells in the region.""" # interpolated field values "aI" # evaluated at quadrature point "p" # for cell "c" return np.einsum( "ca...,apc->...pc", self.values[self.region.mesh.cells], self.region.h )
[docs] def interpolate(self): # extend dimension of in-plane 2d-gradient return np.pad(self._interpolate_2d(), ((0, 1), (0, 0), (0, 0)))
def _grad_2d(self, sym=False): """In-plane 2D gradient as partial derivative of field values at points w.r.t. the undeformed coordinates, evaluated at the integration points of all cells in the region. Optionally, the symmetric part of the gradient is returned. Arguments --------- sym : bool, optional (default is False) Calculate the symmetric part of the gradient. Returns ------- array In-plane 2D-gradient as partial derivative of field values at points w.r.t. undeformed coordinates, evaluated at the integration points of all cells in the region. """ # gradient as partial derivative of field component "I" at point "a" # w.r.t. undeformed coordinate "J" evaluated at quadrature point "p" # for each cell "e" g = np.einsum( "ca...,aJpc->...Jpc", self.values[self.region.mesh.cells], self.region.dhdX, ) if sym: return symmetric(g) else: return g
[docs] def grad(self, sym=False): """3D-gradient as partial derivative of field values at points w.r.t. the undeformed coordinates, evaluated at the integration points of all cells in the region. Optionally, the symmetric part of the gradient is returned. .. code-block:: | dudX(2d) : 0 | dudX(planestrain) = | ..................| | 0 : 0 | Arguments --------- sym : bool, optional (default is False) Calculate the symmetric part of the gradient. Returns ------- array Full 3D-gradient as partial derivative of field values at points w.r.t. undeformed coordinates, evaluated at the integration points of all cells in the region. """ # extend dimension of in-plane 2d-gradient g = np.pad(self._grad_2d(sym=sym), ((0, 1), (0, 1), (0, 0), (0, 0))) return g