Released Dmaps

8/14/2019 Released Dmaps

1/395

NX Nastran 5 DMAP Updates and

Additions


2/395

Proprietary & Restricted Rights Notice

2007 UGS Corp. All Rights Reserved. This software and related documentation are proprietaryto UGS Corp.

NASTRAN is a registered trademark of the National Aeronautics and Space Administration. NXNastran is an enhanced proprietary version developed and maintained by UGS Corp.

MSC is a registered trademark of MSC.Software Corporation. MSC.Nastran and MSC.Patranare trademarks of MSC.Software Corporation.

All other trademarks are the property of their respective owners.

2 NX Nastran 5 DMAP Updates and Additions


3/395

Contents

Updated and New DMAP Modules . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-1

Updated Modules . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-1New Modules . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-92

Updated and New Datablocks . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-1

Updated Datablocks . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-1New Datablocks . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-283

Updated and New Subdmaps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-1

NX Nastran 5 DMAP Updates and Additions 3


4/395


5/395

Chapter

1 Updated and New DMAPModules

1.1 Updated Modules

BCDR

Drives a DMAP loop based on the boundary condition Case Control commands SPC and MPC.

Format:

BCDR CASECC//SEID/SOLAPP/S,N,NSKIP/S,N,NLOADS/S,N,BCFLAG/S,N,SPC/S,N,MPC/S,N,SUPORT/S,N,BCSET/S,N,BGSET/S,N,BOLTPRE/S,N,RIGID/S,N,LOAD/S,N,LSEQ/S,N,STATSUB/S,N,BC/BCLABL $

Input Data Block:

CASECC Table of Case Control command images. Output by IFP1.

Output Data Block:

None.

Parameters:

SEID Input-integer-no default. Superelement identification number.

SOLAPP Input-character-no default. Design optimization analysis type. Currentlynot used.

NSKIP Input/output-integer-no default. The record number in CASECC correspondingto the first subcase of the current boundary condition.

NLOADS Output-integer-default=0. The number of subcase records contiguous withrespect to the MPC and SPC command in the first subcase of the currentboundary condition.

BCFLAG Output-logical-no default. Set to FALSE at the last boundary condition.

NX Nastran 5 DMAP Updates and Additions 1-1


6/395

Chapter 1 Updated and New DMAP Modules

SPC Output-integer-default=0. SPC Case Control command set identificationnumber specified in the third word of the SKIP-th record of CASECC

MPC Output-integer-default=0. MPC Case Control command set identification

number specifi

ed in the second word of the NSKIP-th record of CASECC.

SUPORT Output-integer-default=0. SUPORT Case Control command set identificationnumber specified in the 255-th word of the NSKIP-th record of CASECC.

BCSET Output-integer-default=0. Contact set identification number.

BGSET Output-integer-default=0. Glue set identification number.

BOLTPRE Output-integer-default=0. Bolt preload identification number.

RIGID Output-character-default=linear. Selects the rigid element processing methodfor RBAR, RBE1, RBE2, RROD and RTRPLT elements.

LOAD Output-integer-default=0. LOAD Case Control command set identificationnumber specified in the fourth word of the NSKIP-th record of CASECC.

LSEQ Output-integer-default=0. LOADSET Case Control command set identificationnumber specified in the 205-th word of the NSKIP-th record of CASECC.

STATSUB Output-integer-default=0. STATSUB Case Control command set identificationnumber specified in the 256-th word of the NSKIP-th record of CASECC.

BC Output-integer-default=0. BC Case Control command set identificationnumber specified in the 257-th word of the NSKIP-th record of CASECC.

BCLBL Input-integer-default=0. f06fi

le page header control.-1 Clear page header

0 Initialize page header without page eject

1 Initialize page header with page eject.

Example:

Excerpt from subDMAP PHASE0:

BCFLAG=TRUE $

NSKIP=0 $DO WHILE ( BCFLAG ) $

BCDR CASES//SEID//S,N,NSKIP/S,N,NLOADS/S,N,BCFLAG/S,N,SPC/S,N,MPC/S,N,SUPORT/S,N,BCSET/S,N,BGSET/S,N,BOLTPRE/S,N,RIGID/S,N,LOAD/S,N,LSEQ//S,N,BC $

.

.

.ENDDO $

1-2 NX Nastran 5 DMAP Updates and Additions


7/395

Updated and New DMAP Modules

DDRMM

Performs matrix method data recovery

Computes data recovery items (stress, displacements, forces, strains, forces) directly from the

modal solution in frequency response, transient response, or scaled response spectra analysisusing the matrix method.

Format:

DDRMM CASECC,UH,OL,IUG,IQG,IES,IEF,XYCDB,IUG1,IQG1,IES1,IEF1/OUG,OQG,OES,OEF,UNUSED5/OPTION/NOCOMP/PEXIST/ACOUSTIC/ACOUT/PREFDB/SEID $

Input Data Blocks:

CASECC Table of Case Control command images

EST Element summary table

UH Solution matrix for the h-set (modal degrees-of-freedom). Modal displacementsonly in frequency response. Modal displacements, velocities, and accelerationsin transient response.

OL Transient response time output list or frequency response frequency output list

IUG Table of displacements due to unit modal displacement in SORT1 or SORT2format

IQG Table of single point forces of constraint due to unit modal displacement inSORT1 or SORT2 format

IES Table of element stresses or strains due to unit modal displacement in SORT1or SORT2 format. For strains, NOCOMP must be set to 3.

IEF Table of element forces due to unit modal displacement in SORT1 or SORT2format.

XYCDB Table of x-y plotting commands

IUG1 Table of displacements due to unit modal displacement in SORT1. Requiredonly for the sparse data recovery in SORT2 format.

IQG1 Table of single point forces of constraint due to unit modal displacement inSORT1 format. Required only for the sparse data recovery in SORT2 format.

IES1 Table of element stresses or strains due to unit modal displacement in SORT1format. Required only for the sparse data recovery in SORT2 format.

IEF1 Table of element forces due to unit modal displacement in SORT1 format.Required only for the sparse data recovery in SORT2 format.



8/395


Output Data Blocks:

OUG Table of displacements in SORT1 or SORT2 format

OQG Table of single point forces of constraint SORT1 or SORT2 format

OES Table of element stresses or strains in SORT1 or SORT2 format

OEF Table of element forces in SORT1 or SORT2 format

UNUSED5 Unused and can be purged

Parameters:

OPTION Input-character-default=ABS Response summation method for scaledresponse spectra analysis only. Possible values are:

ABS Absolute

SRSS Square root of the sum of the squares

NRL Naval Research Laboratory (new)

NRLO Naval Research Laboratory (old)

NOCOMP Input-integer-default=0. Set to 3 if IES is element strains.

ACOUSTIC Input-integer-default=0. Fluid-structure analysisflag.

PEXIST Input-logical-default=FALSE. Set to TRUE if p-elements are present.

0 No fluid elements exist.

1 Penalty or fluid acoustic elements exists.

2 Fluid/structure coupling exists.

ACOUT Input-character-default=PEAK Type of acoustic pressure output influid-structural analysis.

RMS Root-mean-square

PEAK Peak

PREFDB Input-real-default=1.0. Peak pressure reference for pressure level in unitsof dB or dBA.

SEID Input-integer-default=0. Superelement identification number.

Remarks:

1. If UD is a real matrix and OL is purged, a scaled response spectra analysis is assumed.



9/395


2. OUG, OQG, OES, and OEF are suitable for printing or punching by the OFP module.

3. SDR2 is used to compute IUG, IQG, IES, and IEF, which are results due to a unit modaldisplacement (eigenvector).

EFFMASS

Computes modal effective mass

Compute the modal effective mass based on the normal modes.

Format:

EFFMASS CASECC,MAA,PHA,LAMA,USET,BGPDT,UNUSED,CSTM,VGQ/TEMF,EMM,DMA,MEMF,MPFEM,MEM,MEW/SEID/WTMASS/S,N,CARDNO/SETNAM/IUNIT/EFOPT $

Input Data Blocks:


MAA Mass matrix in a-set or g-set

PHA Normal modes eigenvector matrix in the a-set or g-set

LAMA Normal modes eigenvalue summary table

USET Degree-of-freedom set membership table for g-set

BGPDT Basic grid point definition table

UNUSED Unused

CSTM Table of coordinate system transformation matrices

VGQ Partitioning vector which is g-set size and contains values of 1.0 at rowscorresponding to degrees-of-freedom in the q-set

Output Data Blocks:

TEMF Total effective mass fraction table

EMM Effective mass matrix

MA Rigid body mass matrix for the a-set

MEMF Modal effective mass fraction table

MPFEM Modal participation factors for effective mass

MEM Modal effective mass matrix



10/395


MEW Modal effective weight matrix

Parameters:

SEID Input-integer-no default. Superelement identification number.

WTMASS Input-real-no default. Scale factor on structural mass matrix. See the NXNastran Quick Reference Guide.

CARDNO Input/output-integer-default=0. Punch file line counter. CARDNO isincremented by one for each line written to the punch file and is also writteninto columns 73-80 of each line.

SETNAM Input-character-default=g Degree-of-freedom set name.

IUNIT Input-integer-no default. IUNIT is the Fortran unit number on which the data

blocks are to be written if the PLOT option is requested.

EFOPT Input-integer-default=0. When set to 1, ignores case control parameters andforces output of all results.

Remarks:

None.

EMG

Computes elemental matricesComputes elemental matrices for stiffness, differential stiffness, mass, damping, heat conduction,or heat capacity.

Format:

EMG EST,CSTM,MPT,DIT,UNUSED5,UG,ETT,EDT,DEQATN,DEQIND,BGPDT,GPSNT,ECTA,EPTA,EHTA,DITID,EBOLT/KELM,KDICT,MELM,MDICT,BELM,BDICT/S,N,NOKGG/S,N,NOMGG/S,N,NOBGG/S,N,NOK4GG/S,N,NONLHT/COUPMASS/TEMPSID/DEFRMSID/PENFAC/NOPNLT/LUMPD/LUMPM/MATCPX/KDGEN/TABS/SIGMA/K6ROT/LANGLE/NOBKGG/ALTSHAPE/PEXIST/FREQTYP/FREQVAL/FREQWA/UNSYMF/S,N,BADMESH/S,N,UNUSED $

