Array Manager

Follow

The Array Manager in OnScale controls and contains all the data arrays used, accessed and modified during the set-up and running of a simulation. Its contents can be viewed directly through the MGR command during a simulation.

These arrays can be manipulated and exported in various formats via the DATA commands. They can also be exported to MATLAB for post-processing.

List of arrays used in OnScale:

Flex arrays:

units and sign conventions:
all normal stresses and strains are positive in tension, negative in compression.
shear strains are tensor strains (epsilon), which are half of the engineering strains (gamma).
all rotation angles are in radians unless otherwise specified.

Nodal arrays:

nodal coordinates
xcrd  i,j,k nodal indices                   original x-coordinate in non-skewed system
ycrd  i,j,k nodal indices                   original y-coordinate in non-skewed system
zcrd  i,j,k nodal indices                   original z-coordinate in non-skewed system
xcrs  i,j,k nodal indices                   updated x-coordinate in skewed system
ycrs  i,j,k nodal indices                   updated y-coordinate in skewed system
zcrs  i,j,k nodal indices                   updated z-coordinate in skewed system
xcro  i,j,k nodal indices                   original x-coordinate for large-deformation problems
ycro  i,j,k nodal indices                   original y-coordinate for large-deformation problems
zcro  i,j,k nodal indices                   original z-coordinate for large-deformation problems

(note: xcro,ycro,zcro do not exist for double-precision versions of the code; to obtain them, you need to subtract xdsp from xcrs etc.)
nodal displacements xdsp i,j,k nodal indices *calc disp* x-displacement ydsp i,j,k nodal indices *calc disp* y-displacement zdsp i,j,k nodal indices *calc disp* z-displacement rxdp i,j,k nodal indices *calc rotn* x-axis rotation (shells/beams only) rydp i,j,k nodal indices *calc rotn* y-axis rotation (shells/beams only) rzdp i,j,k nodal indices *calc rotn* z-axis rotation (shells/beams only)
nodal velocities xvel i,j,k nodal indices x-velocity yvel i,j,k nodal indices y-velocity zvel i,j,k nodal indices z-velocity rxvl i,j,k nodal indices x-axis rot'l velocity (shells/beams only) ryvl i,j,k nodal indices y-axis rot'l velocity (shells/beams only) rzvl i,j,k nodal indices z-axis rot'l velocity (shells/beams only)
nodal forces xfrc i,j,k nodal indices x-force yfrc i,j,k nodal indices y-force zfrc i,j,k nodal indices z-force rxfc i,j,k nodal indices x-axis rotational force (shells/beams only) ryfc i,j,k nodal indices y-axis rotational force (shells/beams only) rzfc i,j,k nodal indices z-axis rotational force (shells/beams only)
other nodal data velm i,j,k nodal indices *calc velm* velocity magnitude tmas i,j,k nodal indices current time step / nodal mass (zero if nodal mass is zero)

Continuum element arrays:

continuum element stresses
sgxx  i,j,k elemental indices  *calc strs*  xx-stress
sgyy  i,j,k elemental indices  *calc strs*  yy-stress
sgzz  i,j,k elemental indices  *calc strs*  zz-stress
sgxy  i,j,k elemental indices  *calc strs*  xy-stress
sgyz  i,j,k elemental indices  *calc strs*  yz-stress
sgxz  i,j,k elemental indices  *calc strs*  xz-stress

continuum element strains epxx i,j,k elemental indices *calc strn* xx-strain epyy i,j,k elemental indices *calc strn* yy-strain epzz i,j,k elemental indices *calc strn* zz-strain epxy i,j,k elemental indices *calc strn* xy-strain epyz i,j,k elemental indices *calc strn* yz-strain epxz i,j,k elemental indices *calc strn* xz-strain
continuum element maximum strains mnxx i,j,k elemental indices *calc maxe* minimum xx-strain mnyy i,j,k elemental indices *calc maxe* minimum yy-strain mnzz i,j,k elemental indices *calc maxe* minimum zz-strain mnxy i,j,k elemental indices *calc maxe* minimum xy-strain mnyz i,j,k elemental indices *calc maxe* minimum yz-strain mnxz i,j,k elemental indices *calc maxe* minimum xz-strain mxxx i,j,k elemental indices *calc maxe* maximum xx-strain mxyy i,j,k elemental indices *calc maxe* maximum yy-strain mxzz i,j,k elemental indices *calc maxe* maximum zz-strain mxxy i,j,k elemental indices *calc maxe* maximum xy-strain mxyz i,j,k elemental indices *calc maxe* maximum yz-strain mxxz i,j,k elemental indices *calc maxe* maximum xz-strain
continuum element misc data pres i,j,k elemental indices *calc pres* element pressure = (sgxx+sgyy+sgzz)/3 aprs i,j,k elemental indices *calc pres* negative of pressure (aprs is positive in compression) note: only calculated if the "aprs" option is input on the "pres" command
p2in i,j,k elemental indices *calc p2in* integral of element pressure sj2p i.j.k elemental indices *calc sj2p* square root of the second invariant of the deviatoric stress tensor = sqrt(j2') el i,j,k elemental indices cap position (units of stress) for most continuum materials note: for plih materials, used for total equivalent plastic strain
vdmg i,j,k elemental indices measure of damage, or fraction of softening, for softening continuum materials, 0 to 1 (how far beyond the peak capacity the element is: 0 = at or before peak, 1 = fully softened) note: for plih materials, used for expended fraction of post-peak plastic strain
loss i.j.k elemental indices *calc loss* accumulated energy dissipation (vdmp, sdmp, nwtn damping only; small defmn only?) volr i,j,k elemental indices *calc volr* volumetric strain rate xcrl i,j,k elemental indices *calc volr* x-component of displacement field curl vector (3d models only) ycrl i,j,k elemental indices *calc volr* y-component of displacement field curl vector (3d models only) zcrl i,j,k elemental indices *calc volr* z-component of displacement field curl vector evol i,j,k elemental indices *calc curl* volumetric strain evbr i,j,k elemental indices plastic component of volumetric strain esbr i,j,k elemental indices plastic component of deviatoric strain (von mises; similar to stress j2') eqte i,j,k elemental indices *calc eqep* equivalent total strain eqpe i,j,k elemental indices *calc eqep* equivalent plastic strain volm i,j,k elemental indices 1 / ( 12 times original element volume ) pflp needed for erosion in plot generation based on data out's sitp needed for erosion in plot generation based on data out's (changes size, so close this array before reading it in) estp i,j,k elemental indices element time step = sqrt(estp) x time step stability safety factor (typically 0.8)
material property data matn i=material# material name (in order defined in input, materials not used by elements left blank; last one is "void")
mats i=dataid,j=material# wave speeds - 1 = elastic cp, 2 = elastic cs, 3 = impedance cp, 4 = impedance cs, 5 = user-specified volumetric strain rate for artificial viscosity calculations matp i=parameter#,j=material# material parameters (varies with material type, but typically the first three are mass density, bulk modulus and shear modulus); for plih materials, 9 = yield stress sitm i,j,k elemental indices material number for each element (can be converted to material name using "matn", above) note: only exists when you are within the "site" command (special compacted form later) matt i=material# number indicating material type, for example matt(*)=20 is plih

Arrays for other element types:

bar element data
(note: bars may be reordered during the process step)
barl  i=bar#                                user-specified bar id
baro  i=bar#                                relates order in which bars were defined with the final bar element number
bart  i=bar#                                current time step = factor(typically 0.8) x sqrt(bart)
barb  i=dataid,j=bar#                       bar geometry: i: 1 = 1/length, 2 = area, 3-5 = direction cosines
                                            (needed for erosion in plot generation based on data out's)
bari  i=2,j=bar#                            global node numbers of end nodes; note that i=1 is negative if element is eroded
                                            (needed for erosion in plot generation based on data out's)

For pout & outp only:

basg  i=barid                               axial stress
baep  i=barid                               axial strain
bapl  i=barid,j=dataid                      j: 1 = equivalent plastic strain, 2 = expended fraction of available fracture energy
                                            (fracture energy calculated only when requested with the "prop" subcommand)

b3sg  i=beamid,j=int pt,k=comp              stress - k: 1 = local xx, 2 = local xy, 3 = local xz
b3ep  i=beamid,j=int pt,k=comp              strain - k: 1 = local xx, 2 = local xy, 3 = local xz
b3pl  i=beamid,j=int pt,k=dataid            for regular beam materials and beam materials based on plih continuum materials,
                                               k: 1 = equivalent plastic strain, 2 = expended fraction of available fracture energy
                                               (similar to vdmg)
                                            note: for regular beam materials, fracture energy is calculated only when requested with
                                               the "prop" subcommand)
                                            for beam materials based on sft1 continuum materials, k: 1 = cap position (el),
                                               2 = measure of damage (vdmg), 3 = plastic component of volumetric strain (evbr)

Review arrays:

time history data (* represents the local name of a time history file)
*/time i=point# time array */data i=point#,j=history#(not incl. time) data arrays */labl i=1,j=history# time history 4-character identifying label (yvel, pize, etc.) */indx i=dataid,j=history# time history identifying indices - i: 1 = history number, 2 = associated i-index, 3 = associated j-index, 4 = associated k-index
*/cord i=dataid,j=history# time history max/min and identifying coordinate information - i: 1 = associated x-coordinate, 2 = associated y-coordinate, 3 = associated z-coordinate, 4 = minimum value in history, 5 = maximum value in history
(note: to plot time history data, use "plot */data", "plot */cord", etc.)
extrapolated pressure data (after extr show command has been used) extr/peak peak amplitude (absolute value) extr/ptop peak-to-peak amplitude fft data (after freq fft command has been used; * represents the local name of an fft file) */time fft frequencies

For direct use of data arrays:

basg  i=1,j=bar#                            axial stress
baep  i=1,j=bar#                            axial strain
bapl  i=dataid,j=bar#                       i: 1 = equivalent plastic strain, 2 = expended fraction of available fracture energy
                                            (fracture energy calculated only when requested with the "prop" subcommand)

b3sg  i=comp,j=int pt,k=beam#               stress - i: 1 = local xx, 2 = local xy, 3 = local xz
b3ep  i=comp,j=int pt,k=beam#               strain - i: 1 = local xx, 2 = local xy, 3 = local xz
b3pl  i=dataid,j=int pt,k=beam#             for regular beam materials, i: 1 = equivalent plastic strain,
                                               2 = expended fraction of available fracture energy
                                            (fracture energy calculated only when requested with the "prop" subcommand)
                                            for beam materials based on sft1 continuum materials, i: 1 = cap position (el),
                                               2 = measure of damage (vdmg), 3 = plastic component of volumetric strain (evbr)

mbsg  i=comp,j=membrane#                    stress - i: 1 = local xx, 2 = local yy, 3 = local xy (3 for 3d only)
mbep  i=comp,j=membrane#                    strain - i: 1 = local xx, 2 = local yy, 3 = local xy (3 for 3d only)
mbpl  i=dataid,j=membrane#                  i: 1 = equivalent plastic strain, 2 = expended fraction of available fracture energy
                                               (similar to vdmg)
                                            (fracture energy calculated only when requested with the "prop" subcommand)
rebar stresses
rbr1  i,j,k elemental indices               rebar axial stress in local 1-direction
rbr2  i,j,k elemental indices               rebar axial stress in local 2-direction
rbr3  i,j,k elemental indices               rebar axial stress in local 3-direction
pore fluid model data
mdem  i,j,k elemental indices               fluid material model # (in order input)
pflu  i,j,k elemental indices               pore fluid pressure
vflu  i,j,k elemental indices               element volume ratio

 For field plots using data arrays:

bm3d_ep11                                   local 11-component normal strain
bm3d_ep22                                   local 22-component normal strain
bm3d_ep12                                   local 12-component shear  strain
bm3d_ep13                                   local 13-component shear  strain
bm3d_ep23                                   local 23-component shear  strain
bm3d_epls                                   equivalent plastic strain (same as array b3pl k=1)
bm3d_sg11                                   local 11-component normal stress
bm3d_sg22                                   local 22-component normal stress
bm3d_sg12                                   local 12-component shear  stress
bm3d_sg13                                   local 13-component shear  stress
bm3d_sg23                                   local 23-component shear  stress

shell element data (also see "special information commands") (note: shells may be reordered during the process step)
shsg i=comp,j=layer# stress - i for 2d: 1 = xx, 2 = yy, 3 = xz (all local) i for 3d: 1 = xx, 2 = yy, 3 = xy, 4 = xz, 5 = yz (all local) shep i=comp,j=layer# strain - i same as above sj2s i=layer# *calc sj2p* square root of the second invariant of the deviatoric stress tensor = sqrt(j2') shpl i=dataid,j=layer# for standard shell materials & plih materials, i: 1 = equivalent plastic strain, 2 = expended frac of fracture energy (use for data out) (fracture energy calculated only when requested with the "prop" subcommand) for shell materials based on sft1 continuum materials, i: 1 = el, 2 = vdmg, 3 = evbr for orthotropic shell materials, i: 1 = failure indicator for direction 1 (fiber), 2 = failure indicator for direction 2 (matrix), 3 = failure indicator for in-plane shear for yeoh type shell materials, i: 1 = first deviatoric strain invariant shlt i=shell# current time step = factor(typically 0.8) x sqrt(shlt) (note: zones may change the index) shli i=dataid,j=shell# for 2d, i: 1-2 = global node number of corner nodes, 3 = layer def number, 4 = pointer to first layer in stress and strain arrays for 3d, i: 1-4 = global node number of corner nodes, 5 = like 3 above, 6 = like 4 above i=1 negative if element eroded (needed for erosion in plot generation based on data out's) shln i=material# list of shell material model names shla i=layer_defn# list of shell layer definition names shlc i=dataid,j=layer_defn# layer data - i: 1 = number of layers, 2 = layer # of first layer shld i=layer# shell material model name for each layer (not just for each layer definition, since layers of a shell can have different material models) shlp i=dataid,j=mat# shell material properties (38 total) by material - i: 1=rho, 2=e, 3=nu, 4=sigyld, 5=ehard, 6=transgfac, 7=b=e/(1-nu**2), 8=g, etc. shll i=dataid,j=layer# shell material properties (38 total) by shell layer - i: 1=rho, 2=b, 3=g, 4=nu, 5=sigyld, 6=ehard, 7=layer thickness, etc. shlz i=dataid,j=zone# i: 1 = pointer to first shell element in zone, 2 = number of shell elements in zone shlo i=shell# order in which shells will be processed orst i=dir,j=layer# *calc orth* orthotropic shell ply direction true stresses, i=1-5 are direction components orep i=dir,j=layer# *calc orth* orthotropic shell ply direction true strains, i=1-5 are direction components 

shel_ep11                                   local 11-component normal strain
shel_ep22                                   local 22-component normal strain
shel_ep12                                   local 12-component shear  strain
shel_ep13                                   local 13-component shear  strain
shel_ep23                                   local 23-component shear  strain
shel_epls                                   equivalent plastic strain (same as array shpl i=1)
shel_sg11                                   local 11-component normal stress
shel_sg22                                   local 22-component normal stress
shel_sg12                                   local 12-component shear  stress
shel_sg13                                   local 13-component shear  stress
shel_sg23                                   local 23-component shear  stress

line element data
lnri  i=dataid,j=line#                      i: 1 = node 1, 2 = node 2, 3 = spring model, 4 = reference node 3, 5 = reference node 4
                                            (note: node 1 is set to negative after a line element is eroded)

lnrb  i=dataid,j=line#                      for type gspr, i: 1 = current distance between nodes
lnrd  i=dataid,j=line#                      i: 1 = axial force, 3 = axial relative motion (except for types eos, cap, gspr & joint);
                                               see flex manual for other i-index definitions (different for different spring types)
lnrp  i=dataid,j=spring model #             spring model properties for bolt springs - i: 2 = area, 3 = length, 4 = material number,
                                            5 = shear length, 6 = shear area
orfd  i=orifice#                            orifice air flow velocity

2d liner element data
lnrd  i=comp,j=liner#                       i: 1 = long'l stress, 2 = lateral stress, 3 = long'l strain, 4 = lateral strain

membrane element data (flex1h)

mbrd  i=comp,j=membrane#                    stress and strain - i: 1 = local xx-stress, 2 = local yy-stress, 3 = local xy-stress,
                                                                   4 = local xx-strain, 5 = local yy-strain, 6 = local xy-strain
mbnt  i=membrane#                           current time step = factor(typically 0.8) x sqrt(mbnt) (note: zones may change the index)

rigd                                        rigid substructures
rigd_xvel                                   x velocity
rigd_yvel                                   y velocity
rigd_zvel                                   z velocity
rigd_rxvel                                  x rotational velocity
rigd_ryvel                                  y rotational velocity
rigd_rzvel                                  z rotational velocity

airblast data & pressure load data
pldt  i=dataid,j=surface#                   conventional explosion or applied pressure data - i: 1 = range (ft), 2 = current pressure
                                            (model units); note that if pressure does not vary between elements, j degenerates to 1

pldw  i=dataid,j=surface#                   conwep loading coefficients - i: 2 = peak pressure (psi), 4 = time of arrival (msec),
                                            5 = positive phase duration (msec), 7 = peak negative phase pressure (psi, <0) if any,
                                            8 = negative phase duration (msec) if any
                                            notes: max i = 6 if positive phase only is requested
                                                   model time units must be msec for conwep to work properly unless the time scale
                                                      factor is used in the "plod" command to perform the conversion

pldm  i=plod defn #                         names of plods (as defined with the pdef command)
pldi  i=corner#,j=surface#                  plod surface corner indices (in global node index format), i = 1-4 corners
pldp  i=dataid,j=surface#                   load direction vector = i: 1-3 = load direction vector multipled by the element face
                                            area, 4 = element face area
pldf  i=dataid,j=plod defn #                various plod input quantities: i=2 is the time shift
pdat  i=record#,j=time point#               user-created pressure time history array for loading the model
plod_sum  i,j,k=(3,8,3)                     plod summary info - i = direction(x/y/z), k = current time, max, and time of max,
                                                 j: 1 = negative force, 2 = positive force, 3 = negative force average pressure,
                                                 4 = positive force average pressure, 5 = neg force area, 6 = pos force area,
                                                 7 = area for all force, 8 = varies (i=1 => # of faces, i=2 => # of deactivated plods) 

plod_pres                                   plod pressure
plod_range                                  plod range (for conwep)

3d interaction surface data (also see "special information commands")
it3s  i=dataid,j=inter'n elem#              elemental spring constants - i: 1 = spring stiffness, 2 = element area
it3n  i=dataid,j=inter'n node#              node and segment data - i: 1 = absolute node #, 2 = nearest surface #, 3 = nearest node #,
                                                 4 = nearest segment #, 5 = node at other segment end (for edges), 6 = surface type,
                                                 7 = previous value of nearest segment #, 8 = edge indicator (1=edge, 2=edge on
                                                 symmetry boundary, 3 = slide-on/slide-off), 9 = back node (for slide-on/slide-off)

it3d  i=dataid,j=inter'n node#              current interaction data - i: 1 = node's tributary area, 2-4 = updated x,y,z coordinates,
                                            5-7 = x,y,z velocities, 8-10 = x,y,z interaction forces, 11 = penetration,
                                            12 = xi (normalized element coordinate), 13 = eta (normalized element coordinate),
                                            14 = total relative slip, 15-17 = x,y,z friction forces, 18 = confining stress,
                                            19 = shear stress, 20-22 = original x,y,z coordinates, 23-25 = total x,y,z displacements,
                                            26-29 = node 1,2,3,4 force distribution factor, 30-32 = x,y,z components of normal vector,
                                            33-35 = x,y,z components of shear force vector, 36-38 = x,y,z components of increment of
                                            sliding displacement, 41 = distance from edge for slid_on logic, 42 = ramp factor for
                                            slid_on logic, 44-46 = x,y,z offsets (slid_on logic), 47 = dilation distance

psi3  i=dataid,j=property#                  nonlinear elastic normal spring property input (defined using the "intr jprp" command):
                                            i=1-10:  10 penetration values (last ones zero if fewer than 10 were input)
                                            i=11-20: 10 corresponding input normal stress values (ditto note above)

it3s                                        slideline spring values
it3d>xx                                     slideline data array value xx

pinball data
pnbl  i,j=part#                             pinball part links: pnbl(i,j)=1 if part i is linked to part j, zero if not linked
pnbi  i=dataid,j=defn#                      part info - i: 1 = #pinballs in part, 2 = index of first pinball in pnbn,
                                            3 = parameter model, 4 = first subsurface, 5 = #subsurfaces
pnbd  i=dataid,j=pb node#                   node data - i: 1 = radius, 2 = spring stiffness, 3-5 = x,y,z forces,
                                            6 = max penetration of any other pinball, 7-9 = x,y,z coordinates
pnbp  i=dataid,j=subsurface#                subsurface info = i: 1 = stiffness, 2 = radius , 3 = max radius
pnbc  i=pb node#                            global node number of each pinball node
pnbwf i=dataid,j=wall#                      wall data - i: 1-3 = total x,y,z forces on wall

it3s                                        slideline spring values
it3d>xx                                     slideline data array value xx

pinball data
pnbl  i,j=part#                             pinball part links: pnbl(i,j)=1 if part i is linked to part j, zero if not linked
pnbi  i=dataid,j=defn#                      part info - i: 1 = #pinballs in part, 2 = index of first pinball in pnbn,
                                            3 = parameter model, 4 = first subsurface, 5 = #subsurfaces
pnbd  i=dataid,j=pb node#                   node data - i: 1 = radius, 2 = spring stiffness, 3-5 = x,y,z forces,
                                            6 = max penetration of any other pinball, 7-9 = x,y,z coordinates
pnbp  i=dataid,j=subsurface#                subsurface info = i: 1 = stiffness, 2 = radius , 3 = max radius
pnbc  i=pb node#                            global node number of each pinball node
pnbwf i=dataid,j=wall#                      wall data - i: 1-3 = total x,y,z forces on wall

pnbl                                        pinballs
pnbl_rad                                    radii
pnbl_stiff                                  stiffness
pnbl_xfrc                                   x force
pnbl_yfrc                                   y force
pnbl_zfrc                                   z force
pnbl_pen                                    penetration
pnbl_xcrd                                   x coordinate
pnbl_ycrd                                   y coordinate
pnbl_zcrd                                   z coordinate

generated fragment data (from frag command)
frgi  i=dataid,j=fragment set#              summary fragment info: i: 1 = number of fragments, 2 = pointer into frgs array for first
                                            fragment in the set (note: typically, there will be only one fragment set)
frgs  i=dataid,j=fragment#                  detailed fragment info: i: 1 = x-coord, 2 = y-coord, 3 = z-coord,
                                            4 = x-velocity, 5 = y-velocity, 6 = z-velocity, 7 = mass

procedure data

exec:t  i="rate" proc #, j=dataid           periodically repeated procedure ("proc rate") data - j: 1 = next time procedure will be
                                            invoked (sec or msec), 2 = rate (sec or msec)

piezoelectric data

potn  i,j,k wndo nodal indices              nodal electric potential field, voltage
qfrc  i,j,k wndo nodal indices              nodal electric charge
ex    i,j,k wndo eleml indices *calc elec*  elemental electric field x-component
ey    i,j,k wndo eleml indices *calc elec*  elemental electric field y-component
ez    i,j,k wndo eleml indices *calc elec*  elemental electric field z-component (3d models only)
pize  i=dataid,j=electrode#                 electrode data - i: 1 = potential voltage, 2 = total charge, 3 = current (d(charge)/dt)

piez                                        electrodes
piez_voltage                                electrode voltage
piez_charge                                 electrode charge
piez_current                                electrode current

circuit data
crtv  i=circuit node #                      voltage at a circuit node (two circuit nodes per circuit element, plus electrode and source)
crtq  i=circuit node #                      charge at a circuit node (two circuit nodes per circuit element, plus electrode and source)

magnetostrictive materials
mage  i=dataid, j=coil#                     coil info - i: 1 = voltage across coil, 2 = total charge, 3 = current
bx    i,j,k local elemental indices         x-component of element magnetic field (must issue "magn calc" command)
by    i,j,k local elemental indices         y-component of element magnetic field (must issue "magn calc" command)
bz    i,j,k local elemental indices         z-component of element magnetic field (must issue "magn calc" command)

special behavior model data (associated with "modl" command group)
rbr1  i,j,k elemental indices               axial rebar stress in 1-direction
rbr2  i,j,k elemental indices               axial rebar stress in 2-direction
rbr3  i,j,k elemental indices               axial rebar stress in 3-direction
dgap  i,j,k elemental indices               gap strains caused by debonding

static analysis data
stat  i=dataid                              conjugate gradient solution variables - i: 1 = |p||k|{p}, 2 = lambda, 3 = [du][q]{du},
                                            4 = alpha

dynamic relaxation data

drlx  i=dataid                              time-varying dynamic relaxation data: 1 = critically damped angular frequency (for "auto"
                                            option only), 7 = mass-weighted velocity norm, 8 = prior peak value of drlx(7),
                                            9 = moving average of drlx(7), 10 = drlx(9)/drlx(8) (velocity norm),
                                            11 = time of peak drlx(7)

pulse/echo extrapolation data
echo  (scalar)                              extrapolated pressure at the point of interest (current value)
echp  i=time step #                         full time history of "echo"

transformed coordinate system quantities
displacements, velocities, stresses, strains: see flex manual, "calc tfrm"

other

timstp (scalar)                             current model timestep
func  (scalar)                              user-defined function
vole  i=volume#                *calc volm*  volume of a group of continuum elements
clvn  i=dataid,j=volume#   *calc volmintg*  element array volume integral input data - i: 1 = user-input name, 2 = element array
clvr  i=volume#            *calc volmintg*  element array volume integral
fragd i=dataid,j=frag set #    *calc frag*  fragment data - i: 1 = # frag nodes, 2 = total mass, 3 = x-momentum, 4 = y-momentum,
                                            5 = z-momentum, 6 = lowest velocity, 7 = highest velocity
enrd  i=comp,j=energy region   *calc enrg*  energy quantities - i if option = totl (default): 1 = x-momentum, 2 = y-momentum,
                                            3 = z-momentum, 4 = kinetic energy, 5-8 = peak values of previous quantities,
                                            9 = total translational mass
                                            i if option = splt: 1 = + x-momentum, 2 = + y-momentum, 3 = + z-momentum,
                                            4 = + kinetic energy, 5 = - x-momentum, 6 = - y-momentum, 7 = - z-momentum,
                                            8 = - kinetic energy, 9 = total translational mass
enrname i=energy region    *calc enrgintg*  name of energy integral region
enrint  i=energy region    *calc enrgintg*  internal energy per unit volume (optionally, volume integral over specified window)
masv  i=added mass #                        added mass (can be used to do nodal field plots of added mass magnitude)
matr  i=1-6,j=1-2                           summary of element types in model & similar; to plot quantities as a function of
                                            coordinate in a model without continuum elements, issue "set matr on 2 2 1 1"

Miscellaneous arrays: 

data array averages
avgn i=dataid,j=request#       *calc avrg*  array average input data - i: 1 = user-assigned name, 2 = array to be averaged
avrg i=dataid,j=request#       *calc avrg*  array average output data - i: 1 = weighted average, 2 = weighted sum, 3 = sum of weights

cutting-plane resultant forces and moments resfrc i=comp,j=plane# *calc resultant* resultant forces on cutting plane - i: 1 = local x, 2 = local y, 3 = local z resmom i=comp,j=plane# *calc resultant* resultant moments on cutting plane - i: 1 = local x, 2 = local y, 3 = local z resnpt i=comp,j=plane# *calc resultant* global coordinates of neutral axis intersection with cutting plane - i: 1-3 = global xyz for in-plane moments, 4-6 = global xyz for out-of-plane moment erosion ("mods") data eroded i=dataid erosion statistics: 1 = total # eroded elements, 2 = # eroded continuum elements, 3 = # eroded shells, 4 = # eroded beams, 5 = # eroded bars, 6 = # eroded membranes, 7 = total # bad (zero timestep) elements, 8 = # bad continuum elements, 9 = # bad shells, 10 = # bad beams, 11 = # bad bars, 12 = # bad membranes
rigid substructure data rgfr i=comp,j=substructure# forces on rigid substructure - i: 1 = x-force, 2 = y-force, 3 = z-force, 4 = x-rot'l force, 5 = y-rot'l force, 6 = z-rot'l force rgvl i=comp,j=substructure# velocities of rigid substructure - i: 1 = x-vel, 2 = y-vel, 3 = z-vel, 4 = x-rot'l vel, 5 = y-rot'l vel, 6 = z-rot'l vel rgds i=comp,j=substructure# displacements of rigid substructure - i: 1 = x-disp, 2 = y-disp, 3 = z-disp, 4 = x-rot'l disp, 5 = y-rot'l disp, 6 = z-rot'l disp (note: components 4 through 6 are based on the integrated velocities, and thus may not be accurate for large rotations) rgnm i=substructure# list of rigid substructure names rgcg i=dataid,j=substructure# rigid substructure center of gravity data - i: 1 = x-cg, 2 = y-cg, 3 = z=cg, 4 = 1.0 if user-specified cg or 0.0 if internally calculated rgms i=dataid,j=substructure# rigid substructure mass data - i: 1-3 = current time step / total x,y,z translation mass, 4-9 = current time step / total moment of inertia components, 10-12 = user-defined additional mass rgdt i=dataid,j=substructure# rigid substructure index data - i: 1 = number of nodes in substructure, 2 = pointer into rigid node arrays for this substructure, 3 = pointer into rigid node radius arrays rgfx i=dof,j=substructure# rigid substructure fixities: i = 1-6 dof's, set to 1 if fixed or 0 if not rgdf i=dof,j=substructure# rigid substructure releases: i = 1-6 dof's, set to 1 if released (free) or 0 if not rgsf i=dof,j=substructure# rigid substructure scale factor array for boundary conditions: i = 1-6 dof's rgbc i=dataid,j=substructure# rigid substructure boundary condition array - i: 1 = type (forc or vel), 2 = definition option ('scale' if same time history applied to all dofs, 'thist' if different histories), 3-8 = time history name for each dof (note: for 'scale', the same time history should be applied to all dofs) rgsl i=index,j-substructure# rigid substructure slave index array, i,j,k of master node if any rgnd i=rigid_node# nodal array (nodal location for each rigid node) 

Special information commands:

shell element information commands
shel
     info indx {n}                          prints corner nodes and associated layers for shell element n
     info layr {l}                          prints shell element number and corner nodes for layer l
     info ijk  {i} {j} {k}                  prints all shell elements connected to node i j k
3d interaction surface information commands (after "prcs" step only)
intr
     info pent                              prints i,j,k,penetration, and slideline table index for node with largest penetration
     info indx {i} {j} {k}                  prints global node index and slideline table index of node i j k
     info ijk  {n}                          prints i,j,k, and global node index of slideline table node n
0 out of 0 found this helpful

Comments

0 comments

Please sign in to leave a comment.