LE3dFrm
Purpose
3d linear frame element under linear or nonlinear geometry
Synopsis
ElemResp = LE3dFrm (action,el_no,xyz,ElemData,ElemState)
Description
LE3dFRM 3d linear frame element under linear or nonlinear geometry ELEMRESP = LE3dFRM (ACTION,EL_NO,XYZ,ELEMDATA,ELEMSTATE) response of 3d linear frame element; 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 'defo': report function handle for deformed shape '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 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 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 = FunHandle for action = 'defo' 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 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; it has the fields Geom = character variable for geometric transformation of node variables (linear, PDelta or corotational) (default=linear) E = modulus of elasticity G = shear modulus rho = mass density A = cross sectional area Iy = moment of inertia about y-axis Iz = moment of inertia about z-axis J = polar moment of inertia yornt = local y-axis orientation in global reference system (column vector) JntOff = rigid joint offsets in global X, Y, Z at element ends; column 1 for node i, column 2 for node j w = uniform element load ( w(1) = longitudinal, w(2),w(3) = transverse in y and z, resp.) e0 = initial deformations ( e(1) = axial strain, e(2),e(3) = curvature about y and z, resp.) Release = axial, torsional and end flexural releases in column vector (0=cont,1=hinge) (default=[0;0;0;0;0;0]) ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 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) ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ELEMMASS is a data structure with element mass information in fields ml = lumped mass vector mc = consistent mass matrix ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ELEMPOST is a data structure with response information for post-processing in fields v = element deformations q = element basic forces
Cross-Reference Information
This function calls:
- ExtrReshu extracts displacements and increments from State and reshapes into array
- Check3dFrmAxes check that y-axis is not co-linear with element chord
- ElmLenOr element length and x-axis orientation (direction cosines)
- GeomTran_3dFrm kinematic matrices and deformations for a 2-node 3d frame element
- kg_3dFrm geometric stiffness matrix for 2-node 3d frame element different options