Input Data Blocks:

EST Element summary table


MPT Table of Bulk Data entry images related to material properties



11/395


DIT Table of TABLEij Bulk Data entry images

UNUSED5 Unused and can be purged

UG Displacement matrix in g-set. Required only for differential stiffnessgeneration.

ETT Element temperature table

EDT Table of Bulk Data entry images related to element deformation. Requiredonly for differential stiffness generation.

DEQATN Table of DEQATN Bulk Data entry images

DEQIND Index table to DEQATN data block


GPSNT Grid point shell normal table

ECTA Secondary element connectivity table

EPTA Secondary table of Bulk Data entry images related to element properties

EHTA Secondary element hierarchical table

DITID Table of identification numbers in DIT

EBOLT Table of element ids of bolt elements

Output Data Blocks:

KELM Table of element matrices for stiffness, heat conduction, differential, orfollower stiffness

KDICT KELM dictionary table

MELM Table of element mass matrices

MDICT MELM dictionary table

BELM Table of element damping or heat capacity matrices

BDICT BELM dictionary table

Parameters:

NOKGG Input/output-integer-default=-1. KELM and KDICT generation flag.

Input

Generate



12/395


-1 Do not generate

Output

>Generated

-1 Not generated

NOMGG Input/output-integer-default=-1. Same as NOKGG except for MELM andMDICT.

NOBGG Input/output-integer-default=-1. Same as NOKGG except for BELM andBDICT.

NOK4GG Input/output-integer-default=1. Structural damping generation flag. Set to-1 if a nonzero damping constant (GE field on MATi Bulk Data entry) isnot found for any element.

NONLHT Output-integer-default=-1. Nonlinear heat transfer or differential stiffnessgeneration flag. Set to 1 if nonlinear heat transfer elements are detected.

On input:

>3 Compute geometric nonlinear effects


13/395


LUMPM Input-real-default=0.0. Lumping factor for electromagnetic mass.

MATCPX Input-integer-default=-1. Complex material properties flag forelectromagnetic elements.

KDGEN Input-integer-default=0. Differential or follower stiffness matrix generationflag. Usually the column number in UG to use in differential stiffnessmatrix generation. If KDGEN is negative, follower stiffness is generated.

TABS Input-real-default=0.0. Absolute temperature conversion. For example, setto 273.16 when specifying temperatures in Celsius or 459.69 in Fahrenheit.

SIGMA Input-real-default=0.0. The Stefan-Boltzmann constant. Used to computeradiant heat flux.

K6ROT Input-real-default=0.0. Normal rotational stiffness factor for CQUAD4and CTRIA3 elements.

LANGLE Input-integer-default=1. Large rotation calculation method:

1 Fimbal angle

2 Rotation vector

NOBKGG Input-integer-default=0. Slideline contact stiffness generation flag. Set to 1to generate slideline contact stiffness.

ALTSHAPE Input-integer-default=0. Specifies set of displacement functions inp-element analysis. ALTSHAPE=0 selects the MacNeal set and 1 selectsthe Full Product Space set.

PEXIST Input-logical-default=FALSE. P-element flag. Set to TRUE if p-elementsare present and to be processed.

FREQTYP Input-character-default= Frequency dependent element processing mode

ESTFCompute frequency dependent stiffness ESTNFCompute nominalfrequency dependent stiffness

FREQVAL Input-real-default=0.0. Frequency value for frequency dependent elementgeneration.

FREQWA Input-real-default=0.0. Parameter for electromagnetic analysis.

UNSYMF Input-character-default=NO Unsymmetric stiffness generation forslideline contact stiffness. If set to YES, stiffness matrix is unsymmetricfor slideline contact.

BADMESH Output-logical-default=FALSE. Bad geometry was detected.

UNUSED Unused parameter..



14/395


Remarks:

1. CSTM can be purged. MPT can be purged only if elements which do not reference anymaterial data are used. The DIT can be purged only if the material properties are nottemperature dependent.

2. If either of a matrix-dictionary data block pair is purged, that particular data block pairis not formed.

FOCOEL

Form contact elements

Format:

FOCOEL CASECC,BGPDT,CSTM,GEOM2,EST,MPT,CONTACT,SIL/CNELMS,GPECTC/S,N,NSKIP/S,N,OPTION/S,N,NLHEAT/S,N,CNTSET/S,N,NCELS/S,N,MAXS/S,N,MAXF/S,N,CTOL/S,N,AITK/S,N,MPER/S,N,RESET$

Input Data Blocks:

CASECC Table of case control command images.

BGPDT Basic grid point definition table.


GEOM2 Table of Bulk Data entry images related to element connectivity and scalarpoints.

EST Element summary table.

MPT Material property table.

CONTACT Contact data table.

SIL Scalar index list.

Output Data Blocks:

CNELM Contact element definition table.

GPECTC Grid point contact element connection table.

Parameters:

NSKIP Input-integer-no default. Record number in CASECC corresponding to the firstsubcase of the current boundary condition.



15/395


OPTION Input-integer-no default.

OPTION=1 to form contact elements.

OPTION=2 to form glue elements.

NLHEAT Input-logical-no default. TRUE if nonlinear heat transfer analysis.

CNTSET Output-integer. Contact set ID.

NCELS Output-integer. Number of contact elements.

MAXS Output-integer. Maximum number of iterations for status loop.

MAXF Output-integer. Maximum number of iterations for force loop.

CTOL Output-real. Contact force convergence tolerance.

AITK Output-integer. Aitken iterationflag.

MPER Output-integer. Minimum contact set percentage.

RESET Output-integer. Reset flag.

Remarks:

NCELS will equal -1 if error is encountered forming the contact element.

FOELCS

Form contact element stiffness.

Format:

FOELCS CNELM,BGPDT,CSTM,USET/ELCNST,ELCTST $S,N,NLHEAT$

Input Data Blocks:

CNELM Contact element definition table.

BGPDT Basic grid point data table.

CSTM Table of coordinate system transformation matrices.

USET Degree-of-freedom set membership table for g-set.

Output Data Blocks:

ELCNST Contact element normal stiffness table.



16/395


ELCTST Contact element tangential stiffness table.

Parameters:

NLHEAT Input-logical-no default. TRUE if nonlinear heat transfer analysis.

FRLG

Generates frequency-dependent loads or time-dependent loads

Generates frequency-dependent loads or time-dependent loads via Fourier transform forfrequency response analysis.

Format:

FRLG CASECC,USETD,DLT,FRL,GMD,GOD,DIT,PHDH/PPF,PSF,PDF,FOL,PHF,YPF/SOLTYP/OPT/S,N,FOURIER/S,N,APP $

Input Data Blocks:


USETD Degree-of-freedom set membership table for p-set

DLT Table of dynamic loads

FRL Frequency response list

GMD Multipoint constraint transformation matrix with extra points, m-set by ne-set

GOD Omitted degree-of-freedom transformation matrix with extra points, o-setby d-set

DIT Table of TABLEij Bulk Data entry images

PHDH Transformation matrix from d-set to modal coordinates

Output Data Blocks:

PPF Frequency response load matrix in the p-set

PSF Frequency response load matrix in the s-set

PDF Frequency response load matrix in the d-set

FOL Frequency response frequency output list

PHF Frequency response load matrix in the h-set (modal)



17/395


YPF Frequency response enforced motion matrix in the p-set

Parameters:

SOLTYP Input-character-no default. Solution method.

MODALModal: that is, compute PH

DIRECTDirect: that is, do not compute PH

OPT Input-integer-default=0. Processing option.0 = generate frequency-or time-dependent loads1 = generate load scale function only (FRL not used)2 = generate frequency-or time-dependent load table values only

FOURIER Output-integer-default=-1. Fourier transform. Set to 1 if TLOADi Bulk Data

entries are referenced by the DLOAD set identification number in CASECC.

APP Output-character-default=FREQ Dynamic load type. Set to FREQ, ifRLOAD1 or RLOAD2 entries are referenced. Set to TRAN, if TLOAD1 orTLOAD2 entries are referenced.

Remarks:

1. CASECC, FRL, and FOL cannot be purged.

2. DLT can be purged if PP, PS, PD, and PH are purged.

3. If USETD is not purged, PS cannot be purged if single-point constraints exist. Also, GMDand GOD cannot be purged if multipoint constraints or omitted degrees-of-freedom exist.

4. PHDH and PH cannot be purged if SOLTYP=MODAL

5. DIT cannot be purged if a dynamic load references TABLEDij Bulk Data entries.

6. PS, PD, and PH can be purged if USETD is also purged.

7. If TLOAD1 or TLOAD2 Bulk entries are referenced, the loads are computed and transformedto the frequency domain.

GP1

Performs basic geometry processing

Format:

GP1 GEOM1,GEOM2,geom3,GDNTAB,MEDGE,SGPDT,DYNAMIC/GPL,EQEXIN,GPDT,CSTM,BGPDT,SIL,VGF,GEOM3B,DYNAMICB/S,N,LUSET/S,N,NOCSTM/S,N,NOPOINTS/UNIT/UPERM/UPRMT/NUFLAG/SEID/NUMLM $



18/395


Input Data Blocks:

GEOM1 Table of Bulk Data entry images related to geometry

GEOM2 Table of Bulk Data entry images related to element connectivity and scalarpoints

GEOM3 Table of Bulk Data entry images related to static and thermal loads

UNUSED Unused and can be purged

GDNTAB Table of grid points generated for p-element analysis

MEDGE Edge table for p-element analysis

SGPDT Superelement basic grid point definition table

DYNAMIC Table of Bulk Data entry images related to dynamics

Output Data Blocks:

GPL External grid/scalar point identification number list

EQEXIN Equivalence table between external and internal grid/scalar identificationnumbers

