Nas101 S+N+B Seminar
Nas101 S+N+B Seminar
Nas101 S+N+B Seminar
Frimley, Camberley
Surrey, United Kingdom
Tel: (+44)-1276-67 1000
Fax: (+44)-1276-69 1111
MSC/NASTRAN Support
Tel: 1-800-336-4858
Fax: 323-259-4999
Tokyo, Japan
Tel: 81-3-3505-0266
Fax: 81-3-3505-0914
Munich, Germany
Tel: (+49)-89-43 19 87 0
Fax: (+49)-89-43 61 716
Paris, France
Tel: (+33)-1-69 36 69 36
Fax: (+33)-1-69 36 45 17
Rome, Italy
Tel: (+39)-6 5 91 64 50
Fax: (+39)-6 5 91 25 05
Moscow, Russia
Tel: (+7095)-236-61-77
Fax: (+7095)-232-3575
Madrid, Spain
Tel: (+34)-915-560919
Fax: (+34)-915-567280
BASIC MSC.NASTRAN
Linear Static and Normal Modes Analysis
SEMINAR NOTES
NA*V70.5*Z*Z*Z*SM-NAS101-NT1
June 1998
DISCLAIMER
The concepts, methods, and examples presented in this text are for educational purposes only and
are not intended to be exhaustive or to apply to any particular engineering problem or design. The
MSC.Software Corporation assumes no liability or responsibility to any person or company for direct
or indirect damages resulting from the use of any information contained herein.
Printed in U.S.A.
1999 by The MacNeal-Schwendler Corporation
All rights reserved.
Seminar
Date
Location
Instructor
Your Name
Title
Company
Your feedback is very important in helping MSC meet the training needs of its clients. Your specific
written comments are especially helpful. Please use the following scale:
1 = poor
Rating
2 = marginal
3 = average
Category
4 = good
5 = very good
Comments
Effectiveness of instructor in
explaining technical concepts.
Organization of material in a
clear and understandable
manner.
(Continued on back)
Rev 5/93
1 = poor
Rating
2 = marginal
3 = average
Category
4 = good
5 = very good
Comments
Impression of hotel
accommodations. At which
hotel did you stay? Any
problems?
Would you recommend this instructor for future classes? Why or why not? ____________________
________________________________________________________________________________
Rev 5/93
TABLE OF CONTENTS
TITLE
1.0
SECTION
INTRODUCTION TO MSC . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-1
COMPANY OVERVIEW . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
WHAT IS MSC/NASTRAN?. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
MSC CLIENT SUPPORT. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
ADDITIONAL SOURCES OF INFORMATION . . . . . . . . . . . . . . . . . .
MSC TRAINING AND EDUCATION . . . . . . . . . . . . . . . . . . . . . . . . . .
2.0
3.0
2-1
2-3
2-4
2-5
2-6
2-7
2-8
2-10
2-14
2-16
2-18
2-20
2-22
4.0
1-1
1-2
1-4
1-6
1-7
3-1
3-2
3-4
3-5
3-6
4-1
4-3
4-4
4-7
SECTION
SMALL FIELD FORMAT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
FREE FIELD FORMAT IN THE BULK DATA . . . . . . . . . . . . . . . . . . .
LARGE FIELD FORMAT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
GENERAL RULES FOR BULK DATA INPUT . . . . . . . . . . . . . . . . . . .
CONTINUATION ENTRIES. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
GENERATION/REPLICATION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
GENERATION/REPLICATION OF
CONTINUATION ENTRIES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
GENERATION/REPLICATION EXAMPLES . . . . . . . . . . . . . . . . . . . .
COMMON ENTRY FORMAT ERRORS . . . . . . . . . . . . . . . . . . . . . . .
THE NASTRAN STATEMENT. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
FILE MANAGEMENT SECTION (FMS) . . . . . . . . . . . . . . . . . . . . . . .
5.0
5-1
5-2
5-3
5-4
5-5
5-6
5-7
5-9
5-10
7.0
4-16
4-17
4-19
4-20
4-21
6.0
4-9
4-10
4-11
4-12
4-13
4-15
6-1
6-2
6-4
6-9
6-11
6-12
6-14
6-17
6-19
7.1-1
7.1-2
7.2-1
7.3-1
SECTION
MATERIAL PROPERTIES. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
MSC/NASTRAN ELEMENTS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
ZERO-DIMENSIONAL (SCALAR) ELEMENT LIBRARY. . . . . . . . . . .
ONE-DIMENSIONAL ELEMENT LIBRARY. . . . . . . . . . . . . . . . . . . . .
CALCULATION OF TORSIONAL CONSTANT J
FOR PROD, PBAR, AND PBEAM ENTRIES. . . . . . . . . . . . . . . . .
CROD EXAMPLE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
CBAR ELEMENT DEFINITION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
CBARAO . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
CBEAM ELEMENT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
TWO-DIMENSIONAL ELEMENT LIBRARY . . . . . . . . . . . . . . . . . . . .
TWO-DIMENSIONAL ELEMENTS . . . . . . . . . . . . . . . . . . . . . . . . . . .
REFERENCES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
CQUAD4 ELEMENT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
CQUAD4 ELEMENT DEFINITION . . . . . . . . . . . . . . . . . . . . . . . . . . .
QUAD4 ELEMENT DEFINITION ALTERNATE PROPERTY ENTRY . . . . . . . . . . . . . . . . . . . . . . . .
THREE-DIMENSIONAL ELEMENT LIBRARY . . . . . . . . . . . . . . . . . .
RECOMMENDED USES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
CHEXA . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
MSC/NASTRAN LOADS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
FORCE AND MOMENT ENTRIES . . . . . . . . . . . . . . . . . . . . . . . . . . .
EXAMPLE INPUT FOR FORCE AND
MOMENT ENTRIES. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
DISTRIBUTED LOADS (PLOAD) IN MSC/NASTRAN . . . . . . . . . . . .
EXAMPLE INPUT FOR PLOAD1 ENTRY . . . . . . . . . . . . . . . . . . . . . .
EXAMPLE INPUT FOR PLOAD4 ENTRY . . . . . . . . . . . . . . . . . . . . . .
GRAV ENTRY . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
THERMAL EFFECTS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
EXAMPLE OF APPLYING THERMAL EFFECTS . . . . . . . . . . . . . . . .
LOAD ENTRY . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
THE LOAD REQUEST TREE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
8.0
7.4-1
7.5-1
7.6-1
7.7-1
7.7-6
7.7-8
7.7-14
7.7-37
7.7-44
7.8-1
7.8-2
7.8-3
7.8-4
7.8-6
7.8-18
7.9-1
7.9-2
7.9-3
7.10-1
7.10-2
7.10-3
7.10-5
7.10-6
7.10-10
7.10-12
7.10-14
7.10-15
7.10-18
7.10-19
8.1-1
8.1-2
8.1-5
8.1-6
8.1-7
8.2-1
8.2-8
8.2-10
PARAMETERS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9-1
INTRODUCTION. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
PARAMETERS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
PRINTOUT GENERATED BY PARAM,PROUT,1. . . . . . . . . . . . . . . .
PRINTOUT GENERATED BY PARAM,GPECT,1 . . . . . . . . . . . . . . . .
PRINTOUT GENERATED BY PARAM,EST,1 . . . . . . . . . . . . . . . . . .
10.0
11.0
10-1
10-3
10-7
10-11
10-12
10-16
10-18
10-20
12.0
9-1
9-2
9-4
9-5
9-6
11-1
11-2
11-3
11-5
11-6
11-7
11-8
11-10
11-16
12-1
12-2
12-3
12-4
12-5
12-6
12-7
12-8
12-9
12-11
12-14
12-15
12-16
12-18
12-21
12-22
12-23
12-24
12-25
12-26
12-27
14.0
15.0
B-1
B-3
B-4
B-6
B-7
B-8
B-9
vii
SECTION 1
INTRODUCTION TO MSC
COMPANY OVERVIEW . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-1
WHAT IS MSC/NASTRAN? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-2
MSC CLIENT SUPPORT. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-4
ADDITIONAL SOURCES OF INFORMATION . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-6
MSC TRAINING AND EDUCATION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-7
COMPANY OVERVIEW
___________________________________________________________________________________________
1-1
WHAT IS MSC/NASTRAN?
MSC/NASTRAN is a general-purpose, finite element analysis
program capable of solving a wide variety of engineering problems
including:
Linear static analysis
Static analysis with geometric and material nonlinearity
Transient analysis with geometric and material nonlinearity
Normal modes and buckling analysis
Direct and modal complex eigenvalue analysis
Direct and modal frequency analysis and random response
Direct and modal transient analysis (including response spectrum
analysis)
Linear static and vibration analysis with cyclic symmetry
Linear and nonlinear steady-state heat transfer
Transient heat transfer
Aeroelasticity
Multilevel superelements
1-2
1-3
1-4
1-5
1-6
1-7
Element
Modeling
1-8
Techniques
Using
SECTION 2
INTRODUCTION TO FINITE ELEMENT THEORY
INTRODUCTION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-1
MATRIX THEORY OF STRUCTURAL ANALYSIS
(HISTORICAL SKETCH) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-3
OVERVIEW OF THE MATRIX METHOD OF STRUCTURAL ANALYSIS . . . . . . 2-4
FUNDAMENTAL STRUCTURAL ANALYSIS REQUIREMENTS . . . . . . . . . . . . . 2-5
BASIC EQUATION OF THE DISPLACEMENT METHOD . . . . . . . . . . . . . . . . . . 2-6
INTERPRETATION OF ELEMENTAL STIFFNESS
MATRIX [K] AND STIFFNESS COEFFICIENTS (kij) . . . . . . . . . . . . . . . . . . . . . 2-7
DISCRETIZATION OF CONTINUOUS STRUCTURE. . . . . . . . . . . . . . . . . . . . . . 2-8
EXAMPLE: DIRECT FORMULATION OF ROD ELEMENT
STIFFNESS MATRIX BY DISPLACEMENT METHOD. . . . . . . . . . . . . . . . . . . . . 2-10
OTHER EXAMPLES OF ELEMENTAL STIFFNESS MATRICES. . . . . . . . . . . . . 2-14
GLOBAL STIFFNESS MATRIX. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-16
CONSTRAINING THE STRUCTURE - RIGID BODY MOTION . . . . . . . . . . . . . . 2-18
OVERVIEW OF FINITE ELEMENT SOLUTION PROCEDURE . . . . . . . . . . . . . . 2-20
REFERENCES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-22
INTRODUCTION
Energy
Boundary Element
Numerical
Methods
Finite Difference
Finite Element
Classical (analytical) methods
Exact and Approximate. Offer a high degree of insight, but the
problems are difficult or impossible to solve for anything but simple
geometry and loadings
2-1
INTRODUCTION (Cont.)
Numerical methods
Energy: Minimizes an expression for the potential energy of the
structure over the whole domain. This approach is feasible only for
simple configurations.
Approximates functions satisfying the
Boundary Element:
governing differential equations not the boundary conditions.
Problem size is reduced because elements represent only the
boundary of the domain. However, the application of this method
relies on knowing the fundamental solution to the governing
equations. This can sometimes be difficult to obtain.
Finite Difference: Replaces governing differential equations and
boundary conditions with corresponding algebraic finite difference
equations.
Problems occur when attempting to represent
irregularly shaped structures or complicated boundary conditions.
Finite Element: Offers considerable problem generality and ease
of solution by permitting irregularly shaped elements. Some prior
knowledge of the overall behavior of the structure is required to
guide model generation and element selection.
2-2
2-3
Displacement
Method
basic
2-4
F = 0 , M = 0
1.
2.
3.
A homogeneous, isotropic material [E] reduces to two indepedent material constants E and .
For such a material under uniaxial load,
= E
These three conditions can be used to generate a system of
equations in which the displacements are unknown (the displacement
method).
2-5
1.
2.
3.
The compatibility condition correlates the external grid point displacements to end deformations of the elements.
The force-displacement relationship is established between the member end forces and displacements, and between the grid point forces
and displacements. The stiffness matrix [ K ] is used to relate the
forces acting on the structure and the displacements resulting from
these forces in the following manner:
{F} = [K]{u}
where
2-7
2-8
uy
x
ux
uz
z
3 translations
{ ux , uy , uz }
3 rotations
{ x , y , z }
{ u } = vector of displacements = { u x u y u z x y z }
Note:
2-9
2
u1
F2
u2
Fx = F1 + F2 = 0
or
(1)
F2 = F1
2-10
(2)
2-11
(3)
(4)
F1
A
Step 5: Relate force to displacement
Substitution of Equations 2 and 3 into Equation 4 yields:
EA
F 1 = x A = E x A = -------- ( u 2 u 1 )
L
or
(5)
EA
EA
F 1 = -------- u 2 -------- u 2
2
L
2-12
(6)
EA 1 1
-------u
L 1 1
1
u2
F1
=
F2
or
(7)
{F} = [K]{u}
where
2-13
x2
T1
T2
L
Two rotational degrees of freedom about x ( x1 , x2 )
GJ
- 1 1
T
1 ------L
1 1
=
T2
K
where G
x1
x2
2-14
y1
y2
L
2 rotations about z ( z1 , z2 )
2 translations in ( y 1 , y 2 )
{P}
3L
3L
3L 2L 3L L
6 3L
3L L
3L
3L 2L
[K]
y1
z1
y2
z2
P y1
M z1
2EI
P y2 = -------3
L
M z2
{u}
2-15
ka
1
u1 , F1
kb
u2 , F2
u3 , F3
ka
0 u1
ka
F1
F2 = ka ( ka + kb ) kb u2
0
kb
kb u3
F3
2-16
ka ka u1
F1
and
F
k
u
k
2
a a 2
F1
F2 =
F3
kb kb u2
F2
F
kb kb u3
3
ka
ka
ka
ka + kb
kb
kb
kb
u1
u2
u3
2-17
{F}
2-18
Note:
2-19
Solve matrix
displacements.
equation
2-20
{F}
[K]{u}
and
stresses
for
from
2-21
REFERENCES
Matrix Structural Analysis
H. C. Martin
Introduction to Matrix Methods of Structural Analysis
McGraw-Hill Book Company
1966
J. S. Przemieniecki
Theory of Matrix Structural Analysis
McGraw-Hill Book Company and Dover Publications
1968
M. F. Rubinstein
Matrix Computer Analysis of Structures
Prentice-Hall, Inc.
1966
M. D. Vanderbilt
Matrix Structural Analysis
Quantum Publishers, Inc.
1974
2-22
REFERENCES (Cont.)
Finite Element Method
K. J. Bathe (theoretical)
Finite Element Procedures in Engineering Analysis
Prentice-Hall, Inc.
1982
R. D. Cook (theoretical, applications)
Concepts and Applications of Finite Element Analysis, 2nd Edition
John Wiley and Sons, Inc.
1981
R. H. Gallagher (applications)
Finite Element Analysis Fundamentals
Prentice-Hall, Inc.
1975
O. C. Zienkiewicz and R. L. Taylor (theoretical)
The Finite Element Method, 4rd Edition
McGraw-Hill Book Company
1989
R. H. MacNeal (theoretical applications)
Finite Elements: Their Design and Performance
Marcel Dekker, Inc.
1994
A. O. Cifuentes (MSC/NASTRAN applications)
Using MSC/NASTRAN: Statics and Dynamics
Springer-Verlag
1989
2-23
2-24
SECTION 3
INTRODUCTION TO FINITE ELEMENT MODELING
GOAL OF FINITE ELEMENT ANALYSIS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-1
INPUT TO THE FINITE ELEMENT MODEL . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-2
OUTPUTS FROM THE FINITE ELEMENT ANALYSIS. . . . . . . . . . . . . . . . . . . . . 3-4
MODELING GUIDELINES. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-5
UNITS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3-6
X
uy
ux
Element
Level
[ K ]e { u }e = { f }e
The behavior of the structure is obtained by considering the collective
behavior of the discrete elements.
3-1
3-2
Geometry
Element Connectivities
Element Properties
Material Properties
Constraints
Loads and
Enforced
Displacements
3-3
3-4
MODELING GUIDELINES
3-5
UNITS
MSC/NASTRAN knows nothing about units.
Example
Input
Units
English
Metric
inches
mm
Elastic Modulus
lb/in2
N/mm2
Applied Moments
inch-lb
mm-N
Applied Forces
lb
Corresponding Output
Displacements
Units
inches
mm
Stresses
lb/in2
N/mm2
3-6
Must Be
Consistent
UNITS (Cont.)
System of
Units
Mass
Length
Time
Force
English
lb sec2/ft
ft
sec
lb
32.174 ft/sec2
SI
kg
sec
Newton
9.81 m/sec2
F = Ma:
Note:
weight
mass(M) = -------------------g
Example:
density
0.283 lb/in
Mass density = weight
----------------------------------------- = -----------------------------------------------------g
ft 12 in
of mild steel
32.174 ------------- --------------
2
sec 1 ft
2
lb sec
= 7.324E 4 -------------------4
in
3-7
3-8
SECTION 4
INTRODUCTION TO MSC/NASTRAN INPUT FILE
MSC/NASTRAN INPUT FILE - OVERALL VIEW . . . . . . . . . . . . . . . . . . . . . . . . . 4-1
DELIMETER ENTRIES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-3
SAMPLE MSC/NASTRAN INPUT FILE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-4
FORMAT OF MSC/NASTRAN BULK DATA ENTRIES. . . . . . . . . . . . . . . . . . . . . 4-7
SMALL FIELD FORMAT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-9
FREE FIELD FORMAT IN THE BULK DATA . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-10
LARGE FIELD FORMAT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-11
GENERAL RULES FOR BULK DATA INPUT. . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-12
CONTINUATION ENTRIES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-13
GENERATION/REPLICATION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-15
GENERATION/REPLICATION OF CONTINUATION ENTRIES . . . . . . . . . . . . . . 4-16
GENERATION/REPLICATION EXAMPLES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-17
COMMON ENTRY FORMAT ERRORS. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-19
THE NASTRAN STATEMENT. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-20
FILE MANAGEMENT SECTION (FMS) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4-21
ENDDATA
BULK DATA SECTION
BEGIN BULK
CASE CONTROL SECTION
CEND
EXECUTIVE CONTROL SECTION
ID A, B $ OPTIONAL
FMS
4-1
2.
3.
4.
4-2
DELIMETER ENTRIES
ID A, B
CEND
BEGIN BULK
ENDDATA
Note:
4-3
E = 30.0 E6 psi
t = 0.25 inches
= 0.33
y
8
4
5
z
1
10 inches
2
x
2
10 inches
= grid ID
= element ID
Note:
4-4
4-5
4-6
CROD
Defines a tension-compression-torsion element.
Format:
1
CROD
EID
PID
G1
G2
13
21
23
10
Example:
CROD
12
Field
Contents
EID
PID
G1, G2
Remarks:
1. Element identification numbers should be unique with respect to all other element
identification numbers.
2. See CONROD for alternative method of rod definition.
3. Only one element may be defined on a single entry.
BCD
.7E1
0.7+1
.70+1
7.+0
70. -1
The first field on the first line of an entry is the entry name.
4-8
10
GRID
10
7.5
8.6
4-9
456
4-10
10A
16
16
16
16
7.5
8.6
*GRID10
GRID*
10
Second Entry
1B
10B
16
16
16
16
*GRID10
9.0
456
4-11
1 0
All real numbers, including zero, must contain a decimal point WARNING - COMMON ERROR!
A blank field will be interpreted as a real number or integer number
as required.
4-12
CONTINUATION ENTRIES
Definition
Continuation entries are used when more than one line is required
to define a single Bulk Data entry.
Free field format example
Parent entry:
PBAR,10,20,1.25,,,,,,+PB10
The plus sign is optional. The first
column of field 10 on the parent entry is
ignored by MSC/NASTRAN.
Continuation entry:
+PB10,2.0,5.0,-2.0,-5.0
The plus sign is required in column 1,
field 1 of a continuation entry.
4-13
4-14
GENERATION/REPLICATION
To avoid the time-consuming input of each Bulk Data entry
individually, repetitive fields can be generated from a single entry
definition. Rules governing this capability are:
Duplication of a field from the preceding entry is accomplished by
coding the symbol = in the associated field.
Duplication of all remaining fields from the preceding entry is
accomplished by coding the symbol == in the first of the fields to
be repeated.
Generation of a incremeted value from the previous entry is
defined by coding *X or *(X) where X is the real or integer value of
the increment.
Note: Parentheses are optional.
Repeated replication is indicated by coding =n or the optional =(n)
in field 1, where n is the number of entry images to be generated
using the values of the increments on the preceding generation
entry.
Generation/replication rules apply to all Bulk Data entries unless
denoted otherwise on specific entry definition pages in
Section 12.0 of the MSC/NASTRAN Handbook for Linear Analysis.
4-15
GENERATION/REPLICATION OF
CONTINUATION ENTRIES
Continuation fields (fields 1 and 10) may be replicated using the
following conventions:
Only letters of the alphabet and integers may be used. They are
coded into a base 36 number. That is, the sequence of numbers is
0, 1, 2,...8, 9, A, B,...Z.
The first character in the field 1 or 10 is not incremented.
MSC/NASTRAN increments continuation fields by +1. User input
values for these fields are ignored.
The number of characters in an incremented field is not increased.
For example, if the field in the first entry is 0, the field in the 37th
entry is also 0 resulting in an illegal duplicate entry. A method to
solve this problem is to start a first entry with 00. This will provide
36 squared unique fields.
See Section 3.5.1 of the MSC/NASTRAN Handbook for Linear
Analysis for examples of continuation entries in small field and
large field formats.
4-16
GENERATION/REPLICATION EXAMPLES
I N P U T
B U L K
D A T A
4-17
.
1 ..
2 ..
3 ..
4 ..
5 ..
6 ..
7
.
GRID,1,,1.,0.,0.
=,*(1),,*(1.),== $GRID GENERATION CARD
=(3)
$GRID REPLICATION CARD
CBEAM,100,1,1,2,0.,1.,0.
=,*10,=,*1,*1,== $CBEAM GENERATION CARD
=2
$CBEAM REPLICATION CARD
PBEAM,1,1,24.,72.,32.
,3.,-2.,3.,2.,-3.,2.,-3.,-2. $PBEAM CONTINUATION CARD
,YESA,0.2 $PBEAM CONTINUATION CARD
,=,*0.2
$GENERATION OF PREVIOUS CARD
=3
$REPLICATION OF PREVIOUS CARD
EIGR,10,MGIV,,30.
,MASS
$EIGR CONTINUATION CARD
$
ENDDATA
INPUT BULK DATA CARD COUNT =
15
D E C K
..
E C H O
..
..
10
4-18
CARD
COUNT
123456789101112131415161718-
.
1 ..
2
CBEAM
100
CBEAM
110
CBEAM
120
CBEAM
130
EIGR
10
++000008MASS
GRID
1
GRID
2
GRID
3
GRID
4
GRID
5
PBEAM
1
++0000013.
++000002YESA
++000003YESA
++000004YESA
++000005YESA
++000006YESA
ENDDATA
TOTAL COUNT=
19
..
3
1
1
1
1
MGIV
1
-2.
0.2
.4
.6
.8
1.0
..
1
2
3
4
1.
2.
3.
4.
5.
24.
3.
B U L K
..
5
2
3
4
5
30.
0.
0.
0.
0.
0.
72.
2.
..
0.
0.
0.
0.
D A T A
..
1.
1.
1.
1.
E C H O
..
0.
0.
0.
0.
..
..
10
+000008
+000009
0.
0.
0.
0.
0.
32.
-3.
2.
-3.
-2.
+000001
+000002
+000003
+000004
+000005
+000006
+000007
4-19
4-20
4-21
4-22
SECTION 5
EXECUTIVE CONTROL SECTION
INTRODUCTION TO THE EXECUTIVE CONTROL SECTION . . . . . . . . . . . . . . 5-1
SELECTED EXECUTIVE CONTROL STATEMENTS . . . . . . . . . . . . . . . . . . . . . . 5-2
DEFINITION OF DMAP . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-3
MSC/NASTRAN SOLUTION SEQUENCE FOR STATIC ANALYSIS . . . . . . . . . . 5-4
A ROAD-MAP FROM MSC/NASTRAN. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-5
F04 OUTPUT FILE FROM SOL 24 STATIC ANALYSIS . . . . . . . . . . . . . . . . . . . . 5-6
EPSILON (RESIDUAL ERROR VECTOR) AND STRAIN ENERGY . . . . . . . . . . 5-7
STRUCTURED SOLUTION SEQUENCES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-9
RIGID FORMAT SOLUTION SEQUENCES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-10
5-1
ID A1,A2
Optional
Second Entry:
SOL K
Last Entry:
CEND
Required entry
Optional Entries
TIME K
Diagnostic:
DIAG K
Some
5-2
DEFINITION OF DMAP
MSC/NASTRANs Executive System uses an internal, data block
oriented language called Direct Matrix Abstraction Programming
(DMAP) to direct MSC/NASTRAN to solve problems.
DMAP
Performs the operations of converting input lists to matrices and/or
tables
Performs the matrix solutions
Converts matrix solutions to output lists
Prints the solution (and/or any intermediate information)
These operations are selected by requesting one of the solution
sequences (SOL K). Each SOL assembles hundreds to thousands
of DMAP commands to perform a specific type of analysis.
User access to DMAP enables:
Altering (modifying) a built-in solution sequence flow
Writing a solution sequence to meet special analysis needs
Note:
5-3
Modules
Preface
MSGMESH
Geometry Processors
GP1,GP2,GP3,GP4
Plot Preparation
(Undeformed)
PLOT
Matrix Assembly
EMG,EMA,GPSP1
Multipoint Constraint
Elimination
MCE1,MCE2
Single-Point Constraint
Elimination
SCE1
Stiffness
Matrix Decomposition
DCMP
Static Solution
SSG1,SSG2,SSG3
Data Recovery
SDR1,SDR2,GPFDR
Plot Preparation
(Deformed)
PLOT
5-4
DAY
DEL_MB
CPU
SEC
DEL_CPU
16:56:39
0:37
2.9
0.0
8.9
SEPREP2
17
GP1
BEGN
16:56:40
0:38
2.9
0.0
9.5
SEPREP2
17
GP1
END
Module Name
DMAP Sequence ID
subDMAP
Elasped Time for Job
(used for "time" limit)
File Operations
Wall Clock - Elasped Seconds
Time of Day
If anything should go wrong during a run, this output is invaluable
as a trouble-shooting aid.
5-5
5-6
DAY TIME
ELAPSED
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:14
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
17:40:15
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:03
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
0:04
I/O SEC
23.4
23.5
23.9
23.9
23.9
24.0
24.0
24.0
24.0
24.0
24.0
24.0
24.1
24.1
24.1
24.1
24.1
24.1
24.1
24.1
24.1
24.1
24.1
24.1
24.1
24.1
24.1
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
24.2
DEL_I/O
2.6
0.2
0.3
0.0
0.0
0.1
0.0
0.0
0.0
0.0
0.0
0.0
0.1
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.1
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
CPU SEC
1.2
1.2
1.6
1.6
1.7
1.7
1.8
1.8
1.9
1.9
1.9
1.9
2.0
2.0
2.0
2.0
2.0
2.0
2.1
2.1
2.1
2.1
2.1
2.1
2.1
2.1
2.1
2.1
2.2
2.2
2.2
2.2
2.2
2.2
2.3
2.3
2.4
2.4
2.4
2.4
2.4
2.4
2.4
2.4
2.5
2.5
2.5
2.5
DEL_CPU
0.1
0.1
0.4
0.0
0.1
0.0
0.1
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.1
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
MODULE/SUBROUTINE MESSAGES
STATICS
STATICS
STATICS
IFP
STATICS
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
SOL24
4
5
6
IFP1
XSORT
IFP
12 SOL24
BEGN
9
GP1
BEGN
11 GP2
BEGN
26 GP3
BEGN
28 TA1
BEGN
31 EMG
BEGN
35 XPURGE BEGN
37 EMA
BEGN
39 XPURGE BEGN
43 ELTPRT BEGN
73 XEQUIV BEGN
85 XEQUIV BEGN
89 XEQUIV BEGN
97 GP4
BEGN
99 GPSP
BEGN
100 XPURGE BEGN
104 XEQUIV BEGN
108 XEQUIV BEGN
110 SCE1
BEGN
112 XEQUIV BEGN
113 XPURGE BEGN
134 XEQUIV BEGN
135 XPURGE BEGN
140 DCMP
BEGN
145 SSG1
BEGN
146 XEQUIV BEGN
148 SSG2
BEGN
150 XEQUIV BEGN
150 XEQUIV BEGN
155 SSG3
BEGN
160 SDR1
BEGN
176 SDR2
BEGN
185 OFP
BEGN
186 SDRX
BEGN
187 XEQUIV BEGN
187 XEQUIV BEGN
188 OFP
BEGN
189 GPFDR
BEGN
190 OFP
BEGN
193 OFP
BEGN
210 SDR2
BEGN
211 XEQUIV BEGN
215 OFP
BEGN
220 EXIT
BEGN
{ P } = { P } [ K ] { u }
The residual load vector exists due to numerical roundoff in the
calculations. If there were no numerical roundoff,
{ P } = 0
It also calculates the residual error vector
5-7
10-6
5-8
Full Name
Abbreviation
101
Statics
SESTATIC
103
Normal modes
SEMODES
105
Buckling
SEBUCKL
106
Nonlinear Statics
NLSTATIC
107
SEDCEIG
108
SEDFREQ
109
SEDTRAN
110
SEMCEIG
111
SEMFREQ
112
SEMTRAN
114
CYCSTATX
115
CYCMODE
116
CYCFREQ
129
NLTRAN
144
AESTAT
145
Aerodynamic Flutter
SEFLUTTER
146
Aeroelastic Response
SEAERO
153
NLHEAT
159
TRHEAT
200
Design Optimization
DESOPT
5-9
Rigid Format
Number
(SOL K)
Full Name
Abbreviation
Static analysis
STATICS1
Normal Modes
Modes
Geometric Nonlinear
GNOLIN
Buckling
BUCKLING
DCEIG
DFREQ
DTRAN
10
MCEIG
11
MFREQ
12
MTRAN
5-10
SECTION 6
CASE CONTROL SECTION
INTRODUCTION TO THE CASE CONTROL SECTION. . . . . . . . . . . . . . . . . . . . 6-1
CASE CONTROL COMMAND - DATA SELECTION . . . . . . . . . . . . . . . . . . . . . . 6-2
CASE CONTROL COMMANDS - OUTPUT SELECTION . . . . . . . . . . . . . . . . . . 6-4
CASE CONTROL SET DEFINITIONS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6-9
CASE CONTROL SORT OPTIONS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6-11
SORT1 FORMAT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6-12
SORT2 FORMAT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6-14
CASE CONTROL COMMANDS - SUBCASE DEFINITION . . . . . . . . . . . . . . . . . 6-17
CASE CONTROL EXAMPLES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6-19
6-1
DEFORM
Example:
LOAD = 10
DEFORM = 20
specifies that the set of initial element deformations defined by all
Bulk Data DEFORM entries with set ID number 20 and all Bulk Data
static load entries having set ID number 10 define the set of grid point
loads for this static analysis
6-2
MPC
6-3
SUBTITLE
LABEL
UNSORT
BOTH
NONE
PUNCH
Example:
ECHO = BOTH
causes both the alphabetically sorted Bulk Data and the Bulk Data in
the order input by the user to be printed in the output file.
6-4
ELSTRESS or STRESS
STRAIN
ESE
DISPLACEMENT(PLOT)
SPCFORCES
OLOAD
GPFORCE
6-6
contraint
Case Control
CEND
TITLE = Use GPFORCE Request
TEMP(LOAD) = 100
SPC = 200
LOAD = 120
DISP = ALL
FORCE = ALL
STRESS = ALL
GPFORCE = ALL
BEGIN BULK
Applied loads
SPC forces
Element elastic forces
6-7
Note:
The applied loads and the element forces include the selfequilibrating forces that arise from the initial thermal strain
(TEMP(LOAD)).
6-8
Example:
CEND
TITLE = OUTPUT SELECTION EXAMPLE
SUBTITLE = ILLUSTRATES USE OF SETS
LOAD = 15
SET 1 = 3,4,7,9,11
SET 5 = 2,9,15 THRU 21,33
DISP = 1 $ DISPLACEMENTS FOR GRIDS 3,4,7,9,11
FORCE = 1 $ FORCES FOR ELEMENTS 3,4,7,9,11
STRESS = 5
OLOAD = ALL
BEGIN BULK
For this example, the following items are computed and printed to
the output file
The displacement for Grid Points 3,4,7,9, and 11
The force in Elements 3,4,7,9, and 11
The stresses in Elements 2,9,15 THRU 21,33
6-9
6-10
SORT2
A request for SORT2 format with any output request results in all
output requests being printed in SORT2 format.
Warning:
6-11
SORT1 FORMAT
EXAMPLE OF SORT1 OUTPUT
SORT1 IS THE DEFAULT OUTPUT FORMAT
AUGUST
16, 1994
D I S P L A C E M E N T
POINT ID.
TYPE
1
2
3
4
5
G
G
G
G
G
T1
0.0
2.758621E-02
5.517241E-02
8.275862E-02
1.103448E-01
T3
6-12
TYPE
1
2
3
4
5
T1
G
G
G
G
G
0.0
0.0
0.0
0.0
0.0
16, 1994
MSC/NASTRAN
0.0
0.0
0.0
0.0
0.0
5/31/94
R3
0.0
0.0
0.0
0.0
0.0
PAGE
11
SUBCASE 2
V E C T O R
T2
0.0
0.0
0.0
0.0
0.0
10
R2
0.0
0.0
0.0
0.0
0.0
D I S P L A C E M E N T
POINT ID.
PAGE
SUBCASE 1
R1
0.0
0.0
0.0
0.0
0.0
AUGUST
5/31/94
V E C T O R
T2
0.0
0.0
0.0
0.0
0.0
MSC/NASTRAN
T3
0.0
0.0
0.0
0.0
0.0
R1
0.0
4.563709E-02
9.127419E-02
1.369113E-01
1.825484E-01
R2
0.0
0.0
0.0
0.0
0.0
R3
0.0
0.0
0.0
0.0
0.0
AUGUST
F O R C E S
16, 1994
ELEMENT
ID.
1
3
AXIAL
FORCE
2.000000E+05
2.000000E+05
TORQUE
0.0
0.0
6-13
AUGUST
F O R C E S
ELEMENT
ID.
1
3
AXIAL
FORCE
0.0
0.0
I N
TORQUE
1.000000E+05
1.000000E+05
5/31/94
PAGE
12
SUBCASE 1
R O D
E L E M E N T S
( C R O D )
ELEMENT
AXIAL
ID.
FORCE
TORQUE
2
2.000000E+05
0.0
4
2.000000E+05
0.0
16, 1994
I N
MSC/NASTRAN
MSC/NASTRAN
5/31/94
PAGE
13
SUBCASE 2
R O D
E L E M E N T S
ELEMENT
AXIAL
ID.
FORCE
2
0.0
4
0.0
( C R O D )
TORQUE
1.000000E+05
1.000000E+05
SORT2 FORMAT
EXAMPLE OF SORT2 OUTPUT
AUGUST 17, 1994
SORT2 SELECTION IS DISP(SORT2)=ALL, FORCE(SORT2)=ALL
POINT-ID =
1
D I S P L A C E M E N T
SUBCASE
1
2
TYPE
G
G
T1
0.0
0.0
T2
0.0
0.0
6-14
TYPE
G
G
T1
2.758621E-02
0.0
T2
0.0
0.0
TYPE
G
G
T1
5.517241E-02
0.0
T3
0.0
0.0
T2
0.0
0.0
T3
0.0
0.0
PAGE
10
V E C T O R
R1
R2
0.0
0.0
5/31/94
T3
0.0
0.0
MSC/NASTRAN
MSC/NASTRAN
R3
0.0
0.0
5/31/94
0.0
0.0
PAGE
11
V E C T O R
R1
0.0
4.563709E-02
MSC/NASTRAN
R2
R3
0.0
0.0
5/31/94
0.0
0.0
PAGE
12
V E C T O R
R1
0.0
9.127419E-02
R2
0.0
0.0
R3
0.0
0.0
TYPE
G
G
T1
8.275862E-02
0.0
T2
0.0
0.0
T3
0.0
0.0
R1
0.0
1.369113E-01
6-15
SUBCASE
1
2
TYPE
G
G
T1
1.103448E-01
0.0
T2
0.0
0.0
T3
0.0
0.0
5/31/94
PAGE
13
R2
R3
0.0
0.0
5/31/94
0.0
0.0
R1
0.0
1.825484E-01
PAGE
14
R2
R3
0.0
0.0
0.0
0.0
SUBCASE
1
AXIAL
FORCE
2.000000E+05
TORQUE
0.0
SUBCASE
2
AXIAL
FORCE
0.0
15
TORQUE
1.000000E+05
SUBCASE
1
AXIAL
FORCE
2.000000E+05
TORQUE
0.0
5/31/94
PAGE
16
( C R O D )
AXIAL
FORCE
SUBCASE
2
TORQUE
1.000000E+05
0.0
6-16
SUBCASE
1
AXIAL
FORCE
2.000000E+05
TORQUE
0.0
SUBCASE
1
AXIAL
FORCE
2.000000E+05
TORQUE
0.0
AXIAL
FORCE
SUBCASE
2
17, 1994
R O D
0.0
MSC/NASTRAN
5/31/94
E L E M E N T S
SUBCASE
2
17
TORQUE
1.000000E+05
PAGE
18
( C R O D )
AXIAL
FORCE
0.0
TORQUE
1.000000E+05
SUBCOM
SUBSEQ
Note:
6-17
SYMCOM
SYMSEQ
MODES
6-18
Case Control
CEND
TITLE = Equipment Rack Analysis
ECHO = BOTH
SPC = 20
SET 1 = 1 THRU 50
DISP = 1
SUBCASE 1
SUBTITLE = Dead Load
LOAD = 10
DISP = ALL
SUBCASE 2
SUBTITLE = NW Wind Load
LOAD = 20
SET 10 = 2,4,6
DISP = 10
SUBCASE 3
SUBTITLE = SW Wind Load
LOAD = 30
SUBCOM 10
SUBTITLE = Load Combination 1
LABEL = Dead Load + NW Wind
SUBSEQ = 1.0,1.0,0.0
SUBCOM 20
SUBTITLE = Load Combination 2
LABEL = Dead Load + (-) 1.5 SW Wind
SUBSEQ = 1.0,0.0,-1.5
STRESS = ALL
BEGIN BULK
6-19
Case Control
CEND
TITLE = Analysis of Half Tank
SPCFORCE = ALL
SET 1 = 1 THRU 100
SET 2 = 3,17,26,31
DISP = 1
ELFORCE = 2
SYM 1
SUBTITLE = Symmetry Boundary,Load
LOAD = 40
SPC = 10
OLOAD = ALL
SYM 2
SUBTITLE = Anti-Symmetry Boundary,Load
SPC = 15
LOAD = 50
SYMCOM 3
SYMSEQ 1., 1.
SYMCOM 4
SYMSEQ 1., -1.
BEGIN BULK
6-20
Must be defined
since element force
requests for sets 2
and 3 cannot be made
within Subcase 1
Note:
Case Control
CEND
TITLE = Analysis of Large Model
SUBTITLE = Use REPCASE to Get
Output
SPC = 22
LOAD = 900
SET 1 = 1 THRU 10
SET 2 = 20 THRU 50
SET 3 = 500 THRU 1100
SUBCASE 1
LABEL = Nodal Data for Whole Model
DISP = ALL
SPCF = ALL
ELFO = 1
REPCASE 2
LABEL = Forces for Set 2
ELFO = 2
REPCASE 3
LABEL = Forces for Set 3
ELFO = 3
BEGIN BULK
6-21
Case Control
CEND
TITLE = CRANKSHAFT ANALYSIS
SET 1 = 1000 THRU 1200
SET 2 = 2000 THRU 2400
i=1
6-22
SECTION 7
BULK DATA SECTION
INTRODUCTION TO THE BULK DATA SECTION . . . . . . . . . . . . . . . . . . . . . . . . 7.1-1
MSC/NASTRAN COORDINATE SYSTEMS . . . . . . . . . . . . . . . . . . . . . . . . . . 7.1-2
GRID POINTS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.2-1
SINGLE-POINT CONSTRAINTS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.3-1
MATERIAL PROPERTIES. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.4-1
MSC/NASTRAN ELEMENTS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.5-1
ZERO-DIMENSIONAL (SCALAR) ELEMENT LIBRARY . . . . . . . . . . . . . . . . . . . 7.6-1
ONE-DIMENSIONAL ELEMENT LIBRARY . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.7-1
CALCULATION OF TORSIONAL CONSTANT J
FOR PROD, PBAR, AND PBEAM ENTRIES . . . . . . . . . . . . . . . . . . . . . . . . . 7.7-6
CROD EXAMPLE. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.7-8
ONE-DIMENSIONAL ELEMENT LIBRARY . . . . . . . . . . . . . . . . . . . . . . . . . . 7.7-12
CBAR ELEMENT DEFINITION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.7-14
CBARAO . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.7-37
CBEAM ELEMENT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.7-44
TWO-DIMENSIONAL ELEMENT LIBRARY . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.8-1
TWO-DIMENSIONAL ELEMENTS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.8-2
REFERENCES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.8-3
CQUAD4 ELEMENT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.8-4
CQUAD4 ELEMENT DEFINITION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.8-6
QUAD4 ELEMENT DEFINITION - ALTERNATE PROPERTY ENTRY . . . . . . 7.8-18
SECTION 7 (Cont.)
BULK DATA SECTION
THREE-DIMENSIONAL ELEMENT LIBRARY . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.9-1
RECOMMENDED USES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.9-2
CHEXA . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.9-3
MSC/NASTRAN LOADS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.10-1
FORCE AND MOMENT ENTRIES. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.10-2
EXAMPLE INPUT FOR FORCE AND
MOMENT ENTRIES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.10-3
DISTRIBUTED LOADS (PLOAD) IN MSC/NASTRAN . . . . . . . . . . . . . . . . . . 7.10-5
EXAMPLE INPUT FOR PLOAD1 ENTRY . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.10-6
EXAMPLE INPUT FOR PLOAD4 ENTRY . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.10-10
GRAV ENTRY . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.10-12
THERMAL EFFECTS. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.10-14
EXAMPLE OF APPLYING THERMAL EFFECTS . . . . . . . . . . . . . . . . . . . . . . 7.10-15
LOAD ENTRY. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.10-18
THE LOAD REQUEST TREE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.10-19
The Bulk Data Section contains all data necessary for describing a
structural model
Bulk Data definitions include:
Geometry
User-definable coordinate systems
Geometric locations of grid points
Constraints
Material properties
Element connections
Element properties
Loads
The Bulk Data Section is not required to be input in any particular
order. It is sorted automatically at the beginning of the analysis.
7.1-1
Rectangular
CORD1R
CORD2R
Cylindrical
CORD1C
CORD2C
Spherical
CORD1S
CORD2S
7.1-2
7.1-3
u2
B
u1
y
Note:
U1 =
x direction
U2 =
y direction
U3 =
z direction
7.1-4
z
B
U
P
Ur
C
z
A
x
Point A
Point P
Point C
Point B
(Ur, U, Uz)
7.1-5
B
P
U
y
Point A
Point P
Point C
Point B
(Ur, U, U)
Note:
cylindrical spherical
7.1-6
50.0
z1
zB
xB
R1
100.
7.1-7
GRID POINTS
Grid points are used to specify:
Structural geometry
Degrees of freedom of the structure
Locations of points at which displacements are constrained or
loads are applied
Locations where output quantities are to be calculated
Each GRID entry refers to 2 coordinate systems. One for locating
the grid point and the other for establishing the grid point
displacement coordinate system that defines for the given grid
point the directions of the nodal displacements, degrees of
freedom, constraints, and solution vectors.
NOTE: constraints, MPCs, R-type elements, and output all use the
displacement coordinate system.
All matrices used inside
MSC/NASTRAN use the displacement coordinate system.
If you specify the displacement coordinate system as a cylindrical
or spherical system, the program creates a right-handed cartesian
coordinate system at the grid point, which is the projection of the
coordinate system used. This used as the displacement system for
that point.
7.2-1
5
2
7.2-2
GRID
ID
CP
X1
X2
X3
CD
PS
SEID
GRID
101
5.0
10.0
2.0
123
Field
10
Contents
ID
CP
X1, X2, X3
CD
PS
SEID
Superelement ID
7.2-3
1
20
YB
XB
7.2-4
10
R
1
20
YB
XB
7.2-5
1
10
7.2-6
SINGLE-POINT CONSTRAINTS
A single-point constraint (SPC) is a constraint applied to one or
more components of motion at selected grid or scalar points. Uses
of SPCs include:
Support a structure
Apply symmetric or antisymmetric boundary conditions by
restraining the DOFs that must have zero values in order to satisfy
symmetry or antisymmetry
Remove degrees of freedom unconnected or weakly coupled to the
structure
Remove degrees of freedom not used in the structural analysis
(e.g., out-of-plane DOFs for a 2-D analysis)
Apply zero or nonzero enforced displacements to grid points
SPCs can be applied as:
Permanent constraints - defined on GRID entry
User-requested in Case Control with SPC=SID.
Bulk Data entry on SPC, SPC1, or SPCD entries
Defined in the
Automatic - PARAM,AUTOSPC,YES
Reaction forces at SPCd grids (termed single-point forces of
constraint), may be obtained by including the Case Control request
SPCFORCES=ALL.
7.3-1
GRID
ID
CP
X1
X2
X3
CD
PS
SEID
GRID
1.0
-2.0
3.0
10
245
7.3-2
SPC
SID
SPC
15
1501
123
1502
456
0.0
Enforced Displacement
(Default = 0.0)
Components
Grid Point ID
Set ID, Selected in Case Control as SPC = 15
7.3-3
10
SPC1
SID
G1
G2
G3
G4
G5
G6
SPC1
10
G7
68
G8
G9
-etc.-
+BC
Alternate Form
SPC1
SID
GID1
THRU
GID2
SPC1
313
12456
THRU
32
7.3-4
10
ABC
SPC Set
SPC1
SID
G1
G2
G3
G4
G5
SPC1
13
14
SPCD
SID
SPCD
11
13
.05
14
.07
G6
7.3-5
7.3-6
MATERIAL PROPERTIES
Definitions
Stress-Strain Curve
(Typical Structural Steel)
Stress (psi) or
(N/mm2)
Linear Region
E
1
Strain
Linear
Youngs Modulus E
Elastic
Poisson Ratio
7.4-1
7.4-2
types
of
materials
can
be
defined
with
Isotropic MAT1
Two-dimensional anisotropic MAT2
Axisymmetric solid orthotropic MAT3
Two-dimensional orthotropic MAT8
Three-dimensional anisotropic MAT9
Temperature-dependent material properties are defined on MATTi
entries.
Note:
MAT2 through MAT9 and MATTi entries are beyond the scope
of this class.
For further information, refer to the
MSC/NASTRAN Reference Manual, Vol. I, Chapter 6.
7.4-3
E
G = --------------------2(1 + )
For line and surface elements
E
Tref
7.4-4
MAT1
MID
NU
RHO
TREF
GE
MAT1
2.9+7
6.5E-6
50.0
0.33
Field
10
+M1
Contents
MID
NU
RHO
TREF
GE
ST
+M1
SC
SS
MCSID
36000.
ST,SC,SS
MCSID
7.4-5
7.4-6
MSC/NASTRAN ELEMENTS
Element Library
Over 50 finite elements
One-dimensional
Two-dimensional
Three-dimensional
Scalar
Axisymmetric
Rigid
Mass and damping
Heat transfer
Genel user-supplied element
Fluid-structure
p-version
Contact
7.5-1
Surface
Elements
Solid
Elements
CROD
CQUAD4
CHEXA
CELASi (i=1,2,3,4)
CONROD
CQUAD8
CPENTA
GENEL
CTUBE
CTRIA3
CTETRA
CONM2
CBAR
CTRIA6
CTRIAX6
CBEAM
CSHEAR
CBEND
7.5-2
Other
Elements
B
C
[ ke ]
A
=
7.5-3
[ ke ]
7.5-4
7.5-5
7.5-6
B
X
The CELASi elements are connected by two degrees of freedom one at each grid/ground connection point
The CBUSH elements connects from 1 to 6 dof between two GRID
points.
Force components: axial force P
or moment M
Displacement components:
axial translation u
or rotation
7.6-1
The CBUSH is the recommended form for scalar springs. It avoids the
potential grounding which may occur when two non-coincident points
are connected. The CELASi elements simply insert terms directly into
the stiffness matrix without considering geometry or displacement
coordinate systems. The CBUSH correctly accounts for the effects of
geometry and displacement coordinate systems.
See page 61 and pp. 121 through 125 of the MSC/NASTRAN Linear
Static Analysis Users Guide and Section 5.6 of the MSC/NASTRAN
Reference Manual for detailed information about scalar elements. The
CBUSH element is documented in the V69 Release Guide.
7.6-2
1202
P
1201
K = 100. lb/in
P = 10. lb
S = 50.
100.
CELAS2
EID
G1
C1
G2
C2
GE
CELAS2
1200
100.
1201
1202
GRID
1201
0.
0.
0.
123456
GRID
1202
100.
0.
0.
23456
50.
Field
Contents
EID
G1,G2
C1,C2
GE
Notes:
10
1.
2.
7.6-3
POINT ID.
1201
1202
TYPE
G
G
T1
0.0
1.000000E-01
T2
V E C T O R
T3
0.0
0.0
R1
0.0
0.0
R2
0.0
0.0
R3
0.0
0.0
0.0
0.0
P
10.
= ---- = ----------- = 0.1 in
K
100.
F O R C E S
7.6-4
ELEMENT
ID.
1200
FORCE
I N
ELEMENT
ID.
S C A L A R
FORCE
S P R I N G S
ELEMENT
ID.
( C E L A S 2 )
FORCE
ELEMENT
ID.
FORCE
-1.000000E+01
STRESS
ELEMENT
ID.
I N
S C A L A R
STRESS
S P R I N G S
ELEMENT
ID.
STRESS
-5.000000E+02
( C E L A S 2 )
ELEMENT
ID.
STRESS
D I S P L A C E M E N T
CBUSH
EID
PID
GA
GB
GO/X1
CBUSH
OCID
S1
S2
S3
7
X2
X3
CID
10
39
100
75
39
39
100
Field
Contents
EID
PID
GA, GB
Xi
GO
CID
7.6-5
Remarks:
CID >_ 0 overrides GO and Xi. Then the element x-axis is along T1, the
element y-axis is along T2, and the element z-axis is along
T3 of the CID coordinate system. If the CID refers to a
cylindrical coordinate system of a shperical coordinate
system, then grid GA is used to locate the system. If for
cylindrical or spherical coordinate, GA falls on the z-axis
used to define them, it is recommended that another CID
be selected to define the element x-axis.
For noncoincident grids (GA =/ GB), when GO or (X1, X2, X3) is given
and no CID is specified, the line AB is the element x-axis
and the orientation vector v lies in the x-y plane (similar to
the CBEAM element).
For noncoincident grids (GA =/ GB), if neither GO or (X1, X2, X3) is
specified and no CID is specified, then the line AB is the
element x-axis. This option is valid only when K1 (or B1)
or K4 (or B4) or both on the PBUSH entry are specified
(but K2, K3, K5, K6 or B2, B3, B5, B6 are not specified). If
K2, K3, K5, or K6 (or B2, B3, B5, or B6) are specified, a
fatal message will be issued.
7.6-6
GA
S*
yelem
(1 - S)*
GB
xelem
zelem
GA
S*
yelem
(1 - S ) *
GB
xelem
Figure 1. CBUSH Element.
7.6-7
GB
GA
Note: 1. The material stiffness and damping properties of the
elastomer are located at (S1, S2, S3).
zelem
yelem
GB
GA
Note: 1. The material stiffness and damping properties of the elastomer are located at (S1, S2, S3).
2. The elastomer itself has zero length; i.e., GA and GB are coincident. It is shown here in an exploded view
7.6-8
PBUSH
PID
"K"
K1
K2
K3
K4
K5
K6
"B"
B1
B2
B3
B4
B5
B6
"GE"
GE1
"RCV"
SA
ST
EA
ET
10
35
4.35
GE
.06
RCV
7.3
2.4
3.1
.03
3.3
35
2.3
Field
Contents
PID
"K"
Ki
"B"
Flag indicating that the next 1 to 6 fields are force-pervelcoity damping. (Character)
Bi
"GE"
GE1
7.6-9
SA
ST
EA
ET
Remarks:
Ki, Bi, or GE1 may be made frequency dependent for both direct and
modal frequency response by use of PBUSHT entry.
The "K", "B", "GE", or "RCV" entries may be specified in any order.
7.6-10
BAR:
BEAM:
BEND:
7.7-1
axial force P
torque T
Displacement components:
ui
i
CONROD
CTUBE
7.7-2
A
T
Xe
The rod element x axis (Xe) is defined along the line connecting End A
to End B.
Note:
7.7-3
CROD
EID
PID
G1
G2
CROD
101
Field
10
Contents
EID
PID
G1,G2
7.7-4
PROD
PID
MID
NSM
PROD
101
201
0.10
Field
10
Contents
PID
MID
Area of rod
1.
2.
7.7-5
1 4
J = --- r
2
2r
J = 2.25 a
2a
7.7-6
J = ab
2b
4
b
b
------ 3.36 --- 1 -------------
4
a
3
12a
3 16
2a
7.7-7
CROD EXAMPLE
y
102
100
101
P
100.
Element Properties
A = 5. in
r = 1.262 in.
4
4
r
J = --------- = 3.984 in
2
Material Properties
E = 29. E + 6 psi
G = 11. E + 6 psi
y = 36000. psi
Applied Loads
P = 2. E + 6 psi
T = 1. E + 5 in lbs
7.7-8
10
CROD
100
101
102
GRID
101
0.
0.
0.
123456
GRID
102
100.
0.
0.
2356
PROD
201
5.
3.984
1.262
MAT1
201
2.9+7
11.+6
+M1
36000.
Note:
+M1
36000.
7.7-9
D I S P L A C E M E N T
POINT ID.
101
102
TYPE
G
G
T1
0.0
1.379310E-01
T2
T3
0.0
0.0
0.0
0.0
R1
0.0
2.281855E-01
R2
0.0
0.0
R3
0.0
0.0
TL
( 1.E + 5 ) ( 100. )
1 = -------- = ----------------------------------------------- = 2.281855E 01
JG ( 3.987 ) ( 11.E + 6 )
7.7-10
F O R C E S
ELEMENT
ID.
100
AXIAL
FORCE
2.000000E+05
S T R E S S E S
I N
R O D
E L E M E N T S
ELEMENT
ID.
TORQUE
1.000000E+05
I N
T H E
( C R O D )
R O D
E L E M E N T S
AXIAL
FORCE
TORQUE
( C R O D )
PL ( 2.E + 5 ) ( 100. )
= -------- = ---------------------------------------- = 1.37931E 01
AE
( 5. ) ( 29.E + 6 )
V E C T O R
ELEMENT
AXIAL
ID.
STRESS
100 4.000000E+04
SAFETY
MARGIN
-1.0E-01
TORSIONAL
STRESS
3.167671E+04
P 2.E + 5
= ---- = ------------------ = 4.0E + 4
A
5.
SAFETY
MARGIN
1.4E-01
ELEMENT
ID.
AXIAL
STRESS
SAFETY
MARGIN
TORSIONAL
STRESS
TC ( 1.E + 5 ) ( 1.262 )
= -------- = -------------------------------------------- = 3.167671E + 4
J
3.984
Stress allow
M.S. = ---------------------------------- 1.0
Stress allow
7.7-11
36000.
M.S. axial = ------------------ 1.0 = 0.1
40000.
36000.
M.S. torsion = ------------------ 1.0 = 0.14
31680.
SAFETY
MARGIN
7.7-12
7.7-13
CBAR
EID
PID
GA
GB
X1,G0
X2
X3
CBAR
102
1.0
1.0
PA
PB
W1A
W2A
W3A
W1B
W2B
0.5
-1.2
+B1
0.5
10
+B1
W3B
0.6
Continuation Entry
Field
Contents
EID
PID
GA,GB
X1,X2,X3
G0
PA,PB
W1A,W2A,W3A
W1B,W2B,W3B
7.7-14
z
End b
Plane 1
ye
xe
y
x
Grid Point b
Plane 2
End a
Offset W a
Offset W b
ze
y
x
Grid Point a
7.7-15
Element Z-axis:
Cross-product Xex V .
Xe Ye plane.
Note:
7.7-16
V V V
V
GO
V
V
V
7.7-17
1
2
7.7-18
Ye
V ( X 1, X 2, X 3 )
GO
WA
WB
GA
GB
7.7-19
Xe GA when V is defined
using GO.
7.7-20
PBAR
PID
MID
I1
I2
NSM
PBAR
102
202
0.30
.0090
.0063
0.0124
C1
C2
D1
D2
E1
E2
F1
F2
0.3
-0.25
0.3
0.25
-0.3
0.25
-0.3
-0.25
K1
K2
I12
+P1
10
+P1
+P2
+P2
Field
Contents
PID
MID
I1,I2,I12
Torsional constant
Ci,Di,Ei,Fi
K1,K2
7.7-21
= beam length
= shear factor
and the term 1/KAG is called the shear flexibility of the beam.
K accounts for the shear distribution over the element cross section,
and its value depends on the geometry of the cross section.
7.7-22
Value of K
Rectangular
5/6
Circular
9/10
1/2
Af 2A
Major Axis
Aw A
References:
1.
Roark and Young, Formulas for Stress and Strain, 5th ed., p.
185.
7.7-23
BH
I 1 = I zz = -----------12
ze
HB
I 2 = I yy = -----------12
C
D
E
F
= .5H, .5B
= .5H, .5B
= .5H, .5B
= .5H, .5B
ze
HB
I 1 = I zz = -----------12
BH
I 2 = I yy = -----------12
C
D
E
F
H
ye
= .5B, .5H
= .5B, .5H
= .5B, .5H
= .5B, .5H
7.7-24
Format:
1
PBARL
PID
MID
DIM1
DIM2
DIM3
DIM9
-etc.-
NSM
10
GROUP TYPE
DIM4
DIM5
DIM6
DIM7
DIM8
Example:
PBARL
Field
39
14..
6.
.5
.5
.5
Contents
PID
MID
DIMi
NSM
7.7-25
yelem
C
DIM 1
DIM 1
zelem
yelem
TYPE = "ROD"
zelem
DIM 2
TYPE = "TUBE"
DIM 4
DIM 3
DIM 2
zelem
E
D
DIM 1
yelem
TYPE = "L"
(PBEAML only)
DIM 3
yelem
DIM 4
C
F
C
zelem
DIM 6
DIM 2
zelem
DIM 4
DIM 3
DIM 1
DIM 5
E
D
DIM 2
D
DIM 1
TYPE = "I"
TYPE = "CHAN"
7.7-26
yelem
DIM 3
DIM 1
zelem
DIM 2
DIM 3
zelem
DIM 2
D
DIM 4
E
yelem
TYPE = "T"
F
DIM 4
DIM 1
C
TYPE = "BOX"
zelem
DIM 2
D
yelem
DIM 1
2
yelem
E
DIM 1
DIM 1
2
DIM 2
TYPE = "BAR"
C
DIM 2
2
C
DIM 3
DIM 3
DIM 4
D
zelem
zelem
DIM 4
E
D
DIM 1
DIM 2
TYPE = "H"
TYPE = "CROSS"
7.7-27
DIM 1
DIM 2
2
DIM 1
2
C
DIM 2
DIM 1
IM 3
zelem
DIM 4
DIM 3
D
TYPE = "T1"
yelem
DIM 2
TYPE = "I1"
DIM 1
DIM 3
zelem
DIM 4
E
yelem
DIM 1
F
TYPE = "CHAN1"
DIM 1
DIM 2
DIM 1
C
yelem
DIM 3
zelem
DIM 3
DIM 4
DIM 2
E
TYPE = Z
DIM 1
TYPE = "CHAN2"
7.7-28
C
zelem
PBARL
yelem
yelem
DIM 1
DIM 4
C
DIM 3
zelem
DIM 3
DIM 2
DIM 2
zelem
E
DIM 1
DIM 4
E
TYPE = T2
DIM 6
DIM 5
TYPE = BOX1
yelem
DIM 4
DIM 4
DIM 3
E
F
DIM 3
DIM 1
zelem
DIM 2
DIM 1
E
DIM 2
TYPE = HAT
TYPE = HEXA
7.7-29
7.7-30
Fx
Xe
ye
ve
M1a
M1b
xe
Fx
a
Plane 1
Fx
v1
ze
ve
M2a
M2b
xe
a
Plane 2
v2
7.7-31
yb
P2
P2 = 5000. lb
T
= 4.0E + 4 in-lb
3400
3401
3402
P1
100.
Element Properties
A = 24. in2
ye
I1 = I = 72. in4
zz
I2 = I = 32. in4
yy
6.
J
= 75.12
ze
in4
7.7-32
4.
xB
= 30.E + 6 psi
= 0.3
= 36000 psi
E
= --------------------- = 11.54E + 6 psi
2(1 + )
10
CBAR
3400
3401
3402
0.
1.
0.
GRID
3401
0.
0.
0.
GRID
3402
100.
0.
0.
PBAR
10
24.
72.
32.
75.12
+PB1
3.
-2.
3.
2.
-3.
2.
MAT1
10
30.+6
0.3
+M1
+M1
36000.
BAR
+PB1
123456
+PB1
-3.
-2.
or
PBAR
10
+PB1
4.
6.
7.7-33
D I S P L A C E M E N T
POINT ID.
3401
3402
TYPE
G
G
T1
0.0
3.333333E-03
T2
0.0
-7.716050E-01
T3
0.0
0.0
V E C T O R
R1
0.0
4.614223E-03
R2
0.0
0.0
7.7-34
P1 L
( 2.4E + 4 ) ( 100. )
1 = ----------- = ------------------------------------------- = 3.33E-3 in.
AE
( .24 ) ( 30.E + 6 )
P2 L
( 5000. ) ( 100. )
2 = -------------- = ------------------------------------------------- = 7.716E-1 in.
3EI 1
( 3 ) ( 30.E + 6 ) ( 72. )
R3
0.0
-1.157407E-02
I N
B A R
E L E M E N T S
30, 1992
( C B A R )
ELEMENT
BEND-MOMENT END-A
BEND-MOMENT END-B
- SHEAR AXIAL
ID.
PLANE 1
PLANE 2
PLANE 1
PLANE 2
PLANE 1
PLANE 2
FORCE
3400 -5.000000E+05 0.0
0.0
0.0
-5.000000E+03 0.0
2.400000E+04
P2 = 5000.
7.7-35
Shear
V = 5000.
V=P
P2 = 5000.
0
Moment
M1 = 5.0E+5
M1 = Px = 5000(100)
M2 = 0
PAGE 6
TORQUE
4.000000E+04
F O R C E S
MARCH
ELEMENT
SA1
ID.
SB1
3400 2.083333E+04
0.0
I N
B A R
30, 1992
E L E M E N T S
SA2
SA3
SA4
AXIAL
SA-MAX
SA-MIN
SB2
SB3
SB4
STRESS
SB-MAX
SB-MIN
2.083333E+04 -2.083333E+04 -2.083333E+04 1.000000E+03 2.183333E+04 -1.983333E+04
0.0
0.0
0.0
1.000000E+03 1.000000E+03
Normal Stress
Due to Axial
Load Only
Combined Axial
and Bending
Stress
My
( 5.0E + 5 ) ( 3. )
bending = -------- = ------------------------------------- = 2.0833E + 4
I1
72.
P
2.4E + 4
axial = ---- = ---------------------- = 1.0E + 3
A
24.
S allow
36000.
M.S. ( T ) = ----------------------- 1.0 = ------------------ 1.0 = 0.65
S actual
21833.
S allow
36000.
M.S. ( C ) = ----------------------- 1.0 = ------------------ 1.0 = 0.82
S actual
19833.
PAGE 7
Optional
M.S.-T
M.S.-C
6.5E-01
8.2E-01
S T R E S S E S
MARCH
CBARAO
The CBARAO Bulk Data entry defines a series of locations along a
CBAR element at which stresses and/or internal element forces can
be recovered. For CBEAM elements, the PBEAM entry incorporates
the capabilities of the CBARAO.
The CBARAO entry is normally used only to evaluate the effects of
applying a nonuniformly distributed load (PLOAD1) to a CBAR
element.
7.7-37
CBARAO (Cont.)
CBARAO Entry Format
1
CBARAO
EID
SCALE
X1
X2
X3
X4
X5
X6
CBARAO
100
FR
.2
.4
.6
.8
Field
10
Contents
EID
SCALE
Scale of X1 values
LE = Xi are actual distances along the element length.
FR = Xi are ratios of actual distances to the bar length.
Xi
7.7-38
CBARAO (Cont.)
CBAR/CBARAO Example
Applied Loading
yB
P2
P1 = 24000. lb.
P2 = 1000. lb.
3402
3401
P1
3400
50.
100.
Element Properties
ye
A = 24 in2
C
ze
in4
7.7-39
4.
xB
CBARAO (Cont.)
Material Properties
E
30.E+6 psi
0.3
CBAR
3400
3401
3402
0.
1.
0.
CBARAO
3400
LE
25.
50.
75.
GRID
3401
0.
0.
0.
12345
GRID
3402
100.
0.
0.
2345
PBAR
24.
72.
32.
75.12
+PB1
3.
-2.
3.
2.
-3.
2.
7.7-40
10
+PB1
-3.
-2.
PAGE 5
D I S P L A C E M E N T
TYPE
G
G
T1
0.0
-3.333333E-03
T2
T3
0.0
0.0
R1
0.0
0.0
0.0
0.0
MARCH 30, 1992
R2
0.0
0.0
R3
-2.893519E-04
2.893519E-04
ELEMENT
ID.
3400
3400
3400
STATION
(PCT)
0.000
0.250
0.500
D I S T R I B U T I O N
BEND-MOMENT
PLANE 1
PLANE 2
0.0
1.250000E+04
2.500000E+04
0.0
0.0
0.0
I N
B A R
E L E M E N T S
SHEAR FORCE
PLANE 1
PLANE 2
-5.000000E+02
-5.000000E+02
-5.000000E+02
0.0
0.0
0.0
( C B A R )
AXIAL
FORCE
TORQUE
-2.400000E+04
-2.400000E+04
-2.400000E+04
0.0
0.0
0.0
PAGE 6
CBARAO (Cont.)
7.7-41
POINT ID.
3401
3402
V E C T O R
3400
3400
0.750
1.000
1.250000E+04
0.0
0.0
0.0
5.000000E+02
5.000000E+02
0.0
0.0
-2.400000E+04
-2.400000E+04
0.0
0.0
V = +500.
CBARAO (Cont.)
0
Shear
-500.
7.7-42
M = +25000.
M = +12500.
Moment
PAGE 8
SXC
D I S T R I B U T I O N
SXD
0.0
0.0
-5.208334E+02 -5.208334E+02
-1.041667E+03 -1.041667E+03
I N
SXE
0.0
5.208334E+02
1.041667E+03
B A R
E L E M E N T S
SXF
0.0
5.208334E+02
1.041667E+03
( C B A R )
AXIAL
-1.000000E+03
-1.000000E+03
-1.000000E+03
S-MAX
S-MIN
-1.000000E+03 -1.000000E+03
-4.791666E+02 -1.520833E+03
4.166675E+01 -2.041667E+03
M.S.
3.5E+01
2.3E+01
1.7E+01
3400
3400
0.750
1.000
-5.208334E+02 -5.208334E+02
0.0
0.0
5.208334E+02
0.0
5.208334E+02
0.0
-1.000000E+03
-1.000000E+03
-4.791666E+02 -1.520833E+03
-1.000000E+03 -1.000000E+03
My
( 25000. ) ( 3. )
bending = -------- = --------------------------------- = 1041.7
I
72.
P
24000.
axial = ---- = ------------------ = 1000.
A
24.
7.7-43
S allow
36000.
M.S. = ----------------------- 1.0 = ------------------ 1.0 = 16.63
S actual
1041.7
where S actual = max { S max , S min }
2.3E+01
3.5E+01
CBEAM ELEMENT
Connected to two grid points
Force components
Axial force P
Total torque T
Warping torque Tw
Bending moments in Planes 1 and 2 Mi
Shears in Planes 1 and 2 Vi
Displacement component
ui
i
(d/dx)i (represented by SPOINTs)
7.7-44
7.7-45
CBEAM
EID
PID
GA
GB
X1,G0
X2
X3
CBEAM
39
13
PA
PB
W1A
W2A
W3A
+23
+34
Field
513
SA
SB
10
123
W1B
W2B
W3B
3.0
234
Contents
EID
PID
GA,GB
X1,X2,X3
G0
PA,PB
7.7-46
Contents
W1A,W2A,W3A
W1B,W2B,W3B
SA,SB
7.7-47
M1(A), M2(A)
zelem
N1A
Nonstructural Mass
Center of Gravity
N2A
yelem
Plane 2
Plane 1
(0, 0, 0)
V
Neutral Axis
xelem
w a Offset
Shear Center
N2B
Grid Point GA
N1B
(xb, 0, 0)
w b Offset
7.7-48
PBEAM
PID
MID
A(A)
I1(A)
I2(A)
I12(A)
J(A)
NSM(A)
PBEAM
39
2.9
3.5
5.97
C1(A)
C2(A)
D1(A)
D2(A)
E1(A)
2.0
-4.0
+23
+34
X/XB
I1
I2
YES
1.0
5.3
56.2
78.6
C1
C2
D1
D2
E1
2.5
-5.0
S1
S2
K1
K2
+56
E2(A)
I12
M2(A)
M1(B)
F1(A)
F2(A)
NSM
345
E2
F1
F2
+56
1.1
M1(A)
0.21
M2(B)
N1(A)
+67
Field
123
+34
S0
+45
10
N2(A)
0.5
Contents
N1(B)
567
N2(B)
0.0
Default
Values
PID
Required
MID
Required
A(A)
Required
I1(A)
7.7-49
Contents
Default
Values
I2(A)
Required
I12(A)
0.0
( I 1 I 2 I 12 > 0 )
J(A)
0.0
NSM(A)
0.0
Ci(A),Di(A),
Ei(A), Fi(A)
yi = zi = 0.0
S0
Required
YES:
X/XB
7.7-50
Required
Contents
Default
Values
Ci,Di,Ei,Fi
The
y,z
locations
in
element
coordinates for the cross section
located at X/XB. The values are fiber
locations for stress data recovery.
K1,K2
1.0, 1.0
S1,S2
0.0, 0.0
0.0, same as
end A
0.0, same as
end A
7.7-51
Field
Contents
Default
Values
7.7-52
MB
hB
hA
MA
P
Tension
Tension
No Shear
Shear
Compression
Compression
The value of the shear coefficient for a tapered beam with heavy
flanges that sustain the entire moment load may then be written as
2 ( hA hB )
S 1 = ----------------------------( hA + hB )
See the MSC/NASTRAN Reference Manual, Section 5.2.1.
7.7-53
d
G ------- J
dx
where E
2
2
d
d
d
------ E ---------- C w ---------- = m
dx
2
2
dx
dx
Cw = warping constant
Note:
= shear modulus
= torsion constant
The
The warping constant Cw has units of (length)6.
development of the above differential equation and methods
for the numerical evaluations of the warping constant are
available in the literature. (See, for example, Timoshenko and
Gere, Theory of Elastic Stability, McGraw Hill Book Company,
1961. Also see Roark & Young, Formulas for Stress and
Strain, for values for different sections.)
7.7-54
_
>
7.7-55
PID
MID
DIM1(A) DIM2(A)
10
GROUP TYPE
-etc.-
DIMn(A)
NSM(A)
SO(1)
X(1)XB
DIM1(1)
DIM1(2)
DIM2(2)
DIM2(1)
-etc.-
DIMn(1)
NSM(1)
SO(2)
X(2)XB
-etc.-
DIMn(2)
-etc.-
NSM(m)
SO(m)
X(m)XB DIM1(m)
DIMn(m) NSM(m)
SO(B)
1.0
DIM1(B) DIM2(B)
-etc.-
-etc.-
DIMn(B)
NO
0.4
6.
0.6
6.
7.8
Example:
PBEAML
99
21
12.
14.8
2.5
7.
1.2
2.6
5.6
2.3
2.6
YES
YES
Field
PID
MID
Group
TYPE
Contents
Property identification number. (Integer > 0)
Material identification number. (Integer > 0)
Cross-section group. (Character; Default="MSCBMLO"
Cross-section shape. See Remark 4.. (Character: "ROD",
"TUBE", "L", "I", "CHAN", "T", "BOX", "BAR", "CROSS", "H",
"T1", "I1", "CHAN1", "Z", CHAN2", "T2", "BOX1", "HEX",
"HAT" for GROUP="MSCBMLO")
DIMi(A), DIMi(B)
Cross-section dimensions at end A and B. (Real > 0.0 for
GROUP="MSCBMLO"
NSM(A), NSM(B)
Nonstructural mass per unit length. (Real _>0.0; Default = 0.0)
7.7-56
2.
3.
4.
5.
7.7-57
6.
7.
8.
9.
7.7-58
yelem
Plane 1
Plane 2
M1
V2
M2
Neutral Axis
V1
Fx
xelem
Shear Axis
Tx
7.7-59
v1
M1a
M1b
Fx
xe
a
Plane 1
v1
v2
ze
M2a
M2b
xe
a
Plane 2
v2
For further information on the BEAM element, see Section 1.3.2 of the
MSC/NASTRAN Users Manual, Section 2.6 of the MSC/NASTRAN
Application Manual and Section 2.3.5.2 of the MSC/NASTRAN Handbook for
Linear Analysis.
7.7-60
Applied Loading
ye
P1 = 2.4E+4
P2
200
P2 = 5000. lb.
T = 4.0E+4
202
xe
201
P1
100.
Element Properties
ye
A = 24.
C
I1 = Izz = 72.
I2 = Iyy = 32.
J = 75.12
6.
ze
4.
7.7-61
30.E+6
0.3
36000.
E
--------------------- = 11.54E+6
2(1 + )
CBEAM
200
201
202
0.
1.
0.
GRID
201
0.
0.
0.
GRID
202
100.
0.
0.
PBEAM
24.
72.
32.
+PB1A
3.
-2.
3.
2.
-3.
+PB1B
YESA
1.0
MAT1
30.+6
+M1
36000.
0.3
10
123456
75.12
2.
-3.
+PB1A
-2.
+PB1B
+M1
or
PBEAMl
BAR
+PB1A
4.
6.
1.
+PB1A
YES
7.7-62
MARCH 30,1992
POINT ID.
201
202
TYPE
G
G
T1
0.0
-3.333333E-03
T2
0.0
-7.734105E-01
T3
V E C T O R
R1
0.0
4.614838E-03
0.0
0.0
7.7-63
F O R C E S
ELEMENT-ID
200
GRID
201
202
STAT DIST/
LENGTH
0.000
1.000
I N
R2
E L E M E N T S
- WEB
PLANE 1
SHEARS PLANE 2
-5.000000E+03 0.0
-5.000000E+03 0.0
R3
0.0
-1.157407E-02
0.0
0.0
B E A M
PAGE 4
PAGE 6
( C B E A M )
AXIAL
FORCE
TOTAL
TORQUE
WARPING
TORQUE
-2.400000E+04
4.000000E+04
0.0
-2.400000E+04
4.000000E+04
0.0
D I S P L A C E M E N T
S T R E S S E S
ELEMENT-ID
200
GRID
201
202
STAT DIST/
LENGTH
0.000
1.000
SXC
SXD
I N
B E A M
E L E M E N T S
SXE
SXF
( C B E A M )
S-MAX
S-MIN
7.7-64
PAGE 7
M.S.-T
M.S.-C
8.2E-01
6.5E-01
Maximum and
Minimum of Combined
Stresses at Points
C,D,E, and F
Applied Loading
P2
P1 = 2.4E+4
P2 = 1000.
202
201
P1
200
50.
100.
ye
Element Properties
C
A = 24.
I1 = Izz = 72.
I2 = Iyy = 32.
ze
6.
J = 75.12
F
7.7-65
4.
E = 30.E+6 psi
y = 36000. psi
CBEAM
200
201
202
0.
1.
0.
GRID
201
0.
0.
GRID
202
100.
0.
0.
PBEAM
24.
72.
32.
+PB1A
3.
-2.
3.
2.
-3.
+PB1B
YESA
.25
+PB1C
+PB1C
YESA
.5
+PB1D
+PB1D
YESA
.75
+PB1E
+PB1E
YESA
1.
MAT1
30.+6
+M1
36000.
0.3
75.12
10
2.
-3.
+PB1A
-2.
+PB1B
+M1
or
PBEAML
BAR
+PB1A
4.
6.
YES
+PB1B
YES
.5
BAR
+PB1C
1.0
+BP1D
+PB1C
YES
7.7-66
+PB1A
.25
+PB1B
D I S P L A C E M E N T
TYPE
G
G
T1
0.0
-3.333333E-03
T2
V E C T O R
T3
0.0
0.0
R1
0.0
0.0
7.7-67
ELEMENT-ID
200
GRID
201
0
0
0
202
STAT DIST/
LENGTH
0.000
0.250
0.500
0.750
1.000
I N
7.812500E-03 0.0
1.250001E+04 0.0
2.500000E+04 0.0
1.250000E+04 0.0
0.0
0.0
E L E M E N T S
- WEB
PLANE 1
SHEARS PLANE 2
-4.999999E+02 0.0
-4.999999E+02 0.0
-4.999999E+02 0.0
5.000001E+02 0.0
5.000001E+02 0.0
R3
-2.835957E-04
2.835955E-04
0.0
0.0
B E A M
R2
0.0
0.0
F O R C E S
PAGE 5
PAGE 7
( C B E A M )
AXIAL
FORCE
TOTAL
TORQUE
-2.400000E+04
0.0
-2.400000E+04
0.0
-2.400000E+04
0.0
-2.400000E+04
0.0
-2.400000E+04
0.0
WARPING
TORQUE
0.0
0.0
0.0
0.0
0.0
POINT ID.
201
202
S T R E S S E S
ELEMENT-ID
200
GRID
7.7-68
201
0
0
0
202
STAT DIST/
LENGTH
0.000
0.250
0.500
0.750
1.000
SXC
SXD
I N
B E A M
E L E M E N T S
SXE
SXF
PAGE 8
( C B E A M )
S-MAX
S-MIN
M.S.-T
M.S.-C
8.6E+02
1.7E+01
Three-noded
isoparametric
flat
plate
element.
Commonly used for mesh
transitions. May have excessive stiffness
particularly for membrane strain.
QUAD4
SHEAR
TRIA6
QUAD8
TRIAR
QUADR
7.8-1
TWO-DIMENSIONAL ELEMENTS
Plates - Background
Definition:
7.8-2
REFERENCES
References on basic plate theory:
1.
Theory of Plates and Shells, by S. Timoshenko and S. Woinowsky-Krieger, 2nd ed., McGraw Hill, 1959
2.
7.8-3
CQUAD4 ELEMENT
Element most commonly used for representing plates, shells, and
membranes.
ye
Lateral, transverse, normal
xe
Membrane, in-plane
Fxy
Mx,My
Mxy
Vx,Vy
ui
x, y (no rotation normal to element)
7.8-4
Ze
Fx
Vy
Ye
Fy
Mxy
My
Fyx Vx
Mxy
Fyx
Mxy
Fxy
Fyx Mx
My
Fx
Xe
Mxy
Fy
Vy
Forces and moments: Calculated at element centroid by default for
CQUAD4 and CTRIA3 (may be calculated at
corners if desired)
Calculated at centroid and grid points for
CQUAD8, CTRIA6, CQUADR, and CTRIAR
Ye y
xy
xy
x
x
xy
Xe
xy
y
Stresses:
7.8-5
CQUAD4
EID
PID
G1
G2
G3
G4
ZOFFS
CQUAD4
111
203
31
74
75
32
2.6
T1
T2
T3
T4
1.77
2.04
2.09
1.80
+BC
Field
10
ABC
Contents
EID
PID
G1,G2,
G3,G4
T1,T2,
T3,T4
ZOFFS
7.8-6
G4
xelement
xmaterial
G1
G2
7.8-7
7.8-8
PSHELL
PID
MID1
MID2
12/T3
MID3
TS/T
NSM
PSHELL
204
.025
204
Z1
Z2
MID4
204
10
+PS1
+PS1
Field
Contents
PID
MID1
MID2
12I/T3
MID3
7.8-9
Contents
TS/T
NSM
Z1,Z2
MID4
The QUAD4 element can have in-plane, bending, and transverse shear
behavior. The element mechanical behavior is specified by the
presence or absence of a material ID number in the appropriate
field(s) on the PSHELL entry.
7.8-10
PID
MID1
PSHELL
204
.025
MID2
12I/T3
MID3
TS/T
MID3
TS/T
PID
PSHELL
MID1
MID2
.025
204
12I/T3
PID
PSHELL
Note:
MID1
MID2
.025
204
12I/T3
7.8-11
MID3
204
TS/T
t =1
20 50
MID3
t
Note:
MID4:
MID2
MID3
MID4
No bending-membrane coupling
7.8-12
3303
P1
10.0
3300
3302 P2
3301
P1
10.0
P2
0.1
P1
E = 30.E+6
P 1 = 12000.
= 0.3
P 2 = 5000.
CQUAD4
3300
3301
3302
3303
3304
GRID
3301
0.
0.
0.
123456
GRID
3302
10.
0.
0.
GRID
3303
10.
10.
0.
GRID
3304
0.
10.
0.
123456
0.1
PSHELL
MAT1
30.+6
0.3
7.8-13
10
D I S P L A C E M E N T
7.8-14
POINT ID.
3301
3302
3303
3304
TYPE
G
G
G
G
T1
0.0
7.806971E-03
7.806971E-03
0.0
T2
0.0
1.756568E-03
-1.756568E-03
0.0
T3
0.0
-1.261228E+02
-1.261228E+02
0.0
V E C T O R
R1
0.0
-3.192983E+00
3.192983E+00
0.0
R2
0.0
1.915789E+01
1.915789E+01
0.0
R3
0.0
0.0
0.0
0.0
F O R C E S
ELEMENT
ID
3300
I N
Q U A D R I L A T E R A L
E L E M E N T S
( Q U A D 4 )
7.8-15
- MEMBRANE FORCES FX
FY
FXY
2.400000E+03 1.930295E+02 7.105427E-15
- BENDING MOMENTS MX
MY
MXY
-5.000000E+03 -7.017544E+02 1.776357E-12
Transverse
Shear Forces
(Force/Length)
P2
M = 100000.
MX = 50000.
in-lb/10 in
= 5000.
in-lb/in
ELEMENT
ID.
3300
FIBRE
DISTANCE
-5.000000E-02
5.000000E-02
I N
Q U A D R I L A T E R A L
E L E M E N T S
( Q U A D 4 )
VON MISES
2.790149E+06
2.836263E+06
P
24000.
axial = ---- = ------------------ = 2.4E + 4
A
1
7.8-16
My
( 50000. ) ( .05 )
bending = -------- = ------------------------------------ = 3.0E + 6
3
I
10(.1) 12
+ = a + b = 3.024E + 6
- = a + b = 2.976E + 6
2
HVM = ( x x y + y + 3 xy )
2
12
2 12
= 2.836E + 6
S T R E S S E S
S T R A I N S
ELEMENT
ID.
3300
FIBRE
DISTANCE
-5.000000E-02
5.000000E-02
I N
Q U A D R I L A T E R A L
E L E M E N T S
( Q U A D 4 )
7.8-17
This output was obtained with the Case Control request: STRAIN(FIBER) = ALL
VON MISES
6.920528E-02
7.037997E-02
7.8-18
PENTA
(6-15 nodes)
HEXA
(8-20 nodes)
TETRA
(4-10 nodes)
HEXA
PENTA
Any or all edge midside
nodes may be deleted. Corner
nodes may not be deleted.
TETRA
7.9-1
RECOMMENDED USES
HEXA
PENTA
TETRA
7.9-2
CHEXA
Connected to 8 to 20 grid points (8 or 20 nodes recommended for
best results)
Stress components: x, y, z, xy, yz, zx (at center and corner
points)
Displacement components: ui
Can reference either an isotropic or anisotropic material entry
7.9-3
CHEXA (Cont.)
Geometry
G18
G7
G6
G19
G14
G17
G15
G20
G2
G8
G10
G5
G3
G16
G13
G9
G11
G12
G4
G1
CHEXA
EID
PID
G1
G2
G3
G4
G5
G6
CHEXA
71
G7
G8
G9
G10
G11
G12
G13
G14
10
30
31
53
54
G15
G16
G17
G18
G19
G20
55
56
57
58
59
60
+BC
+EF
Field
10
ABC
DEF
Contents
EID
PID
G1,...,G20
7.9-4
CHEXA (Cont.)
Element Coordinate System
R vector
of
faces
G4-G1-G5-G8
and
S vector
of
faces
G1-G2-G6-G5
and
T vector
of
faces
G1-G2-G3-G4
and
XYZ axes
G7
T
G6
G8
G5
Centroid
Centroid
G3
G4
G1
7.9-5
G2
CHEXA (Cont.)
Properties
1
PSOLID
PID
MID
CORDM
IN
STRESS
ISOP
PSOLID
100
TWO
GRID
REDUCED
Field
10
FCTN
Contents
PID
MID
FCTN
7.9-6
CHEXA (Cont.)
Properties
Remarks:
1.
PSOLID entries should have unique identification numbers with respect to all other property entries.
2.
3.
See the CHEXA, CPENTA, or CTETRA entry for the defintion of the
element coordinate system.
The material coordinate system
(CORDM) may be the basic system (0 or blank), any defined
system (Integer>0), or the element coordinate system (-1). The
default value for CORDM is zero unless it is overridden by the
NASTRAN statement with the CORDM keyword. See Chaper 1 of
the MSC/NASTRAN Quick Reference Guide.
4.
5.
b.
c.
7.9-7
CHEXA (Cont.)
Properties
6.
7.
8.
9.
10. The gauss point locations for the solid elements are documented in
Section 15.3 of the MSC/NASTRAN Reference Manual.
7.9-8
CHEXA (Cont.)
Properties
CHEXA
Integration
BUBBLE or
2x2x2 Reduced
Shear with Bubble Bland or 0
(default)
Function (default)
2x2x2 Reduced
Shear Only
8 Node
2x2x2 Standard
Isoparametric
2x2x2 Reduced
Shear Only
2x2x2 Standard
Isoparametric
2x2x2 Reduced
Shear Only
9-20
NODE
3x3x3 Standard
Isoparametric
Blank or REDUCED
(Default*)
Yes
FULL or 1
Blank or REDUCED
THREE or
3
Blank or GRID
No
FULL or 1
Blank or REDUCED
TWO or 2
2x2x2 Standard
Isoparametric
3x3x3 Reduced
Shear Only
(default)
p-elements
TWO or 2
Blank or GRID
or
GAUSS or 1
Blank or GRID
0 or 1
Bubble, P+ISOP
Integration
No Bubble,
P+ISOP Integration
2 or 3
No
Blank or REDUCED
(Default*)
Blank or
THREE or
3
(Default)
Reduced (p-order)
Bubble
FULL or 1
FULL or 1
Not
applicable
No
-10 < ISOP < 10
_
_
7.9-9
CHEXA (Cont.)
Interpretation of Element Output
Stress output consists of the six components of stress measured
in the material coordinate system defined in the CORDM field of the
PSOLID entry (default=BASIC system).
Additional output includes the magnitude and direction of the three
principal stresses, the mean pressure, and the octahedral stresses.
These stresses are provided at the corner grid points and at the
center of each element.
See Section 5.4 of the MSC/NASTRAN Reference Manual and section
4.3 of the MSC/NASTRAN Linear Static Analysis Users Guide for
more detailed information about solid elements.
7.9-10
CHEXA (Cont.)
Example Input
z
6711
6714
10.
6712 6701 6713
y
6704
6702
6703
x
10.
E = 30.E+6 psi
= 0.3
p = 8. psi
7.9-11
10.
S O R T E D
.
1 ..
2
CHEXA
6700
+CH1
6713
GRID
6701
GRID
6702
GRID
6703
GRID
6704
GRID
6711
GRID
6712
GRID
6713
GRID
6714
MAT1
1
PLOAD4 1
PSOLID 1
ENDDATA
TOTAL COUNT=
..
3
1
6714
D A T A
..
5
6702
..
6
6703
0.
10.
10.
0.
0.
10.
10.
0.
0.
0.
10.
10.
0.
0.
10.
10.
.3
8.
0.
0.
0.
0.
10.
10.
10.
10.
30.E6
6700
1
..
4
6701
8.
8.
..
7
6704
E C H O
..
8
6711
..
9
6712
.. 10
+CH1
123456
23456
3456
3456
456
456
456
456
8.
6711
6713
14
EPSILON
1.3239410E-17
EXTERNAL WORK
1.0666667E-03
CHEXA (Cont.)
7.9-12
CARD
COUNT
12345678910111213-
B U L K
S T R E S S E S
ELEMENT-ID
6700
I N
CORNER
------CENTER
GRID-ID
NORMAL
-1GRID CS 8 GP
CENTER X
1.665335E-16
Y
2.775558E-17
Z -8.000000E+00
H E X A H E D R O N
S O L I D
E L E M E N T S
( H E X A )
DIR. COSINES
-A- -B- -C-
MEAN
PRESSURE
VON MISES
XY
YZ
ZX
4.163336E-17
5.551115E-17
5.551115E-17
A
B
C
-8.000000E+00
1.110223E-16
1.110223E-16
LX 0.00 0.0
LY 0.00 0.0
LZ 1.00 0.0
0.0
0.0
0.0
2.666667E+00
8.000000E+00
X
Y
Z
3.747003E-16
5.273559E-16
-8.000000E+00
XY
YZ
ZX
1.249001E-16
2.775558E-17
5.150794E-17
A
B
C
-8.000000E+00
4.996004E-16
4.996004E-16
LX 0.00 0.0
LY 0.00 0.0
LZ 1.00 0.0
0.0
0.0
0.0
2.666667E+00
8.000000E+00
6702
X
Y
Z
3.122502E-16
-1.249001E-16
-8.000000E+00
XY
YZ
ZX
1.249001E-16
8.593553E-17
5.150794E-17
A
B
C
-8.000000E+00
1.110223E-16
1.110223E-16
LX 0.00 0.0
LY 0.00 0.0
LZ 1.00 0.0
0.0
0.0
0.0
2.666667E+00
8.000000E+00
6703
X
Y
Z
-2.775558E-17
-3.538836E-16
-8.000000E+00
XY
YZ
ZX
1.249001E-16
8.593553E-17
2.775558E-17
A
B
C
-8.000000E+00
-2.220446E-16
-2.220446E-16
LX 0.00 0.0
LY 0.00 0.0
LZ 1.00 0.0
0.0
0.0
0.0
2.666667E+00
8.000000E+00
6704
X
Y
Z
-8.326673E-17
3.608225E-16
-8.000000E+00
XY
YZ
ZX
1.249001E-16
2.775558E-17
2.775558E-17
A
B
C
-8.000000E+00
1.665335E-16
1.665335E-16
LX 0.00 0.0
LY 0.00 0.0
LZ 1.00 0.0
0.0
0.0
0.0
2.666667E+00
8.000000E+00
6711
X
Y
Z
-4.163336E-17
-1.942890E-16
-8.000000E+00
XY
YZ
ZX
-4.163336E-17
2.775558E-17
5.150794E-17
A
B
C
-8.000000E+00
-1.110223E-16
-1.110223E-16
LX 0.00 0.0
LY 0.00 0.0
LZ 1.00 0.0
0.0
0.0
0.0
2.666667E+00
8.000000E+00
6712
X
Y
Z
-6.938894E-18
3.330669E-16
-8.000000E+00
XY
YZ
ZX
-4.163336E-17
8.593553E-17
5.150794E-17
A
B
C
-8.000000E+00
1.110223E-16
1.110223E-16
LX 0.00 0.0
LY 0.00 0.0
LZ 1.00 0.0
0.0
0.0
0.0
2.666667E+00
8.000000E+00
6713
X
Y
Z
3.053113E-16
3.747003E-16
-8.000000E+00
XY
YZ
ZX
-4.163336E-17
8.593553E-17
2.775558E-17
A
B
C
-8.000000E+00
3.885781E-16
3.885781E-16
LX 0.00 0.0
LY 0.00 0.0
LZ 1.00 0.0
0.0
0.0
0.0
2.666667E+00
8.000000E+00
6714
X
Y
Z
8.326673E-17
-3.330669E-16
-8.000000E+00
XY
YZ
ZX
-4.163336E-17
2.775558E-17
2.775558E-17
A
B
C
-8.000000E+00
-1.110223E-16
-1.110223E-16
LX 0.00 0.0
LY 0.00 0.0
LZ 1.00 0.0
0.0
0.0
0.0
2.666667E+00
8.000000E+00
CHEXA (Cont.)
7.9-13
6701
Example Output
CHEXA (Cont.)
Example Stress Calculations
For the output on the preceding page,
p = pressure = 8.0 psi
p =
net force on element face = -80 lb/in2 (10 in x 10 in) = -800.0 lbs
Thus,
z = principal stress = P= 8.0 psi
po
1
mean pressure = --- ( x + y + z )
3
1
--- ( 8.0 ) = 2.667 psi
3
Note:
2
2
2
2
2 12
1
------- [ ( x y ) + ( y z ) + ( z x ) + 6 yz + 6 xy ]
2
12
1
------- ( 64 + 64 )
2
= 8.000 psi
7.9-14
MSC/NASTRAN LOADS
Summary of Static Loads in MSC/NASTRAN
Type of Load
SLOAD
PLOAD1
PLOAD,
PLOADX
Gravity
Centrifugal Force
Thermal Expansion
DEFORM
SPCD, SPC
LOAD
7.10-1
PLOAD2,
PLOAD4,
7.10-2
x
2
P1
T1
FORCE
SID
CID
FORCE
10
FORCE
N1
N2
N3
5000.
0.
-1.
0.
10
24000.
1.
MOMENT
SID
N1
N2
N3
MOMENT
10
40000.
1.
CID
Field
10
Contents
SID
CID
F or M
N1,N2,N3
7.10-3
f = FN
m = MN
where N = ( N1 , N2 , N3 )
7.10-4
PLOAD1
PLOAD2
PLOAD4
PLOADX
Element
PLOAD1
BAR
BEAM
BEND
(Restricted)
QUAD4
PLOAD2
PLOAD4
QUAD8
TRIA3
TRIA6
SHEAR
PLOADX
PLOAD
HEXA
PENTA
TETRA
TRIAX6
See pages 12-115 - 12-124 of the Handbook for Linear Analysis and
pages 2.4-255 - 2.4-264 of the MSC/NASTRAN Users Manual for
detailed information about the PLOAD type entries.
7.10-5
XB
100
Length = 4.0
PLOAD1
SID
EID
PLOAD1
100
TYPE SCALE
FY
Field
FR
X1
P1
X2
P2
0.0
-15.0
1.0
-15.0
10
Contents
SID
EID
Element ID number
TYPE
P1,P2
7.10-6
YB
w2 = 20.0
w1 = 8.0
XB
100
1.0
PLOAD1
SID
EID
PLOAD1
100
2.0
1.0
TYPE SCALE
FY
LE
X1
P1
X2
P2
1.0
-8.0
3.0
-20.0
7.10-7
10
YB
P = 1000.
XB
100
2.0
PLOAD1
SID
EID
PLOAD1
100
Note:
2.0
TYPE SCALE
FY
X1
P1
X2
P2
.5
-1000.
FR
10
7.10-8
YB
144
100
XB
R1
1
PLOAD1
SID
EID
PLOAD1
100
Note:
100
5
TYPE SCALE
FY
LEPR
R2
X1
P1
X2
P2
0.0
-50.
144.
-50.
7.10-9
10
11
14
12
2
1
PLOAD
PLOAD4
13
1
4
3
SID
EID
P1
P2
P3
P4
G1
G3 or G4
CID
N1
N2
N3
10
100
10.
11
13
Field
10
Contents
SID
EID
Element ID number
P1,P2,P3,P4
7.10-10
Contents
G1
G3
Note:
7.10-11
GRAV ENTRY
The GRAV entry is used to define the direction and magnitude of a
uniform linear (gravity) acceleration vector in any defined
coordinate system.
The GRAV entry may be used to apply accelerations to a model.
The resulting load is computed by using the gravity vector and
mass matrix.
Cannot be used at scalar points
1
GRAV
SID
CID
N1
N2
N3
GRAV
14
13.0
1.0
Field
10
2.0
Contents
SID
CID
N1,N2,N3
7.10-12
2.
3.
7.10-13
THERMAL EFFECTS
Several temperature definition entries are available for including
thermal effects in an analysis. Selection of the appropriate entry is
based on the component the temperature will be applied to.
Use TEMP, TEMPD for grid points.
Use TEMPRB for ROD, BAR, BEAM, BEND, CONROD, TUBE.
Use TEMPP1, 3 for 2-D plates.
For the thermal effects defined on any of these temperature entries
to be included in the analysis, the user must define the reference
temperature (TREF) and coefficient of thermal expansion () on the
material entries. Also, the Case Control request TEMP=SID must
be included.
If thermal effects are requested, all elements must have a
temperature field defined. If only a portion of the model requires
thermal effects, the remaining portion can:
Reference a material entry having the same material properties, a
different MID, and = 0.0
Reference a material entry having the same material properties, a
different MID, and TREF = the applied temperature (i.e., T = 0)
See pages 12-160 - 12-167 of the MSC/NASTRAN Handbook for
Linear Analysis and the MSC/NASTRAN Users Manual for detailed
information about each temperature entry.
7.10-14
E
= 1.73E9 psf
Solution:
(1)
L
Substituting = ------- = ---- into Equation 1 gives
L E
-------- = T
E
(2)
7.10-15
(3)
TEMPRB
SID
EID1
TA
TB
TEMPRB
400
100
-275.
-275.
10
MAT1
MID
NU
RHO
MAT1
1000
1.728E9
0.3
ALPHA TREF
6.56E-6
9
GE
0.
7.10-16
10
7.10-17
LOAD ENTRY
Defines a static load as a linear combination of load sets
consisting of:
Direct grid point loads
Pressures
Element dependent loads
Gravity loads (mandatory if gravity loads are used in combination
with other types of applied loads)
Selected with Case Control request LOAD=SID
1
LOAD
SID
S1
L1
S2
L2
S3
L3
LOAD
16
1.0
2.0
11
1.0
13
0.5
15
10
LOAD = [ S ( S1 L1 ) + ( S2 L2 ) + ( S3 L3 ) ] . . .
where L1, L2, L3
S1, S1, S3
7.10-18
LOAD = SID1
LOADS
FORCE
FORCE1
FORCE2
MEMENT
MOMENT1
MOMENT2
PLOAD
PLOAD2
PLOAD4
PLOADX
RFORCE
SLOAD
Bulk
Data
Entries
Notes:
GRAV
SPCD
DEFORM = SID3
Bulk
Data
Entries
SID SID1
SID = SID1
Bulk
Data
Entries
TEMP(LOAD) = SID2
SID = SID2
TEMP
TEMPD
TEMPP1
TEMPP3
TEMPRB
SID = SID3
DEFORMS
1.
2.
An SPCD entry cannot be the only requested load. Thus, if the user desires
to enforce displacements only, he must also request a legally filled out FORCE,
MOMENT, etc., entry with zero or small magnitude and the same SID as the
SPCD entry.
7.10-19
7.10-20
SECTION 8
INTERMEDIATE MODELING PRACTICES
SET NOTATIONS
Displacement Sets
Consider all grid point and scalar point degrees of freedom in a finite
element model as the members of a single displacement set. This
global set is called the g-set; the displacement set is known as ug.
The mathematical expression for the static equilibrium of a finite
element model is
[ K gg ] { u g } = { P g }
[Kgg] must be nonsingular in order to solve this equation. To achieve
a nonsingular stiffness matrix, the user can specify the independent
subsets of {ug} to be partitioned out during matrix reduction.
For example:
um
us
8.1-1
MULTIPOINT CONSTRAINTS
A multipoint constraint (MPC) is a user-imposed linear equation that
relates displacement degrees of freedom.
MPCs are useful to
Define the relative motion between two or more grid points as a
degree of freedom
Join dissimilar elements; for example, to join elements with
rotational degrees of freedom to elements which have only
translational degrees of freedom (e.g., to join shell elements to
solid elements)
Distribute loads to several points in a structure
Model rigid connections between grid points
8.1-2
Aj uj = 0
j
where A
= constraint coefficient
1
MPC
+MPC1
SID
+MPC1
+MPC2
Field
10
Contents
SID
Constraint coefficient
8.1-3
u101
SPOINT
ID
ID
ID
ID
ID
ID
ID
ID
SPOINT
103
MPC
SID
MPC
103
1.
101
-1.
102
1.
8.1-4
10
Description
RROD
RBAR
RTRPLT
RBE1
RBE2
RBE3
RSPLINE
RSSCON
8.1-5
SINGULARITIES
A singularity is caused by a degree of freedom lacking any structural
stiffness. Some examples of singularities are
Plate Elements - Rotation normal to plate
Nonsupported
Weakly Supported
3
1
8.1-6
PARAM, AUTOSPC
When PARAM,AUTOSPC,YES is included, the following table is
generated by MSC/NASTRAN:
G R I D
POINT
TYPE
ID
P O I N T
S I N G U L A R I T Y
FAILED
STIFFNESS
DIRECTION
RATIO
OLD USET
EXCLUSIVE
T A B L E
NEW USET
UNION
EXCLUSIVE
UNION
0.00E+00
SB
SB
0.00E+00
SB
SB
0.00E+00
SB
SB
0.00E+00
SB
SB
Default stiffness
ratio = 1.0E-8
USET membership is
changed only if
PARAM,AUTOSPC,YES
is included in Bulk Data
Section
8.1-7
8.1-8
STRUCTURAL SYMMETRY
Reflective symmetry can often be employed in the modeling process
to reduce the cost of the analysis.
z, z
y, y
x, x
z, z
y, y
x, x
Reflected
Coordinate
System
Note:
Right-hand
Coordinate
System
Kright = Krefl
Mright = Mrefl
etc.
Symmetric boundary constraints for
coordinate system shown:
ux = 0 , y = 0 , z = 0
Antisymmetric boundary constraints
for coordinate system shown:
uy = 0 , uz = 0 , x = 0
8.2-1
4
A
8.2-2
2500 lb
8.2-3
P
2
8.2-4
P
2
8.2-5
A
+
SYM 1
SYM 2
SYMCOM 3
SYM 1
SYM 2
8.2-6
SYMCOM 4
8.2-7
MESH TRANSITIONS
Mesh transitions are most frequently needed to transition between a
fine mesh in areas of stress concentrations to a coarser mesh away
from these areas.
All mesh transitions sacrifice some degree of accuracy in the model.
Some examples of mesh transitions are:
Belt of Triangular Elements
8.2-8
MESHTRANSITIONS(Cont.)
Higher-Order Elements with Deleted Mid-Side Nodes
Q4
Q4
Q8
Q8
Q4
Q4
Q4
Q4
Q4
Q8
Q4
Q8
Q4
Q4
Q4
Spline Interpolation
RSPLINE
8.2-9
Aspect ratio
b
a
< 4:1 *
b
a
Aspect ratio should be less than about 4:1 (much less in regions
where stress levels change rapidly). In cases of nearly-uniaxial
stress fields, larger aspect ratios are acceptable.
Skew
Angle < 20 - 30 *
Taper (2 directions)
8.2-10
h
a
Warp
8.2-11
15 arc or less
8.2-12
MODEL DEBUGGING
Structural plots are useful primarily to visually verify model geometry.
Other tools must be employed to assess the numerical accuracy of a
finite element model. These tools include automatic error checks
performed by MSC/NASTRAN and user-supplied diagnostic requests
in the form of DIAG, PARAM,DMAP Alters, and Case Control requests.
MSC/NASTRAN performs numerous error checks during an analysis
to ensure that all input data is in a proper and usable format. If an
error is detected during data processing, an error message is
generated. If the error is fatal, the analysis terminates. Following is
an example of a fatal error message and ways to determine the cause
of error.
13
3
100 lb
12
11
Y
X
Analyze this frame for the 100 lb force. Determine displacements at
all grid points.
8.3-1
N A S T R A N
E X E C U T I V E
C O N T R O L
D E C K
E C H O
ID ERROR,EX
SOL 101
CEND
*** USER WARNING MESSAGE 4519, A TIME CARD IS MISSING. DEFAULT TIME LIMIT IS ONE MINUTE.
C A S E
C O N T R O L
D E C K
E C H O
8.3-2
CARD
COUNT
1
2
3
4
5
6
S O R T E D
CARD
COUNT
1-
B U L K
.
1
CBAR
..
1
D A T A
..
10
E C H O
..
1
..
2
..
0.
..
1.
..
0.
..
..
10
8.3-3
23456789101112131415-
CBAR
CBAR
CBAR
CBAR
CBAR
FORCE
GRID
GRID
GRID
GRID
GRID
GRID
MAT1
PBAR
ENDDATA
2
5
6
11
12
200
1
2
3
11
12
13
10
10
TOTAL COUNT=
***
***
***
***
***
***
***
***
USER
USER
USER
USER
USER
USER
USER
USER
FATAL
FATAL
FATAL
FATAL
FATAL
FATAL
FATAL
FATAL
MESSAGE
MESSAGE
MESSAGE
MESSAGE
MESSAGE
MESSAGE
MESSAGE
MESSAGE
FATAL ERROR
* * * END OF JOB * * *
10
10
10
10
10
3
2
2
3
11
12
3
12
13
12
13
100.
0.
10.
20.
0.
10.
20.
0.3
0.
0.
0.
10.
10.
10.
10.5E6
10
0.
0.
0.
0.
0.
1.
0.
0.
0.
0.
0.
0.
1.
1.
1.
1.
1.
0.
0.
0.
0.
0.
0.
0.
123456
123456
.25
16
2026,ELEMENT
2026,ELEMENT
2026,ELEMENT
2026,ELEMENT
2026,ELEMENT
2026,ELEMENT
2026,ELEMENT
2026,ELEMENT
1
1
2
2
11
11
12
12
GEOMETRY
GEOMETRY
GEOMETRY
GEOMETRY
GEOMETRY
GEOMETRY
GEOMETRY
GEOMETRY
YIELDS
YIELDS
YIELDS
YIELDS
YIELDS
YIELDS
YIELDS
YIELDS
UNREASONABLE
UNREASONABLE
UNREASONABLE
UNREASONABLE
UNREASONABLE
UNREASONABLE
UNREASONABLE
UNREASONABLE
MATRIX
MATRIX
MATRIX
MATRIX
MATRIX
MATRIX
MATRIX
MATRIX
Item (3) listed here is the cause of this error. If this message did not
provide enough explanation, the user should also check the .F04 file
to determine at what point during data processing the analysis
terminated.
8.3-4
DAY TIME
ELAPSED
20:03:23
0:09
20:03:23
0:09
20:03:23
0:09
20:03:23
0:09
**** LINK 2****
20:03:27
0:13
20:03:27
0:13
20:03:29
0:15
**** LINK 10****
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
**** LINK 2****
20:03:29
0:15
20:03:29
0:15
20:03:29
0:15
20:03:34
0:20
20:03:34
0:20
**** LINK 10****
20:03:34
0:20
**** LINK 2****
20:03:34
0:20
**** LINK 8****
20:03:35
0:21
**** LINK 10****
I/O MB
DEL_MB
V A X
E X E C U T I O N
CPU SEC
DEL_CPU
S U M M A R Y
SUB_DMAP/DMAP_MODULE MESSAGES
9.3
9.4
9.4
9.4
0.3
0.0
0.0
0.0
4.8
4.9
4.9
4.9
0.2
0.1
0.0
0.1
SESTATIC
SUPER1
IFPL
IFPL
16
16
15
22
SUPER1
IFPL
IFP1
XSORT
BEGN
BEGN
BEGN
BEGN
9.4
9.4
9.4
0.0
0.0
0.0
7.0
7.0
8.6
2.0
0.0
1.6
IFPL
IFP
IFPL
24
IFP
BEGN
25
MODEPT
BEGN
9.4
9.4
9.6
9.6
9.6
9.7
9.8
9.8
9.8
9.8
10.0
10.1
10.2
10.3
10.3
0.0
0.0
0.1
0.0
0.0
0.1
0.1
0.0
0.0
0.0
0.1
0.1
0.1
0.1
0.0
8.6
8.7
8.8
8.8
8.8
8.9
9.0
9.0
9.0
9.0
9.1
9.2
9.3
9.4
9.4
0.0
0.1
0.1
0.0
0.0
0.1
0.1
0.0
0.0
0.0
0.1
0.1
0.1
0.1
0.0
IFPL
IFPL
IFPL
IFPL
IFPL
IFPL
IFPL
IFPL
IFPL
IFPL
IFPL
IFPL
IFPL
IFPL
IFPL
26
32
33
34
35
36
37
38
39
41
42
43
45
46
58
PVT
XEQUIV
XEQUIV
XEQUIV
XEQUIV
XEQUIV
XEQUIV
XEQUIV
XEQUIV
XEQUIV
XEQUIV
XEQUIV
XEQUIV
XEQUIV
PVT
BEGN
BEGN
BEGN
BEGN
BEGN
BEGN
BEGN
BEGN
BEGN
BEGN
BEGN
BEGN
BEGN
BEGN
BEGN
10.3
10.3
10.3
10.5
10.5
0.0
0.0
0.0
0.1#
0.0
9.4
9.4
9.4
12.1
12.1
0.1
0.0
0.0
2.6#
0.0
IFPL
IFPL
IFPL
IFPL
SUPER1
71
72
97
97
28
DELETE
RESTART
SEQP
SEQP
PHASE0
BEGN
BEGN
BEGN
END
BEGN
10.5
0.0
12.1
0.0
PHASE0
26
DTIIN
BEGN
10.5
0.0
12.1
0.0
PHASE0
27
GP1
BEGN
10.6
0.0
12.4
0.3
PHASE0
30
SEP1
BEGN
8.3-5
*
*
*
*
*
*
8.3-6
SECTION 9
PARAMETERS
INTRODUCTION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9-1
PARAMETERS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9-2
PRINTOUT GENERATED BY PARAM,PROUT,1 . . . . . . . . . . . . . . . . . . . . . 9-4
PRINTOUT GENERATED BY PARAM,GPECT,1 . . . . . . . . . . . . . . . . . . . . . 9-5
PRINTOUT GENERATED BY PARAM,EST,1 . . . . . . . . . . . . . . . . . . . . . . . . 9-6
INTRODUCTION
10
9-1
PARAMETERS
If PARAM entries are used, they may be placed in either the Case
Control or the Bulk Data section
Parameters may apear in the Case Control Section and Bulk Data
Section in the unstructured and structured solutions.
Some parameters frequently used are:
PARAM,AUTOSPC,YES
PARAM,GRDPNT,i
PARAM,NOGPF,1
9-2
PARAMETERS (Cont.)
PARAM,CHECKOUT,YES
PARAM,POST,-1
PARAM,WTMASS,X.X
9-3
6 HEXA
9-4
Note:
1 LAST EID =
ID
3
TYPE
HEXA
ID
4
TYPE
HEXA
ID
5
TYPE
HEXA
9-5
HEXA
HEXA
HEXA
HEXA
HEXA
HEXA
HEXA
HEXA
HEXA
HEXA
HEXA
HEXA
11
HEXA
12
HEXA
HEXA
13
HEXA
HEXA
14
HEXA
HEXA
15
HEXA
HEXA
16
HEXA
HEXA
17
HEXA
ELEMENT ID VOLUME
1
.02
ELEMENT ID VOLUME
2
.02
ELEMENT ID VOLUME
3
.02
ELEMENT ID VOLUME
4
.02
ELEMENT ID VOLUME
5
.02
*** USER WARNING MESSAGE 3103, EMGCOR OF EMG MODULE FINDS EITHER OF DATA BLOCKS 201 OR
STIF. MATRIX WILL NOT BE FORMED.
ELEMENT ID VOLUME
6
.02
9-6
SECTION 10
NATURAL FREQUENCIES AND NORMAL MODES
GOVERNING EQUATIONS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-1
MASS MATRIX . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-3
THEORETICAL RESULTS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-7
REASONS TO COMPUTE NATURAL FREQUENCIES
AND NORMAL MODES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-11
IMPORTANT FACTS AND RESULTS REGARDING
NORMAL MODES AND NATURAL FREQUENCIES . . . . . . . . . . . . . . . . . . . . . 10-12
METHODS OF COMPUTATION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-16
COMPARISON OF EIGENVALUE
EXTRACTION METHODS. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-18
NORMAL MODES ANALYSIS ENTRIES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10-20
11
GOVERNING EQUATIONS
12
x
where m
k
= mass
= stiffness
The equation of motion for free vibrations (i.e., without external load
or damping) is:
mx = kx
or
mx + kx = 0
10-1
10-2
MASS MATRIX
The mass matrix represents the inertia properties of the structure.
MSC/NASTRAN provides the user with two choices:
1.
2.
10-3
3
L
0
0
0
0
0
0
12
0
0
0
0
0
1 12
0
5 12
0
0
0
0
0
10-4
0
0
0
0
10-5
10-6
THEORETICAL RESULTS
Consider
[ M ] { x} + [ K ] { x } = 0
(1)
it
(2)
{ x} = { }e
it
(3)
[ M ] [ ]e
it
+ [ K ] [ ]e
it
=0
which simplifies to
2
([K] [M]{}) = 0
This is an eigenvalue problem.
10-7
(4)
1.
2.
det ( [ K ] [ M ] ) = 0
or
det ( [ K ] [ M ] ) = 0
where = 2
10-8
10-9
Mode 1
Mode 2
Mode 3
etc.
10-10
10-11
x1
m
k
1 = 0 { 1 } = 1
1
10-12
10-13
x1
x2
{ 1 } = 1 , { 1 } = 300
0.5
150
and
{ 1 } = .66
.33
represent the same mode of vibration.
10-14
det ( [ K ] [ M ] ) = 0
is a difficult problem. The solution to this problem must be
determined using a numerical approach.
10-15
METHODS OF COMPUTATION
MSC/NASTRAN provides the user with the following three types of
methods for eigenvalue extraction.
Tracking Methods
Eigenvalues (or natural frequencies) are determined one at a
time using an iterative technique. Two variations of the inverse
power method are provided INV and SINV. This approach is
more convenient when few natural frequencies are to be
determined. In general, SINV is more reliable than INV.
Transformation Methods
The original eigenvalue problem
([K] [M]){} = 0
is transformed to the form
[A]{} = {}
Then, the matrix [ A ] is transformed into a tridiagonal matrix
using either the Givens technique or the Householder
technique. Finally, all the eigenvalues are extracted at once
using the QR Algorithm.
Two variations of the Givens
technique and two variations of the Householder technique are
provided: GIV, MGIV, HOU, and MHOU. These methods are
more efficient for small models when a large proportion of
eigenvalues are needed.
10-16
10-17
COMPARISON OF EIGENVALUE
EXTRACTION METHODS
Method
Tridiagonal
Method
Tracking
Method with
Shifts
Obtain Eigenvalues
in Order
All at Once
Nearest to
Shift Point
Takes Advantage of
Bandwidth
No
Yes
Yes
Number of
Calculations
(Order of)
N3
NB2
NB2
Characteristic
Lanczos
Nearest
to Shift
Point
HOU,
GIV
MHOU,
MGIV
INV
SINV
Lanczos
Singular M Allowed?
No
Yes
Yes
Yes
Yes
Yes
Yes
No
Yes
Yes
Be
Careful
Be
Careful
Yes
Yes
Yes
Limitations
where N
B
Small, Dense
Systems.
Many Vectors.
Large,
Sparse
Systems.
Few Vectors
Very Large
The AGIV and AHOU methods inspect the mass matrix and
automatically select MGIV or MHOU if needed. Otherwise, GIV or HOU
is used.
10-18
COMPARISON OF EIGENVALUE
EXTRACTION METHODS (Cont.)
Number of
Eigenvalues
Desired
1000
Reduction Plus
Transformation Methods
(GDR + GIV or HOU)
or
Lanczos Method
100
GIV
MGIV
HOU
MHOU
AGIV
AHOU
10
10
INV
SINV
100
1000
10,000
10-19
Order of
Problem
Bulk Data
EIGR entry - Eigenvalue extraction entry
or
EIGRL entry for Lanczos Method
Mass properties are required.
10-20
MAT1
MID
NU
RHO
MAT1
10.+7
0.3
0.1
Nonstructural Mass
10
Mass per unit dimension (mass per unit area in this case)
PSHELL
PID
MID1
MID2
PSHELL
0.1
12I/T3
MID3
1
TS/T
NSM
0.15
Concentrated Mass
Explicit mass properties at a point (CONM2) (i.e., center of gravity
of the concentrated mass offset from the grid point, moments, and
products of inertia
10-21
10-22
F R O M
G R I D
P O I N T
REFERENCE POINT =
G E N E R A T O R
M O
0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 -3.820500E-01 *
2.547000E-01 0.000000E+00 0.000000E+00 0.000000E+00 3.820500E-01 *
0.000000E+00 2.547000E-01 3.820500E-01 -3.820500E-01 0.000000E+00 *
0.000000E+00 3.820500E-01 1.146150E+00 -5.730750E-01 0.000000E+00 *
0.000000E+00 -3.820500E-01 -5.730750E-01 1.146150E+00 0.000000E+00 *
3.820500E-01 0.000000E+00 0.000000E+00 0.000000E+00 2.292300E+00 *
S
*
*
*
1.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.000000E+00
0.000000E+00
0.000000E+00 *
0.000000E+00 *
1.000000E+00 *
DIRECTION
MASS AXIS SYSTEM (S)
10-23
X
Y
Z
MASS
2.547000E-01
2.547000E-01
2.547000E-01
X-C.G.
Y-C.G.
0.000000E+00
1.500000E+00
1.500000E+00
1.500000E+00
0.000000E+00
1.500000E+00
Z-C.G.
0.000000E+00
0.000000E+00
0.000000E+00
I(S)
* 5.730749E-01 -5.960464E-08
* -5.960464E-08 5.730749E-01
* 0.000000E+00 0.000000E+00
0.000000E+00 *
0.000000E+00 *
1.146150E+00 *
*
*
1.146150E+00 *
0.000000E+00 *
0.000000E+00 *
1.000000E+00 *
I(Q)
*
*
*
5.730749E-01
5.730749E-01
Q
*
*
*
1.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.000000E+00
0.000000E+00
1.
This is standard Grid Point Weight Generator output. It is obtained by setting the parameter GRDPNT
to an integer value which defines a grid point to be used as a reference point. If the integer is zero or
is not a defined grid point, the reference point is taken as the origin of the basic coordinate system.
* 2.547000E-01
* 0.000000E+00
* 0.000000E+00
* 0.000000E+00
* 0.000000E+00
* -3.820500E-01
W E I G H T
O U T P U T
Cost
SUPORT
ID
SUPORT
16
125
Notes:
ID
ID
ID
1.
2.
3.
10-24
EIGRL
SID
V1
V2
ND
EIGRL
0.1
3.2
10
Field
10
NORM
Contents
SID
V1, V2
ND
MSGLVL
MAXSET
10-25
NORM
MAX
10-26
EIGR
SID
EIGR
+3C
Field
SID
METHOD
F1
F2
NE
ND
13
SINV
1.9
15.6
NORM
POINT
32
12
10
ABC
Contents
Set identification number (integer > 0)
When METHOD = INV
or SINV
F1,F2
NE
ND
Desired
number
of
roots.
(Default is 3 NE for INV only,
integer > 0). If blank, all roots
between F1 and F2 are searched
for (SINV only).
10-27
Frequency
range
of
interest (real 0.0, F1 < F2).
If ND is not blank, F1 and
F2 are ignored. If ND is
blank, eigenvectors are
found
whose
natural
frequencies lie in the range
between F1 and F2.
Desired
number
of
eigenvectors (integer > 0).
If ND is blank or zero, the
number of eigenvectors is
determined from F1 and F2.
(Default = 0)
Contents
Method for normalizing eigenvectors, one of the BCD values,
MASS, MAX, or POINT
MASS
MAX
POINT
10-28
SECTION 11
LINEAR BUCKLING ANALYSIS
THEORY OF BUCKLING . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-1
SOLUTION OF THE EIGENVALUE PROBLEM . . . . . . . . . . . . . . . . . . . . . . . . . 11-2
SOLUTION SEQUENCES FOR BUCKLING AND
STABILITY PROBLEMS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-3
EXAMPLES OF NONLINEAR BUCKLING. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-5
RULES FOR SOL 105 BUCKLING ANALYSIS . . . . . . . . . . . . . . . . . . . . . . . . . 11-6
DATA ENTRIES FOR LINEAR BUCKLING . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-7
EIGB ENTRY . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-8
EXAMPLE - SIMPLE EULER COLUMN . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-10
REFERENCE FOR BUCKLING AND
STABILITY ANALYSIS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11-16
13
THEORY OF BUCKLING
14
(1)
(2)
11-1
(3)
(1)
i = 1, 2, ..., E
each
eigenvalue
i ,
there
is
corresponding
distinct
eigenvector { i } .
11-2
Linear buckling
SOL 106
Nonlinear buckling
11-3
Slender
Intermediate
Short
Fails by elastic
buckling. Prebuckled
deflections are small
and critical load is
reached before
material yields. This
is a Euler column.
Note:
Fails by combination of
yielding and buckling.
Prebuckled deflections are
small, but some prebuckled
stresses are inelastic.
11-4
11-5
11-6
SUBCASE 2
METHOD = N
STATSUB = i
or
EIGRL
11-7
EIGB ENTRY
1
EIGB
SID
METHOD
L1
L2
NEP
NDP
NDN
EIGB
13
SINV
0.1
2.5
NORM
+BC
10
ABC
MAX
Field
Contents
SID
METHOD
L1,L2
NEP
NDP, NDN
NORM
Point
11-8
Contents
11-9
21 in.
diameter
E
I
A
=
=
=
=
0.25 inches
30 x 106 psi
1.917E-4 in4
4.909E-2 in2
Fixed
Theoretical Solution
where Leff
11-10
z
MSC/NASTRAN Solution
Load Value on Force Entry
Pcr = 32.18 x 1.0 = 32.18 lbs
Eigenvalue
11-11
E X E C U T I V E
C O N T R O L
D E C K
PAGE 1
E C H O
ID BASIC,EXAMPLE
TIME 10
SOL 105
CEND
EULER BUCKLING, FIXED-FREE UNIFORM COLUMN
C A S E
11-12
CARD
COUNT
1
2
3
4
5
6
7
8
9
10
11
C O N T R O L
E C H O
PAGE 2
S O R T E D
11-13
CARD
COUNT
1234567891011-
.
1 ..
CBEAM
1
CBEAM
2
CBEAM
3
EIGB
14
FORCE
5
GRID
1
GRID
2
GRID
3
GRID
4
MAT1
2
PBEAM
1
ENDDATA
TOTAL COUNT=
..
3
1
1
1
SINV
4
30.E6
2
12
..
1
2
3
0.
B U L K
D A T A
..
5 ..
6 ..
2
0.
0.
3
0.
0.
4
0.
0.
1.E4
1
1
1.
-1.
0.
0.
0.
0.
7.
0.
0.
14.
0.
0.
21.
0.
0.
0.33
4.909E-21.917E-41.917E-4
PAGE 4
E C H O
..
1.
1.
1.
..
0.
123456
345
345
345
3.835E-4
..
10
E I G E N V A L U E
A N A L Y S I S
S U M M A R Y
16
PAGE 9
11-14
PAGE 10
MODE
NO.
EXTRACTION
ORDER
1
1
2
2
EIGENVALUE
3.217839E+01
3.769160E+03
RADIANS
5.672600E+00
6.139349E+01
E I G E N V A L U E S
CYCLES
9.028223E-01
9.771076E+00
GENERALIZED
MASS
5.873542E-02
1.644818E+00
GENERALIZED
STIFFNESS
1.890011E+00
6.199583E+03
PAGE 12
SUBCASE 2
3.217839E+01
R E A L
11-15
POINT ID.
1
2
3
4
TYPE
G
G
G
G
T1
0.0
0.0
0.0
0.0
T2
0.0
1.339744E-01
4.999996E-01
1.000000E+00
E I G E N V E C T O R
T3
0.0
0.0
0.0
0.0
N O .
R1
0.0
0.0
0.0
0.0
R2
0.0
0.0
0.0
0.0
R3
0.0
3.739759E-02
6.477462E-02
7.479534E-02
2.
3.
MSC/NASTRAN Verification
January 1986 Edition):
Problem
Manual
(Version
64,
5.
September 1979
February 1982
November 1985
11-16
SECTION 12
FILE MANAGEMENT SECTION
FMS OEVERVIEW. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-1
DBSETS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-2
REQUIRED DBSETS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-3
TYPICAL FILES GENERATED . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-4
RESTARTS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-5
AUTOMATIC RESTARTS (SOLDS 101-200) . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-6
RESTART EXAMPLE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-7
RESTART EXAMPLE - COLD START . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-8
COLD START F04 FILE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-9
COLD START F06 FILE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-11
RESTART EXAMPLE - DATA RECOVERY RESTART . . . . . . . . . . . . . . . . . . . . 12-14
RESTART EXAMPLE - PREFERRED METHOD . . . . . . . . . . . . . . . . . . . . . . . . 12-15
RESTART F04 FILE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-16
RESTART F06 FILE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-18
RESTART EXAMPLE SUMMARY. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-21
MANIPULATING THE DATABASE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-22
ASSIGN AND INIT STATEMENTS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-23
EXAMPLES - ASSIGN AND INIT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-24
EXPAND STATEMENT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-25
EXAMPLE - EXPAND STATEMENT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-26
EXAMPLE ASSIGNING A FORTRAN FILE . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12-27
15
FMS OVERVIEW
Provides user dynamic file allocation
Creation/attachment of database files
Creation/attachment of FORTRAN files
Provides data manipulation
Restarting
Initializing and expanding the database
Management of data on the database
Database directory prints
Provides means of attaching external databases
12-1
16
DBSETS
A DBSET is a collection of files which are used for a purpose by
the program
The MSC/PATRAN database consists of several DBSETs which are
automatically created by the program:
MASTER - Master directory - contains the description of the
database, a list of all files contained in all DBSETs used by the run,
and a "table of contents" for each DBSETs.
DBALL - Permanent data is stored here by default
USROBJ - User DMAP source - no longer automatically created
USRSOU - Complied user DMAPs - no longer automatically created
OBJSCR - Temporary scratch for DMAPs - automatically deleted at
the end of the run
SCRATCH - Temporary working files for scratch data blocks automatically deleted at the end of the run
SCR300 - Temporary working files for the modules - automatically
deleted at the end of the run
12-2
REQUIRED DBSETS
MASTER and DBALL must be retained if the database is to be used
in a subsequent restart run. If no restart is planned, they may be
deleted upon completion of the run.
USRSOU and USROBJ are only used for user DMAP functions.
They have no effect on the database or future restarts, unless you
are modifying the solutions provided by MSC to create your own
solutions.
OBJSCR, SCRATCH, and SCR300 are used for temporary file
storage and are deleted automatically upon completion of the run.
They are typically created in a scratch directory and are not seen
by the user.
For automatic deletion of MASTER, DBALL, USRSCR, and USROBJ,
use PARAM, DBALL, SCRATCH or submit with scr=yes.
12-3
12-4
RESTARTS
Rigid Format Solutions (SOLs 1 thru 12)
These solutions no longer support restarts.
Structured Solutions (SOLs 101,103,etc)
Include automatic restart capabilities
Minimal user input required for restarts
RECOMMENDED - Let the program decide what has changed in
the model and process accordingly.
12-5
12-6
RESTART EXAMPLE
A cantilver beam composed of 1000 QUAD4 elements is subjected to
gravity load of 1G downward as shown.
Fixed
Edge
100
10
1 G Gravity
Load
E = 30 106
= .3
T = .25
Mesh size 10 100
Total elements = 1,000
12-7
12-8
1****
1
MACHINE
MODEL
SUN
Sun4/75
===
M S C / N A S T R A N
DEL_CPU
SUB_DMAP/DMAP_MODULE MESSAGES
OPERATING SYSTEM
SunOS
E X E C U T I O N
MSC/NASTRAN
VERSION 67
S U M M A R Y
12-9
20:03:14
0:00
0.0
0.0
0.0
0.0
DBINIT
BGN
** CURRENT PROJECT ID = BLANK ** CURRENT VERSION ID =
0
S U M M A R Y
O F
F I L E
A S S I G N M E N T
F O R
T H E
0
ASSIGNED PHYSICAL FILE NAME (/ORIGINAL)
LOGICAL NAME
DBSET
-----------------------------------------------------./cold.MASTER
MASTER
MASTER
./cold.USROBJ
USROBJ
USROBJ
./cold.USRSOU
USRSOU
USRSOU
./cold.DBALL
DBALL
DBALL
/msc/scratch/N2123.OBJSCR
OBJSCR
OBJSCR
**** MEM FILE ****
* N/A *
SCRATCH
/msc/scratch/N2123.SCRATCH
SCRATCH
SCRATCH
/msc/scratch/N2123.SCR300
SCR300
SCRATCH
20:03:15
0:01
1.1
1.1
0.8
0.8
DBINIT
END
20:03:15
0:01
1.1
0.0
0.8
0.0
XCSA
BGN
0
S U M M A R Y
O F
F I L E
A S S I G N M E N T
F O R
T H E
0
ASSIGNED PHYSICAL FILE NAME (/ORIGINAL)
LOGICAL NAME
DBSET
-----------------------------------------------------/msc/nast67/nast67/del/SSS.MASTERA
MASTERA
MASTER
/SSS.MASTERA
/msc/nast67/nast67/del/SSS.MSCOBJ
MSCOBJ
MSCOBJ
/SSS.MSCOBJ
/msc/nast67/nast67/del/SSS.MSCSOU
MSCSOU
MSCSOU
/SSS.MSCSOU
20:03:21
0:07
7.8
6.7
3.6
2.8
XCSA
END
20:03:21
0:07
7.8
0.0
3.6
0.0
CGPI
BGN
20:03:21
0:07
7.8
0.0
3.7
0.0
CGPI
END
20:03:21
0:07
7.8
0.0
3.7
0.0
LINKER
BGN
BUILD DATE
OCT 2, 1991
=== DAY TIME ELAPSED
RUN DATE
MAR 11, 1992
I/O SEC DEL_I/O
1
P R I M A R Y
D A T A B A S E
( DBSNO
STATUS
BUFFSIZE
CLUSTER SIZE
-----------------------NEW
2049
1
NEW
2049
1
NEW
2049
1
NEW
2049
1
NEW
2049
1
NEW
NEW
2049
2049
1
1
CPU SEC
1, V67
)
TIME STAMP
-----------920311200314
920311200315
920311200316
920311200317
920311200318
920311200319
920311200320
D E L I V E R Y
D A T A B A S E
( DBSNO
2, V67
)
STATUS
BUFFSIZE
CLUSTER SIZE
TIME STAMP
----------------------------------OLD
2049
1
910905163143
OLD
2049
910905163147
OLD
2049
910905163148
1:57
1:57
33.7
33.8
0.0
0.0
69.0
69.1
0.0
0.1
SEDRCVR
SEDRCVR
218 OFP
282 DBC
BEGN
BEGN
0.0
0.0
0.0
69.2
69.3
69.3
0.2
0.0
0.0
SEDRCVR
SEDRCVR
SEDRCVR
283 DBC
284 DBC
285 DBC
BEGN
BEGN
BEGN
.
.
.
20:05:12
20:05:12
20:05:12
1:58
1:58
1:58
33.8
33.8
33.8
12-10
.
.
.
20:05:12
1:58
34.0
0.2
69.4
0.2
SESTATIC
40
EXIT
BEGN *
0
0
N A S T R A N
ID
SOL
TIME
CEND
1
E X E C U T I V E
C O N T R O L
D E C K
E C H O
BASIC, COLD
101
100
MARCH
11, 1992
0
0
12-11
C A S E
C O N T R O L
D E C K
CARD
COUNT
1
2
3
4
5
6
7
8
9
10
0
0
ECHO = NONE
TITLE = PLATE MODEL USING 1000 QUAD4S - COLD START
SUBTITLE = ONE G GRAVITY LOAD
SET 10 = 991 THRU 1000 $ ELEM AT THE FIXED END
SET 20 = 1 THRU 11 $ GRID AT THE TIP
STRESS = 10
DISP = 20
SPC = 3
LOAD = 4
BEGIN BULK
INPUT BULK DATA CARD COUNT =
TOTAL COUNT=
2119
2126
E C H O
PAGE
MARCH
11, 1992
PAGE
PAGE
POINT ID.
1
2
TYPE
G
G
T1
0.0
0.0
D I S P L A C E M E N T
T2
T3
0.0
-2.231765E+01
0.0
-2.231768E+01
V E C T O R
R1
-3.415722E-05
-3.220524E-05
R2
-2.985874E-01
-2.985862E-01
R3
0.0
0.0
.
.
.
1
0
11
G
0.0
0.0
PLATE MODEL USING 1000 QUAD4S - COLD START
ONE G GRAVITY LOAD
-2.231765E+01
3.415669E-05 -2.985874E-01
0.0
MARCH 11, 1992 MSC/NASTRAN 10/ 2/91
12-12
S T R E S S E S
ID.
991
992
I N
Q U A D R I L A T E R A L
ELEMENT
FIBRE
STRESSES IN ELEMENT COORD SYSTEM
DISTANCE
NORMAL-X
NORMAL-Y
SHEAR-XY
ANGLE
-1.250000E-01
5.023495E+03
2.443366E+04 -3.849042E+03
-79.1833
1.250000E-01
-5.023495E+03 -2.443366E+04
3.849042E+03
10.8167
-1.250000E-01
9.018255E+03
3.238159E+04 -1.627886E+03
-86.0333
1.250000E-01
-9.018255E+03 -3.238159E+04
1.627886E+03
3.9667
E L E M E N T S
( Q U A D 4 )
.
.
.
12-13
0
1
1000
-1.250000E-01
1.250000E-01
5.023495E+03
-5.023495E+03
2.443366E+04
-2.443366E+04
0
* * * END OF JOB * * *
3.849042E+03
-3.849042E+03
79.1833
-10.8167
MARCH
2.516907E+04
-4.288091E+03
11, 1992
4.288091E+03
-2.516907E+04
2.332257E+04
2.332257E+04
PAGE
11
Submittal
nastran run2a dbs=cold
12-14
Submittal
nastran run2b
12-15
MODEL
Sun4/75
===
M S C / N A S T R A N
SUB_DMAP/DMAP_MODULE MESSAGES
OPERATING SYSTEM
MSC/NASTRAN
SunOS
VERSION 67
E X E C U T I O N
S U M M A R Y
BUILD DATE
OCT 2, 1991
=== DAY TIME
12-16
20:05:43
0:00
0.0
0.0
0.0
0.0
DBINIT
BGN
** CURRENT PROJECT ID = BLANK ** CURRENT VERSION ID =
2
0
S U M M A R Y
O F
F I L E
A S S I G N M E N T
F O R
T H E
P R I M A R Y
0
ASSIGNED PHYSICAL FILE NAME (/ORIGINAL)
LOGICAL NAME
DBSET
STATUS
----------------------------------------------------------cold.MASTER
MASTER
MASTER
OLD
/./cold.MASTER
cold.USROBJ
USROBJ
USROBJ
OLD
/./cold.USROBJ
cold.USRSOU
USRSOU
USRSOU
OLD
/./cold.USRSOU
cold.DBALL
DBALL
DBALL
OLD
/./cold.DBALL
/msc/scratch/N2169.OBJSCR
OBJSCR
OBJSCR
NEW
**** MEM FILE ****
* N/A *
SCRATCH
/msc/scratch/N2169.SCRATCH
SCRATCH
SCRATCH
NEW
/msc/scratch/N2169.SCR300
SCR300
SCRATCH
NEW
20:05:45
0:02
0.5
0.5
0.4
0.4
DBINIT
END
20:05:45
0:02
0.5
0.0
0.4
0.0
XCSA
BGN
0
S U M M A R Y
O F
F I L E
A S S I G N M E N T
F O R
T H E
D E L I V E R
0
ASSIGNED PHYSICAL FILE NAME (/ORIGINAL)
LOGICAL NAME
DBSET
STATUS
----------------------------------------------------------/msc/nast67/nast67/del/SSS.MASTERA
MASTERA
MASTER
OLD
/SSS.MASTERA
/msc/nast67/nast67/del/SSS.MSCOBJ
MSCOBJ
MSCOBJ
OLD
/SSS.MSCOBJ
/msc/nast67/nast67/del/SSS.MSCSOU
MSCSOU
MSCSOU
OLD
/SSS.MSCSOU
RUN DATE
MAR 11, 1992
ELAPSED I/O SEC
DEL_I/O
CPU SEC
D A T A B A S E
( DBSNO
1, V67
)
BUFFSIZE
CLUSTER SIZE
TIME STAMP
-----------------------------2049
1
920311200314
2049
920311200315
2049
920311200316
2049
920311200317
2049
920311200544
2049
2049
1
1
920311200545
920311200546
Y
D A T A B A S E
( DBSNO
2, V67
)
BUFFSIZE
CLUSTER SIZE
TIME STAMP
-----------------------------2049
1
910905163143
2049
910905163147
2049
910905163148
20:05:48
0:05
6.5
6.0
20:05:48
0:05
6.5
0.0
20:05:48
0:05
6.5
0.0
20:05:48
0:05
6.5
0.0
20:05:50
0:07
7.7
1.2
20:05:50
0:07
7.7
0.0
20:05:52
0:09
12.9
5.2
0
** MASTER DIRECTORIES ARE LOADED IN
20:05:52
20:05:52
20:05:52
0:09
0:09
0:09
13.2
13.2
13.4
2.8
2.3
XCSA
END
2.8
0.0
CGPI
BGN
2.8
0.1
CGPI
END
2.8
0.0
LINKER
BGN
3.8
1.0
LINKER
END
3.8
0.0
DBERST
BGN
4.8
1.0
DBERST
END
MEMORY. DAY TIME ELAPSED I/O SEC
0.3
0.1
0.1
5.0
5.0
5.1
0.2
0.0
0.1
SESTATIC
SUPER1
IFPL
0.0
0.2
17.1
17.2
0.0
0.1
SEDRCVR
SESTATIC
16
16
15
DEL_I/O
SUPER1
IFPL
IFP1
CPU SEC
DEL_CPU
SUB_DMAP/DMAP_MODULE MESSAGES
BEGN
BEGN
BEGN *
.
.
.
12-17
20:06:14
20:06:14
0:31
0:31
21.0
21.2
285 DBC
40 EXIT
BEGN
BEGN *
MARCH
0
0
N A S T R A N
F I L E
M A N A G E M E N T
S E C T I O N
RESTART
0*** USER INFORMATION MESSAGE 736 (RDREST)
THE RESTART VERSION ID IS NOT DEFINED ON THE RESTART COMMAND.
A RESTART OF LAST VERSION ID IS ASSUMED.
0*** USER INFORMATION MESSAGE 1144 (RSPRVR)
LAST VERSION CREATED IS BEING USED FOR THIS RESTART JOB.
LAST VERSION =
1
PROJECT = BLANK
1
12-18
0
0
N A S T R A N
ID
SOL
TIME
CEND
1
BASIC, COLD
101
100
E X E C U T I V E
C O N T R O L
D E C K
11, 1992
PAGE
11, 1992
PAGE
E C H O
MARCH
E C H O
C A S E
12-19
CARD
COUNT
1
2
3
4
5
6
7
8
9
10
C O N T R O L
MARCH
D E C K
11, 1992
E C H O
ECHO = NONE
TITLE = PLATE MODEL USING 1000 QUAD4S - RESTART FOR ADD OUTPUT
SUBTITLE = ONE G GRAVITY LOAD
SET 10 = 501 THRU 511 $ ELEM AT THE FIXED END
SET 20 = 551 THRU 561 $ GRIDS AT THE CENTER
STRESS = 10
DISP = 20
SPC = 3
LOAD = 4
BEGIN BULK
INPUT BULK DATA CARD COUNT =
1
TOTAL COUNT= 2119
0
0
RESTART F06 FILE (Cont.)
.
.
.
PAGE
TYPE
G
T1
0.0
T2
0.0
V E C T O R
T3
-7.857040E+00
R1
-3.428952E-03
R2
-2.609812E-01
R3
0.0
.
.
1
0
1
561
G
0.0
0.0
-7.857040E+00
PLATE MODEL USING 1000 QUAD4S - RESTART FOR ADD OUTPUT
ONE G GRAVITY LOAD
3.428952E-03 -2.609812E-01
0.0
MARCH 11, 1992 MSC/NASTRAN 10/ 2/91
MARCH
11, 1992
PAGE
PAGE
12-20
ELEMENT
ID.
0
501
S T R E S S E S
I N
Q U A D R I L A T E R A L
STRESSES IN ELEMENT COORD SYSTEM
NORMAL-X
NORMAL-Y
SHEAR-XY
-7.436593E+00
8.650283E+03
3.532764E+02
7.436593E+00 -8.650283E+03 -3.532764E+02
FIBRE
DISTANCE
-1.250000E-01
1.250000E-01
E L E M E N T S
( Q U A D 4 )
PRINCIPAL STRESSES (ZERO SHEAR)
ANGLE
MAJOR
MINOR
87.6672
8.664674E+03
-2.182804E+01
-2.3328
2.182804E+01
-8.664674E+03
VON MISES
8.675609E+03
8.675609E+03
87.7124
9.010879E+03
-2.182634E+01
-2.2876
2.182634E+01
-9.010879E+03
MARCH 11, 1992 MSC/NASTRAN 10/ 2/91
9.021812E+03
9.021812E+03
PAGE
10
.
.
0
1
0
1
511
-1.250000E-01
-7.434482E+00
8.996487E+03
3.602642E+02
1.250000E-01
7.434482E+00 -8.996487E+03 -3.602642E+02
PLATE MODEL USING 1000 QUAD4S - RESTART FOR ADD OUTPUT
ONE G GRAVITY LOAD
* * * END OF JOB * * *
12-21
12-22
12-23
12-24
EXPAND STATEMENT
Allows the user to add a new files to an existing DBset
The ASSIGN statement should be used in conjuction with the
EXPAND statement.
Format:
ASSIGN dbset_member_name=phyiscal_file_name
EXPAND dbset_name LOGICAL=(dbset_member_name(p),...)
where p = number of blocks allocated for the member
12-25
12-26
12-27
12
2KGG
1
1
7
3.000000142E-01 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00
0.000000000E+00-3.000000142E-01
2
2
11
3.600000171E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 1.800000085E+00
0.000000000E+00-3.600000171E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00
1.800000085E+00
3
3
9
3.600000171E+00 0.000000000E+00-1.800000085E+00 0.000000000E+00 0.000000000E+00
0.000000000E+00-3.600000171E+00 0.000000000E+00-1.800000085E+00
4
4
7
2.307692370E-01 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00
0.000000000E+00-2.307692370E-01
5
3
9
-1.800000085E+00 0.000000000E+00 1.200000057E+00 0.000000000E+00 0.000000000E+00
0.000000000E+00 1.800000085E+00 0.000000000E+00 6.000000285E-01
6
2
11
1.800000085E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 1.200000057E+00
0.000000000E+00-1.800000085E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00
6.000000285E-01
7
1
7
-3.000000142E-01 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00
0.000000000E+00 3.000000142E-01
8
2
11
-3.600000171E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00-1.800000085E+00
0.000000000E+00 3.600000171E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00
-1.800000085E+00
9
3
9
-3.600000171E+00 0.000000000E+00 1.800000085E+00 0.000000000E+00 0.000000000E+00
0.000000000E+00 3.600000171E+00 0.000000000E+00 1.800000085E+00
10
4
7
-2.307692370E-01 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00
0.000000000E+00 2.307692370E-01
11
3
9
-1.800000085E+00 0.000000000E+00 6.000000285E-01 0.000000000E+00 0.000000000E+00
0.000000000E+00 1.800000085E+00 0.000000000E+00 1.200000057E+00
12
2
11
1.800000085E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 6.000000285E-01
0.000000000E+00-1.800000085E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00
1.200000057E+00
13
1
1
1.974999905E+00
12
12
6
2MGG
1
1
1
3.660000232E-07
2
2
1
3.660000232E-07
3
3
1
3.660000232E-07
7
7
1
3.660000232E-07
8
8
1
3.660000232E-07
9
9
1
3.660000232E-07
13
1
1
2.979723513E-01
12-28
APPENDIX A
MSC/NASTRAN DOCUMENTATION
DOCUMENTATION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-1
ADDITIONAL DOCUMENTATION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . A-4
17
DOCUMENTATION
18
A-1
installing,
DOCUMENTATION
19
A-2
DOCUMENTATION
20
MSC/NASTRAN Documentation
Reference Manual
DMAP Module Dictionary
Quick Reference Guide
Users Guides
Online Encyclopedia CD-ROM
Bibliography
Common Questions & Answers
Configuration and Operations Guide
Release Guide
A-3
pertaining
to
ADDITIONAL DOCUMENTATION
Older User Documentation Available
MSC/NASTRAN Handbook for Superelement Analysis
MSC/NASTRAN Handbook for Nonlinear Analysis
MSC/NASTRAN Demonstration Problem Manual
Shows illustrative examples with MSC/NASTRAN input and output.
MSC/NASTRAN Verification Probblem Manual
Provides examples for which there are known solutions and
compares MSC/NASTRAN output with theory. Examples using
most of the Unstructured Solution Sequences are included.
The NASTRAN Theoretical Manual
Explains the theoretical basis of the elements and algorithms. The
numerical algorithm material is superseded by the MSC/NASTRAN
Handbook for Numerical Methods and the MSC/NASTRAN
Numerical Methods Users Guide.
MSC/NASTRAN Programmers Manual
Describes the input and output for the modules. This manual is not
recommended for general use but is useful for advanced DMAP
development.
All of the documents listed above can be ordered from your local
MSC office or representative.
A-4
ADDITIONAL DOCUMENTATION
MSC Corporate Web Site
The MacNeal-Schwendler Corporations Web site provides several
sources of information that can assist you in running
MSC/NASTRAN and our other products. MSCs Web address is
http://www.msc.com
Here you can find out what is new with MSC, read white papers on
the use of MSC products, download technical papers from previous
Users Conference, review the minutes of the Technical Forum,
order documentation from the MSC Bookstore, obtain a schedule
of training courses, share feedback and suggestions interactively
with other users, subscribe to MSCs corporate newsletter, and
even download software patches and utilities. The sssalter library
and error list will also be posted on the Web in the future.
A-5
A-6
APPENDIX B
STRUCTURE PLOTS
INTRODUCTION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . B-1
PLOT ENTRIES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . B-3
VIEWING OPTIONS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . B-4
EXAMPLE - PLOT INPUT FILE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . B-6
PLOT SET 1 LABEL BOTH . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . B-7
PLOT SET 1 SHRINK, LABEL GSPC. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . B-8
PLOT STATIC DEFORMATION 0 SET . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . B-9
21
INTRODUCTION
22
These plots can be viewed using the batch plotter routine supplied
with MSC/NASTRAN: PLOTPS.
PLOTPS converts the file createed by MSC/NASTRAN into a
postscript format which may be printed on any postscript
printer or viewed in any postscript viewer.
Instructions for using PLOTPS are given in Section 6.7 of the
Configuration and Operations Guide
B-1
INTRODUCTION (Cont.)
23
B-2
PLOT ENTRIES
A request for structure plots is made in the Case Control beginning
with the OUTPUT(PLOT) commands and ending with the BEGIN
BULK command.
Only elements can be plotted, and grid points not associated with
elements cannot be plotted. Grid points can be connected with
PLOTEL elements for plotting purposes.
Required commands for structure plots:
OUTPUT(PLOT)
SET
FIND
PLOT
Plot request entries are free-field format, and subject to the rules
given in Section 4 of the MSC/NASTRAN Quick Refernce Guide.
Order of plot package data entries:
Set entries can be defined anywhere in the plot package.
Parameters describing plot characteristics are evaluated every
time a PLOT or FIND entry is encountered. To minimize mistakes,
pay attention to the order of the plot entries.
B-3
VIEWING OPTIONS
The plot frame can be visualized as occupying the S-T plane of the
viewing coordinate system R, S, T. The AXES entry associates the R,
S, T with X, Y, Z. Three angles (, , and ) define the angular
relationship between R, S, T and the X, Y, Z axes of the object (in the
basic coordinate system).
T
B-4
, ,
(default
for
orthographic
B-5
Plot titles
Up to four lines of title information are printed in the lower left-hand
corner of each plot.
TITLE
SUBTITLE Case Control
LABEL
PLOT TYPE (e.g., deformed plot)
B-6
11
12
504
6
13
604
104
14
704
204
2
304
804
3
404
4
5
B-7
10
15
11
12
13
7
1
14
8
15
2
9
3
10
4
5
B-8
B-9
B-10
APPENDIX C
MSC/NASTRAN MODEL CHECKOUT
Presented at the
MSC/NASTRAN USERS CONFERENCE
MARCH 20-21, 1986