UNAVCO Logo
 
 
Software Help with Software Modeling 3D~Def Guide User's Manual Definitions Introduction How to drive the deformation Input file and example 1 Output More examples Pitfalls and advice Get and run the code Selected references

Output | 3D~Def Guide


  1. Specification of output
  2. Output files format
  3. Output filenames and descriptions

1. Specification of output

The bottom part of the input files describes the types of output requested. Here is an example:

bottom part of input file explanation
* Run everything else interactively comment
Non-interactive grid specification specify grid to which output is referred below
plane let this grid be a plane
vector request vector output (note: xyz format is always output)
* output coordinate system specify output coodinate system
global choose global rather than in-plane
* output files to create select whichever output is required (there are many choices, read on)
elem relative displacements on elements
displacements absolute displacements on inspection grid
rbr rigid-body rotations on inspection grid
* output file suffix 4 characters or less
t output file suffix (format allows no more than 4 characters)
* Xo,Yo,Zo (Ref. corner),strike,dip,length,width: specify inspection plane
0 55 0 90 0 59 59 see below
* number of grid points in strike and dip directions: actually the number of equal divisions along strike and dip
59 59 see below

In this example, results will be written with reference to coordinates that define a horizontal plane at the surface whose origin in global coordinates lies 55 km north of the coordinate origin, and that is 59 km along strike and dip. The plane is then divided into 59 intervals along both strike and dip, so that quantities (in the case, components of displacements and rbr) will be written at 1 km spacings. This will yield 60 results along both strike and dip, for a total of 3,600 x 3 displacement results. Displacements of all subelements will also be written to elem files.

return to top of page

2. Output files format

Output files formats can be 'xyz' (column) or 'vector' (matrix) format, and they are referenced to either the inspection grid or the elements themselves.

'xyz' format lists each component of the calculated quantity in column format next to its location in either global or in-plane coordinates.

The first three columns are the x, y, and z locations in global coordinates (it's usually better or easier to use global coordinates), the next columns are deformation quantities calculated at this location.

The equivalent vector output comprises three matrices, one for each component of displacement (or whatever quantity is being output): disp1, disp2, and disp3, where 1, 2, and 3 are cryptic mnemonics for x, y, and z. Locations are implicit in the position of the matrix elements. Vector output is useful for plotting purposes (using something like MATLAB or Spyglass) and for quick visual inspection of the results.

Remember that the vector files elem.. are written in matrix format (with (1,1) at top-left, etc.), in contrast to their order in the input file. There are some applications in which you might want to put the results of one model run back into another model. A few MATLAB commands to do this one elem block at a time are:

  • load elemds01...
  • new_elem = reshape(flipud(elemds01), rows x cols, 1)

which produces a column of, in this case, elemds01 components that can be pasted directly into the input file. A more complete matlab m-file to do this for all elem-files can be found in the m-file directory.

return to top of page

3. Output filenames and descriptions

A description of the some of the output files and the type of output is given below. Global coordinates are referenced if the quantities have x, y, or z as part of their output name (e.g., Exy, Exx, etc.) In-plane coordinates are referenced if s, d, or n (strike, dip, or normal) are part of the output file names (e.g., Esd, Ess, etc.). To avoid a cluttered look, we refer always to the global coordinates in the list below.

The underlined word or portion of a word is the name of the output file as it should be entered in the input file. For example, if you want the displacement gradient tensor, enter grad, or for rigid-body rotations, enter rbr.

Each filename, in italics below, whether output in xyz or in vector format carries a suffix (e.g., strainten.<suffix>). For clarity, this is not indicated below.

Stress tensor:

xyz: stressten (x y z sxx txy txz syy tyz szz tmax tmax's rake, tstrike, tdip)
Note: tmax is the max shear stress, tstrike and tdip are shear stresses along strike and dip of the inspection plane.
vector: Sxx, Sxy, Sxz, Syy, Syz, Szz, Syy, tmax, tmaxa, tmx, tmy

Strain tensor:

xyz: strainten (x y z Exx Exy Exz Eyy Eyz Ezz Evol)
Note: Evol is the trace of the tensor and is equal to the volume change.
vector: Exx Exy Exz Eyy Eyz Ezz Evol
Although we label these quantities as simple strains, they ARE indeed strain tensor components, e.g., Exy = (dUx/dy + dUy/dx)/2, etc.

Absolute displacements:

xyz: displacements (x y z Ux Uy Uz)
vector: disp1, disp2, disp3 (1, 2, and 3 refer to x, y, and z)

Displacement gradients:

xyz: dgten (x y z dUx/dx, dUx/dy, dUx/dz, dUy/dx, dUy/dy, dUy/dz, dUz/dx, dUz/dy, dUz/dz)
vector: dU/dx, dU/dy, etc.

Magnitudes and orientations of principal strains:

xyz: strainorient (x y z Emax, Pmax, Tmax, Eint, Pint, Tint, Emin, Pmin, Tmin
Note: For each principal strain axis, in descending order, given is the magnitude, trend, and plunge of the strain.
vector: prince_max, plunge_max, trend_max, prince_int, plunge_int, trend_int, prince_min, plunge_min, trend_min.

Invariants:

xyz: invariants (x y z Evol, failure stresses, octahedral shear stress, strain energy density)
vector: delv, aftershocks, oss, work

Rigid body rotations (rbr):

xyz: rbr (x y z rbr_x rbr_y rbr_z)
vector: rbr_deg_x, rbr_deg_y, rbr_deg_z

Optimal failure planes:

xyz: failure (x y z strike dip sv_east, sv_north, sv_up)
Note: sv_east, etc. are direction cosines of the slip vector w/r to east (global x-axis), etc.
vector: strike1, strike2, ce1, cn1, cz1, ce2, cn2, cz2
Note: ce1, etc. are the direction cosines of the slip vector for plane 1w/r to east (global x-axis), etc.

Relative sub-element displacements:

xyz: elements
vector: elemds01, elemdd01, elemdn01, etc.
Note: ds, dd, and dn refer to displacements of the hanging wall relative to the footwall in the strike, dip, and normal directions, respectively. Thus, negative values of elemds, for example, mean that the displacement is dextral or right-lateral; negative values of elemdd refers to a normal fault, etc. The numeral 01 refers to the first element specified in the input file; 02 would refer to the second, and so on.

return to top of page

 

Last modified: 2019-12-24  02:29:41  America/Denver