GPDT Grid point definition table



SIL Scalar index list

VGF Fluid/structure partitioning vector with ones at the rows corresponding tofluid degrees-of-freedom

GEOM3B Table of Bulk Data entry images related to static and thermal loads withDAREA entry images converted to equivalent FORCE and MOMENT entryimages

DYNAMICB Table of Bulk Data entry images related to dynamics without DAREA entryimages

Parameters:

LUSET Output-integer-default=0. The number of degrees-of-freedom in the g-set.

NOCSTM Output-integer-no default. Number of coordinate systems found in GEOM1.Set to -1 if none are found.



19/395


NOPOINTS Output-integer-no default. Grid point flag. Set to -1 if none are found.Otherwise, set to 1.

UNIT Input-real-default=1.0. AUNIT record factor for electromagnetic analysis.

UPERM Input-real-default=1.2566E-06. Permeability for electromagnetic analysis.

UPRMT Input-real-default=8.8542E-12. Permittivity for electromagnetic analysis.

NUFLAG Input-integer-default=10. Unit type for electromagnetic analysis.

SEID Input-integer-default=-1. Superelement identification number.

NUMLM Input-integer-default=0. Number of Lagrange Multiplies from Lagrange rigidelements.

Remarks:

1. GP1 assembles a list of all grid and scalar points and places them in internal order, computescoordinate system transformation matrices, and transforms all grid points to the basiccoordinate system.

2. No output data block, except VGF, can be purged.

GPSTR2

Computes grid point stresses or strains

Computes grid point stresses or strains interpolated from element centroid stresses or strains.

Format:

GPSTR2 CASECC,EGPSF,BGPDT,OES1,OESNLXR/OGS1,EGPSTR,OGSR1/S,N,NOOGS1/S,N,NOEGPSTR/APP/NLSTRAIN $

Input Data Blocks:


EGPSF Table of element to grid point interpolation factors


OES1 Table of element stresses or strains in SORT1 format

OESNLXR Table of nonlinear element stresses in SORT1 format and appended for allsubcases



20/395


Output Data Blocks:

OGS1 Table of grid point stresses or strains in SORT1 format

EGPSTR Table of grid point stresses or strains for post-processing in the DBC module

OGSR1 Table of raw grid point stresses or strains in SORT1 format.

Parameters:

NOOGS1 Output-integer-default=-1. OGS1 creation flag. Set to 0 if OGS1 is created.

NOEGPSTR Output-integer-default=-1. EGPSTR creation flag. Set to 0 if EGPSTR iscreated.

APP Input-character-default=STATICS Analysis type. Allowable values are:

STATICS: statics

REIGEN: normal modes

TRANRESP: transient response

NLSTRAIN Logical-input-default=FALSE. Nonlinear strain data recovery, otherwise theflag at word 11 of OES1 takes precedence. Set to TRUE if nonlinear strainsare to be processed.

Remarks:

1. The GPSTRESS Case Control command controls the contents of OGS1.

2. The STRFIELD Case Control command controls the contents of EGPSTR.

MATMOD

Matrix modification

Transforms matrix or table data blocks according to one of many options into output matrixor table data blocks.

Format:

MATMOD I1,I2,I3,I4,I5,I6,I7,I8,I9/O1,O2/P1/P2/P3/P4/P5/P6/P7/P8/P9/P10/P11/P12/p13/p14/p15 $

Input Data Blocks:

Ii Input data blocks. I1 is required; I2 through I9 may not be necessarydepending on the value P1.



21/395


Output Data Blocks:

Oi Output data blocks

Parameters:

P1 Input-integer-no default. Option selection described in the table that follows.

P2, P3, P4 Input/output-integer-default=0. Parametric data depending on P1.

P5, P6 Input/output-real-default=0. Parametric data depending on P1.

P7 throughP11

Input/output-integer-default=0. Parametric data depending on P1.

P12 Input/output-character-default=blank. Parametric data depending on P1.

P13 Input/output-character-default= Parametric data depending on P1.



Remarks:

Each option corresponds to a different value of the first parameter, P1. The following summarytable provides descriptions of the options.

Option P1 = 1

Extract a block(s) of columns from a matrix.

Format:

MATMOD I1,,,,,,,,/O1,/1/STARTCOL/ENDCOL/COLINC $

Input Data Block:

I1 Any matrix (real or complex)

Output Data Block:

O1 Column vector containing column P2 of I1

Parameter:

STARTCOL Input-integer. Starting column number to extract from I1.

ENDCOL Input-integer. Ending column number to extract from I1



22/395


COLINC Input-integer. Column increment. Extract every COLINCth column betweenSTARTCOL and ENDCOL.

Remarks:1. If ENDCOL is zero, ENDCOL=STARTCOL.

2. If COLINC is zero, every column between STARTCOL and ENDCOL is extracted.

Examples:

1. Extract the seventh column from A and call it A7.

MATMOD A,,,,,,,,/A7, /1/7 $

2. Extract the third, fourth, and fifth columns from A.

MATMOD A,,,,,,,,/A345, /1/3/3 $

Option P1 = 2

Filter small magnitude terms of a matrix.

Format:

MATMOD I1,,,,,,,,/O1,/2/PURGE///FILTER $

Input Data Block:


Output Data Block:

O1 A copy of I1 with terms smaller in magnitude than P5 set to 0.0

Parameter:

PURGE Input-integer-default=0. If PURGE=0, and the input matrix has no nonzeroterms, the output matrix is purged. If PURGE 0, and the input matrix has no

nonzero terms, the output matrix is null.FILTER Input-real-default=0.0. Value offilter. Terms in I1 with an absolute magnitude

less than the absolute value of FILTER are set to zero.

Example:

Print terms in A smaller in magnitude than 100.0.

MATMOD A,,,,,,,,/AFILTER,/2////100.0 $ ADD A,AFILTER/ASMALL//-1 $MATPRN ASMALL// $



23/395


Option P1 = 3

Zeros out rows and columns of a matrix according to degree of freedom component number.

Format:

MATMOD I1,,,,,,,,/O1,/3/CODE $

Input Data Block:


Output Data Block:

O1 I1 with rows and columns according to DOFs described by P2.

Parameter:

CODE Input-integer-default=0. Packed DOF code that identifies rows and columns ofI1 to be made null (for example, 136 means that degrees of freedom 1, 3, and 6for each grid point are set to zero).

Remarks:

I1 is assumed to consist only of grid point degrees of freedom. A code of 345 simply zeros rowsand columns 3, 4, 5, 9, 10, 11, 15, 16, 17, and so on, of matrix I1. You should exercise caution

when selecting this option on a resequenced matrix.

Example:

Zero out degrees of freedom 1, 2, and 6 in stiffness matrix KGG.

MATMOD KGG,,,,,,,,/KGG1,/3/126 $EQUIVX KGGQ/KGG/ALWAYS $

Option P1 = 4

Replicates a matrix six rows by N columns row-wise to a g-row by N-column matrix. The inputmatrix is replicated for each grid point.

Format:

Form 1

MATMOD I1,SIL,,,,,,,/O1,/4 $

Form 2

MATMOD I1,,,,,,,,/O1,/4/LUSET $



24/395


Input Data Blocks:

I1 Any six-row by N-column matrix (real or complex)

SIL Scalar index list (SIL) table generated by the GP1 module

Output Data Block:

O1 g-row by N-column matrix containing I1 at every grid point

Parameter:

LUSET Integer-input-default=0. Used to supply the length of the g-set when SIL ispurged.

Remarks:

1. If SIL is purged, MATMOD uses LUSET for the size of the g-set. The assumption is madethat only grids exist in the g-set. LUSET must not be zero if SIL is purged.

2. If SIL is not purged, LUSET is ignored. I1 is inserted at the rows of every grid point. Scalarand extra points are ignored.

Option P1 = 5

Accepts a DMI matrix six rows by six columns and outputs either a g-row by g-column matrixwhere the input matrix is inserted at the diagonal 6x6 of each grid point, or a g-row by g-column

transform matrix.

Format:

Form 1 (Inserts 6x6 matrix along diagonal)

MATMOD I1,SIL,,,,,,,/O1,/5/LUSET/0 $

Form 2 (generates specified coordinate system to basic coordinate system transformation matrix)

MATMOD CSTM,SIL,,,,,,,/01,/5/LUSET/P3 $

Form 3 (generates a global-to-basic coordinate system transformation matrix)

MATMOD CSTM,SIL,BGPDT,,,,,,/01,/5//-1 $

Input Data Blocks:

I1 Any six-row by six-column matrix (real or complex)

SIL Scalar index list table output from the GP1 module



25/395



Output Data Block:

O1 g-row by N-column matrix containing I1 at every grid point

Parameter:

LUSET Integer-input-default=0. Used to supply the length of the g-set when SIL ispurged. This parameter is valid for Forms 1 and 2 only.

P3 Integer-input-default=0. Coordinate system identification number. Thisparameter is valid for Form 3 only.

Remarks:1. If P3 = 0, this option accepts a six row by six column matrix and the SIL table. The output is

a g-size square matrix containing the 6 x 6 input matrix along the diagonal at every gridpoint. Scalar and extra points contain 0.0 values.

If P3 > 0, the 6 x 6 single-precision matrix is the transformation matrix from the coordinatesystem with the coordinate ID = P3 to the basic system. If P3 points to a cylindrical orspherical coordinate system, the transformation location is at the origin of the P3 system.Scalar and extra points contain 1.0.

If P3 = -1, the 6 x 6 single-precision matrix is the global-to-basic transformation for each gridpoint. Scalar and extra points contain 1.0.

2. If SIL is purged, the MATMOD uses LUSET for the size of the g-set. The assumption is madethat only grids exist in the g-set. LUSET must not be zero if SIL is purged.

Examples:

Transform KGG to another coordinate system.

1. Assume TRANS to be a 6 x 6 transformation matrix and KGG was formed using only onecoordinate system (global coordinate system is the same at each grid point). TransformKGG using TRANS.

