chrono::fea::ChElasticityReissnerOrthotropic Class Reference

Description

Elasticity of 6-field Reissner-Mindlin shells (kinematically-exact shell theory as in Witkowski et al.) to be used in a ChMaterialShellReissner.

This class implements material properties for a layer from the Reissner theory, for the case of orthotropic linear elastic material. This is useful for laminated shells. One direction can be made softer than the other. Note that the angle and the thickness are defined when adding a material with this elasticity to a shell finite element (ex. ChElementShellReissner4) as a layer. Previously: ChMaterialShellReissnerOrthotropic

#include <ChMaterialShellReissner.h>

Inheritance diagram for chrono::fea::ChElasticityReissnerOrthotropic:
Collaboration diagram for chrono::fea::ChElasticityReissnerOrthotropic:

Public Member Functions

 ChElasticityReissnerOrthotropic (double m_E_x, double m_E_y, double m_nu_xy, double m_G_xy, double m_G_xz, double m_G_yz, double m_alpha=1.0, double m_beta=0.1)
 Construct an orthotropic material. More...
 
 ChElasticityReissnerOrthotropic (double m_E, double m_nu, double m_alpha=1.0, double m_beta=0.1)
 Construct an orthotropic material as sub case isotropic. More...
 
double Get_E_x () const
 Return the elasticity moduli, on x.
 
double Get_E_y () const
 Return the elasticity moduli, on y.
 
double Get_nu_xy () const
 Return the Poisson ratio, for xy.
 
double Get_nu_yx () const
 Return the Poisson ratio, for yx (follows xy as it must be nu_yx*E_x = nu_xy*E_y)
 
double Get_G_xy () const
 Return the shear mod, in plane.
 
double Get_G_xz () const
 Return the shear mod, transverse.
 
double Get_G_yz () const
 Return the shear mod, transverse.
 
double Get_alpha () const
 Return the shear factor.
 
double Get_beta () const
 Return the torque factor.
 
virtual void ComputeStress (ChVector<> &n_u, ChVector<> &n_v, ChVector<> &m_u, ChVector<> &m_v, const ChVector<> &eps_u, const ChVector<> &eps_v, const ChVector<> &kur_u, const ChVector<> &kur_v, const double z_inf, const double z_sup, const double angle)
 The FE code will evaluate this function to compute per-unit-length forces/torques given the u,v strains/curvatures. More...
 
virtual void ComputeStiffnessMatrix (ChMatrixRef mC, const ChVector<> &eps_u, const ChVector<> &eps_v, const ChVector<> &kur_u, const ChVector<> &kur_v, const double z_inf, const double z_sup, const double angle)
 /// Compute the 12x12 stiffness matrix [Km] , that is [ds/de], the tangent of the constitutive relation stresses/strains. More...
 

Additional Inherited Members

- Public Attributes inherited from chrono::fea::ChElasticityReissner
ChMaterialShellReissnersection
 

Constructor & Destructor Documentation

◆ ChElasticityReissnerOrthotropic() [1/2]

chrono::fea::ChElasticityReissnerOrthotropic::ChElasticityReissnerOrthotropic ( double  m_E_x,
double  m_E_y,
double  m_nu_xy,
double  m_G_xy,
double  m_G_xz,
double  m_G_yz,
double  m_alpha = 1.0,
double  m_beta = 0.1 
)

Construct an orthotropic material.

Parameters
m_E_xYoung's modulus on x
m_E_yYoung's modulus on y
m_nu_xyPoisson ratio xy (for yx it holds: nu_yx*E_x = nu_xy*E_y)
m_G_xyShear modulus, in plane
m_G_xzShear modulus, transverse
m_G_yzShear modulus, transverse
m_alphashear factor
m_betatorque factor

◆ ChElasticityReissnerOrthotropic() [2/2]

chrono::fea::ChElasticityReissnerOrthotropic::ChElasticityReissnerOrthotropic ( double  m_E,
double  m_nu,
double  m_alpha = 1.0,
double  m_beta = 0.1 
)

Construct an orthotropic material as sub case isotropic.

Parameters
m_EYoung's modulus on x
m_nuPoisson ratio
m_alphashear factor
m_betatorque factor

Member Function Documentation

◆ ComputeStiffnessMatrix()

void chrono::fea::ChElasticityReissnerOrthotropic::ComputeStiffnessMatrix ( ChMatrixRef  mC,
const ChVector<> &  eps_u,
const ChVector<> &  eps_v,
const ChVector<> &  kur_u,
const ChVector<> &  kur_v,
const double  z_inf,
const double  z_sup,
const double  angle 
)
virtual

/// Compute the 12x12 stiffness matrix [Km] , that is [ds/de], the tangent of the constitutive relation stresses/strains.


Parameters
mCtangent matrix
eps_ustrains along u direction
eps_vstrains along v direction
kur_ucurvature along u direction
kur_vcurvature along v direction
z_inflayer lower z value (along thickness coord)
z_suplayer upper z value (along thickness coord)
anglelayer angle respect to x (if needed) -not used in this, isotropic

Reimplemented from chrono::fea::ChElasticityReissner.

◆ ComputeStress()

void chrono::fea::ChElasticityReissnerOrthotropic::ComputeStress ( ChVector<> &  n_u,
ChVector<> &  n_v,
ChVector<> &  m_u,
ChVector<> &  m_v,
const ChVector<> &  eps_u,
const ChVector<> &  eps_v,
const ChVector<> &  kur_u,
const ChVector<> &  kur_v,
const double  z_inf,
const double  z_sup,
const double  angle 
)
virtual

The FE code will evaluate this function to compute per-unit-length forces/torques given the u,v strains/curvatures.

Parameters
n_uforces along u direction (per unit length)
n_vforces along v direction (per unit length)
m_utorques along u direction (per unit length)
m_vtorques along v direction (per unit length)
eps_ustrains along u direction
eps_vstrains along v direction
kur_ucurvature along u direction
kur_vcurvature along v direction
z_inflayer lower z value (along thickness coord)
z_suplayer upper z value (along thickness coord)
anglelayer angle respect to x (if needed) -not used in this, isotropic

Implements chrono::fea::ChElasticityReissner.


The documentation for this class was generated from the following files:
  • /builds/uwsbel/chrono/src/chrono/fea/ChMaterialShellReissner.h
  • /builds/uwsbel/chrono/src/chrono/fea/ChMaterialShellReissner.cpp