nasread

Contents Functions

nasread

Purpose

Read results from outputs of the MSC/NASTRAN finite element code. This function is part of FEMLink.

Syntax

 out = nasread('FileName','Command')

Description

nasread reads bulk data deck (NASTRAN input), direct reading of model and result information in OUTPUT2 and OUTPUT4 files generated using NASTRAN PARAM,POST,-i cards. This is the most efficient and accurate method to import NASTRAN results for post-processing (visualization with feplot, normal model handling with nor2ss, ...) or parameterized model handling with upcom. Results in the .f06 text file (no longer supported).

Available commands are

Bulk file

model=nasread('FileName','bulk') reads NASTRAN bulk files for nodes (grid points), element description matrix, material and element properties, and coordinate transformations, MPC, SPC, DMIG, SETS, ...

Use 'BulkNo' for a file with no BEGIN BULK card. Unsupported cards are displayed to let you know what was not read. You can omit the 'bulk' command when the file name has the .dat or .bdf extension.

Each row of the bas.bas output argument contains the description of a coordinate system.

The following table gives a partial conversion list. For an up to date table use nas2up('convlist')

NASTRAN	SDT
CELAS1, CELAS2, RBAR	celas
RBE2	rigid
RBE3	rbe3 in Case
CONROD	bar1
CBAR, CBEAM, CROD	beam1
CBUSH	cbush
CSHEAR	quad4
CONM1, CONM2	mass2
CHEXA	hexa8, hexa20
CPENTA	penta6, penta15
CTETRA	tetra4, tetra10
CTRIA3, CTRIAR	tria3
CTRIA6	tria6
CQUAD4, CQUADR	quad4
CQUAD8	quadb

Details on properties are given under naswrite WritePLIL. NASTRAN Scalar points are treated as standard SDT nodes with the scalar DOF being set to DOF .01 (this has been tested for nodes, DMIG and MPC).

OUTPUT2 binary

model=nasread('FileName','output2') reads output2 binary output format for tables, matrices and labels. You can omit the output2 command if the file names end with 2. The output model is a model data structure described in section 7.6. If deformations are present in the binary file, the are saved OUG(i) entries in the stack (see section 7.8). With no output argument, the result is shown in feplot.

Warning: do not use the FORM = FORMATTED in the eventual ASSIGN OUTPUT2 statement.

The optional out argument is a cell array with fields the following fields

.name	Header data block name (table, matrix) or label (label)
.dname	Data block name (table, matrix) or NASTRAN header (label)
.data	cell array with logical records (tables), matrix (matrix), empty (label)
.trl	Trailer (7 integers) followed by record 3 data if any (for table and matrix), date (for label)

Translation is provided for the following tables

GEOM1	nodes with support for local coordinates and output of nodes in global coordinates
GEOM2	elements with translation to SDT model description matrix (see bulk command).
GEOM4	translates constraints (MPC, OMIT, SPC) and rigid links (RBAR, RBE1, RBE2, RBE3, RROD, ...) to SDT model description matrix
GPDT	with use of GPL and CSTM to obtain nodes in global coordinates
KDICT	reading of element mass (MDICT, MELM) and stiffness (KDICT, KELM) matrix dictionaries and transformation of a type 3 superelement handled by upcom. This is typically obtained from NASTRAN with PARAM,POST,-4. To choose the file name use Up.file='FileName';Up=nasread(Up,'Output2.op2');
MPT	material property tables
OUG	transformation of shapes (modes, time response, static response, ...) as curve entries in the stack (possibly multiple if various outputs are requested).
	Note : by default deformations are in the SDT global coordinate system (basic in NASTRAN terminology). You may switch to output in the local (global in NASTRAN terminology) using PARAM,OUGCORD,GLOBAL.
	To avoid Out of Memory errors when reading deformations, you can set use a smaller buffer sdtdef('OutOfCoreBufferSize',10) (in MB). When too large, def.def is left in the file and read as a v_handle object that lets you access deformations with standard indexing commands. Use def.def=def.def(:,:) to load all.
	To get the deformation in the stack use calls of the form def=stack_get(model,'curve','OUG(1)','get')
OEE	tables of element energy
OES	tables of element stresses or strains.

This translation allows direct reading/translation of output generated with NASTRAN PARAM,POST commands simply using out=nasread('FileName.op2'). For model and modeshapes, use PARAM,POST,-1. For model and element matrices use PARAM,POST,-4 or PARAM,POST,-5 (see BuildUp command below).

BuildUp,BuildOrLoad

A standard use of FEMLink is to import a model including element matrices to be used later with upcom. You must first run NASTRAN SOL103 with PARAM,POST,-4 to generate the appropriate .op2 file (note that you must include the geometry in the file, that is not use PARAM,OGEOM,NO). Assuming that you have saved the bulk file and the .op2 result in the same directory with the same name (different extension), then

 Up=nasread('FileName.blk','buildup')

reads the bulk and .op2 file to generate a superelement saved in FileName.mat.

It is necessary to read the bulk because linear constraints are not saved in the .op2 file during the NASTRAN run. If you have no such constraints, you can read the .op2 only with Up=upcom('load FileName);Up=nasread(Up,'FileName.op2').

The BuildOrLoad command is used to generate the upcom file on the first run and simply load it if it already exists.

 nasread('FileName.blk','BuildOrLoad') % result in global variable Up

OUTPUT4 binary

out=nasread('FileName','output4') reads output4 binary output format for matrices. The result out is a cell array containing matrix names and values stored as MATLAB sparse matrices.

All double precision matrix types are now supported. If you encounter any problem, ask for a patch which will be provided promptly.

Output4 text files are also supported with less performance and no support for non sequential access to data with the SDT v_handle object.

Supported options

-full : assumes that the matrix to be read should be stored as full (default sparse).
-transpose : transpose data while reading.
-hdf : save data in a hdf file. Reading is performed using buffer (sdtdef('OutOfCoreBufferSize',100) for a 100MB buffer). It is useful to overcome the 2GB limit on 32 bit Matlab: see sdthdffor details about how to build v_handle on hdf file.

.f06 output (obsolete)

ASCII reading in .f06 files is slow and often generates round-off errors. You should thus consider reading binary OUTPUT2 and OUTPUT4 files, which is now the only supported format. You may try reading matrices with nasread('FileName','matprt'), tables with nasread('F','tabpt') and real modes with

 [vector,mdof]=nasread('filename','vectortype')

Supported vectors are displacement (displacement), applied load vector (oload) and grid point stress (gpstress).