MATMOD TRANS,SIL,,,,,,,/TRANSG,/5 $TRNSP TRANS/TRANSGT $SMPYAD TRANSG,KGG,TRANSGT,,,/KGGPRIME/3 $

2. Assume KGG was formed using coordinate system 10 as the global coordinate system for allgrid points. Transform KGG to the basic coordinate system.

MATMOD CSTM,SIL,,,,,,,/TRANSG,/5//10 $TRNSP TRANSG/TRANSGT $SMPYAD TRANSG,KGG,TRANSGT,,,/KGGBASIC/3 $

3. Assume KGG was formed using arbitrary coordinate systems as the global coordinate systemfor each grid point. Transform KGG to the basic coordinate system.

MATMOD CSTM,SIL,BGPDT,,,,,,/TRANSG,/5//-1 $TRNSP TRANSG/TRANSGT $SMPYAD TRANSG,KGG,TRANSGT,,,/KGGBASIC/3 $



26/395


Transform KGGBASIC coordinate system 10.

MATMOD CSTM,SIL,,,,,,,/TRANS10,/5//10 $TRNSP TRANS10/TRANST10T $SMPYAD TRANS10T,KGGBASIC,TRANS10,,,/KGG10/3 $

Option P1 = 6

Find the maximum absolute value for each row over all columns of a matrix.

Format:

MATMOD I1,,,,,,,,/O1,/6 $

Input Data Block:

I1 Any matrix (real only)

Output Data Block:

O1 Column vector with terms that represent the absolute maximum over allcolumns of I1 for each row

Example:

Find the maximum displacements over all loading conditions.

MATMOD UG,,,,,,,,/UGMAX,/6 $

Option P1 = 7

Find the maximum absolute value for each column over all the rows of a matrix.

Format:

MATMOD I1,,,,,,,,/O1,/7 $

Input Data Block:

I1 Any matrix (real only)

Output Data Block:

O1 Column vector with a term that represents the absolute maximum over allrows of I1 for each column

Example:

Find the maximum displacement for each loading condition.



27/395


MATMOD UG,,,,,,,,/MAXDISP,/7 $

Option P1 = 8

Normalize matrix.

Format:

MATMOD I1,,,,,,,,/O1,/8////S,N,NORMREAL/S,N,NORMIMAG $

Input Data Block:


Output Data Block:

O1 Matrix shaped like I1 with every term divided by the term in I1 with thelargest absolute value

Parameters:

NORMREAL Output-real single precision. Set to the real part of the normalizing factor.

NORMIMAG Output-real single precision. Set to the imaginary part of the normalizing

factor if I1 is complex.

Option P1 = 9

Find the maximum (absolute magnitude) value of each of the three columns of the UHT-transientresponse solution matrix. (The columns of UHT represent displacement, velocity, andacceleration for each output time step.)

Format:

MATMOD I1,,,,,,,,/O1,/9 $

Input Data Block:

I1 Transient response solution matrix consisting of H rows by three columnmatrices (which represent displacement, velocity, and acceleration for eachoutput time step) appended to form a matrix H rows by three times thenumber of output time steps columns. (Real only).

