This section explains the content of the reduced displacement file (jobname.rdsp).
See The Standard Header for Mechanical APDL Binary Files for a description of this set. File number (Item 1) is 10.
*comdeck,fdrdsp
c *** Copyright ANSYS. All Rights Reserved.
c *** ansys, inc
c ********** description of reduced displacement file **********
character*8 RDSPNM
parameter (RDSPNM='rdsp ')
integer RDSPHDLEN
parameter (RDSPHDLEN=80)
LONGINT rdspfpL, rdspfp
integer rdspbk, rdsput
common /fdrdsp/ rdspfpL, rdspbk, rdsput
equivalence (rdspfp,rdspfpL)
c write: lnfrcl,lnfrin,lnfrwr
c read: rdtrs
c ********** common variable descriptions ***********
co rdspfpL file position on file rdsp
co rdspbk block number for file rdsp
co rdsput file unit for file rdsp
c See fddesc for documentation of how binary files are stored.
c
c ********** file format **********
c recid tells the identifier for this record. Not all records will have
c identifiers -- they are only indicated for those records whose
c record pointers are stored in the second file header.
c type tells what kind of information is stored in this record:
c i - integer
c dp - double precision
c cmp - complex
c nrec tells how many records of this description are found here
c lrec tells how long the records are (how many items are stored)
c recid type nrec lrec contents
c --- i 1 100 standard ANSYS file header (see binhed8 for
c details of header contents)
c --- i 1 80 .RDSP FILE HEADER
c
c fun10, nmrow, nmatrx, nmode, numdof,
c maxn, wfmax, lenbac, ngaps, ncumit, (10)
c kan, nres, ndva, nvect, DSPfmt,
c minmod, 0, modlstp, ndefdval, nEnfDof, (20)
c ptrDOF, ptrDAMP, ptrDAMPh, 0, 0,
c ptrFRQ, ptrDSP, 0, 0, 0, (30)
c ptrFRQh, ptrDSPh, ptrDVA, ptrDVAh,nrkeyPert,
c kPerturb, keyVA,Glblenbac, 0, 0, (40)
c 0, 0, 0, 0, 0,
c 0, 0, 0, 0, 0, (50)
c 0, 0, 0, 0, 0,
c 0, 0, 0, 0, 0, (60)
c 0, 0, 0, 0, 0,
c 0, 0, 0, 0, 0, (70)
c 0, 0, 0, 0, 0,
c 0, 0, 0, 0, 0 (80)
c each item in header is described below:
c fun10 - unit number (rdsp file is 10)
c nmrow - number of rows/columns in matrices
c nmatrx - number of reduced matrices on the
c file
c nmode - number of modes extracted during
c modal analysis
c numdof - number of dofs per node
c maxn - maximum node number
c wfmax - maximum wavefront
c lenbac - number of nodes
c ngaps - number of gaps
c ncumit - total number of iterations done
c during analysis
c kan - analysis type
c = 5 for MSUP transient analysis
c nres - number of residual vectors used
c modlstp- multiple load step key
c ndva - length of DVA (for restart)
c nvect - number of available load vectors
c DSPfmt - 0,physical disps .ne.0,modal coords
c minmod - smallest mode number used
c ndefdval- number of defined enforced motion
c maxEnf - maximum enforced values
c ptrDOF - pointer to degree of freedom set
c ptrDAMP - pointer to damping values
c ptrDAMPh - High part of DAMP pointer
c ptrDNC - pointer to nodal constraints
c ptrFRQ - pointer to the frequencies
c ptrDSP - pointer to the step solution data:
c - calculated displacements
c - load vector scale factors
c - gap restoring forces
c - calculated velocities
c - calculated accelerations
c ptrFRQh- High part of frequency ptr
c ptrDSPh- High part of displacement ptr
c ptrDVA - pointer to modal disp, velo and acc (for restart)
c ptrDVAh- High part of modal disp, velo and acc (for restart)
c nrkeyPert - nrkey setting of base analysis (Linear Perturbation)
c kPerturb - Linear Perturbation key
c keyVA - Key for velocities/accelerations on file
c Glblenbac - global number of nodes (== 0 unless using Distributed Ansys)
c 0 - position not used
c --- i 1 numdof Degrees of freedom per node
c (curdof(i),i=1,numdof)
c dof reference numbers are:
c UX = 1, UY = 2, UZ = 3, ROTX= 4, ROTY= 5, ROTZ= 6, AX = 7, AY = 8
c AZ = 9, VX =10, VY =11, VZ =12, GFV1=13, GFV2=14, GFV3=15, WARP=16
c CONC=17, HDSP=18, PRES=19, TEMP=20, VOLT=21, MAG =22, ENKE=23, ENDS=24
c EMF =25, CURR=26, SP01=27, SP02=28, SP03=29, SP04=30, SP05=31, SP06=32
c --- i 1 lenbac Nodal equivalence table. This table equates
c the number used for storage to the actual
c node number
c (Back(i),i=1,lenbac)
c --- i 1 Glblenbac Global nodal equivalence table. This
c table equates the number used for storage
c to the actual node number. Only written
c by the master process in Distributed Ansys
c (GlbBack(i),i=1,Glblenbac)
c --- dp 1 10 Time information:
c dtime, 0.0, 0.0, 0.0, 0.0,
c 0.0, 0.0, 0.0, 0.0, timend
c each item is described below:
c
c dtime - the time increment
c timend - the final time of the analysis
c 0.0 - position not used
c DOF i 1 nmrow Degree of freedom set used
c The DOFs are calculated as (N-1)*numdof+DOF,
c where N is the position number of the node in
c the nodal equivalence table and DOF is the
c DOF reference number given above.
c --- i 1 nmrow+1 Original reduced set of DOFs used.
c The DOFs are calculated as (N-1)*numdof+DOF,
c where N is the position number of the node in
c the nodal equivalence table and DOF is the
c DOF reference number given above.
c DAMP dp 1 nmode+10 Damping values.
c This record is present only for analysis using
c the mode superposition method with MCkey
c activated - TRNOPT cmd.
c There are nmode+10 entries:
c 1 : nmode - effective damping ratios
c nmode+1 : nmode+5 - alphad, betad, dmprat, dmpstr, dsmpFreq
c nmode+6 : nmode+10 - additional entries - zeroed out for now
c FRQ dp 1 nmrow/ Frequencies extracted from the modal analysis.
c ndva/ This record is present only for analyses using
c ndva-nEnfDof the mode superposition method.
c If DSPfmt is 0, the first nmrow values are the
c frequencies extracted from the modal analysis.
c The remaining values have no meaning. Otherwise,
c there are ndva frequencies if there are no
c enforced motions and ndva less the number of enforcced
c motions otherwise.
c If DSPfmt=0 :(freq(i),i=1,nmrow)
c If DSPfmt!=0 and ndefdval=0:
c (freq(i),i=1,ndva)
c If DSPfmt!=0 and ndefdval>0:
c (freq(i),i=1,ndva-nEnfDof)
c *** The next 6 to 7 records are repeated (as a group) until the time value
c *** equals the value of timend. The number of iterations is stored as
c *** ncumit. (see above records that deal with time)
c DSP dp 1 nmrow+7/ Calculated displacements
c ndva+7 The first nmrow entries are the displacements
c in the same order as the original set of DOFs
c (see record AFTER ptrDOF). If DSPfmt=0, these
c are physical displacements, If DSPfmt!=0,
c these are the ndva modal coordinates instead
c of the nmrow entries.
c
c For the last six entries:
c 1. Time for these displacements
c 2. Load step number
c 3. Substep number
c 4. Cumulative iteration number
c 5. Scale factor
c 6. numdeflvs - number of scale factors
c 7. kwrval - key to control writing velocities
c and accelerations at each substep (write if keyVA = 1 and kwrval = 1)
c
c --- i 1 Note: If, upon reading of this record, there
c is less than nmrow+5 items in the record,
c then only a selected set of nodes were
c output. Another record follows (integer, less
c than lenbac long) which contains the list of
c nodes for which DOF solutions are available.
c --- i 1 numdeflvs lvscal table scale factor IDs
c (ilvscID(i),i=1,numdeflvs)
c --- dp 1 numdeflvs lvscal table scale factor values
c (dlvscVal(i),i=1,numdeflvs)
c --- dp 1 ngaps Gap restoring forces. This record is
c present only if ngaps > 0.
c (fgaps(i),i=1,ngaps)
c --- dp 1 nmrow/ndva Calculated velocities (present if keyVA = 1).
c nmrow entries if DSPfmt=0, ndva entries if DSPfmt!=0
c --- dp 1 nmrow/ndva Calculated accelerations (present if keyVA = 1)
c nmrow entries if DSPfmt=0, ndva entries if DSPfmt!=0
c *** The next 3 records are kept for possible restart in mode superposition
c *** transient. They are overwritten upon restarting. They are written once (last
c *** loadstep).
c DVA dp 1 ndva+6 Calculated modal displacements
c The first ndva entries are the modal
c displacements. For the last six entries:
c 1. Time for these displacements
c 2. Load step number
c 3. Substep number
c 4. Cumulative iteration number
c 5. Scale factor
c 6. numdeflvs - number of scale factors
c --- dp 1 ndva Calculated modal velocities
c --- dp 1 ndva Calculated modal accelerations