Inel4to9nodeMindlinShell

Purpose

isoparametric 4-9 node quad layer shell element with inelastic material

Synopsis

ElemResp = Inel4to9nodeMindlinShell (action,el_no,xyz,ElemData,ElemState)

Description

INEL4to9NODEMINDLINSHELL isoparametric 4-9 node quad layer shell element with inelastic material
 ELEMRESP = INEL4to9NODEMINDLINSHELL (ACTION,EL_NO,XYZ,ELEMDATA,ELEMSTATE)
  response of isoparametric 4-9 node quadrilateral layer shell element with inelastic material;
  the element accounts for linear and nonlinear geometry for the nodal dof transformations; 
  depending on the value of the character variable ACTION the function returns information
  in data structure ELEMRESP for the element with number EL_NO, end node coordinates XYZ,
  and material and loading properties in the data structure ELEMDATA.
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  When the character variable ACTION has one of the following values,
  the function performs the listed operations and returns the results in ELEMRESP:
  ACTION = 'size': report size of element arrays
           'chec': check element property data for omissions and assign default values
           'init': initialize element history variables
           'forc': report element resisting forces
           'stif': report element stiffness matrix and resisting forces
           'mass': report lumped mass vector and consistent mass matrix
           'post': report post-processing information
           'stre': stress recovery to element nodes
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  The data structure ELEMRESP stands for the following data object(s) for each ACTION:
  ELEMRESP = ARSZ        for action = 'size' 
  ELEMRESP = ELEMDATA    for action = 'chec'
  ELEMRESP = ELEMSTATE   for action = 'init'
  ELEMRESP = ELEMSTATE   for action = 'stif'
  ELEMRESP = ELEMSTATE   for action = 'forc'
  ELEMRESP = ELEMMASS    for action = 'mass'
  ELEMRESP = ELEMPOST    for action = 'post'
  ELEMRESP = STRSREC     for action = 'stre'
  ELEMRESP is empty      for unsupported keywords
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  ARSZ   is an Boolean array of size NDF x NEN,
         where NDF = number of DOFs/node, NEN = number of nodes,
         with unit values corresponding to the active element DOFs
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  ELEMDATA is a data structure with element property information in fields
         Geom    = character variable for geometric transformation of node variables
                   (linear or corotational) (default=linear)
         t       = element thickness (default t=1)
         w       = transverse element load (default w=0)
         LdIdw   = load ID for transverse element load
         b       = element body forces (default b=[0;0;0] )
         LdIdx   = load ID for body forces b(1)
         LdIdy   = load ID for body forces b(2)
         LdIdz   = load ID for body forces b(3)
         nIP     = number of integration points for flexure (default nIP =2)
         nl      = number of layers for across depth integration (default nl=5)
         MatName = function name for material model
         MatData = properties for element material model
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  ELEMSTATE is a data structure with the current element state; it has the fields
         u     = vector of total element displacements in global reference
         Du    = vector of element displacement increments from last convergence
         DDu   = vector of element displacement increments from last iteration
         ke    = element stiffness matrix in global reference; updated under ACTION = 'stif'
         p     = element resisting force vector in global reference; updated under ACTION = 'stif' or 'forc'
         Past  = element history variables at last converged state
         Pres  = current element history variables
         lamda = row vector of current load factor(s)
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  ELEMPOST is a data structure with element response information for post-processing in fields
         Tr0         = transformation matrix from local to global coordinates
         Mat{i}.xyz  = coordinates           of membrane integration point i
         Mat{i}.eps  = strain tensor         at membrane integration point i
         Mat{i}.sig  = stress tensor         at membrane integration point i
         Mat{i}.Fxyz = coordinates           of flexural integration point i
         Mat{i}.M    = bending moment tensor at flexural integration point i
         Mat{i}.Sxyz = coordinates           of shear    integration point i
         Mat{i}.V    = shear forces          at shear    integration point i
  ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  STRSREC is a data structure with nodal stress information in field(s)
         sigNd = membrane forces
         MomNd = nodal moments
         ShrNd = nodal shear forces

Cross-Reference Information

This function calls:

• ElementLoading determines current distributed element load value
• ExtrReshu extracts displacements and increments from State and reshapes into array
• Extract_El2MatState extract material state from element state
• ADOFTran_PlateShell transforms active dof array from local to global reference for plate and shell elements
• RectPatch2Layer integration points and weights for 1d-integration of rectangular patch
• TransformMtensor2GL transformation of moment tensor from plane to global reference system
• TransformStr2GL transformation of strain and stress tensors from local to global reference system
• shape2d shape functions for 4-9 node quadrilateral element
• DefGeom_Quad determines local coordinates and corotational triad of quadrilateral element
• GeomTran_QuadPlate kinematic matrices and local displacements for a quadrilateral plate element
• kg_QuadPlate geometric stiffness matrix for quadrilateral plate element for different options
• Gauss2d Gauss integration rule in two dimensions