[[{u1} {1} {a1}] [{u2} {2} {a2}]...[{ui} {#i} {ai}]]



28/395


Output Data Block:

O1 H-row by three column matrix of peak displacements, velocities, andaccelerations

Example:

Find and output maximum transient response.

MATMOD UHT,,,,,,,,/UHTMAX,/9 $DDRMM CASEXX,UHTMAX,PHIP1,,,,/OUPMAX,,,,/ $OFP OUPMAX,,,,// $

Option P1 = 10

Convert matrix I1 into its complex conjugate.

Format:

MATMOD I1,,,,,,,,/O1,/10 $

Input Data Block:


Output Data Block:

O1 Matrix shaped like I1 with every term converted to its complex conjugate

Example:

Find the magnitude of terms of a complex vector.

MATMOD CMPLX,,,,,,,,/CMPLXC,/10 $ ADD CMPLX,CMPLXC/CMPLXSQ///1 $DIAGONAL CMPLXSQ/CMLPXMAG/WHOLE/0.5 $

Option P1 = 11

Form a new Basic Grid Point Definition Table (BGPDT) with grid locations that are given by I1.

Format:MATMOD LOCVEC,BGPDT,,,,,,,/BGPDTN,/11 $

Input Data Blocks:

LOCVEC G-size vector with values that represent grid locations in the basic coordinatesystem




29/395


Output Data Block:

BGPDTN New BGPDT table with grid point locations that are displaced by LOCVEC

Example:

Build new BGPDT table based on the deformed state.

$ Convert displacement vector to basic$ coordinate systemPARAML CSTM//PRESENCE////S,N,NOCSTMIF (NOCSTM > -1) THEN $

MATMOD CSTM,SIL,BGPDT,,,,,,/TRANS,/5//-1 $MPYAD TRANS,UG/UGBASIC $

ELSE $EQUIVX UG/UGBASIC/ALWAYS $

ENDIF $$ Form vector containing new grid locations in$ basic coordinate system

VECPLOT UGBASIC,BGPDT,SCSTM,CSTM,,,,/LOCVEC/0/0/3 $$ Generate new BGPDTMATMOD LOCVEC,BGPDT,,,,,,,/BGPDTNEW,/11 $

Option P1 = 12

Perform simultaneous null column search on up to three matrices.

Format:

MATMOD I1,I2,I3,,,,,,/O1,O2/12/S,N,NONULL/NMATRIX $

Input Data Blocks:

I1, I2, I3 Matrices to search for simultaneous null columns (real or complex)

Output Data Blocks:

O1 Column vector which has 1.0 at those rows where all matrices selected forsearching have null columns

O2 Square symmetric matrix which has 1.0 on the diagonal of those columnswhere all matrices selected for searching have null columns

Parameters:

NONULL Output-integer. Set to -1 if no simultaneous null columns found; otherwise, itis set to the number of simultaneous null columns.

NMATRIX Input-integer-default=0. Number of matrices to be included in null columnsearch.



30/395


Remarks:

1. I2 and I3 can be purged.

2. O2 can be purged.

Example:

Search for simultaneous null columns in the g-size mass, damping, and stiffness matrices andremove rows and columns corresponding to these columns.

MATMOD MGG,BGG,KGG,,,,,,/PARTNULL,/12/S,N,NONULL/3 $IF (NONULL > 0) THEN $

PARTN MGG,PARTNULL,/MGGNEW,,,/-1 $EQUIVX MGGNEW/MGG/ALWAYS $PARTN BGG,PARTNULL,/BGGNEW,,,/-1 $EQUIVX BGGNEW/BGG/ALWAYS $PARTN KGG,PARTNULL,/KGGNEW,,,/-1 $EQUIVX KGGNEW/KGG/ALWAYS $

ENDIF $

Option P1 = 13

Copies any data block.

Format:

MATMOD I1,,,,,,,,/O1,/13 $

Input Data Block:

I1 Any data block (table or matrix)

Output Data Block:

O1 Copy of I1

Remarks:

COPY module is preferred over this option.

Option P1 = 14

Filter small magnitude terms from a matrix; more capabilities than option 2.

Format:

MATMOD I1,,,,,,,,/O1,/14/PURGE/UPLOW/STRTR/FILTER/RELFLT/TRUNC $



31/395


Input Data Block:

I1 Matrix to be filtered (real or complex)

Output Data Block:

O1 I1 is modified according to specifications set by parameters.

Parameters:

PURGE Input-integer-default=0. If PURGE = 0, and the input matrix has no nonzeroterms, the output matrix is purged. If PURGE 0, and the input matrix has nononzero terms, the output matrix is null.

UPLOW Input-integer-default=0. If UPLOW < 0, all lower triangular terms are set tozero. If UPLOW > 0, all upper triangular terms are set to zero. If UPLOW = 0,the action of this parameter is ignored.

STRTR Input-integer-default=0. If STRTR=0, string trailer is written.

FILTER Input-real single precision-default=0.0. Terms in I1 with an absolutemagnitude less than the absolute value of FILTER are set to zero.

RELFLT Input-real single precision-default=0.0. If RELFLT 0.0, terms of I1 are setto zero when

Equation 1-1.

I1 must be square for this option

TRUNC Input-integer-default=0. If TRUNC 0, truncate terms of I1 accordingly

Equation 1-2.

Remarks:

1. If FILTER = 0.0, O1 is a copy of I1.

2. If relative filtering is desired, FILTER must be zero (default).

3. If FILTER 0 or RELFLT < 0, the absolute filter technique is used.



32/395


4. If I1 is not square and the relative filtering option is selected, FILTER is set to RELFLT andthe absolute filtering technique is used. A user warning message is also issued.

Option P1 = 15

Not implemented.

Option P1 = 16

Put matrix into MATPOOL format, optional DMIG punched output.

Format:

MATMOD I1,I2,,,,,,,/O1,/16/PNDMIG//TYPOUT////////CCHAR $

Input Data Blocks:

I1 Any matrix of g-rows arranged in external (ascending by GRID ID) sequence.(Real or complex).

I2 EQEXIN table from module GP1

Output Data Block:

O1 Table data block in MATPOOL format containing I1

Parameters:

PNDMIG Input-integer-default=0. If PNDMIG 0, I1 is printed in DMIG format tothe punch output file (.pch).

TYPOUT Input-integer-default=0. Default is to set DMIG precision to machineprecision. The default can be overridden by specifying:

1 Real single-precision output

2 Real double-precision output

3 Complex single-precision output

4 Complex double-precision output

CCHAR Input-character-default = blank. Continuation characters to be usedfor DMIG output. If nonblank continuation characters are specified, amaximum of 9999 DMIG entries can be printed for any single matrix.Only the first two characters of the nonblank mnemonic are used for thecontinuation string.



33/395


Remarks:

1. I1 must be a g-row size matrix arranged in external sequence (see PARAM, OLDSEQ, andMATGEN Option 9). The rows are always labeled with the external sequence (grid or scalarIDs and component numbers). If the input matrix form is 1 (square) or 6 (symmetric) the

columns are also labeled with the external sequence, and the IFO entry on the generatedDMIG entry is set to 1 or 6. If the form is 6, only the terms in one triangle are output. If theinput matrix form is 2 (rectangular form), the columns are labeled sequentially, starting withunity. The IFO is set to 9. If the form is not 1, 2, or 6, the module returns with no output.

2. EQEXIN table must not be purged.

Example:

Output KGG and PG matrices in MATPOOL formatted table and punch to DMIG Bulk Dataentries (MATGEN Option 9 is being used to resequence them from internal to external sort).

MATGEN EQEXIN/INTEXT/9/0/LUSET $MPYAD INTEXT,KGG,/KGGE/1 $MPYAD KGGE,INTEXT,/KGGEXT $MATMOD KGGEXT,EQEXIN,,,,,,,/MATPOOL1,/16/1 $MPYAD INTEXT,PG,/PGEXT/1 $MATMOD PGEXT,EQEXIN,,,,,,,/MATPOOL2,/16/1 $

Option P1 = 17

Generate a g-size partitioning vector from a user-defined set of grid and/or scalar points or froma user-supplied bit position that designates one of the USET sets.

Format:

MATMOD EQEXIN,USET,SIL,CASECC,,,,,/CP/17/UBIT/SETFLG/S,N,NOCP////////SETSTR1/SETSTR2/SETSTR3/SETSTR4 $

Input Data Blocks:


USET Degree-of-freedom set membership table for g-set

SIL Scalar index list


Output Data Block:

CP Column partitioning vector



34/395


Parameters:

UBIT Input-integer-default=0. Obsolete method for set selection. For a moreuser-friendly method, use SETSTRi. In order to select specific sets for UBIT,

add the corresponding decimal equivalent numbers from the table below. Forexample, sets R, O, and M, UBIT=8+4+1=13. For supersets, add the decimalequivalent numbers of the mutually exclusive sets which are contained inthe superset. For example, set S combines the SB and SG set and thereforeUBIT=1024+512=1536. The presence of any grid point degree of freedom in theassociated sets causes all degrees of freedom associated with that grid point to begiven a value of 1.0 in the output vector.

SetName DecimalEquivalentNumber

Q 4194304

BE 2097152

C 1048576

K 262144

SA 131072

E 2048

SB 1024

SG 512

R 8

O 4

BF 2

M 1

SETFLG Input-integer-default=0. If SETFLG 0, SETFLG selects a set of grid pointidentification numbers of which all degrees of freedom associated with eachpoint are assigned a value of 1.0 in the corresponding row of CP. If no SETcommand is found, the UBIT parameter is used.

If SETFLG>0, the PARTN=SID Case Control command selects the SETcommand.

If SETFLG


35/395


SETSTRi Input-character-default= Set name string. SETSTR1 through SETSTR4form a single string of set name(s) and is 32 characters in length.The set names must be separated by a plus sign, "+". For example,SETSTR1=M+R+N+SGand SETSTR2=+A+Qspecifies the m, r, n, sg, a,

and q sets.

Remarks:

1. None of the data blocks can be purged.

2. UBIT entry must be a legitimate value.

Example:

Generate a partitioning vector from a set of grid points defined in the Case Control Section.

MATMODEQEXIN,USET,SIL,CASECC,,,,,/VECX,/17/128/1 $

If no set had been selected in the Case Control Section, the vector VECX would have beengenerated using the a-set degrees of freedom since UBIT=128.

The Case Control Section contains:

SET 10 = 1 THRU 50PARTN = 10

Option P1 = 18

Insert or modify a GEOM3 table temperature record.

Format:

MATMOD GEOM,GPL,UG,,,,,,/GEOM3T,/18/SID $

Input Data Blocks:


GPL External grid/scalar point identification number list

UG Temperature matrix in g-set

Output Data Block:

GEOM3T GEOM3 table with new or modified temperatures

Parameter:

SID Input-integer. Temperature set to be modified or added.



36/395


Remarks:

1. This option should only be used in heat transfer analysis.


3. Only grid temperature records (not elements) are modified.

Example:

Put data from the temperature vector UG into a record for SID = 30.

MATMOD GEOM3,GPL,UG,,,,,,/GEOM3NEW,/18/30 $

Option P1 = 19

Extract a temperature vector from a GEOM3 table.

Format:

MATMOD GEOM3,EQEXIN,,,,,,,/UGT,/19/SID

Input Data Blocks:



Output Data Block:

UGT Updated temperature matrix in g-set

Parameter:

SID Input-integer. Temperature set to extract.

Remarks:

1. This option should only be used in heat transfer analysis.


3. Only grid temperature records (not elements) are extracted.

Example:

Extract the temperature vector UGN for SID = 40.

MATMOD GEOM3,EQEXIN,,,,,,,/UGN,/19/40 $



37/395


Option P1 = 20

Print the magnitude of the largest terms of up to six matrices.

Format:

MATMOD I1,I2,I3,I4,I5,I6,,,/O1,/20////S,N,SUM $

Input Data Blocks:

I1 though I6 Any matrix (real or complex)

Output Data Block:

O1 Dummy output data block

Parameter:

SUM Output-real. Sum of the absolute values of the largest terms in the inputmatrices.

Remarks:

1. Any input matrix can be purged.

2. All computations are performed in single or double precision depending on the matrix type.The magnitudes of the largest terms and SUM are converted to single precision for output.

Option P1 = 21

Extracts the components of a factor matrix and converts them to a standard form suitable forinput to any matrix module.

Format:

MATMOD LD,,,,,,,,/T,LP/21 $

Input Data Block:

LD Lower triangular factor/diagonal matrix

Output Data Blocks:

T Diagonal from symmetric decomposition

LP Lower triangular [L] and permutation matrix appended together



38/395


Remarks:

Symmetric decomposition forms the equivalent matrix representation of a symmetric matrix.

[A] =[P}T [L] [D] [L]T[P]

where P is a permutation matrix (row and column interchange used to improve efficiency), L,a lower triangular matrix, and D, a tridiagonal matrix. Option 21 extracts P, L, and D andconverts them to a standard form, suitable for input to any matrix module. L and P are appendedcolumn-wise in output D.

Example:

Extract components of factor matrix LLL.

MATMOD LLL,,,,,,,,/TT,LP/21 $PARAML KLL//TRAILER/1/S,N,NL $TYPE PARM,,I,N,TUNL $T U N L = 2 * N L $

MATGEN ,/V21/6/TUNL/NL/NL $PARTN LP,V21,/LL,,PP,/0 $

Option P1 = 22

Generate special aeroelasticity matrix with modified trailers.

Format:

MATMOD MKLIST,Qij,,,,,,,/QijL/22 $

Input Data Blocks:

MKLIST Table of Mach number and reduced frequency pairs

Qij Aerodynamic matrix

Output Data Block:

QijL Special aerodynamic matrix with modified trailers

Option P1 = 23

Determines type of eigenvalue analysis requested and optionally extracts values from theselected EIGR or EIGRL Bulk Data entry.

Format:

MATMOD CASECC,DYNAMIC,,,,,,,/,/23/ S,N,METHTYP/S,N,LANCZOS/S,N,EIGRVALI/S,N,EIS,N,NFOUND/ICASE////S,N,EIGRFLD $



39/395


Input Data Blocks:


DYNAMIC Table of Bulk Data entry images related to dynamics

Parameters:

METHTYP Output-integer. Set to 1 if Inverse Power or Lanczos method selected;otherwise, set to zero.

LANCZOS Output-integer. Set to -1 if Lanczos method selected; otherwise, set to zero.

EIGRVALI Output-integer-default=0. Extracted integer value from the EIGR or EIGRLentry.

EIGRVALR Output-real-default=0.0. Extracted real value from the EIGR or EIGRL entryNFOUND Output-integer-default=0. EIGR* entry found flag; 0 if entry was found and -1

if entry was not found.

ICASE Input-integer-default=1. Case Control record number which contains theMETHOD command.

EIGRFLD EIGRFLD Input/output-character-default= Field name of EIGR or EIGRLentry. EIGRFLD is also an output if the field value is a character string.

Examples:

1. Determine method type on the EIGR entry.

MATMOD CASECC,DYNAMIC,,,,,,,/,/23/S,N,METHTYP $IF ( METHTYP=1 ) THEN $

MESSAGE //LANCZOS OR SINV IS SELECTED.$ELSE $

MESSAGE //GIVENS OR HOUSEHOLDER IS SELECTED.$ENDIF $

2. Extract the F2 field value from the EIGR entry:

MATMOD CASECC,DYNAMIC,,,,,,,/,/23////S,N,F2///////F2$

3. Extract the NORM field value from the EIGR entry:

NORM=NORM$ initialize and will change on output

MATMOD CASECC,DYNAMIC,,,,,,,/,/23///////////S,N,NORM $

Option P1 = 24

Generate a square matrix that has a 1.0 at the intersection of every null row and null column ofI1, I2, and I3 simultaneously.

Format:

MATMOD I1,I2,I3,,,,,,/O1,/24/S,N,NOOUT/NMATX/S,N,NRNENC $



40/395


Input Data Blocks:

I1, I2, I3 Square, commensurate matrices (real or complex)

Output Data Block:

O1 Square matrix that has a 1.0 at the intersection of every null row and nullcolumn of I1, I2, and I3 simultaneously

Parameters:

NOOUT Output-integer. Set to -1 if O1 is null.

NMATX Input-integer-default=1. Number of input matrices to be used for search,

starting from the first input.

NRNENC Output-integer. Set to -1 if the number of null rows does not equal thenumber of null columns.

Remarks:

Any two input matrices can be purged.

Example:

Add a unit value to the stiffness matrix for degrees of freedom that have no associated mass,damping, or stiffness. This is usually done to prevent potential singularities during directtransient and frequency analyses.

MATMOD MAA,BAA,KAA,,,,,,/KAAX,/24/S,N,NOADD/3/S,N,NRNENC $IF (NRNENC < 0) THEN $

MESSAGE //ERROR: MATRICES ARE NOT/ SYMMETRIC $

EXIT $ENDIF $IF (NOADD < -1) THEN

ADD KAA,KAAX/KAANEW $EQUIVX KAANEW/KAA/ALWAYS $

ENDIF $

Option P1 = 25

Generate vectors that have 1.0 corresponding to each null row and null column in I1, I2, and I3simultaneously.

Format:

MATMOD I1,I2,I3,,,,,,/O1,O2/25/S,N,NOOUT/NMATX/S,N,NRNENC///S,N,SYM $



41/395


Input Data Blocks:

I1, I2, I3 Square, commensurate matrices (real or complex)

Output Data Blocks:

O1 Vector that has 1.0 corresponding to each null row in I1, I2, and I3simultaneously. See Remark 1.

O2 Vector that has 1.0 corresponding to each null column in I1, I2, and I3simultaneously.

Parameters:

NOOUT Output-integer. Set to -1 if both output vectors are null, set to zero otherwise.

NMATX Input-integer-default=1. Number of input matrices to be used for search,starting from the first input.

NRNENC Output-integer. Set to -1 if number of null rows does not equal the number ofnull columns; otherwise zero.

SYM Output-integer. Set to -1 if I1, I2, and I3 are symmetric; otherwise zero. SeeRemark 2.

Remarks:

1. Any two input matrices can be purged.

2. If I1, I2, and I3 are symmetric, O2 is purged.

Example:

Remove null rows and columns from matrix A.

MATMOD A,,,,,,,,/RPARTN,CPARTN/25/S,N,NOOUT/1////S,N,SYM $IF (NOOUT > -1) THEN $

PARTN A,CPARTN,RPARTN/ANEW/SYM $EQUIVX ANEW/A/ALWAYS $

ENDIF $

Option P1 = 26Used internally for development testing.

Option P1 = 27

Convert a diagonal matrix (form 3) to a symmetric matrix (form 6).

Format:

MATMOD I1,,,,,,,,/O1,/27 $



42/395


Input Data Block:

I1 Diagonal matrix of form 3 (real or complex)

Output Data Block:

O1 Symmetric matrix of form 6 containing diagonal terms of I1

Remarks:

1. Form 3 matrices are not output by any module. They are only allowed as input by theINPUTT2, INPUTT4, and DMIIN modules.

2. The SMPYAD, MPYAD, and ADD modules do not accept form 3 matrices. The matricesshould now be converted to form 6 before use in these modules.

Example:

DMIIN DMI,DMINDX/A3,,,,,,,,, $MATMOD A3,,,,,,,,/A6,/27 $

where A3 is the DMI matrix defined by the Bulk Data entries

DMI,A3,0,3,1,1,,4,1DMI,A3,1,2,2.0,3.0,4.0

and A6 is the matrix

Option P1 = 28

Convert the first column of a matrix to a symmetric matrix (form 6) with the terms of the firstcolumn along the diagonal and off-diagonal terms set to zero.

Format:

MATMOD I1,,,,,,,,/O1,/28 $

Input Data Block:

I1 Any matrix of form 1, 2, or 6 (real or complex)



43/395


Output Data Block:

O1 Symmetric matrix (form 6) with terms of the first column of I1 along thediagonal and off-diagonal terms set to zero.

Example:

MATMOD B,,,,,,,,/BDIAG,/28 $

If B is the matrix

Figure 1-1.

BDIAG is

Figure 1-2.

Option P1 = 29

Used internally for development testing.

Option P1 = 30

Print data blocks or a portion of data blocks as a table of hexadecimal values.

Format:

MATMOD I1,I2,I3,I4,I5,I6,,,/,/30/BBLK/EBLK $

Input Data Block:

Ii Any data block (matrix or table)

Parameters:

BBLK Input-integer-default=1. Beginning GINO block number.

EBLK Input-integer-default=-1. Ending GINO block number. Default value impliesthe total number of blocks.



44/395


Option P1=31

Writes the bit map of a matrix to the punch file.

Format:

MATMOD MAT,,,,,,,,/,/31/MAXSIZ $

Input Data Block:

MAT Any matrix

Output Data Block:

None.

Parameter:

MAXSIZ Input-integer-default=0. Maximum size of the bit map matrix (row and/orcolumn)

Option P1=32

Convert tables created by DRMH1 into DTI Bulk Data entry format and write to the punch file.Also converts DRMH1 directory tables in DTI Bulk Data entry format into DRMH1 outputtable format.

Format:

MATMOD TXY,,,,,,,,/TOUT,/32/CONVERT $

Input Data Block:

TXY DRMH1 directory table in DTI or table data block format

Output Data Block:

TOUT DRMH1 directory table in table data block format or DTI format

Parameter:

CONVERT Input-integer-default=0. Convert option.

0 Table data block format to DTI format

1 DDTI to table data block format



45/395


Remarks:

Table record 3 is all character and reading DTI entries produces all numbers. Therefore,CONVERT=1 converts the DTI numbers to character values.

Option P1=33

Create a single column matrix from the frequency response output list table, FOL. Thefrequencies are also converted to radian units.

Format:

MATMOD FOL,,,,,,,,/FOLMAT,/33 $

Input Data Block:

FOL Frequency response frequency output list

Output Data Block:

FOLMAT Matrix of frequencies in radian units

Parameters:

None.

Option P1=34

Extract the real and imaginary parts of complex matrix into two real matrices.

Format:

MATMOD CMAT,,,,,,,,/RMAT,IMAT/34//PREC $

Input Data Block:

CMAT Complex matrix

Output Data Blocks:

RMAT Matrix containing real part of CMAT

IMAT Matrix containing imaginary part of CMAT

Parameter:

PREC Input-integer-default=0. Precision of output matrices.



46/395


47/395


2 Absolute value in ascending order

NKEYS Input-integer-default=1. Duplicate value sort option specification.

1 Single key sort

2 Double key sort to maintain original order of terms in case of duplicateterms

Remarks:

For complex matrices, only SORTOPT=2 or -2 is allowed.

Example:

Given the input matrix, IM, generate an algebraic ascending order sort. The input matrix andits sorted order (algebraically ascending) are:

IM =

The MATMOD call looks like this:

P1=35 $MATMOD

IM,,,,,,,,/LIST,BOOL/S,N,P1//S,N,SORTED////2 $

and the output matrix LIST for NKEYS=2 contains

LIST =

For NKEYS=1, it is equally likely that the indices for equal values are in a different order. Forexample, if the MATMOD call statement were

P1=35 $MATMOD

IM,,,,,,,,/LIST,BOOL/S,N,P1//S,N,SORTED////1 $

the output matrix LIST for NKEYS=1 contain either



48/395


LIST =

or

LIST =

since there are duplicate terms (0.0s) in the input matrix column and a single key sort was used.

The Boolean square matrix contains unit values in the appropriate positions so that it can beused to create the sorted input matrix by means of a simple matrix multiply as in:

MPYAD BOOL,IM,/IMS/1 $

producing the sorted IM matrix, IMS, as

IMS =

Remarks:

1. If the matrix input into this option contains more than one column and SORTBOOL is usedsubsequently to operate on this matrix, all columns have their rows re-ordered according tothe sort obtained from the column processed by the MATMOD operation. In this case, thecolumn that was selected during the MATMOD operation to produce the sorted ordering isguaranteed in sort. Other columns may or may not have their rows in sorted order.

2. The output data blocks are in machine precision, regardless of the precision of the inputmatrix.



49/395


3. NKEYS=2 provides a more repeatable sort in the presence of equal values in the input, atthe cost of longer run times. A test on a typical vector showed a difference of a factor ofapproximately ten. If repeatability is not essential, NKEYS=1 is the preferred choice.

Option P1=36Reduce the GRID record in the GEOM1 table to the entries corresponding to grid identificationnumbers specified in a Case Control set.

Format:

MATMOD GEOM1,CASECC,,,,,,,/GEOM1R,/36/GRIDSET/S,N,NOGEOM1 $

Input Data Blocks:


CASECC Table of Case Control command selections

Output Data Block:

GEOM1R GEOM1 table with reduced GRID record

Parameters:

GRIDSET Input-integer. SET Case Control command identification number which

contains a list grid point identifi

cation numbers.

NOGEOM1 Output-integer. Processing status flag.

+1 No grid data found matching gridset.

0 GRIDSET found and contents match some GRIDs in GEOM1.

-1 GRIDSET found and contents matches all GRIDs in GEOM1.

Remarks:

Only the GRID record is processed and all other GEOM1 records are copied as is to GEOM1R.

Option P1=37

Reduce the element and SPOINT records in the GEOM2 table to the entries corresponding toelement or SPOINT identification numbers specified in a Case Control set.

Format:

MATMOD GEOM2,CASECC,,,,,,,/GEOM2R,/37/ELEMSET/GRIDSET/S,N,NOGEOM2 $



50/395


Input Data Blocks:


CASECC Table of Case Control command selections

Output Data Block:

GEOM2R GEOM2 table with reduced element record

Parameters:

ELEMSET Input-integer. SET Case Control command identification number whichcontains a list element point identification numbers.

GRIDSET Input-integer. SET Case Control command identification number whichcontains a list SPOINT identification numbers.

NOGEOM2 Output-integer. Processing status flag.

+1 No element and SPOINTs found matching ELEMSET and GRIDSET.

0 ELEMSET and GRIDSET found and contents match some elements andSPOINTs in GEOM2.

-1 ELEMSET and GRIDSET found and contents match all elements andSPOINTs in GEOM2.

Option P1=38

Reduce the records in the EST table to the entries corresponding to element numbers specified ina Case Control set.

Format:

MATMOD EST,CASECC,,,,,,,/ESTR,/38/ELEMSET/GRIDSET/S,N,NOEST $

Input Data Blocks:

EST Table of Bulk Data entry images related to geometry

CASECC Table of Case Control command selections.

Output Data Block:

ESTR EST table with reduced records



51/395


Parameters:

ELEMSET Input-integer. SET Case Control command identification number thatcontains a list element point identification numbers.

GRIDSET Input-integer. SET Case Control command identification number whichcontains a list grid point identification numbers.

NOEST Output-integer. Processing status flag.

+1 No element found matching contents of ELEMSET

0 ELEMSET found and contents match some elements in EST

-1 ELEMSET and contents match all elements and SPOINTs in EST

Option P1=39

Remove and identify explicit zero terms in a matrix.

Format:

MATMOD I1,,,,,,,,/O1,O2/39/S,N,NOXPLZER $

Input Data Block:

I1 Any matrix

Output Data Blocks:

O1 Matrix I1 with explicit zero terms removed

O2 Matrix containing a 1.0 at the row and column, where an explicit zero wasfound in I1

Parameter:

NOXPLZER Output-integer. Explicit zero existence flag. Set to -1 if no explicit zeros arefound.

Option P1=40

Remove unused q-set degrees-of-freedom from USET.

Format:

MATMOD USET,VAXW,,,,,,,/USETN,/40 $



52/395


Input Data Block:


VAXW Column vector for a-set degrees of freedom which has a 1.0 at those locationswhere a-set structural matrices are null.

Output Data Blocks:

USETN New degree-of-freedom set membership table for g-set (with unused q-setdegrees of freedom removed).

Option P1=41

Generate a partitioning vector from a set of grid points defined in the Bulk Data Section on aROTORG card.

Format:

MATMOD EQEXIN,USET,SIL,CASECC,DYNAMICS,,,,/PARTG,/41/0/ROTGID/NOVEC $

Input Data Block:

EQEXIN Equivalence table between external and internal grid/scalaridentification numbers.



CASECC Table of Case Control command images.

DYNAMICS Table of Bulk Data entry images related to dynamics.

Output Data Blocks:

PARTG Column partitioning vector.



53/395


Parameters:

Parameter 2 not used at this time.

ROTGID Input-integer. ROTGID must be less than zero. The absolute value ofROTGID selects the ROTORG Bulk Data card containing a set of grid pointidentification numbers. Six degrees of freedom will be associated to each pointby assigning six values of 1.0 in the corresponding row of PARTG.

NOVEC Output-integer. NOVEC will be set to -1 if the partitioning vector is null orcannot be generated. Otherwise, it will be set to zero.

Remarks:

1. None of the data blocks may be purged.

2. ROTGID parameter must be a negative integer value less than zero.

Example:

MATMOD EQEXIN,USET,SIL,CASECC,DYNAMICS,,,,/PARTG,/41/0/-11/NOVEC $

The Bulk Data Section contains:

ROTORG 11 1 3 THRU 5 7

The resulting partitioning vector PARTG:

1 0 1 1 1 0 1

Option P1=42

Generate a partitioning vector from a set of data defined in the Bulk Data Section on SOL 200.

Format:

MATMOD BGPDT,EQEXIN,SIL,CASECC,EDOM,EPT,EPTABF,GEOM2,GEOM2A/RGPV,/42/NEPT/NOSE $

Input Data Block:

BGPDT Basic grid point definition table.






54/395


EDOM Table of Bulk Data entries related to design sensitivity and optimization.

EPT Table of Bulk Data entry related to element properties.

EPTABF Family of tables of designed property attributes.

GEOM2 Table of Bulk Data entry images related to element connectivity andscalar point.

GEOM2A Table of secondary Bulk Data entry images related to elementconnectivity and updated for the current p-level.

Output Data Blocks:

RGPV Row partitioning vector.

Parameters:

NEPT Input-integer-default=0. The number of EPOINT used in FE model.

NOSE Input-integer default=0. NOSE should be set to 1 if superelement is used.

Remarks:

1. The partitioning vector RGPV is generated from the following data defined in the Bulk Data:

- grids and elements defined in set

- properties defined in DESVAR

- grids, elements, properties defined in DRESP.

2. Option 42 does not support shape optimization.

3. Option 42 does not support material optimization.

4. Option 42 does not support XYPLOT, PLOT, POST output request.

Example:

The partitioning vector RGPV includes 0 or 1:

1 0 1 1 1 0 1



55/395


Option P1=43

Generate a partitioning vector from a set of modal indices defined in the Case Control Section ona SET card referenced to by a MODSEL card.

Format:

MATMOD EQEXIN,USET,SIL,CASECC,DYNAMICS,,,,/PARTMS1,/43/NEIGV/MODSEL/NOVEC $

Input Data Block:





DYNAMICS Table of Bulk Data entry images related to dynamics.

Output Data Blocks:

PARTMS1 Column partitioning vector.

Parameters:

Parameter 2 not used at this time.

NEIGV Input-integer. When the Case Control MODSEL card references a CaseControl "SET = ALL" card the value of NEIGV will be used to determine theset of grid point identification numbers. Thus the set will contain grids 1 to nwhere n equals the value of NEIGV.

MODSEL Input-integer. MODSEL must be less than zero. The absolute value ofMODSEL selects a Case Control MODSEL card which points to a Case ControlSET card containing a set of grid point identification numbers. Each point inthe set will be assigned a value of 1.0 in the corresponding row of PARTMS1.

NOVEC Output-integer. NOVEC will be set to -1 if the partitioning vector is null orcannot be generated. Otherwise, it will be set to zero.



56/395


Remarks:

1. None of the data blocks may be purged.

2. MODSEL parameter must be a negative integer value less than zero.

Example:

MATMOD EQEXIN,USET,SIL,CASECC,DYNAMICS,,,,/PARTMS1,/43/0/-10/NOVEC $

The Case Control Section contains:

SET 10 = 1, 3 THRU 5, 7MODSEL = 10The resulting partitioning vector PARTMS1:

1 0 1 1 1 0 1

Option P1=44

Create a g-set partitioning vector for sparse data recovery based on the output requests inthe case control section. This option is only valid for frequency and transient response outputrequests.

Format:

MATMOD BGPDTS,SILS,CASES,USETS,EPTS,GEOM2S,EQEXINS,XYCDB,GMTG/PARTV,/44/NOEPT/NOSE/IWHO $

Input Data Blocks:

BGPDTS Basic grid point definition table.

SILS Scalar index list.

CASES Table of Case Control command images.

USETS Table of degree-of-freedom sets.

EPTS Element property table.

GEOM2S Table of bulk data entries related to element connectivity.

EQEXINS Table of external and internal grid/scalar numbers.

XYCDB Table of XY plotting commands.



57/395


GMTG Transpose of GM (multipoint constraint transformation matrix; m-set byn-set) where the n-set dimension has been expanded to g-set size.

Output Data Block:

PARTV Sparse data recovery partitioning vector (length is g-set). Values of 1 indicateretained dof; values of 0 indicate ignored dof.

Parameters:

NOEPT Input-integer default=0. Number of EPOINTs to process.

NOSE Input-integer default=0. NOSE should be set to 1 if superelement is used.

IWHO Input-integer default=0. IWHO should be set to 0 for enforced motion; set to 1for SOL 111 or SOL 112 with superelements.

Remarks:

1. All data blocks must exist except for GMTG. If any of the required data blocks are purged,the module will not produce PARTV.

2. GMTG must exist in order to process m-set degrees-of-freedom.

3. If the user requests ALL for any of the supported case control output requests, PARTV willcontain all values of 1 and the SPARSEDR system cell (421) will be reset to 0, indicatingthat sparse data recovery is OFF.

4. The supported case control output requests are:

a. DISPLACEMENT

b. VELOCITY

c. ACCELERATION

d. OLOAD

e. STRESS/ELSTRESS

f. STRAIN/ELSTRAINg. FORCE/ELFORCE

h. EDE

I. EKE

j. ESE

k. XYPEAK

l. XYPRINT

m. XYPLOT

n. XYPUNCH



58/395


Option P1=45

Construct XCASECC by selecting the desired records of CASECC.

Format:

MATMOD CASECC,,,,,,,,/XCASECC,/45/P2/P3/P4 $

Input Data Blocks:


Output Data Block:

XCASECC Table of Case Control command images for each processor in DMP environment

Parameters:

P2 Input-integer default=0. Skip P2 records from the beginning of CASECC.

P3 Input-integer default=0. Copy the subsequent P3 records to XCASECC after

P2 records.

P4 Input-integer default=0. The remaining number of records after P3 records iscopied to XCASECC.

Remarks:

Example:

Extract 2, 3, and 4th record from CASECC that has 6 records.

MATMOD CASECC,,,,,,,,/XCASECC,/45/1/3/2 $

MTRXIN

Converts matrices input on DMIG Bulk Data entries to matrix data blocks.

Format:

Form 1 Simplified (CASECC is purged)



59/395


MTRXIN

Form 2 Case Control Command Selection of stiffness, mass, and damping (or square) matrices(IOPT=1 for K2GG, and so on, and IOPT=0 for K2PP, and so on, and TF)

MTRXIN

Form 3 Case Control Command selection of load (or rectangular) matrix (IOPT=2)

MTRXIN

Form 4 Selection of DMIK, DMIJ and DMIJI by data block names MATKi, MATJi, and MATJIi.

MTRXIN

Form 5 Selection of stiffness, mass, damping, and loads (or square) matrices by K2PNAM, andso on, input parameter values (IOPT=10 through 12).



60/395


MTRXIN

Form 6 - Selection of DMIK, DMIJ, and DMIJI matrices by the MATNAMi input parametervalues (IOPT=13 through 15)

MTRXIN

Input Data Blocks:


MATPOOL Table of Bulk Data entry images related to hydroelastic boundary, heattransfer radiation, virtual mass, DMIG, and DMIAX entries


EQDYN Equivalence table between external and internal grid/scalar/extra pointidentification numbers. (EQEXIN appended with extra point data)

TFPOOL Table of TF Bulk Data entry images

AEBGPDTK Basic grid point definition table for the aerodynamic ks-set degrees of freedom.

AEBGPDTJ Basic grid point definition table for the aerodynamic js-set degrees of freedom



61/395


AEBGPDTI Basic grid point definition table for the aerodynamic js-set interferencedegrees of freedom

Output Data Blocks:

NAMEi Matrices defined on DMIG Bulk Data entries

K2GG, etc. Matrices defined on DMIG Bulk Data entries and referenced by the K2GG,M2GG, B2GG, K42GG, K2PP, M2PP, B2PP, or P2G Case Control commands

MATPi Matrices defined on DMIG Bulk Data entries and intended for the p-set

MATGi Matrices defined on DMIG Bulk Data entries and intended for the g-set

RMATG Rectangular matrix defined on DMIG Bulk Data entries and can have anarbitrary number of columns but g-set rows, similar to P2G

MATKi Matrices defined on DMIK Bulk Data entries

MATJi Matrices defined on DMIJ Bulk Data entries

MATIi Matrices defined on DMIJI Bulk Data entries

Parameters:

LUSET Input-integer-no default. The number of degrees-of-freedom in the g-set.

LUSETD Input-integer-no default. The number of degrees-of-freedom in the p-set.

NONAMEi Output-integer-default=-1. NAMEi generation flag. Set to +1 if NAMEi isgenerated; -1 otherwise.

NOK2, etc. Output-integer-default=-1. K2GG, and so on, generation flag. Set to +1 ifK2GG, and so on is generated; -1 otherwise.

IOPT Input-integer-default=0. Case Control command selection flag.

0 No Case Control command selection (see Form 1) or K2GG, and so on,and TFL Case Control command selection (see Form 2)

1 K2GG, and so on, Case Control command selection (see Form 2)

2 P2G Case Control command selection (see Form 3)

3 DMIK selection by output data block name (see Form 4)

4 DMIJ selection by output data block name (see Form 4)

5 DMIJI selection by output data block name (see Form 4)

10 K2PP, M2PP, and B2PP selection by input parameter value (see Form 5)



62/395


11 K2GG, M2GG, B2GG, and K42GG selection by input parameter value(see Form 5)

12 P2G selection by input parameter value (see Form 5)

13 DMIK selection by input parameter value (see Form 6)

14 DMIJ selection by input parameter value (see Form 6)

15 DMIJI selection by input parameter value (see Form 6)

LKSETLJSET,LISET

Output-integer-default=0. Size of ks-set, js-set, and inteference js-setextracted from the AEBGPTK, AEBGPDTJ and AEBGPDTI tables.

NOMATi Output-integer-default=1. Generation flag. Set to +1 if MAT* isgenerated; 1 otherwise.

MATNAMi Input-character-default= Matrix name found on DMIG, DMIJ, DMIK,and DMIJI Bulk Data entries.

TFLID Input-integer-default=0. Transfer function set identification number.TFLID is ignored if IOPT=3, 4, 5, 13, 14, or 15.

NFEXIT Input-logical-default=TRUE. Termination flag. If FALSE, do not issueUser Fatal Message 2070 and do not terminate the module if the matrixis not found.

Remarks:

1. Any output data block can be purged.

2. Form 1 is used to input matrices from DMIG entries named in the DMAP statement outputsection. No Case Control commands are required.

3. Forms 2 and 3 are used to select the matrices with Case Control commands: K2GG, M2GG,B2GG, K42GG, K2PP, M2PP, B2PP, or P2G. -2GG matrices are of dimension g by g.-2PP matrices are of dimension p by p. The P2G matrix has g-rows, with the number ofcolumns determined by the several methods used to input rectangular matrices describedon the DMIG entry.

4. if the output data blocks are specified on a CALL statement and the DMIIN module isspecified in the subDMAP referenced by the CALL statement, the data block name specifiedon the CALL statement must be the same as the name specified on the DMIIN module.

Examples:

1. Assume the Bulk Data contains two DMIG matrices, named M1 and M2, which referencegrid and/or scalar points only. The following set of DMAP instructions generate these twomatrices in matrix format, multiply them, and print the result.

MTRXIN ,,MATPOOL,EQEXIN,,/Ml,M2,,/LUSET/S,N,NOMl/S,N,NOM2$IF (NOM1 > -1 AND NOM2 > -1) THEN $

MPYAD M1,M2,/PRODUCT $MATPRN PRODUCT//$



63/395


ENDIF $

2. Assume the Bulk Data contains two DMIG matrices, MASS and STIFF, which reference gridand/or scalar points only. The following Case Control and DMAP instructions generate thesetwo matrices in matrix format and add them to the structural mass and stiffness.

Case Control:

M2GG = MASSK2GG = STIFF

DMAP instructions:

MTRXIN CASECC,MATPOOL,EQEXIN,,/STIFF,MASS,,/LUSET/S,N,NOSTIFF/S,N,NOMASS///1 $

IF (NOSTIFF > -1) THEN $ ADD KGG,STIFF/KGGNEW $EQUIVX KGGNEW/KGG/ALWAYS $

ENDIF $IF (NOMASS > -1) THEN $

ADD MGG,MASS/MGGNEW $

EQUIVX MGGNEW/MGG/ALWAYS $ENDIF

NXNADAMS

Creates an ADAMS MNF for a superelement.

Creates an ADAMS Modal Neutral File (MNF) for a superelement. The output is based on thedefinitions in the ADAMSMNF case control command and the DTI,UNITS bulk data entry.

Format:

NXNADAMSUNITS,CASES,BGPDTS,GEOM2,GEOM4,USET,LAMA,PHIG,MGGDIAG,PCDB,OGPWG,OGSR1,OGSTRR1,CSTMS//SEID/FLXONL/FLXERR/WTMASS/GRDPNT $

Input Data Blocks:

UNITS UNITS data block from the DTI,UNITS bulk data entry

CASES Case Control table associated with superelement

BGPDTS Basic Grid Point Definition Table associated with superelement

GEOM2 Table of Bulk Data entries related to element connectivity

GEOM4 Table of Bulk Data entries related to constraints

USET Table of degree-of-freedom sets

LAMA Eigenvalue summary table for superelement

PHIG Matrix of eigenvectors (g-set size) corresponding to LAMA

MGGDIAG Column vector of diagonal values from superelement mass matrix (g-set size)



64/395


PCDB Plot Control Data Block from OUTPUT(PLOT) case control

OGPWG Output table of Grid Point Weight Generator of superelement

OGSR1 Output table of grid point stresses of superelement

OGSTRR1 Output table of grid point strains of superelement

CSTMS Coordinate System Transformation Matrices for superelement

Parameters:

SEID Input-integer-default=0. Superelement ID number (0 for residual-onlyanalysis).

FLXONL Output-integer-no default. Value of FLEXONLY keyword from ADAMSMNF

case control command. Options are:

0: Continue with solution of residual structure

1: Do not perform solution of residual structure

FLXERR Output-integer-no default. Error flag. Options are:

0: No error

1: Error occurred creating MNF. Terminate processing.

WTMASS Input-real-default=1.0. Value of WTMASS parameter from

PARAM,WTMASS,value.

GRDPNT Input-integer-default=-1. Value of GRDPNT parameter fromPARAM,GRDPNT,value.

Remarks:

The MNF naming convention is as follows: jid_SEID.mnf, where jid is the job ID of the run(that is, the name of the job input file) and SEID is the superelement ID number (that is, theSEID parameter). The location of the created MNF is the same as the jid.f06 file.

PARAML

Sets parameters from a data block



65/395


Format:

PARAML DB/DBNAME/P1/S,N,P2/S,N,P3/S,N,P4/S,N,P5/S,N,P6/S,N,SET1/S,N,F1/S,N,SET2/S,N,F2/S,N,SET3/S,N,F3/S,N,SET4/S,N,F4/

S,N,SET5/S,N,F5/S,N,SET6/S,N,F6/S,N,SET7/S,N,F7/S,N,SET8/S,N,F8/S,N,SET9/S,N,F9/S,N,SET10/S,N,F10/S,N,SET11/S,N,F11/S,N,SET12/S,N,F12 $

Input Data Block:

DB Any matrix or table

Output Data Block:

DBNAME Any data block. Used only when P1 = NAME; otherwise DBname must notbe specified.

Parameters:

P1 Input-character-no default. Only the first 4 characters are required. Forexample, PRES and PRESENCE are equivalent.

P2 Input/output-integer-default = 1

P3 Input/output-integer-default = 1

P4 Output-real-default = 0.0

P5 Output-integer-default = 0.0

P6 Output-real-default = 0.0

SETi Input/output-character-default =

Fi Output-integer-default = 0

The following describes the various options and their formats. The meaning and usage of

parameters P2 through P6, SETi, and Fi depend on the value of P1. Under all options, P5 is setto -1, if the input data block does not exist, and no other parameters are set.

Option P1 = BULK

Check for the presence of Bulk Data entry records by examining the trailer bits of its IFPmodule related table.



66/395


Format:

PARAML IFPDB//BULK//////BULKNM1/S,N,BULKFG1/BULKNM2/S,N,BULKFG2/BULKNM3/S,N,BULKFG3/BULKNM4/S,N,BULKFG4/

BULKNM5/S,N,BULKFG5/BULKNM6/S,N,BULKFG6/BULKNM7/S,N,BULKFG7/BULKNM8/S,N,BULKFG8/BULKNM9/S,N,BULKFG9/BULKNM10/S,N,BULKFG10/BULKNM11/S,N,BULKFG11/BULKNM12/S,N,BULKFG12 $

Input Data Block:

IFPDB Table with trailer bits indicating existence of Bulk Data entry records

Parameters:

BULKNMi Input-character. Name of Bulk Data entry.

BULKFGi Output-integer. Set to -1 if Bulk Data entry exists.

Remarks:

To determine which table contains BULKNMi, see the data block description.

Example:

Check for the presence of the rigid elements.

PARAML GEOM4//BULK//////RBE1/S,N,RBE1/RBE2/S,N,RBE2/

RBE3/S,N,RBE3/RROD/S,N,RROD/RBAR/S,N,RBAR/RTRPLT/S,N,RTRPLT/RSPLINE/S,N,RSPLINE $

Released Dmaps

Documents