solidmechanics/plate/mindlin_reissner/model.hh File Reference

Mindlin-Reissner plate model. More...

#include <dumux/common/properties.hh>
#include <dumux/common/properties/model.hh>
#include <dumux/solidmechanics/plate/kirchhoff_love/model.hh>
#include <dumux/solidmechanics/plate/kirchhoff_love/volumevariables.hh>
#include <dumux/solidmechanics/plate/kirchhoff_love/couplingmanager.hh>
#include "localresidual.hh"

Go to the source code of this file.

Description

Mindlin-Reissner plate model for small strains and small rotations. The stress resultants are defined by integration over the plate thickness \( [-\tfrac{t}{2}, \tfrac{t}{2}] \):

\[ N_{\alpha\beta} := \int_{-t/2}^{t/2}\sigma_{\alpha\beta}\,dx_3,\quad M_{\alpha\beta} := \int_{-t/2}^{t/2} x_3\,\sigma_{\alpha\beta}\,dx_3,\quad Q_\alpha := \int_{-t/2}^{t/2}\sigma_{\alpha 3}\,dx_3. \]

Ignoring in-plane stresses, the equilibrium equations are

\begin{align} \nabla\cdot\mathbf{Q} &= F,\\ -\nabla\cdot\mathbf{M} + \mathbf{Q} &= \mathbf{0}, \end{align}

where \( F \) is the out-of-plane load. Using a plane-stress constitutive law, the moment resultants for an isotropic material are

\[ \mathbf{M}(\boldsymbol{\theta}) = -D\left\{(1-\nu)\,\boldsymbol{\varepsilon}(\boldsymbol{\theta}) + \nu\operatorname{tr}(\boldsymbol{\varepsilon}(\boldsymbol{\theta}))\,\mathbf{I}\right\}, \]

with \( \boldsymbol{\varepsilon}(\boldsymbol{\theta}) = \frac{1}{2}(\nabla\boldsymbol{\theta}+(\nabla\boldsymbol{\theta})^T) \), bending modulus \( D = Et^3/(12(1-\nu^2)) \), and the shear resultants

\[ \mathbf{Q}(\boldsymbol{\theta}, w) = \kappa G t\,(\nabla w - \boldsymbol{\theta}), \]

where \( G = E/(2(1+\nu)) \) is the shear modulus and \( \kappa \) the shear correction factor.

Mixed form using potentials

The shear resultants admit the Helmholtz decomposition \( \mathbf{Q} = \nabla\varphi + \mathbf{J}\nabla\psi \) in terms of scalar potentials \( \varphi \) and \( \psi \). With the rotation matrix \( \mathbf{J} = \begin{bmatrix}0&1\\-1&0\end{bmatrix} \), the system becomes

\begin{align} \nabla\cdot\nabla\varphi &= F,\\ -\nabla\cdot(\nabla w - \boldsymbol{\theta}) &= -\nabla\cdot((\kappa Gt)^{-1}\nabla\varphi),\\ -\nabla\cdot(\mathbf{J}\boldsymbol{\theta}) &= \nabla\cdot((\kappa Gt)^{-1}\mathbf{J} \mathbf{J}\nabla\psi),\\ -\nabla\cdot(\mathbf{M}(\boldsymbol{\theta}) - \mathbf{I}\varphi - \mathbf{J}\psi) &= \mathbf{0}. \end{align}

Equations (1)-(3) are the deformation-and-potentials sub-problem according to the implemented order for the primary variables \( (\varphi, w, \psi) \). All three are scalar second-order equations in \( \varphi \), \( w \), and \( \psi \), respectively. Equation (4) is a vector second-order equation for the rotation field \( \boldsymbol{\theta} \).

Primary variables

The deformation sub-problem has three primary variables per DOF:

• shear gradient potential \( \varphi \)
• vertical deformation \( w \)
• shear curl potential \( \psi \)

The rotation sub-problem has two primary variables per DOF:

• rotation component \( \theta_x \)
• rotation component \( \theta_y \)

Classes
class	Dumux::MindlinReissnerPlateCouplingManager< MDTraits >
struct	Dumux::MindlinReissnerPlateVolumeVariablesTraits< PV, MT >
struct	Dumux::Properties::TTag::MindlinReissnerPlateDeformation
struct	Dumux::Properties::TTag::MindlinReissnerPlateRotation
struct	Dumux::Properties::ModelTraits< TypeTag, TTag::MindlinReissnerPlateDeformation >
struct	Dumux::Properties::LocalResidual< TypeTag, TTag::MindlinReissnerPlateDeformation >
struct	Dumux::Properties::VolumeVariables< TypeTag, TTag::MindlinReissnerPlateDeformation >
	Set the volume variables property. More...
struct	Dumux::Properties::ModelTraits< TypeTag, TTag::MindlinReissnerPlateRotation >
struct	Dumux::Properties::LocalResidual< TypeTag, TTag::MindlinReissnerPlateRotation >
struct	Dumux::Properties::VolumeVariables< TypeTag, TTag::MindlinReissnerPlateRotation >
	Set the volume variables property. More...

Namespaces
namespace	Dumux
namespace	Dumux::Properties
	The energy balance equation for a porous solid.
namespace	Dumux::Properties::TTag
	Type tag for numeric models.

Typedefs
using	Dumux::MindlinReissnerPlateTraits = KirchhoffLovePlateTraits
using	Dumux::MindlinReissnerPlateRotationModelTraits = KirchhoffLovePlateRotationModelTraits
template<class T >
using	Dumux::MindlinReissnerPlateDeformationVolumeVariables = KirchhoffLovePlateDeformationVolumeVariables< T >
template<class T >
using	Dumux::MindlinReissnerPlateRotationVolumeVariables = KirchhoffLovePlateRotationVolumeVariables< T >

Include dependency graph for solidmechanics/plate/mindlin_reissner/model.hh: