Inel4nodeMindlinReShellwITC
Purpose
isoparametric inelastic 4 node quad Mindlin RC layer shell element with ITC
Synopsis
ElemResp = Inel4nodeMindlinReShellwITC (action,el_no,xyz,ElemData,ElemState)
Description
INEL4NODEMINDLINRESHELLwITC isoparametric inelastic 4 node quad Mindlin RC layer shell element with ITC
ELEMRESP = INEL4NODEMINDLINRESHELLwITC (ACTION,EL_NO,XYZ,ELEMDATA,ELEMSTATE)
response of inelastic isoparametric 4 node quad Mindlin layer shell element with reinforcing layers with ITC;
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)
nrl = number of reinforcing layers (default nrl=0)
rhol = reinforcing ratio of each reinforcing layer
langle = angle of each reinforcing layer relative to x-axis
MatName{1} = function name for material model of element
MatName{2} = function name for material model of reinforcement
MatData{1} = properties for element material model
MatData{2} = properties for reinforcing 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
- shapeBhat shape function Bhat
- 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